Search

Protein and solid-state nanopores are used for DNA/RNA sequencing as well as for single molecule analysis. We proposed that selective labeling/tagging may improve base-to-base resolution of nucleic acids via nanopores. We have explored one specific tag, the Osmium tetroxide 2,2′-bipyridine (OsBp), which conjugates to pyrimidines and leaves purines intact. Earlier reports using OsBp-tagged oligodeoxyribonucleotides demonstrated proof-of-principle during unassisted voltage-driven translocation via either alpha-Hemolysin or a solid-state nanopore. Here we extend this work to RNA oligos and a third nanopore by employing the MinION, a commercially available device from Oxford Nanopore Technologies (ONT). Conductance measurements demonstrate that the MinION visibly discriminates oligoriboadenylates with sequence A15PyA15, where Py is an OsBp-tagged pyrimidine. Such resolution rivals traditional chromatography, suggesting that nanopore devices could be exploited for the characterization of RNA oligos and microRNAs enhanced by selective labeling. The data also reveal marked discrimination between a single pyrimidine and two consecutive pyrimidines in OsBp-tagged AnPyAn and AnPyPyAn. This observation leads to the conjecture that the MinION/OsBp platform senses a 2-nucleotide sequence, in contrast to the reported 5-nucleotide sequence with native nucleic acids. Such improvement in sensing, enabled by the presence of OsBp, may enhance base-calling accuracy in enzyme-assisted DNA/RNA sequencing.

Nanopores can serve as single molecule sensors. We exploited the MinION, a portable nanopore device from Oxford Nanopore Technologies (ONT), and repurposed it to detect any DNA/RNA oligo (target) in a complex mixture by conducting voltage-driven ion-channel measurements. The detection and quantitation of the target is enabled by the use of a unique complementary probe. Using a validated labeling technology, probes are tagged with a bulky Osmium tag (Osmium tetroxide 2,2’-bipyridine), in a way that preserves strong hybridization between probe and target. Intact oligos traverse the MinION’s nanopore relatively quickly compared to the device’s acquisition rate, and exhibit count of events comparable to the baseline. Counts are reported by a publicly available software,OsBp_detect. Due to the presence of the bulky Osmium tag, probes traverse more slowly, produce multiple counts over the baseline, and are even detected at single digit attomole (amole) range. In the presence of the target the probe is “silenced”. Silencing is attributed to a 1:1 double stranded (ds) complex that doesn’t fit and can’t traverse this nanopore. This ready-to-use platform can be tailored as a diagnostic test to meet the requirements for point-of-care cell-free tumor DNA (ctDNA) and microRNA (miRNA) detection and quantitation in body fluids.

The influence of an electric field on an isolated channel or nanopore separating two compartments filled with electrolytes produces a constant ion flux through the pore. Nucleic acids added to one compartment traverse the pore, and modulate the current in a sequence-dependent manner. While translocation is faster than detection, the α-hemolysin nanopore (α-HL) successfully senses base modifications in ssDNA immobilized within the pore. With the assistance of a processing enzyme to slow down translocation, nanopore-based DNA sequencing is now a commercially available platform. However, accurate base calling is challenging because α-HL senses a sequence, and not a single nucleotide. Osmylated DNA was recently proposed as a surrogate for nanopore-based sequencing. Osmylation is the addition of osmium tetroxide 2,2’-bipyridine (OsBp) to the C5–C6 pyrimidine double bond. The process is simple, selective for deoxythymidine (dT) over deoxycytidine (dC), unreactive towards the purines, practically 100% effective, and strikingly independent of length, sequence, and composition. Translocation of an oligodeoxynucleotide (oligo) dA10XdA9 via α-HL is relatively slow, and exhibits distinct duration as well as distinct residual current when X = dA, dT(OsBp), or dC(OsBp). The data indicate that the α-HL constriction zone/β-barrel interacts strongly with both OsBp and the base. A 23 nucleotide long oligo with four dT(OsBp) traverses 18-times slower, and the same oligo with nine (dT+dC)(OsBp) moieties traverses 84-times slower compared to dA20, suggesting an average rate of 40 or 180 μs/base, respectively. These translocation speeds are well above detection limits, may be further optimized, and clear the way for nanopore-based sequencing using osmylated DNA.

Osmium tetroxide 2,2’-bipyridine (OsBp) is known to react with pyrimidines in ssDNA and preferentially label deoxythymine (T) over deoxycytosine (C). The product, osmylated DNA, was proposed as a surrogate for nanopore-based DNA sequencing due to OsBp’s “perfect” label attributes. Osmylated deoxyoligos translocate unassisted and measurably slow via sub-2 nm SiN solid-state nanopores, as well as via the alpha-hemolysin (α-HL) pore. Both nanopores discriminate clearly between osmylated and intact nucleobase; α-HL was also shown to discriminate between osmylated T and osmylated C. Experiments presented here confirm that the kinetics of osmylation are comparable for short oligos and long ssDNA and show that pyrimidine osmylation is practically complete in two hours at room temperature with less than 15 mM OsBp. Under the proposed labeling conditions: deoxyoligo backbone degradation measures less than 1/1,000,000; false positives such as osmylated deoxyadenine (A) and osmylated deoxyguanine (G) measure less than 1/100,000; false negatives, i.e., unosmylated C measure less than 1/10,000; and unosmylated T must measure substantially lower than 1/10,000 due to the 27-fold higher reactivity of T compared to C. However, osmylated C undergoes degradation that amounts to about 1–2% for the duration of the labeling protocol. This degradation may be further characterized, possibly suppressed, and the properties of the degradation products via nanopore translocation can be evaluated to assure base calling quality in a DNA sequencing effort.

