RNA hydration: a detailed look

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

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

Inosine-Induced Base Pairing Diversity during Reverse Transcription

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

On Water Arrangements in Right- and Left-Handed DNA Structures

DNA requires hydration to maintain its structural integrity. Crystallographic analyses have enabled patterns of water arrangements to be visualized. We survey these water motifs in this review, focusing on left- and right-handed duplex and quadruplex DNAs, together with the i-motif. Common patterns of linear spines of water organization in grooves have been identified and are widely prevalent in right-handed duplexes and quadruplexes. By contrast, a left-handed quadruplex has a distinctive wheel of hydration populating the almost completely circular single groove in this structure.

Structural basis for water modulating RNA duplex formation in the CUG repeats of myotonic dystrophy type 1

Secondary structures formed by expanded CUG RNA are involved in the pathobiology of myotonic dystrophy type 1. Understanding the molecular basis of toxic RNA structures can provide insights into the mechanism of disease pathogenesis and accelerate the drug discovery process. Here, we report the crystal structure of CUG repeat RNA containing three U-U mismatches between C-G and G-C base pairs. The CUG RNA crystallizes as an A-form duplex, with the first and third U-U mismatches adopting a water-mediated asymmetric mirror isoform geometry. We found for the first time that a symmetric, water-bridged U-H₂O-U mismatch is well tolerated within the CUG RNA duplex, which was previously suspected but not observed. The new water-bridged U-U mismatch resulted in high base-pair opening and single-sided cross-strand stacking interactions, which in turn dominate the CUG RNA structure. Furthermore, we performed molecular dynamics (MD) simulations that complemented the structural findings and proposed that the first and third U-U mismatches are interchangeable conformations, while the central water-bridged U-U mismatch represents an intermediate state that modulates the RNA duplex conformation. Collectively, the new structural features provided in this work are important for understanding the recognition of U-U mismatches in CUG repeats by external ligands such as proteins or small molecules.

Modified Nucleotides for Chemical and Enzymatic Synthesis of Therapeutic RNA

In recent years, RNA has emerged as a medium with a broad spectrum of therapeutic potential, however, for years, a group of short RNA fragments was studied and considered therapeutic molecules. In nature, RNA plays both functions, with coding and non-coding potential. For RNA, like any other therapeutic, to be used clinically, certain barriers must be crossed. Among them, there are biocompatibility, relatively low toxicity, bioavailability, increased stability, target efficiency and low off-target effects. In the case of RNA, most of these obstacles can be overcome by incorporating modified nucleotides into its structure. This may be achieved by both, in vitro and in vivo biosynthetic methods, as well as chemical synthesis. Some advantages and disadvantages of each approach are summarized here. The wide range of nucleotide analogues has been tested for their utility as monomers for RNA synthesis. Many of them have been successfully implemented, and a lot of pre-clinical and clinical studies involving modified RNA have been carried out. Some of these medications have already been introduced into clinics. After the huge success of RNA-based vaccines that were introduced into widespread use in 2020, and the introduction to the market of some RNA-based drugs, RNA therapeutics containing modified nucleotides appear to be the future of medicine.

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.

Measuring RNA UNCG Tetraloop Refolding Dynamics Using Temperature-Jump/Drop Infrared Spectroscopy

Determining the structural dynamics of RNA and DNA is essential to understanding their cellular function, but direct measurement of strand association or folding remains experimentally challenging. Here we illustrate a temperature-jump/drop method able to reveal refolding dynamics. Time-resolved temperature-jump/drop infrared spectroscopy is used to measure the melting and refolding dynamics of a 12-nucleotide RNA sequence comprising a UACG tetraloop and a four-base-pair double-stranded GC stem, comparing them to an equivalent DNA (TACG) sequence. Stem-loop melting occurred an order of magnitude more slowly in RNA than DNA (6.0 ± 0.1 μs versus 0.8 ± 0.1 μs at 70 °C). In contrast, the refolding dynamics of both sequences occurred on similar time scales (200 μs). While the melting and refolding dynamics of RNA and DNA hairpins both followed Arrhenius temperature dependences, refolding was characterized by an apparent negative activation energy, consistent with a mechanism involving multiple misfolded intermediates prior to zipping of the stem base pairs.

Cryo neutron crystallography demonstrates influence of RNA 2'-OH orientation on conformation, sugar pucker and water structure

The ribose 2′-hydroxyl is the key chemical difference between RNA and DNA and primary source of their divergent structural and functional characteristics. Macromolecular X-ray diffraction experiments typically do not reveal the positions of hydrogen atoms. Thus, standard crystallography cannot determine 2′-OH orientation (H2′-C2′-O2′-HO2′ torsion angle) and its potential roles in sculpting the RNA backbone and the expansive fold space. Here, we report the first neutron crystal structure of an RNA, the Escherichia coli rRNA Sarcin-Ricin Loop (SRL). 2′-OD orientations were established for all 27 residues and revealed O-D bonds pointing toward backbone (O3′, 13 observations), nucleobase (11) or sugar (3). Most riboses in the SRL stem region show a 2′-OD backbone-orientation. GAGA-tetraloop riboses display a 2′-OD base-orientation. An atypical C2′-endo sugar pucker is strictly correlated with a 2′-OD sugar-orientation. Neutrons reveal the strong preference of the 2′-OH to donate in H-bonds and that 2′-OH orientation affects both backbone geometry and ribose pucker. We discuss 2′-OH and water molecule orientations in the SRL neutron structure and compare with results from a solution phase 10 μs MD simulation. We demonstrate that joint cryo-neutron/X-ray crystallography offers an all-in-one approach to determine the complete structural properties of RNA, i.e. geometry, conformation, protonation state and hydration structure.

Data-informed reparameterization of modified RNA and the effect of explicit water models: application to pseudouridine and derivatives

