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

Modulation of CRISPR-Cas9 Cleavage with an Oligo-Ribonucleoprotein Design

Clustered regularly interspaced short palindromic repeats (CRISPR) associated protein Cas9 system has been widely used for genome editing. However, the editing or cleavage specificity of CRISPR Cas9 remains a major concern due to the off-target effects. The existing approaches to control or modulate CRISPR Cas9 cleavage include engineering Cas9 protein and development of anti-CRISPR proteins. There are also attempts on direct modification of sgRNA, for example, structural modification via truncation or hairpin design, or chemical modification on sgRNA such as partially replacing RNA with DNA. The above-mentioned strategies rely on extensive protein engineering and direct chemical or structural modification of sgRNA. In this study, we proposed an indirect method to modulate CRISPR Cas9 cleavage without modification on Cas9 protein or sgRNA. An oligonucleotide was used to form an RNA-DNA hybrid structure with the sgRNA spacer, creating steric hindrance during the Cas9 mediated DNA cleavage process. We first introduced a simple and robust method to assemble the oligo-ribonucleoprotein (oligo-RNP). Next, the cleavage efficiency of the assembled oligo-RNP was examined using different oligo lengths in vitro. Lastly, we showed that the oligo-RNP directly delivered into cells could also modulate Cas9 activity inside cells using three model gene targets with reduced off-target effects.

Unraveling DNA Triplex Assembly: Mass Spectrometric Investigation of Modified Triplex Forming Oligonucleotides for Enhanced Gene Targeting

Deoxyribonucleic acid triplexes have potential roles in a range of biological processes involving gene and transcriptional regulation. A major challenge in exploiting the formation of these higher-order structures to target genes in vivo is their low stability, which is dependent on many factors including the length and composition of bases in the sequence. Here, different DNA base modifications have been explored, primarily using native mass spectrometry, in efforts to enable stronger binding between the triplex forming oligonucleotide (TFO) and duplex target sites. These modifications can also be used to overcome pyrimidine interruptions in the duplex sequence in promoter regions of genomes, to expand triplex target sequences for antigene therapies. Using model sequences with a single pyrimidine interruption, triplex forming oligonucleotides containing locked nucleic acid base modifications were shown to have a higher triplex binding propensity than DNA-only and dSpacer-containing TFOs. However, the triplex forming ability of these systems was limited by the competitive formation of multiple higher order assemblies. Triplex forming sequences that correspond to specific gene targets from the Pseudomonas aeruginosa genome were also investigated, with LNA-containing TFOs the only variant able to form triplex using these sequences. This work indicates the advantages of utilizing synthetically modified TFOs to form triplex assemblies in vivo for potential therapeutic applications and highlights the advantages of native mass spectrometry for the study of their formation.

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.

Iso-FRET: an isothermal competition assay to analyze quadruplex formation in vitro

Algorithms have been widely used to predict G-quadruplexes (G4s)-prone sequences. However, an experimental validation of these predictions is generally required. We previously reported a high-throughput technique to evidence G4 formation in vitro called FRET-MC. This method, while convenient and reproducible, has one known weakness: its inability to pin point G4 motifs of low thermal stability. As such quadruplexes may still be biologically relevant if formed at physiological temperature, we wanted to develop an independent assay to overcome this limitation. To this aim, we introduced an isothermal version of the competition assay, called iso-FRET, based on a duplex-quadruplex competition and a well-characterized bis-quinolinium G4 ligand, PhenDC3. G4-forming competitors act as decoys for PhenDC3, lowering its ability to stabilize the G4-forming motif reporter oligonucleotide conjugated to a fluorescence quencher (37Q). The decrease in available G4 ligand concentration restores the ability of 37Q to hybridize to its FAM-labeled short complementary C-rich strand (F22), leading to a decrease in fluorescence signal. In contrast, when no G4-forming competitor is present, PhenDC3 remains available to stabilize the 37Q quadruplex, preventing the formation of the F22 + 37Q complex. Iso-FRET was first applied to a reference panel of 70 sequences, and then used to investigate 23 different viral sequences.

Biomolecules of 2-Thiouracil, 4-Thiouracil and 2,4-Dithiouracil: A DFT Study of the Hydration, Molecular Docking and Effect in DNA:RNAMicrohelixes

