Role of unsatisfied hydrogen bond acceptors in RNA energetics and specificity

Unusual Base Pair between Two 2-Thiouridines and Its Implication for Nonenzymatic RNA Copying

2-Thiouridine (s2U) is a nucleobase modification that confers enhanced efficiency and fidelity both on modern tRNA codon translation and on nonenzymatic and ribozyme-catalyzed RNA copying. We have discovered an unusual base pair between two 2-thiouridines that stabilizes an RNA duplex to a degree that is comparable to that of a native A:U base pair. High-resolution crystal structures indicate similar base-pairing geometry and stacking interactions in duplexes containing s2U:s2U compared to those with U:U pairs. Notably, the C═O···H-N hydrogen bond in the U:U pair is replaced with a C═S···H-N hydrogen bond in the s2U:s2U base pair. The thermodynamic stability of the s2U:s2U base pair suggested that this self-pairing might lead to an increased error frequency during nonenzymatic RNA copying. However, competition experiments show that s2U:s2U base-pairing induces only a low level of misincorporation during nonenzymatic RNA template copying because the correct A:s2U base pair outcompetes the slightly weaker s2U:s2U base pair. In addition, even if an s2U is incorrectly incorporated, the addition of the next base is greatly hindered. This strong stalling effect would further increase the effective fidelity of nonenzymatic RNA copying with s2U. Our findings suggest that s2U may enhance the rate and extent of nonenzymatic copying with only a minimal cost in fidelity.

A unified computational view of DNA duplex, triplex, quadruplex and their donor-acceptor interactions

DNA can assume various structures as a result of interactions at atomic and molecular levels (e.g., hydrogen bonds, π–π stacking interactions, and electrostatic potentials), so understanding of the consequences of these interactions could guide development of ways to produce elaborate programmable DNA for applications in bio- and nanotechnology. We conducted advanced ab initio calculations to investigate nucleobase model structures by componentizing their donor-acceptor interactions. By unifying computational conditions, we compared the independent interactions of DNA duplexes, triplexes, and quadruplexes, which led us to evaluate a stability trend among Watson–Crick and Hoogsteen base pairing, stacking, and even ion binding. For a realistic solution-like environment, the influence of water molecules was carefully considered, and the potassium-ion preference of G-quadruplex was first analyzed at an ab initio level by considering both base-base and ion-water interactions. We devised new structure factors including hydrogen bond length, glycosidic vector angle, and twist angle, which were highly effective for comparison between computationally-predicted and experimentally-determined structures; we clarified the function of phosphate backbone during nucleobase ordering. The simulated tendency of net interaction energies agreed well with that of real world, and this agreement validates the potential of ab initio study to guide programming of complicated DNA constructs.

The Emergence of RNA from the Heterogeneous Products of Prebiotic Nucleotide Synthesis

Recent advances in prebiotic chemistry are beginning to outline plausible pathways for the synthesis of the canonical ribonucleotides and their assembly into oligoribonucleotides. However, these reaction pathways suggest that many noncanonical nucleotides are likely to have been generated alongside the standard ribonucleotides. Thus, the oligomerization of prebiotically synthesized nucleotides is likely to have led to a highly heterogeneous collection of oligonucleotides comprised of a wide range of types of nucleotides connected by a variety of backbone linkages. How then did relatively homogeneous RNA emerge from this primordial heterogeneity? Here we focus on nonenzymatic template-directed primer extension as a process that would have strongly enriched for homogeneous RNA over the course of multiple cycles of replication. We review the effects on copying the kinetics of nucleotides with altered nucleobase and sugar moieties, when they are present as activated monomers and when they are incorporated into primer and template oligonucleotides. We also discuss three variations in backbone connectivity, all of which are nonheritable and regenerate native RNA upon being copied. The kinetic superiority of RNA synthesis suggests that nonenzymatic copying served as a chemical selection mechanism that allowed relatively homogeneous RNA to emerge from a complex mixture of prebiotically synthesized nucleotides and oligonucleotides.

