Hydrophobic Groups Stabilize the Hydration Shell of 2'-O-Methylated RNA Duplexes

Regulation of the epigenome through RNA modifications

Chemical modifications of nucleotides expand the complexity and functional properties of genomes and transcriptomes. A handful of modifications in DNA bases are part of the epigenome, wherein DNA methylation regulates chromatin structure, transcription, and co-transcriptional RNA processing. In contrast, more than 150 chemical modifications of RNA constitute the epitranscriptome. Ribonucleoside modifications comprise a diverse repertoire of chemical groups, including methylation, acetylation, deamination, isomerization, and oxidation. Such RNA modifications regulate all steps of RNA metabolism, including folding, processing, stability, transport, translation, and RNA's intermolecular interactions. Initially thought to influence all aspects of the post-transcriptional regulation of gene expression exclusively, recent findings uncovered a crosstalk between the epitranscriptome and the epigenome. In other words, RNA modifications feedback to the epigenome to transcriptionally regulate gene expression. The epitranscriptome achieves this feat by directly or indirectly affecting chromatin structure and nuclear organization. This review highlights how chemical modifications in chromatin-associated RNAs (caRNAs) and messenger RNAs (mRNAs) encoding factors involved in transcription, chromatin structure, histone modifications, and nuclear organization affect gene expression transcriptionally.

Influence of Sugar Modifications on the Nucleoside Conformation and Oligonucleotide Stability: A Critical Review

Ribofuranose sugar conformation plays an important role in the structure and dynamics of functional nucleic acids such as siRNAs, AONs, aptamers, miRNAs, etc. To improve their therapeutic potential, several chemical modifications have been introduced into the sugar moiety over the years. The stability of the oligonucleotide duplexes as well as the formation of stable and functional protein-oligonucleotide complexes are dictated by the conformation and dynamics of the sugar moiety. In this review, we systematically categorise various ribofuranose sugar modifications employed in DNAs and RNAs so far. We discuss different stereoelectronic effects imparted by different substituents on the sugar ring and how these effects control sugar puckering. Using this data, it would be possible to predict the precise use of chemical modifications and design novel sugar-modified nucleosides for therapeutic oligonucleotides that can improve their physicochemical properties.

Roles of non-canonical structures of nucleic acids in cancer and neurodegenerative diseases

Cancer and neurodegenerative diseases are caused by genetic and environmental factors. Expression of tumour suppressor genes is suppressed by mutations or epigenetic silencing, whereas for neurodegenerative disease-related genes, nucleic acid-based effects may be presented through loss of protein function due to erroneous protein sequences or gain of toxic function from extended repeat transcripts or toxic peptide production. These diseases are triggered by damaged genes and proteins due to lifestyle and exposure to radiation. Recent studies have indicated that transient, non-canonical structural changes in nucleic acids in response to the environment can regulate the expression of disease-related genes. Non-canonical structures are involved in many cellular functions, such as regulation of gene expression through transcription and translation, epigenetic regulation of chromatin, and DNA recombination. Transcripts generated from repeat sequences of neurodegenerative disease-related genes form non-canonical structures that are involved in protein transport and toxic aggregate formation. Intracellular phase separation promotes transcription and protein assembly, which are controlled by the nucleic acid structure and can influence cancer and neurodegenerative disease progression. These findings may aid in elucidating the underlying disease mechanisms. Here, we review the influence of non-canonical nucleic acid structures in disease-related genes on disease onset and progression.

New Insights into the Functions of Nucleic Acids Controlled by Cellular Microenvironments

The right-handed double-helical B-form structure (B-form duplex) has been widely recognized as the canonical structure of nucleic acids since it was first proposed by James Watson and Francis Crick in 1953. This B-form duplex model has a monochronic and static structure and codes genetic information within a sequence. Interestingly, DNA and RNA can form various non-canonical structures, such as hairpin loops, left-handed helices, triplexes, tetraplexes of G-quadruplex and i-motif, and branched junctions, in addition to the canonical structure. The formation of non-canonical structures depends not only on sequence but also on the surrounding environment. Importantly, these non-canonical structures may exhibit a wide variety of biological roles by changing their structures and stabilities in response to the surrounding environments, which undergo vast changes at specific locations and at specific times in cells. Here, we review recent progress regarding the interesting behaviors and functions of nucleic acids controlled by molecularly crowded cellular conditions. New insights gained from recent studies suggest that nucleic acids not only code genetic information in sequences but also have unknown functions regarding their structures and stabilities through drastic structural changes in cellular environments.

The 'hidden side' of spin labelled oligonucleotides: Molecular dynamics study focusing on the EPR-silent components of base pairing