The molecular structure of 2-thiouracil, 4-thiouracil and 2,4-dithiouracil was analyzed under the effect of the first and second hydration shell by using the B3LYP density functional (DFT) method, and the results were compared to those obtained for the uracil molecule. A slight difference in the water distribution appears in these molecules. On the hydration of these molecules several trends in bond lengths and atomic charges were established. The ring in uracil molecule appears easier to be deformed and adapted to different environments as compared to that when it is thio-substituted. Molecular docking calculations of 2-thiouracil against three different pathogens: Bacillus subtilis, Escherichia coli and Candida albicans were carried out. Docking calculations of 2,4-dithiouracil ligand with various targeted proteins were also performed. Different DNA: RNA hybrid microhelixes with uridine, 2-thiouridine, 4-thiouridine and 2,4-dithiouridine nucleosides were optimized in a simple model with three nucleotide base pairs. Two main types of microhelixes were analyzed in detail depending on the intramolecular H-bond of the 2'-OH group. The weaker Watson-Crick (WC) base pair formed with thio-substituted uracil than with unsubstituted ones slightly deforms the helical and backbone parameters, especially with 2,4-dithiouridine. However, the thio-substitution significantly increases the dipole moment of the A-type microhelixes, as well as the rise and propeller twist parameters.

Effect of the sulfur atom on S2 and S4 positions of the uracil ring in different DNA:RNA hybrid microhelixes with three nucleotide base pairs

The effect of the sulphur atom on the uracil ring was analyzed in different DNA:RNA microhelixes with three nucleotide base-pairs, including uridine, 2-thiouridine, 4-thiouridine, 2,4-dithiouridine, cytidine, adenosine and guanosine. Distinct backbone and helical parameters were optimized at different density functional (DFT) levels. The Watson-Crick pair with 2-thiouridine appears weaker than with uridine, but its interaction with water molecules appears easier. Two types of microhelixes were found, depending on the H-bond of H2' hydroxyl atom: A-type appears with the ribose ring in ³ E-envelope C_3' -endo, and B-type in ² E-envelope C_2' -endo. B-type is less common but it is more stable and with higher dipole-moment. The sulphur atoms significantly increase the dipole-moment of the microhelix, as well as the rise and propeller twist parameters. Simulations with four Na atoms H-bonded to the phosphate groups, and further hydration with explicit water molecules were carried out. A re-definition of the numerical value calculation of several base-pair and base-stacking parameters is suggested.

Partial DNA-guided Cas9 enables genome editing with reduced off-target activity

CRISPR-Cas9 is a versatile RNA-guided genome editing tool. Here we demonstrate that partial replacement of RNA nucleotides with DNA nucleotides in CRISPR RNA (crRNA) enables efficient gene editing in human cells. This strategy of partial DNA replacement retains on-target activity when used with both crRNA and sgRNA, as well as with multiple guide sequences. Partial DNA replacement also works for crRNA of Cpf1, another CRISPR system. We find that partial DNA replacement in the guide sequence significantly reduces off-target genome editing through focused analysis of off-target cleavage, measurement of mismatch tolerance and genome-wide profiling of off-target sites. Using the structure of the Cas9-sgRNA complex as a guide, the majority of the 3' end of crRNA can be replaced with DNA nucleotide, and the 5 - and 3'-DNA-replaced crRNA enables efficient genome editing. Cas9 guided by a DNA-RNA chimera may provide a generalized strategy to reduce both the cost and the off-target genome editing in human cells.

RNA-based boronate internucleosidic linkages: an entry into reversible templated ligation and loop formation

The synthesis of a 5′-boronoribonucleotidic phosphoramidite building block has been achieved and incorporated at the 5′ extremities of RNA sequences for the templated assembly of RNA shortmers.

Direct Observation of Carbohydrate Hydroxyl Protons in Hydrogen Bonds with a Protein

Hydroxyl proton resonances of uniformly ¹³C-labeled Manα(1-2)Manα(1-2)ManαOMe (Man₃) bound to cyanovirin-N (CV-N) were detected at ambient temperature in aqueous solution by NMR spectroscopy. The directions of the hydroxyl groups were determined on the basis of NOEs, and a previously unknown hydrogen-bonding network between Man₃ and CV-N was discovered. This is the first report on detecting hydroxyl protons of a protein-bound carbohydrate in aqueous solution by NMR. Approaches such as those presented here may open the door for accurately determining intermolecular hydrogen bonds in carbohydrate-protein complexes.

RNA Structures as Mediators of Neurological Diseases and as Drug Targets