Pseudouridine is one of the most abundant post-transcriptional modifications in RNA. We have previously shown that the FF99-derived parameters for pseudouridine and some of its naturally occurring derivatives in the AMBER distribution either alone or in combination with the revised γ torsion parameters (parmbsc0) failed to reproduce their conformational characteristics observed experimentally (Deb et al. in J Chem Inf Model 54:1129–1142, 2014; Deb et al. in J Comput Chem 37:1576–1588, 2016; Dutta et al. in J Chem Inf Model 60:4995–5002, 2020). However, the application of the recommended bsc0 correction did lead to an improvement in the description not only of the distribution in the γ torsional space but also of the sugar pucker distributions. In an earlier study, we examined the transferability of the revised glycosidic torsion parameters (χ_IDRP) for Ψ to its derivatives. We noticed that although these parameters in combination with the AMBER FF99-derived parameters and the revised γ torsional parameters resulted in conformational properties of these residues that were in better agreement with experimental observations, the sugar pucker distributions were still not reproduced accurately. Here we report a new set of partial atomic charges for pseudouridine, 1-methylpseudouridine, 3-methylpseudouridine and 2′-O-methylpseudouridine and a new set of glycosidic torsional parameters (χ_ND) based on chosen glycosidic torsional profiles that most closely corresponded to the NMR data for conformational propensities and studied their effect on the conformational distributions using REMD simulations at the individual nucleoside level. We have also studied the effect of the choice of water model on the conformational characteristics of these modified nucleosides. Our observations suggest that the current revised set of parameters and partial atomic charges describe the sugar pucker distributions for these residues more accurately and that the choice of a suitable water model is important for the accurate description of their conformational properties. We have further validated the revised sets of parameters by studying the effect of substitution of uridine with pseudouridine within single stranded RNA oligonucleotides on their conformational and hydration characteristics.

Backbone Hydrocarbon-Constrained Nucleic Acids Modulate Hybridization Kinetics for RNA

The binding affinity of therapeutic oligonucleotides (ONs) for their cognate RNA is determined by the rates of association (k_a) and dissociation (k_d). Single-stranded ONs are highly flexible and can adopt multiple conformations in solution, some of which may not be conducive for hybridization. We investigated if restricting rotation around the sugar-phosphate backbone, by tethering two adjacent backbone phosphonate esters using hydrocarbon bridges, can modulate hybridization kinetics of the modified ONs for complementary RNA. Given the large number of possible analogues with different tether lengths and configurations at the phosphorus atoms, we employed molecular dynamic simulations to optimize the size of the hydrocarbon bridge to guide the synthetic efforts. The backbone-constrained nucleotide trimers with stereodefined configurations at the contiguous backbone phosphorus atoms were assembled using a ring-closing metathesis reaction, then incorporated into oligonucleotides by an in situ synthesis of the phosphoramidites followed by coupling to solid supports. Evaluation of the modified oligonucleotides revealed that 15-membered macrocyclic-constrained analogues displayed similar or slightly improved on-rates but significantly increased off-rates compared to unmodified DNA ONs, resulting in reduced duplex stability. In contrast, LNA ONs with conformationally preorganized furanose rings showed similar on-rates to DNA ONs but very slow off-rates, resulting in net improvement in duplex stability. Furthermore, the experimental data generally supported the molecular dynamics simulation results, suggesting that this strategy can be used as a predictive tool for designing the next generation of constrained backbone ON analogues with improved hybridization properties.

Contact pairs of RNA with magnesium ions-electrostatics beyond the Poisson-Boltzmann equation

The electrostatic interaction of RNA with its aqueous environment is most relevant for defining macromolecular structure and biological function. The attractive interaction of phosphate groups in the RNA backbone with ions in the water environment leads to the accumulation of positively charged ions in the first few hydration layers around RNA. Electrostatics of this ion atmosphere and the resulting ion concentration profiles have been described by solutions of the nonlinear Poisson-Boltzmann equation and atomistic molecular dynamics (MD) simulations. Much less is known on contact pairs of RNA phosphate groups with ions at the RNA surface, regarding their abundance, molecular geometry, and role in defining RNA structure. Here, we present a combined theoretical and experimental study of interactions of a short RNA duplex with magnesium (Mg2+) ions. MD simulations covering a microsecond time range give detailed hydration geometries as well as electrostatics and spatial arrangements of phosphate-Mg2+ pairs, including both pairs in direct contact and separated by a single water layer. The theoretical predictions are benchmarked by linear infrared absorption and nonlinear two-dimensional infrared spectra of the asymmetric phosphate stretch vibration which probes both local interaction geometries and electric fields. Contact pairs of phosphate groups and Mg2+ ions are identified via their impact on the vibrational frequency position and line shape. A quantitative analysis of infrared spectra for a range of Mg2+-excess concentrations and comparison with fluorescence titration measurements shows that on average 20-30% of the Mg2+ ions interacting with the RNA duplex form contact pairs. The experimental and MD results are in good agreement. In contrast, calculations based on the nonlinear Poisson-Boltzmann equation fail in describing the ion arrangement, molecular electrostatic potential, and local electric field strengths correctly. Our results underline the importance of local electric field mapping and molecular-level simulations to correctly account for the electrostatics at the RNA-water interface.

The mutual interactions of RNA, counterions and water - quantifying the electrostatics at the phosphate-water interface

The structure and dynamics of polyanionic biomolecules, like RNA, are decisively determined by their electric interactions with the water molecules and the counterions in the environment. The solvation dynamics of the biomolecules involves a subtle balance of non-covalent and many-body interactions with structural fluctuations due to thermal motion occurring in a femto- to subnanosecond time range. This complex fluctuating many particle scenario is crucial in defining the properties of biological interfaces with far reaching significance for the folding of RNA structures and for facilitating RNA-protein interactions. Given the inherent complexity, suited model systems, carefully calibrated and benchmarked by experiments, are required to quantify the relevant interactions of RNA with the aqueous environment. In this feature article we summarize our recent progress in the understanding of the electrostatics at the biological interface of double stranded RNA (dsRNA) and transfer RNA (tRNA). Dimethyl phosphate (DMP) is introduced as a viable and rigorously accessible model system allowing the interaction strength with water molecules and counterions, their relevant fluctuation timescales and the spatial reach of interactions to be established. We find strong (up to ≈90 MV cm-1) interfacial electric fields with fluctuations extending up to ≈20 THz and demonstrate how the asymmetric stretching vibration νAS(PO2)- of the polarizable phosphate group can serve as the most sensitive probe for interfacial interactions, establishing a rigorous link between simulations and experiment. The approach allows for the direct interfacial observation of interactions of biologically relevant Mg2+ counterions with phosphate groups in contact pair geometries via the rise of a new absorption band imposed by exchange repulsion interactions at short interatomic distances. The systematic extension to RNA provides microscopic insights into the changes of the hydration structure that accompany the temperature induced melting of the dsRNA double helix and quantify the ionic interactions in the folded tRNA. The results show that pairs of negatively charged phosphate groups and Mg2+ ions represent a key structural feature of RNA embedded in water. They highlight the importance of binding motifs made of contact pairs in the electrostatic stabilization of RNA structures that have a strong impact on the surface potential and enable the fine tuning of the local electrostatic properties which are expected to be relevant for mediating the interactions between biomolecules.

