Hydration of the RNA duplex r(CGCAAAUUUGCG)2 determined by NMR

Molecular Origin of Distinct Hydration Dynamics in Double Helical DNA and RNA Sequences

Water molecules are essential to determine the structure of nucleic acids and mediate their interactions with other biomolecules. Here, we characterize the hydration dynamics of analogous DNA and RNA double helices with unprecedented resolution and elucidate the molecular origin of their differences: first, the localization of the slowest hydration water molecules─in the minor groove in DNA, next to phosphates in RNA─and second, the markedly distinct hydration dynamics of the two phosphate oxygen atoms OR and OS in RNA. Using our Extended Jump Model for water reorientation, we assess the relative importance of previously proposed factors, including the local topography, water bridges, and the presence of ions. We show that the slow hydration dynamics at RNA OR sites is not due to bridging water molecules but is caused by both the larger excluded volume and the stronger initial H-bond next to OR, due to the different phosphate orientations in A-form double helical RNA.

Inosine-Induced Base Pairing Diversity during Reverse Transcription

A-to-I editing catalyzed by adenosine deaminase acting on RNAs impacts numerous physiological and biochemical processes that are essential for cellular functions and is a big contributor to the infectivity of certain RNA viruses. The outcome of this deamination leads to changes in the eukaryotic transcriptome functionally resembling A-G transitions since inosine preferentially pairs with cytosine. Moreover, hyper-editing or multiple A to G transitions in clusters were detected in measles virus. Inosine modifications either directly on viral RNA or on cellular RNA can have antiviral or pro-viral repercussions. While many of the significant roles of inosine in cellular RNAs are well understood, the effects of hyper-editing of A to I on viral polymerase activity during RNA replication remain elusive. Moreover, biological strategies such as molecular cloning and RNA-seq for transcriptomic interrogation rely on RT-polymerase chain reaction with little to no emphasis placed on the first step, reverse transcription, which may reshape the sequencing results when hypermodification is present. In this study, we systematically explore the influence of inosine modification, varying the number and position of inosines, on decoding outcomes using three different reverse transcriptases (RTs) followed by standard Sanger sequencing. We find that inosine alone or in clusters can differentially affect the RT activity. To gain structural insights into the accommodation of inosine in the polymerase site of HIV-1 reverse transcriptase (HIV-1-RT) and how this structural context affects the base pairing rules for inosine, we performed molecular dynamics simulations of the HIV-1-RT. The simulations highlight the importance of the protein-nucleotide interaction as a critical factor in deciphering the base pairing behavior of inosine clusters. This effort sets the groundwork for decrypting the physiological significance of inosine and linking the fidelity of reverse transcriptase and the possible diverse transcription outcomes of cellular RNAs and/or viral RNAs where hyper-edited inosines are present in the transcripts.

Characterizing micro-to-millisecond chemical exchange in nucleic acids using off-resonance R1ρ relaxation dispersion

This review describes off-resonance R1ρ relaxation dispersion NMR methods for characterizing microsecond-to-millisecond chemical exchange in uniformly 13C/15N labeled nucleic acids in solution. The review opens with a historical account of key developments that formed the basis for modern R1ρ techniques used to study chemical exchange in biomolecules. A vector model is then used to describe the R1ρ relaxation dispersion experiment, and how the exchange contribution to relaxation varies with the amplitude and frequency offset of an applied spin-locking field, as well as the population, exchange rate, and differences in chemical shifts of two exchanging species. Mathematical treatment of chemical exchange based on the Bloch-McConnell equations is then presented and used to examine relaxation dispersion profiles for more complex exchange scenarios including three-state exchange. Pulse sequences that employ selective Hartmann-Hahn cross-polarization transfers to excite individual 13C or 15N spins are then described for measuring off-resonance R1ρ(13C) and R1ρ(15N) in uniformly 13C/15N labeled DNA and RNA samples prepared using commercially available 13C/15N labeled nucleotide triphosphates. Approaches for analyzing R1ρ data measured at a single static magnetic field to extract a full set of exchange parameters are then presented that rely on numerical integration of the Bloch-McConnell equations or the use of algebraic expressions. Methods for determining structures of nucleic acid excited states are then reviewed that rely on mutations and chemical modifications to bias conformational equilibria, as well as structure-based approaches to calculate chemical shifts. Applications of the methodology to the study of DNA and RNA conformational dynamics are reviewed and the biological significance of the exchange processes is briefly discussed.

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.

Linker dependent intercalation of bisbenzimidazole-aminosugars in an RNA duplex; selectivity in RNA vs. DNA binding

Neomycin and Hoechst 33258 are two well-known nucleic acid binders that interact with RNA and DNA duplexes with high affinities respectively. In this manuscript, we report that covalent attachment of bisbenzimidazole unit derived from Hoechst 33258 to neomycin leads to intercalative binding of the bisbenzimidazole unit (oriented at 64-74° with respected to the RNA helical axis) in a linker length dependent manner. The dual binding and intercalation of conjugates were supported by thermal denaturation, CD, LD and UV-Vis absorption experiments. These studies highlight the importance of linker length in dual recognition by conjugates, for effective RNA recognition, which can lead to novel ways of recognizing RNA structures. Additionally, the ligand library screens also identify DNA and RNA selective compounds, with compound 9, containing a long linker, showing a 20.3°C change in RNA duplex Tm with only a 13.0°C change in Tm for the corresponding DNA duplex. Significantly, the shorter linker in compound 3 shows almost the reverse trend, a 23.8°C change in DNA Tm, with only a 9.1°C change in Tm for the corresponding RNA duplex.

