Necessary conditions for the emergence of homochirality via autocatalytic self-replication

Linear stability analysis of chemical mechanism, Listanalchem: A tool for the search of spontaneous mirror symmetry breaking

Homochirality, the phenomenon by which one of two virtually identical (non-superimposable mirror images) compounds is favored over the other in the chemistry of life, has been regarded as a requisite for the emergence of all living things on earth. Spontaneous mirror symmetry breaking has been proposed to produce the imbalance. Under this framework, Frank presented, in his foundational article "On spontaneous asymmetric synthesis", a simple chemical reaction network that displays spontaneous symmetry breaking for a specific set of reaction rates. Research has since focused on finding more complex and plausible models, each one with its advantages and disadvantages. Nevertheless, finding reaction rate values that make a model exhibit spontaneous symmetry breaking is a complex task, even for specially crafted models. LInear STability ANALysis of CHEmical Mechanism, Listanalchem, is a method and software for the search for appropriate reaction rates. It includes four different algorithms inspired by the analysis of Frank's network. Two classical algorithms are also included in Listanalchem: the Trace-Determinant plane and the Stoichiometric Network Analysis by Bruce Clarke. Listanalchem reads a chemical reaction network from plain text and runs one or more of the available algorithms according to the user selection. Listanalchem is tested and verified by studying classical, modified, and recent models proposed to explain the origin of biological homochirality.•Listanalchem allows a fast and reliable search for instability behavior in chemical mechanisms that pretend to explain spontaneous mirror symmetry breaking.•Listanalchem contains several model examples, including the most cited in the related literature.•Listanalchem is a tool that tests models that pretend to explain the origin of biological homochirality, helps find errors, and aids in designing new models.

The protometabolic nature of prebiotic chemistry

The field of prebiotic chemistry has been dedicated over decades to finding abiotic routes towards the molecular components of life. There is nowadays a handful of prebiotically plausible scenarios that enable the laboratory synthesis of most amino acids, fatty acids, simple sugars, nucleotides and core metabolites of extant living organisms. The major bottleneck then seems to be the self-organization of those building blocks into systems that can self-sustain. The purpose of this tutorial review is having a close look, guided by experimental research, into the main synthetic pathways of prebiotic chemistry, suggesting how they could be wired through common intermediates and catalytic cycles, as well as how recursively changing conditions could help them engage in self-organized and dissipative networks/assemblies (i.e., systems that consume chemical or physical energy from their environment to maintain their internal organization in a dynamic steady state out of equilibrium). In the article we also pay attention to the implications of this view for the emergence of homochirality. The revealed connectivity between those prebiotic routes should constitute the basis for a robust research program towards the bottom-up implementation of protometabolic systems, taken as a central part of the origins-of-life problem. In addition, this approach should foster further exploration of control mechanisms to tame the combinatorial explosion that typically occurs in mixtures of various reactive precursors, thus regulating the functional integration of their respective chemistries into self-sustaining protocellular assemblies.

Inducing Homochirality Through Intermediary Catalytic Species: A Stochastic Approach

A new chiral amplification mechanism based on a stochastic approach is proposed. The mechanism includes five different chemical species, an achiral substrate (A), two chiral forms (L, D), and two intermediary species (LA, DA). The process occurs within a small, semipermeable compartment that can be diffusively coupled with the outside environment. The study considers two alternative primary sources for chiral species within the compartment, one chemical and the other diffusive. As a remarkable fact, the chiral amplification process occurs due to stochastic fluctuations of an intermediary catalytic species (LA, DA) produced in situ, given the interaction of the chiral species with the achiral substrate. The net process includes two different steps: the synthesis of the catalyst (LA and DA) and the catalytic production of new chiral species from the substrate. Stochastic simulations show that proper parameterization can induce a robust chiral state within the compartment regardless of whether the system is open or closed. We also show how an increase in the non-catalytic production of chiral species tends to negatively impact the homochirality degree of the system. By its conception, the proposed mechanism suggests a deeper connection with how most biochemical processes occur in living beings, a fact that could open new avenues for studying this fascinating phenomenon.

Spontaneous emergence of enantioenriched chiral aldol reaction products from Achiral precursors in solution and origin of biological homochirality of sugars: a first-principles study

