Journal of the American Chemical Society
Home A Model for the Emergence of RNA from a Prebiotically Plausible Mixture of Ribonucleotides, Arabinonucleotides, and 2′-Deoxynucleotides
A Model
for the Emergence of RNA from a Prebiotically
Plausible Mixture of Ribonucleotides, Arabinonucleotides, and 2′-Deoxynucleotides
A Model for the Emergence of RNA from a Prebiotically Plausible Mixture of Ribonucleotides, Arabinonucleotides, and 2′-Deoxynucleotides

Article Type: research-article Article History
Abstract

The abiotic synthesis of ribonucleotides is thought to have been an essential step toward the emergence of the RNA world. However, it is likely that the prebiotic synthesis of ribonucleotides was accompanied by the simultaneous synthesis of arabinonucleotides, 2′-deoxyribonucleotides, and other variations on the canonical nucleotides. In order to understand how relatively homogeneous RNA could have emerged from such complex mixtures, we have examined the properties of arabinonucleotides and 2′-deoxyribonucleotides in nonenzymatic template-directed primer extension reactions. We show that nonenzymatic primer extension with activated arabinonucleotides is much less efficient than with activated ribonucleotides, and furthermore that once an arabinonucleotide is incorporated, continued primer extension is strongly inhibited. As previously shown, 2′-deoxyribonucleotides are also less efficiently incorporated in primer extension reactions, but the difference is more modest. Experiments with mixtures of nucleotides suggest that the coexistence of ribo- and arabinonucleotides does not impede the copying of RNA templates. Moreover, chimeric oligoribonucleotides containing 2′-deoxy- or arabinonucleotides are effective templates for RNA synthesis. We propose that the initial genetic polymers were random sequence chimeric oligonucleotides formed by untemplated polymerization, but that template copying chemistry favored RNA synthesis; multiple rounds of replication may have led to pools of oligomers composed mainly of RNA.

Kim,Zhou,Zhang,O’Flaherty,Rondo-Brovetto,and Szostak: A Model for the Emergence of RNA from a Prebiotically Plausible Mixture of Ribonucleotides, Arabinonucleotides, and 2′-Deoxynucleotides

Introduction

The RNA world hypothesis proposes that RNA oligonucleotides served dual roles as both genetic and functional polymers during the origin and early evolution of life.13 One of the essential requirements for the emergence of the RNA world is the prebiotic availability of ribonucleotides via abiotic synthesis.46 There has been significant progress in understanding the prebiotic synthesis of the pyrimidine ribonucleotides, led by Sutherland who has demonstrated prebiotically viable syntheses of pyrimidine ribonucleotides (U and C) and their 2-thio variants.7 Despite recent progress, equivalently direct and high yielding routes to the purine ribonucleotides remain unknown.8,9 A recent proposal for nucleotide synthesis from the Carrell group is unconstrained with respect to sugar chemistry,10 suggesting the possibility of a great diversity of nucleotides with differing sugar moieties. Thus, the continued exploration of plausible routes to the canonical ribonucleotides suggests that a diverse set of noncanonical nucleotides may have been present together with ribonucleotides on the early Earth.11,12 While UV photolysis has been elegantly shown to selectively destroy many undesired side-products in prebiotic pyrimidine nucleotide synthesis,7 the pool of available nucleotides still remains diverse,13 with nucleobase variations such as 8-oxo-purine,8 inosine,8,14 and the 2-thio-pyrimidines,15 as well as sugar variants including arabino-,9 2′-deoxyribo-,16 and threonucleotides.17,18

Intriguingly, some of these variant ribonucleotides have been shown to be compatible with or even beneficial for the copying of RNA templates. For instance, 2-thiouridine greatly improves the rate and fidelity of nonenzymatic primer extension compared with uridine,19 due to the 2-thio substituent increasing the binding affinity to adenosine and destabilizing wobble pairing with guanosine.2022 Inosine, which can arise from the deamination of adenosine, improves the fidelity of nonenzymatic primer extension relative to guanosine.14 In contrast, the 8-oxo-purine ribonucleotides are extremely poor substrates for nonenzymatic primer extension.14 Long before recent suggestions of prebiotic routes to the 2′-deoxyribonucleotides,16,23 these DNA building blocks were known to be less effective than ribonucleotides in nonenzymatic template copying reactions.24 However, recent work showing that mixed ribo/deoxy oligonucleotide duplexes have much lower melting temperatures than either homogeneous RNA or DNA duplexes25,26 suggests that this heterogeneity might facilitate replication by enhancing strand separation. Given the wide range of effects of noncanonical nucleotides on nonenzymatic template directed copying chemistry, a systematic study of the effects of prebiotically plausible variant nucleotides on copying chemistry is needed in order to understand if and how more homogeneous modern RNA could have arisen from complex prebiotic mixtures of nucleotides.

Among the plausible candidates for primordial nucleosides, arabinonucleosides, the building blocks of arabinonucleic acid (ANA), are of particular interest (Figure 1a). Powner et al. have recently demonstrated the first potentially prebiotic synthesis of all four canonical arabinonucleosides.27 In this approach, arabino- and ribonucleotides derive from the same synthetic intermediate, where external nucleophiles, such as hydroxide or sulfide, furnish the corresponding arabinonucleotides, while an intramolecular 3′ phosphate inverts the stereocenter at the 2′ position to furnish ribonucleotides (Figure 1b). The divergent syntheses of arabino- and ribonucleotides from a common intermediate suggests that the first genetic polymer may have been assembled from a heterogeneous mixture of ANA and RNA. Although homogeneous ANA alone cannot form stable duplexes,28 ANA can form a stable Watson−Crick paired duplex with the corresponding complementary RNA oligonucleotide.29,30 Recent studies have shown that ANA oligonucleotides can have RNA cleaving catalytic activity when paired across from an RNA oligonucleotide substrate.31 These results have been taken to suggest that ANA oligonucleotides can serve as both genetic and functional polymers in the presence of RNA oligonucleotides.

Divergent prebiotic syntheses of nucleosides.
(a) Schematic representation
of ribo-, arabino-, and 2′-deoxynucleosides. NB: nucleobase.
(b) Reported divergent syntheses of RNA (left), ANA (middle), and
DNA (right) from a common anhydronucleoside intermediate.
Figure 1

Divergent prebiotic syntheses of nucleosides. (a) Schematic representation of ribo-, arabino-, and 2′-deoxynucleosides. NB: nucleobase. (b) Reported divergent syntheses of RNA (left), ANA (middle), and DNA (right) from a common anhydronucleoside intermediate.

