Stability and structure of RNA duplexes containing isoguanosine and isocytidine.

Isoguanosine (iG) and isocytidine (iC) differ from guanosine (G) and cytidine (C), respectively, in that the amino and carbonyl groups are transposed. The thermodynamic properties of a set of iG, iC containing RNA duplexes have been measured by UV optical melting. It is found that iG-iC replacements usually stabilize duplexes, and the stabilization per iG-iC pair is sequence-dependent. The sequence dependence can be fit to a nearest-neighbor model in which the stabilities of iG--iC pairs depend on the adjacent iG--iC or G--C pairs. For 5'-CG-3'/3'-GC-5' and 5'-GG-3'/3'-CC-5' nearest neighbors, the free energy differences upon iG-iC replacement are smaller than 0.2 kcal/mol at 37 degrees C, regardless of the number of replacements. For 5'-GC-3'/3'-CG-5', however, each iG--iC replacement adds 0.6 kcal/mol stabilizing free energy at 37 degrees C. Stacking propensities of iG and iC as unpaired nucleotides at the end of a duplex are similar to those of G and C. An NMR structure is reported for r(CiGCGiCG)(2) and found to belong to the A-form family. The structure has substantial deviations from standard A-form but is similar to published NMR and/or crystal structures for r(CGCGCG)(2) and 2'-O-methyl (CGCGCG)(2). These results provide benchmarks for theoretical calculations aimed at understanding the fundamental physical basis for the thermodynamic stabilities of nucleic acid duplexes.     

Understanding the role of base stacking in nucleic acids. MD and QM analysis of tandem GA base pairs in RNA duplexes

Preceding NMR experiments show that the conformation of tandem GA base pairs, an important recurrent non-canonical building block in RNA duplexes, is context dependent. The GA base pairs adopt "sheared" N3(G)-N6(A), N2(G)-N7(A) geometry in the r(CGAG)(2) and r(iGGAiC)(2) contexts while switching to "imino" N1(G)-N1(A), O6(G)-N6(A) geometry in the r(GGAC)(2) and r(iCGAiG)(2) contexts (iC and iG stand for isocytosine and isoguanine, respectively). As base stacking is likely to be one of the key sources of the context dependence of the conformation of GA base pairs, we calculated base stacking energies in duplexes containing such base pairs, to see if this dependence can be predicted by stacking energy calculations. When investigating the context dependence of the GA geometry two different conformations of the same duplex were compared (imino vs. sheared). The geometries were generated via explicit solvent MD simulations of the respective RNA duplexes, while the subsequent QM energy calculations focused on base stacking interactions of the four internal base pairs. Geometrical relaxation of nucleobase atoms prior to the stacking energy computations has a non-negligible effect on the results. The stacking energies were derived at the DFT-D/6-311++G(3df,3pd) level. We show a rather good correspondence between the intrinsic gas-phase stacking energies and the NMR-determined GA geometries. The conformation with more favorable gas-phase stacking is in most cases the one observed in experiments. This correlation is not improved when including solvent effects via the COSMO method. On the other side, the stacking calculations do not predict the relative thermodynamic stability of duplex formation for different sequences.     

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.     

Crystal structures, reactivity and inferred acylation transition states for 2'-amine substituted RNA

Ribose 2'-amine substitutions are broadly useful as structural probes in nucleic acids. In addition, structure-selective chemical reaction at 2'-amine groups is a robust technology for interrogating local nucleotide flexibility and conformational changes in RNA and DNA. We analyzed crystal structures for several RNA duplexes containing 2'-amino cytidine (C(N)) residues that form either C(N)-G base pairs or C(N)-A mismatches. The 2'-amine substitution is readily accommodated in an A-form RNA helix and thus differs from the C2'-endo conformation observed for free nucleosides. The 2'-amide product structure was visualized directly by acylating a C(N)-A mismatch in intact crystals and is also compatible with A-form geometry. To visualize conformations able to facilitate formation of the amide-forming transition state, in which the amine nucleophile carries a positive partial charge, we analyzed crystals of the C(N)-A duplex at pH 5, where the 2'-amine is protonated. The protonated amine moves to form a strong electrostatic interaction with the 3'-phosphodiester. Taken together with solution-phase experiments, 2'-amine acylation is likely facilitated by either of two transition states, both involving precise positioning of the adjacent 3'-phosphodiester group.     