Terpene Chain Length Affects the Base Pairing Discrimination of S-geranyl-2-thiouridine in RNA Duplex

Geranylation is a hydrophobic modification discovered in several bacteria tRNAs with the function of promoting codon bias during translation. However, why nature selects this C10-geranyl group remains a question. We conduct synthesis, UV-thermal denaturation, and molecular simulation studies in RNA duplexes and reveal possible reasons behind this natural selection. Among methyl-(C1), dimethylallyl-(C5), geranyl-(C10), and farnesyl-(C15) modified 2-thiouridines, only geranyl-group promotes U:G over U:A pair. Molecular simulation shows all the modified terpene groups point to the minor groove of RNA duplexes. The discrimination between U:G and U:A pairs derives from the difference in hydrogen bonding and interactions of the chain with the hydrophobic area in the minor groove. Geranyl group has perfect length to discriminate U:G and U:A pairs, whereas the others are either too long or too short to achieve the same behavior. This work indicates that geranyl group cannot be replaced by other terpene groups in promoting codon-specificity.

•

Geranyl group on ges²U has perfect length to enhance binding to G over A in RNA

•

Other ges²U analogues has no pairing preference between G or A

•

MD shows the molecular dynamic behavior of different terpene chains in RNA duplex

Incorporating a Thiophosphate Modification into a Common RNA Tetraloop Motif Causes an Unanticipated Stability Boost

GNRA (N = A, C, G, or U; R = A or G) tetraloops are common RNA secondary structural motifs and feature a phosphate stacked atop a nucleobase. The rRNA sarcin/ricin loop (SRL) is capped by GApGA, and the phosphate p stacks on G. We recently found that regiospecific incorporation of a single dithiophosphate (PS2) but not a monothiophosphate (PSO) instead of phosphate in the backbone of RNA aptamers dramatically increases the binding affinity for their targets. In the RNA:thrombin complex, the key contribution to the 1000-fold tighter binding stems from an edge-on contact between PS2 and a phenylalanine ring. Here we investigated the consequences of replacing the SRL phosphate engaged in a face-on interaction with guanine with either PS2 or PSO for stability. We found that PS2···G and Rp-PSO···G contacts stabilize modified SRLs compared to the parent loop to unexpected levels: up to 6.3 °C in melting temperature T_m and -4.7 kcal/mol in ΔΔG°. Crystal structures demonstrate that the vertical distance to guanine for the closest sulfur is just 0.05 Å longer on average compared to that of oxygen despite the larger van der Waals radius of the former (1.80 Å for S vs 1.52 Å for O). The higher stability is enthalpy-based, and the negative charge as assessed by a neutral methylphosphonate modification plays only a minor role. Quantum mechanical/molecular mechanical calculations are supportive of favorable dispersion attraction interactions by sulfur making the dominant contribution. A stacking interaction between phosphate and guanine (SRL) or uracil (U-turn) is also found in newly classified RNA tetraloop families besides GNRA.

A systematic study on the influence of thermodynamic asymmetry of 5'-ends of siRNA duplexes in relation to their silencing potency

siRNA molecules possess high potential as molecular tools and can be used as effective therapeutics in humans. One of the key steps in the action of these molecules is the choice of antisense strand by the RNA-induced silencing complex (RISC). To explain this process, we verified the theory which states that antisense strand selection is based on the thermodynamically less stable 5′ end of siRNA. Based on the studies presented herein, we observed that for the tested siRNA duplexes, the difference in the thermodynamic stability of the terminal, penultimate and pre-penultimate pairs in the duplex siRNA is not the dominant factor in antisense strand selection. We found that both strands in each tested siRNA molecule are used as an antisense strand. The introduction of modified nucleotides, whose impact on the thermodynamic stability of siRNA duplexes was studied, results in changes in antisense strand selection by the RISC complex. The presence of a modified residue often caused predominant selection of only one antisense strand which is at variance with the theory of siRNA strand bias.

