Viedma deracemization mechanisms in self-assembly processes

Simulations on an ODE-based model shows that there are many common points between Viedma deracemization and chiral self-assemblies of achiral building blocks towards chiral nanoparticles. Both systems occur in a closed system with energy exchange but no matter exchange with the surroundings and show parallel reversible growth mechanisms which coexist with an irreversible cluster breaking (grinding). The various mechanisms of growth give rise to the formation of polymerization/depolymerization cycles while the consecutive transformation of achiral monomer into chiral cluster results into an indirect enantioselective autocatalysis. Deracemization occurs by the destabilization of the racemic non-equilibrium stationary state likely because of the excess of entropy production generated by the coupling of the reversible cluster growth mechanisms with grinding. Results show that the SMSB bias from the racemic composition occurs already at the oligomeric level of polymerization. Our model goes beyond the scope of the effect of grinding by the stirring of solutions which is thoroughly reported in supramolecular chirality. For instance, some unique characteristics, as those of a SMSB in closed systems, the simultaneous presence of different coupled reversible growth mechanisms, the activation by a depolymerization agent and the reincorporation of oligomers to the polymer growth reactions, could be adapted to replicator selectivity and to the emergence of biological homochirality scenarios.

Physical Chemistry Models for Chemical Research in the XXth and XXIst Centuries

Thermodynamic hypotheses and models are the touchstone for chemical results, but the actual models based on time-invariance, which have performed efficiently in the development of chemistry, are nowadays invalid for the interpretation of the behavior of complex systems exhibiting nonlinear kinetics and with matter and energy exchange flows with the surroundings. Such fields of research will necessarily foment and drive the use of thermodynamic models based on the description of irreversibility at the macroscopic level, instead of the current models which are strongly anchored in microreversibility.

Mathematical Models of Chiral Symmetry-breaking – A Review of General Theories, and Adiabatic Approximations of the APED System

We review the literature surrounding chiral symmetry-breaking in chemical systems, with a focus on understanding the mathematical models underlying these chemical processes. We comment in particular on the toy model of Sandars, Viedma’s crystal grinding systems and the APED model. We include a few new results based on asymptotic analysis of the APED system.

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.

Spontaneous Emergence of Transient Chirality in Closed, Reversible Frank-like Deterministic Models

To explore abiotic theories related to the origin of biomolecular homochirality, we analyze two entirely reversible kinetic models composed of an enantioselective autocatalysis with limited stereoselectivity that is coupled to an enantiomeric mutual inhibition (Frank-like models). The two models differ in their autocatalytic steps in respect to the formation of monomer species in one model and of dimer species in the other. While fully reversible and running in a closed system, spontaneous mirror symmetry breaking (SMSB) gives rise to transient chiral excursions, even when starting from a strictly achiral situation. Before the SMSB, the two models differ in the main dissipative processes. At the SMSB, the entropy production rate reaches its maximum in both models. Here it is the enantioselective autocatalysis with retention of the winner enantiomer that dominates. During the terminal phase, the enantioselective autocatalysis with inversion prevails, while the entropy production rate vanishes, thus fulfilling the conditions of microscopic reversibility. SMSB does not occur if the autocatalytic rate constant is too strong or too weak. However, when the autocatalysis is relatively weak, the temporary chiral excursions last for long periods of time and could be the starting point of a cascade of asymmetric reactions. The realism of such Frank-like models is discussed from the viewpoint of their relevance to prebiotic chemistry.

Stoichiometric network analysis in reaction networks yielding spontaneous mirror symmetry breaking in a prebiotic atmosphere