Experimental reports about observation of spontaneous mirror symmetry breaking and chiral amplification in stereoselective Mannich and aldol reactions, run under fully achiral initial conditions, have drawn a lot of attention, fuelled partly by the role these reactions could have played in chemical evolution as a cause for still puzzling observed homochirality of biomolecules, often considered a prerequisite for the origin of life. We have now revisited this still unresolved problem, using DFT computation of all combinatorially possible transition states and numerical solution of complete set of resulting coupled kinetic rate equations to model the aldol reaction rigorously "from the first principles" and without making any a priori assumptions. Spontaneous mirror symmetry breaking in this autocatalytic, reversible, closed and homogenous system is explained by a supercritical pitchfork bifurcation, occurring in concentrations of enantiomers due to time-delayed kinetic instability of racemic composition of reaction mixture, when reactants are initially provided in non-stoichiometric quantities. Same process, taking place under similar conditions in primordial "soup" of chemicals, might conceivably explain origin of biological homochirality of sugar molecules on early earth billions of years ago. Our results suggest that seemingly innocuous chemical reactions could exhibit unexpected and counter-intuitive emergent behaviour, when initial conditions are appropriately chosen. Chiral amplification in self-catalyzed aldol reaction occurs during approach of thermodynamic equilibrium in accord with principle of microscopic reversibility and second law of thermodynamics.

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.

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.

Equilibrium and non-equilibrium furanose selection in the ribose isomerisation network

The exclusive presence of β-D-ribofuranose in nucleic acids is still a conundrum in prebiotic chemistry, given that pyranose species are substantially more stable at equilibrium. However, a precise characterisation of the relative furanose/pyranose fraction at temperatures higher than about 50 °C is still lacking. Here, we employ a combination of NMR measurements and statistical mechanics modelling to predict a population inversion between furanose and pyranose at equilibrium at high temperatures. More importantly, we show that a steady temperature gradient may steer an open isomerisation network into a non-equilibrium steady state where furanose is boosted beyond the limits set by equilibrium thermodynamics. Moreover, we demonstrate that nonequilibrium selection of furanose is maximum at optimal dissipation, as gauged by the temperature gradient and energy barriers for isomerisation. The predicted optimum is compatible with temperature drops found in hydrothermal vents associated with extremely fresh lava flows on the seafloor.

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.

Magnetic circular dichroism in Archean atmosphere and asymmetric photolysis of biomolecules: enantiomeric excess of prebiotic sugar

In the terrestrial dipolar magnetic field, magnetic circular dichroism (MCD) of UV sunlight by paramagnetic O₂ in an Archean atmosphere (mostly CO₂ and N₂) results in circular polarization anisotropy (~ 10^-10). This is used to calculate enantiomeric excess (EE~10^-13) of glyceraldehyde (3-carbon sugar) with a model that includes racemic production and asymmetric photolysis of its enantiomers. The sign and magnitude of enantiomeric excess (EE) vary with the Earth's latitude. Unlike random noise fluctuation in spontaneous mirror symmetry breaking (SMSB) models, the sign of EE is deterministic and constant over large areas of prebiotic Earth. The magnitude is several orders greater than the mean amplitude of stochastically fluctuating EE. MCD could provide the initial EE for growth of homochirality by asymmetric autocatalysis.

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.

In search of a primitive signaling code

Cells must have preceded by simpler chemical systems (protocells) that had the capacity of a spontaneous self-assembly process and the ability to confine chemical reaction networks together with a form of information. The presence of lipid molecules in the early Earth conditions is sufficient to ensure the occurrence of spontaneous self-assembly processes, not defined by genetic information, but related to their chemical amphiphilic nature. Ribozymes are plausible molecules for early life, being the first small polynucleotides made up of random oligomers or formed by non-enzymatic template copying. Compartmentalization represents a strategy for the evolution of ribozymes; the attachment of ribozymes to surfaces, such as formed by lipid micellar aggregates may be particular relevant if the surface itself catalyzes RNA polymerization.It is conceivable that the transition from pre-biotic molecular aggregates to cellular life required the coevolution of the RNA world, capable of synthesizing specific, instead of statistical proteins, and of the Lipid world, with a transition from micellar aggregates to semipermeable vesicles. Small molecules available in the prebiotic inventory might promote RNA stability and the evolution of hydrophobic micellar aggregates into membrane-delimited vesicles. The transition from ribozymes catalyzing the assembly of statistical polypeptides to the synthesis of proteins, required the appearance of the genetic code; the transition from hydrophobic platforms favoring the stability of ribozymes and of nascent polypeptides to the selective transport of reagents through a membrane, required the appearance of the signal transduction code.A further integration between the RNA and Lipid worlds can be advanced, taking into account the emerging roles of phospholipid aggregates not only in ensuring stability to ribozymes by compartmentalization, but also in a crucial step of evolution through natural selection mechanisms, based on signal transduction pathways that convert environmental changes into biochemical responses that could vary according to the context. Here I present evidences on the presence of traces of the evolution of a signal transduction system in extant cells, which utilize a phosphoinositide signaling system located both at nucleoplasmic level as well as at the plasma membrane, based on the very same molecules but responding to different rules. The model herewith proposed is based on the following assumptions on the biomolecules of extant organisms: i) amphiphils can be converted into structured aggregates by hydrophobic forces thus giving rise to functional platforms for the interaction of other biomolecules and to their compartmentalization; ii) fundamental biochemical pathways, including protein synthesis, can be sustained by natural ribozymes of ancient origin; iii) ribozymes and nucleotide-derived coenzymes could have existed long before protein enzymes emerged; iv) signaling molecules, both derived from phospholipids and from RNAs could have guided the evolution of complex metabolic processes before the emergence of proteins.