Evaluation of Phosphorus and Non-Phosphorus Neutral Oligonucleotide Backbones for Enhancing Therapeutic Index of Gapmer Antisense Oligonucleotides

The phosphorothioate (PS) linkage in an essential component of therapeutic oligonucleotides. PS in the DNA region of gapmer antisense oligonucleotides (ASOs) supports RNaseH1 activity and enhances nuclease stability. PS also promotes binding to plasma, cell surface, and intracellular proteins, which facilitates tissue distribution, cellular uptake, and endosomal escape of PS ASOs. We recently showed that site-specific replacement of PS in the DNA gap with methoxylpropyl phosphonate (MOP) linkages can enhance the therapeutic index of gapmer ASOs. In this article, we explored 18 phosphorus- and non-phosphorus-based neutral backbone modifications to determine the structure-activity relationship of neutral linkages for enhancing therapeutic index. Replacing MOP with other alkyl phosphonate and phosphotriester linkages enhanced therapeutic index, but these linkages were susceptible to chemical degradation during oligonucleotide deprotection from solid supports following synthesis. Replacing MOP with non-phosphorus linkages resulted in improved chemical stability, but these linkages were introduced into ASOs as nucleotide dimers, which limits their versatility. Overall, linkages such as isopropyl and isobutyl phosphonates and O-isopropyl and O-tetrahydrofuranosyl phosphotriesters, formacetal, and C3-amide showed improved activity in mice relative to MOP. Our data suggest that site-specific incorporation of any neutral backbone linkage can improve therapeutic index, but the size, hydrophobicity, and RNA-binding affinity of the linkage influence ASO activity.

Water structure around a left-handed Z-DNA fragment analyzed by cryo neutron crystallography

Even in high-quality X-ray crystal structures of oligonucleotides determined at a resolution of 1 Å or higher, the orientations of first-shell water molecules remain unclear. We used cryo neutron crystallography to gain insight into the H-bonding patterns of water molecules around the left-handed Z-DNA duplex [d(CGCGCG)]₂. The neutron density visualized at 1.5 Å resolution for the first time allows us to pinpoint the orientations of most of the water molecules directly contacting the DNA and of many second-shell waters. In particular, H-bond acceptor and donor patterns for water participating in prominent hydration motifs inside the minor groove, on the convex surface or bridging nucleobase and phosphate oxygen atoms are finally revealed. Several water molecules display entirely unexpected orientations. For example, a water molecule located at H-bonding distance from O6 keto oxygen atoms of two adjacent guanines directs both its deuterium atoms away from the keto groups. Exocyclic amino groups of guanine (N2) and cytosine (N4) unexpectedly stabilize waters H-bonded to O2 keto oxygens from adjacent cytosines and O6 keto oxygens from adjacent guanines, respectively. Our structure offers the most detailed view to date of DNA solvation in the solid-state undistorted by metal ions or polyamines.

Phosphate Vibrations Probe Electric Fields in Hydrated Biomolecules: Spectroscopy, Dynamics, and Interactions

Electric interactions have a strong impact on the structure and dynamics of biomolecules in their native water environment. Given the variety of water arrangements in hydration shells and the femto- to subnanosecond time range of structural fluctuations, there is a strong quest for sensitive noninvasive probes of local electric fields. The stretching vibrations of phosphate groups, in particular the asymmetric (PO₂)⁻ stretching vibration ν_AS(PO₂)⁻, allow for a quantitative mapping of dynamic electric fields in aqueous environments via a field-induced redshift of their transition frequencies and concomitant changes of vibrational line shapes. We present a systematic study of ν_AS(PO₂)⁻ excitations in molecular systems of increasing complexity, including dimethyl phosphate (DMP), short DNA and RNA duplex structures, and transfer RNA (tRNA) in water. A combination of linear infrared absorption, two-dimensional infrared (2D-IR) spectroscopy, and molecular dynamics (MD) simulations gives quantitative insight in electric-field tuning rates of vibrational frequencies, electric field and fluctuation amplitudes, and molecular interaction geometries. Beyond neat water environments, the formation of contact ion pairs of phosphate groups with Mg²⁺ ions is demonstrated via frequency upshifts of the ν_AS(PO₂)⁻ vibration, resulting in a distinct vibrational band. The frequency positions of contact geometries are determined by an interplay of attractive electric and repulsive exchange interactions.

Magnesium Contact Ions Stabilize the Tertiary Structure of Transfer RNA: Electrostatics Mapped by Two-Dimensional Infrared Spectra and Theoretical Simulations

Ions interacting with hydrated RNA play a central role in defining its secondary and tertiary structure. While spatial arrangements of ions, water molecules, and phosphate groups have been inferred from X-ray studies, the role of electrostatic and other noncovalent interactions in stabilizing compact folded RNA structures is not fully understood at the molecular level. Here, we demonstrate that contact ion pairs of magnesium (Mg²⁺) and phosphate groups embedded in local water shells stabilize the tertiary equilibrium structure of transfer RNA (tRNA). Employing dialyzed tRNA^Phe from yeast and tRNA from Escherichia coli, we follow the population of Mg²⁺ sites close to phosphate groups of the ribose-phosphodiester backbone step by step, combining linear and nonlinear infrared spectroscopy of phosphate vibrations with molecular dynamics simulations and ab initio vibrational frequency calculations. The formation of up to six Mg²⁺/phosphate contact pairs per tRNA and local field-induced reorientations of water molecules balance the phosphate–phosphate repulsion in nonhelical parts of tRNA, thus stabilizing the folded structure electrostatically. Such geometries display limited sub-picosecond fluctuations in the arrangement of water molecules and ion residence times longer than 1 μs. At higher Mg²⁺ excess, the number of contact ion pairs per tRNA saturates around 6 and weakly interacting ions prevail. Our results suggest a predominance of contact ion pairs over long-range coupling of the ion atmosphere and the biomolecule in defining and stabilizing the tertiary structure of tRNA.

Microsecond Equilibrium Dynamics of Hairpin-Forming Oligonucleotides Quantified by Two-Color Two-Dimensional Fluorescence Lifetime Correlation Spectroscopy