The generation of amino acid homochirality under prebiotic atmosphere conditions is a relevant issue in the study of the origin of life. This research is based on the production of amino acids via Strecker synthesis and how it is adjusted to the Kondepudi-Nelson autocatalytic model. The spontaneous mirror symmetry breaking (SMSB) of the new Kondepudi-Nelson-Strecker model, subject to two modifications (with Limited Enantioselective and Cross Inhibition), and also their combination were studied using the stoichiometric network analysis (SNA). In the calculations, the values obtained from the literature for alanine were considered. A total production of alanine of 7.56 × 10⁹ mol year^-1 was determined under prebiotic atmosphere conditions and starting from that value, the reaction rates for the models studied were estimated. Only the model with cross inhibition or achiral dimer formation is driven by stochastic fluctuations during SMSB. The stochastic fluctuation was estimated for a value of 2.619 × 10^-15 mol L^-1. This perturbation was sufficient to trigger SMSB. Finally, the results of SMSB were used to calculate the entropy production for the cross inhibition model.

A robust transition to homochirality in complex chemical reaction networks

Homochirality, i.e. the dominance across all living matter of one enantiomer over the other among chiral molecules, is thought to be a key step in the emergence of life. Building on ideas put forward by Frank and many others, we proposed recently one such mechanism in Laurent et al. (Laurent, 2021 Proc. Natl Acad. Sci. USA 118 , e2012741118. ( doi:10.1073/pnas.2012741118 )) based on the properties of large out of equilibrium chemical networks. We showed that in such networks, a phase transition towards a homochiral state is likely to occur as the number of chiral species in the system becomes large or as the amount of free energy injected into the system increases. This paper aims at clarifying some important points in that scenario, not covered by our previous work. We first analyse the various conventions used to measure chirality, introduce the notion of chiral symmetry of a network and study its implications regarding the relative chiral signs adopted by different groups of molecules. We then propose a generalization of Frank’s model for large chemical networks, which we characterize completely using methods of random matrices. This analysis is extended to sparse networks, which shows that the emergence of homochirality is a robust transition.

Symmetry of Post-Translational Modifications in a Human Enzyme

Paraoxonase 2 (PON2) is a member of a small family of human lactonases. Recently, post-translational modifications (PTMs) of PON2 were highlighted, one of which involved the modulation of the enzyme activity. Furthermore, two important single nucleotide polymorphisms (SNPs) involved in type 2 diabetes and its consequences, were found to modulate the enzyme activity as well. The position on the PON2 structural model of both residues corresponding to SNPs and PTMs suggested a symmetry of the molecule. By sequence and structure superposition we were able to confirm this finding. The result will be discussed in light of the evolution of symmetry in biological molecules and their function.

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.

The Coordinate Reaction Model: An Obstacle to Interpreting the Emergence of Chemical Complexity

The way chemical transformations are described by models based on microscopic reversibility does not take into account the irreversibility of natural processes, and therefore, in complex chemical networks working in open systems, misunderstandings may arise about the origin and causes of the stability of non-equilibrium stationary states, and general constraints on evolution in systems that are far from equilibrium. In order to be correctly simulated and understood, the chemical behavior of complex systems requires time-dependent models, otherwise the irreversibility of natural phenomena is overlooked. Micro reversible models based on the reaction-coordinate model are time invariant and are therefore unable to explain the evolution of open dissipative systems. The important points necessary for improving the modeling and simulations of complex chemical systems are: a) understanding the physical potential related to the entropy production rate, which is in general an inexact differential of a state function, and b) the interpretation and application of the so-called general evolution criterion (GEC), which is the general thermodynamic constraint for the evolution of dissipative chemical systems.

Effect of configuration-dependent multi-body forces on interconversion kinetics of a chiral tetramer model