Nitroxide labels are combined with nucleic acid structures and are studied using electron paramagnetic resonance experiments (EPR). As X-ray/NMR structures are unavailable with the nitroxide labels, detailed residue level information, down to atomic resolution, about the effect of these nitroxide labels on local RNA structures is currently lacking. This information is critical to evaluate the choice of spin label. In this study, we compare and contrast the effect of TEMPO-based (N^T) and rigid spin (Ç) labels (in both 2'-O methylated and not-methylated forms) on RNA duplexes. We also investigate sequence- dependent effects of N^T label on RNA duplex along with the more complex G-quadruplex RNA. Distances measured from molecular dynamics simulations between the two spin labels are in agreement with the EPR experimental data. To understand the effect of labelled oligonucleotides on the structure, we studied the local base pair geometries and global structure in comparison with the unlabelled structures. Based on the structural analysis, we can conclude that TEMPO-based and Ç labels do not significantly perturb the base pair arrangements of the native oligonucleotide. When experimental structures for the spin labelled DNA/RNA molecules are not available, general framework offered by the current study can be used to provide information critical to the choice of spin labels to facilitate future EPR studies.

2'-O-Methylation can increase the abundance and lifetime of alternative RNA conformational states

2'-O-Methyl (Nm) is a highly abundant post-transcriptional RNA modification that plays important biological roles through mechanisms that are not entirely understood. There is evidence that Nm can alter the biological activities of RNAs by biasing the ribose sugar pucker equilibrium toward the C3'-endo conformation formed in canonical duplexes. However, little is known about how Nm might more broadly alter the dynamic ensembles of flexible RNAs containing bulges and internal loops. Here, using NMR and the HIV-1 transactivation response (TAR) element as a model system, we show that Nm preferentially stabilizes alternative secondary structures in which the Nm-modified nucleotides are paired, increasing both the abundance and lifetime of low-populated short-lived excited states by up to 10-fold. The extent of stabilization increased with number of Nm modifications and was also dependent on Mg2+. Through phi-value analysis, the Nm modification also provided rare insights into the structure of the transition state for conformational exchange. Our results suggest that Nm could alter the biological activities of Nm-modified RNAs by modulating their secondary structural ensembles as well as establish the utility of Nm as a tool for the discovery and characterization of RNA excited state conformations.

Comparative patterns of modified nucleotides in individual tRNA species from a mesophilic and two thermophilic archaea

To improve and complete our knowledge of archaeal tRNA modification patterns, we have identified and compared the modification pattern (type and location) in tRNAs of three very different archaeal species, Methanococcus maripaludis (a mesophilic methanogen), Pyrococcus furiosus (a hyperthermophile thermococcale), and Sulfolobus acidocaldarius (an acidophilic thermophilic sulfolobale). Most abundant isoacceptor tRNAs (79 in total) for each of the 20 amino acids were isolated by two-dimensional gel electrophoresis followed by in-gel RNase digestions. The resulting oligonucleotide fragments were separated by nanoLC and their nucleotide content analyzed by mass spectrometry (MS/MS). Analysis of total modified nucleosides obtained from complete digestion of bulk tRNAs was also performed. Distinct base- and/or ribose-methylations, cytidine acetylations, and thiolated pyrimidines were identified, some at new positions in tRNAs. Novel, some tentatively identified, modifications were also found. The least diversified modification landscape is observed in the mesophilic Methanococcus maripaludis and the most complex one in Sulfolobus acidocaldarius Notable observations are the frequent occurrence of ac⁴C nucleotides in thermophilic archaeal tRNAs, the presence of m⁷G at positions 1 and 10 in Pyrococcus furiosus tRNAs, and the use of wyosine derivatives at position 37 of tRNAs, especially those decoding U1- and C1-starting codons. These results complete those already obtained by others with sets of archaeal tRNAs from Methanocaldococcus jannaschii and Haloferax volcanii.

Aberrant regulation of RNA methylation during spermatogenesis

Natural modifications of cellular RNA include various chemical modifications, such as N6-methyladenosine (m⁶ A), which enable the orderly metabolism and function of RNA structural diversity, thereby affecting gene expression. Spermatogenesis is a complex differentiating developmental process, which includes the proliferation of spermatogonial stem cells, spermatocyte meiosis and sperm maturation. Emerging evidence has shown that RNA methylation can influence RNA splicing, exportation and translation, which are controlled in the male germline in order to ensure coordinated gene expression. In this review, we summarize the typical characteristics of different types of RNA methylation during the process of spermatogenesis. In particular, we emphasize the functions of the RNA methylation effectors during the male germ cell development.

Chemical biology of non-canonical structures of nucleic acids for therapeutic applications

DNA forms not only the canonical duplex structure but also non-canonical structures. Most potential sequences that induce the formation of non-canonical structures are present in disease-related genes. Interestingly, biological reactions are inhibited or dysregulated by non-canonical structure formation in disease-related genes. To control biological reactions, methods for inducing the formation of non-canonical structures have been developed using small molecules and oligonucleotides. In this feature article, we review biological reactions such as replication, transcription, and reverse transcription controlled by non-canonical DNA structures formed by disease-related genes. Furthermore, we discuss recent studies aimed at developing methods for regulating these biological reactions using drugs targeting the DNA structure.

Same fold, different properties: polarizable molecular dynamics simulations of telomeric and TERRA G-quadruplexes

DNA and RNA sequences rich in guanine can fold into noncanonical structures called G-quadruplexes (GQs), which exhibit a common stem structure of Hoogsteen hydrogen-bonded guanine tetrads and diverse loop structures. GQ sequence motifs are overrepresented in promoters, origins of replication, telomeres, and untranslated regions in mRNA, suggesting roles in modulating gene expression and preserving genomic integrity. Given these roles and unique aspects of different structures, GQs are attractive targets for drug design, but greater insight into GQ folding pathways and the interactions stabilizing them is required. Here, we performed molecular dynamics simulations to study two bimolecular GQs, a telomeric DNA GQ and the analogous telomeric repeat-containing RNA (TERRA) GQ. We applied the Drude polarizable force field, which we show outperforms the additive CHARMM36 force field in both ion retention and maintenance of the GQ folds. The polarizable simulations reveal that the GQs bind bulk K+ ions differently, and that the TERRA GQ accumulates more K+ ions, suggesting different ion interactions stabilize these structures. Nucleobase dipole moments vary as a function of position and also contribute to ion binding. Finally, we show that the TERRA GQ is more sensitive than the telomeric DNA GQ to water-mediated modulation of ion-induced dipole-dipole interactions.

Influence of nucleotide modifications at the C2' position on the Hoogsteen base-paired parallel-stranded duplex of poly(A) RNA

Polyadenylate (poly(A)) has the ability to form a parallel duplex with Hoogsteen adenine:adenine base pairs at low pH or in the presence of ammonium ions. In order to evaluate the potential of this structural motif for nucleic acid-based nanodevices, we characterized the effects on duplex stability of substitutions of the ribose sugar with 2'-deoxyribose, 2'-O-methyl-ribose, 2'-deoxy-2'-fluoro-ribose, arabinose and 2'-deoxy-2'-fluoro-arabinose. Deoxyribose substitutions destabilized the poly(A) duplex both at low pH and in the presence of ammonium ions: no duplex formation could be detected with poly(A) DNA oligomers. Other sugar C2' modifications gave a variety of effects. Arabinose and 2'-deoxy-2'-fluoro-arabinose nucleotides strongly destabilized poly(A) duplex formation. In contrast, 2'-O-methyl and 2'-deoxy-2'-fluoro-ribo modifications were stabilizing either at pH 4 or in the presence of ammonium ions. The differential effect suggests they could be used to design molecules selectively responsive to pH or ammonium ions. To understand the destabilization by deoxyribose, we determined the structures of poly(A) duplexes with a single DNA residue by nuclear magnetic resonance spectroscopy and X-ray crystallography. The structures revealed minor structural perturbations suggesting that the combination of sugar pucker propensity, hydrogen bonding, pKa shifts and changes in hydration determine duplex stability.

Additive CHARMM force field for naturally occurring modified ribonucleotides

More than 100 naturally occurring modified nucleotides have been found in RNA molecules, in particular in tRNAs. We have determined molecular mechanics force field parameters compatible with the CHARMM36 all‐atom additive force field for all these modifications using the CHARMM force field parametrization strategy. Emphasis was placed on fine tuning of the partial atomic charges and torsion angle parameters. Quantum mechanics calculations on model compounds provided the initial set of target data, and extensive molecular dynamics simulations of nucleotides and oligonucleotides in aqueous solutions were used for further refinement against experimental data. The presented parameters will allow for computational studies of a wide range of RNAs containing modified nucleotides, including the ribosome and transfer RNAs. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Synthesis and properties of oligonucleotides modified with 2'-O-(2-carboxyethyl)nucleotides and their carbamoyl derivatives

2'-O-Methyl oligoribonucleotides with four kinds of 2'-O-modified uridine derivatives were synthesised. Their duplex stability, hydration behavior and exonuclease resistance were studied by spectroscopic analyses and molecular dynamics simulations. Consequently, 2'-O-modification of the uridine residue with 2-carbamoylethyl or 2-(N-methylcarbamoyl)ethyl groups resulted in a significant improvement of the exonuclease resistance without the loss of duplex stability.

A new modified cytosine base capable of base pairing with guanine using four hydrogen bonds

Oligonucleotides, containing 4-N-(1H-pyrrol-2-ylcarbonyl)deoxycytidine (dC(Pyc)) and related derivatives, were synthesized via deprotection using 1.5 M NaOMe/MeOH. Among them, oligodeoxynucleotides containing dC(Pyc) exhibited a higher hybridization affinity for DNA and RNA than the unmodified oligodeoxynucleotides. Comparative analysis between dC(Pyc) and its derivatives by molecular dynamic simulation indicated that the C(Pyc) residue could form four hydrogen bonds with the opposite G nucleobase keeping a more planar structure than the C(Inc) residue where the Pyc group was replaced with a 1H-indol-2-ylcarbonyl group.

Structure, stability and behaviour of nucleic acids in ionic liquids

Nucleic acids have become a powerful tool in nanotechnology because of their conformational polymorphism. However, lack of a medium in which nucleic acid structures exhibit long-term stability has been a bottleneck. Ionic liquids (ILs) are potential solvents in the nanotechnology field. Hydrated ILs, such as choline dihydrogen phosphate (choline dhp) and deep eutectic solvent (DES) prepared from choline chloride and urea, are 'green' solvents that ensure long-term stability of biomolecules. An understanding of the behaviour of nucleic acids in hydrated ILs is necessary for developing DNA materials. We here review current knowledge about the structures and stabilities of nucleic acids in choline dhp and DES. Interestingly, in choline dhp, A-T base pairs are more stable than G-C base pairs, the reverse of the situation in buffered NaCl solution. Moreover, DNA triplex formation is markedly stabilized in hydrated ILs compared with aqueous solution. In choline dhp, the stability of Hoogsteen base pairs is comparable to that of Watson-Crick base pairs. Moreover, the parallel form of the G-quadruplex is stabilized in DES compared with aqueous solution. The behaviours of various DNA molecules in ILs detailed here should be useful for designing oligonucleotides for the development of nanomaterials and nanodevices.

Control of stability and structure of nucleic acids using cosolutes

The stabilities, structures, and functions of nucleic acids are responsive to surrounding conditions. Living cells contain biomolecules, including nucleic acids, proteins, polysaccharides, and other soluble and insoluble low-molecular weight components, that occupy a significant fraction of the cellular volume (up to 40%), resulting in a highly crowded intracellular environment. We have proven that conditions that mimic features of this intra-cellular environment alter the physical properties affect the stability, structure, and function of nucleic acids. The ability to control structure of nucleic acids by mimicking intra-cellular conditions will be useful in nanotechnology applications of nucleic acids. This paper describes methods that can be used to analyze quantitatively the intra-cellular environment effects caused by cosolutes on nucleic acid structures and to regulate properties of nucleic acids using cosolutes.

Synthesis, gene silencing, and molecular modeling studies of 4'-C-aminomethyl-2'-O-methyl modified small interfering RNAs

The linear syntheses of 4'-C-aminomethyl-2'-O-methyl uridine and cytidine nucleoside phosphoramidites were achieved using glucose as the starting material. The modified RNA building blocks were incorporated into small interfering RNAs (siRNAs) by employing solid phase RNA synthesis. Thermal melting studies showed that the modified siRNA duplexes exhibited slightly lower T(m) (∼1 °C/modification) compared to the unmodified duplex. Molecular dynamics simulations revealed that the 4'-C-aminomethyl-2'-O-methyl modified nucleotides adopt South-type conformation in a siRNA duplex, thereby altering the stacking and hydrogen-bonding interactions. These modified siRNAs were also evaluated for their gene silencing efficiency in HeLa cells using a luciferase-based reporter assay. The results indicate that the modifications are well tolerated in various positions of the passenger strand and at the 3' end of the guide strand but are less tolerated in the seed region of the guide strand. The modified siRNAs exhibited prolonged stability in human serum compared to unmodified siRNA. This work has implications for the use of 4'-C-aminomethyl-2'-O-methyl modified nucleotides to overcome some of the challenges associated with the therapeutic utilities of siRNAs.

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.

Short-RNA selective binding of oligonucleotides modified using adenosine and guanosine derivatives that possess cyclohexyl phosphates as substituents

We have developed new artificial oligonucleotides which distinguish short RNA targets from long ones. The modification of the 5' termini of oligonucleotides by using adenosine derivatives that possess a bulky cyclohexyl phosphate moiety at their base moiety and a phosphate group at the position of their 5'-hydroxyl group maximized their short RNA selectivity. The 2'-O-methyl-RNA (5'-XC(m)A(m)A(m)C(m)C(m)U(m)A(m)C(m)U(m)) having these modifications exhibits ca. 10 °C higher T(m) in the duplexes with the complementary short RNA (3'-GUUGGAUGA-5') than with the long RNA (3'-AUUAUAUGUUGGAUGAUGGUUA-5'). The oligodeoxynucleotides having the same modification exhibited similar selectivity. Such short-RNA selective binding of terminally modified oligonucleotides can be employed to distinguish between mature microRNAs and pre-microRNAs.