RNA and DNA play distinct roles in biological systems. However, the underlying physicochemical difference has been poorly understood, in particular, that in dynamical aspects. In this paper, we report on a comparative study of the formation-dissociation dynamics of a hairpin structure of RNA and DNA with development of two-color two-dimensional fluorescence lifetime correlation spectroscopy (two-color 2D FLCS). In this extension of 2D FLCS, we newly introduce the two-color detection scheme to analyze not only donor fluorescence photons but also acceptor fluorescence photons from a doubly labeled Förster resonance energy transfer (FRET) pair. This new 2D FLCS is utilized to resolve multiple species present in an equilibrated condition with a microsecond time resolution and enhanced sensitivity, and the combined use with the filtered fluorescence correlation spectroscopy (FCS) method enables a quantitative discussion on microsecond structural dynamics occurring in the equilibrium. This integrated approach is applied to FRET-labeled RNA/DNA oligonucleotides having analogous hairpin-forming sequences, and it was revealed that the hairpin dissociation rate of RNA is an order of magnitude slower than that of DNA while their hairpin-forming rates are comparable. This marked difference is attributable to the distinct duplex structure of RNA and DNA. The present study demonstrates that the integrated approach combining two-color 2D FLCS and filtered FCS has a high potential for quantifying microsecond kinetics at the single-molecule level, which allows us to experimentally construct a free energy landscape.

Synthesis, chirality-dependent conformational and biological properties of siRNAs containing 5'-(R)- and 5'-(S)-C-methyl-guanosine

Various chemical modifications have been identified that enhance potency of small interfering RNAs (siRNAs) and that reduce off-target effects, immune stimulation, and toxicities of metabolites of these therapeutic agents. We previously described 5′-C-methyl pyrimidine nucleotides also modified at the 2′ position of the sugar. Here, we describe the synthesis of 2′-position unmodified 5′-(R)- and 5′-(S)-C-methyl guanosine and evaluation of these nucleotides in the context of siRNA. The (R) isomer provided protection from 5′ exonuclease and the (S) isomer provided protection from 3′ exonuclease in the context of a terminally modified oligonucleotide. siRNA potency was maintained when these modifications were incorporated at the tested positions of sense and antisense strands. Moreover, the corresponding 5′ triphosphates were not substrates for mitochondrial DNA polymerase. Models generated based on crystal structures of 5′ and 3′ exonuclease oligonucleotide complexes with 5′-(R)- and 5′-(S)-C-methyl substituents attached to the 5′- and 3′-terminal nucleotides, respectively, provided insight into the origins of the observed protections. Structural properties of 5′-(R)-C-methyl guanosine incorporated into an RNA octamer were analysed by X-ray crystallography, and the structure explains the loss in duplex thermal stability for the (R) isomer compared with the (S) isomer. Finally, the effect of 5′-C-methylation on endoribonuclease activity has been explained.

Intrinsic Contributions of 2'-Hydroxyl to the Hydration of Nucleosides at the Monomeric Level