We describe a reformulation of the four-site molecular model for chiral phenomena introduced by Latinwo et al. ["Molecular model for chirality phenomena," J. Chem. Phys. 145, 154503 (2016)]. The reformulation includes an additional eight-body force that arises from an explicit configuration-dependent term in the potential energy function, resulting in a coarse-grained energy-conserving force field for molecular dynamics simulations of chirality phenomena. In this model, the coarse-grained interaction energy between two tetramers depends on their respective chiralities and is controlled by a parameter λ, where λ < 0 favors local configurations involving tetramers of opposite chirality and λ > 0 gives energetic preference to configurations involving tetramers of the same chirality. We compute the autocorrelation function for a quantitative chirality metric and demonstrate that the multi-body force modifies the interconversion kinetics such that λ ≠ 0 increases the effective barrier for enantiomer inversion. Our simulations reveal that for λ > 0 and temperatures below a sharply defined threshold value, this effect is dramatic, giving rise to spontaneous chiral symmetry breaking and locking molecules into their chiral identity.

Symmetry breaking meets multisite modification

Multisite modification is a basic way of conferring functionality to proteins and a key component of post-translational modification networks. Additional interest in multisite modification stems from its capability of acting as complex information processors. In this paper, we connect two seemingly disparate themes: symmetry and multisite modification. We examine different classes of random modification networks of substrates involving separate or common enzymes. We demonstrate that under different instances of symmetry of the modification network (invoked explicitly or implicitly and discussed in the literature), the biochemistry of multisite modification can lead to the symmetry being broken. This is shown computationally and consolidated analytically, revealing parameter regions where this can (and in fact does) happen, and characteristics of the symmetry-broken state. We discuss the relevance of these results in situations where exact symmetry is not present. Overall, through our study we show how symmetry breaking (i) can confer new capabilities to protein networks, including concentration robustness of different combinations of species (in conjunction with multiple steady states); (ii) could have been the basis for ordering of multisite modification, which is widely observed in cells; (iii) can significantly impact information processing in multisite modification and in cell signalling networks/pathways where multisite modification is present; and (iv) can be a fruitful new angle for engineering in synthetic biology and chemistry. All in all, the emerging conceptual synthesis provides a new vantage point for the elucidation and the engineering of molecular systems at the junction of chemical and biological systems.

Proteins help our cells perform the chemical reactions necessary for life. Once proteins are made, they can also be modified in different ways. This can simply change their activity, or otherwise make them better suited for their specific jobs within the cell. Biological ‘catalysts’ called enzymes carry out protein modifications by reversibly adding (or removing) chemical groups, such as phosphate groups.

‘Multisite modifications’ occur when a protein has two or more modifications in different areas, which can be added randomly or in a specific sequence. The combination of all the modifications attached to a protein acts like a chemical barcode and confers a specific function to the protein.

Modification networks add levels of complexity above individual proteins. These encompass not only the proteins in a cell or tissue, but also the different enzymes that can modify them, and how they all interact with each other. Although our knowledge of these networks is substantial, basic aspects, such as how the ordering of multisite modification systems emerges, is still not well understood.

Using a simple set of multisite modifications, Ramesh and Krishnan set out to study the potential mechanisms allowing the creation of order in this context. Symmetry is a pervasive theme across the sciences. In biology, symmetry and how it may be broken, is important to understand, for example, how organism develop. Ramesh and Krishnan used the perspective of symmetry in protein networks to uncover the origins of ordering.

First, mathematical models of simple modification networks were created based on their basic descriptions. This system centred on proteins that could have phosphate modifications at two possible sites. The network was ‘symmetric’, meaning that the rate of different sets of chemical reactions was identical, as were the amounts of all the enzymes involved.

Dissecting the simulated network using a variety of mathematical approaches showed that its initial symmetry could break, giving rise to sets of ordered multisite modifications. Breaking symmetry did not require any additional features or factors; the basic chemical ‘ingredients’ of protein modification were all that was needed. The prism of symmetry also revealed other aspects of these multisite modification networks, such as robustness and oscillations.

This study sheds new light on the mechanism behind ordering of protein modifications. In the future, Ramesh and Krishnan hope that this approach can be applied to the study of not just proteins but also a wider range of biochemical networks.

Emergence of homochirality in large molecular systems

