The RNA backbone plays a crucial role in mediating the intrinsic stability of the GpU dinucleotide platform and the GpUpA/GpA miniduplex

Dissecting the energetic architecture within an RNA tertiary structural motif via high-throughput thermodynamic measurements

Structured RNAs and RNA/protein complexes perform critical cellular functions. They often contain structurally conserved tertiary contact "motifs," whose occurrence simplifies the RNA folding landscape. Prior studies have focused on the conformational and energetic modularity of intact motifs. Here, we turn to the dissection of one common motif, the 11nt receptor (11ntR), using quantitative analysis of RNA on a massively parallel array to measure the binding of all single and double 11ntR mutants to GAAA and GUAA tetraloops, thereby probing the energetic architecture of the motif. While the 11ntR behaves as a motif, its cooperativity is not absolute. Instead, we uncovered a gradient from high cooperativity amongst base-paired and neighboring residues to additivity between distant residues. As expected, substitutions at residues in direct contact with the GAAA tetraloop resulted in the largest decreases to binding, and energetic penalties of mutations were substantially smaller for binding to the alternate GUAA tetraloop, which lacks tertiary contacts present with the canonical GAAA tetraloop. However, we found that the energetic consequences of base partner substitutions are not, in general, simply described by base pair type or isostericity. We also found exceptions to the previously established stability-abundance relationship for 11ntR sequence variants. These findings of "exceptions to the rule" highlight the power of systematic high-throughput approaches to uncover novel variants for future study in addition to providing an energetic map of a functional RNA.

Single molecule force spectroscopy of a streptomycin-binding RNA aptamer: An out-of-equilibrium molecular dynamics study

Here, we investigate the unfolding behavior of a streptomycin-binding ribonucleic acid (RNA) aptamer under application of force in shear geometry. Using Langevin out-of-equilibrium simulations to emulate the single-molecule force spectroscopy (SMFS) experiment, we were able to understand the hierarchical unfolding process that occurs in the RNA molecule under application of stretching force and the influence of streptomycin modifying this unfolding. Subsequently, the application of the Jarzynski equality to the force profiles obtained in the pulling simulations shows that the free energies for individual systems and the difference of unfolding free energy upon streptomycin binding to the RNA free aptamer are in fair agreement with the experimental values, obtained through SMFS by Nick et al. [J. Phys. Chem. B 120, 6479 (2016)].

In silico survey of the central conserved regions in viroids of the Pospiviroidae family for conserved asymmetric loop structures

Viroids are the smallest replicative pathogens, consisting of RNA circles (∼300 nucleotides) that require host machinery to replicate. Structural RNA elements recruit these host factors. Currently, many of these structural elements and the nature of their interactions are unknown. All Pospiviroidae have homology in the central conserved region (CCR). The CCR of potato spindle tuber viroid (PSTVd) contains a sarcin/ricin domain (SRD), the only viroid structural element with an unequivocal replication role. We assumed that every member of this family uses this region to recruit host factors, and that each CCR has an SRD-like asymmetric loop within it. Potential SRD or SRD-like motifs were sought in the CCR of each Pospiviroidae member as follows. Motif location in each CCR was predicted with MUSCLE alignment and Vienna RNAfold. Viroid-specific models of SRD-like motifs were built by superimposing noncanonical base pairs and nucleotides on a model of an SRD. The RNA geometry search engine FR3D was then used to find nucleotide groups close to the geometry suggested by this superimposition. Atomic resolution structures were assembled using the molecular visualization program Chimera, and the stability of each motif was assessed with molecular dynamics (MD). Some models required a protonated cytosine. To be stable within a cell, the pK_a of that cytosine must be shifted up. Constant pH-replica exchange MD analysis showed such a shift in the proposed structures. These data show that every Pospiviroidae member could form a motif that resembles an SRD in its CCR, and imply there could be undiscovered mimics of other RNA domains.

Effects of Noncanonical Base Pairing on RNA Folding: Structural Context and Spatial Arrangements of G·A Pairs