Although many reports have revealed structural differences between DNA and RNA at the polymeric level, there are no comparative studies with 2'-deoxyribonucleoside and ribonucleoside to explore the role of the 2'-OH group at the monomeric level under the same conditions. Inspired by this, herein, the intrinsic contributions of the 2'-OH group in the nucleoside have been systematically investigated by directly solving the single-crystal structures of 2'-deoxy-2-aminoadenosine (1), 2-aminoadenosine (2), and 2-aminoarabinofuranosyladenine (3) in water. The 2'-OH group not only influenced the conformation and base-pair pattern of the single-nucleoside molecule, but also played a fundamental role in the entire supramolecular structure. Interestingly, compound 1, which did not contain the 2'-OH group, displayed strong hydration, whereas 2 and 3 (with the 2'-OH group in the opposite direction) exhibited no hydration, which was completely different from that observed in nucleic acids. Meanwhile, compound 1 trapped water molecules to form unique trihydrol moieties, which further served as the backbone to construct the simplest double-chain DNA-like structures. To this end, to investigate the effect of the biological environment on these unique structures, the solvent was changed from water to phosphate-buffered saline (PBS). Surprisingly, such a subtle adjustment led to entirely different superstructures, consisting of 2D lamellar structures in water and 3D porous structures in PBS. These large morphological differences could be attributed to delicate ion hydration, which was also confirmed through variable-temperature X-ray analysis, SEM, and intermolecular interaction energy calculations. In summary, this study comprehensively investigated the intrinsic contributions of 2'-hydroxyl to the hydration of nucleosides at the monomeric level; this is helpful to further understand the differences in DNA/RNA and the impact of their surrounding environment.

Effect of Molecular Crowding on the Stability of RNA G-Quadruplexes with Various Numbers of Quartets and Lengths of Loops

G-Quadruplexes are noncanonical structures formed by guanine-rich regions of not only DNA but also RNA. RNA G-quadruplexes are widely present in the transcriptome as mRNAs and noncoding RNAs and take part in various essential functions in cells. Furthermore, stable RNA G-quadruplexes control the extent of biological functions, such as mRNA translation and antigen presentation. To understand and regulate the functions controlled by RNA G-quadruplexes in cellular environments, which are molecularly crowded, we would be required to investigate the stability of G-quadruplexes in molecular crowding. Here, we systematically investigated the thermodynamic stability of RNA G-quadruplexes with different numbers of G-quartets and lengths of loops. The molecular crowding conditions of polyethylene glycol with an average molecular weight of 200 (PEG200) were found to stabilize RNA G-quadruplexes with three and four G-quartets, while G-quadruplexes with two G-quartets did not exhibit any stabilization upon addition of PEG200. On the other hand, no difference in stabilization by PEG200 was observed among the G-quadruplexes with different loop lengths. Thermodynamic analysis of the RNA G-quadruplexes revealed more appropriate motifs for identifying G-quadruplex-forming sequences. The informatics analysis with new motifs demonstrated that the distributions of G-quadruplexes in human noncoding RNAs differed depending on the number of G-quartets. Therefore, RNA G-quadruplexes with different numbers of G-quartets may play different roles in response to environmental changes in cells.

Origins of the Increased Affinity of Phosphorothioate-Modified Therapeutic Nucleic Acids for Proteins

The phosphorothioate backbone modification (PS) is one of the most widely used chemical modifications for enhancing the drug-like properties of nucleic acid-based drugs, including antisense oligonucleotides (ASOs). PS-modified nucleic acid therapeutics show improved metabolic stability from nuclease-mediated degradation and exhibit enhanced interactions with plasma, cell-surface, and intracellular proteins, which facilitates their tissue distribution and cellular uptake in animals. However, little is known about the structural basis of the interactions of PS nucleic acids with proteins. Here, we report a crystal structure of the DNA-binding domain of a model ASO-binding protein PC4, in complex with a full PS 2'-OMe DNA gapmer ASO. To our knowledge this is the first structure of a complex between a protein and fully PS nucleic acid. Each PC4 dimer comprises two DNA-binding interfaces. In the structure one interface binds the 5'-terminal 2'-OMe PS flank of the ASO, while the other interface binds the regular PS DNA central part in the opposite polarity. As a result, the ASO forms a hairpin-like structure. ASO binding also induces the formation of a dimer of dimers of PC4, which is stabilized by base pairing between homologous regions of the ASOs bound by each dimer of PC4. The protein interacts with the PS nucleic acid through a network of electrostatic and hydrophobic interactions, which provides insights into the origins for the enhanced affinity of PS for proteins. The importance of these contacts was further confirmed in a NanoBRET binding assay using a Nano luciferase tagged PC4 acting as the BRET donor, to a fluorescently conjugated ASO acting as the BRET acceptor. Overall, our results provide insights into the molecular forces that govern the interactions of PS ASOs with cellular proteins and provide a potential model for how these interactions can template protein-protein interactions causative of cellular toxicity.

Change of Hydration Patterns upon RNA Melting Probed by Excitations of Phosphate Backbone Vibrations

The water hydration shell has a decisive impact on the structural and functional properties of RNA. Changes of the RNA structure upon melting and in biochemical processes are accompanied by a change of hydration patterns, a process that is barely characterized. To discern hydration geometries around the backbone phosphate groups of an RNA double helix at the molecular level, we combine two-dimensional infrared spectroscopy of phosphate vibrations with theoretical simulations. There are three distinct coexisting hydration motifs of the RNA A-helix: an ordered chain-like arrangement of water molecules with links between neighboring phosphate groups, separate local hydration shells of up to six water molecules, and hydrated phosphate/counterion contact pairs. RNA disordering upon melting is connected with a transition from predominantly ordered water structures to local hydration shells around phosphate units. Structural fluctuations are dominated by librational water motions occurring on a 300 fs time scale, without exchange between hydration motifs.

Contact Ion Pairs of Phosphate Groups in Water: Two-Dimensional Infrared Spectroscopy of Dimethyl Phosphate and ab Initio Simulations

The interaction of phosphate groups with ions in an aqueous environment has a strong impact on the structure and folding processes of DNA and RNA. The dynamic variety of ionic arrangements, including both contact pairs and water separated ions, and the molecular coupling mechanisms are far from being understood. In a combined experimental and theoretical approach, we address the properties of contact ion pairs of the prototypical system dimethyl phosphate with Na⁺, Ca²⁺, and Mg²⁺ ions in water. Linear and femtosecond two-dimensional infrared (2D-IR) spectroscopy of the asymmetric (PO₂)^- stretching vibration separates and characterizes the different species via their blue-shifted vibrational signatures and 2D-IR line shapes. Phosphate-magnesium contact pairs stand out as the most compact geometry while the contact pairs with Ca²⁺ and Na⁺ display a wider structural variation. Microscopic density functional theory simulations rationalize the observed frequency shifts and reveal distinct differences between the contact geometries.

Why are Hoogsteen base pairs energetically disfavored in A-RNA compared to B-DNA?

A(syn)-U/T and G(syn)-C⁺ Hoogsteen (HG) base pairs (bps) are energetically more disfavored relative to Watson–Crick (WC) bps in A-RNA as compared to B-DNA by >1 kcal/mol for reasons that are not fully understood. Here, we used NMR spectroscopy, optical melting experiments, molecular dynamics simulations and modified nucleotides to identify factors that contribute to this destabilization of HG bps in A-RNA. Removing the 2′-hydroxyl at single purine nucleotides in A-RNA duplexes did not stabilize HG bps relative to WC. In contrast, loosening the A-form geometry using a bulge in A-RNA reduced the energy cost of forming HG bps at the flanking sites to B-DNA levels. A structural and thermodynamic analysis of purine-purine HG mismatches reveals that compared to B-DNA, the A-form geometry disfavors syn purines by 1.5–4 kcal/mol due to sugar-backbone rearrangements needed to sterically accommodate the syn base. Based on MD simulations, an additional penalty of 3–4 kcal/mol applies for purine-pyrimidine HG bps due to the higher energetic cost associated with moving the bases to form hydrogen bonds in A-RNA versus B-DNA. These results provide insights into a fundamental difference between A-RNA and B-DNA duplexes with important implications for how they respond to damage and post-transcriptional modifications.

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.

Interactions of RNA and Water probed by 2D-IR Spectroscopy

Combined experimental-theoretical investigation of ultrafast hydration dynamics of an A-form RNA double helix in water reveals an ordered arrangement of water molecules and provides boundary conditions for the ion atmosphere around the polyanionic RNA.

Polarizable force field for RNA based on the classical drude oscillator

RNA molecules are highly dynamic and capable of adopting a wide range of complex, folded structures. The factors driving the folding and dynamics of these structures are dependent on a balance of base pairing, hydration, base stacking, ion interactions, and the conformational sampling of the 2'-hydroxyl group in the ribose sugar. The representation of these features is a challenge for empirical force fields used in molecular dynamics simulations. Toward meeting this challenge, the inclusion of explicit electronic polarization is important in accurately modeling RNA structure. In this work, we present a polarizable force field for RNA based on the classical Drude oscillator model, which represents electronic degrees of freedom via negatively charged particles attached to their parent atoms by harmonic springs. Beginning with parametrization against quantum mechanical base stacking interaction energy and conformational energy data, we have extended the Drude-2017 nucleic acid force field to include RNA. The conformational sampling of a range of RNA sequences were used to validate the force field, including canonical A-form RNA duplexes, stem-loops, and complex tertiary folds that bind multiple Mg2+ ions. Overall, the Drude-2017 RNA force field reproduces important properties of these structures, including the conformational sampling of the 2'-hydroxyl and key interactions with Mg2+ ions. © 2018 Wiley Periodicals, Inc.

Cryo-neutron crystallographic data collection and preliminary refinement of left-handed Z-DNA d(CGCGCG)

Crystals of left-handed Z-DNA [d(CGCGCG)]2 diffract X-rays to beyond 1 Å resolution, feature a small unit cell (∼18 × 31 × 44 Å) and are well hydrated, with around 90 water molecules surrounding the duplex in the asymmetric unit. The duplex shows regular hydration patterns in the narrow minor groove, on the convex surface and around sugar-phosphate backbones. Therefore, Z-DNA offers an ideal case to test the benefits of low-temperature neutron diffraction data collection to potentially determine the donor-acceptor patterns of first- and second-shell water molecules. Nucleic acid fragments pose challenges for neutron crystallography because water molecules are located on the surface rather than inside sequestered spaces such as protein active sites or channels. Water molecules can be expected to display dynamic behavior, particularly in cases where water is not part of an inner shell and directly coordinated to DNA atoms. Thus, nuclear density maps based on room-temperature diffraction data with a resolution of 1.6 Å did not allow an unequivocal determination of the orientations of water molecules. Here, cryo-neutron diffraction data collection for a Z-DNA crystal on the Macromolecular Neutron Diffractometer at the Spallation Neutron Source at Oak Ridge National Laboratory and the outcome of an initial refinement of the structure are reported. A total of 12 diffraction images were recorded with an exposure time of 3.5 h per image, whereby the crystal was static for each diffraction image with a 10° ϕ rotation between images. Initial refinements using these neutron data indicated the positions and orientations of 30 water molecules within the first hydration shell of the DNA molecule. This experiment constitutes a state-of-the-art approach and is the first attempt to our knowledge to determine the low-temperature neutron structure of a DNA crystal.

RNA as a Complex Polymer with Coupled Dynamics of Ions and Water in the Outer Solvation Sphere

We unravel the internal and collective modes of a widely studied 58-nucleotide rRNA fragment in solvent using atomically detailed molecular dynamics simulations. The variation of lifetimes for water hydrogen bonds with nucleotide groups indicates heterogeneity of water dynamics on the RNA surface. The time scales of interactions of the discrete water molecules with RNA nucleotides extend from several hundred picoseconds to a few nanoseconds. We determine all of the association sites and the spatial distribution of residence times for Mg²⁺, K⁺, and water molecules in those sites. We provide insights into the population of Mg²⁺ and K⁺ ions and water molecules in the outer sphere and how their fluctuations are intricately linked with the kinetics of the 58-mer. We find that many of the long-lived Mg²⁺ sites identified from the simulations agree with the locations of ions in the X-ray structure. We determine the excess ion atmosphere around the rRNA and compare it with experimental data. We investigate the collective behavior of RNA, ions, and water, by performing a joint principle component analysis for the Cartesian coordinates of the RNA phosphorus atoms and for the occupation counts of the association sites. Our results indicate that the 58-mer system is a complex polymer, composed of RNA that is encased by a fluctuating network of associated counterions, co-ions, and water.

Rapid access to RNA resonances by proton-detected solid-state NMR at >100 kHz MAS

Fast (>100 kHz) magic angle spinning solid-state NMR allows combining high-sensitive proton detection with the absence of an intrinsic molecular weight limit. Here we apply this technique to RNA and assign nucleotide spin systems through highly sensitive multidimensional experiments.

Fast (>100 kHz) magic angle spinning solid-state NMR allows combining high-sensitive proton detection with the absence of an intrinsic molecular weight limit. Using this technique we observe for the first time narrow ¹H RNA resonances and assign nucleotide spin systems with only 200 μg of uniformly ¹³C,¹⁵N-labelled RNA.

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

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

Physics-based all-atom modeling of RNA energetics and structure

The database of RNA sequences is exploding, but knowledge of energetics, structures, and dynamics lags behind. All-atom computational methods, such as molecular dynamics, hold promise for closing this gap. New algorithms and faster computers have accelerated progress in improving the reliability and accuracy of predictions. Currently, the methods can facilitate refinement of experimentally determined nuclear magnetic resonance and x-ray structures, but are 'unreliable' for predictions based only on sequence. Much remains to be discovered, however, about the many molecular interactions driving RNA folding and the best way to approximate them quantitatively. The large number of parameters required means that a wide variety of experimental results will be required to benchmark force fields and different approaches. As computational methods become more reliable and accessible, they will be used by an increasing number of biologists, much as x-ray crystallography has expanded. Thus, many fundamental physical principles underlying the computational methods are described. This review presents a summary of the current state of molecular dynamics as applied to RNA. It is designed to be helpful to students, postdoctoral fellows, and faculty who are considering or starting computational studies of RNA. WIREs RNA 2017, 8:e1422. doi: 10.1002/wrna.1422.

Vibrational Dynamics and Couplings of the Hydrated RNA Backbone: A Two-Dimensional Infrared Study

The equilibrium structure of the RNA sugar-phosphate backbone and its hydration shell is distinctly different from hydrated DNA. Applying femtosecond two-dimensional infrared (2D-IR) spectroscopy in a range from 950 to 1300 cm^-1, we elucidate the character, dynamics, and couplings of backbone modes of a double-stranded RNA A-helix geometry in its aqueous environment. The 2D-IR spectra display a greater number of backbone modes than for DNA, with distinctly different lineshapes of diagonal peaks. Phosphate-ribose interactions and local hydration structures are reflected in the complex coupling pattern of RNA modes. Interactions with the fluctuating water shell give rise to spectral diffusion on a 300 fs time scale, leading to a quasi-homogeneous line shape of the symmetric (PO₂)^- stretching mode of the strongly hydrated phosphate groups. The RNA results are benchmarked by 2D-IR spectra of DNA oligomers in water and analyzed by molecular dynamics and quantum mechanical molecular mechanics simulations.

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.

Impact of Oligonucleotide Structure, Chemistry, and Delivery Method on In Vitro Cytotoxicity

Single-stranded (ss) 2'-fluoro (2'-F)-modified oligonucleotides (ONs) with a full phosphorothioate (PS) backbone have been reported to be cytotoxic and cause DNA double-strand breaks (DSBs) when transfected into HeLa cells. However, the molecular determinants of these effects have not been fully explored. In this study, we investigated the impact of ON structure, chemistry, delivery method, and cell type on in vitro cytotoxicity and DSBs. We found that ss PS-ONs were more cytotoxic than double-stranded (ds) PS-ONs, irrespective of the 2'-ribose chemistry, inclusive of the 2'-F modification. Cytotoxicity of ss ONs was most affected by the total PS content, with an additional contribution of 2'-F substitutions in HeLa, but not HepG2, cells. The relatively mild cytotoxicity of ds ONs was most impacted by long contiguous PS stretches combined with 2'-F substitutions. None of the tested ds 2'-F-modified PS-ONs caused DSBs, while the previously reported DSBs caused by ss 2'-F-modified PS-ONs were PS dependent. HeLa cells were more sensitive to ON-mediated toxicity when transfected with Lipofectamine 2000 versus Lipofectamine RNAiMax. Importantly, asialoglycoprotein receptor-mediated uptake of N-acetylgalactosamine-conjugated ss or ds PS-ONs, even those with long PS stretches and high 2'-F content, was neither cytotoxic nor caused DSBs at transfection-equivalent exposures. These results suggest that in vitro cytotoxicity and DSBs associated with ONs are delivery method dependent and primarily determined by single-stranded nature and PS content of ONs.