Inosine, but none of the 8-oxo-purines, is a plausible component of a primordial version of RNA

The RNA world hypothesis assumes the abiotic synthesis of nucleotides, as well as their participation in nonenzymatic RNA replication. Whereas prebiotic syntheses of the canonical purine nucleotides remain inefficient, a prebiotically plausible route to the 8-oxo-purines has been reported. Although these noncanonical purine nucleotides are known to engage in non-Watson–Crick pairing with their canonical purine counterparts, their behavior in nonenzymatic RNA copying has not been evaluated. Our study indicates that none of the 8-oxo-purines behaves as a suitable substrate for nonenzymatic RNA copying. However, inosine turns out to exhibit reasonable rates and fidelities in RNA copying reactions. We propose that inosine could have served as a surrogate for guanosine in the early emergence of life.

The emergence of primordial RNA-based life would have required the abiotic synthesis of nucleotides, and their participation in nonenzymatic RNA replication. Although considerable progress has been made toward potentially prebiotic syntheses of the pyrimidine nucleotides (C and U) and their 2-thio variants, efficient routes to the canonical purine nucleotides (A and G) remain elusive. Reported syntheses are low yielding and generate a large number of undesired side products. Recently, a potentially prebiotic pathway to 8-oxo-adenosine and 8-oxo-inosine has been demonstrated, raising the question of the suitability of the 8-oxo-purines as substrates for prebiotic RNA replication. Here we show that the 8-oxo-purine nucleotides are poor substrates for nonenzymatic RNA primer extension, both as activated monomers and when present in the template strand; their presence at the end of a primer also strongly reduces the rate and fidelity of primer extension. To provide a proper comparison with 8-oxo-inosine, we also examined primer extension reactions with inosine, and found that inosine exhibits surprisingly rapid and accurate nonenzymatic RNA copying. We propose that inosine, which can be derived from adenosine by deamination, could have acted as a surrogate for G in the earliest stages of the emergence of life.

An intricate balance of hydrogen bonding, ion atmosphere and dynamics facilitates a seamless uracil to cytosine substitution in the U-turn of the neomycin-sensing riboswitch

The neomycin sensing riboswitch is the smallest biologically functional RNA riboswitch, forming a hairpin capped with a U-turn loop-a well-known RNA motif containing a conserved uracil. It was shown previously that a U→C substitution of the eponymous conserved uracil does not alter the riboswitch structure due to C protonation at N3. Furthermore, cytosine is evolutionary permitted to replace uracil in other U-turns. Here, we use molecular dynamics simulations to study the molecular basis of this substitution in the neomycin sensing riboswitch and show that a structure-stabilizing monovalent cation-binding site in the wild-type RNA is the main reason for its negligible structural effect. We then use NMR spectroscopy to confirm the existence of this cation-binding site and to demonstrate its effects on RNA stability. Lastly, using quantum chemical calculations, we show that the cation-binding site is altering the electronic environment of the wild-type U-turn so that it is more similar to the cytosine mutant. The study reveals an amazingly complex and delicate interplay between various energy contributions shaping up the 3D structure and evolution of nucleic acids.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Thiolated uridine substrates and templates improve the rate and fidelity of ribozyme-catalyzed RNA copying

Ribozyme-catalyzed RNA polymerization is inefficient and error prone. Here we demonstrate that two alternative bases, 2-thio-uridine (s(2)U) and 2-thio-ribo-thymidine (s(2)T), improve the rate and fidelity of ribozyme catalyzed nucleotide addition as NTP substrates and as template bases. We also demonstrate the functionality of s(2)U and s(2)T-containing ribozymes.

Base pairing and structural insights into the 5-formylcytosine in RNA duplex