Herein, we evaluate the effects of arabinonucleotides on nonenzymatic primer extension, both as activated monomers and when incorporated into primers and templates. We have also revisited the effects of 2′-deoxyribonucleotides in parallel experiments, using more modern and prebiotically plausible activation chemistry than in earlier studies. In addition, we have examined the behavior of heterogeneous chimeric, predominantly ribonucleotide oligomers, by systematically introducing single variant residues to the RNA system.

Results and Discussion

Primer Extension with Monomers

To examine arabino- and 2′-deoxyribonucleotides as substrates in nonenzymatic primer extension, we measured extension rates using the 2-aminoimidazole activated monomers,32 2AI-araG, 2AI-araA, 2AI-dG, and 2AI-dA across from their respective Watson−Crick pairing partners C and s2U in RNA templates (Figure 2a). We used 2-thio-uridine (s2U) instead of uridine in the template because s2U is prebiotically plausible and exhibits stronger pairing with A than U. Previous reports from our lab have demonstrated that activated downstream oligonucleotides can greatly accelerate nonenzymatic primer extensions33 by forming highly reactive imidazolium-bridged intermediates34 and because of the structural preorganization afforded by an extended helical geometry.35 Therefore, we prepared 2-aminoimidazole activated RNA trimer helpers and utilized these activated downstream helpers in the evaluation of arabinonucleotides, 2′-deoxyribonucleotides, and ribonucleotides in nonenzymatic primer extension. The rates of primer extension using activated ribonucleotides were the fastest, at 11 h−1 for 2AI-rG and 16 h−1 for 2AI-rA. Primer extension rates with activated 2′-deoxynucleotides were slightly slower than with ribonucleotide monomers (8.2 and 13 h−1 for 2AI-dG and 2AI-dA respectively). In contrast, primer extension with arabinonucleotide monomers was much slower, at 0.88 h−1 for 2AI-araG and 0.52 h−1 for 2AI-araA (Figure 2b).

Evaluation of arabino- and 2′-deoxynucleotides in nonenzymatic
primer extension. (a) Schematic representation of a primer extension
reaction using an RNA template. 2AIpN represents 2-aminoimidazole
activated monomers and 2AIpNNN represents 2-aminoimidazole activated
RNA trimer helpers. (b) Gel electrophoresis images and rates of primer
extension for 2-aminoimidazole activated ribonucleotides, 2′-deoxynucleotides,
and arabinonucleotides. All reactions were performed at pH 8.0, 200
mM HEPES, 50 mM Mg2+, 20 mM 2AIpP1, 0.5 mM of
activated RNA trimer. See SI for kinetic
analysis of primer extension reactions (SI, Figures S8−S13). Values are reported as the mean ± SD in
brackets, with the last digit reported being the last significant
figure and the one in which the error arises from triplicate experiments.
Figure 2

Evaluation of arabino- and 2′-deoxynucleotides in nonenzymatic primer extension. (a) Schematic representation of a primer extension reaction using an RNA template. 2AIpN represents 2-aminoimidazole activated monomers and 2AIpNNN represents 2-aminoimidazole activated RNA trimer helpers. (b) Gel electrophoresis images and rates of primer extension for 2-aminoimidazole activated ribonucleotides, 2′-deoxynucleotides, and arabinonucleotides. All reactions were performed at pH 8.0, 200 mM HEPES, 50 mM Mg2+, 20 mM 2AIpP1, 0.5 mM of activated RNA trimer. See SI for kinetic analysis of primer extension reactions (SI, Figures S8−S13). Values are reported as the mean ± SD in brackets, with the last digit reported being the last significant figure and the one in which the error arises from triplicate experiments.

In order to distinguish between effects due to the formation, degradation, or steady state levels of the imidazolium-bridged intermediate, as opposed to the reactivity of that intermediate once formed, we prepared the 2-aminoimidazolium-bridged dinucleotides from ribo-, 2′-deoxy- and arabino-adenosine. We measured the rates of hydrolysis of the three imidazolium bridged dimers under primer extension conditions (50 mM Mg2+, pH 8 in 10% D2O; for details, see Supporting Information (SI), Figures S1−S4). Measured rates of hydrolysis revealed no significant difference in the stability of activated dimer species with k = 0.31, 0.33, and 0.30 h−1 for rA dimer, araA dimer, and dA dimer respectively, showing that the stereochemical information at the 2′ position does not impact the decomposition of the dimer intermediate under relevant conditions. In our evaluation of the rate of dimer intermediate formation, ribonucleotides, arabinonucleotides, and deoxynucleotides show no meaningful difference with k = 0.0032, 0.0031, and 0.0028 h−1 mM−1 for rA, araA, and dA respectively (see SI, Figures S5−S7).

We then evaluated the Michaelis−Menten parameters for primer extension for all three imidazolium-bridged dinucleotides. For these experiments, we prepared RNA templates containing two consecutive s2U residues by solid phase synthesis; we also used downstream oligonucleotides to further promote the nonenzymatic copying (Figure 3a). The Km values of the imidazolium-bridged dimer intermediates are very low, at 0.11, 0.065, and 0.64 mM for rA, dA, and araA respectively (Figure 3c). While these Km values span an order of magnitude, all values, including Km for araA, are low, in line with the ability to form ANA:RNA duplexes.25 Interestingly, Vmax for the three dimer intermediates differ by approximately 70 fold, at 22, 1.4, and 0.32 mM h−1 for rA, dA, and araA, respectively (Figure 3c). The large differences in Vmax are surprising, since the 2′-substituent is remote from the site of the reaction between the primer 3′-hydroxyl and the phosphate of the adjacent incoming nucleotide. Furthermore, simple nucleophiles such as hydroxide or water attack the phosphate of all three imidazolium-bridged dimers with similar rates, indicating that the 2′-substituent influences the rate of phosphate substitution only with stereochemically rich nucleophiles, such as the ribose at the end of the primer. Arabinonucleotides have been shown to adopt a C2′ endo pucker with intramolecular hydrogen bonding between the 2′-hydroxyl and 5′-oxygen,36 while 2′-deoxynucleotides show a more modest preference for a C2′ endo pucker and ribonucleotides adopt a C3′ endo pucker. This difference in sugar conformation could affect the position of the 5′-phosphate of the nucleotide adjacent to the 3′-end of the primer, potentially explaining the observed rate differences.