B-DNA Helix Stability in a Solvent-Free Environment

B-DNA is the most common DNA helix conformation under physiological conditions. However, when the amount of water in a DNA solution is decreased, B-to-A helix transitions have been observed. To understand what type of helix conformations exist in a solvent-free environment, a series of poly d(CG)(n) and mixed sequence DNA duplexes from 18 to 30 bp were examined with circular dichroism (CD), ESI-MS, ion mobility, and molecular dynamics. From the CD spectra, it was observed that all sequences had B-form helices in solution. However, the solvent-free results were more complex. For the poly d(CG)(n) series, the 18 bp duplex had an A-form helix conformation, both A- and B-helices were present for the 22 bp duplex, and only B-helices were observed for the 26 and 30 bp duplexes. Since these sequences were all present as B-DNA in solution, the observed solvent-free structures illustrate that smaller helices with fewer base pairs convert to A-DNA more easily than larger helices in the absence of solvent. A similar trend was observed for the mixed sequence duplexes where both an A- and B-helix were present for the 18 bp duplex, while only B-helices occur for the larger 22, 26, and 30 bp duplexes. Since the solvent-free B-helices appear at smaller sizes for the mixed sequences than for the pure d(CG)(n) duplexes, the pure d(CG)(n) duplexes have a greater A-philicity.     

Pyrene aromatic arrays on RNA duplexes as helical templates

RNA oligomers having multiple (2 to 4) pyrenylmethyl substituents at the 2'-O-sugar residues were synthesized. UV-melting studies showed that the pyrene-modified RNAs could form duplexes with complementary RNA sequences without loss of thermal stability. Absorption, fluorescence, and circular dichroism (CD) spectra revealed that the incorporated pyrenes projected toward the outside of A-form RNA duplexes and assembled in helical aromatic arrays along the minor grooves of the RNA duplexes. Results of computer simulations agreed with the assembled structures of the pyrenes. The helical pyrene arrays exhibited remarkably strong excimer fluorescence, which was dependent on the sequence contexts of RNA duplexes.     

Evaluating nonpolarizable nucleic acid force fields: A systematic comparison of the nucleobases hydration free energies and chloroform‐to‐water partition coefficients

Nucleic acid force fields have been shown to reproduce structural properties of DNA and RNA very well, but comparative studies with respect to thermodynamic properties are rare. As a test for thermodynamic properties, we have computed hydration free energies and chloroform‐to‐water partition coefficients of nucleobases using the AMBER‐99, AMBER‐gaff, CHARMM‐27, GROMOS‐45a4/53a6 and OPLS‐AA force fields. A mutual force field comparison showed a very large spread in the calculated thermodynamic properties, demonstrating that some of the parameter sets require further optimization. The choice of solvent model used in the simulation does not have a significant effect on the results. Comparing the hydration free energies obtained by the various force fields to the adenine and thymine experimental values showed a very large deviation for the GROMOS and AMBER parameter sets. Validation against experimental partition coefficients showed good agreement for the CHARMM‐27 parameter set. In view of mutation studies, differences in partition coefficient between two bases were also compared, and good agreement between experiments and calculations was found for the AMBER‐99 parameter set. Overall, the CHARMM‐27 parameter set performs best with respect to the thermodynamic properties tested here. © 2012 Wiley Periodicals, Inc.     

Structural studies of LNA:RNA duplexes by NMR: conformations and implications for RNase H activity