How Mg2+ ion and water network affect the stability and structure of non-Watson-Crick base pairs in E. coli loop E of 5S rRNA: a molecular dynamics and reference interaction site model (RISM) study

The non-Watson-Crick (non-WC) base pairs of Escherichia coli loop E of 5S rRNA are stabilized by Mg²⁺ ions through water-mediated interaction. It is important to know the synergic role of Mg²⁺ and the water network surrounding Mg²⁺ in stabilizing the non-WC base pairs of RNA. For this purpose, free energy change of the system is calculated using molecular dynamics (MD) simulation as Mg²⁺ is pulled from RNA, which causes disturbance of the water network. It was found that Mg²⁺ remains hexahydrated unless it is close to or far from RNA. In the pentahydrated form, Mg²⁺ interacts directly with RNA. Water network has been identified by two complimentary methods; MD followed by a density-based clustering algorithm and three-dimensional-reference interaction site model. These two methods gave similar results. Identification of water network around Mg²⁺ and non-WC base pairs gives a clue to the strong effect of water network on the stability of this RNA. Based on sequence analysis of all Eubacteria 5s rRNA, we propose that hexahydrated Mg²⁺ is an integral part of this RNA and geometry of base pairs surrounding it adjust to accommodate the [Formula: see text]. Overall the findings from this work can help in understanding the basis of the complex structure and stability of RNA with non-WC base pairs.