Synthesis and hybridization properties of 2'-O-methylated oligoribonucleotides incorporating 2'-O-naphthyluridines

2'-O-(1-Naphthyl)uridine and 2'-O-(2-naphthyl)uridine were synthesized by a microwave-mediated reaction of 2,2'-anhydrouridine with naphthols. Using the 3'-phosphoramidite building blocks, these 2'-O-aryluridine derivatives were incorporated into 2'-O-methylated oligoribonucleotides. Incorporation of five 2'-O-(2-naphthyl)uridines into a 2'-O-methylated RNA sense strand significantly increased the thermostability of the duplex with a 2'-O-methylated RNA antisense strand. Circular dichroism spectroscopy and molecular dynamic simulation of the duplexes formed between the modified RNAs and 2'-O-methyl RNAs suggested that there are π-π interactions between two neighboring naphthyl groups in a sequence of the five consecutively modified nucleosides.

Linear relationship between deformability and thermal stability of 2'-O-modified RNA hetero duplexes

We describe the relationship between the experimentally determined melting temperatures of 2′-O-modified-RNA/RNA duplexes and their deformability estimated from molecular dynamics simulations. To clarify this relationship, we synthesized several fully modified oligoribonucleotides such as 2′-O-cyanoethyl RNAs and 2′-O-methoxyethyl RNAs and compared the actual melting temperatures of the duplexes with their calculated deformabilities. An increase of the melting temperatures by 2′-O-modifications was found to correlate strongly with an increase of the helical elastic constants in U₁₄/A₁₄, (CU)₇/(AG)₇, and (GACU)₃/(AGUC)₃ sequences. Linear regression analyses could be used to estimate the melting temperature with an accuracy of ±2.0 °C in our model case. Although the strong correlation was observed in the same base sequence, the linear regression functions were different from each base sequence. Our results indicated the possibility of predicting the thermal stability of 2′-O-modified duplexes at the computer-aided molecular design stage.

