Organization of Nucleotides in Different Environments and the Formation of Pre-Polymers

(R)evolution of Viruses: Introduction to biothermodynamics of viruses

As of 26 April 2024, the International Committee on Taxonomy of Viruses has registered 14690 virus species. Of these, only several dozen have been chemically and thermodynamically characterized. Every virus species is characterized by a specific empirical formula and thermodynamic properties - enthalpy, entropy and Gibbs energy. These physical properties are used in a mechanistic model of virus-host interactions at the cell membrane and in the cytoplasm. This review article presents empirical formulas and Gibbs energies for all major variants of SARS-CoV-2. This article also reports and suggests a mechanistic model of evolutionary changes, with the example of time evolution of SARS-CoV-2 from 2019 to 2024.

Wet-dry cycles cause nucleic acid monomers to polymerize into long chains

The key first step in the oligomerization of monomers is to find an initiator, which is usually done by thermolysis or photolysis. We present a markedly different approach that initiates acid-catalyzed polymerization at the surface of water films or water droplets, which is the reactive phase during a wet-dry cycle in freshwater hot springs associated with subaerial volcanic landmasses. We apply this method to the oligomerization of different nucleic acids, a topic relevant to how it might be possible to go from simple nucleic acid monomers to long-chain polymers, a key step in forming the building blocks of life. It has long been known that dehydration at elevated temperatures can drive the synthesis of ester and peptide bonds, but this reaction has typically been carried out by incubating dry monomers at elevated temperatures. We report that single or multiple cycles of wetting and drying link mononucleotides by forming phosphodiester bonds. Mass spectrometric analysis reveals uridine monophosphate oligomers up to 53 nucleotides, with an abundance of 35 and 43 nt in length. Long-chain oligomers are also observed for thymidine monophosphate, adenosine monophosphate, and deoxyadenosine monophosphate after exposure to a few wet-dry cycles. Nanopore sequencing confirms that long linear chains are formed. Enzyme digestion shows that the linkage is the phosphodiester bond, which is further confirmed by 31P NMR and Fourier transform infrared spectroscopy. This suggests that nucleic acid oligomers were likely to be present on early Earth in a steady state of synthesis and hydrolysis.

Membrane Catalyzed Formation of Nucleotide Clusters and Their Role in the Origins of Life: Insights from Molecular Simulations and Lattice Modeling

One of the mysteries in studying the molecular "Origin of Life" is the emergence of RNA and RNA-based life forms, where nonenzymatic polymerization of nucleotides is a crucial hypothesis in formation of large RNA chains. The nonenzymatic polymerization can be mediated by various environmental settings, such as cycles of hydration and dehydration, temperature variations, and proximity to a variety of organizing matrices, such as clay, salt, fatty acids, lipid membrane, and mineral surface. In this work, we explore the influence of different phases of the lipid membrane toward nucleotide organization and polymerization in a simulated prebiotic setting. Our molecular simulations quantify the localization propensity of a mononucleotide, uridine monophosphate (UMP), in distinct membrane settings. We perform all-atom molecular dynamics (MD) simulations to estimate the role of the monophasic and biphasic membranes in modifying the behavior of UMPs localization and their clustering mechanism. Based on the interaction energy of mononucleotides with the membrane and their diffusion profile from our MD calculations, we developed a lattice-based model to explore the thermodynamic limits of the observations made from the MD simulations. The mathematical model substantiates our hypothesis that the lipid layers can act as unique substrates for "catalyzing" polymerization of mononucleotides due to the inherent spatiotemporal heterogeneity and phase change behavior.

Crystallization as a selection force at the polymerization of nucleotides in a prebiotic context

Accumulation and selection of nucleotides is one of the most challenging problems surrounding the origin of the first RNA molecules on our planet. In the current work we propose that guanosine 3',5' cyclic monophosphate could selectively crystallize upon evaporation of an acidic prebiotic pool containing various other nucleotides. The conditions of the evaporative crystallization are fully compatible with the subsequent acid catalyzed polymerization of this cyclic nucleotide reported in earlier studies and may be relevant in a broad range of possible prebiotic environments. Albeit cytidine 3',5' cyclic monophosphate has the ability to selectively accumulate under the same conditions, its crystal structure is not likely to support polymer formation.