Michaelis−Menten study of nonenzymatic
copying using rA,
dA, and araA imidazolium bridged dimers. All reactions were performed
at pH 8.0, 200 mM HEPES, 50 mM Mg2+, and various concentrations
of homodimers ranging from 0.05 to 1 mM. (a) Schematic representation
of a primer extension reaction with imidazolium bridged dimers. (b)
Gel electrophoresis images of primer extension for imidazolium bridged
dimers at 1 mM concentration. (c) Plot of k (h−1) against the concentration of rA, dA, and araA imidazolium
bridged dimers. See SI for kinetic analysis
of primer extension reactions (SI, Figures S14−S28). Values are reported as the mean ± SD in brackets, with the
last digit reported being the last significant figure and the one
in which the error arises from triplicate experiments.
Figure 3

Michaelis−Menten study of nonenzymatic copying using rA, dA, and araA imidazolium bridged dimers. All reactions were performed at pH 8.0, 200 mM HEPES, 50 mM Mg2+, and various concentrations of homodimers ranging from 0.05 to 1 mM. (a) Schematic representation of a primer extension reaction with imidazolium bridged dimers. (b) Gel electrophoresis images of primer extension for imidazolium bridged dimers at 1 mM concentration. (c) Plot of k (h−1) against the concentration of rA, dA, and araA imidazolium bridged dimers. See SI for kinetic analysis of primer extension reactions (SI, Figures S14−S28). Values are reported as the mean ± SD in brackets, with the last digit reported being the last significant figure and the one in which the error arises from triplicate experiments.

Effect of Arabino- and 2′-Deoxynucleotides at the 3′-End of the Primer

Previous studies indicate that the nature of the 3′-end of the primer can strongly influence the rate of subsequent nonenzymatic primer extension;37 up to a 300 fold decrease in the rate of primer extension has been observed for a mismatched base pair vs a Watson−Crick base pair at the 3′-end of the primer.14 We therefore evaluated the effect on nonenzymatic primer extension of replacing the ribonucleotide at the 3′-end of the primer with either a 2′-deoxynucleotide or an arabinonucleotide. We prepared the three primers by solid phase synthesis, and then measured the rate of primer extension in the presence of 2-aminoimidazole activated guanosine monomers using templates containing the Watson−Crick paired primer binding region followed by a C3 templating region (Figure 4a). Strikingly, primers containing a single arabinonucleotide at the 3′-end showed no measurable extension even after 4 h. In contrast, primers with ribo- or 2′-deoxy-nucleotides at the 3′-end were rapidly extended under the same reaction conditions (kobs = 4.6 and 2.8 h−1 for primers containing a 3′ terminal rG and rA vs 2.1 and 1.3 h−1 for dG and dA, respectively), showing that the change in the configuration of the 2′-hydroxyl group on the pentose at the 3′ end of primer strongly affects nonenzymatic primer extension.

Evaluation
of nonenzymatic primer extension with primers containing
alternative nucleotides at the 3′ end. (a) Schematic representation
of nonenzymatic primer extension of 2AIrG with a primer that contains
either a ribo-, a 2′-deoxy-, or an arabinonucleotide at the
3′ end. (b) Gel electrophoresis images and rates of primer
extension for primers containing ribo-, 2′-deoxy-, and arabinonucleotides
at the 3′ end. See SI for kinetic
analysis of primer extension reactions (for further discussion see SI, Figures S29−S34). Values are reported
as the mean ± SD in brackets, with the last digit reported being
the last significant figure and the one in which the error arises
from triplicate experiments.
Figure 4

Evaluation of nonenzymatic primer extension with primers containing alternative nucleotides at the 3′ end. (a) Schematic representation of nonenzymatic primer extension of 2AIrG with a primer that contains either a ribo-, a 2′-deoxy-, or an arabinonucleotide at the 3′ end. (b) Gel electrophoresis images and rates of primer extension for primers containing ribo-, 2′-deoxy-, and arabinonucleotides at the 3′ end. See SI for kinetic analysis of primer extension reactions (for further discussion see SI, Figures S29−S34). Values are reported as the mean ± SD in brackets, with the last digit reported being the last significant figure and the one in which the error arises from triplicate experiments.

Given the dramatic stalling effect of an arabinonucleotide at the 3′ end of the primer, we sought to understand the reason for this effect at the atomic level. We therefore crystallized a chimeric oligonucleotide with a 3′-arabinonucleotide, complexed with RNA-monomers. The rGMP monomer was cocrystallized with a self-complementary 14mer oligonucleotide 5′-mCmCmCGACUUAAGUCaG-3′, which contains 4 locked nucleotides at the 5′ end (in bold) to rigidify the template and speed up crystallization (Figure 5A). The residue at the 3′ end is arabino-G (aG), to mimic the nonenzymatic primer extension shown in Figure 4. At each end, two GMP nucleotides bind to the templating locked 5-methyl Cs, forming a 16mer RNA-monomer helical structure. The assembled complex crystallizes in the rhombohedral R32 space group, and the structure was determined to 2.0-Å resolution. In one asymmetric unit, one strand of RNA is complexed with 2 GMP monomers. Except for the crystallographic symmetry, this arabino-modified RNA complex structure resembles our previously determined RNA-GMP complex structures,38,39 including the A-form double helices, C3′-endo conformation of the RNA sugars and slip-stacking of the neighboring duplexes (Figure 5B). The crystal structure reveals two key features of the terminal arabino-G residue: it is in the C2′-endo conformation, and forms a Watson−Crick base pair with the templating mC. The sugar of the terminal arabino-G is well-ordered, and its C2′-endo conformation contrasts with reported structures of an RNA/ANA heteroduplex, in which the arabinonucleotides are in the O4′-endo (Eastern) conformation25,29,40 (Figure 5C). The nucleobase of the arabino-G residue is Watson−Crick base paired with the templating mC, with hydrogen bonds ranging from 2.8 to 3.1 Å.

Crystal structure of
RNA primer/template complex with an arabinonucleotide
at the 3′-end of the primer (PDB: 6OWL). Red: primer; green: GMP monomers. The
wheat mesh indicates the 2Fo − Fc omit
map contoured at 1.5σ. (A) Schematic of the primer/template
complex. Italic letters: locked nucleotides; aG at the
3′ end: arabino-guanosine. (B) Overall structure of the 3′-arabino
RNA-GMP complex, stacking with neighboring complexes. (C) The local
structure of the arabino-G-terminated primer. The corresponding omit
map indicates that the primer arabinose is in the C2′-endo
conformation, and the guanine is Watson−Crick base paired with
the template. (D) The GMP monomers are bound through both canonical
and noncanonical G:C base pairs. The primer 2′−OH and
N7 of the +1 GMP monomer are within hydrogen bonding distance.
Figure 5

