Chemical Memory with Discrete Turing Patterns Appearing in the Glycolytic Reaction

Memory is an essential element in information processing devices. We investigated a network formed by just three interacting nodes representing continuously stirred tank reactors (CSTRs) in which the glycolytic reaction proceeds as a potential realization of a chemical memory unit. Our study is based on the 2-variable computational model of the reaction. The model parameters were selected such that the system has a stable limit cycle and several distinct, discrete Turing patterns characterized by stationary concentrations at the nodes. In our interpretation, oscillations represent a blank memory unit, and Turing patterns code information. The considered memory can preserve information on one of six different symbols. The time evolution of the nodes was individually controlled by the inflow of ATP. We demonstrate that information can be written with a simple and short perturbation of the inflow. The perturbation applies to only one or two nodes, and it is symbol specific. The memory can be erased with identical inflow perturbation applied to all nodes. The presented idea of pattern-coded memory applies to other reaction networks that allow for discrete Turing patterns. Moreover, it hints at the experimental realization of memory in a simple system with the glycolytic reaction.

Identification of the best medium for experiments on chemical computation with Belousov–Zhabotinsky reaction and ferroin-loaded Dowex beads

Our study is focused on identification of the best medium for future experiments on information processing with Belousov–Zhabotinsky reaction proceeding in Dowex beads with immobilized catalyst inside. The optimum medium should be characterized by long and stable nonlinear behavior, mechanical stability and should allow for control with electric potential. We considered different types of Dowex ion-exchange resins, bead distributions and various initial concentrations of substrates: malonic acid and 1,4-cyclohexanedione. The electric potential on platinum electrodes, stabilized by a potentiostat is used to control medium evolution. A negative electric potential generates activator species HBrO2 on the working electrode according to the reaction: BrO3−  + 2e−  + 3H+  → HBrO2 + H2O, while positive electric potential attracts inhibitor species Br− to the proximity of it. We study oscillation amplitude and period stability in systems with ferroin loaded Dowex 50W-X2 and Dowex 50W-X8 beads during experiments exceeding 16 h. It has been observed, that the above mentioned resins generate a smaller number of CO2 bubbles close to the beads than Dowex 50W-X4, which makes Dowex 50W-X2 and Dowex 50W-X8 more suitable for applications in chemical computing. We report amplitude stability, oscillation frequency, merging and annihilation of travelling waves in a lattice of Dowex 50W-X8 beads (mesh size 50–100) in over 19 h long experiments with equimolar solution of malonic acid and 1,4-cyclohexanedione. This system looks as a promising candidate for chemical computing devices that can operate for a day.

Advanced Chemical Computing Using Discrete Turing Patterns in Arrays of Coupled Cells

We examine dynamical switching among discrete Turing patterns that enable chemical computing performed by mass-coupled reaction cells arranged as arrays with various topological configurations: three coupled cells in a cyclic array, four coupled cells in a linear array, four coupled cells in a cyclic array, and four coupled cells in a branched array. Each cell is operating as a continuous stirred tank reactor, within which the glycolytic reaction takes place, represented by a skeleton inhibitor-activator model where ADP plays the role of activator and ATP is the inhibitor. The mass coupling between cells is assumed to be operating in three possible transport regimes: (i) equal transport coefficients of the inhibitor and activator (ii) slightly faster transport of the activator, and (iii) strongly faster transport of the inhibitor. Each cellular array is characterized by two pairs of tunable parameters, the rate coefficients of the autocatalytic and inhibitory steps, and the transport coefficients of the coupling. Using stability and bifurcation analysis we identified conditions for occurrence of discrete Turing patterns associated with non-uniform stationary states. We found stable symmetric and/or asymmetric discrete Turing patterns coexisting with stable uniform periodic oscillations. To switch from one of the coexisting stable regimes to another we use carefully targeted perturbations, which allows us to build systems of logic gates specific to each topological type of the array, which in turn enables to perform advanced modes of chemical computing. By combining chemical computing techniques in the arrays with glycolytic excitable channels, we propose a cellular assemblage design for advanced chemical computing.

Applications of Information Theory Methods for Evolutionary Optimization of Chemical Computers

It is commonly believed that information processing in living organisms is based on chemical reactions. However, the human achievements in constructing chemical information processing devices demonstrate that it is difficult to design such devices using the bottom-up strategy. Here I discuss the alternative top-down design of a network of chemical oscillators that performs a selected computing task. As an example, I consider a simple network of interacting chemical oscillators that operates as a comparator of two real numbers. The information on which of the two numbers is larger is coded in the number of excitations observed on oscillators forming the network. The parameters of the network are optimized to perform this function with the maximum accuracy. I discuss how information theory methods can be applied to obtain the optimum computing structure.

Reaction fronts of the autocatalytic hydrogenase reaction

We have built a model to describe the hydrogenase catalyzed, autocatalytic, reversible hydrogen oxidation reaction where one of the enzyme forms is the autocatalyst. The model not only reproduces the experimentally observed front properties, but also explains the found hydrogen ion dependence. Furthermore, by linear stability analysis, two different front types are found in good agreement with the experiments.

Chemical memory with states coded in light controlled oscillations of interacting Belousov–Zhabotinsky droplets

The information storing potential of droplets, in which an oscillatory, photosensitive Belousov–Zhabotinsky (BZ) reaction proceeds, is investigated experimentally.