Noncanonical base pairs play important roles in assembling the three-dimensional structures critical to the diverse functions of RNA. These associations contribute to the looped segments that intersperse the canonical double-helical elements within folded, globular RNA molecules. They stitch together various structural elements, serve as recognition elements for other molecules, and act as sites of intrinsic stiffness or deformability. This work takes advantage of new software (DSSR) designed to streamline the analysis and annotation of RNA three-dimensional structures. The multiscale structural information gathered for individual molecules, combined with the growing number of unique, well-resolved RNA structures, makes it possible to examine the collective features deeply and to uncover previously unrecognized patterns of chain organization. Here we focus on a subset of noncanonical base pairs involving guanine and adenine and the links between their modes of association, secondary structural context, and contributions to tertiary folding. The rigorous descriptions of base-pair geometry that we employ facilitate characterization of recurrent geometric motifs and the structural settings in which these arrangements occur. Moreover, the numerical parameters hint at the natural motions of the interacting bases and the pathways likely to connect different spatial forms. We draw attention to higher-order multiplexes involving two or more G·A pairs and the roles these associations appear to play in bridging different secondary structural units. The collective data reveal pairing propensities in base organization, secondary structural context, and deformability and serve as a starting point for further multiscale investigations and/or simulations of RNA folding.

Going beyond base-pairs: topology-based characterization of base-multiplets in RNA

Identification and characterization of base-multiplets, which are essentially mediated by base-pairing interactions, can provide insights into the diversity in the structure and dynamics of complex functional RNAs, and thus facilitate hypothesis driven biological research. The necessary nomenclature scheme, an extension of the geometric classification scheme for base-pairs by Leontis and Westhof, is however available only for base-triplets. In the absence of information on topology, this scheme is not applicable to quartets and higher order multiplets. Here we propose a topology-based classification scheme which, in conjunction with a graph-based algorithm, can be used for the automated identification and characterization of higher order base-multiplets in RNA structures. Here, the RNA structure is represented as a graph, where nodes represent nucleotides and edges represent base-pairing connectivity. Sets of connected components (of n nodes) within these graphs constitute subgraphs representing multiplets of "n" nucleotides. The different topological variants of the RNA multiplets thus correspond to different nonisomorphic forms of these subgraphs. To annotate RNA base-multiplets unambiguously, we propose a set of topology-based nomenclature rules for quartets, which are extendable to higher multiplets. We also demonstrate the utility of our approach toward the identification and annotation of higher order RNA multiplets, by investigating the occurrence contexts of selected examples in order to gain insights regarding their probable functional roles.

Identification of D Modification Sites by Integrating Heterogeneous Features in Saccharomyces cerevisiae

As an abundant post-transcriptional modification, dihydrouridine (D) has been found in transfer RNA (tRNA) from bacteria, eukaryotes, and archaea. Nonetheless, knowledge of the exact biochemical roles of dihydrouridine in mediating tRNA function is still limited. Accurate identification of the position of D sites is essential for understanding their functions. Therefore, it is desirable to develop novel methods to identify D sites. In this study, an ensemble classifier was proposed for the detection of D modification sites in the Saccharomyces cerevisiae transcriptome by using heterogeneous features. The jackknife test results demonstrate that the proposed predictor is promising for the identification of D modification sites. It is anticipated that the proposed method can be widely used for identifying D modification sites in tRNA.

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.

All-atom MD indicates ion-dependent behavior of therapeutic DNA polymer

Understanding the efficacy of and creating delivery mechanisms for therapeutic nucleic acids requires understanding structural and kinetic properties which allow these polymers to promote the death of cancerous cells. One molecule of interest is a 10 mer of FdUMP (5-fluoro-2'-deoxyuridine-5'-O-monophosphate) - also called F10. Here we investigate the structural and kinetic behavior of F10 in intracellular and extracellular solvent conditions along with non-biological conditions that may be efficacious in in vitro preparations of F10 delivery systems. From our all-atom molecular dynamics simulations totaling 80 microseconds, we predict that F10's phosphate groups form close-range interactions with calcium and zinc ions, with calcium having the highest affinity of the five ions investigated. We also predict that F10's interactions with magnesium, potassium and sodium are almost exclusively long-range interactions. In terms of intramolecular interactions, we find that F10 is least structured (in terms of hydrogen bonds among bases) in the 150 mM NaCl (extracellular-like solvent conditions) and most structured in 150 mM ZnCl2. Kinetically, we see that F10 is unstable in the presence of magnesium, sodium or potassium, finding stable kinetic traps in the presence of calcium or zinc.