5-Formylcytidine (f⁵C), a previously discovered natural nucleotide in the mitochondrial tRNA of many species including human, has been recently detected as the oxidative product of 5-methylcytidine (m⁵C) through 5-hydroxymethylcytidine (hm⁵C) in total RNA of mammalian cells. The discovery indicated that these cytosine derivatives in RNA might also play important epigenetic roles similar as in DNA, which has been intensively investigated in the past few years. In this paper, we studied the base pairing specificity of f⁵C in different RNA duplex contexts. We found that the 5-formyl group could increase duplex thermal stability and enhance base pairing specificity. We present three high-resolution crystal structures of an octamer RNA duplex [5′-GUA(f⁵C)GUAC-3′]₂ that have been solved under three crystallization conditions with different buffers and pH values. Our results showed that the 5-formyl group is located in the same plane as the cytosine base and forms an intra-residue hydrogen bond with the amino group in the N4 position. In addition, this modification increases the base stacking between the f⁵C and the neighboring bases while not causing significant global and local structure perturbations. This work provides insights into the effects of 5-formylcytosine on RNA duplex.

Polymerase recognition of 2-thio-iso-guanine·5-methyl-4-pyrimidinone (iGs·P)--A new DD/AA base pair

Polymerase specificity is reported for a previously unknown base pair with a non-standard DD/AA hydrogen bonding pattern: 2-thio-iso-guanine·5-methyl-4-pyrimidinone. Our findings suggest that atomic substitution may provide a solution for low fidelity previously associated with enzymatic copying of iso-guanine.

Watson-Crick hydrogen bonding of unlocked nucleic acids

We herein describe the synthesis of two new unlocked nucleic acid building blocks containing hypoxanthine and 2,6-diaminopurine as nucleobase moieties and their incorporation into oligonucleotides. The modified oligonucleotides were used to examine the thermodynamic properties of UNA against unmodified oligonucleotides and the resulting thermodynamic data support that the hydrogen bonding face of UNA is Watson-Crick like.

Thermodynamic insights into 2-thiouridine-enhanced RNA hybridization

Nucleobase modifications dramatically alter nucleic acid structure and thermodynamics. 2-thiouridine (s(2)U) is a modified nucleobase found in tRNAs and known to stabilize U:A base pairs and destabilize U:G wobble pairs. The recently reported crystal structures of s(2)U-containing RNA duplexes do not entirely explain the mechanisms responsible for the stabilizing effect of s(2)U or whether this effect is entropic or enthalpic in origin. We present here thermodynamic evaluations of duplex formation using ITC and UV thermal denaturation with RNA duplexes containing internal s(2)U:A and s(2)U:U pairs and their native counterparts. These results indicate that s(2)U stabilizes both duplexes. The stabilizing effect is entropic in origin and likely results from the s(2)U-induced preorganization of the single-stranded RNA prior to hybridization. The same preorganizing effect is likely responsible for structurally resolving the s(2)U:U pair-containing duplex into a single conformation with a well-defined H-bond geometry. We also evaluate the effect of s(2)U on single strand conformation using UV- and CD-monitored thermal denaturation and on nucleoside conformation using (1)H NMR spectroscopy, MD and umbrella sampling. These results provide insights into the effects that nucleobase modification has on RNA structure and thermodynamics and inform efforts toward improving both ribozyme-catalyzed and nonenzymatic RNA copying.

The role of N7 protonation of guanine in determining the structure, stability and function of RNA base pairs

Ab initio computations and bioinformatics studies reveal that stabilization of some important RNA structural motifs might involve N7 protonation of guanine.

Crystal structure studies of RNA duplexes containing s(2)U:A and s(2)U:U base pairs