Screening for Primordial RNA–Peptide Interactions Using High-Density Peptide Arrays

RNA–peptide interactions are an important factor in the origin of the modern mechanism of translation and the genetic code. Despite great progress in the bioinformatics of RNA–peptide interactions due to the rapid growth in the number of known RNA–protein complexes, there is no comprehensive experimental method to take into account the influence of individual amino acids on non-covalent RNA–peptide bonds. First, we designed the combinatorial libraries of primordial peptides according to the combinatorial fusion rules based on Watson–Crick mutations. Next, we used high-density peptide arrays to investigate the interaction of primordial peptides with their cognate homo-oligonucleotides. We calculated the interaction scores of individual peptide fragments and evaluated the influence of the peptide length and its composition on the strength of RNA binding. The analysis shows that the amino acids phenylalanine, tyrosine, and proline contribute significantly to the strong binding between peptides and homo-oligonucleotides, while the sum charge of the peptide does not have a significant effect. We discuss the physicochemical implications of the combinatorial fusion cascade, a hypothesis that follows from the amino acid partition used in the work.

The Formation of RNA Pre-Polymers in the Presence of Different Prebiotic Mineral Surfaces Studied by Molecular Dynamics Simulations

We used all-atom Molecular Dynamics (MD) computer simulations to study the formation of pre-polymers between the four nucleotides in RNA (AMP, UMP, CMP, GMP) in the presence of different substrates that could have been present in a prebiotic environment. Pre-polymers are C3'-C5' hydrogen-bonded nucleotides that have been suggested to be the precursors of phosphodiester-bonded RNA polymers. We simulated wet-dry cycles by successively removing water molecules from the simulations, from ~60 to 3 water molecules per nucleotide. The nine substrates in this study include three clay minerals, one mica, one phosphate mineral, one silica, and two metal oxides. The substrates differ in their surface charge and ability to form hydrogen bonds with the nucleotides. From the MD simulations, we quantify the interactions between different nucleotides, and between nucleotides and substrates. For comparison, we included graphite as an inert substrate, which is not charged and cannot form hydrogen bonds. We also simulated the dehydration of a nucleotide-only system, which mimics the drying of small droplets. The number of hydrogen bonds between nucleotides and nucleotides and substrates was found to increase significantly when water molecules were removed from the systems. The largest number of C3'-C5' hydrogen bonds between nucleotides occurred in the graphite and nucleotide-only systems. While the surface of the substrates led to an organization and periodic arrangement of the nucleotides, none of the substrates was found to be a catalyst for pre-polymer formation, neither at full hydration, nor when dehydrated. While confinement and dehydration seem to be the main drivers for hydrogen bond formation, substrate interactions reduced the interactions between nucleotides in all cases. Our findings suggest that small supersaturated water droplets that could have been produced by geysers or springs on the primitive Earth may play an important role in non-enzymatic RNA polymerization.

Visualizing RNA polymers produced by hot wet-dry cycling

It is possible that the transition from abiotic systems to life relied on RNA polymers that served as ribozyme-like catalysts and for storing genetic information. The source of such polymers is uncertain, but previous investigations reported that wet-dry cycles simulating prebiotic hot springs provide sufficient energy to drive condensation reactions of mononucleotides to form oligomers and polymers. The aim of the study reported here was to verify this claim and visualize the products prepared from solutions composed of single mononucleotides and 1:1 mixture of two mononucleotides. Therefore, we designed experiments that allowed comparisons of all such mixtures representing six combinations of the four mononucleotides of RNA. We observed irregular stringy patches and crystal strands when wet-dry cycling was performed at room temperature (20 °C). However, when the same solutions were exposed to wet-dry cycles at 80 °C, we observed what appeared to be true polymers. Their thickness was consistent with RNA-like products composed of covalently bonded monomers, while irregular strings and crystal segments of mononucleotides dried or cycled at room temperature were consistent with structures assembled and stabilized by weak hydrogen bonds. In a few instances we observed rings with short polymer attachments. These observations are consistent with previous claims of polymerization during wet-dry cycling. We conclude that RNA-like polymers and rings could have been synthesized non-enzymatically in freshwater hot springs on the prebiotic Earth with sizes sufficient to fold into ribozymes and genetic molecules required for life to begin.