All-Atom MD Predicts Magnesium-Induced Hairpin in Chemically Perturbed RNA Analog of F10 Therapeutic

Given their increasingly frequent usage, understanding the chemical and structural properties which allow therapeutic nucleic acids to promote the death of cancer cells is critical for medical advancement. One molecule of interest is a 10-mer of FdUMP (5-fluoro-2'-deoxyuridine-5'-O-monophosphate) also called F10. To investigate causes of structural stability, we have computationally restored the 2' oxygen on each ribose sugar of the phosphodiester backbone, creating FUMP[10]. Microsecond time-scale, all-atom, simulations of FUMP[10] in the presence of 150 mM MgCl₂ predict that the strand has a 45% probability of folding into a stable hairpin-like secondary structure. Analysis of 16 μs of data reveals phosphate interactions as likely contributors to the stability of this folded state. Comparison with polydT and polyU simulations predicts that FUMP[10]'s lowest order structures last for one to 2 orders of magnitude longer than similar nucleic acid strands. Here we provide a brief structural and conformational analysis of the predicted structures of FUMP[10], and suggest insights into its stability via comparison to F10, polydT, and polyU.

All-Atom Molecular Dynamics Reveals Mechanism of Zinc Complexation with Therapeutic F10

Advancing the use of therapeutic nucleic acids requires understanding the chemical and structural properties that allow these polymers to promote the death of malignant cells. Here we explore Zn²⁺ complexation by the fluoropyrimidine polymer F10, which has strong activities in multiple preclinical models of cancer. Delivery of fluoropyrimidine FdUMP in the 10-residue polymer F10 rather than the nucleobase (5-fluorouracil) allows consideration of metal ion binding effects on drug delivery. The differences in metal ion interactions with fluoropyrimidine compared to normal DNA results in conformation changes that affect protein binding, cell uptake, and codelivery of metals such as zinc, and the cytoxicity thereof. Microsecond-time-scale, all-atom simulations of F10 predict that zinc selectively stabilizes the polymer via interactions with backbone phosphate groups and suggest a mechanism of complexation for the zinc-base interactions shown in previous experimental work. The positive zinc ions are attracted to the negatively charged phosphate groups. Once the Zn²⁺ ions are near F10, they cause the base's N3 nitrogen to deprotonate. Subsequently, magnesium atoms displace zinc from their interactions with phosphate, freeing the zinc ions to interact with the FdU bases by forming weak interactions with the O4 oxygen and the fluorine attached to C5. These interactions of magnesium with phosphate groups and zinc with nucleobases agree with previous experimental results and are seen in MD simulations only when magnesium is introduced after N3 deprotonation, indicating a specific order of metal binding events. Additionally, we predict interactions between zinc and F10's O2 atoms, which were not previously observed. By comparison to 10mers of polyU and polydT, we also predict that the presence of fluorine increases the binding affinity of zinc to F10 relative to analogous strands of RNA and DNA consisting of only native nucleotides.

Quantum chemical benchmark study on 46 RNA backbone families using a dinucleotide unit

We have created a benchmark set of quantum chemical structure-energy data denoted as UpU46, which consists of 46 uracil dinucleotides (UpU), representing all known 46 RNA backbone conformational families. Penalty-function-based restrained optimizations with COSMO TPSS-D3/def2-TZVP ensure a balance between keeping the target conformation and geometry relaxation. The backbone geometries are close to the clustering-means of their respective RNA bioinformatics family classification. High-level wave function methods (DLPNO-CCSD(T) as reference) and a wide-range of dispersion-corrected or inclusive DFT methods (DFT-D3, VV10, LC-BOP-LRD, M06-2X, M11, and more) are used to evaluate the conformational energies. The results are compared to the Amber RNA bsc0χOL3 force field. Most dispersion-corrected DFT methods surpass the Amber force field significantly in accuracy and yield mean absolute deviations (MADs) for relative conformational energies of ∼0.4-0.6 kcal/mol. Double-hybrid density functionals represent the most accurate class of density functionals. Low-cost quantum chemical methods such as PM6-D3H+, HF-3c, DFTB3-D3, as well as small basis set calculations corrected for basis set superposition errors (BSSEs) by the gCP procedure are also tested. Unfortunately, the presently available low-cost methods are struggling to describe the UpU conformational energies with satisfactory accuracy. The UpU46 benchmark is an ideal test for benchmarking and development of fast methods to describe nucleic acids, including force fields.