Uncovering the thermodynamics of monomer binding for RNA replication

The nonenzymatic replication of primordial RNA is thought to have been a critical step in the origin of life. However, despite decades of effort, the poor rate and fidelity of model template copying reactions have thus far prevented an experimental demonstration of nonenzymatic RNA replication. The overall rate and fidelity of template copying depend, in part, on the affinity of free ribonucleotides to the RNA primer-template complex. We have now used (1)H NMR spectroscopy to directly measure the thermodynamic association constants, Kas, of the standard ribonucleotide monophosphates (rNMPs) to native RNA primer-template complexes. The binding affinities of rNMPs to duplexes with a complementary single-nucleotide overhang follow the order C > G > A > U. Notably, these monomers bind more strongly to RNA primer-template complexes than to the analogous DNA complexes. The relative binding affinities of the rNMPs for complementary RNA primer-template complexes are in good quantitative agreement with the predictions of a nearest-neighbor analysis. With respect to G:U wobble base-pairing, we find that the binding of rGMP to a primer-template complex with a 5'-U overhang is approximately 10-fold weaker than to the complementary 5'-C overhang. We also find that the binding of rGMP is only about 2-fold weaker than the binding of rAMP to 5'-U, consistent with the poor fidelity observed in the nonenzymatic copying of U residues in RNA templates. The accurate Ka measurements for ribonucleotides obtained in this study will be useful for designing higher fidelity, more effective RNA replication systems.

Impact of 2'-hydroxyl sampling on the conformational properties of RNA: update of the CHARMM all-atom additive force field for RNA

Here, we present an update of the CHARMM27 all-atom additive force field for nucleic acids that improves the treatment of RNA molecules. The original CHARMM27 force field parameters exhibit enhanced Watson-Crick base pair opening which is not consistent with experiment, whereas analysis of molecular dynamics (MD) simulations show the 2'-hydroxyl moiety to almost exclusively sample the O3' orientation. Quantum mechanical (QM) studies of RNA related model compounds indicate the energy minimum associated with the O3' orientation to be too favorable, consistent with the MD results. Optimization of the dihedral parameters dictating the energy of the 2'-hydroxyl proton targeting the QM data yielded several parameter sets, which sample both the base and O3' orientations of the 2'-hydroxyl to varying degrees. Selection of the final dihedral parameters was based on reproduction of hydration behavior as related to a survey of crystallographic data and better agreement with experimental NMR J-coupling values. Application of the model, designated CHARMM36, to a collection of canonical and noncanonical RNA molecules reveals overall improved agreement with a range of experimental observables as compared to CHARMM27. The results also indicate the sensitivity of the conformational heterogeneity of RNA to the orientation of the 2'-hydroxyl moiety and support a model whereby the 2'-hydroxyl can enhance the probability of conformational transitions in RNA.

Molecular dynamics of potential rRNA binders: single-stranded nucleic acids and some analogues

By hindering or "silencing" protein translation in vivo, antisense nucleic acid analogues that hybridize to bacterial rRNA could serve as a promising class of antibacterial compounds. Thus, we performed a comparative analysis of the dynamical properties of modified oligonucleotides based upon a sequence (5')r(UGUUACGACU)(3') that is complementary to bacterial ribosomal A-site RNA. In particular, 25 ns explicit solvent molecular dynamics simulations were computed for the following six single-stranded decamers: (1) the above RNA in unmodified form; (2) the 2'-O-methyl-modified RNA; (3) peptide nucleic acid (PNA) analogues of the above sequence, containing either (a) T or (b) U; and (4) two serine-substituted PNAs. Our results show that 2'-O-methylation attenuates RNA backbone dynamics, thereby preventing interconversion between stacked and unstacked conformations. The PNA analogue is rendered less flexible by replacing uracil with thymine; in addition, we found that derivatizing the PNA backbone with serine leads to enhanced base-stacking interactions. Consistent with known solubility properties of these classes of molecules, both RNAs exhibited greater localization of water molecules than did PNA. In terms of counterions, the initially helical conformation of the 2'-O-methyl RNA exhibits the highest Na(+) density among all the simulated decamers, while Na(+) build-up was most negligible for the neutral PNA systems. Further studies of the conformational and physicochemical properties of such modified single-stranded oligomers may facilitate better design of nucleic acid analogues, particularly those capable of serving as specific, high-affinity ribosomal A-site binders.

Hydration and conformational transitions in DNA, RNA, and mixed DNA-RNA triplexes studied by gravimetry and FTIR spectroscopy