Structural studies of modified nucleobases in RNA duplexes are critical for developing a full understanding of the stability and specificity of RNA base pairing. 2-Thio-uridine (s²U) is a modified nucleobase found in certain tRNAs. Thermodynamic studies have evaluated the effects of s²U on base pairing in RNA, where it has been shown to stabilize U:A pairs and destabilize U:G wobble pairs. Surprisingly, no high-resolution crystal structures of s²U-containing RNA duplexes have yet been reported. We present here two high-resolution crystal structures of heptamer RNA duplexes (5′-uagcs²Ucc-3′ paired with 3′-aucgAgg-5′ and with 3′-aucgUgg-5′) containing s²U:A and s²U:U pairs, respectively. For comparison, we also present the structures of their native counterparts solved under identical conditions. We found that replacing O2 with S2 stabilizes the U:A base pair without any detectable structural perturbation. In contrast, an s²U:U base pair is strongly stabilized in one specific U:U pairing conformation out of four observed for the native U:U base pair. This s²U:U stabilization appears to be due at least in part to an unexpected sulfur-mediated hydrogen bond. This work provides additional insights into the effects of 2-thio-uridine on RNA base pairing.

Insights into gas-phase structural conformers of hydrated rubidium and cesium cations, M(+)(H2O)(n)Ar (M = Rb, Cs; n = 3-5), using infrared photodissociation spectroscopy

Infrared photodissociation (IRPD) spectra of M(+)(H2O)nAr (M = Rb, Cs; n = 3-5) with simultaneous monitoring of [Ar] and [Ar+H2O] fragmentation channels are reported. The comparison between the spectral features in the two channels and corresponding energy analysis provide spectral assignments of the stable structural conformers and insight into the competition between ion-water electrostatic and water-water hydrogen bonding interactions. Results show that as the level of hydration increases, the water-water interaction exhibits the tendency to dominate over the ion-water interaction. Cyclic water tetramer and water pentamer substructures appear in Cs(+)(H2O)4Ar and Cs(+)(H2O)5Ar systems, respectively. However, cyclic water tetramer and pentamer structures were not observed for Rb(+)(H2O)4Ar and Rb(+)(H2O)5Ar systems, respectively, due to the stronger influence of the rubidium ion-water electrostatic interaction. The energy analysis, including the available internal energy and the IR photon energy, helped provide an experimental estimate of water binding energies.

Incorporation of thio-pseudoisocytosine into triplex-forming peptide nucleic acids for enhanced recognition of RNA duplexes

Peptide nucleic acids (PNAs) have been developed for applications in biotechnology and therapeutics. There is great potential in the development of chemically modified PNAs or other triplex-forming ligands that selectively bind to RNA duplexes, but not single-stranded regions, at near-physiological conditions. Here, we report on a convenient synthesis route to a modified PNA monomer, thio-pseudoisocytosine (L), and binding studies of PNAs incorporating the monomer L. Thermal melting and gel electrophoresis studies reveal that L-incorporated 8-mer PNAs have superior affinity and specificity in recognizing the duplex region of a model RNA hairpin to form a pyrimidine motif major-groove RNA₂–PNA triplex, without appreciable binding to single-stranded regions to form an RNA–PNA duplex or, via strand invasion, forming an RNA–PNA₂ triplex at near-physiological buffer condition. In addition, an L-incorporated 8-mer PNA shows essentially no binding to single-stranded or double-stranded DNA. Furthermore, an L-modified 6-mer PNA, but not pseudoisocytosine (J) modified or unmodified PNA, binds to the HIV-1 programmed −1 ribosomal frameshift stimulatory RNA hairpin at near-physiological buffer conditions. The stabilization of an RNA₂–PNA triplex by L modification is facilitated by enhanced van der Waals contacts, base stacking, hydrogen bonding and reduced dehydration energy. The destabilization of RNA–PNA and DNA–PNA duplexes by L modification is due to the steric clash and loss of two hydrogen bonds in a Watson–Crick-like G–L pair. An RNA₂–PNA triplex is significantly more stable than a DNA₂–PNA triplex, probably because the RNA duplex major groove provides geometry compatibility and favorable backbone–backbone interactions with PNA. Thus, L-modified triplex-forming PNAs may be utilized for sequence-specifically targeting duplex regions in RNAs for biological and therapeutic applications.