Identification and analysis of the N(6)-methyladenosine in the Saccharomyces cerevisiae transcriptome

Knowledge of the distribution of N⁶-methyladenosine (m⁶A) is invaluable for understanding RNA biological functions. However, limitation in experimental methods impedes the progress towards the identification of m⁶A site. As a complement of experimental methods, a support vector machine based-method is proposed to identify m⁶A sites in Saccharomyces cerevisiae genome. In this model, RNA sequences are encoded by their nucleotide chemical property and accumulated nucleotide frequency information. It is observed in the jackknife test that the accuracy achieved by the proposed model in identifying the m⁶A site was 78.15%. For the convenience of experimental scientists, a web-server for the proposed model is provided at http://lin.uestc.edu.cn/server/m6Apred.php.

Prediction of hydrogen and carbon chemical shifts from RNA using database mining and support vector regression

The Biological Magnetic Resonance Data Bank (BMRB) contains NMR chemical shift depositions for over 200 RNAs and RNA-containing complexes. We have analyzed the (1)H NMR and (13)C chemical shifts reported for non-exchangeable protons of 187 of these RNAs. Software was developed that downloads BMRB datasets and corresponding PDB structure files, and then generates residue-specific attributes based on the calculated secondary structure. Attributes represent properties present in each sequential stretch of five adjacent residues and include variables such as nucleotide type, base-pair presence and type, and tetraloop types. Attributes and (1)H and (13)C NMR chemical shifts of the central nucleotide are then used as input to train a predictive model using support vector regression. These models can then be used to predict shifts for new sequences. The new software tools, available as stand-alone scripts or integrated into the NMR visualization and analysis program NMRViewJ, should facilitate NMR assignment and/or validation of RNA (1)H and (13)C chemical shifts. In addition, our findings enabled the re-calibration a ring-current shift model using published NMR chemical shifts and high-resolution X-ray structural data as guides.

DSSR: an integrated software tool for dissecting the spatial structure of RNA

Insight into the three-dimensional architecture of RNA is essential for understanding its cellular functions. However, even the classic transfer RNA structure contains features that are overlooked by existing bioinformatics tools. Here we present DSSR (Dissecting the Spatial Structure of RNA), an integrated and automated tool for analyzing and annotating RNA tertiary structures. The software identifies canonical and noncanonical base pairs, including those with modified nucleotides, in any tautomeric or protonation state. DSSR detects higher-order coplanar base associations, termed multiplets. It finds arrays of stacked pairs, classifies them by base-pair identity and backbone connectivity, and distinguishes a stem of covalently connected canonical pairs from a helix of stacked pairs of arbitrary type/linkage. DSSR identifies coaxial stacking of multiple stems within a single helix and lists isolated canonical pairs that lie outside of a stem. The program characterizes 'closed' loops of various types (hairpin, bulge, internal, and junction loops) and pseudoknots of arbitrary complexity. Notably, DSSR employs isolated pairs and the ends of stems, whether pseudoknotted or not, to define junction loops. This new, inclusive definition provides a novel perspective on the spatial organization of RNA. Tests on all nucleic acid structures in the Protein Data Bank confirm the efficiency and robustness of the software, and applications to representative RNA molecules illustrate its unique features. DSSR and related materials are freely available at http://x3dna.org/.

An atlas of RNA base pairs involving modified nucleobases with optimal geometries and accurate energies