We have studied by gravimetric measurements and FTIR spectroscopy the hydration of duplexes and triplexes formed by combinations of dA(n), dT(n), rA(n), and rU(n) strands. Results obtained on hydrated films show important differences in their hydration and in the structural transitions which can be induced by varying the water content of the samples. The number of water molecules per nucleotide (w/n) measured at high relative humidity (98% R.H.) is found to be 21 for dA(n).dT(n) and 15 for rA(n).rU(n). Addition of a third rU(n) strand does not change the number of water molecules per nucleotide: w/n=21 for rU(n)*dA(n).dT(n) and w/n=15 for rU(n)*rA(n).rU(n). On the contrary, the addition of a third dT(n) strand changes the water content but in a different way, depending whether the duplex is DNA or RNA. Thus, a loss of four water molecules per nucleotide is measured for dT(n)*dA(n).dT(n) while an increase of two water molecules per nucleotide is observed for dT(n)*rA(n).rU(n). The final hydration is the same for both triplexes (w/n=17). The desorption profiles obtained by gravimetry and FTIR spectroscopy are similar for the rA(n).rU(n) duplex and the rU(n)*rA(n).rU(n) triplex. On the contrary, the desorption profiles of the dA(n).dT(n) duplex and the triplexes formed with it (rU(n)*dA(n).dT(n) and dT(n)*dA(n).dT(n)) are different from each other. This is correlated with conformational transitions induced by varying the hydration content of the different structures, as shown by FTIR spectroscopy. Modifications of the phosphate group hydration and of the sugar conformation (S to N type repuckering) induced by decrease of the water content are observed in the case of triplexes formed on the dA(n).dT(n) duplex.

Influence of the 2'-hydroxyl group conformation on the stability of A-form helices in RNA

The 2'-hydroxyl group plays fundamental roles in both the structure and the function of RNA, and is the major determinant of the conformational and thermodynamic differences between RNA and DNA. Here, we report a conformational analysis of 2'-OH groups of the HIV-2 TAR RNA by means of NMR scalar coupling measurements in solution. Our analysis supports the existence of a network of water molecules spanning the minor groove of an RNA A-form helix, as has been suggested on the basis of a high-resolution X-ray study of an RNA duplex. The 2'-OH protons of the lower stem nucleotides of the TAR RNA project either towards the O3' or towards the base, where the 2'-OH group can favorably participate in H-bonding interactions with a water molecule situated in the nucleotide base plane. We observe that the k(ex) rate of the 2'-OH proton with the bulk solvent anti-correlates with the base-pair stability, confirming the involvement of the 2'-OH group in a collective network of H-bonds, which requires the presence of canonical helical secondary structure. The methodology and conformational analysis presented here are broadly applicable and facilitate future studies aimed to correlate the conformation of the 2'-OH group with both the structure and the function of RNA and RNA-ligand complexes.

Assignment and NOE analysis of 2'-hydroxyl protons in RNA: implications for stabilization of RNA A-form duplexes

The ribose 2'-OH hydroxyl group distinguishes RNA from DNA. The 2'-OH hydroxyl protons are responsible for differences in conformation, hydration, and thermodynamic stability of RNA and DNA oligonucleotides. Additionally, the 2'-OH group plays a central role in RNA catalysis. This important group lies in the shallow groove of RNA, where it is involved in a network of hydrogen bonds with water molecules stabilizing RNA A-form duplexes. Structural and dynamical information on 2'-OH hydroxyl protons is essential to understand their respective roles. Here we report the 2'-OH hydroxyl proton assignments for a 30mer RNA, the HIV-2 transactivation region, in water using solution NMR techniques. We provide structural information on 2'-OH hydroxyl groups in the form of orientational preferences contradicting the paradigm that the 2'-OH hydroxyl typically points away from the ribose H1' proton.

The solution structure of a DNA*RNA duplex containing 5-propynyl U and C; comparison with 5-Me modifications

The addition of the propynyl group at the 5 position of pyrimidine nucleotides is highly stabilising. We have determined the thermodynamic stability of the DNA.RNA hybrid r(GAAGAGAAGC)*d(GC(p)U(p)U(p)C(p)U(p) C(p)U(p)U(p)C) where p is the propynyl group at the 5 position and compared it with that of the unmodified duplex and the effects of methyl substitutions. The incorporation of the propyne group at the 5 position gives rise to a very large stabilisation of the hybrid duplex compared with the analogous 5-Me modification. The duplexes have been characterised by gel electrophoresis and NMR spectroscopy, which indicate that methyl substitutions have a smaller influence on local and global conformation than the propynyl groups. The increased NMR spectral dispersion of the propyne-modified duplex allowed a larger number of experimental restraints to be measured. Restrained molecular dynamics in a fully solvated system showed that the propyne modification leads to substantial conformational rearrangements stabilising a more A-like structure. The propynyl groups occupy a large part of the major groove and make favourable van der Waals interactions with their nearest neighbours and the atoms of the rings. This enhanced overlap may account at least in part for the increased thermodynamic stability. Furthermore, the simulations show a spine of hydration in the major groove as well as in the minor groove involving the RNA hydroxyl groups.

Detection of scalar couplings involving 2'-hydroxyl protons across hydrogen bonds in a frameshifting mRNA pseudoknot