We have used NMR and CD spectroscopy to study the conformations of modified oligonucleotides (locked nucleic acid, LNA) containing a conformationally restricted nucleotide (T(L)) with a 2'-O,4'-C-methylene bridge. We have investigated two LNA:RNA duplexes, d(CTGAT(L)ATGC):r(GCAUAUCAG) and d(CT(L)GAT(L)AT(L)GC):r(GCAUAUCAG), along with the unmodified DNA:RNA reference duplex. Increases in the melting temperatures of +9.6 degrees C and +8.1 degrees C per modification relative to the unmodified duplex were observed for these two LNA:RNA sequences. The three duplexes all adopt right-handed helix conformations and form normal Watson-Crick base pairs with all the bases in the anti conformation. Sugar conformations were determined from measurements of scalar coupling constants in the sugar rings and distance information derived from 1H-1H NOE measurements; all the sugars in the RNA strands of the three duplexes adopt an N-type conformation (A-type structure), whereas the sugars in the DNA strands change from an equilibrium between S- and N-type conformations in the unmodified duplex towards more of the N-type conformation when modified nucleotides are introduced. The presence of three modified T(L) nucleotides induces drastic conformational shifts of the remaining unmodified nucleotides of the DNA strand, changing all the sugar conformations except those of the terminal sugars to the N type. The CD spectra of the three duplexes confirm the structural changes described above. On the basis of the results reported herein, we suggest that the observed conformational changes can be used to tune LNA:RNA duplexes into substrates for RNase H: Partly modified LNA:RNA duplexes may adopt a duplex structure between the standard A and B types, thereby making the RNA strand amenable to RNase H-mediated degradation.     

2-Aminopurine: a probe of structural dynamics and charge transfer in DNA and DNA:RNA hybrids

Spectroscopic techniques are employed to probe relationships between structural dynamics and charge transfer (CT) efficiency in DNA duplexes and DNA:RNA hybrids containing photoexcited 2-aminopurine (Ap). To better understand the variety of interactions and reactions, including CT, between Ap and DNA, the fluorescence behavior of Ap is investigated in a full series of redox-inactive as well as redox-active assemblies. Thus, Ap is developed as a dual reporter of structural dynamics and base-base CT reactions in nucleic acid duplexes. CD, NMR, and thermal denaturation profiles are consistent with the family of DNA duplexes adopting a distinct conformation versus the DNA:RNA hybrids. Fluorescence measurements establish that the d(A)-r(U) tract of the DNA:RNA hybrid exhibits enhanced structural flexibility relative to that of the d(A)-d(T) tract of the DNA duplexes. The yield of CT from either G or 7-deazaguanine (Z) to Ap in the assemblies was determined by comparing Ap emission in redox-active G- or Z-containing duplexes to otherwise identical duplexes in which the G or Z is replaced by inosine (I), the redox-inactive nucleoside analogue. Investigations of CT not only demonstrate efficient intrastrand base-base CT in the DNA:RNA hybrids but also reveal a distance dependence of CT yield that is more shallow through the d(A)-r(U) bridge of the A-form DNA:RNA hybrids than through the d(A)-d(T) bridge of the B-form DNA duplexes. The shallow distance dependence of intrastrand CT in DNA:RNA hybrids correlates with the increased conformational flexibility of bases within the hybrid duplexes. Measurements of interstrand base-base CT provide another means to distinguish between the A- and B-form helices. Significantly, in the A-form DNA:RNA hybrids, a similar distance dependence is obtained for inter- and intrastrand reactions, while, in B-DNA, a more shallow distance dependence is evident with interstrand CT reactions. These observations are consistent with evaluations of intra- and interstrand base overlap in A- versus B-form duplexes. Overall, these data underscore the sensitivity of CT chemistry to nucleic acid structure and structural dynamics.     

C-F...H-C hydrogen bonds in ribonucleic acids