Posttranscriptional modifications greatly enhance the chemical information of RNA molecules, contributing to explain the diversity of their structures and functions. A significant fraction of RNA experimental structures available to date present modified nucleobases, with half of them being involved in H-bonding interactions with other bases, i.e. ‘modified base pairs’. Herein we present a systematic investigation of modified base pairs, in the context of experimental RNA structures. To this end, we first compiled an atlas of experimentally observed modified base pairs, for which we recorded occurrences and structural context. Then, for each base pair, we selected a representative for subsequent quantum mechanics calculations, to find out its optimal geometry and interaction energy. Our structural analyses show that most of the modified base pairs are non Watson–Crick like and are involved in RNA tertiary structure motifs. In addition, quantum mechanics calculations quantify and provide a rationale for the impact of the different modifications on the geometry and stability of the base pairs they participate in.

MINT: software to identify motifs and short-range interactions in trajectories of nucleic acids

Structural biology experiments and structure prediction tools have provided many high-resolution three-dimensional structures of nucleic acids. Also, molecular dynamics force field parameters have been adapted to simulating charged and flexible nucleic acid structures on microsecond time scales. Therefore, we can generate the dynamics of DNA or RNA molecules, but we still lack adequate tools for the analysis of the resulting huge amounts of data. We present MINT (Motif Identifier for Nucleic acids Trajectory) — an automatic tool for analyzing three-dimensional structures of RNA and DNA, and their full-atom molecular dynamics trajectories or other conformation sets (e.g. X-ray or nuclear magnetic resonance-derived structures). For each RNA or DNA conformation MINT determines the hydrogen bonding network resolving the base pairing patterns, identifies secondary structure motifs (helices, junctions, loops, etc.) and pseudoknots. MINT also estimates the energy of stacking and phosphate anion-base interactions. For many conformations, as in a molecular dynamics trajectory, MINT provides averages of the above structural and energetic features and their evolution. We show MINT functionality based on all-atom explicit solvent molecular dynamics trajectory of the 30S ribosomal subunit.

Chemical probing of RNA with the hydroxyl radical at single-atom resolution

While hydroxyl radical cleavage is widely used to map RNA tertiary structure, lack of mechanistic understanding of strand break formation limits the degree of structural insight that can be obtained from this experiment. Here, we determine how individual ribose hydrogens of sarcin/ricin loop RNA participate in strand cleavage. We find that substituting deuterium for hydrogen at a ribose 5'-carbon produces a kinetic isotope effect on cleavage; the major cleavage product is an RNA strand terminated by a 5'-aldehyde. We conclude that hydroxyl radical abstracts a 5'-hydrogen atom, leading to RNA strand cleavage. We used this approach to obtain structural information for a GUA base triple, a common tertiary structural feature of RNA. Cleavage at U exhibits a large 5' deuterium kinetic isotope effect, a potential signature of a base triple. Others had noted a ribose-phosphate hydrogen bond involving the G 2'-OH and the U phosphate of the GUA triple, and suggested that this hydrogen bond contributes to backbone rigidity. Substituting deoxyguanosine for G, to eliminate this hydrogen bond, results in a substantial decrease in cleavage at G and U of the triple. We conclude that this hydrogen bond is a linchpin of backbone structure around the triple.

Structure of the complete bacterial SRP Alu domain

The Alu domain of the signal recognition particle (SRP) arrests protein biosynthesis by competition with elongation factor binding on the ribosome. The mammalian Alu domain is a protein-RNA complex, while prokaryotic Alu domains are protein-free with significant extensions of the RNA. Here we report the crystal structure of the complete Alu domain of Bacillus subtilis SRP RNA at 2.5 Å resolution. The bacterial Alu RNA reveals a compact fold, which is stabilized by prokaryote-specific extensions and interactions. In this 'closed' conformation, the 5' and 3' regions are clamped together by the additional helix 1, the connecting 3-way junction and a novel minor groove interaction, which we term the 'minor-saddle motif' (MSM). The 5' region includes an extended loop-loop pseudoknot made of five consecutive Watson-Crick base pairs. Homology modeling with the human Alu domain in context of the ribosome shows that an additional lobe in the pseudoknot approaches the large subunit, while the absence of protein results in the detachment from the small subunit. Our findings provide the structural basis for purely RNA-driven elongation arrest in prokaryotes, and give insights into the structural adaption of SRP RNA during evolution.