With the recent advances in electron microscopy (EM), computation, and nanofabrication, the original idea of reading DNA sequence directly from an image can now be tested. One approach is to develop heavy atom labels that can provide the contrast required for EM imaging. While evaluating tentative labels for the respective nucleobases in synthetic oligodeoxynucleotides (oligos), we developed a streamlined CE protocol to assess the label stability, reactivity, and selectivity. We report our protocol using osmium tetroxide 2,2'-bipyridine (Osbipy) as a thymidine (T) specific label. The observed rates show that the labeling process is kinetically independent of both the oligo length, and the base composition. The conditions, i.e. temperature, optimal Osbipy concentration, and molar ratio of reagents, to promote 100% conversion of the starting oligo to labeled product were established. Hence, the optimized conditions developed with the oligos could be leveraged to allow osmylation of effectively all Ts in ssDNA, while achieving minimal mislabeling. In addition, the approach and methods employed here may be adapted to the evaluation of other prospective contrasting agents/labels to facilitate next-generation DNA sequencing by EM.

Saenger sequencing has led the advances in molecular biology, while faster and cheaper next generation technologies are urgently needed. A newer approach exploits nanopores, natural or solid-state, set in an electrical field, and obtains base sequence information from current variations due to the passage of a ssDNA molecule through the pore. A hurdle in this approach is the fact that the four bases are chemically comparable to each other which leads to small differences in current obstruction. 'Base calling' becomes even more challenging because most nanopores sense a short sequence and not individual bases. Perhaps sequencing DNA via nanopores would be more manageable, if only the bases were two, and chemically very different from each other; a sequence of 1s and 0s comes to mind. Osmylated DNA comes close to such a sequence of 1s and 0s. Osmylation is the addition of osmium tetroxide bipyridine across the C5-C6 double bond of the pyrimidines. Osmylation adds almost 400% mass to the reactive base, creates a sterically and electronically notably different molecule, labeled 1, compared to the unreactive purines, labeled 0. If osmylated DNA were successfully sequenced, the result would be a sequence of osmylated pyrimidines (1), and purines (0), and not of the actual nucleobases. To solve this problem we studied the osmylation reaction with short oligos and with M13mp18, a long ssDNA, developed a UV-vis assay to measure extent of osmylation, and designed two protocols. Protocol A uses mild conditions and yields osmylated thymidines (1), while leaving the other three bases (0) practically intact. Protocol B uses harsher conditions and effectively osmylates both pyrimidines, but not the purines. Applying these two protocols also to the complementary of the target polynucleotide yields a total of four osmylated strands that collectively could define the actual base sequence of the target DNA.

Nanopores are a promising platform in next generation DNA sequencing. In this platform, an individual DNA strand is threaded into nanopore using an electric field, and enzyme-based ratcheting is used to move the strand through the detector. During this process the residual ion current through the pore is measured, which exhibits unique levels for different base combinations inside the pore. While this approach has shown great promise, accuracy is not optimal because the four bases are chemically comparable to one another, leading to small differences in current obstruction. Nucleobase-specific chemical tagging can be a viable approach to enhancing the contrast between different bases in the sequence. Herein we show that covalent modification of one or both of the pyrimidine bases by an osmium bipyridine complex leads to measureable differences in the blockade amplitudes of DNA molecules. We qualitatively determine the degree of osmylation of a DNA strand by passing it through a solid-state nanopore, and are thus able to gauge T and C base content. In addition, we show that osmium bipyridine reacts with dsDNA, leading to substantially different current blockade levels than exhibited for bare dsDNA. This work serves as a proof of principle for nanopore sequencing and mapping via base-specific DNA osmylation.

Polymeric compounds similar to oligonucleotides are relevant to the origin of life and particularly to the concept of an RNA world. Although short oligomers of RNA can be synthesized nonenzymatically under laboratory conditions by second-order reactions in concentrated solutions, there is no consensus on how these polymers could have been synthesized de novo on the early Earth from dilute solutions of monomers. To address this question in the context of an RNA world, we have explored ice eutectic phases as a reaction medium. When an aqueous solution freezes, the solutes become concentrated in the spaces between the ice crystals. The increased concentration offsets the effect of the lower temperature and accelerates the reaction. Here we show that in the presence of metal ions in dilute solutions, frozen samples of phosphoimidazolide-activated uridine react within days at -18 degrees C to form oligouridylates up to 11 bases long. Product yields typically exceed 90%, and approximately 30% of the oligomers include one or more 3'-5' linkages. These conditions facilitate not only the notoriously difficult oligouridylate synthesis, but also the oligomerization of activated cytidylate, adenylate, and guanylate. To our knowledge, this represents the first report to indicate that ice matrices on the early Earth may have accelerated certain prebiotic polymerization reactions.