Tracing the Primordial Chemical Life of Glycine: A Review from Quantum Chemical Simulations

Glycine (Gly), NH2CH2COOH, is the simplest amino acid. Although it has not been directly detected in the interstellar gas-phase medium, it has been identified in comets and meteorites, and its synthesis in these environments has been simulated in terrestrial laboratory experiments. Likewise, condensation of Gly to form peptides in scenarios resembling those present in a primordial Earth has been demonstrated experimentally. Thus, Gly is a paradigmatic system for biomolecular building blocks to investigate how they can be synthesized in astrophysical environments, transported and delivered by fragments of asteroids (meteorites, once they land on Earth) and comets (interplanetary dust particles that land on Earth) to the primitive Earth, and there react to form biopolymers as a step towards the emergence of life. Quantum chemical investigations addressing these Gly-related events have been performed, providing fundamental atomic-scale information and quantitative energetic data. However, they are spread in the literature and difficult to harmonize in a consistent way due to different computational chemistry methodologies and model systems. This review aims to collect the work done so far to characterize, at a quantum mechanical level, the chemical life of Gly, i.e., from its synthesis in the interstellar medium up to its polymerization on Earth.

Ensemble Switching Unveils a Kinetic Rheostat Mechanism of the Eukaryotic Thiamine Pyrophosphate Riboswitch

Thiamine pyrophosphate (TPP) riboswitches regulate thiamine metabolism by inhibiting the translation of enzymes essential to thiamine synthesis pathways upon binding to thiamine pyrophosphate in cells across all domains of life. Recent work on the Arabidopsis thaliana TPP riboswitch suggests a multi-step TPP binding process involving multiple riboswitch configurational ensembles and that Mg2+ dependence underlies the mechanism of TPP recognition and subsequent transition to the expression-inhibiting state of the aptamer domain followed by changes in the expression platform. However, details of the relationship between TPP riboswitch conformational changes and interactions with TPP and Mg2+ ¬¬in the aptamer domain constituting this mechanism are unknown. Therefore, we integrated single-molecule multiparameter fluorescence and force spectroscopy with atomistic molecular dynamics simulations and found that conformational transitions within the aptamer domain's sensor helices associated with TPP and Mg2+ ligand binding occurred between at least five different ensembles on timescales ranging from µs to ms. These dynamics are orders of magnitude faster than the 10 second-timescale folding kinetics associated with expression-state switching in the switch sequence. Together, our results show that a TPP and Mg2+ dependent mechanism determines dynamic configurational state ensemble switching of the aptamer domain's sensor helices that regulates the stability of the switch helix, which ultimately may lead to the expression-inhibiting state of the riboswitch. Additionally, we propose that two pathways exist for ligand recognition and that this mechanism underlies a kinetic rheostat-like behavior of the Arabidopsis thaliana TPP riboswitch.

AFM Images of Viroid-Sized Rings That Self-Assemble from Mononucleotides through Wet-Dry Cycling: Implications for the Origin of Life

It is possible that early life relied on RNA polymers that served as ribozyme-like catalysts and for storing genetic information. The source of such polymers is uncertain, but previous investigations reported that wet–dry cycles simulating prebiotic hot springs provide sufficient energy to drive condensation reactions of mononucleotides to form oligomers. The aim of the study reported here was to visualize the products by atomic force microscopy. In addition to globular oligomers, ring-like structures ranging from 10–200 nm in diameter, with an average around 30–40 nm, were abundant, particularly when nucleotides capable of base pairing were present. The thickness of the rings was consistent with single stranded products, but some had thicknesses indicating base pair stacking. Others had more complex structures in the form of short polymer attachments and pairing of rings. These observations suggest the possibility that base-pairing may promote polymerization during wet–dry cycling followed by solvation of the rings. We conclude that RNA-like rings and structures could have been synthesized non-enzymatically on the prebiotic Earth, with sizes sufficient to fold into ribozymes and genetic molecules required for life to begin.