Crystallographic studies of chemically modified nucleic acids: a backward glance

Chemically modified nucleic acids (CNAs) are widely explored as antisense oligonucleotide or small interfering RNA (siRNA) candidates for therapeutic applications. CNAs are also of interest in diagnostics, high-throughput genomics and target validation, nanotechnology and as model systems in investigations directed at a better understanding of the etiology of nucleic acid structure, as well as the physicochemical and pairing properties of DNA and RNA, and for probing protein-nucleic acid interactions. In this article, we review research conducted in our laboratory over the past two decades with a focus on crystal-structure analyses of CNAs and artificial pairing systems. We highlight key insights into issues ranging from conformational distortions as a consequence of modification to the modulation of pairing strength, and RNA affinity by stereoelectronic effects and hydration. Although crystal structures have only been determined for a subset of the large number of modifications that were synthesized and analyzed in the oligonucleotide context to date, they have yielded guiding principles for the design of new analogs with tailor-made properties, including pairing specificity, nuclease resistance, and cellular uptake. And, perhaps less obviously, crystallographic studies of CNAs and synthetic pairing systems have shed light on fundamental aspects of DNA and RNA structure and function that would not have been disclosed by investigations solely focused on the natural nucleic acids.

Synthesis and hybridization of 2'-O-methyl-RNAs incorporating 2'-O-carbamoyluridine and unique participation of the carbamoyl group in U-G base pair