How to understand quantum chemical computations on DNA and RNA systems? A practical guide for non-specialists

In this review primarily written for non-experts we explain basic methodological aspects and interpretation of modern quantum chemical (QM) computations applied to nucleic acids. We introduce current reference QM computations on small model systems consisting of dozens of atoms. Then we comment on recent advance of fast and accurate dispersion-corrected density functional theory methods, which will allow computations of small but complete nucleic acids building blocks in the near future. The qualitative difference between QM and molecular mechanics (MM, force field) computations is discussed. We also explain relation of QM and molecular simulation computations to experiments.

On the Origin of the Canonical Nucleobases: An Assessment of Selection Pressures across Chemical and Early Biological Evolution

The native bases of RNA and DNA are prominent examples of the narrow selection of organic molecules upon which life is based. How did nature "decide" upon these specific heterocycles? Evidence suggests that many types of heterocycles could have been present on the early Earth. It is therefore likely that the contemporary composition of nucleobases is a result of multiple selection pressures that operated during early chemical and biological evolution. The persistence of the fittest heterocycles in the prebiotic environment towards, for example, hydrolytic and photochemical assaults, may have given some nucleobases a selective advantage for incorporation into the first informational polymers. The prebiotic formation of polymeric nucleic acids employing the native bases remains, however, a challenging problem to reconcile. Hypotheses have proposed that the emerging RNA world may have included many types of nucleobases. This is supported by the extensive utilization of non-canonical nucleobases in extant RNA and the resemblance of many of the modified bases to heterocycles generated in simulated prebiotic chemistry experiments. Selection pressures in the RNA world could have therefore narrowed the composition of the nucleic acid bases. Two such selection pressures may have been related to genetic fidelity and duplex stability. Considering these possible selection criteria, the native bases along with other related heterocycles seem to exhibit a certain level of fitness. We end by discussing the strength of the N-glycosidic bond as a potential fitness parameter in the early DNA world, which may have played a part in the refinement of the alphabetic bases.

Recognition of RNA duplexes by chemically modified triplex-forming oligonucleotides

Triplex is emerging as an important RNA tertiary structure motif, in which consecutive non-canonical base pairs form between a duplex and a third strand. RNA duplex region is also often functionally important site for protein binding. Thus, triplex-forming oligonucleotides (TFOs) may be developed to regulate various biological functions involving RNA, such as viral ribosomal frameshifting and reverse transcription. How chemical modification in TFOs affects RNA triplex stability, however, is not well understood. Here, we incorporated locked nucleic acid, 2-thio U- and 2′-O methyl-modified residues in a series of all pyrimidine RNA TFOs, and we studied the binding to two RNA hairpin structures. The 12-base-triple major-groove pyrimidine–purine–pyrimidine triplex structures form between the duplex regions of RNA/DNA hairpins and the complementary RNA TFOs. Ultraviolet-absorbance-detected thermal melting studies reveal that the locked nucleic acid and 2-thio U modifications in TFOs strongly enhance triplex formation with both parental RNA and DNA duplex regions. In addition, we found that incorporation of 2′-O methyl-modified residues in a TFO destabilizes and stabilizes triplex formation with RNA and DNA duplex regions, respectively. The (de)stabilization of RNA triplex formation may be facilitated through modulation of van der Waals contact, base stacking, hydrogen bonding, backbone pre-organization, geometric compatibility and/or dehydration energy. Better understanding of the molecular determinants of RNA triplex structure stability lays the foundation for designing and discovering novel sequence-specific duplex-binding ligands as diagnostic and therapeutic agents targeting RNA.

Hydration sites of unpaired RNA bases: a statistical analysis of the PDB structures

BACKGROUND

Hydration is crucial for RNA structure and function. X-ray crystallography is the most commonly used method to determine RNA structures and hydration and, therefore, statistical surveys are based on crystallographic results, the number of which is quickly increasing.