Emergence of ribozyme and tRNA-like structures from mineral-rich muddy pools on prebiotic earth

The RNA world hypothesis, although a viable one regarding the origin of life on earth, has so far failed to provide a compelling explanation for the synthesis of RNA enzymes from free nucleotides via abiotic processes. To tackle this long-standing problem, we develop a realistic model for the onset of the RNA world, using experimentally determined rates for polymerization reactions. We start with minimal assumptions about the initial state that only requires the presence of short oligomers or just free nucleotides and consider the effects of environmental cycling by dividing a day into a dry, semi-wet and wet phases that are distinguished by the nature of reactions they support. Long polymers, with maximum lengths sometimes exceeding 100 nucleotides, spontaneously emerge due to a combination of non-enzymatic, non-templated polymer extension and template-directed primer extension processes. The former helps in increasing the lengths of RNA strands, whereas the later helps in producing complementary copies of the strands. Strands also undergo hydrolysis in a structure-dependent manner that favour breaking of bonds connecting unpaired nucleotides. We identify the most favourable conditions needed for the emergence of ribozyme and tRNA-like structures and double stranded RNA molecules, classify all RNA strands on the basis of their secondary structures and determine their abundance in the population. Our results indicate that under suitable environmental conditions, non-enzymatic processes would have been sufficient to lead to the emergence of a variety of ribozyme-like molecules with complex secondary structures and potential catalytic functions.

NMR analysis of nucleotide π-stacking in prebiotically relevant crowded environment

The prebiotic soup of a putative 'RNA World' would have been replete with a plethora of molecules resulting from complex chemical syntheses and exogeneous delivery. The presence of background molecules could lead to molecular crowding, potentially affecting the course of the reactions facilitated therein. Using NMR spectroscopy, we have analyzed the effect of crowding on the stacking ability of RNA monomers. Our findings corroborate that the purines stack more efficiently than the pyrimidine ribonucleotides. This competence is further enhanced in the presence of a crowding agent. This enhanced stacking could result in greater sequestration of the purine monomers, putting their ready availability for relevant nonenzymatic reactions into question. Thus, this study demonstrates the need for systematic characterization of molecular crowding in the context of prebiotically pertinent processes. Unraveling such phenomena is essential for our understanding of the transition from abiotic to biotic, during the origin of life.

The Hot Spring Hypothesis for an Origin of Life

We present a testable hypothesis related to an origin of life on land in which fluctuating volcanic hot spring pools play a central role. The hypothesis is based on experimental evidence that lipid-encapsulated polymers can be synthesized by cycles of hydration and dehydration to form protocells. Drawing on metaphors from the bootstrapping of a simple computer operating system, we show how protocells cycling through wet, dry, and moist phases will subject polymers to combinatorial selection and draw structural and catalytic functions out of initially random sequences, including structural stabilization, pore formation, and primitive metabolic activity. We propose that protocells aggregating into a hydrogel in the intermediate moist phase of wet-dry cycles represent a primitive progenote system. Progenote populations can undergo selection and distribution, construct niches in new environments, and enable a sharing network effect that can collectively evolve them into the first microbial communities. Laboratory and field experiments testing the first steps of the scenario are summarized. The scenario is then placed in a geological setting on the early Earth to suggest a plausible pathway from life's origin in chemically optimal freshwater hot spring pools to the emergence of microbial communities tolerant to more extreme conditions in dilute lakes and salty conditions in marine environments. A continuity is observed for biogenesis beginning with simple protocell aggregates, through the transitional form of the progenote, to robust microbial mats that leave the fossil imprints of stromatolites so representative in the rock record. A roadmap to future testing of the hypothesis is presented. We compare the oceanic vent with land-based pool scenarios for an origin of life and explore their implications for subsequent evolution to multicellular life such as plants. We conclude by utilizing the hypothesis to posit where life might also have emerged in habitats such as Mars or Saturn's icy moon Enceladus. "To postulate one fortuitously catalyzed reaction, perhaps catalyzed by a metal ion, might be reasonable, but to postulate a suite of them is to appeal to magic." -Leslie Orgel.