Are Waters around RNA More than Just a Solvent? - An Insight from Molecular Dynamics Simulations

Hydrating water molecules are believed to be an inherent part of the RNA structure and have a considerable impact on RNA conformation. However, the magnitude and mechanism of the interplay between water molecules and the RNA structure are still poorly understood. In principle, such hydration effects can be studied by molecular dynamics (MD) simulations. In our recent MD studies, we observed that the choice of water model has a visible impact on the predicted structure and structural dynamics of RNA and, in particular, has a larger effect than type, parametrization, and concentration of the ions. Furthermore, the water model effect is sequence dependent and modulates the sequence dependence of A-RNA helical parameters. Clearly, the sensitivity of A-RNA structural dynamics to the water model parametrization is a rather spurious effect that complicates MD studies of RNA molecules. These results nevertheless suggest that the sequence dependence of the A-RNA structure, usually attributed to base stacking, might be driven by the structural dynamics of specific hydration. Here, we present a systematic MD study that aimed to (i) clarify the atomistic mechanism of the water model sensitivity and (ii) discover whether and to what extent specific hydration modulates the A-RNA structural variability. We carried out an extended set of MD simulations of canonical A-RNA duplexes with TIP3P, TIP4P/2005, TIP5P, and SPC/E explicit water models and found that different water models provided a different extent of water bridging between 2'-OH groups across the minor groove, which in turn influences their distance and consequently also inclination, roll, and slide parameters. Minor groove hydration is also responsible for the sequence dependence of these helical parameters. Our simulations suggest that TIP5P is not optimal for RNA simulations.

Biophysical Characterization of G-Quadruplex Recognition in the PITX1 mRNA by the Specificity Domain of the Helicase RHAU

Nucleic acids rich in guanine are able to fold into unique structures known as G-quadruplexes. G-quadruplexes consist of four tracts of guanylates arranged in parallel or antiparallel strands that are aligned in stacked G-quartet planes. The structure is further stabilized by Hoogsteen hydrogen bonds and monovalent cations centered between the planes. RHAU (RNA helicase associated with AU-rich element) is a member of the ATP-dependent DExH/D family of RNA helicases and can bind and resolve G-quadruplexes. RHAU contains a core helicase domain with an N-terminal extension that enables recognition and full binding affinity to RNA and DNA G-quadruplexes. PITX1, a member of the bicoid class of homeobox proteins, is a transcriptional activator active during development of vertebrates, chiefly in the anterior pituitary gland and several other organs. We have previously demonstrated that RHAU regulates PITX1 levels through interaction with G-quadruplexes at the 3’-end of the PITX1 mRNA. To understand the structural basis of G-quadruplex recognition by RHAU, we characterize a purified minimal PITX1 G-quadruplex using a variety of biophysical techniques including electrophoretic mobility shift assays, UV-VIS spectroscopy, circular dichroism, dynamic light scattering, small angle X-ray scattering and nuclear magnetic resonance spectroscopy. Our biophysical analysis provides evidence that the RNA G-quadruplex, but not its DNA counterpart, can adopt a parallel orientation, and that only the RNA can interact with N-terminal domain of RHAU via the tetrad face of the G-quadruplex. This work extends our insight into how the N-terminal region of RHAU recognizes parallel G-quadruplexes.

Locked 2'-Deoxy-2',4'-Difluororibo Modified Nucleic Acids: Thermal Stability, Structural Studies, and siRNA Activity

2'-Deoxy-2',4'-difluorouridine (2',4'-diF-rU) was readily incorporated into DNA and RNA oligonucleotides via standard solid phase synthesis protocols. NMR and thermal denaturation (Tm) data of duplexes was consistent with the 2',4'-diF-rU nucleotides adopting a rigid North (RNA-like) sugar conformation, as previously observed for the nucleoside monomer. The impact of this modification on Tm is neutral when incorporated within RNA:RNA duplexes, mildly destabilizing when located in the RNA strand of a DNA:RNA duplex, and highly destabilizing when inserted in the DNA strand of DNA:RNA and DNA:DNA duplexes. Molecular dynamics calculations suggest that the destabilization effect in DNA:DNA and DNA:RNA duplexes is the result of structural distortions created by A/B junctions within the helical structures. Quantum mechanics calculations suggest that the "neutral" effect imparted to A-form duplexes is caused by alterations in charge distribution that compensate the stabilizing effect expected for a pure North-puckered furanose sugar. 2',4'-diF-RNA modified siRNAs were able to trigger RNA interference with excellent efficiency. Of note, incorporation of a few 2',4'-diF-rU residues in the middle of the guide (antisense) strand afforded siRNAs that were more potent than the corresponding siRNAs containing LNA and 2'-F-ANA modifications, and as active as the 2'-F-RNA modified siRNAs.