The selection of a single molecular handedness, or homochirality across all living matter, is a mystery in the origin of life. Frank's seminal model showed in the '50s how chiral symmetry breaking can occur in nonequilibrium chemical networks. However, an important shortcoming in this classic model is that it considers a small number of species, while there is no reason for the prebiotic system, in which homochirality first appeared, to have had such a simple composition. Furthermore, this model does not provide information on what could have been the size of the molecules involved in this homochiral prebiotic system. Here, we show that large molecular systems are likely to undergo a phase transition toward a homochiral state, as a consequence of the fact that they contain a large number of chiral species. Using chemoinformatics tools, we quantify how abundant chiral species are in the chemical universe of all possible molecules of a given length. Then, we propose that Frank's model should be extended to include a large number of species, in order to possess the transition toward homochirality, as confirmed by numerical simulations. Finally, using random matrix theory, we prove that large nonequilibrium reaction networks possess a generic and robust phase transition toward a homochiral state.

Entropic analysis of bistability and the general evolution criterion

An open flow bistable Schlögl model in a well-stirred isothermal reaction volume V showing the balance of the entropy production and the exchange entropy on the non-equilibrium stationary states.

Chirality: The Backbone of Chemistry as a Natural Science

Chemistry as a natural science occupies the length and temporal scales ranging between the formation of atoms and molecules as quasi-classical objects, and the formation of proto-life systems showing catalytic synthesis, replication, and the capacity for Darwinian evolution. The role of chiral dissymmetry in the chemical evolution toward life is manifested in how the increase of chemical complexity, from atoms and molecules to complex open systems, accompanies the emergence of biological homochirality toward life. Chemistry should express chirality not only as molecular structural dissymmetry that at the present is described in chemical curricula by quite effective pedagogical arguments, but also as a cosmological phenomenon. This relates to a necessarily better understanding of the boundaries of chemistry with physics and biology.

Spontaneous mirror symmetry breaking: an entropy production survey of the racemate instability and the emergence of stable scalemic stationary states

Stability of non-equilibrium stationary states and spontaneous mirror symmetry breaking, provoked by the destabilization of the racemic thermodynamic branch, is studied for enantioselective autocatalysis in an open flow system, and for a continuous range n of autocatalytic orders.

Computing the Parameter Values for the Emergence of Homochirality in Complex Networks

The goal of our research is the development of algorithmic tools for the analysis of chemical reaction networks proposed as models of biological homochirality. We focus on two algorithmic problems: detecting whether or not a chemical mechanism admits mirror symmetry-breaking; and, given one of those networks as input, sampling the set of racemic steady states that can produce mirror symmetry-breaking. Algorithmic solutions to those two problems will allow us to compute the parameter values for the emergence of homochirality. We found a mathematical criterion for the occurrence of mirror symmetry-breaking. This criterion allows us to compute semialgebraic definitions of the sets of racemic steady states that produce homochirality. Although those semialgebraic definitions can be processed algorithmically, the algorithmic analysis of them becomes unfeasible in most cases, given the nonlinear character of those definitions. We use Clarke’s system of convex coordinates to linearize, as much as possible, those semialgebraic definitions. As a result of this work, we get an efficient algorithm that solves both algorithmic problems for networks containing only one enantiomeric pair and a heuristic algorithm that can be used in the general case, with two or more enantiomeric pairs.

Chemical Basis of Biological Homochirality during the Abiotic Evolution Stages on Earth