Interaction of the mononucleotide UMP with a fluid phospholipid bilayer

Interaction between mononucleotides and lipid membranes is believed to have played an important role in the origin of life on Earth. Studies on mononucleotide-lipid systems hitherto have focused on the influence of the lipid environment on the organization of the mononucleotide molecules, and the effect of the latter on the confining medium has not been investigated in detail. We have probed the interaction of the mononucleotide, uridine 5'-monophosphate (UMP), and its disodium salt (UMPDSS) with fluid dimyristoylphosphatidylcholine (DMPC) membranes, using small-angle X-ray scattering (SAXS), cryogenic scanning electron microscopy (cryo-SEM) and computer simulations. UMP adsorbs and charges the lipid membrane, resulting in the formation of unilamellar vesicles in dilute solutions. Adsorption of UMP reduces the bilayer thickness of DMPC. UMPDSS has a much weaker effect on interbilayer interactions. These observations are in very good agreement with the results of an all-atom molecular dynamics simulation of these systems. In the presence of counterions, such as Na⁺, UMP forms small aggregates in water, which bind to the bilayer without significantly perturbing it. The phosphate moiety in the lipid headgroup is found to bind to the hydrogens from the sugar ring of UMP, while the choline group tends to bind to the two oxygens from the nucleotide base. These studies provide important insights into lipid-nucleotide interactions and the effect of the nucleotide on lipid membranes.

Formation of Abasic Oligomers in Nonenzymatic Polymerization of Canonical Nucleotides

Polymerization of nucleotides under prebiotically plausible conditions has been a focus of several origins of life studies. Non-activated nucleotides have been shown to undergo polymerization under geothermal conditions when subjected to dry-wet cycles. They do so by a mechanism similar to acid-catalyzed ester-bond formation. However, one study showed that the low pH of these reactions resulted in predominantly depurination, thereby resulting in the formation of abasic sites in the oligomers. In this study, we aimed to systematically characterize the nature of the oligomers that resulted in reactions that involved one or more of the canonical ribonucleotides. All the reactions analyzed showed the presence of abasic oligomers, with purine nucleotides being affected the most due to deglycosylation. Even in the reactions that contained nucleotide mixtures, the presence of abasic oligomers was detected, which suggested that information transfer would be severely hampered due to losing the capacity to base pair via H-bonds. Importantly, the stability of the N-glycosidic linkage, under conditions used for dry-wet cycling, was also determined. Results from this study further strengthen the hypothesis that chemical evolution in a pre-RNA World would have been vital for the evolution of informational molecules of an RNA World. This is evident in the high degree of instability displayed by N-glycosidic bonds of canonical purine ribonucleotides under the same geothermal conditions that otherwise readily favors polymerization. Significantly, the resultant product characterization in the reactions concerned underscores the difficulty associated with analyzing complex prebiotically relevant reactions due to inherent limitation of current analytical methods.

Organization and Compartmentalization by Lipid Membranes Promote Reactions Related to the Origin of Cellular Life

Liquid crystals have certain physical properties that promote chemical reactions which cannot occur in bulk phase media. These properties are displayed, among other molecules, by amphiphilic compounds which assemble into membrane structures then concentrate and organize biologically relevant monomers within their confined spaces. When mixtures of lipids and nucleotides are cycled multiple times between hydrated and anhydrous conditions, the monomers polymerize in the dry phase into oligonucleotides. Upon rehydration, mixtures of the polymers are encapsulated in lipid-bounded compartments called protocells. Reactions in liquid crystalline organizing matrices represent a promising approach for future research on how primitive cells could emerge on the early Earth and other habitable planets.

How do Nucleotides Adsorb Onto Clays?