Crystal structure of RNA primer/template complex with an arabinonucleotide at the 3′-end of the primer (PDB: 6OWL). Red: primer; green: GMP monomers. The wheat mesh indicates the 2FoFc omit map contoured at 1.5σ. (A) Schematic of the primer/template complex. Italic letters: locked nucleotides; aG at the 3′ end: arabino-guanosine. (B) Overall structure of the 3′-arabino RNA-GMP complex, stacking with neighboring complexes. (C) The local structure of the arabino-G-terminated primer. The corresponding omit map indicates that the primer arabinose is in the C2′-endo conformation, and the guanine is Watson−Crick base paired with the template. (D) The GMP monomers are bound through both canonical and noncanonical G:C base pairs. The primer 2′−OH and N7 of the +1 GMP monomer are within hydrogen bonding distance.

In the evaluation of various nucleotides in nonenzymatic primer extension, several prior reports have shown a strong correlation between primer extension rates and the conformation of the 3′-terminal nucleotide. Primers composed of or ending in nucleotides that prefer a C3′-endo conformation (A-form) engage in faster primer extension than primers composed of or ending in nucleotides that do not adopt a C3′-endo conformation (e.g., DNA).24,41,42 Additionally, pseudoaxial/axial nucleophiles generally exhibit lower reactivities than pseudoequatorial/equatorial nucleophiles.43 This is due to steric encumbrance of axial/pseudoaxial substituents, especially by 1,3-diaxial interactions, which are exacerbated by a more crowded steric environment in the product and in the transition state. One plausible explanation for the significant stalling effect we observe in nonenzymatic polymerization with an arabinonucleotide terminated primer is that the preferred C2′-endo conformation of the arabinonucleotide places the 3′-hydroxyl too far away from the phosphate of the incoming monomer.

The base-pairing of the unactivated GMP monomer bound downstream of the primer is significantly different from that observed in our previously determined RNA-GMP structures. A hydrogen bonding interaction between the N7 atom of the +1 GMP and the 2′−OH group of the arabinonucleotide at the end of the primer is observed (2.9 Å, Figure 5D). This flips the nucleobase of the +1 GMP into a noncanonical G:C base pair in which N3 of the +1 GMP forms a hydrogen bond with the exocyclic amine of the templating mC, while the exocyclic amine of the guanine hydrogen bonds with N3 of the mC. This noncanonical base pairing positions the 5′-phosphate of the monomer distant from the 3′−OH of the primer. The phosphate of the +1 GMP is highly disordered, and the primer 3′−OH of the arabinose sugar is pointing to the minor groove instead of the incoming monomer. In previously determined RNA-GMP structures, the first bound GMP monomer was Watson−Crick paired with the template, with its binding affinity and the distance between the primer 3′−OH and the incoming phosphorus improved by the downstream bound monomer at the +2 position. The structural changes observed in the case of the arabino-terminated primer would be detrimental to nonenzymatic primer extension, in line with our experimental data; however, the crystal structure may not reflect the conformation of the template bound imidazolium-bridged dinucleotide.

Competition Experiment between Activated rA and araA

The common intermediate in a recently proposed prebiotic synthesis of RNA, DNA, and ANA raises the possibility that these nucleotides may have coexisted on early earth.27 Given the strong stalling effect when the 3′ end of the primer is an arabinonucleotide, the incorporation of an arabinonucleotide into a growing primer will essentially act as a chain terminator, potentially preventing the synthesis of oligonucleotides long enough to have catalytic function. Fortunately, arabinonucleotide incorporation rates were 15−30 fold slower than ribonucleotide incorporation rates in the case of trimer assisted nonenzymatic copying, suggesting that ribonucleotides may outcompete arabinonucleotides well enough to minimize this concern. To test this idea, we performed a nonenzymatic copying reaction in a one pot competition experiment with a 1:1 ratio of activated araA and rA monomers, utilizing the reaction conditions in Figure 2. The initial primer, P1, was exposed to activated rA and araA monomers with activated trimer helper (2AI-rGAC) for 1 h, to allow for incorporation of either rA (generating primer P2) or araA (generating primer P3). A subsequent nonenzymatic copying step was then performed by adding the activated RNA monomer 2AI-rG and the trimer helper 2AI-rACA. Primers that initially incorporated araA to generate P3 cannot be extended in the second step (Lane 7, Figure 6b) while primers that initially incorporated rA to generate primer P2 can be further extended (Lane 5, Figure 6b). After 4 h, the reaction was quenched and PAGE gel analysis was performed to analyze the product ratios (Figure 6b). Only a small fraction (approximately 7%) fails to extend after the initial addition (Lane 3, Figure 6b). Generally, iterative nonenzymatic copying of oligonucleotides is less than 100% complete,33,44 presumably due to the presence of multiple trimer helpers, which have high binding affinity to the template and are thus inhibitory to subsequent copying. Our results therefore represent a lower bound on the true ratio of rA vs araA incorporation, and indicate that activated rA monomers substantially outcompete activated araA in the nonenzymatic copying of oligonucleotides. Thus, the coexistence of arabino- and ribonucleotides should not significantly interfere with oligonucleotide copying.

Competition experiment
between 2AI-araA and 2AI-rA. (a) Schematic
representation of nonenzymatic primer extension with both ribonucleotide
and arabinonucleotide monomers followed by a nonenzymatic copying
step, in which we expect the ribo-terminated primer [P2] to extend
and the arabino-terminated primer [P3] not to extend (shown in lanes
5 and 7 in b). (b) Gel electrophoresis images of the competition experiment.
Primers used for the experiment are shown on the top ([P1] for lanes
2/3, [P2] for lanes 4/5, and [P3] for lanes 6/7) and the activated
nucleotide monomers and trimer helpers added are shown below. Synthetic
primers [P2] in lane 4 and [P3] in lane 6 were used as controls. Ctrl
lane 1 includes synthetic primers [P1] and [P2].
Figure 6

Competition experiment between 2AI-araA and 2AI-rA. (a) Schematic representation of nonenzymatic primer extension with both ribonucleotide and arabinonucleotide monomers followed by a nonenzymatic copying step, in which we expect the ribo-terminated primer [P2] to extend and the arabino-terminated primer [P3] not to extend (shown in lanes 5 and 7 in b). (b) Gel electrophoresis images of the competition experiment. Primers used for the experiment are shown on the top ([P1] for lanes 2/3, [P2] for lanes 4/5, and [P3] for lanes 6/7) and the activated nucleotide monomers and trimer helpers added are shown below. Synthetic primers [P2] in lane 4 and [P3] in lane 6 were used as controls. Ctrl lane 1 includes synthetic primers [P1] and [P2].