Spontaneous mirror symmetry breaking (SMSB), a phenomenon leading to non-equilibrium stationary states (NESS) that exhibits biases away from the racemic composition is discussed here in the framework of dissipative reaction networks. Such networks may lead to a metastable racemic non-equilibrium stationary state that transforms into one of two degenerate but stable enantiomeric NESSs. In such a bifurcation scenario, the type of the reaction network, as well the boundary conditions, are similar to those characterizing the currently accepted stages of emergence of replicators and autocatalytic systems. Simple asymmetric inductions by physical chiral forces during previous stages of chemical evolution, for example in astrophysical scenarios, must involve unavoidable racemization processes during the time scales associated with the different stages of chemical evolution. However, residual enantiomeric excesses of such asymmetric inductions suffice to drive the SMSB stochastic distribution of chiral signs into a deterministic distribution. According to these features, we propose that a basic model of the chiral machinery of proto-life would emerge during the formation of proto-cell systems by the convergence of the former enantioselective scenarios.

Entropic Analysis of Mirror Symmetry Breaking in Chiral Hypercycles

Replicators are fundamental to the origin of life and evolvability. Biology exhibits homochirality: only one of two enantiomers is used in proteins and nucleic acids. Thermodynamic studies of chemical replicators able to lead to homochirality shed valuable light on the origin of homochirality and life in conformity with the underlying mechanisms and constraints. In line with this framework, enantioselective hypercyclic replicators may lead to spontaneous mirror symmetry breaking (SMSB) without the need for additional heterochiral inhibition reactions, which can be an obstacle for the emergence of evolutionary selection properties. We analyze the entropy production of a two-replicator system subject to homochiral cross-catalysis which can undergo SMSB in an open-flow reactor. The entropy exchange with the environment is provided by the input and output matter flows, and is essential for balancing the entropy production at the non-equilibrium stationary states. The partial entropy contributions, associated with the individual elementary flux modes, as defined by stoichiometric network analysis (SNA), describe how the system's internal currents evolve, maintaining the balance between entropy production and exchange, while minimizing the entropy production after the symmetry breaking transition. We validate the General Evolution Criterion, stating that the change in the chemical affinities proceeds in a way as to lower the value of the entropy production.

Stoichiometric network analysis of a reaction system with conservation constraints

Stoichiometric Network Analysis (SNA) is a powerful method that can be used to examine instabilities in modelling a broad range of reaction systems without knowing the explicit values of reaction rate constants. Due to a lack of understanding, SNA is rarely used and its full potential is not yet fulfilled. Using the oscillatory carbonylation of a polymeric substrate [poly(ethylene glycol)methyl ether acetylene] as a case study, in this work, we consider two mathematical methods for the application of SNA to the reaction models when conservation constraints between species have an important role. The first method takes conservation constraints into account and uses only independent intermediate species, while the second method applies to the full set of intermediate species, without the separation of independent and dependent variables. Both methods are used for examination of steady state stability by means of a characteristic polynomial and related Jacobian matrix. It was shown that both methods give the same results. Therefore, as the second method is simpler, we suggest it as a more straightforward method for the applications.

Stoichiometric network analysis of entropy production in chemical reactions

SNA extreme currents allow for the evaluation and understanding of entropy production of NESS in open system reaction networks.

Spontaneous mirror symmetry breaking and origin of biological homochirality

Recent reports on both theoretical simulations and on the physical chemistry basis of spontaneous mirror symmetry breaking (SMSB), that is, asymmetric synthesis in the absence of any chiral polarizations other than those arising from the chiral recognition between enantiomers, strongly suggest that the same nonlinear dynamics acting during the crucial stages of abiotic chemical evolution leading to the formation and selection of instructed polymers and replicators, would have led to the homochirality of instructed polymers. We review, in the first instance, which reaction networks lead to the nonlinear kinetics necessary for SMSB, and the thermodynamic features of the systems where this potentiality may be realized. This could aid not only in the understanding of SMSB, but also the design of reliable scenarios in abiotic evolution where biological homochirality could have taken place. Furthermore, when the emergence of biological chirality is assumed to occur during the stages of chemical evolution leading to the selection of polymeric species, one may hypothesize on a tandem track of the decrease of symmetry order towards biological homochirality, and the transition from the simple chemistry of astrophysical scenarios to the complexity of systems chemistry yielding Darwinian evolution.