Adsorption of prebiotic building blocks is proposed to have played a role in the emergence of life on Earth. The experimental and theoretical study of this phenomenon should be guided by our knowledge of the geochemistry of the habitable early Earth environments, which could have spanned a large range of settings. Adsorption being an interfacial phenomenon, experiments can be built around the minerals that probably exhibited the largest specific surface areas and were the most abundant, i.e., phyllosilicates. Our current work aims at understanding how nucleotides, the building blocks of RNA and DNA, might have interacted with phyllosilicates under various physico-chemical conditions. We carried out and refined batch adsorption studies to explore parameters such as temperature, pH, salinity, etc. We built a comprehensive, generalized model of the adsorption mechanisms of nucleotides onto phyllosilicate particles, mainly governed by phosphate reactivity. More recently, we used surface chemistry and geochemistry techniques, such as vibrational spectroscopy, low pressure gas adsorption, X-ray microscopy, and theoretical simulations, in order to acquire direct data on the adsorption configurations and localization of nucleotides on mineral surfaces. Although some of these techniques proved to be challenging, questioning our ability to easily detect biosignatures, they confirmed and complemented our pre-established model.

Templated replication (or lack thereof) under prebiotically pertinent conditions

Accurate replication of encoded information would have been crucial for the formation and propagation of functional ribozymes during the early evolution of life. Studies aimed at understanding prebiotically pertinent nonenzymatic reactions have predominantly used activated nucleotides. However, the existence of concentrated pools of activated monomers on prebiotic Earth is debatable. In this study, we explored the feasibility of nonenzymatic copying reactions using the more prebiotically relevant 5′-nucleoside monophosphates (5′-NMP). These reactions, involving a 20-mer primer, were performed in the presence of amphiphiles, under volcanic geothermal conditions. Interestingly, the extended primer was not comparable to the expected full length 21-mer product. Our results suggest loss of the nitrogenous base in the extended primer. This phenomenon persisted even after lowering the temperature and when different rehydration solutions were used. We envisage that the loss of the informational moiety on the incoming 5′-NMP, might be occurring during addition of this monomer to the pre-existing oligomer. Significantly, when 5′-ribose monophosphate was used, multiple additions to the aforementioned primer were observed that resulted in hybrid polymers. Such hybrid oligomers could have been important for exploring a vast chemical space of plausible alternate nucleobases, thus having important implications for the origin of primitive informational polymers.

Mobility of a Mononucleotide within a Lipid Matrix: A Neutron Scattering Study

An essential question in studies on the origins of life is how nucleic acids were first synthesized and then incorporated into compartments about 4 billion years ago. A recent discovery is that guided polymerization within organizing matrices could promote a non-enzymatic condensation reaction allowing the formation of RNA-like polymers, followed by encapsulation in lipid membranes. Here, we used neutron scattering and deuterium labelling to investigate 5′-adenosine monophosphate (AMP) molecules captured in a multilamellar phospholipid matrix. The aim of the research was to determine and compare how mononucleotides are captured and differently organized within matrices and multilamellar phospholipid structures and to explore the role of water in organizing the system to determine at which level the system becomes sufficiently anhydrous to lock the AMP molecules into an organized structure and initiate ester bond synthesis. Elastic incoherent neutron scattering experiments were thus employed to investigate the changes of the dynamic properties of AMP induced by embedding the molecules within the lipid matrix. The influence of AMP addition to the lipid membrane organization was determined through diffraction measurement, which also helped us to define the best working Q range for dynamical data analysis with respect to specific hydration. The use of different complementary instruments allowed coverage of a wide time-scale domain, from ns to ps, of atomic mean square fluctuations, providing evidence of a well-defined dependence of the AMP dynamics on the hydration level.

Co-operation between Polymerases and Nucleotide Synthetases in the RNA World