Primer Extension Across an Arabinonucleotide in the Template Strand

Although the above results suggest that arabinonucleotide-containing oligonucleotides would only rarely be generated by nonenzymatic primer extension, untemplated polymerization of activated monomers could lead to the synthesis of oligonucleotide templates containing arabinonucleotides in internal positions.45,46 We have evaluated the effects of internal arabinonucleotides in template strands on the rate of primer extension with standard activated ribonucleotides. We prepared RNA templates containing isolated arabino- or 2′-deoxynucleotides by solid phase synthesis, and then performed nonenzymatic primer extensions with activated RNA monomers and trimer helpers, and compared the observed rates on all-RNA templates and RNA templates containing single arabino- or 2′-deoxynucleotide residues (Figure 7a). Primer extension with 2-aminoimidazole activated cytidine exhibited a minimal rate difference when copying across rG, dG, and araG in the template strand, with kobs = 10, 8.0, and 6.8 h−1 (Figure 7b). Primer extension with activated 2-thiouridine also exhibited minor rate differences when copying across rA, dA and araA in the template strand, with kobs = 12, 11, and 2.1 h−1, respectively. In all cases, we observe complete conversion of the primer to extended products in the presence of the corresponding complementary Watson−Crick RNA monomers.

Evaluation of primer extension for 2-aminoimidazole
activated pyrimidine
nucleotides across from templating ribo-, 2′-deoxy-, and arabinonucleotides.
(a) Schematic representation of a primer extension reaction with template
containing a ribo-, 2′-deoxy-, or arabinonucleotide. 2AIpN
represents the corresponding 2-aminoimidazole activated ribonucleotide
monomers and 2AIpNNN represents 2-aminoimidazole activated RNA trimer
helper. (b) Gel electrophoresis images and rates of primer extension
for 2-aminoimidazole activated pyrimidine ribose across from a ribo-,
2′-deoxy- and arabinonucleotide (for further discussion see SI, Figures S35−S40). Values are reported
as the mean ± SD in brackets, with the last digit reported being
the last significant figure and the one in which the error arises
from triplicate experiments.
Figure 7

Evaluation of primer extension for 2-aminoimidazole activated pyrimidine nucleotides across from templating ribo-, 2′-deoxy-, and arabinonucleotides. (a) Schematic representation of a primer extension reaction with template containing a ribo-, 2′-deoxy-, or arabinonucleotide. 2AIpN represents the corresponding 2-aminoimidazole activated ribonucleotide monomers and 2AIpNNN represents 2-aminoimidazole activated RNA trimer helper. (b) Gel electrophoresis images and rates of primer extension for 2-aminoimidazole activated pyrimidine ribose across from a ribo-, 2′-deoxy- and arabinonucleotide (for further discussion see SI, Figures S35−S40). Values are reported as the mean ± SD in brackets, with the last digit reported being the last significant figure and the one in which the error arises from triplicate experiments.

Conclusions

The discovery of potentially prebiotic synthetic routes to the canonical ribonucleotides has given rise to an improved understanding of the most likely alternative nucleotide byproducts of these pathways. The likely presence of these byproducts leads to a significant problem with regard to the emergence of the RNA world, since the initially synthesized oligonucleotides would be expected to be quite heterogeneous in composition. How could such a heterogeneous mixture of oligonucleotides give rise to the relatively homogeneous RNAs that are thought to be required for the evolution of functional RNAs such as ribozymes?

In addressing the problem of the emergence of RNA, we have, in this paper, focused on the consequences of coexisting activated arabino- and 2′-deoxy-nucleotides for nonenzymatic template directed primer extension. We find that both primer extension with activated arabinonucleotides, and the extension of an arabinonucleotide-terminated primer with ribonucleotides, are very slow. Our one pot competition experiment between ribo- and arabinonucleotides directly demonstrates that highly RNA enriched oligonucleotides could arise by template copying with a heterogeneous mixture of nucleotides. Our kinetic and structural studies are consistent with this rate difference being due to the 2′-endo sugar conformation of arabinonucleotides. In parallel experiments, we find, as expected,24 that primer extension with 2′-deoxynucleotides, and extension of a 2′-deoxynucleotide terminated primer, are slower than with ribonucleotides, but the differences are much smaller than in the case of arabinonucleotides. In contrast, the copying of templates containing either an arabino- or a 2′-deoxynucleotide is only moderately slower than the copying of an all RNA template. Lastly, additional experiments with 2′-deoxynucleotide or arabinonucleotide terminated primers with activated deoxynucleotide and arabinonucleotide monomers show much slower, if any, extension across both homogeneous RNA templates or templates containing 2′-deoxynucleotide or arabinonucleotide overhangs (SI, Figures S41 and S42). Our results suggest that an initial set of highly heterogeneous oligonucleotides, generated by nontemplated polymerization, could act as templates for the synthesis of complementary strands that would be enriched in RNA. Multiple cycles of template copying by primer extension could therefore provide a pathway for a transition from heterogeneity to homogeneity. We have seen similar effects in previous work with 8-oxo-purine ribonucleotides, where primer extension with these activated nucleotides is extremely slow, but they can be copied when present in the template strand, albeit at lower rates.14

The Krishnamurthy group has recently proposed a similar pathway leading from an initially heterogeneous set of templates to an increasingly homogeneous set of replication products.26,47,48 In their case the physical mechanism leading to homogeneity is different, and is based on the surprising observation that homogeneous all-RNA or all-DNA oligonucleotides bind more strongly than do otherwise identical but heterogeneous oligonucleotides to complementary but heterogeneous templates. As a result of the stronger binding, the template-directed ligation of all RNA or all DNA oligonucleotides is favored over the ligation of mixed oligonucleotides, such that over repeated cycles of ligation-mediated replication, relatively homogeneous replication products should emerge. Because these experiments relied on nonprebiotic ligation chemistry and the nucleotide variation was restricted to internal positions, effects on ligation rates could not be evaluated. However, variations in chemical reactivity might also be expected to contribute to the emergence of homogeneous products. We have previously suggested that both primer extension with activated monomers and ligation mediated template copying may have acted together to mediate primitive replication cycles,49 in which case both ligation and primer extension effects may have contributed to the emergence of homogeneous RNA oligonucleotides as a result of cycles of replication. We note that both primer extension with monomers and oligonucleotide ligation are driven by the energy of nucleotide (or oligonucleotide) activation. The chemical energy of activated substrates drives template copying and ultimately replication in a far from equilibrium reaction network, which in turn powers the emergence of homogeneity from heterogeneity. Thus, local order in the form of homogeneous oligonucleotides emerges at the expense of chemical energy that is dissipated through an overall increase in disorder.