We report the synthesis of 1'-deoxy-1'-(benzimidazol-1-yl)-beta-D-ribofuranose 7 and 1'-deoxy-1'-phenyl-beta-D-ribofuranose 2. With these two ribonucleoside analogues we have a set of nine different RNA building blocks in hand, which are isostere to the natural bases. Now it is possible to investigate their duplex stabilizing forces. These forces are hydrogen bonds, base stacking, and solvation. The phosphoramidites of all building blocks were incorporated into a 12mer RNA, and the resulting RNA duplexes were investigated by UV- and CD-spectroscopy. We found that some of the RNA analogues are universal bases. The best universal bases with the lowest destabilization and the smallest discrimination between the natural bases are 1 (B) and 9 (E). On the basis of UV measurements we determined the melting points and the thermodynamic data. We were able to show that there are no hydrogen bonds between the natural bases and the RNA analogues. From thermodynamic data we calculated the contributions for base stacking and solvation of all modified building blocks. Comparison of calculated and measured data of double modified base pairs in 12mer RNA duplexes showed a further duplex stabilizing force in base pairs containing fluorine atoms at the Watson-Crick binding site. This stabilizing force can be defined as C-F.H-C hydrogen bond as is observed in crystal structures of 1'-deoxy-1'-(4-fluorophenyl)-beta-D-ribofuranose.     

Stacking effects on local structure in RNA: changes in the structure of tandem GA pairs when flanking GC pairs are replaced by isoG-isoC pairs

The Watson-Crick-like isoG-isoC (iGiC) pair, with the amino and carbonyl groups transposed relative to the Watson-Crick GC pair, provides an expanded alphabet for understanding interactions that shape nucleic acid structure. Here, thermodynamic stabilities of tandem GA pairs flanked by iGiC pairs are reported along with the NMR structures of the RNA self-complementary duplexes (GCiGGAiCGCA)2 and (GGiCGAiGCCA)2. A sheared GA pairing forms in (GCiGGAiCGCA)2, and an imino GA pairing forms in (GGiCGAiGCCA)2. The structures contrast with the formation of tandem imino and sheared GA pairs flanked by GC pairs in the RNA self-complementary duplexes (GCGGACGC)2 and (GGCGAGCC)2, respectively. In both iGiC duplexes, Watson-Crick-like hydrogen bonds are formed between iG and iC, and iGiC substitutions result in less favorable loop stability. The results provide benchmarks for testing computations of molecular interactions that shape RNA three-dimensional structure.     

NMR solution structure of an oligonucleotide hairpin with a 2'F-ANA/RNA stem: implications for RNase H specificity toward DNA/RNA hybrid duplexes.