It is believed that life passed through an RNA World stage in which replication was sustained by catalytic RNAs (ribozymes). The two most obvious types of ribozymes are a polymerase, which uses a neighbouring strand as a template to make a complementary sequence to the template, and a nucleotide synthetase, which synthesizes monomers for use by the polymerase. When a chemical source of monomers is available, the polymerase can survive on its own. When the chemical supply of monomers is too low, nucleotide production by the synthetase is essential and the two ribozymes can only survive when they are together. Here we consider a computational model to investigate conditions under which coexistence and cooperation of these two types of ribozymes is possible. The model considers six types of strands: the two functional sequences, the complementary strands to these sequences (which are required as templates), and non-functional mutants of the two sequences (which act as parasites). Strands are distributed on a two-dimensional lattice. Polymerases replicate strands on neighbouring sites and synthetases produce monomers that diffuse in the local neighbourhood. We show that coexistence of unlinked polymerases and synthetases is possible in this spatial model under conditions in which neither sequence could survive alone; hence, there is a selective force for increasing complexity. Coexistence is dependent on the relative lengths of the two functional strands, the strand diffusion rate, the monomer diffusion rate, and the rate of deleterious mutations. The sensitivity of this two-ribozyme system suggests that evolution of a system of many types of ribozymes would be difficult in a purely spatial model with unlinked genes. We therefore speculate that linkage of genes onto mini-chromosomes and encapsulation of strands in protocells would have been important fairly early in the history of life as a means of enabling more complex systems to evolve.

Trans-acting polymerases are cooperative, because they copy neighbouring strands, and do not copy themselves directly. Inaccurate replication creates parasitic strands that act as templates but not ribozymes. It is known that in spatially distributed models with slow strand diffusion, clusters of cooperating polymerases arise that can survive in the presence of parasites provided that the error rate is less than a maximum limit (the error threshold). In the RNA World, we envisage multiple types of ribozymes working together. We would like to understand how a multi-ribozyme system could evolve from a system with a single type of polymerase ribozyme. As a first step in increasing complexity, we consider a two-ribozyme system in which there is one polymerase and one nucleotide synthetase that produces monomers for use by the polymerase. We are particularly interested to find conditions in which the chemical supply of monomers is too low for the polymerase to survive alone, but the additional monomers created by the synthetase allow the two-ribozyme system to survive where the single-ribozyme system could not. There is then a selective force for increasing the complexity of the system. Here we show that spatial clustering is sufficient to allow cooperation and survival of systems of unlinked ribozymes with different functions. Clusters form in which synthetases form fringes around the polymerases. Survival of the two-ribozyme system depends on several factors. The strand diffusion rate must be slow enough for cooperative clusters to emerge. The replication rate of the polymerase must be comparable to that of the synthetase. The diffusion rate of the monomers must be neither too slow nor too fast. The model considers the most difficult case for cooperation–unlinked genes with no compartments. The sensitivity of the two-ribozyme system that we study here suggests that evolution of a spatial system with multiple unlinked ribozymes would become increasingly more difficult as the number of components increased, and suggests that linkage and protocells would need to evolve relatively early in the history of life.

The Effect of Limited Diffusion and Wet-Dry Cycling on Reversible Polymerization Reactions: Implications for Prebiotic Synthesis of Nucleic Acids

A long-standing problem for the origins of life is that polymerization of many biopolymers, including nucleic acids and peptides, is thermodynamically unfavourable in aqueous solution. If bond making and breaking is reversible, monomers and very short oligomers predominate. Recent experiments have shown that wetting and drying cycles can overcome this problem and drive the formation of longer polymers. In the dry phase, bond formation is favourable, but diffusion is restricted, and bonds only form between monomers that are initially close together. In the wet phase, some of the bonds are hydrolyzed. However, repositioning of the molecules allows new bonds to form in the next dry phase, leading to an increase in mean polymer length. Here, we consider a simple theoretical model that explains the effect of cycling. There is an equilibrium length distribution with a high mean length that could be achieved if diffusion occurred freely in the dry phase. This equilibrium is inaccessible without diffusion. A single dry cycle without diffusion leads to mean lengths much shorter than this. Repeated cycling leads to a significant increase in polymerization relative to a single cycle. In the most favourable case, cycling leads to the same equilibrium length distribution as would be achieved if free diffusion were possible in the dry phase. These results support the RNA World scenario by explaining a potential route to synthesis of long RNAs; however, they also imply that cycling would be beneficial to the synthesis of other kinds of polymers, including peptides, where bond formation involves a condensation reaction.