Our observations show that arabino-terminated products of primer extension cannot be further extended at significant rates. In the past, such dead-end products of incomplete replication would have been considered to be highly detrimental to the overall process of genomic replication in protocells. However, recent results suggest that such oligonucleotides could in fact have multiple positive or even essential roles in primordial genomic replication. For example, when suitably activated, short oligonucleotides act as catalysts of template copying by forming imidazolium-bridged intermediates with activated monomers.33,34 Slightly longer oligonucleotides may also act as catalysts of strand displacement synthesis through a toehold/branch migration mechanism that assists in opening up a template and allowing primer extension to proceed.50 Similar oligonucleotides could also act as splint templates for nonenzymatic ligation; since they could not be elongated by primer extension they could be much easier to dissociate from the longer ligated product, which could facilitate the construction of functional oligonucleotides.51 Short oligonucleotides may also contribute to homeostatic regulation of ribozyme activity in protocells.52

Multiple distinct processes, in addition to primer extension and ligation, are likely to have contributed to the transition from heterogeneous primordial nucleic acids to the relatively homogeneous RNA genomes of the first cells. For example, photochemical reactions preferentially degrade the noncanonical nucleobases (and the corresponding nucleosides and nucleotides) and also preferentially degrade the alpha-anomeric byproducts of nucleotide synthesis.7 Steric constraints and variations in chemical reactivity may have influenced the composition of the first oligonucleotides formed through nontemplated polymerization; for example, nucleotides with acyclic sugars rapidly cyclize to unreactive products following phosphate activation. The synthesis of standard 3′−5′ phosphodiester bonds when copying templates with 2′−5′ linkages53 or 3′−5′ pyrophosphate linkages54 may also have contributed to the gradual elimination of the variability in nucleic acid structure that is the inevitable consequence of nonenzymatic polymerization. It has also been shown that while mixed RNA/DNA oligonucleotides can exhibit functional properties such as molecular recognition, the homogeneous systems exhibit better function and, thus, the emergence of homogeneous genetic polymers could be the result of selective pressures for superior function.55 Taken together these and other mechanisms may explain the transition from a heterogeneous mixture of prebiotically synthesized nucleotides and oligonucleotides to a relatively homogeneous set of RNAs. Considerable additional experimental work must be done to extend this model, as only a fraction of the likely prebiotic variability in nucleotide and nucleic acid structure has been explored to date.

Supporting Information Available

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/jacs.9b11239.

    Methods and materials; oligonucleotide synthesis; synthesis of activated nucleotides; hydrolysis of activated araA and rA imidazolium bridged dimers; formation of activated araA and rA imidazolium bridged dimers; data collection statistics; PAGE analysis; kinetic analysis; sequences of oligonucleotides used in this study; additional references (PDF)

Author Present Address

§ Alnylam Pharmaceuticals, 300 Third Street, Cambridge, Massachusetts 02142, United States.

Notes

The authors declare no competing financial interest.

Acknowledgments

J.W.S. is an Investigator of the Howard Hughes Medical Institute. This work was supported in part by a grant (290363) from the Simons Foundation to J.W.S. and by a grant from the NSF (CHE-1607034) to J.W.S. D.K.O. is a recipient of a Postdoctoral Research Scholarship (B3) from the Fonds de recherche du Québec − Nature et technologies (FRQNT), Quebec, Canada, and a Postdoctoral Fellowship from Canadian Institutes of Health Research (CIHR) from Canada. The authors would like to thank Andrew J. Bendelsmith, Li Li, and Victor S. Lelyveld for their helpful discussions. The authors also thank the entire Szostak group for feedback.

References