Polycytidylate, poly(C), serves as a scaffold or template to direct and catalyze the synthesis of long oligoguanylates from guanosine 5'-phosphate 2-methylimidazolide, 2-MeImpG. In the absence of poly(C), small amounts of three isomeric dimers, i.e., the 2'-5'-, the 3'-5'-, and the pyrophosphate-linked, are formed slowly. In the presence of poly(C) oligomers that are primarily 3'-5'-linked are formed quickly and in high yield. Product analysis suggests that the oligomers are elongation products of the 3'-5'-linked dimer, abbreviated D. Assuming that D is formed slowly from two molecules of 2-MeImpG (Scheme 1) and elongates relatively fast, the initial rate of dimerization, d[D]/dt in M h(-1), was determined using two independent methods. The first method is based on the approximation that at the onset of the reaction the substrate is consumed only via hydrolysis and dimerization, and thus elongation can be neglected. The second, more accurate, method exploits the assertion that every oligomer was once a 3'-5'-linked dimer. Hence the concentration of D was obtained indirectly from the concentration of the oligomer products. These two methods gave comparable results. Experiments were run in aqueous solution in the presence of 1.0 M NaCl, 0.2 M MgCl(2) at pH 7.9 +/- 0.1 and 23 degrees C. Controls were run in the absence of poly(C) and in the presence of other polynucleotides. The kinetics were determined as a function of both monomer and polymer concentration the latter expressed in C equivalents. The kinetic data obtained in the presence of poly(C) confirmed an earlier conclusion regarding the remarkable effect of poly(C) on the formation of the 3'-5'-linked diguanylate. Initial dimerization rates were quantitatively correlated using a simple template-directed (TD) model that presumes cooperative binding (two association constants) of 2-MeImpG on poly(C) and reaction between adjacent template-bound molecules. The model allows for the estimation of the association constants and the intrinsic rate constant of dimerization, k(2). Insights into the detailed mechanism are also gained from this analysis. The fact that the proposed model can successfully correlate kinetic data that vary by more than 5000-fold between the slowest and the fastest reaction adds confidence and suggests the suitability of this model for describing TD reactions in general. It is anticipated that similar analysis of other known TD reactions may lead to clues that will facilitate the design of more efficient polynucleotide-synthesizing systems.

Aqueous solutions of deoxyguanosine 5'-monophosphate 2-methylimidazolide, 2-MeImpdG, yield primarily deoxyguanosine 5'-monophosphate, 5'dGMP, and pyrophosphate-linked dideoxyguanylate, dG5'ppdG, abbreviated G2p (see Chart 1). The initial rate of G2p formation, d[G2p]/dt in M h-1, determined at 23 degrees C, pH 7.8, 1.0 M NaCl and 0.2 M Mg2+ by timed high-performance liquid chromatography (HPLC) analysis, exhibits a second-order dependence on 2-MeImpdG concentration, [G]o, indicating a bimolecular mechanism of dimerization in the range 0.02 M < or = [G]o < or = 0.09 M. In the presence of polycytidylate, poly(C), G2p synthesis is accelerated and oligodeoxyguanylate products are formed by incorporation of 2-MeImpdG molecules. The kinetics of G2p formation as a function of both monomer and polymer concentration, expressed in C equivalents, were also determined under the above conditions and exhibited a complex behavior. Specifically, at a constant [poly(C)], values of d[G2p]/dt typically increased with [G]o with a parabolic upward curvature. At a constant [G]o, values of d[G2p]/dt increase with [poly(C)], but level off at the higher poly(C) concentrations. As [G]o increases this saturation occurs at a higher poly(C) concentration, a result opposite to expectation for a simple complexation of two reacting monomers with the catalyst prior to reaction. Nevertheless, these results are shown to be quantitatively consistent with a template-directed (TD) mechanism of dimerization where poly(C) acts as the template to bind 2-MeImpdG in a cooperative manner and lead, for the first time, to the formulation of principles that govern template-directed chemistry. Analysis of the kinetic data via a proposed TD cooperative model provides association constants for the affinity between polymer and monomer and the intrinsic reactivity of 2-MeImpdG toward pyrophosphate synthesis. To the best of our knowledge, poly(C)/2-MeImpdG is the first system that could serve as a textbook example of a TD reaction under conditions such that the template is fully saturated by monomers and under conditions that it is not.