Sequence dependent variations in RNA duplex are related to non-canonical hydrogen bond interactions in dinucleotide steps

Background

Sequence determines the three-dimensional structure of RNAs, and thereby plays an important role in carrying out various biological functions. RNA duplexes containing Watson-Crick (WC) basepairs, interspersed with non-Watson-Crick basepairs, are the dominant structural unit and form the scaffold for the 3-dimensional structure of RNA. It is therefore crucial to understand the geometric variation in the dinucleotide steps that form the helices. We have carried out a detailed analysis of the dinucleotide steps formed by AU and GC Watson-Crick basepairs in RNA structures (both free and protein bound) and compared the results to that seen in DNA. Further, the effect of protein binding on these steps was examined by comparing steps in free RNA structures with protein bound RNA structures.

Results

Characteristic sequence dependent geometries are observed for the RR, RY and YR type of dinucleotide steps in RNA. Their geometric parameters show correlated variations that are different from those observed in B-DNA helices. Subtle, but statistically significant differences are seen in roll, slide and average propeller-twist values, between the dinucleotide steps of free RNA and protein bound RNA structures. Many non-canonical cross-strand and intra-strand hydrogen bonds were identified that can stabilise the RNA dinucleotide steps, among which YR steps show presence of many new unreported interactions.

Conclusions

Our work provides for the first time a detailed analysis of the conformational preferences exhibited by Watson-Crick basepair containing steps in RNA double helices. Overall, the WC dinucleotide steps show considerable conformational variability. Furthermore, we have identified hydrogen bond interactions in several of the dinucleotide steps that could play a role in determining the preferred geometry, in addition to the intra-basepair hydrogen bonds and stacking interactions. Protein binding affects the conformation of the steps that are in direct contact, as well as allosterically affect the steps that are not in direct physical contact.

Crystal structure of metallo DNA duplex containing consecutive Watson-Crick-like T-Hg(II)-T base pairs

The metallo DNA duplex containing mercury-mediated T-T base pairs is an attractive biomacromolecular nanomaterial which can be applied to nanodevices such as ion sensors. Reported herein is the first crystal structure of a B-form DNA duplex containing two consecutive T-Hg(II)-T base pairs. The Hg(II) ion occupies the center between two T residues. The N3-Hg(II) bond distance is 2.0 Å. The relatively short Hg(II)-Hg(II) distance (3.3 Å) observed in consecutive T-Hg(II)-T base pairs suggests that the metallophilic attraction could exist between them and may stabilize the B-form double helix. To support this, the DNA duplex is largely distorted and adopts an unusual nonhelical conformation in the absence of Hg(II). The structure of the metallo DNA duplex itself and the Hg(II)-induced structural switching from the nonhelical form to the B-form provide the basis for structure-based design of metal-conjugated nucleic acid nanomaterials.

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

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

Isosteric and nonisosteric base pairs in RNA motifs: molecular dynamics and bioinformatics study of the sarcin-ricin internal loop

The sarcin-ricin RNA motif (SR motif) is one of the most prominent recurrent RNA building blocks that occurs in many different RNA contexts and folds autonomously, that is, in a context-independent manner. In this study, we combined bioinformatics analysis with explicit-solvent molecular dynamics (MD) simulations to better understand the relation between the RNA sequence and the evolutionary patterns of the SR motif. A SHAPE probing experiment was also performed to confirm the fidelity of the MD simulations. We identified 57 instances of the SR motif in a nonredundant subset of the RNA X-ray structure database and analyzed their base pairing, base-phosphate, and backbone-backbone interactions. We extracted sequences aligned to these instances from large rRNA alignments to determine the frequency of occurrence for different sequence variants. We then used a simple scoring scheme based on isostericity to suggest 10 sequence variants with a highly variable expected degree of compatibility with the SR motif 3D structure. We carried out MD simulations of SR motifs with these base substitutions. Nonisosteric base substitutions led to unstable structures, but so did isosteric substitutions which were unable to make key base-phosphate interactions. The MD technique explains why some potentially isosteric SR motifs are not realized during evolution. We also found that the inability to form stable cWW geometry is an important factor in the case of the first base pair of the flexible region of the SR motif. A comparison of structural, bioinformatics, SHAPE probing, and MD simulation data reveals that explicit solvent MD simulations neatly reflect the viability of different sequence variants of the SR motif. Thus, MD simulations can efficiently complement bioinformatics tools in studies of conservation patterns of RNA motifs and provide atomistic insight into the role of their different signature interactions.