Mechanistic insights into temperature-dependent regulation of the simple cyanobacterial hsp17 RNA thermometer at base-pair resolution

The cyanobacterial hsp17 ribonucleicacid thermometer (RNAT) is one of the smallest naturally occurring RNAT. It forms a single hairpin with an internal 1×3-bulge separating the start codon in stem I from the ribosome binding site (RBS) in stem II. We investigated the temperature-dependent regulation of hsp17 by mapping individual base-pair stabilities from solvent exchange nuclear magnetic resonance (NMR) spectroscopy. The wild-type RNAT was found to be stabilized by two critical CG base pairs (C14-G27 and C13-G28). Replacing the internal 1×3 bulge by a stable CG base pair in hsp17^rep significantly increased the global stability and unfolding cooperativity as evidenced by circular dichroism spectroscopy. From the NMR analysis, remote stabilization and non-nearest neighbour effects exist at the base-pair level, in particular for nucleotide G28 (five nucleotides apart from the side of mutation). Individual base-pair stabilities are coupled to the stability of the entire thermometer within both the natural and the stabilized RNATs by enthalpy–entropy compensation presumably mediated by the hydration shell. At the melting point the Gibbs energies of the individual nucleobases are equalized suggesting a consecutive zipper-type unfolding mechanism of the RBS leading to a dimmer-like function of hsp17 and switch-like regulation behaviour of hsp17^rep. The data show how minor changes in the nucleotide sequence not only offset the melting temperature but also alter the mode of temperature sensing. The cyanobacterial thermosensor demonstrates the remarkable adjustment of natural RNATs to execute precise temperature control.

The cellular environment of cancerous human tissue. Interfacial and dangling water as a "hydration fingerprint"

Despite a large number of publications, the role of water in the cellular environment of biological tissue has not been clarified. Characterizing the biological interface is a key challenge in understanding the interactions of water in the tissue. Although we often assume that the properties of the bulk water can be translated to the crowded biological environment, this approach must be considerably revised when considering the biological interface. To our knowledge, few studies have directly monitored the interactions and accumulation of water in the restricted environments of the biological tissue upon realistic crowding conditions. The present study focuses on a molecular picture of water molecules at the biological interface, or specifically, water molecules adjacent to the hydrophobic and hydrophilic surfaces of normal and cancerous tissues. We recorded and analyzed the IR and Raman spectra of the νs(OH) stretching modes of water at the biological interfaces of the human breast and neck tissues. The results revealed dramatic changes in the water content in the tissue and are potentially relevant to both the fundamental problems of interfacial water modeling and the molecular diagnostics of cancer as a 'hydration fingerprint'. Herein, we will discuss the origin of the vibrational substructures observed for the νs(OH) stretching modes of water, showing that the interfacial water interacting via H-bond with other water molecules and biomolecules at the biological surface and free OH vibration of the dangling water are sensitive indicators of the pathology between the normal (noncancerous) and cancerous tissue and cancer types.

Hydration of protein-RNA recognition sites

We investigate the role of water molecules in 89 protein–RNA complexes taken from the Protein Data Bank. Those with tRNA and single-stranded RNA are less hydrated than with duplex or ribosomal proteins. Protein–RNA interfaces are hydrated less than protein–DNA interfaces, but more than protein–protein interfaces. Majority of the waters at protein–RNA interfaces makes multiple H-bonds; however, a fraction do not make any. Those making H-bonds have preferences for the polar groups of RNA than its partner protein. The spatial distribution of waters makes interfaces with ribosomal proteins and single-stranded RNA relatively ‘dry’ than interfaces with tRNA and duplex RNA. In contrast to protein–DNA interfaces, mainly due to the presence of the 2′OH, the ribose in protein–RNA interfaces is hydrated more than the phosphate or the bases. The minor groove in protein–RNA interfaces is hydrated more than the major groove, while in protein–DNA interfaces it is reverse. The strands make the highest number of water-mediated H-bonds per unit interface area followed by the helices and the non-regular structures. The preserved waters at protein–RNA interfaces make higher number of H-bonds than the other waters. Preserved waters contribute toward the affinity in protein–RNA recognition and should be carefully treated while engineering protein–RNA interfaces.

Amides are excellent mimics of phosphate internucleoside linkages and are well tolerated in short interfering RNAs

RNA interference (RNAi) has become an important tool in functional genomics and has an intriguing therapeutic potential. However, the current design of short interfering RNAs (siRNAs) is not optimal for in vivo applications. Non-ionic phosphate backbone modifications may have the potential to improve the properties of siRNAs, but are little explored in RNAi technologies. Using X-ray crystallography and RNAi activity assays, the present study demonstrates that 3'-CH2-CO-NH-5' amides are excellent replacements for phosphodiester internucleoside linkages in RNA. The crystal structure shows that amide-modified RNA forms a typical A-form duplex. The amide carbonyl group points into the major groove and assumes an orientation that is similar to the P-OP2 bond in the phosphate linkage. Amide linkages are well hydrated by tandem waters linking the carbonyl group and adjacent phosphate oxygens. Amides are tolerated at internal positions of both the guide and passenger strand of siRNAs and may increase the silencing activity when placed near the 5'-end of the passenger strand. As a result, an siRNA containing eight amide linkages is more active than the unmodified control. The results suggest that RNAi may tolerate even more extensive amide modification, which may be useful for optimization of siRNAs for in vivo applications.

Solid-phase synthesis and hybrization behavior of partially 2'/3'-O-acetylated RNA oligonucleotides

Synthesis of partially 2'/3'-O-acetylated oligoribonucleotides has been accomplished by using a 2'/3'-O-acetyl orthogonal protecting group strategy in which non-nucleophilic strong-base (DBU) labile nucleobase protecting groups and a UV-light cleavable linker were used. Strong-base stability of the photolabile linker allowed on-column nucleobase and phosphate deprotection, followed by a mild cleavage of the acetylated oligonucleotides from the solid support with UV light. Two 17nt oligonucleotides, which were synthesized possessing one specific internal 2'- or 3'-acetyl group, were used as synthetic standards in a recent report from this laboratory detailing the prebiotically plausible ligation of RNA oligonucleotides. In order to further investigate the effect of 2'/3'-O-acetyl groups on the stability of RNA duplex structure, two complementary bis-acetylated RNA oligonucleotides were also expediently obtained with the newly developed protocols. UV melting curves of 2'-O-acetylated RNA duplexes showed a consistent ~3.1 °C decrease in Tm per 2'-O-acetyl group.