Actin droplet machine
Andrew Adamatzky, Florian Huber, Joerg Schnauss

TL;DR
This paper introduces a computer model of an actin droplet network that processes binary information through traveling excitations, potentially serving as a basis for future biopolymer-based parallel computing systems.
Contribution
It presents a novel design and analysis of an actin droplet machine capable of implementing binary string mappings and state transitions for biopolymer computing.
Findings
The actin droplet network can generate and record impulses as binary responses.
The machine's state transitions can be mapped and analyzed through graph models.
Potential application as an elementary processor in biopolymer-based parallel computers.
Abstract
The actin droplet machine is a computer model of a three-dimensional network of actin bundles developed in a droplet of a physiological solution, which implements mappings of sets of binary strings. The actin bundle network is conductive to travelling excitations, i.e. impulses. The machine is interfaced with an arbitrary selected set of k electrodes through which stimuli, binary strings of length k represented by impulses generated on the electrodes, are applied and responses are recorded. The responses are recorded in a form of impulses and then converted to binary strings. The machine's state is a binary string of length k: if there is an impulse recorded on the th electrode, there is a `1' in the i-th position of the string, and `0' otherwise. We present a design of the machine and analyse its state transition graphs. We envisage that actin droplet machines could form an…
| 1 | 369 | 567 | 6 |
|---|---|---|---|
| 2 | 509 | 580 | 10 |
| 3 | 631 | 590 | 10 |
| 4 | 382 | 322 | 12 |
| 5 | 533 | 331 | 23 |
| 6 | 626 | 463 | 7 |
| 7 | 358 | 676 | 22 |
| 8 | 369 | 424 | 7 |
| 9 | 572 | 691 | 17 |
| 10 | 705 | 394 | 17 |
| 1 | 3 | 8, 9, 1, |
| 2 | 3 | 16, 32, 8, |
| 3 | 3 | 1, 16, 32, |
| 4 | 3 | 8, 1, 16, |
| 5 | 3 | 1, 8, 16, |
| 6 | 3 | 16, 8, 1, |
| 7 | 4 | 8, 1, 16, 4, |
| 8 | 4 | 1, 16, 8, 5, |
| 9 | 5 | 16, 1, 8, 4, 5, |
| 10 | 4 | 16, 1, 8, 4, |
| 11 | 5 | 8, 1, 16, 20, 4, |
| 12 | 4 | 1, 16, 8, 20, |
| 13 | 6 | 16, 8, 1, 17, 4, 20, |
| 14 | 8 | 8, 16, 17, 4, 20, 1, 32, 2, |
| 15 | 8 | 1, 16, 8, 4, 2, 10, 20, 32, |
| 16 | 6 | 16, 4, 8, 1, 10, 32, |
| 17 | 5 | 16, 1, 4, 8, 9, |
| 18 | 7 | 8, 16, 4, 1, 17, 10, 9, |
| 19 | 6 | 1, 8, 16, 17, 4, 10, |
| 20 | 8 | 16, 1, 8, 17, 4, 24, 10, 2, |
| 21 | 9 | 8, 16, 1, 17, 32, 24, 9, 4, 10, |
| 22 | 6 | 16, 1, 8, 32, 9, 4, |
| 23 | 7 | 8, 1, 16, 4, 32, 9, 17, |
| 24 | 6 | 1, 16, 17, 4, 32, 8, |
| 25 | 7 | 16, 1, 8, 4, 17, 32, 9, |
| 26 | 6 | 8, 16, 4, 12, 1, 17, |
| 27 | 6 | 1, 8, 16, 4, 17, 32, |
| 28 | 6 | 16, 8, 4, 1, 24, 32, |
| 29 | 7 | 8, 1, 4, 16, 12, 24, 32, |
| 30 | 7 | 16, 1, 8, 4, 17, 2, 32, |
| 31 | 9 | 8, 1, 24, 16, 12, 4, 2, 17, 32, |
| 32 | 7 | 1, 16, 8, 24, 17, 2, 40, |
| 33 | 9 | 16, 8, 1, 4, 40, 17, 24, 32, 2, |
| 34 | 7 | 8, 1, 16, 24, 40, 4, 32, |
| 35 | 6 | 1, 16, 8, 4, 24, 2, |
| 36 | 6 | 16, 8, 1, 17, 4, 32, |
| 37 | 7 | 8, 16, 17, 4, 1, 40, 2, |
| 38 | 7 | 1, 8, 16, 17, 4, 24, 2, |
| 39 | 7 | 16, 1, 8, 17, 9, 4, 2, |
| 40 | 7 | 8, 16, 4, 1, 24, 40, 2, |
| 41 | 10 | 1, 8, 16, 9, 17, 4, 18, 24, 40, 2, |
| 42 | 8 | 16, 1, 8, 4, 18, 33, 40, 24, |
| 43 | 9 | 8, 1, 16, 4, 24, 33, 18, 32, 34, |
| 44 | 9 | 1, 16, 8, 4, 17, 33, 24, 32, 40, |
| 45 | 7 | 16, 8, 4, 1, 12, 24, 34, |
| 46 | 7 | 8, 1, 16, 4, 24, 18, 34, |
| 47 | 5 | 1, 16, 8, 4, 33, |
| 48 | 5 | 16, 8, 1, 4, 17, |
| 49 | 8 | 8, 1, 16, 4, 20, 32, 24, 19, |
| 50 | 6 | 1, 16, 8, 4, 17, 32, |
| 51 | 8 | 16, 8, 1, 4, 17, 32, 41, 19, |
| 52 | 9 | 8, 16, 4, 1, 32, 33, 41, 2, 19, |
| 53 | 10 | 1, 8, 16, 4, 20, 10, 2, 41, 32, 19, |
| 54 | 9 | 16, 1, 8, 5, 17, 4, 2, 32, 19, |
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Actin droplet machine
Andrew Adamatzky
Unconventional Computing Laboratory, Department of Computer Science, University of the West of England, Bristol, UK
Florian Huber
Netherlands eScience Center, Science Park 140, 1098 XG Amsterdam, The Netherlands
Jörg Schnauß
Soft Matter Physics Division, Peter Debye Institute for Soft Matter Physics, Faculty of Physics and Earth Science, Leipzig University, Germany & Fraunhofer Institute for Cell Therapy and Immunology (IZI), DNA Nanodevices Group, Leipzig, Germany
Abstract
The actin droplet machine is a computer model of a three-dimensional network of actin bundles developed in a droplet of a physiological solution, which implements mappings of sets of binary strings. The actin bundle network is conductive to travelling excitations, i.e. impulses. The machine is interfaced with an arbitrary selected set of electrodes through which stimuli, binary strings of length represented by impulses generated on the electrodes, are applied and responses are recorded. The responses are recorded in a form of impulses and then converted to binary strings. The machine’s state is a binary string of length : if there is an impulse recorded on the th electrode, there is a ‘1’ in the th position of the string, and ‘0’ otherwise. We present a design of the machine and analyse its state transition graphs. We envisage that actin droplet machines could form an elementary processor of future massive parallel computers made from biopolymers.
Keywords: actin network, computing, waves, logical gates, finite state machine, automata
1 Introduction
Actin is a protein presented in forms of monomeric, globular actin (G-actin) and filamentous actin (F-actin) [34, 18, 38]. G-actin polymerises into filamentous actin forming a double helical structure [35, 19, 13]. The filaments can be further arranged into bundles by various different mechanisms such as crowding effects, cross-linking or counterion condensation [46, 45, 7, 29, 30, 37, 36, 14, 15]. The bundles are conductive to travelling localisations — defects, ionic waves, solitons [41, 42, 43, 40, 22, 44, 25, 27, 26, 17]. By interpreting presence or absence of a travelling localisation at a given site of the network at a given time step, we can implement a logical function. This approach was comprehensively developed and successfully tested on chemical systems in the framework of collision-based computing [1, 16, 8, 39, 49, 2].
Our approach — computing with excitation waves propagating on overall ‘density’ of the conductive material — has previously been presented by us in [5]. As conductive material we looked at networks of actin bundles which were arranged by crowding effects without the need of additional accessory proteins [29, 30]. We demonstrated how to discover logical gates on a two-dimensional slice of the actin bundle network by representing Boolean inputs and outputs as spikes of the network activity.
In the present paper we develop a novel concept and computer modelling implementation of the actin network machine, which implements a mapping , where is a number of electrodes, and ‘1’ signifies a presence of an impulse on the electrode and ‘0’ the absence. At a higher level, the machine acts as a finite state machine, at the lower level a structure of the mapping is determined by interactions of impulses propagating on the three-dimensional network of actin bundles.
We also offer an alternative to a numerical integration used in [5]: an automaton model of a three-dimensional actin network. There is a substantial body of evidence confirming that automaton models are sufficient and appropriate discrete tools for modelling dynamics of spatially extended non-linear excitable media [21, 10, 47], propagation [20], action potential [48, 6], electrical pulses in the heart [28, 9, 33]. A major advantage of automata is that they require less computational resources than typical numerical integration approaches.
The paper is structured as follows. Our modelling approach is described in detail in Sect. 2. This includes a representation of a three-dimensional actin bundle network (Subsect. 2.1), a structure of an automaton model to simulate propagation of impulses on the actin bundle network (Subsect. 2.2), and an interface with the actin network (Subsect. 2.3). In Sect. 3 we analyse dependencies of a number of Boolean gates implemented in the network on an excitation threshold and refractory period. Thus, we justify the selection of these parameters for the construction of the actin machine. The actin droplet machine is designed and analysed in Sect. 4. Section 5 discusses the results in a context of cytoskeleton computing and outlines directions for future research.
2 Methods
The overall approach is the following: we simulate the actin bundle network using three-dimensional arrays of finite-state machines, cellular automata. We select several domains of the network and assign them as inputs and outputs. We represent Boolean logic values with spikes of electrical activity, which are schematically represented as a virtual experiment in Fig. 1. We stimulate the network with all possible configurations of input strings and record spikes on the outputs. Based on the mapping of configurations of input spikes to output spikes, we reconstruct logical functions implemented by the network. In our design of the actin droplet machine we consider outputs recorded on all electrodes at a given time step as a binary string and then represent the actin droplet machine as a finite-state machine whose states are binary strings of a given length.
2.1 Three-dimensional actin network
As a template for our actin droplet machine we used an actual three-dimensional actin bundle network produced in laboratory experiments with purified proteins (Fig. 2). The underlying experimental method was shown to reliably produce regularly spaced bundle networks from homogeneous filament solutions inside small isolated droplets in the absence of molecular motor-driven processes or other accessory proteins [15]. These structures effectively form very stable and long-living three-dimensional networks, which can be readily imaged with confocal microscopy resulting in stacks of optical two-dimensional slices (Fig. 2). Dimensions of the network are the following: size along coordinate is 225 m (width), along coordinate is 222 m (height), along coordinate is 112 m (depth), voxel width is 0.22 m, height 0.22 m and depth 4 m.
Original image: , , where , , are RGB values of the element at , was converted to a conductive matrix as follows: if , and . The conductive matrices are shown in Fig. 3. The 3D conductive matrix is compressed along -axis to reduce consumption of computational resources, scenario of the non-compressed matrix will be considered in future papers.
2.2 Automaton model
To model activity of an actin bundle network we represent it as an automaton . is a set of voxels, or a conductive matrix defined in Sect. 2.1. Each voxel takes states from the set , excited (), refractory (), resting () and is complemented by a counter to handle the temporal decay of the refractory state. Following discrete time steps, each voxel updates its state depending on its current state and the states of its neighbourhood , where is an Euclidean distance between voxels and ; is a neighbourhood radius. is an excitation threshold and is refractory delay. All voxels update their states in parallel and by the same rule:
[TABLE]
[TABLE]
Every resting () voxel of excites () at the moment if a number of its excited neighbours at the moment , , exceeds a threshold . An excited voxel takes the refractory state at the next time step and at the same moment a counter of refractory state is set to the refractory delay . The counter is decremented, at each iteration until it becomes 0. When the counter becomes zero the voxel returns to the resting state . For all results shown in this manuscript, the neighbourhood radius was set to . Choices of and are considered in Sect. 3.
2.3 Interfacing with the network
To stimulate the network and to record activity of the network we assigned several domains of as electrodes. We calculated a potential at an electrode location as , where is an Euclidean distance between sites and in 3D space. We have chosen an electrode radius of voxels and conducted two families of experiments with two configurations of electrodes.
In the first family of experiments we studied frequencies of two-input-one-output Boolean functions implementable in the network. We used ten electrodes, their coordinates are listed in Tab. 2 and a configuration is shown in Fig. 4. Electrodes representing input and representing input are the input electrodes, all others are output electrodes representing outputs . Results are presented in Sect. 3. In the second family of experiments we used six electrodes (Tab. 2 and Fig. 4). All electrodes were considered as inputs during stimulation and outputs during recording of the network activity.
Exemplary snapshots of excitation dynamics on the network are shown in Fig. 5. Domains corresponding to the two electrodes and (Tab. 2 and Fig. 4) have been excited (Fig. 5(a)). The excitation wave fronts propagates away from and (Fig. 5(b)). The fronts traverse the whole breadth of the network (Fig. 5(c)). Due to the presence of circular conductive paths in the network, the repetitive patterns of activity emerge (Fig. 5(d)). Videos of the experiments can be found in http://doi.org/10.5281/zenodo.2649293.
3 Frequencies of gates
To map dynamics of the network onto sets of gates, we undertook the following trials of stimulation
fixed refractory delay and excitation threshold , 2. 2.
fixed excitation threshold , and refractory delay .
For each combination we counted numbers of gates or, and, xor, not-and, and-not and select. We found that in overall a total number of gates realised by the network decreases with increase of (Fig. 6). The function is non-linear and could be adequately described by a five degree polynomial. The function reaches its maximal value at (Fig. 6). or gates are most commonly realised at , and gates at and gates at as well as (Fig. 6. A number of and-not gates implemented by the network reaches its highest value at then drops sharply after (Fig. 6). not-and gates are more common at , while select has its peak at and select at (Fig. 6). A total number of gates realised in the network with the excitability threshold fixed to decreases with the increase of . Oscillations of are visible at (Fig. 6). The three highest values of are achieved at and . Let us look now at the dependence of the numbers of or, and and xor gates of the refractory delay in Fig. 6. The number of or gates increases with increasing from 10 to 15, but then drops substantially at to reach its maximum at . Numbers of gates and and xor behave similarly to each other. They both have a pronounced peak at (Fig. 6). Thus, to maximise a number of logical gates produced and their diversity we selected and for our construction of the actin droplet machine.
4 Actin droplet machine
An actin droplet machine is defined as a tuple , where is an actin network automaton, defined in Sect. 2.2, is a number of electrodes, is a configuration of electrodes, , is a state-transition function that implements a mapping between sets of all possible configurations of binary strings of length . In the experiments reported here .
In our experiments we have chosen six electrodes, their locations are shown in Fig. 4 and exact coordinates in Tab, 2. Thus, and the machine has 64 states. We represent the inputs and the machine states in decimal encoding. Spikes detected in response to every input from are shown in Fig. 7.
Global transition graphs of for selected inputs are shown in Fig. 8. Nodes of the graphs are states of , edges show transitions between the states. These directed graphs are defined as follows. There is an edge from node to node if there is such that and .
Let us now define a weighted global transition graph , where is a set of nodes (isomorphic to the ), and is a set of edges, and weighting function assigning a number of a unit interval to each edge. Let and then a normalised weight is calculated as , with takes value ‘1’ when the conditions are true and ‘0’ otherwise. In words, is a number of transitions from to observed in the evolution of for all possible inputs from during time interval normalised by a total number of transition from to all other nodes. The graph is visualised in Fig. 9. Nodes which have predecessors are 1–6, 8–10, 12, 16–21, 24, 25, 28, 32–34, 36–38, 40, 41, 44, 48–50, 52, 53, 56. Nodes without predecessors are 7, 11, 13–15, 22, 23, 26, 27, 29–31, 35, 39, 42, 43, 45–47, 51, 54, 55, 57–63.
Let us convert to an acyclic non-weighted graph of more likely transitions , where if . That is for each node we select an outgoing edge with maximum weight. The graph is a tree, see Fig. 9. Most states apart of 1, 2, 4, 8, 16, 20, 32 are Garden-of-Eden configurations, which have no predecessors. Indegrees of not-Garden-of-Eden nodes are . There is one fixed point, the state 1, corresponding to the situation when a spike is recorded only on electrode ; it has no successors.
By analysing we can characterise a richness of ’s responses to input stimuli. We define a richness as a number of different states over all inputs, as shown in Tab. 3, and distribution in Fig. 11. A number of states produced increases from under five for beginning of evolution and then reaches circa seven states in average. Oscillations around this value are seen in (Fig. 11). Figure 11 shows a number of different nodes, generated in evolution of , stimulated by a given input. There is below fifteen different states found in the evolution in responses to inputs 1 to 21 (21 corresponds to binary input string 010101); then a number of different nodes stay around 25. The diagram Fig. 11 shows how many inputs might lead to a given state/node of . Some of the states/nodes are seen to be Garden-of-Eden configurations (nodes without predecessors) and thus could not be generated by stimulating by sequences from .
Assume is a set of temporal moments when the machine responded at least to one input string with a non-zero state. Configurations at each transition can be considered as outputs representing the function . As we can see in Tab. 3, transitions at and correspond to the highest number of different binary strings . The graph corresponding to at is shown in Fig. 10 and is not connected. The small component consists of fixed point 40 (string ‘101000’) with two leafs 39 (‘100111’) and 38 (‘100110’). The largest component has a tree structure at large, with cycle 2 (‘000010’) – 1 (‘000001’) as a root. Other nodes with most predecessors are 8 (‘001000’), 16 (‘010000’), and 18 (‘010010’).
From the transitions we can reconstruct Boolean functions realised at each of six electrodes (the functions are minimised and represented in a disjunctive normal form):
:
:
:
:
:
:
5 Discussion
Early concepts of sub-cellular computing on cytoskeleton networks as microtubule automata [12, 24, 11] and information processing in actin-tubulin networks [23] did not specify what type of ‘computation’ or ‘information processing’ the cytoskeleton networks could execute and how exactly they do this. We implemented several concrete implementations of logical gates and functions on a single actin filament [32] and on an intersection of several actin filaments [31] via collisions between solitons. We also used a reservoir-computing-like approach to discover functions on a single actin unit [3] and filament [4]. Later, we realised that it might be unrealistic to expect someone to initiate and record a travelling localisations (solitons, impulses) on a single actin filament. Therefore, we developed a numerical model of spikes propagating on a network of actin filament bundles and demonstrated that such a network can implement Boolean gates [5].
In present paper, we reconsidered the whole idea of the information processing on actin networks and designed an actin droplet machine. The machine is a model of a three-dimensional network, based on an experimental network developed in a droplet, which executes mapping of a space of binary strings of length on itself. The machine acts as a finite state machine, which behaviour at a low level is governed by localisations travelling along the networks and interacting with each other. By focusing on a single element of a string, i.e. a single location of an electrode, we can reconstruct functions with arguments, as we have exemplified at the end of the Sect. 4. Exact structure of each -ary function is determined by , which, in turn, is determined by the exact architecture of a three-dimensional actin network and a configuration of electrodes.
Thus, potential future directions could be in detailed analysis of possible architectures of actin networks developed in laboratory experiments and evaluation on how far an exact configuration of electrodes affects a structure of mapping and corresponding distribution of functions implementable by the actin droplet machine. The ultimate goal would be to implement actin droplet machines in laboratory experiments and to cascade several machines into a multi-processors computing architecture.
Author contributions statements
A.A., F.H., J.S. undertook the research and wrote the manuscript.
Competing interests
The authors declare that they have no competing interests.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] A. Adamatzky. Collision-based computing in biopolymers and their automata models. International Journal of Modern Physics C , 11:1321–1346, 2000.
- 2[2] Andrew Adamatzky, editor. Collision-based Computing . Springer, 2002.
- 3[3] Andrew Adamatzky. Logical gates in actin monomer. Scientific reports , 7(1):11755, 2017.
- 4[4] Andrew Adamatzky. On discovering functions in actin filament automata. ar Xiv preprint ar Xiv:1807.06352 , 2018.
- 5[5] Andrew Adamatzky, Florian Huber, and Jörg Schnauß. Computing on actin filament bundles. ar Xiv preprint ar Xiv:1903.10186 , 2019.
- 6[6] Felipe Alonso Atienza, Jesús Requena Carrión, Arcadi García Alberola, José L Rojo Álvarez, Juan J Sánchez Muñoz, Juan Martínez Sánchez, and Mariano Valdés Chávarri. A probabilistic model of cardiac electrical activity based on a cellular automata system. Revista Española de Cardiología (English Edition) , 58(1):41–47, 2005.
- 7[7] James R Bartles. Parallel actin bundles and their multiple actin-bundling proteins. Current opinion in cell biology , 12(1):72–78, 2000.
- 8[8] Steve Blair and Kelvin Wagner. Gated logic with optical solitons. In Collision-based computing , pages 355–380. Springer, 2002.