Analyzing and building nucleic acid structures with 3DNA

The 3DNA software package is a popular and versatile bioinformatics tool with capabilities to analyze, construct, and visualize three-dimensional nucleic acid structures. This article presents detailed protocols for a subset of new and popular features available in 3DNA, applicable to both individual structures and ensembles of related structures. Protocol 1 lists the set of instructions needed to download and install the software. This is followed, in Protocol 2, by the analysis of a nucleic acid structure, including the assignment of base pairs and the determination of rigid-body parameters that describe the structure and, in Protocol 3, by a description of the reconstruction of an atomic model of a structure from its rigid-body parameters. The most recent version of 3DNA, version 2.1, has new features for the analysis and manipulation of ensembles of structures, such as those deduced from nuclear magnetic resonance (NMR) measurements and molecular dynamic (MD) simulations; these features are presented in Protocols 4 and 5. In addition to the 3DNA stand-alone software package, the w3DNA web server, located at http://w3dna.rutgers.edu, provides a user-friendly interface to selected features of the software. Protocol 6 demonstrates a novel feature of the site for building models of long DNA molecules decorated with bound proteins at user-specified locations.

Molecular mechanism of preQ1 riboswitch action: a molecular dynamics study

Riboswitches often occur in the 5'-untranslated regions of bacterial mRNA where they regulate gene expression. The preQ(1) riboswitch controls the biosynthesis of a hypermodified nucleoside queuosine in response to binding the queuosine metabolic intermediate. Structures of the ligand-bound and ligand-free states of the preQ(1) riboswitch from Thermoanaerobacter tengcongensis were determined recently by X-ray crystallography. We used multiple, microsecond-long molecular dynamics simulations (29 μs in total) to characterize the structural dynamics of preQ(1) riboswitches in both states. We observed different stabilities of the stem in the bound and free states, resulting in different accessibilities of the ribosome-binding site. These differences are related to different stacking interactions between nucleotides of the stem and the associated loop, which itself adopts different conformations in the bound and free states. We suggest that the loop not only serves to bind preQ(1) but also transmits information about ligand binding from the ligand-binding pocket to the stem, which has implications for mRNA accessibility to the ribosome. We explain functional results obscured by a high salt crystallization medium and help to refine regions of disordered electron density, which demonstrates the predictive power of our approach. Besides investigating the functional dynamics of the riboswitch, we have also utilized this unique small folded RNA system for analysis of performance of the RNA force field on the μs time scale. The latest AMBER parmbsc0χ(OL3) RNA force field is capable of providing stable trajectories of the folded molecule on the μs time scale. On the other hand, force fields that are not properly balanced lead to significant structural perturbations on the sub-μs time scale, which could easily lead to inappropriate interpretation of the simulation data.

The DNA and RNA sugar-phosphate backbone emerges as the key player. An overview of quantum-chemical, structural biology and simulation studies

Knowledge of geometrical and physico-chemical properties of the sugar-phosphate backbone substantially contributes to the comprehension of the structural dynamics, function and evolution of nucleic acids. We provide a side by side overview of structural biology/bioinformatics, quantum chemical and molecular mechanical/simulation studies of the nucleic acids backbone. We highlight main features, advantages and limitations of these techniques, with a special emphasis given to their synergy. The present status of the research is then illustrated by selected examples which include classification of DNA and RNA backbone families, benchmark structure-energy quantum chemical calculations, parameterization of the dihedral space of simulation force fields, incorporation of arsenate into DNA, sugar-phosphate backbone self-cleavage in small RNA enzymes, and intricate geometries of the backbone in recurrent RNA building blocks. Although not apparent from the current literature showing limited overlaps between the QM, simulation and bioinformatics studies of the nucleic acids backbone, there in fact should be a major cooperative interaction between these three approaches in studies of the sugar-phosphate backbone.