2'-O-Carbamoyluridine (U(cm)) was synthesized and incorporated into DNAs and 2'-O-Me-RNAs. The oligonucleotides incorporating U(cm) formed less stable duplexes with their complementary and U(cm)-U, U(cm)-C single-base mismatched DNAs and RNAs in comparison with those without the carbamoyl group. On the contrary, the T(m) analyses revealed that the duplexes with a mismatched U(cm)-G base pair showed almost the same thermostability as the corresponding unmodified duplexes. Molecular dynamics (MD) simulations of the U(cm)-modified 2'-O-Me-RNA/RNA duplexes with U(cm)-G mismatched base pair suggested that the carbamoyl group could participate in the U(cm)-G base pair by an additional intermolecular hydrogen bond between the carbamoyl oxygen and the H2 of the guanine base.

Fine tuning of electrostatics around the internucleotidic phosphate through incorporation of modified 2',4'-carbocyclic-LNAs and -ENAs leads to significant modulation of antisense properties

In the antisense (AS) and RNA interference (RNAi) technologies, the native single-stranded 2'-deoxyoligonucleotides (for AS) or double-stranded RNA (for RNAi) are chemically modified to bind to the target RNA in order to give improved downregulation of gene expression through inhibition of RNA translation. It is shown here how the fine adjustment of the electrostatic interaction by alteration of the substituents as well as their stereochemical environment around the internucleotidic phosphodiester moiety near the edge of the minor grove of the antisense oligonucleotides (AON)-RNA heteroduplex can lead to the modulation of the antisense properties. This was demonstrated through the synthesis of various modified carbocyclic-locked nucleic acids (LNAs) and -ethylene-bridged nucleic acids (ENAs) with hydroxyl and/or methyl substituents attached at the carbocyclic part and their integration into AONs by solid-phase DNA synthesis. The target affinity toward the complementary RNA and DNA, nuclease resistance, and RNase H elicitation by these modified AONs showed that both the nature of the modification (-OH versus -CH(3)) and their respective stereochemical orientations vis-a-vis vicinal phosphate play a very important role in modulating the AON properties. Whereas the affinity to the target RNA and the enzymatic stability of AONs were not favored by the hydrophobic and sterically bulky modifications in the center of the minor groove, their positioning at the edge of the minor groove near the phosphate linkage resulted in significantly improved nuclease resistance without loss of target affinity. On the other hand, hydrophilic modification, such as a hydroxyl group, close to the phosphate linkage made the internucleotidic phosphodiester especially nucleolytically unstable, and hence was not recommended. The substitutions on the carbocyclic moiety of the carba-LNA and -ENA did not affect significantly the choice of the cleavage sites of RNase H mediated RNA cleavage in the AON/RNA hybrid duplex, but the cleavage rate depended on the modification site in the AON sequence. If the original preferred cleavage site by RNase H was included in the 4-5nt stretch from the 3'-end of the modification site in the AON, decreased cleavage rate was observed. Upon screening of 52 modified AONs, containing 13 differently modified derivatives at C6' and C7' (or C8') of the carba-LNAs and -ENAs, two excellent modifications in the carba-LNA series were identified, which synergistically gave outstanding antisense properties such as the target RNA affinity, nuclease resistance, and RNase H activity and were deemed to be ideal candidates as potential antisense or siRNA therapeutic agents.

Synthesis and characterization of 2'-modified-4'-thioRNA: a comprehensive comparison of nuclease stability

We report herein the synthesis and physical and physiological characterization of fully modified 2′-modified-4′-thioRNAs, i.e. 2′-fluoro-4′-thioRNA (F-SRNA) and 2′-O-Me-4′-thioRNA (Me-SRNA), which can be considered as a hybrid chemical modification based on 2′-modified oligonucleotides (ONs) and 4′-thioRNA (SRNA). In its hybridization with a complementary RNA, F-SRNA (15mer) showed the highest T_m value (+16°C relative to the natural RNA duplex). In addition, both F-SRNA and Me-SRNA preferred RNA as a complementary partner rather than DNA in duplex formation. The results of a comprehensive comparison of nuclease stability of single-stranded F-SRNA and Me-SRNA along with 2′-fluoroRNA (FRNA), 2′-O-MeRNA (MeRNA), SRNA, and natural RNA and DNA, revealed that Me-SRNA had the highest stability with t_1/2 values of > 24 h against S1 nuclease (an endonuclease) and 79.2 min against SVPD (a 3′-exonuclease). Moreover, the stability of Me-SRNA was significantly improved in 50% human plasma (t_1/2 = 1631 min) compared with FRNA (t_1/2 = 53.2 min) and MeRNA (t_1/2 = 187 min), whose modifications are currently used as components of therapeutic aptamers. The results presented in this article will, it is hoped, contribute to the development of 2′-modified-4′-thioRNAs, especially Me-SRNA, as a new RNA molecule for therapeutic applications.

Analysis of the stability and flexibility of RNA complexes containing bulge loops of different sizes

Molecular dynamics simulations of RNA molecules consisting of an antisense oligonucleotide forming a complex with a target strand thereby creating an internal bulge-loop with 3, 4, or 5 nucleotides have been performed with and without O2' methylation of the antisense strand. The methylation influences minor groove hydration, in particular near guanines but also around the methylated O2', and it also reduces the flexibility of both RNA strands. A G.U wobble pair adjacent to the bulge-loop is also found to increase the flexibility of the bulge nucleotides, compared to the situation with a standard Watson-Crick G-C base-pair in the same position.

Ribozyme catalysis revisited: is water involved?

Enzymatic catalysis by RNA was discovered 25 years ago, yet mechanistic insights are emerging only slowly. Thought to be metalloenzymes at first, some ribozymes proved more versatile than anticipated when shown to utilize their own functional groups for catalysis. Recent evidence suggests that some may also judiciously place structural water molecules to shuttle protons in acid-base catalyzed reactions.

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.

Insights from crystallographic studies into the structural and pairing properties of nucleic acid analogs and chemically modified DNA and RNA oligonucleotides

Chemically modified nucleic acids function as model systems for native DNA and RNA; as chemical probes in diagnostics or the analysis of protein-nucleic acid interactions and in high-throughput genomics and drug target validation; as potential antigene-, antisense-, or RNAi-based drugs; and as tools for structure determination (i.e., crystallographic phasing), just to name a few. Biophysical and structural investigations of chemically modified DNAs and RNAs, particularly of nucleic acid analogs with more significant alterations to the well-known base-sugar-phosphate framework (i.e., peptide or hexopyranose nucleic acids), can also provide insights into the properties of the natural nucleic acids that are beyond the reach of studies focusing on DNA or RNA alone. In this review we summarize results from crystallographic analyses of chemically modified DNAs and RNAs that are primarily of interest in the context of the discovery and development of oligonucleotide-based therapeutics. In addition, we re-examine recent structural data on nucleic acid analogs that are investigated as part of a systematic effort to rationalize nature's choice of pentose in the genetic system.

Oligoribonucleotide analogues containing a mixed backbone of phosphodiester and formacetal internucleoside linkages, together with vicinal 2'-O-methyl groups

Oligoribonucleotides containing formacetal internucleoside linkages have been prepared and studied by UV melting experiments. In RNA duplexes, the formacetal substitution is stabilizing (Deltat(m)=0 to +0.9 degrees C per modification) at physiological salt concentrations (0.1 M) but destabilizing (Deltat(m)=-0.4 to -0.8 degrees C per modification) at high salt concentrations (1 M); this suggests that reduction of electrostatic repulsion contributes substantially to the stabilization. The presence of 2'-O-Me substituents increases the stabilities of the duplexes (Deltat(m)=+0.5 to +1.1 degrees C per modification). The positive effects of formacetals and 2'-O-Me groups were independent and additive. (1)H NMR studies on monomeric model compounds containing 3'-(ethyl phosphate) or 3'-O-ethoxymethyl groups showed that the formacetal and 2'-O-Me substitutions shift the conformational equilibria of the ribose residues towards the North conformers by 5 to 12 %. Although the preference for the North conformers qualitatively correlates with increased duplex stabilities, changes in thermodynamic parameters (DeltaH degrees and TDeltaS degrees ) for formation of oligonucleotide duplexes and differences in dependence on concentrations of sodium acetate, sodium chloride and sodium perchlorate suggest that solvation effects are also important for the duplex stabilities. Overall the formacetal linkages fit well in A-type RNA duplexes, making them potentially interesting modifications for RNA-based gene-control strategies (e.g., antisense and RNA interference).

Simulations of RNA base pairs in a nanodroplet reveal solvation-dependent stability

We show that RNA base pairs have variable stability depending on their degree of solvation. This finding has far-reaching biological implications for nucleic acid structure in a partially solvated cellular environment such as inside RNA-protein complexes. Molecular dynamics simulations of partially solvated Watson-Crick RNA base pairs show that whereas water serves to destabilize a base pair by competing for and disrupting base-base hydrogen bonds, when sufficient water molecules are present, fewer hydrogen bonds are available to disrupt the base pairs and the destabilization effect is reduced. The result is that base pairs exist at a stability minimum when solvated in between 20 and 100 water molecules, the upper limit of which corresponds to the approximate number of water molecules contained in the first hydration shell.

Binary hammerhead ribozymes with improved catalytic activity

A new design of binary hammerhead ribozymes displaying high catalytic activity and nucleolytic stability is described. These catalytic structures consist of two partially complementary oligoribonucleotides, capable of assembling into the hammerhead-like structure without tetraloop II on binding to the RNA target. A series of these binary ribozymes targeting the translation initiation region of multiple drug resistance gene mdr1 mRNA was synthesized and assessed in terms of catalytic activity under single and multiple reaction turnover conditions. Enhanced nuclease resistance of the binary ribozymes was achieved by incorporation of 2'-modified nucleotides at selected positions, along with addition of a 3'-3'-linked thymidine cap. The new binary ribozymes exhibit higher RNA cleavage activity than their full-length analogs because of faster dissociation of cleavage products. Furthermore, an excess of one of the ribozyme strands provides the possibility to unfold structured regions of the target RNA and facilitate productive complex formation.

Recognition of RNA by amide modified backbone nucleic acids: molecular dynamics simulations of DNA-RNA hybrids in aqueous solution

Thermodynamic and structural properties of a chemically modified DNA-RNA hybrid in which a phosphodiester linkage is replaced by a neutral amide-3 linkage (3'-CH(2)-CONH-5') were investigated using UV melting experiments, molecular dynamics simulations in explicit water, and continuum solvent models. van't Hoff analysis of the experimental UV melting curves suggests that the significant increase of the thermodynamic stability of a 15-mer DNA-RNA with seven alternated amide-3 modifications (+11 degrees C) is mainly due to an increased binding enthalpy. To further evaluate the origin in the observed affinities differences, the electrostatic contribution to the binding free energy was calculated by solving the Poisson-Boltzmann equation numerically. The nonelectrostatic contribution was estimated as the product of a hydrophobic surface tension coefficient and the surface area that is buried upon double strand formation. Structures were taken from 10 ns molecular dynamics simulations computed in a consistent fashion using explicit solvent, counterions, and the particle-mesh Ewald procedure. The present preliminary thermodynamic study suggests that the favorable binding free energy of the amide-3 DNA single strand to the complementary RNA is equally driven by electrostatic and nonpolar contributions to the binding compared to their natural analogues. In addition, molecular dynamics simulations in explicit water were performed on an amide-3 DNA single strand and the corresponding natural DNA. Results from the conformations cluster analysis of the simulated amide-3 DNA single strand ensembles suggest that the 25% of the population sampled within 10 ns has a pre-organized conformation where the sugar C3' endo pucker is favored at the 3'-flanking nucleotides. These structural and thermodynamic features contribute to the understanding of the observed increased affinities of the amide-3 DNA-RNA hybrids at the microscopic level.

Simulation and modeling of nucleic acid structure, dynamics and interactions

In moving towards the simulation of larger nucleic acid assemblies over longer timescales that include more accurate representations of the environment, we are nearing the end of an era characterized by single nanosecond molecular dynamics simulation of nucleic acids. We are excited by the promise and predictability of the modeling methods, yet remain prudently cautious of sampling and force field limitations. Highlights include the accurate representation of subtle drug-DNA interactions, the detailed study of modified and unusual nucleic acid structures, insight into the influence of dynamics on the structure of DNA, and exploration of the interaction of solvent and ions with nucleic acids.

Pronounced instability of tandem IU base pairs in RNA

Optical melting was used to determine the stabilities of three series of RNA oligomers containing tandem XU base pairs, GGCXUGCC (5'XU3'), GGCUXGCC (5'UX3') and GGCXXGGC/CCGUUCCG (5'XX3'), where X is either A, G or I (inosine). The helices containing tandem AU base pairs were the most stable in the first two series (5'XU3' and 5'UX3'), with an average melting temperature approximately 11 degrees C higher than the helices with tandem 5'GU3' base pairs and 25 degrees C higher than the helices with tandem 5'IU3' base pairs. For the third series (5'XX3'), the helix containing tandem GG is the most stable, with an average melting temperature approximately 2 degrees C higher than the helix with tandem AA base pairs and approximately 24 degrees C higher than the helix with tandem II base pairs. The thermodynamic stability of the oligomers with tandem IU base pairs was also investigated as a function of magnesium ion concentration. As with normal A-U or G-U tandem duplexes, the data could best be interpreted as non-specific binding of magnesium ions to the inosine-containing RNA oligonucleotides.

Hydration of short DNA, RNA and 2'-OMe oligonucleotides determined by osmotic stressing

Studies on hydration are important for better understanding of structure and function of nucleic acids. We compared the hydration of self-complementary DNA, RNA and 2'-O-methyl (2'-OMe) oligonucleotides GCGAAUUCGC, (UA)6 and (CG)3 using the osmotic stressing method. The number of water molecules released upon melting of oligonucleotide duplexes, Delta(n)W, was calculated from the dependence of melting temperature on water activity and the enthalpy, both measured with UV thermal melting experiments. The water activity was changed by addition of ethylene glycol, glycerol and acetamide as small organic co-solutes. The Delta(n)W was 3-4 for RNA duplexes and 2-3 for DNA and 2'-OMe duplexes. Thus, the RNA duplexes were hydrated more than the DNA and the 2'-OMe oligonucleotide duplexes by approximately one to two water molecules depending on the sequence. Consistent with previous studies, GC base pairs were hydrated more than AU pairs in RNA, whereas in DNA and 2'-OMe oligonucleotides the difference in hydration between these two base pairs was relatively small. Our data suggest that the better hydration of RNA contributes to the increased enthalpic stability of RNA duplexes compared with DNA duplexes.