RNAs adopt diverse folded structures that are essential for function and thus play critical roles in cellular biology. A striking example of this is the ribosome, a complex, three-dimensionally folded macromolecular machine that orchestrates protein synthesis. Advances in RNA biochemistry, structural and molecular biology, and bioinformatics have revealed other non-coding RNAs whose functions are dictated by their structure. It is not surprising that aberrantly folded RNA structures contribute to disease. In this Review, we provide a brief introduction into RNA structural biology and then describe how RNA structures function in cells and cause or contribute to neurological disease. Finally, we highlight successful applications of rational design principles to provide chemical probes and lead compounds targeting structured RNAs. Based on several examples of well-characterized RNA-driven neurological disorders, we demonstrate how designed small molecules can facilitate the study of RNA dysfunction, elucidating previously unknown roles for RNA in disease, and provide lead therapeutics.

Energies and 2'-Hydroxyl Group Orientations of RNA Backbone Conformations. Benchmark CCSD(T)/CBS Database, Electronic Analysis, and Assessment of DFT Methods and MD Simulations

Sugar-phosphate backbone is an electronically complex molecular segment imparting RNA molecules high flexibility and architectonic heterogeneity necessary for their biological functions. The structural variability of RNA molecules is amplified by the presence of the 2'-hydroxyl group, capable of forming multitude of intra- and intermolecular interactions. Bioinformatics studies based on X-ray structure database revealed that RNA backbone samples at least 46 substates known as rotameric families. The present study provides a comprehensive analysis of RNA backbone conformational preferences and 2'-hydroxyl group orientations. First, we create a benchmark database of estimated CCSD(T)/CBS relative energies of all rotameric families and test performance of dispersion-corrected DFT-D3 methods and molecular mechanics in vacuum and in continuum solvent. The performance of the DFT-D3 methods is in general quite satisfactory. The B-LYP-D3 method provides the best trade-off between accuracy and computational demands. B3-LYP-D3 slightly outperforms the new PW6B95-D3 and MPW1B95-D3 and is the second most accurate density functional of the study. The best agreement with CCSD(T)/CBS is provided by DSD-B-LYP-D3 double-hybrid functional, although its large-scale applications may be limited by high computational costs. Molecular mechanics does not reproduce the fine energy differences between the RNA backbone substates. We also demonstrate that the differences in the magnitude of the hyperconjugation effect do not correlate with the energy ranking of the backbone conformations. Further, we investigated the 2'-hydroxyl group orientation preferences. For all families, we conducted a QM and MM hydroxyl group rigid scan in gas phase and solvent. We then carried out set of explicit solvent MD simulations of folded RNAs and analyze 2'-hydroxyl group orientations of different backbone families in MD. The solvent energy profiles determined primarily by the sugar pucker match well with the distribution data derived from the simulations. The QM and MM energy profiles predict the same 2'-hydroxyl group orientation preferences. Finally, we demonstrate that the high energy of unfavorable and rarely sampled 2'-hydroxyl group orientations can be attributed to clashes between occupied orbitals.

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.

High-resolution NMR structure of an RNA model system: the 14-mer cUUCGg tetraloop hairpin RNA

We present a high-resolution nuclear magnetic resonance (NMR) solution structure of a 14-mer RNA hairpin capped by cUUCGg tetraloop. This short and very stable RNA presents an important model system for the study of RNA structure and dynamics using NMR spectroscopy, molecular dynamics (MD) simulations and RNA force-field development. The extraordinary high precision of the structure (root mean square deviation of 0.3 A) could be achieved by measuring and incorporating all currently accessible NMR parameters, including distances derived from nuclear Overhauser effect (NOE) intensities, torsion-angle dependent homonuclear and heteronuclear scalar coupling constants, projection-angle-dependent cross-correlated relaxation rates and residual dipolar couplings. The structure calculations were performed with the program CNS using the ARIA setup and protocols. The structure quality was further improved by a final refinement in explicit water using OPLS force field parameters for non-bonded interactions and charges. In addition, the 2'-hydroxyl groups have been assigned and their conformation has been analyzed based on NOE contacts. The structure currently defines a benchmark for the precision and accuracy amenable to RNA structure determination by NMR spectroscopy. Here, we discuss the impact of various NMR restraints on structure quality and discuss in detail the dynamics of this system as previously determined.

Effect of a water molecule on the sugar puckering of uridine, 2'-deoxyuridine, and 2'-O-methyl uridine inserted in duplexes