The -1 frameshift-stimulating mRNA pseudoknot from pea enation mosaic virus-1 (PEMV-1) is composed nearly entirely of RNA triple helix.(4) The 2'-OH hydroxyl protons of riboses C15 and C16 are hydrogen bond donors to the N1 atoms of adenosines A27 and A25, respectively, positioned in the minor groove of pseudoknot stem S1. In this paper, a nonrefocused (1)H,(15)N CPMG HSQC of uniformly (13)C,(15)N-labeled 33-mer PEMV-1 RNA has been tailored to reveal a correlation of the 2'-OH hydroxyl proton of C15 to the N1 nitrogen resonance of A27 mediated by a cross hydrogen bond scalar coupling. The (1h)J(2'OH,N) cross hydrogen bond scalar coupling constant determined from a quantitative 1D (15N) spin-echo difference experiment for the C15/A27 interaction is 1.7 +/- 0.1 Hz, while that for the C16/A25 interaction appears larger, 3.5 +/- 0.3 Hz, despite the fact that the corresponding direct correlation between the 2'-OH hydroxyl proton of C16 and the N1 of A25 is missing due to unfavorable solvent exchange properties. These findings reveal a detailed picture of critical noncanonical O-H.N hydrogen-bonding loop-stem interactions in an RNA triple helical structure.

Detecting ligand binding to a small RNA target via saturation transfer difference NMR experiments in D(2)O and H(2)O

A small RNA motif is used as a target for ligand-based NMR-screening by saturation transfer difference (STD) NMR experiments. The prerequisites for using a small RNA target in STD experiments, such as saturation time, frequency, and pulses, are discussed. We also show that it is of advantage to use D2O as solvent instead of H2O due to the reduced R1 relaxation rate in D2O. The 27-nucleotide model of the ribosomal A-site was known to bind the aminoglycoside paromomycin with high affinity. This binding interaction could be detected easily, proving the effectiveness of STD NMR experiments as a screening tool for RNA-ligand interactions.

Melting of the solvent structure around a RNA duplex: a molecular dynamics simulation study

From three 2.4-ns molecular dynamics simulations of the r(CpG)(12) duplex conducted at 5, 25 and 37 degrees C, a strong temperature dependence of the dynamics of the water molecules and ions located in the first nucleic acid coordination shell is observed. At 5 degrees C, the highest residence times of bound water molecules exceed 1 ns while, at 37 degrees C, they decrease to 0.5 ns in agreement with available NMR data. Similar temperature dependencies are observed for the potassium ions bound to the duplex. In this temperature range, the structure of the RNA helix remains essentially unchanged. Thus, the observed alterations correspond to a 'premelting' of the solvent structure around the duplex. It is proposed that, before the nucleic acid structure melts, the entropy of the solvent increases to a point where it is no longer compensated by the enthalpic contribution of solute-solute and solute-solvent interactions. At this stage, the weakest structural elements start to melt. In other terms, the experimentally observed melting processes are preceded by a melting of the more labile solvent structure.

The solution structure of [d(CGC)r(amamam)d(TTTGCG)]2

The solution structure and hydration of a DNA.RNA hybrid chimeric duplex [d(CGC)r(amamam)d(TTTGCG)]2 in which the RNA adenines were substituted by 2'-O-methylated riboadenines was determined using two-dimensional NMR, simulated annealing, and restrained molecular dynamics. Only DNA residue 7T in the 2'-OMe-RNA.DNA junction adopted an O4'-endo sugar conformation, while the other DNA residues including 3C in the DNA.2'-OMe-RNA junction, adopted C1'-exo or C2'-endo conformations. The observed NOE intensity of 2'-O-methyl group to H1' proton of 4am at the DNA.2'-OMe-RNA junction is much weaker than those of 5am and 6am. The 2'-O-methyl group of 4am was found to orient towards the minor groove in the trans domain while the 2'-O-methyl groups of 5am and 6am were found to be in the gauche (+) domain. In contrast to the long-lived water molecules found close to the RNA adenine H2 and H1' protons and the methyl group of 7T in the RNA-DNA junction of [d(CGC)r(aaa)d(TTTGCG)]2, there were no long-lived water molecules found in [d(CGC)r(amamam)d(TTTGCG)]2. This is probably due to the hydrophobic enviroment created by the 2'-O-methylated riboadenines in the minor groove or due to the wider minor groove width in the middle of the structure. In addition, the 2'-O-methylation of riboadenines in pure chimeric duplex increses its melting temperature from 48.5 degrees C to 51.9 degrees C. The characteristic structural features and hydration patterns of this chimeric duplex provide a molecular basis for further therapeutic applications of DNA.RNA hybrid and chimeric duplexes with 2'-modified RNA residues.

Water and ion binding around r(UpA)12 and d(TpA)12 oligomers--comparison with RNA and DNA (CpG)12 duplexes

The structural and dynamic properties of the water and ion first coordination shell of the r(A-U) and d(A-T) base-pairs embedded within the r(UpA)12 and d(TpA)12 duplexes are described on the basis of two 2.4 ns molecular dynamics simulations performed in a neutralizing aqueous environment with 0.25 M added KCl. The results are compared to previous molecular dynamics simulations of the r(CpG)12 and d(CpG)12 structures performed under similar conditions. It can be concluded that: (i) RNA helices are more rigid than DNA helices of identical sequence, as reflected by the fact that RNA duplexes keep their initial A-form shape while DNA duplexes adopt more sequence-specific shapes. (ii) Around these base-pairs, the water molecules occupy 21 to 22 well-defined hydration sites, some of which are partially occupied by potassium ions. (iii) These hydration sites are occupied by an average of 21.9, 21.0, 20.1, and 19.8 solvent molecules (water and ions) around the r(G=C), r(A-U), d(G=C), and d(A-T) pairs, respectively. (iv) From a dynamic point of view, the stability of the hydration shell is the strongest for the r(G=C) pairs and the weakest for the d(A-T) pairs. (v) For RNA, the observed long-lived hydration patterns are essentially non-sequence dependent and involve water bridges located in the deep groove and linking OR atoms of adjacent phosphate groups. Maximum lifetimes are close to 400 ps. (vi) In contrast, for DNA, long-lived hydration patterns are sequence dependent and located in the minor groove. For d(CpG)12, water bridges linking the (G)N3 and (C)O2 with the O4' atoms of adjacent nucleotides with 400 ps maximum lifetimes are characterized while no such bridges are observed for d(TpA)12. (vii) Potassium ions are observed to bind preferentially to deep/major groove atoms at RpY steps, essentially d(GpC), r(GpC), and r(ApU), by forming ion-bridges between electronegative atoms of adjacent base-pairs. On average, about half an ion is observed per base-pair. Positive ion-binding determinants are related to the proximity of two or more electronegative atoms. Negative binding determinants are associated with the electrostatic and steric hindrance due to the proximity of electropositive amino groups and neutral methyl groups. Potassium ions form only transient contacts with phosphate groups.

Water and ion binding around RNA and DNA (C,G) oligomers

The dynamics, hydration, and ion-binding features of two duplexes, the A(r(CG)(12)) and the B(d(CG)(12)), in a neutralizing aqueous environment with 0.25 M added KCl have been investigated by molecular dynamics (MD) simulations. The regular repeats of the same C=G base-pair motif have been exploited as a statistical alternative to long MD simulations in order to extend the sampling of the conformational space. The trajectories demonstrate the larger flexibility of DNA compared to RNA helices. This flexibility results in less well defined hydration patterns around the DNA than around the RNA backbone atoms. Yet, 22 hydration sites are clearly characterized around both nucleic acid structures. With additional results from MD simulations, the following hydration scale for C=G pairs can be deduced: A-DNA<RNA (+3 H(2)O) and B-DNA<RNA (+2 H(2)O). The calculated residence times of water molecules in the first hydration shell of the helices range from 0.5 to 1 ns, in good agreement with available experimental data. Such water molecules are essentially found in the vicinity of the phosphate groups and in the DNA minor groove. The calculated number of ions that break into the first hydration shell of the nucleic acids is close to 0.5 per base-pair for both RNA and DNA. These ions form contacts essentially with the oxygen atoms of the phosphate groups and with the guanine N7 and O6 atoms; they display residence times in the deep/major groove approaching 500 ps. Further, a significant sequence-dependent effect on ion binding has been noted. Despite slight structural differences, K(+) binds essentially to GpC and not to CpG steps. These results may be of importance for understanding various sequence-dependent binding affinities. Additionally, the data help to rationalize the experimentally observed differences in gel electrophoretic mobility between RNA and DNA as due to the difference in hydration (two water molecules in favor of RNA) rather than to strong ion-binding features, which are largely similar for both nucleic acid structures.

The solution structure of [d(CGC)r(aaa)d(TTTGCG)](2): hybrid junctions flanked by DNA duplexes

The solution structure and hydration of the chimeric duplex [d(CGC)r(aaa)d(TTTGCG)](2), in which the central hybrid segment is flanked by DNA duplexes at both ends, was determined using two-dimensional NMR, simulated annealing and restrained molecular dynamics. The solution structure of this chimeric duplex differs from the previously determined X-ray structure of the analogous B-DNA duplex [d(CGCAAATTTGCG)](2)as well as NMR structure of the analogous A-RNA duplex [r(cgcaaauuugcg)](2). Long-lived water molecules with correlation time tau(c)longer than 0.3 ns were found close to the RNA adenine H2 and H1' protons in the hybrid segment. A possible long-lived water molecule was also detected close to the methyl group of 7T in the RNA-DNA junction but not with the other two thymines (8T and 9T). This result correlates with the structural studies that only DNA residue 7T in the RNA-DNA junction adopts an O4'-endo sugar conformation, while the other DNA residues including 3C in the DNA-RNA junction, adopt C1'-exo or C2'-endo conformations. The exchange rates for RNA C2'-OH were found to be approximately 5-20 s(-1). This slow exchange rate may be due to the narrow minor groove width of [d(CGC)r(aaa)d(TTTGCG)](2), which may trap the water molecules and restrict the dynamic motion of hydroxyl protons. The minor groove width of [d(CGC)r(aaa)d(TTTGCG)](2)is wider than its B-DNA analog but narrower than that of the A-RNA analog. It was further confirmed by its titration with the minor groove binding drug distamycin. A possible 2:1 binding mode was found by the titration experiments, suggesting that this chimeric duplex contains a wider minor groove than its B-DNA analog but still narrow enough to hold two distamycin molecules. These distinct structural features and hydration patterns of this chimeric duplex provide a molecular basis for further understanding the structure and recognition of DNA. RNA hybrid and chimeric duplexes.

Hydration of [d(CGC)r(aaa)d(TTTGCG)](2)

We have studied the hydration and dynamics of RNA C2'-OH in a DNA. RNA hybrid chimeric duplex [d(CGC)r(aaa)d(TTTGCG)](2). Long-lived water molecules with correlation time tau(c) larger than 0.3 ns were found close to the RNA adenine H2 and H1' protons in the hybrid segment. A possible long-lived water molecule was also detected close to the methyl group of 7T in the RNA-DNA junction but not to the other two thymine bases (8T and 9T). This result correlates with the structural studies that only DNA residue 7T in the RNA-DNA junction adopts an O4'-endo sugar conformation (intermediate between B-form and A-form), while the other DNA residues including 3C in the DNA-RNA junction, adopt C1'-exo or C2'-endo conformations (in the B-form domain). Based on the NOE cross-peak patterns, we have found that RNA C2'-OH tends to orient toward the O3' direction, forming a possible hydrogen bond with the 3'-phosphate group. The exchange rates for RNA C2'-OH were found to be around 5-20 s(-1), compared to 26.7(+/-13.8) s(-1) reported previously for the other DNA.RNA hybrid duplex. This slow exchange rate may be due to the narrow minor groove width of [d(CGC)r(aaa)d(TTTGCG)](2), which may trap the water molecules and restrict the dynamic motion of hydroxyl protons. The distinct hydration patterns of the RNA adenine H2 and H1' protons and the DNA 7T methyl group in the hybrid segment, as well as the orientation and dynamics of the RNA C2'-OH protons, may provide a molecular basis for further understanding the structure and recognition of DNA.RNA hybrid and chimeric duplexes.

Structure, dynamics and hydration of the nogalamycin-d(ATGCAT)2Complex determined by NMR and molecular dynamics simulations in solution

The structure of the 1:1 nogalamycin:d(ATGCAT)2 complex has been determined in solution from high-resolution NMR data and restrained molecular dynamics (rMD) simulations using an explicit solvation model. The antibiotic intercalates at the 5'-TpG step with the nogalose lying along the minor groove towards the centre of the duplex. Many drug-DNA nuclear Overhauser enhancements (NOEs) in the minor groove are indicative of hydrophobic interactions over the TGCA sequence. Steric occlusion prevents a second nogalamycin molecule from binding at the symmetry-related 5'-CpA site, leading to the conclusion that the observed binding orientation in this complex is the preferred orientation free of the complication of end-effects (drug molecules occupy terminal intercalation sites in all X-ray structures) or steric interactions between drug molecules (other NMR structures have two drug molecules bound in close proximity), as previously suggested. Fluctuations in key structural parameters such as rise, helical twist, slide, shift, buckle and sugar pucker have been examined from an analysis of the final 500 ps of a 1 ns rMD simulation, and reveal that many sequence-dependent structural features previously identified by comparison of different X-ray structures lie within the range of dynamic fluctuations observed in the MD simulations. Water density calculations on MD simulation data reveal a time-averaged pattern of hydration in both the major and minor groove, in good agreement with the extensive hydration observed in two related X-ray structures in which nogalamycin is bound at terminal 5'-TpG sites. However, the pattern of hydration determined from the sign and magnitude of NOE and ROE cross-peaks to water identified in 2D NOESY and ROESY experiments identifies only a few "bound" water molecules with long residence times. These solvate the charged bicycloaminoglucose sugar ring, suggesting an important role for water molecules in mediating drug-DNA electrostatic interactions within the major groove. The high density of water molecules found in the minor groove in X-ray structures and MD simulations is found to be associated with only weakly bound solvent in solution.

Unrestrained stochastic dynamics simulations of the UUCG tetraloop using an implicit solvation model

Three unrestrained stochastic dynamics simulations have been carried out on the RNA hairpin GGAC[UUCG] GUCC, using the AMBER94 force field (Cornell et al., 1995. J. Am. Chem. Soc. 117:5179-5197) in MacroModel 5.5 (Mohamadi et al., 1990. J. Comp. Chem. 11:440-467) and either the GB/SA continuum solvation model (Still et al., 1990. J. Am. Chem. Soc. 112:6127-6129) or a linear distance-dependent dielectric (1/R) treatment. The linear distance-dependent treatment results in severe distortion of the nucleic acid structure, restriction of all hydroxyl dihedrals, and collapse of the counterion atmosphere over the course of a 5-ns simulation. An additional vacuum simulation without counterions shows somewhat improved behavior. In contrast, the two GB/SA simulations (1.149 and 3.060 ns in length) give average structures within 1.2 A of the initial NMR structure and in excellent agreement with results of an earlier explicit solvent simulation (Miller and Kollman, 1997. J. Mol. Biol. 270:436-450). In a 3-ns GB/SA simulation starting with the incorrect UUCG tetraloop structure (Cheong et al., 1990. Nature. 346:680-682), this loop conformation converts to the correct loop geometry (Allain and Varani, 1995. J. Mol. Biol. 250:333-353), suggesting enhanced sampling relative to the previous explicit solvent simulation. Thermodynamic effects of 2'-deoxyribose substitutions of loop nucleotides were experimentally determined and are found to correlate with the fraction of time the ribose 2'-OH is hydrogen bonded and the distribution of the hydroxyl dihedral is observed in the GB/SA simulations. The GB/SA simulations thus appear to faithfully represent structural features of the RNA without the computational expense of explicit solvent.

The residence time of the bound water in the hydrophobic minor groove of the carbocyclic-nucleoside analogs of Dickerson-Drew dodecamers

The residence time of the bound water molecules in the antisense oligodeoxyribonucleotides containing 7'-alpha-methyl (TMe) carbocyclic thymidines in duplex (I), d5'(1C2G3C4G5A6A7TMc8TMe9C10G11C12G)23', and 6'-alpha-hydroxy (TOH) carbocyclic thymidines in duplex (II), d5'(1C2G3C4G5AOH6AOH7TOH8 TOH9C10G11C12G)2(3), have been investigated using a combination of NOESY and ROESY experiments. Because of the presence of 7'-alpha-methyl groups of TMe in the centre of the minor groove in duplex (I), the residence time of the bound water molecule is shorter than 0.3 ns. The dramatic reduction of the residence time of the water molecule in the minor groove in duplex (I) compared with the natural counterpart has been attributed to the replacement of second shell of hydration and disruption of hydrogen-bonding with 04' in the minor groove by hydrophobic alpha-methyl groups, as originally observed in the X-ray study. This effect could not be attributed to the change of the width of the minor groove because a comparative NMR study of the duplex (I) and its natural counterpart showed that the widths of their minor grooves are more or less unchanged (r.m.s.d change in the core part is <0.63A). For duplex (II) with polar 6'-alpha-hydroxyl groups pointed to the minor groove, the correlation time is much longer than 0.36 ns as a result of the stabilising hydrogen-bonding interaction with N3 or 02 of the neighbouring nucleotides.

The orientation and dynamics of the C2'-OH and hydration of RNA and DNA.RNA hybrids

The stereochemical and dynamic properties of the C2' hydroxyl group in several DNA.RNA hybrids have been measured by NMR and compared with the homologous RNA duplex. The C2'-OH NMR signals of the RNA strands were identified, and numerous specific assignments were made. The rate constants for exchange of the hydroxyl protons with water were determined at 5 degrees C, and were found to depend on both the position within a particular sequence and the nature of the duplex. On average, the exchange rate constants were slowest for the hybrids of composition rR.dY, and fastest for the RNA duplex, with an overall range of approximately 10-50/s. In the DNA.RNA hybrids, strong NOEs and ROEs were observed between the OH and the H1' of the same sugar, unambiguously showing that the OH proton points toward the H1' most of the time, and not toward the O3' of the same sugar. Evidence for significant hydration in both grooves of the DNA.RNA hybrids and the DNA duplex was found in ROESY and NOESY experiments. On average, the minor groove of the DNA.RNA hybrids showed more kinetically significant hydration than the DNA, which can be attributed to the hydrophilic lining of hydroxyl groups in RNA.

Conformational and dynamic properties of a 14 residue antifreeze glycopeptide from Antarctic cod

The 1H and 13C NMR spectra of a 14-residue antifreeze glycopeptide from Antarctic cod (Tetramatomnus borchgrevinki) containing two proline residues have been assigned. 13C NMR relaxation experiments indicate motional anisotropy of the peptide, with a tumbling time in water at 5 degrees C of 3-4 ns. The relaxation data and lack of long-range NOEs are consistent with a linear peptide undergoing significant segmental motion. However, extreme values of some coupling constants and strong sequential NOEs indicate regions of local order, which are most evident at the two ATPA subsequences. Similar spectroscopic properties were observed in the 16-residue analogue containing an Arg-Ala dipeptide added to the C-terminus. Molecular modeling also showed no evidence of long-range order, but the two ATPA subsequences were relatively well determined by the experimental data. These motifs were quite distinct from helical structures or beta turns commonly found in proteins, but rather resemble sections of an extended polyproline helix. Thus, the NMR data provide a description of the local order, which is of relevance to the mechanism of action of the antifreeze activity of the antifreeze glycopeptides as well as their ability to protect cells during hypothermic storage.

Rules governing the orientation of the 2'-hydroxyl group in RNA

Molecular dynamics simulations reveal that, in C3'-endo sugar puckers, only three orientations are accessible to the 2'-hydroxyl groups distinctive of RNA molecules: towards (i) the O3', (ii) the O4' of the same sugar, and (iii) the shallow groove base atoms. In the rarer C2'-endo sugar puckers, orientations towards the O3' atom of the same sugar are strongly favoured. Surprisingly, in helical regions, the frequently suggested intra-strand O2'-H(n)...O4'(n+1) interaction is not found. This observation led to the detection of an axial C-H...O interaction between the C2'-H2'(n) group and the O4'(n+1) atom contributing to the stabilization of RNA helical regions. Subsequent analysis of crystallographic structures of both RNA and A-DNA helices fully supports this finding. Specific hydration patterns are also thought to play a significant role in the stabilization of RNA structures. In the shallow groove of RNA, known as a favourable RNA or protein-binding region, three well-defined hydration sites are located around the O2' atom. These hydration sites, occupied by water molecules exchanging with the bulk, constitute, after dehydration, anchor points for specific interactions between RNA and nucleic acids, proteins or drugs. Therefore, the fact that the 2'-hydroxyl group is not monopolised by axial stabilization, together with its water-like behaviour, facilitates complex formation involving RNA helical regions.

Comparison of the electrophoretic and hydrodynamic properties of DNA and RNA oligonucleotide duplexes

The electrophoretic behavior of defined DNA and RNA oligonucleotide duplexes from 10 to 20 bp in length has been investigated as a function of salt conditions, gel concentration, and temperature. The RNA oligomers migrated much more slowly than the DNA oligomers of the same sequence under all conditions. From sedimentation equilibrium and velocity measurements, the apparent partial specific volume in 0.1 M KCI, 20 mM NaPi, pH 7, was determined as 0.56 +/- 0.015 ml g(-1) for DNA and 0.508 ml g(-1) for RNA. The translational friction coefficients were determined and compared with the values calculated for cylinders. Taking into account the shape factors, the solution density, and partial specific volumes, the effective degree of hydration was estimated as 0.8-1 g g(-1) DNA. There was no significant difference in the frictional coefficients of the DNA and RNA oligomers, indicating that the effective sizes of DNA and RNA are very similar in solution. The differential electrophoretic mobility of DNA and RNA must arise from the differences in interaction with counterions, which is probably a global property of the oligonucleotides.

Conformational properties and thermodynamics of the RNA duplex r(CGCAAAUUUGCG)2: comparison with the DNA analogue d(CGCAAATTTGCG)2

The thermodynamic stability of nine dodecamers (four DNA and five RNA) of the same base composition has been compared by UV-melting. TheDeltaG of stabilisation were in the order: r(GACUGAUCAGUC)2>r(CGCAAATTTGCG)2 approximately r(CGCAUAUAUGCG)2>d(CGCAAATTTGCG)2 approximately r(CGCAAAUUUGCG)2>d(CGCATATATGCG)2 approximately d(GACTGATCAGTC)2>r(CGCUUUAAAGCG)2 approximately d(CGCTTTAAAGCG)2. Compared with the mixed sequences, both r(AAAUUU) and r(UUUAAA) are greatly destablising in RNA, whereas in DNA, d(TTTAAA) is destabilising but d(AAATTT) is stabilising, which has been attributed to the formation of a special B'structure involving large propeller twists of the A-T base pairs. The solution structure of the RNA dodecamer r(CGCAAAUUUGCG)2has been determined using NMR and restrained molecular dynamics calculations to assess the conformational reasons for its stability in comparison with d(CGCAAATTTGCG)2. The structures refined to a mean pairwise r.m.s.d. of 0.89+/-0.29 A. The nucleotide conformations are typical of the A family of structures. However, although the helix axis displacement is approximately 4.6 A into the major groove, the rise (3.0 A) and base inclination ( approximately 6 degrees ) are different from standard A form RNA. The extensive base-stacking found in the AAATTT tract of the DNA homologue that is largely responsible for the higher thermodynamic stability of the DNA duplex is reduced in the RNA structure, which may account for its low relative stability.

RNA hydration: three nanoseconds of multiple molecular dynamics simulations of the solvated tRNA(Asp) anticodon hairpin

The hydration of the tRNA(Asp) anticodon hairpin was investigated through the analysis of six 500 ps multiple molecular dynamics (MMD) trajectories generated by using the particle mesh Ewald method for the treatment of the long-range electrostatic interactions. Although similar in their dynamical characteristics, these six trajectories display different local hydration patterns reflecting the landscape of the "theoretical" conformational space being explored. The statistical view gained through the MMD strategy allowed us to characterize the hydration patterns around important RNA structural motifs such as a G-U base-pair, the anticodon U-turn, and two modified bases: pseudouridine and 1-methylguanine. The binding of ammonium counterions to the hairpin has also been investigated. No long-lived hydrogen bond between water and a 2'-hydroxyl has been observed. Water molecules with long-residence times are found bridging adjacent pro-Rp phosphate atoms. The conformation of the pseudouridine is stiffened by a water-mediated base-backbone interaction and the 1-methylguanine is additionally stabilized by long-lived hydration patterns. Such long-lived hydration patterns are essential to ensure the structural integrity of this hairpin motif. Consequently, our simulations confirm the conclusion reached from an analysis of X-ray crystal structures according to which water molecules form an integral part of nucleic acid structure. The fact that the same conclusion is reached from a static and a dynamic point of view suggests that RNA and water together constitute the biologically relevant functional entity.

Solution structure of r(gaggacug):d(CAGTCCTC) hybrid: implications for the initiation of HIV-1 (+)-strand synthesis

The three-dimensional solution structure of the hybrid duplex r(gaggacug):d(CAGTCCTC) has been determined by two-dimensional NMR, distance geometry (DG), restrained molecular dynamics (rMD) and NOE back-calculation methods. This hybrid, consisting of a purine-rich RNA strand and a pyrimidine-rich DNA strand, is related to the polypurine (+)-strand primer formed after (-)-strand DNA synthesis and RNase H degradation of the viral RNA strand and contains the site of a specific cleavage by reverse transcription (RT) RNase H at the end of the HIV-1 polypurine tract. This polypurine primer is an important intermediate in the formation of virally encoded double-stranded DNA prior to HIV-1 retrovirus integration. The correct processing of this primer is vital in the life cycle of the human immunodeficiency virus type (HIV-1) retrovirus. The structure of the r(gaggacug):d(CAGTCCTC) hybrid, as determined in solution by NMR, is intermediate between canonical A-type and B-type double helices, and has mixed structural characteristics. It is quantitatively different from the previously determined solution structures of other RNA-DNA hybrids, particularly in the width and shape of the major groove, which is wider than the major groove of other hybrids and is close to the dimension of the major groove of B-type DNA duplexes. The structure of this hybrid duplex contains a prominent bend in the double helix with a magnitude and direction similar to the bend in Okazaki fragments. The structural features of the present duplex may explain the unique interactions of this sequence with HIV-1 RT during both (-)-strand and (+)-strand DNA synthesis.