Is pre-Darwinian evolution plausible?

Background

This essay highlights critical aspects of the plausibility of pre-Darwinian evolution. It is based on a critical review of some better-known open, far-from-equilibrium system-based scenarios supposed to explain processes that took place before Darwinian evolution had emerged and that resulted in the origin of the first systems capable of Darwinian evolution. The researchers’ responses to eight crucial questions are reviewed. The majority of the researchers claim that there would have been an evolutionary continuity between chemistry and “biology”. A key question is how did this evolution begin before Darwinian evolution had begun? In other words the question is whether pre-Darwinian evolution is plausible.

Results

Strengths and weaknesses of the reviewed scenarios are presented. They are distinguished between metabolism-first, replicator-first and combined metabolism-replicator models. The metabolism-first scenarios show major issues, the worst concerns heredity and chirality. Although the replicator-first scenarios answer the heredity question they have their own problems, notably chirality. Among the reviewed combined metabolism-replicator models, one shows the fewest issues. In particular, it seems to answer the chiral question, and eventually implies Darwinian evolution from the very beginning. Its main hypothesis needs to be validated with experimental data.

Conclusion

From this critical review it is that the concept of “pre-Darwinian evolution” appears questionable, in particular because it is unlikely if not impossible that any evolution in complexity over time may work without multiplication and heritability allowing the emergence of genetically and ecologically diverse lineages on which natural selection may operate. Only Darwinian evolution could have led to such an evolution. Thus, Pre-Darwinian evolution is not plausible according to the author. Surely, the answer to the question posed in the title is a prerequisite to the understanding of the origin of Darwinian evolution.

Reviewers

This article was reviewed by Purificacion Lopez-Garcia, Anthony Poole, Doron Lancet, and Thomas Dandekar.

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.

Toy trains, loaded dice and the origin of life: dimerization on mineral surfaces under periodic drive with Gaussian inputs

In a major extension of previous work, we model the putative hydrothermal rock pore setting for the origin of life on Earth as a series of coupled continuous flow units (the toy train). Perfusing through this train are reactants that set up thermochemical and pH oscillations, and an activated nucleotide that produces monomer and dimer monophosphates. The dynamical equations that model this system are thermally self-consistent. In an innovative step that breaks some new ground, we build stochasticity of the inputs into the model. The computational results infer various constraints and conditions on, and insights into, chemical evolution and the origin of life and its physical setting: long, interconnected porous structures with longitudinal non-uniformity would have been favourable, and the ubiquitous pH dependences of biology may have been established in the prebiotic era. We demonstrate the important role of Gaussian fluctuations of the inputs in driving polymerization, evolution and diversification. In particular, we find that the probability distribution of the resulting output fluctuations is left-skewed and right-weighted (the loaded dice), which could favour chemical evolution towards a living RNA world. We tentatively name this distribution 'Goldilocks'. These results also vindicate the general approach of constructing and running a simple model to learn important new information about a complex system.

Mechanically Induced Homochirality in Nucleated Enantioselective Polymerization

Understanding how biological homochirality may have emerged during chemical evolution remains a challenge for origin of life research. In keeping with this goal, we introduce and solve numerically a kinetic rate equation model of nucleated cooperative enantioselective polymerization in closed systems. The microreversible scheme includes (i) solution-phase racemization of the monomers, (ii) linear chain growth by stepwise monomer attachment, in both nucleation and elongation phases, and (iii) annealing or fusion of homochiral chains. Mechanically induced breakage of the longest chains maintains the system out of equilibrium and drives a breakage-fusion recycling mechanism. Spontaneous mirror symmetry breaking can be achieved starting from small initial enantiomeric excesses due to the intrinsic statistical fluctuations about the idealized racemic composition. The subsequent chiral amplification confirms the model's capacity for absolute asymmetric synthesis, without chiral cross-inhibition and without explicit autocatalysis.

Stoichiometric network analysis of spontaneous mirror symmetry breaking in chemical reactions

Stoichiometric network analysis (SNA) is used to study spontaneous mirror symmetry breaking in chemical reaction schemes.