Structure and energy of non-canonical basepairs: comparison of various computational chemistry methods with crystallographic ensembles

Different types of non-canonical basepairs, in addition to the Watson-Crick ones, are observed quite frequently in RNA. Their importance in the three dimensional structure is not fully understood, but their various roles have been proposed by different groups. We have analyzed the energetics and geometry of 32 most frequently observed basepairs in the functional RNA crystal structures using different popular empirical, semi-empirical and ab initio quantum chemical methods and compared their optimized geometry with the crystal data. These basepairs are classified into three categories: polar, non-polar and sugar-mediated, depending on the types of atoms involved in hydrogen bonding. In case of polar basepairs, most of the methods give rise to optimized structures close to their initial geometry. The interaction energies also follow similar trends, with the polar ones having more attractive interaction energies. Some of the C-H...O/N hydrogen bond mediated non-polar basepairs are also found to be significantly stable in terms of their interaction energy values. Few polar basepairs, having amino or carboxyl groups not hydrogen bonded to anything, such as G:G H:W C, show large flexibility. Most of the non-polar basepairs, except A:G s:s T and A:G w:s C, are found to be stable; indicating C-H...O/N interaction also plays a prominent role in stabilizing the basepairs. The sugar mediated basepairs show variability in their structures, due to the involvement of flexible ribose sugar. These presumably indicate that the most of the polar basepairs along with few non-polar ones act as seed for RNA folding while few may act as some conformational switch in the RNA.

Understanding the Sequence Preference of Recurrent RNA Building Blocks using Quantum Chemistry: The Intrastrand RNA Dinucleotide Platform

Folded RNA molecules are shaped by an astonishing variety of highly conserved noncanonical molecular interactions and backbone topologies. The dinucleotide platform is a widespread recurrent RNA modular building submotif formed by the side-by-side pairing of bases from two consecutive nucleotides within a single strand, with highly specific sequence preferences. This unique arrangement of bases is cemented by an intricate network of noncanonical hydrogen bonds and facilitated by a distinctive backbone topology. The present study investigates the gas-phase intrinsic stabilities of the three most common RNA dinucleotide platforms - 5'-GpU-3', ApA, and UpC - via state-of-the-art quantum-chemical (QM) techniques. The mean stability of base-base interactions decreases with sequence in the order GpU > ApA > UpC. Bader's atoms-in-molecules analysis reveals that the N2(G)…O4(U) hydrogen bond of the GpU platform is stronger than the corresponding hydrogen bonds in the other two platforms. The mixed-pucker sugar-phosphate backbone conformation found in most GpU platforms, in which the 5'-ribose sugar (G) is in the C2'-endo form and the 3'-sugar (U) in the C3'-endo form, is intrinsically more stable than the standard A-RNA backbone arrangement, partially as a result of a favorable O2'…O2P intra-platform interaction. Our results thus validate the hypothesis of Lu et al. (Lu Xiang-Jun, et al. Nucleic Acids Res. 2010, 38, 4868-4876), that the superior stability of GpU platforms is partially mediated by the strong O2'…O2P hydrogen bond. In contrast, ApA and especially UpC platform-compatible backbone conformations are rather diverse and do not display any characteristic structural features. The average stabilities of ApA and UpC derived backbone conformers are also lower than those of GpU platforms. Thus, the observed structural and evolutionary patterns of the dinucleotide platforms can be accounted for, to a large extent, by their intrinsic properties as described by modern QM calculations. In contrast, we show that the dinucleotide platform is not properly described in the course of atomistic explicit-solvent simulations. Our work also gives methodological insights into QM calculations of experimental RNA backbone geometries. Such calculations are inherently complicated by rather large data and refinement uncertainties in the available RNA experimental structures, which often preclude reliable energy computations.