Phosphoimidazolide activated ribomononucleotides (*pN, see structure) are useful substrates for the non-enzymatic synthesis of oligonucleotides. In the presence of metal ions, aqueous solutions of *pN yield primarily the two internucleotide-linked (pN2' pN and pN3' pN) and the pyrophosphate-linked (N5' ppN) dimers. Small amounts of cyclic dimers and higher oligomers are also produced. In this study the relative reactivity of 2'-OH vs. 3'-OH was determined from the ratio of the yields of pN2' pN vs. pN3' pN. Experiments were performed at 23 degrees C in the range 7.2 < or = pH < or = 8.4 with substrates that differ in nucleobase (guanosine (G), cytidine (C), uridine (U), and adenosine (A)) and leaving group (imidazole (Im), 2-methylimidazole (2-MeIm) and 2,4-dimethylimidazole (2,4-diMeIm)). Two metal ions (Mg2+ or Mn2+) were employed as catalysts. The conditions used here, i.e. a substrate concentration in the range 0.1 M to 1.0 M and metal ion concentration in the range 0.05 M to 0.2 M, favor base-stacking interactions. The ratio pN2' pN: pN3' pN = 2'-5': 3'-5' was found independent of nucleobase and typically varied between 2 to 3 indicating that the 2'-OH is about 2 to 3 times more reactive than the 3'-OH. *pN with Im, compared to 2-MeIm and 2,4-diMeIm leaving group, produce lower yields of internucleotide linked dimers, and a higher pN2' pN: pN3' pN ratio. Trends in the data, observed with all three leaving groups, suggest an increase in pN2' pN: pN3' pN ratio with decreasing substrate concentration (up to 5.47 with 0.051 M ImpG). The observations are in accord with earlier studies reporting a relative reactivity 2'-5': 3'-5' = 6 to 9 obtained with Im as the leaving group, in dilute nucleotide solutions and under conditions that disfavor stacking. It is speculated that the concentration induced change in the relative reactivity is the result of self-association via base-stacking that enhances selectively the proximity of the 3'-OH of one molecule to the reactive P-N bond of an other molecule. The implication of these conclusions for oligomerization/ligation reactions is discussed.

To identify key parameters which influence the efficiency of nonenzymatic template-directed oligonucleotide reactions, a kinetic study of oligoguanylate synthesis on a polycytidylate (poly(C)) template has been performed. This oligomerization is satisfactorily described by three kinetic processes: (i) dimerization to form the first primer (k2), (ii) elongation of a preformed primer (ki, i ≥ 3), and (iii) hydrolysis of the monomer to form deactivated material (kh), with kh < k2 < ki. This is the first reported study that includes rate determinations of ki as a function of the concentration of both poly(C) template and the activated monomer, guanosine 5‘-monophosphate-2-methylimidazolide (2-MeImpG), in the range 2 mM ≤ [poly(C)] ≤ 50 mM and 5 mM ≤ [2-MeImpG] ≤ 50 mM. ki values determined under conditions where the template is fully saturated with monomer are practically independent of both monomer and polymer concentration and thus strongly support a template-directed elongation model. Values of ki determi...

Phosphoimidazolide activated ribomononucleotides (*pN) are useful substrates for the non-enzymatic synthesis of polynucleotides. However, dilute neutral aqueous solutions of *pN typically yield small amounts of dimers and traces of polymers; most of *pN hydrolyzes to yield nucleoside 5'-monophosphate. Here we report the self-condensation of nucleoside 5'-phosphate 2- methylimidazolide (2-MeImpN with N = cytidine, uridine or guanosine) in the presence of Mg2(+) in concentrated solutions, such as might have been found in an evaporating lagoon on prebiotic Earth. The product distribution indicates that oligomerization is favored at the expense of hydrolysis. At 1.0 M, 2-MelmpU and 2-MelmpC produce about 65% of oligomers including 4% of the 3',5'-Iinked dimer. Examination of the product distribution of the three isomeric dimers in a self-condensation allows identification of reaction pathways that lead to dimer formation. Condensations in a concentrated mixture of all three nucleotides (U,C,G mixtures) is made possible by the enhanced solubility of 2-MeImpG in such mixtures. Although percent yield of intemucleotide linked dimers is enhanced as a function of initial monomer concentration, pyrophosphate dimer yields remain practically unchanged at about 20% for 2-MelmpU, 16% for 2-MeImpC and 25% of the total pyrophosphate in the U,C,G mixtures. The efficiency by which oligomers are produced in these concentrated solutions makes the evaporating lagoon scenario a potentially interesting medium for the prebiotic synthesis of dimers and short RNAs.

Guanosine 5'-phosphate 2-methylimidazolide (2-MeImpG), a labile phosphoimidazolide analog of guanosine triphosphate, was used to test the reactivity of the natural polyamines (PAs), spermine (spm) and spermidine (spd). The products are the guanosine 5'-phosphate-polyamine derivatives (PA-pG: spd-pG and spm-pG) which are quite stable in the range 4 < pH < 11. Our study is the first of which we are aware that reports on the nucleophilicity of these amines. The main findings are as follows. (i) HPLC analysis of the products indicates the formation of only two of the three possible spd products and only one of the two possible spm products. These results can be explained if only the primary amino groups of the two polyamines are reactive, while the secondary amino groups are rendered unreactive by a steric effect. The reactions of 2-MeImpG and other phosphoimidazolide derivatives of nucleosides (ImpNs) with primary and secondary monoamines support this interpretation (Kanavarioti et al. J. Org. Chem. 1995, 60, 632). (ii) The product ratio of the two spd-pG adducts derived from the primary amino groups varies between 2.40 and 0.71 in the range 6.1 < or equal to pH < or equal to 11.9. Such small variation in the product ratio can only be rationalized by the similar, but not identical, basicity of the two primary amino groups and provides strong support for a previously reported model for polyamine ionization (Onasch et. al. Biophys. Chem. 1984, 19, 245). (iii) On the basis of our kinetic determinations conditions at which the nucleophilicity of these amines is at a minimum and at which other interactions with ImpNs could be tested can be chosen.

Symmetrical dinucleoside 5'-pyrophosphates have been synthesized from the corresponding nucleoside 5'-phosphate free acid in high yield. The one-pot procedure is carried out in DMF or DMSO using triphenylphosphine and 2,2'-dipyridyldisulfide as the coupling agents, and 1-methylimidazole as the catalyst.

The formation of the internucleotide bond in diguanylate synthesis was studied in aqueous solution at pH 8 and 0.2 M Mg2+ in the presence and absence of polycytidylate, poly(C). The investigation was simplified by using guanosine 5'-phosphorylmorpholinamide, mor-pG, which can act only as a nucleophile, and deoxyguanosine 5'-phosphoryl-2-methylimidazolide, 2-MeImpdG, which can act only as an electrophile. The time-dependent product distribution was monitored by high-performance liquid chromatography (HPLC) and liquid chromatography mass spectrometry (LC/MS). In the absence of poly(C) the reaction between mor-pG and 2-MeImpdG yielded small amounts of the dimer mor-pGpdG with a regioselectivity of 2'-5':3'-5' = 3.5. In the presence of poly(C) dimer yields increased and a reversal in regioselectivity occurred; both effects were in proportion to the concentration of the polymer. The results can be quantitatively explained with the proposition that poly(C), acting as the template, catalyzes the reaction between template-bound monomers by about a factor of 4-5 over the reaction in solution and yields dimers with a regioselectivity of 2'-5':3'-5' approximately 0.33. These findings illustrate the intrinsic preference of guanosine monomers to correctly self-assemble on the appropriate template.

An earlier study of the reaction of phosphoimidazolide activated nucleosides (ImpN) in aqueous phosphate buffers indicated two modes of reaction of the phosphate monoanion and dianion. The first mode is catalysis of the hydrolysis of the P-N bond in ImpN's which leads to imidazole and nucleoside 5'-monophosphate. The second represents a nucleophilic substitution of the imidazole to yield the nucleoside 5'-diphosphate. This earlier study thus served as a model for the reaction of ImpN with nucleoside monophosphates (pN) because the latter can be regarded as phosphate derivatives. In the present study we investigated the reaction of guanosine 5'-phosphate-2-methylimidazolide, 2-MeImpG, in the presence of pN (N = guanosine, adenosine and uridine) in the range 6.9 less than or equal to pH less than or equal to 7.7. We observed that pN's do act as nucleophiles to form NppG, and as general base to enhance the hydrolysis of the P-N bond in 2-MeImpG, i.e. pN show the same behavior as inorganic phosphate. The kinetic analysis yields the following rate constants for the dianion pN2-: knpN = 0.17 +/- 0.02 M-1 h-1 for nucleophilic attack and khpN = 0.11 +/- 0.07 M-1 h-1 for general base catalysis of the hydrolysis. These rate constants which are independent of the nucleobase compare with kp.2 = 0.415 M-1 h-1 and khp2. = 0.217 M-1 h-1 for the reactions of HPO4(2-). In addition, this study shows that under conditions where pN presumably form stacks, the reaction mechanism remains unchanged although in quantitative terms stacked pN are somewhat less reactive. Attack by the 2'-OH and 3'-OH groups of the ribose moiety in amounts greater than or equal to 1% is not observed; this is attributed to the large difference in nucleophilicity in the neutral pH range between the phosphate group and the ribose hydroxyls. This nucleophilicity rank is not altered by stacking.

Phosphoimidazolide-activated derivatives of guanosine and cytidine 5'-monophosphates, henceforth called ImpN's, exhibit enhanced rates of degradation in the presence of aqueous inorganic phosphate in the range 4.0 < or = pH < or = 8.6. This degradation is been attributed to (i) nucleophilic substitution of the imidazolide and (ii) catalysis of the P-N bond hydrolysis by phosphate. The first reaction results in the formation of nucleoside 5'-diphosphate and the second in nucleoside 5'-monophosphate. Analysis of the observed rates as well as the product ratios as a function of pH and phosphate concentration allow distinction between various mechanistic possibilities. The results show that both H2PO4- and HPO4(2-) participate in both hydrolysis and nucleophilic substitution. Statistically corrected biomolecular rate constants indicate that the dianion is 4 times more effective as a general base than the monoanion, and 8 times more effective as nucleophile. The low Bronsted value beta = 0.15 calculated for these phosphate species, presumed to act as general bases in facilitating water attack, is consistent with the fact that catalysis of the hydrolysis of the P-N bond in ImpN's has not been detected before. The beta nuc = 0.35 calculated for water, H2PO4-, HPO4(2-), and hydroxide acting as nucleophiles indicates a more associative transition state for nucleotidyl (O2POR- with R = nucleoside) transfers than that observed for phosphoryl (PO3(2-)) transfers (beta nuc = 0.25). With respect to the stability/reactivity of ImpN's under prebiotic conditions, our study shows that these materials would not suffer additional degradation due to inorganic phosphate, assuming the concentrations of phosphate, Pi, on prebiotic Earth were similar to those in the present oceans ([Pi] approximately 2.25 micromoles).

Selected imidazolide-activated nucleotides have been subjected to hydrolysis under conditions similar to those that favor their template-directed oligomerization. Rate constants of hydrolysis of the P-N bond in guanosine 5'-monophosphate 2-methylimidazolide (2-MeImpG) and in guanosine 5'-monophosphate imidazolide (ImpG), kh, have been determined in the presence/absence of magnesium ion as a function of temperature and polycytidylate [poly(C)] concentration. Using the rate constant of hydrolysis of 2-MeImpG and the rate constant of elongation, i.e., the reaction of an oligoguanylate with 2-MeImpG in the presence of poly(C) acting as template, the limiting concentration of 2-MeImpG necessary for oligonucleotide elongation to compete with hydrolysis can be calculated. The limiting concentration is defined as the initial concentration of monomer that results in its equal consumption by hydrolysis and by elongation. These limiting concentrations of 2-MeImpG are found to be 1.7 mM at 37 degrees C and 0.36 mM at 1 degrees C. Boundary conditions in the form of limiting concentration of activated nucleotide may be used to evaluate a prebiotic model for chemical synthesis of biopolymers. For instance, the limiting concentration of monomer can be used as a basis of comparison among catalytic, but nonenzymatic, RNA-type systems. We also determined the rate constant of dimerization of 2-MeImpG, k2 = 0.45 +/- 0.06 M-1 h-1 in the absence of poly(C), and 0.45 +/- 0.06 less than or equal to k2 less than or equal to 0.97 +/- 0.13 M-1 h-1 in its presence at 37 degrees C and pH 7.95.(ABSTRACT TRUNCATED AT 250 WORDS)

The RNA world hypothesis requires a plausible mechanism by which RNA itself (or precursor RNA-like polymers) can be synthesized nonenzymatically from the corresponding building blocks. Simulation experiments have exploited chemically reactive mononucleotides as monomers. Solutions of such monomers in the prebiotic environment were likely to be very dilute, but in experimental simulations of polymerization reactions dilute solutions of activated mononucleotides in the millimolar range hydrolyze extensively, and only trace amounts of dimers and trimers are formed. We report here that random medium-size RNA analogues with mixed sequences (5- to 17-mers with traces of longer products) can be synthesized in ice eutectic phases that are produced when dilute solutions of activated monomers and catalysts (Mg(II) and Pb(II)) are frozen and maintained at -18 degrees C for periods up to 38 days. Under these conditions, the monomers are concentrated as eutectics in an ice matrix. Hydrolysis of the activated mononucleotides was suppressed at low-temperature ranges, and polymerization was enhanced with yields up to 90%. Analysis of the mixed oligomers established that incorporation of both purine and pyrimidine bases proceeded at comparable rates and yields. These results suggest that ice deposits on the early Earth could have facilitated the synthesis of short- and medium-size random sequence RNA analogues and thereby provided a microenvironment suitable for the formation of biopolymers or their precursors.

A commonly accepted view is that life began in a marine environment, which would imply the presence of inorganic ions such as Na+, Cl-, Mg2+, Ca2+, and Fe2+. We have investigated two processes relevant to the origin of life--membrane self-assembly and RNA polymerization--and established that both are adversely affected by ionic solute concentrations much lower than those of contemporary oceans. In particular, monocarboxylic acid vesicles, which are plausible models of primitive membrane systems, are completely disrupted by low concentrations of divalent cations, such as magnesium and calcium, and by high sodium chloride concentrations as well. Similarly, a nonenzymatic, nontemplated polymerization of activated RNA monomers in ice/eutectic phases (in a solution of low initial ionic strength) yields oligomers with80% of the original monomers incorporated, but polymerization in initially higher ionic strength aqueous solutions is markedly inhibited. These observations suggest that cellular life may not have begun in a marine environment because the abundance of ionic inorganic solutes would have significantly inhibited the chemical and physical processes that lead to self-assembly of more complex molecular systems.

The formation of the internucleotide bond in diguanylate synthesis was studied in aqueous solution at pH 8 and 0.2 M Mg2+ in the presence and absence of polycytidylate, poly(C). The investigation was simplified by using guanosine 5'-phosphorylmorpholinamide, mor-pG, which can act only as a nucleophile, and deoxyguanosine 5'-phosphoryl-2-methylimidazolide, 2-MeImpdG, which can act only as an electrophile. The time-dependent product distribution was monitored by high-performance liquid chromatography (HPLC) and liquid chromatography mass spectrometry (LC/MS). In the absence of poly(C) the reaction between mor-pG and 2-MeImpdG yielded small amounts of the dimer mor-pGpdG with a regioselectivity of 2'-5':3'-5' = 3.5. In the presence of poly(C) dimer yields increased and a reversal in regioselectivity occurred; both effects were in proportion to the concentration of the polymer. The results can be quantitatively explained with the proposition that poly(C), acting as the template, catalyzes the reaction between template-bound monomers by about a factor of 4-5 over the reaction in solution and yields dimers with a regioselectivity of 2'-5':3'-5' approximately 0.33. These findings illustrate the intrinsic preference of guanosine monomers to correctly self-assemble on the appropriate template.

Phosphoimidazolide-activated ribomononucleotides (*pN; see Scheme I) are useful substrates for the nonenzymatic synthesis of oligonucleotides. In the presence of metal ions dilute neutral aqueous solutions of *pN (0.01 M) typically yield only small amounts of dimers and traces of oligomers; most of *pN hydrolyzes to yield nucleoside 5'-monophosphate (5'NMP). An earlier investigation of *pN reactions in highly concentrated aqueous solutions (up to 1.4 M) showed, as expected, that the percentage yield of the condensation products increases and the yield of the hydrolysis product correspondingly decreases with *pN concentration (Kanavarioti 1997). Here we report product distributions in reactions with one, two, or three reactive components at the same total nucleotide concentration. *pN used as substrates were the nucleoside 5'-phosphate 2-methylimidazolides, 2-MeImpN, with N = cytidine (C), uridine (U), or guanosine (G). Reactions were conducted as self-condensations, i. e., one nucleotide only, with two components in the three binary U,C, U,G, and C,G mixtures, and with three components in the ternary U,C, G mixture. The products are 5'NMP, 5',5'-pyrophosphate-, 2',5'-, 3', 5'-linked dimers, cyclic dimers, and a small percentage of longer oligomers. The surprising finding was that, under identical conditions, including the same total monomer concentration, the product distribution differs substantially from one reaction to another, most likely due to changing intermolecular interactions depending on the constituents. Even more unexpected was the observed trend according to which reactions of the U,C,G mixture produce the highest yield of internucleotide-linked dimers, whereas the self-condensations produce the least and the reactions with the binary mixtures produce yields that fall in between. What is remarkable is that the approximately two-fold increase in the percentage yield of internucleotide-linked dimers is not due to a concentration effect or a catalyst, but to the increased complexity of the system from a single to two and three components. These observations, perhaps, provide an example of how increased complexity in relatively simple chemical systems leads to organization of the material and consequently to chemical evolution. A possible link between prebiotic chemistry and the postulated RNA world is discussed.

