Constrained stoichiometric network analysis

Biological homochirality and stoichiometric network analysis: Variations on Frank’s model

Biological homochirality is modelled using chemical reaction mechanisms that include autocatalytic and inhibition reactions as well as input and output flows. From the mathematical point of view, the differential equations associated with those mechanisms have to exhibit bistability.  The search for those bifurcations can be carried out using stoichiometric network analysis. This algorithm simplifies the mathematical analysis and can be implemented in a computer programme, which can help us to analyse chemical networks. However, regardless of the reduction to linear polynomials, which is made possible by this algorithm, in some cases, the complexity and length of the polynomials involved make the analysis unfeasible. This problem has been partially solved by extending the stoichiometric matrix with rows that code the duality relations between the different reactions occurring in the network given as input. All these facts allow us to analyse 28 different network models, highlighting the basic requirements needed by a chemical mechanism to have spontaneous mirror symmetry breaking.

Origins of oscillatory dynamics in the model of reactive oxygen species in the rhizosphere

Oscillatory processes are essential for normal functioning and survival of biological systems, and reactive oxygen species have a prominent role in many of them. A mechanism representing the dynamics of these species in the rhizosphere is analyzed using stoichiometric network analysis with the aim to determine its capabilities to simulate various dynamical states, including oscillations. A detailed analysis has shown that unstable steady states result from four destabilizing feedback cycles, among which the cycle involving hydroquinone, an electron acceptor, and its semi-reduced form is the dominant one responsible for the existence of saddle-node and Andronov-Hopf bifurcations. This requires a higher steady-state concentration for the reduced electron acceptor compared to that of the remaining species, where the level of oxygen steady-state concentration determines whether the Andronov-Hopf or saddle-node bifurcation will occur.

Automated Exploration of Prebiotic Chemical Reaction Space: Progress and Perspectives

Prebiotic chemistry often involves the study of complex systems of chemical reactions that form large networks with a large number of diverse species. Such complex systems may have given rise to emergent phenomena that ultimately led to the origin of life on Earth. The environmental conditions and processes involved in this emergence may not be fully recapitulable, making it difficult for experimentalists to study prebiotic systems in laboratory simulations. Computational chemistry offers efficient ways to study such chemical systems and identify the ones most likely to display complex properties associated with life. Here, we review tools and techniques for modelling prebiotic chemical reaction networks and outline possible ways to identify self-replicating features that are central to many origin-of-life models.