Woese C. R.The Genetic Code: The Molecular Basis for Genetic Expression; Harper & Row: New York, 1967.
L. E.Orgel Evolution of the Genetic Apparatus. J. Mol. Biol. 1968, 38, 381393. 10.1016/0022-2836(68)90393-8.
W.Gilbert Origin of Life: The RNA World. Nature 1986, 319, 61810.1038/319618a0.
L. E.Orgel Prebiotic Chemistry and the Origin of the RNA World. Crit. Rev. Biochem. Mol. Biol. 2004, 39, 99123. 10.1080/10409230490460765.
G. F.Joyce The Antiquity of RNA-Based Evolution. Nature 2002, 418, 214221. 10.1038/418214a.
M. W.Powner; J. D.Sutherland; J. W.Szostak The Origins of Nucleotides. Synlett 2011, 14, 19561964. 10.1055/s-0030-1261177.
M. W.Powner; B.Gerland; J. D.Sutherland Synthesis of Activated Pyrimidine Ribonucleotides in Prebiotically Plausible Conditions. Nature 2009, 459, 239242. 10.1038/nature08013.
S.Stairs; A.Nikmal; D.-K.Bucar; S.-L.Zheng; J. W.Szostak; M. W.Powner Divergent Prebiotic Synthesis of Pyrimidine and 8-Oxo-Purine Ribonucleotides. Nat. Commun. 2017, 8, 1527010.1038/ncomms15270.
S. J.Roberts; R.Szabla; Z. R.Todd; S.Stairs; D. K.Bučar; J.Šponer; D. D.Sasselov; M. W.Powner Selective Prebiotic Conversion of Pyrimidine and Purine Anhydronucleosides into Watson-Crick Base-Pairing Arabino-Furanosyl Nucleosides in Water. Nat. Commun. 2018, 9, 407310.1038/s41467-018-06374-z.
S.Becker; J.Feldmann; S.Wiedemann; H.Okamura; C.Schneider; K.Iwan; A.Crisp; M.Rossa; T.Amatov; T.Carell Unified Prebiotically Plausible Synthesis of Pyrimidine and Purine RNA Ribonucleotides. Science 2019, 366, 7682. 10.1126/science.aax2747.
S.Becker; C.Schneider; A.Crisp; T.Carell Non-Canonical Nucleosides and Chemistry of the Emergence of Life. Nat. Commun. 2018, 9, 517410.1038/s41467-018-07222-w.
D.Whitaker; M. W.Powner Prebiotic Nucleic Acids Need Space to Grow. Nat. Commun. 2018, 9, 517210.1038/s41467-018-07221-x.
D. M.Fialho; K. C.Clarke; M. K.Moore; G. B.Schuster; R.Krishnamurthy; N. V.Hud Glycosylation of a Model Proto-RNA Nucleobase with Non-Ribose Sugars: Implications for the Prebiotic Synthesis of Nucleosides. Org. Biomol. Chem. 2018, 16, 12631271. 10.1039/C7OB03017G.
S. C.Kim; D. K.O’Flaherty; L.Zhou; V. S.Lelyveld; J. W.Szostak Inosine, but None of the 8-Oxo-Purines, Is a Plausible Component of a Primordial Version of RNA. Proc. Natl. Acad. Sci. U. S. A. 2018, 115, 1331813323. 10.1073/pnas.1814367115.
J.Xu; M.Tsanakopoulou; C. J.Magnani; R.Szabla; J. E.Šponer; J.Šponer; R. W.Góra; J. D.Sutherland A Prebiotically Plausible Synthesis of Pyrimidine β-Ribonucleosides and Their Phosphate Derivatives Involving Photoanomerization. Nat. Chem. 2017, 9, 303309. 10.1038/nchem.2664.
J.Xu; N. J.Green; C.Gibard; R.Krishnamurthy; J. D.Sutherland Prebiotic Phosphorylation of 2-Thiouridine Provides Either Nucleotides or DNA Building Blocks via Photoreduction. Nat. Chem. 2019, 11, 457462. 10.1038/s41557-019-0225-x.
H. J.Kim; S. A.Benner Prebiotic Stereoselective Synthesis of Purine and Noncanonical Pyrimidine Nucleotide from Nucleobases and Phosphorylated Carbohydrates. Proc. Natl. Acad. Sci. U. S. A. 2017, 114 (43), 1131511320. 10.1073/pnas.1710778114.
S.Islam; J. A.Aguilar; M. W.Powner; M.Nilsson; G. A.Morris; J. D.Sutherland Detection of Potential TNA and RNA Nucleoside Precursors in a Prebiotic Mixture by Pure Shift Diffusion-Ordered NMR Spectroscopy. Chem. - Eur. J. 2013, 19 (14), 45864595. 10.1002/chem.201202649.
B. D.Heuberger; A.Pal; F.Del Frate; V. V.Topkar; J. W.Szostak Replacing Uridine with 2-Thiouridine Enhances the Rate and Fidelity of Nonenzymatic RNA Primer Extension. J. Am. Chem. Soc. 2015, 137 (7), 27692775. 10.1021/jacs.5b00445.
R. K.Kumar; D. R.Davis Synthesis and Studies on the Effect of 2-Thiouridine and 4-Thiouridine on Sugar Conformation and RNA Duplex Stability. Nucleic Acids Res. 1997, 25, 127210.1093/nar/25.6.1272.
S. M.Testa; M. D.Disney; D. H.Turner; R.Kierzek Thermodynamics of RNA-RNA Duplexes with 2- or 4-Thiouridines: Implications for Antisense Design and Targeting a Group I Intron. Biochemistry 1999, 38 (50), 1665516662. 10.1021/bi991187d.
A. T.Larsen; A. C.Fahrenbach; J.Sheng; J.Pian; J. W.Szostak Thermodynamic Insights into 2-Thiouridine-Enhanced RNA Hybridization. Nucleic Acids Res. 2015, 43 (16), 76757687. 10.1093/nar/gkv761.
J. S.Teichert; F. M.Kruse; O.Trapp Direct Prebiotic Pathway to DNA Nucleosides. Angew. Chem., Int. Ed. 2019, 58, 99449947. 10.1002/anie.201903400.
I. A.Kozlov; P. K.Politis; A.Van Aerschot; R.Busson; P.Herdewijn; L. E.Orgel Nonenzymatic Synthesis of RNA and DNA Oligomers on Hexitol Nucleic Acid Templates: The Importance of the a Structure. J. Am. Chem. Soc. 1999, 121 (12), 26532656. 10.1021/ja983958r.
J. K.Watts; N.Martín-Pintado; I.Gómez-Pinto; J.Schwartzentruber; G.Portella; M.Orozco; C.González; M. J.Damha Differential Stability of 2′F-ANA•RNA and ANA•RNA Hybrid Duplexes: Roles of Structure, Pseudohydrogen Bonding, Hydration, Ion Uptake and Flexibility. Nucleic Acids Res. 2010, 38 (7), 24982511. 10.1093/nar/gkp1225.
J. V.Gavette; M.Stoop; N. V.Hud; R.Krishnamurthy RNA−DNA Chimeras in the Context of an RNA World Transition to an RNA/DNA World. Angew. Chem., Int. Ed. 2016, 55, 1320413209. 10.1002/anie.201607919.
S. J.Roberts; R.Szabla; Z. R.Todd; S.Stairs; D.-K.Bučar; J.Šponer; D. D.Sasselov; M. W.Powner Selective Prebiotic Conversion of Pyrimidine and Purine Anhydronucleosides into Watson-Crick Base-Pairing Arabino-Furanosyl Nucleosides in Water. Nat. Commun. 2018, 9 (1), 407310.1038/s41467-018-06374-z.
N.Martín-Pintado; M.Yahyaee-Anzahaee; R.Campos-Olivas; A. M.Noronha; C. J.Wilds; M. J.Damha; C.González The Solution Structure of Double Helical Arabino Nucleic Acids (ANA and 2′F-ANA): Effect of Arabinoses in Duplex-Hairpin Interconversion. Nucleic Acids Res. 2012, 40 (18), 93299339. 10.1093/nar/gks672.
A. M.Noronha; C. J.Wilds; C. N.Lok; K.Viazovkina; D.Arion; M. A.Parniak; M. J.Damha Synthesis and Biophysical Properties of Arabinonucleic Acids (ANA): Circular Dichroic Spectra, Melting Temperatures, and Ribonuclease H Susceptibility of ANA·RNA Hybrid Duplexes. Biochemistry 2000, 39 (24), 70507062. 10.1021/bi000280v.
F.Li; S.Sarkhel; C. J.Wilds; Z.Wawrzak; T. P.Prakash; M.Manoharan; M.Egli 2′-Fluoroarabino- and Arabinonucleic Acid Show Different Conformations, Resulting in Deviating RNA Affinities and Processing of Their Heteroduplexes with RNA by RNase H. Biochemistry 2006, 45 (13), 41414152. 10.1021/bi052322r.
A. I.Taylor; V. B.Pinheiro; M. J.Smola; A. S.Morgunov; S.Peak-Chew; C.Cozens; K. M.Weeks; P.Herdewijn; P.Holliger Catalysts from Synthetic Genetic Polymers. Nature 2015, 518, 427430. 10.1038/nature13982.
A. C.Fahrenbach; C.Giurgiu; C. P.Tam; L.Li; Y.Hongo; M.Aono; J. W.Szostak Common and Potentially Prebiotic Origin for Precursors of Nucleotide Synthesis and Activation. J. Am. Chem. Soc. 2017, 139 (26), 87808783. 10.1021/jacs.7b01562.
N.Prywes; J. C.Blain; F.Del Frate; J. W.Szostak Nonenzymatic Copying of RNA Templates Containing All Four Letters Is Catalyzed by Activated Oligonucleotides. eLife 2016, 5, 1775610.7554/eLife.17756.
T.Walton; J. W.Szostak A Highly Reactive Imidazolium-Bridged Dinucleotide Intermediate in Nonenzymatic RNA Primer Extension. J. Am. Chem. Soc. 2016, 138 (36), 1199612002. 10.1021/jacs.6b07977.
W.Zhang; C. P.Tam; L.Zhou; S. S.Oh; J.Wang; J. W.Szostak Structural Rationale for the Enhanced Catalysis of Nonenzymatic RNA Primer Extension by a Downstream Oligonucleotide. J. Am. Chem. Soc. 2018, 140 (8), 28292840. 10.1021/jacs.7b11750.
D.Venkateswarlu; D. M.Ferguson Effects of C2’-Substitution on Arabinonucleic Acid Structure and Conformation. J. Am. Chem. Soc. 1999, 121 (23), 56095610. 10.1021/ja990870k.
A.Pal; R. S.Das; W.Zhang; M.Lang; L. W.McLaughlin; J. W.Szostak Effect of Terminal 3′-Hydroxymethyl Modification of an RNA Primer on Nonenzymatic Primer Extension. Chem. Commun. 2016, 52, 1190511907. 10.1039/C6CC06925H.
W.Zhang; C. P.Tam; J.Wang; J. W.Szostak Unusual Base-Pairing Interactions in Monomer-Template Complexes. ACS Cent. Sci. 2016, 2 (12), 916926. 10.1021/acscentsci.6b00278.
W.Zhang; T.Walton; L.Li; J. W.Szostak Crystallographic Observation of Nonenzymatic RNA Primer Extension. eLife 2018, 7, 3642210.7554/eLife.36422.
M. J.Damha; A. M.Noronha; C. J.Wilds; J. F.Trempe; A.Denisov; R. T.Pon; K.Gehring Properties of Arabinonucleic Acids (ANA & 20′F-ANA): Implications for the Design of Antisense Therapeutics That Invoke RNase H Cleavage of RNA. Nucleosides, Nucleotides Nucleic Acids 2001, 20, 429440. 10.1081/NCN-100002317.
J. P.Schrum; A.Ricardo; M.Krishnamurthy; J. C.Blain; J. W.Szostak Efficient and Rapid Template-Directed Nucleic Acid Copying Using 2′-Amino-2′,3′-Dideoxyribonucleoside-5′- Phosphorimidazolide Monomers. J. Am. Chem. Soc. 2009, 131 (40), 1456014570. 10.1021/ja906557v.
S.Zhang; N.Zhang; J. C.Blain; J. W.Szostak Synthesis of N3′-P5′-Linked Phosphoramidate DNA by Nonenzymatic Template-Directed Primer Extension. J. Am. Chem. Soc. 2013, 135 (2), 924932. 10.1021/ja311164j.
S.Winstein; N. J.Holness Neighboring Carbon and Hydrogen. XIX. t-Butylcyclohexyl Derivatives. Quantitative Conformational Analysis. J. Am. Chem. Soc. 1955, 77 (21), 55625578. 10.1021/ja01626a037.
L.Li; N.Prywes; C. P.Tam; D. K.O’Flaherty; V. S.Lelyveld; E. C.Izgu; A.Pal; J. W.Szostak Enhanced Nonenzymatic RNA Copying with 2-Aminoimidazole Activated Nucleotides. J. Am. Chem. Soc. 2017, 139 (5), 18101813. 10.1021/jacs.6b13148.
H.Sawai; L. E.Orgel Oligonucleotide Synthesis Catalyzed by the Zn2+ Ion. J. Am. Chem. Soc. 1975, 97 (12), 35323533. 10.1021/ja00845a050.
A.Kanavarioti; P.-A.Monnard; D. W.Deamer Eutectic Phases in Ice Facilitate Nonenzymatic Nucleic Acid Synthesis. Astrobiology 2001, 1 (3), 271281. 10.1089/15311070152757465.
R.Krishnamurthy On the Emergence of RNA. Isr. J. Chem. 2015, 55, 837850. 10.1002/ijch.201400180.
S.Bhowmik; R.Krishnamurthy The Role of Sugar-Backbone Heterogeneity and Chimeras in the Simultaneous Emergence of RNA and DNA. Nat. Chem. 2019, 11, 10091018. 10.1038/s41557-019-0322-x.
J. W.Szostak The Eightfold Path to Non-Enzymatic RNA Replication. J. Syst. Chem. 2012, 3, 210.1186/1759-2208-3-2.
L.Zhou; S. C.Kim; K. H.Ho; D. K.O’Flaherty; C.Giurgiu; T. H.Wright; J. W.Szostak Non-Enzymatic Primer Extension with Strand Displacement. eLife 2019, 8, 5188810.7554/eLife.51888.
L.Zhou; D. K.O’Flaherty; J. W.Szostak Assembly of a Functional Ribozyme from Short Oligomers by Enhanced Non-Enzymatic Ligation. bioRxiv 2019, 700229.
A. E.Engelhart; K. P.Adamala; J. W.Szostak A Simple Physical Mechanism Enables Homeostasis in Primitive Cells. Nat. Chem. 2016, 8, 448453. 10.1038/nchem.2475.
C.Giurgiu; L.Li; D. K.O’Flaherty; C. P.Tam; J. W.Szostak A Mechanistic Explanation for the Regioselectivity of Nonenzymatic RNA Primer Extension. J. Am. Chem. Soc. 2017, 139 (46), 1674116747. 10.1021/jacs.7b08784.
T. H.Wright; C.Giurgiu; W.Zhang; A.Radakovic; D. K.O’Flaherty; L.Zhou; J. W.Szostak Prebiotically Plausible “Patching” of RNA Backbone Cleavage Through a 3′-5′ Pyrophosphate Linkage. J. Am. Chem. Soc. 2019, 141 (45), 1810418112. 10.1021/jacs.9b08237.
S. G.Trevino; N.Zhang; M. P.Elenko; A.Lupták; J. W.Szostak Evolution of Functional Nucleic Acids in the Presence of Nonheritable Backbone Heterogeneity. Proc. Natl. Acad. Sci. U. S. A. 2011, 108, 1349213497. 10.1073/pnas.1107113108.