We used high-level quantum mechanical calculations to determine the pucker (north type or south type) of various compounds: uridine, 2'-deoxyuridine, and 2'-O-methyl uridine. Although the dihedrals of the backbone are set close to their experimental values in double-stranded nucleic acids, calculations using density functional theory show that, in vacuo or in a continuum mimicking the dielectric properties of water, the south puckering conformations of uridine is favored. This contrasts with experimental data: most ribonucleosides inserted into a duplex have the north puckering. We show here that the north puckering is favored when an explicit water molecule is introduced into the calculation. The orientations of the 2' group and of the water molecule have implications for the prevalence of the north puckering. We studied several orientations of the water molecule binding uracil O2 and the 2' group and estimated the energy barriers in the path between the north-to-south conformations. The north puckering is more favored in 2'-OH than in 2'-OCH3 compounds in the presence of the explicit water molecule.

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.

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.

NMR structure of the 3' stem-loop from human U4 snRNA

The NMR structure of the 3' stem-loop (3'SL) from human U4 snRNA was determined to gain insight into the structural basis for conservation of this stem-loop sequence from vertebrates. 3'SL sequences from human, rat, mouse and chicken U4 snRNA each consist of a 7 bp stem capped by a UACG tetraloop. No high resolution structure has previously been reported for a UACG tetraloop. The UACG tetraloop portion of the 3'SL was especially well defined by the NMR data, with a total of 92 NOE-derived restraints (about 15 per residue), including 48 inter-residue restraints (about 8 per residue) for the tetraloop and closing C-G base pair. Distance restraints were derived from NOESY spectra using MARDIGRAS with random error analysis. Refinement of the 20mer RNA hairpin structure was carried out using the programs DYANA and miniCarlo. In the UACG tetraloop, U and G formed a base pair stabilized by two hydrogen bonds, one between the 2'-hydroxyl proton of U and carbonyl oxygen of G, another between the imino proton of G and carbonyl oxygen O2 of U. In addition, the amino group of C formed a hydrogen bond with the phosphate oxygen of A. G adopted a syn orientation about the glycosidic bond, while the sugar puckers of A and C were either C2'-endo or flexible. The conformation of the UACG tetraloop was, overall, similar to that previously reported for UUCG tetraloops, another member of the UNCG class of tetraloops. The presence of an A, rather than a U, at the variable position, however, presents a distinct surface for interaction of the 3'SL tetraloop with either RNA or protein residues that may stabilize interactions important for active spliceosome formation. Such tertiary interactions may explain the conservation of the UACG tetraloop motif in 3'SL sequences from U4 snRNA in vertebrates.

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 RNA i-motif

Oligodeoxynucleotides with stretches of cytidine residues associate into a four-stranded structure, the i-motif, in which two head-to-tail, intercalated, parallel-stranded duplexes are held together by hemiprotonated C.C+ pairs. We have investigated the possibility of forming an i-motif structure with C-rich ribonucleic acids. The four C-rich RNAs studied, r(UC5), r(C5), r(C5U) and r(UC3), associate into multiple intercalated structures at acidic pH. r(UC5) forms two i-motif structures that differ by their intercalation topologies. We report on a structural study of the main form and we analyze the small conformational differences found by comparison with the DNA i-motif. The stacking topology of the main structure avoids one of the six 2'-OH/2'-OH repulsive contacts expected in a fully intercalated structure. The C3'-endo pucker of the RNA sugars and the orientation of the intercalated C.C+ pairs result in a modest widening of the narrow grooves at the steps where the hydroxyl groups are in close contact. The free energy of the RNA i-motif, on average -4 kJ mol(-1) per C.C+ pair, is half of the value found in DNA i-motif structures.

Kinetics of formation of hypoxanthine containing base pairs by HIV-RT: RNA template effects on the base substitution frequencies