The first structure of a 2'-deoxy-2'-fluoro-D-arabinose nucleic acid (2'F-ANA)/RNA duplex is presented. We report the structural characterization by NMR spectroscopy of a small hybrid hairpin, r(GGAC)d(TTCG)2'F-a(GTCC), containing a 2'F-ANA/RNA stem and a four-residue DNA loop. Complete (1)H, (13)C, (19)F, and (31)P resonance assignments, scalar coupling constants, and NOE constraints were obtained from homonuclear and heteronuclear 2D spectra. In the chimeric duplex, the RNA strand adopts a classic A-form structure having C3' endo sugar puckers. The 2'F-ANA strand is neither A-form nor B-form and contains O4' endo sugar puckers. This contrasts strongly with the dynamic sugar conformations previously observed in the DNA strands of DNA/RNA hybrid duplexes. Structural parameters for the duplex, such as minor groove width, x-displacement, and inclination, were intermediate between those of A-form and B-form duplexes and similar to those of DNA/RNA duplexes. These results rationalize the enhanced stability of 2'F-ANA/RNA duplexes and their ability to elicit RNase H activity. The results are relevant for the design of new antisense drugs based on sugar-modified nucleic acids.     

Neomycin-neomycin dimer: an all-carbohydrate scaffold with high affinity for AT-rich DNA duplexes

A dimeric neomycin-neomycin conjugate 3 with a flexible linker, 2,2'-(ethylenedioxy)bis(ethylamine), has been synthesized and characterized. Dimer 3 can selectively bind to AT-rich DNA duplexes with high affinity. Biophysical studies have been performed between 3 and different nucleic acids with varying base composition and conformation by using ITC (isothermal calorimetry), CD (circular dichroism), FID (fluorescent intercalator displacement), and UV (ultraviolet) thermal denaturation experiments. A few conclusions can be drawn from this study: (1) FID assay with 3 and polynucleotides demonstrates the preference of 3 toward AT-rich sequences over GC-rich sequences. (2) FID assay and UV thermal denaturation experiments show that 3 has a higher affinity for the poly(dA)·poly(dT) DNA duplex than for the poly(dA)·2poly(dT) DNA triplex. Contrary to neomycin, 3 destabilizes poly(dA)·2poly(dT) triplex but stabilizes poly(dA)·poly(dT) duplex, suggesting the major groove as the binding site. (3) UV thermal denaturation studies and ITC experiments show that 3 stabilizes continuous AT-tract DNA better than DNA duplexes with alternating AT bases. (4) CD and FID titration studies show a DNA binding site size of 10-12 base pairs/drug, depending upon the structure/sequence of the duplex for AT-rich DNA duplexes. (5) FID and ITC titration between 3 and an intramolecular DNA duplex [d(5'-A(12)-x-T(12)-3'), x = hexaethylene glycol linker] results in a binding stoichiometry of 1:1 with a binding constant ～10(8) M(-1) at 100 mM KCl. (6) FID assay using 3 and 512 hairpin DNA sequences that vary in their AT base content and placement also show a higher binding selectivity of 3 toward continuous AT-rich than toward DNA duplexes with alternate AT base pairs. (7) Salt-dependent studies indicate the formation of three ion pairs during binding of the DNA duplex d[5'-A(12)-x-T(12)-3'] and 3. (8) ITC-derived binding constants between 3 and DNA duplexes have the following order: AT continuous, d[5'-G(3)A(5)T(5)C(3)-3'] > AT alternate, d[5'-G(3)(AT)(5)C(3)-3'] > GC-rich d[5'-A(3)G(5)C(5)T(3)-3']. (9) 3 binds to the AT-tract-containing DNA duplex (B* DNA, d[5'-G(3)A(5)T(5)C(3)-3']) with 1 order of magnitude higher affinity than to a DNA duplex with alternating AT base pairs (B DNA, d[5'-G(3)(AT)(5)C(3)-3']) and with almost 3 orders of magnitude higher affinity than a GC-rich DNA (A-form, d[5'-A(3)G(5)C(5)T(3)-3']).     

Molecular Dynamisc Simulation of Polyelectrolites

We carried out simulations of a polymer chain using molecular dynamics algorythm. As a model we used a three-dimensional set monomers (electrically charged material points) connected with its nearest neighbours by harmonic potential. Additionally all pairs of segments interacts by the Lennard-Jones (LJ) and Coulomb forces. The aim of the simulation was to determine chain conformation and other basic properties like radius of gyration and moment of inertia for various polymer length and electric charge distribution.Presented model could be alternative tool for structure prediction to typically used ones based on AMBER 99 [1] or another advanced force field.     

On the accuracy of force fields for predicting the physical properties of dimethylnitramine

The accuracy of three force fields for predicting the physical properties of dimethylnitramine (DMNA) has been investigated by using molecular dynamics simulations. The Sorescu, Rice, and Thompson (SRT) (J. Phys. Chem. B 1997, 101, 798) rigid-molecule, flexible generalized AMBER (J. Comput. Chem. 2004, 25, 1157), and Smith et al. flexible force fields (J. Phys. Chem. B 1999, 103, 705) were tested. The density, lattice parameters, isotherm, and melting point of DMNA are calculated using classical molecular dynamics. Except for the melting point, the predictions of the three force fields are in reasonable agreement with experimental values. The calculated thermodynamic melting points (Tmp) for the SRT, AMBER, and Smith et al. force fields are 380, 360, and 260 K, respectively. The experimental value is 331 K. Modifications of the torsional barriers in the AMBER force field resulted in Tmp = 346 K, in good agreement with the experimental value of 331 K. The calculated lattice parameters and bulk modulus are also improved with the modifications of the AMBER potential. The results indicate that, although not sufficiently accurate without modifications, the general force fields such as AMBER provide the basis for developing force fields that correctly predict the physical properties of nitramines.     

Structural rationalization of a large difference in RNA affinity despite a small difference in chemistry between two 2'-O-modified nucleic acid analogues

Chemical modification of nucleic acids at the 2'-position of ribose has generated antisense oligonucleotides (AONs) with a range of desirable properties. Electron-withdrawing substituents such as 2'-O-[2-(methoxy)ethyl] (MOE) confer enhanced RNA affinity relative to that of DNA by conformationally preorganizing an AON for pairing with the RNA target and by improving backbone hydration. 2'-Substitution of the ribose has also been shown to increase nuclease resistance and cellular uptake via changes in lipophilicity. Interestingly, incorporation of either 2'-O-[2-(methylamino)-2-oxoethyl]- (NMA) or 2'-O-(N-methylcarbamate)-modified (NMC) residues into AONs has divergent effects on RNA affinity. Incorporation of 2'-O-NMA-T considerably improves RNA affinity while incorporation of 2'-O-NMC-T drastically reduces RNA affinity. Crystal structures at high resolution of A-form DNA duplexes containing either 2'-O-NMA-T or 2'-O-NMC-T shed light on the structural origins of the surprisingly large difference in stability given the relatively minor difference in chemistry between NMA and NMC. NMA substituents adopt an extended conformation and use either their carbonyl oxygen or amino nitrogen to trap water molecules between phosphate group and sugar. The conformational properties of NMA and the observed hydration patterns are reminiscent of those found in the structures of 2'-O-MOE-modified RNA. Conversely, the carbonyl oxygen of NMC and O2 of T are in close contact, providing evidence that an unfavorable electrostatic interaction and the absence of a stable water structure are the main reasons for the loss in thermodynamic stability as a result of incorporation of 2'-O-NMC-modified residues.     

Molecular dynamics simulation of the A-DNA to B-DNA transition in aqueous RbCl solution

Unrestrained molecular dynamics (MD) simulations have been carried out to characterize the stability of DNA conformations and the dynamics of A-DNA→B-DNA conformational transitions in aqueous RbCl solutions. The PARM99 force field in the AMBER8 package was used to investigate the effect of RbCl concentration on the dynamics of the A→B conformational transition in the DNA duplex d(CGCGAATTCGCG)2 . Canonical Aand B-form DNA were assumed for the initial conformation and the final conformation had a length per complete turn that matched the canonical B-DNA. The DNA structure was monitored for 3.0 ns and the distances between the C5′ atoms were obtained from the simulations. It was found that all of the double stranded DNA strands of A-DNA converged to the structure of B-form DNA within 1.0 ns during the unrestrained MD simulations. In addition, increasing the RbCl concentration in aqueous solution hindered the A→B conformational transition and the transition in aqueous RbCl solution was faster than that in aqueous NaCl solution for the same electrolyte strength. The effects of the types and concentrations of counterions on the dynamics of the A→B conformational transition can be understood in terms of the variation in water activity and the number of accumulated counterions in the major grooves of A-DNA. The rubidium ion distributions around both fixed A-DNA and B-DNA were obtained using the restrained MD simulations to help explain the effect of RbCl concentration on the dynamics of the A→B conformational transition.     

Study of Tamiflu sensitivity to variants of A/H5N1 virus using different force fields

An accurate estimation of binding free energy of a ligand to receptor ΔG(bind) is one of the most important problems in drug design. The success of solution of this problem is expected to depend on force fields used for modeling a ligand-receptor complex. In this paper, we consider the impact of four main force fields, AMBER99SB, CHARMM27, GROMOS96 43a1, and OPLS-AA/L, on the binding affinity of Oseltamivir carboxylate to the wild-type and Y252H, N294S, and H274Y mutants of glycoprotein neuraminidase from the pandemic A/H5N1 virus. Having used the molecular mechanic-Poisson-Boltzmann surface area method, we have shown that ΔG(bind), obtained by AMBER99SB, OPLS-AA/L, and CHARMM27, shows the high correlation with the available experimental data. They correctly capture the binding ranking Y252H → WT → N294S → H274Y observed in experiments (Collins, P. J. et al. Nature 2008, 453, 1258). In terms of absolute values of binding scores, results obtained by AMBER99SB are in the nearest range with experiments, while OPLS-AA/L, which is applied to study binding of Oseltamivir to the influenza virus for the first time, gives rather big negative values for ΔG(bind). GROMOS96 43a1 provides a lower correlation as it supports Oseltamivir to be more resistant to N294S than H274Y. Our study suggests that force fields have pronounced influence on theoretical estimations of binding free energy of a ligand to receptor. The effect of all-atom models on dynamics of the binding pocket as well as on the hydrogen-bond network between Oseltamivir and receptors is studied in detail. The hydrogen network, obtained by GROMOS, is weakest among four studied force fields.     

Solution structure of a HNA-RNA hybrid

BACKGROUND: Synthetic nucleic acid analogues with a conformationally restricted sugar-phosphate backbone are widely used in antisense strategies for biomedical and biochemical applications. The modified backbone protects the oligonucleotides against degradation within the living cell, which allows them to form stable duplexes with sequences in target mRNAs with the aim of arresting their translation. The biologically most active antisense oligonucleotides also trigger cleavage of the target RNA through activation of endogenous RNase H. Systematic studies of synthetic oligonucleotides have also been conducted to delineate the origin of the chirality of DNA and RNA that are both composed of D-nucleosides. RESULTS: Hexitol nucleic acids (HNA) are the first example of oligonucleotides with a six-membered carbohydrate moiety that can bind strongly and selectively to complementary RNA oligomers. We present the first high resolution nuclear magnetic resonance structure of a HNA oligomer bound to a complementary RNA strand. The HNA-RNA complex forms an anti-parallel heteroduplex and adopts a helical conformation that belongs to the A-type family. Possibly, due to the rigidity of the rigid chair conformation of the six-membered ring both the HNA and RNA strand in the duplex are well defined. The observed absence of end-fraying effects also indicate a reduced conformational flexibility of the HNA-RNA duplex compared to canonical dsRNA or an RNA-DNA duplex. CONCLUSIONS: The P-P distance across the minor groove, which is close to A-form, and the rigid conformation of the HNA-RNA complex, explain its resistance towards degradation by Rnase H. The A-form character of the HNA-RNA duplex and the reduced flexibility of the HNA strand is possibly responsible for the stereoselectivity of HNA templates in non-enzymatic replication of oligonucleotides, supporting the theory that nucleosides with six-membered rings could have existed at some stage in molecular evolution.     

Intrinsic contribution of the 2'-hydroxyl to RNA conformational heterogeneity

Canonical duplex RNA assumes only the A-form conformation at the secondary structure level while, in contrast, a wide range of noncanonical, tertiary conformations of RNA occur. Here, we show how the 2'-hydroxyl controls RNA conformational properties. Quantum mechanical calculations reveal that the orientation of the 2'-hydroxyl significantly alters the intrinsic flexibility of the phosphodiester backbone, favoring the A-form in duplex RNA when it is in the base orientation and facilitating sampling of a wide range of noncanonical, tertiary structures when it is in the O3' orientation. Influencing the orientation of the 2'-hydroxyl are interactions with the environment, as evidenced by crystallographic survey data, indicating the 2'-hydroxyl to sample more of the O3' orientation in noncanonical RNA structures. These results indicate that the 2'-hydroxyl acts as a "switch", both limiting the conformation of RNA to the A-form at the secondary structure level and allowing RNA to sample a wide range of noncanonical tertiary conformations.