Stabilisation of the transition state of phosphodiester bond cleavage within linear single-stranded oligoribonucleotides

The effect of base sequence on the stability of the transition state (TS) of phosphodiester bond cleavage within linear single-stranded oligoribonucleotides has been studied in order to better understand why the reactivity of some phosphodiester bonds is enhanced compared to an unconstrained linkage. Molecular dynamics simulations of 3.0 ns were carried out for 14 oligonucleotides that contain in the place of the scissile phosphodiester bond a phosphorane structure mimicking the TS of the bond cleavage. The hydrolytic stability of the same oligonucleotides had previously been reported. Both the non-bridging oxyanions and the leaving 5[prime or minute]-oxygen of the pentacoordinated phosphorane moiety were observed to form hydrogen bonds with solvent water molecules in a similar way with all the compounds studied. In addition, water mediated hydrogen bonds between the phosphorane non-bridging oxyanions and the bases of the 3[prime or minute]-flanking sequence were detected with some of the compounds, but not with the most labile ones. Hence, it seems that the enhanced cleavage of some internucleosidic linkages does not result from the TS stabilisation by hydrogen bonding. With heterooligomers, the stacking of bases next to the cleavage site was observed to be enhanced on going from the initial state to the TS, whereas within uracil homooligomer, having initially negligible stacking, no change in the magnitude of stacking was seen. Accordingly, while strong stacking in the initial state is known to retard the phosphodiester bond cleavage, it may in the TS accelerate the reaction. Therefore, enhanced stacking on going from the initial state to transition state appears to be a factor that markedly contributes to the hydrolytic stability of phosphodiester bonds within oligonucleotides and may, at least partly, explain accelerated cleavage compared to fully unconstrained bonds, such as those in polyuridylic acid.

The base sequence dependent flexibility of linear single-stranded oligoribonucleotides correlates with the reactivity of the phosphodiester bond

The effect of base sequence on the structure and flexibility of linear single-stranded RNA molecules and the influence of the base sequence on phosphodiester bond reactivity have been studied. Molecular dynamics simulations of 2.1 ns were carried out for nine chimeric oligonucleotides containing only one unsubstituted ribo unit, all the rest of sugars being 2'-O-methylated. The base sequence has recently been reported to make a big contribution to the reactivity of these compounds. A detailed examination of the interaction energies between the base moieties shows that base stacking is strongly context-dependent and cooperative. The strength of stacking at the site susceptible to chain cleavage by intramolecular transesterification was observed to be dependent on both the flanking bases of the cleavage site and those further apart in the molecule. The interaction energies between the bases in the vicinity of the scissile linkage were found to correlate well with the experimental phosphodiester bond cleavage rates: the stronger the bases close to the cleavage site are stacked, the slower the cleavage rate is.

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.