Hypoxanthine (H), the deamination product of adenine, has been implicated in the high frequency of A to G transitions observed in retroviral and other RNA genomes. Although H.C base pairs are thermodynamically more stable than other H.N pairs, polymerase selection may be determined in part by kinetic factors. Therefore, the hypoxanthine induced substitution pattern resulting from replication by viral polymerases may be more complex than that predicted from thermodynamics. We have examined the steady-state kinetics of formation of base pairs opposite template H in RNA by HIV-RT, and for the incorporation of dITP during first- and second-strand synthesis. Hypoxanthine in an RNA template enhances the k(2app) for pairing with standard dNTPs by factors of 10-1000 relative to adenine at the same sequence position. The order of base pairing preferences for H in RNA was observed to be H.C >> H.T > H.A > H.G. Steady-state kinetics of insertion for all possible mispairs formed with dITP were examined on RNA and DNA templates of identical sequence. Insertion of dITP opposite all bases occurs 2-20 times more frequently on RNA templates. This bias for higher insertion frequencies on RNA relative to DNA templates is also observed for formation of mispairs at template A. This kinetic advantage afforded by RNA templates for mismatches and pairing involving H suggests a higher induction of mutations at adenines during first-strand synthesis by HIV-RT.

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.

Comparison of the solution conformation and dynamics of antifreeze glycoproteins from Antarctic fish

The (1)H- and (13)C-NMR spectra of antifreeze glycoprotein fractions 1-5 from Antarctic cod have been assigned, and the dynamics have been measured using (13)C relaxation at two temperatures. The chemical shifts and absence of non-sequential (1)H-(1)H NOEs are inconsistent with a folded, compact structure. (13)C relaxation measurements show that the protein has no significant long-range order, and that the local correlation times are adequately described by a random coil model. Hydroxyl protons of the sugar residues were observed at low temperature, and the presence of exchange-mediated ROEs to the sugar indicate extensive hydration. The conformational properties of AFGP1-5 are compared with those of the previously examined 14-mer analog AFGP8, which contains proline residues in place of some alanine residues (Lane, A. N., L. M. Hays, R. E. Feeney, L. M. Crowe, and J. H. Crowe. 1998. Protein Sci. 7:1555-1563). The infrared (IR) spectra of AFGP8 and AFGP1-5 in the amide I region are quite different. The presence of a wide distribution of backbone torsion angles in AFGP1-5 leads to a rich spectrum of frequencies in the IR spectrum, as interconversion among conformational states is slow on the IR frequency time scale. However, these transitions are fast on the NMR chemical shift time scales. The restricted motions for AFGP8 may imply a narrower distribution of possible o, psi angles, as is observed in the IR spectrum. This has significance for attempts to quantify secondary structures of proteins by IR in the presence of extensive loops.

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.

A molecular dynamics simulation of the flavin mononucleotide-RNA aptamer complex

We report on an unrestrained molecular dynamics simulation of the flavin mononucleotide (FMN)-RNA aptamer. The simulated average structure maintains both cross-strand and intermolecular FMN-RNA nuclear Overhauser effects from the nmr experiments and has all qualitative features of the nmr structure including the G10-U12-A25 base triple and the A13-G24, A8-G28, and G9-G27 mismatches. However, the relative orientation of the hairpin loop to the remaining part of the molecule differs from the nmr structure. The simulation predicts that the flexible phosphoglycerol part of FMN moves toward G27 and forms hydrogen bonds. There are structurally long-lived water molecules in the FMN binding pocket forming hydrogen bonds within FMN and between FMN and RNA. In addition, long-lived water is found bridging primarily RNA backbone atoms. A general feature of the environment of long-lived "structural" water is at least two and in most cases three or four potential acceptor atoms. The 2'-OH group of RNA usually acts as an acceptor in interactions with the solvent. There are almost no intrastrand O2'H(n) vertical ...O4'(n + 1) hydrogen bonds within the RNA backbone. In the standard case the preferred orientation of the 2'-OH hydrogen atoms is approximately toward O3' of the same nucleotide. However, a relatively large number of conformations with the backbone torsional angle gamma in the trans orientation is found. A survey of all experimental RNA x-ray structures shows that this backbone conformation occurs but is less frequent than found in the simulation. Experimental nmr RNA aptamer structures have a higher fraction of this conformation as compared to the x-ray structures. The backbone conformation of nucleotide n + 1 with the torsional angle gamma in the trans orientation leads to a relatively short distance between 2'-OH(n) and O5'(n + 1), enabling hydrogen-bond formation. In this case the preferred orientation of the 2'-OH hydrogen atom is approximately toward O5'(n + 1). We find two relatively short and dynamically stable types of backbone-backbone next-neighbor contacts, namely C2'(H)(n) vertical ...O4'(n + 1) and C5'(H)(n + 1) vertical ...O2'(n). These interactions may affect both backbone rigidity and thermodynamic stability of RNA helical structures.