We have recently shown that the polycytidylic acid-directed polymerization of guanosine 5'-monophosphate 2-methylimidazolide (2-MeImpG) is amenable to kinetic study and that rate determinations as a function of 2-MeImpG concentration can reveal much mechanistic detail (Kanavarioti et al. 1993). Here we report kinetic data which show that, once the reaction has been initiated by the formation of dimers, the elongation of dimers to form longer oligomers is accelerated by decreasing polycytidylate (poly(C)) concentration from 0.05 to 0.002 M. This result is consistent with the previously proposed mechanism. The increase in the observed pseudo-first order rate constant for formation of the trimer, k3', and the corresponding constant for formation of oligomers longer than the trimer, ki' (ki' is independent of oligomer length for i > or = 4), with decreasing template concentration for a given monomer concentration is attributed to an increase in template occupancy as template concentration is reduced.

Evidence is presented for complexation of guanosine 5'-monophosphate 2-methylimidazolide (2-MeImpG) with polycytidylate (poly(C)) at pH 8.0 and 23 degrees C in the presence of 1.0 M NaCl2 and 0.2 M MgCl2 in water. The association of 2-MeImpG with poly(C) was investigated using UV-vis spectroscopy as well as by monitoring the kinetics of the nucleophilic substitution reaction of the imidazole moiety by amines. The results of both methods are consistent with moderately strong poly(C) 2-MeImpG complexation and the spectrophotometric measurements allowed the construction of a binding isotherm with a concentration of 2-MeImpG equal to 5.55 +/- 0.15 mM at half occupancy. UV spectroscopy was employed to establish the binding of other guanosine derivatives on poly(C). These derivatives are guanosine 5'-monophosphate (5'GMP), guanosine 5'-monophosphate imidazolide (ImpG), and guanosine 5'-monophosphate morpholidate (morpG). Within experimental error these guanosine derivatives exhibit the same affinity for poly(C) as 2-MeImpG.

A kinetic study of oligoguanylate synthesis on a polycytidylate template, poly(C), as a function of the concentration of the activated monomer, guanosine 5'-monophosphate 2-methylimidazolide, 2-MeImpG, is reported. Reactions were run with 0.005-0.045 M 2-MeImpG in the presence of 0.05 M poly(C) at 23 degrees C. The kinetic results are consistent with a reaction scheme (eq 1) that consists of a series of consecutive steps, each step representing the addition of one molecule of 2-MeImpG to the growing oligomer. This scheme allows the calculation of second-order rate constants for every step by analyzing the time-dependent growth of each oligomer. Computer simulations of the course of reaction based on the determined rate constants and eq 1 are in excellent agreement with the product distributions seen in the HPLC profiles. In accord with an earlier study (Fakhrai, H.; Inoue, T.; Orgel, L. E. Tetrahedron 1984, 40, 39), rate constants, ki, for the formation of the tetramer and longer oligomers up to the 16-mer were found to be independent of length and somewhat higher than k3 (formation of trimer), which in turn is much higher than k2 (formation of dimer). The ki (i > or = 4), k3, and k2 values are not true second-order rate constants but vary with monomer concentration. Mechanistic models for the dimerization (Scheme I) and elongation reactions (Scheme II) are proposed that are consistent with our results. These models take into account that the monomer associates with the template in a cooperative manner. Our kinetic analysis allowed the determination of rate constants for the elementary processes of covalent bond formation between two monomers (dimerization) and between an oligomer and a monomer (elongation) on the template. A major conclusion from our study is that bond formation between two monomer units or between a primer and a monomer is assisted by the presence of additional next-neighbor monomer units. This is consistent with recent findings with hairpin oligonucleotides (Wu, T.; Orgel, L. E. J. Am. Chem. Soc. 1992, 114, 317). Our study is the first of its kind that shows the feasibility of a thorough kinetic analysis of a template-directed oligomerization and provides a detailed mechanistic model of these reactions.

Symmetrical dinucleoside 5'-pyrophosphates have been synthesized from the corresponding nucleoside 5'-phosphate free acid in high yield. The one-pot procedure is carried out in DMF or DMSO using triphenylphosphine and 2,2'-dipyridyldisulfide as the coupling agents, and 1-methylimidazole as the catalyst.

A computer simulation (KINSIM) modeling up to 33 competing reactions was used in order to investigate the product distribution in a template-directed oligonucleotide synthesis as a function of time and concentration of the reactants. The study is focused on the poly(C)-directed elongation reaction of an oligoguanylate (a 7-mer is chosen) with guanosine 5′-monophosphate-2-methylimidazolide (2-MeImpG), the activated monomer. It is known that theelongation of oligoguanylates to form oligomeric products such as 8-mer, 9-mer, 10-mer, etc., is in competition with (1) thedimerization and further oligomerization reaction of 2-MeImpG that leads to the formation of dimers and short oligomers, and (2) thehydrolysis of 2-MeImpG that forms inactive guanosine 5′-monophosphate, 5′-GMP. Experimentally determined rate constants for the above three processes at 37°C and pH 7.95 were used in the simulation; the initial concentrations of 2-MeImpG, [M]o, and of the oligoguanylate primer, [7-mer]o, were varied, and KINSIM calculated the distribution of products as a function of time until equilibration was reached, i.e., when all the activated monomer has been consumed. In order to sort out how strongly the elongation reaction may be affected by the competing hydrolysis and dimerization, we also simulated the idealized situation in which these competing reactions do not occur. Simulation of the idealized system suggests that (1) the fraction of [7-mer]o that has reacted as well as the product distribution after equilibration do not depend on the absolute concentrations of the reactants, but only on their ratio, [M]o/[7-mer]o; (2) the rate of elongation is proportional to [7-mer]o and not to [M]o; and (3) as the [M]o/[7-mer]o ratio increases longer oligomers are formed. The results of the computer simulation with the experimental system, i.e., elongation in the presence of both hydrolysis and dimerization, are similar to the ones obtained with the idealized system as long as dimerization and hydrolysis are not responsible for consuming a substantial fraction of 2-MeImpG.

Results are presented from a detailed study of the P-N bond hydrolysis in guanosine 5-prime-monophosphate 2-methylimidazolide (2-MeImpG) and in guanosine 5-prime-imidazolide (ImpG) in the presence of 0-0.50 M Mg(2+). Pseudo-first-order rate constants of these compounds were obtained as a function of Mg(2+) concentration, for pH values between 6 and 10 and 37 C. It was found that Mg(2+) catalysis was most effective at pH 10, where a 15-fold increase in hydrolysis was achieved in 0.02 M Mg; at 0.2 M, a 115-fold increase was observed. Implication of these results for the mechanism of template-directed oligomerization is discussed.

We have undertaken a complete kinetic analysis of the template-directed oligoguanylate synthesis originated in Orgel's laboratory (Inoue and Orgel, 1982). The reaction of guanosine 5'-phospho-2-methylimidazolide, 2-MelmpG, with ribooligoguanylates all 3'-5' linked, designated n3 with n = 7-12, was studied in the presence/absence of the complementary template polycytidylic acid, poly(C). Conditions were chosen where poly(C) and 2-MelmpG are in large excess over the oligoguanylate. In the absence of the template at 37 degrees C the reaction leads to three isomeric oligomers that are elongated by one monomer unit. They are the 3'-5' linked, (n + 1)3, the 2'-5' linked, (n + 1)2, and the pyrophosphate product, (n + 1)p, formed in an approximate ratio 1:2:5. In the presence of the template the reaction is 20-fold faster and yields products n + 1, n + 2, n + 3 etc. as long as 2-MelmpG is available. Most importantly the formation of the natural, 3'-5' linked isomer, is enhanced selectively by 140-fold at 37 degrees C. Qualitative observations allow the conclusion that this enhancement is temperature dependent and increases with decreasing temperature. For example, at 1 degree C only the 3'-5' linked isomers were detected. Initial rates for the disappearance of the n3 oligoguanylate were determined at 1, 23, and 37 degrees C. It was found that the pseudo-first order rate constant for oligoguanylate elongation was linearly proportional to the 2-MelmpG concentration. This implies that the reaction complex poly(C).n3.2-MelmpG does not accumulate under the reaction conditions, a conclusion which is also supported by infrared data (Miles and Frazier, 1982). The implication of the above results with respect to chemical evolution is that lower temperatures, i.e., close to freezing, enhance the regioselectivity of these template-directed reactions and that one way to improve replication models may be sought in finding conditions that favor stable reaction complexes.

We have studied the hydrolysis of guanosine 5'-phospho-2-methylimidazolide, 2-MeImpG, in aqueous buffered solutions of various pH's at 75 degrees C and 37 degrees C. At 75 degrees C and pHor = 1.0, two kinetic processes were observed spectrophotometrically: the first and more rapid one is attributed to the hydrolysis of the phosphoimidazolide P-N bond; the second and much slower one, to the cleavage of the glycosidic bond. At 37 degrees C, pH 2.0, the spectrophotometrically determined rate constant of P-N bond hydrolysis was confirmed by using high pressure liquid chromatography, HPLC. With the latter technique it was possible to separate reactants and products and also to extend the pH-rate profile into the neutral region where rates are slower and, therefore, difficult to measure spectrophotometrically. The pH-rate profiles at both temperatures exhibit similar behavior. At pH2 the pseudo-first-order rate constant increases with decreasing pH; in the region 2pH7 there is a plateau followed by a decrease for pH7. These data are consistent with a reactivity order zwitterionanion for P-N bond hydrolysis. It is noteworthy that P-N bond hydrolysis in phosphoimidazolides is very slow compared to other phosphoramidates. This may be one of the reasons why this compound showed extraordinary ability in forming long oligomers under template-directed conditions.

Results are presented from a detailed study of the P-N bond hydrolysis in guanosine 5-prime-monophosphate 2-methylimidazolide (2-MeImpG) and in guanosine 5-prime-imidazolide (ImpG) in the presence of 0-0.50 M Mg(2+). Pseudo-first-order rate constants of these compounds were obtained as a function of Mg(2+) concentration, for pH values between 6 and 10 and 37 C. It was found that Mg(2+) catalysis was most effective at pH 10, where a 15-fold increase in hydrolysis was achieved in 0.02 M Mg; at 0.2 M, a 115-fold increase was observed. Implication of these results for the mechanism of template-directed oligomerization is discussed.