RESULTS

A statistical analysis of the water molecule distribution in high-resolution X-ray structures of unpaired RNA nucleotides showed that: different bases have the same penchant to be surrounded by water molecules; clusters of water molecules indicate possible hydration sites, which, in some cases, match those of the major and minor grooves of RNA and DNA double helices; complex hydrogen bond networks characterize the solvation of the nucleotides, resulting in a significant rigidity of the base and its surrounding water molecules. Interestingly, the hydration sites around unpaired RNA bases do not match, in general, the positions that are occupied by the second nucleotide when the base-pair is formed.

CONCLUSIONS

The hydration sites around unpaired RNA bases were found. They do not replicate the atom positions of complementary bases in the Watson-Crick pairs.

NMR structure of a 4 x 4 nucleotide RNA internal loop from an R2 retrotransposon: identification of a three purine-purine sheared pair motif and comparison to MC-SYM predictions

The NMR solution structure is reported of a duplex, 5'GUGAAGCCCGU/3'UCACAGGAGGC, containing a 4 × 4 nucleotide internal loop from an R2 retrotransposon RNA. The loop contains three sheared purine-purine pairs and reveals a structural element found in other RNAs, which we refer to as the 3RRs motif. Optical melting measurements of the thermodynamics of the duplex indicate that the internal loop is 1.6 kcal/mol more stable at 37°C than predicted. The results identify the 3RRs motif as a common structural element that can facilitate prediction of 3D structure. Known examples include internal loops having the pairings: 5'GAA/3'AGG, 5'GAG/3'AGG, 5'GAA/3'AAG, and 5'AAG/3'AGG. The structural information is compared with predictions made with the MC-Sym program.

Tightening of active site interactions en route to the transition state revealed by single-atom substitution in the guanosine-binding site of the Tetrahymena group I ribozyme

Protein enzymes establish intricate networks of interactions to bind and position substrates and catalytic groups within active sites, enabling stabilization of the chemical transition state. Crystal structures of several RNA enzymes also suggest extensive interaction networks, despite RNA's structural limitations, but there is little information on the functional and the energetic properties of these inferred networks. We used double mutant cycles and presteady-state kinetic analyses to probe the putative interaction between the exocyclic amino group of the guanosine nucleophile and the N7 atom of residue G264 of the Tetrahymena group I ribozyme. As expected, the results supported the presence of this interaction, but remarkably, the energetic penalty for introducing a CH group at the 7-position of residue G264 accumulates as the reaction proceeds toward the chemical transition state to a total of 6.2 kcal/mol. Functional tests of neighboring interactions revealed that the presence of the CH group compromises multiple contacts within the interaction network that encompass the reactive elements, apparently forcing the nucleophile to bind and attack from an altered, suboptimal orientation. The energetic consequences of this indirect disruption of neighboring interactions as the reaction proceeds demonstrate that linkage between binding interactions and catalysis hinges critically on the precise structural integrity of a network of interacting groups.

Stability of DNA containing a structural water mimic in an A-T rich sequence

We describe here the synthesis and properties of A-T rich DNA containing covalently bound water mimics located in the DNA minor groove.

Quantum chemical studies of nucleic acids: can we construct a bridge to the RNA structural biology and bioinformatics communities?

In this feature article, we provide a side-by-side introduction for two research fields: quantum chemical calculations of molecular interaction in nucleic acids and RNA structural bioinformatics. Our main aim is to demonstrate that these research areas, while largely separated in contemporary literature, have substantial potential to complement each other that could significantly contribute to our understanding of the exciting world of nucleic acids. We identify research questions amenable to the combined application of modern ab initio methods and bioinformatics analysis of experimental structures while also assessing the limitations of these approaches. The ultimate aim is to attain valuable physicochemical insights regarding the nature of the fundamental molecular interactions and how they shape RNA structures, dynamics, function, and evolution.