RNA conformational classes

Application of constrained real-space refinement of flexible molecular fragments to automatic model building of RNA structures

New methods have been developed for locating phosphate groups and nucleic acid bases in the electron density of RNA structures. These methods utilize a constrained real-space refinement of molecular fragments and a phased rotation–conformation–translation function. Real-space refinement has also contributed to the improvement of the bone/base method of RNA model building and to redesigning the method of building double helices in nucleic acid structures. This improvement is reflected in the increased accuracy of the model building and the ability to better distinguish between correct and false solutions. A program,RSR, was created, and the programsNUT,HELandDHLwere upgraded and organized into a program system, which is CCP4 oriented. Source codes will also be released.

Discrete RNA libraries from pseudo-torsional space

The discovery that RNA molecules can fold into complex structures and carry out diverse cellular roles has led to interest in developing tools for modeling RNA tertiary structure. While significant progress has been made in establishing that the RNA backbone is rotameric, few libraries of discrete conformations specifically for use in RNA modeling have been validated. Here, we present six libraries of discrete RNA conformations based on a simplified pseudo-torsional notation of the RNA backbone, comparable to phi and psi in the protein backbone. We evaluate the ability of each library to represent single nucleotide backbone conformations, and we show how individual library fragments can be assembled into dinucleotides that are consistent with established RNA backbone descriptors spanning from sugar to sugar. We then use each library to build all-atom models of 20 test folds, and we show how the composition of a fragment library can limit model quality. Despite the limitations inherent in using discretized libraries, we find that several hundred discrete fragments can rebuild RNA folds up to 174 nucleotides in length with atomic-level accuracy (<1.5 Å RMSD). We anticipate that the libraries presented here could easily be incorporated into RNA structural modeling, analysis, or refinement tools.

Correlating the 31P NMR chemical shielding tensor and the 2J(P,C) spin-spin coupling constants with torsion angles ζ and α in the backbone of nucleic acids

Determination of nucleic acid (NA) structure with NMR spectroscopy is limited by the lack of restraints on conformation of NA phosphate. In this work, the (31)P chemical shielding tensor, the Γ(P,C5'H5'1) and Γ(P,C5'H5'2) cross-correlated relaxation rates, and the (2)J(P,C3'), (2)J(P,C5'), and (3)J(P,C4') coupling constants were calculated in dependence on NA backbone torsion angles ζ and α. While the orientation of the (31)P chemical shielding tensor was almost independent of the NA phosphate conformation, the principal tensor components varied by up to ~40 ppm. This variation and the dependence of the phosphate geometry on torsion angles ζ and α had only a minor influence on the calculated Γ(P,C5'H5'1) and Γ(P,C5'H5'2) cross-correlated relaxation rates, and therefore, the so-called rigid tensor approximation was here validated. For the first time, the (2)J(P,C) spin-spin coupling constants were correlated with the conformation of NA phosphate. Although each of the two J-couplings was significantly modulated by both torsions ζ and α, the (2)J(P,C3') coupling could be structurally assigned to torsion ζ and the (2)J(P,C5') coupling to torsion α. We propose qualitative rules for their structural interpretation as loose restraints on torsion angles ζ and α. The (3)J(P,C4') coupling assigned to torsion angle β was found dependent also on torsions ζ and α, implying that the uncertainty in determination of β with standard Karplus curves could be as large as ~25°. The calculations provided a unified picture of NMR parameters applicable for the determination of NA phosphate conformation.

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.

Revision of AMBER Torsional Parameters for RNA Improves Free Energy Predictions for Tetramer Duplexes with GC and iGiC Base Pairs

All-atom force fields are important for predicting thermodynamic, structural, and dynamic properties of RNA. In this paper, results are reported for thermodynamic integration calculations of free energy differences of duplex formation when CG pairs in the RNA duplexes r(CCGG)(2), r(GGCC)(2), r(GCGC)(2), and r(CGCG)(2) are replaced by isocytidine-isoguanosine (iCiG) pairs. Agreement with experiment was improved when ε/ζ, α/γ, β, and χ torsional parameters in the AMBER99 force field were revised on the basis of quantum mechanical calculations. The revised force field, AMBER99TOR, brings free energy difference predictions to within 1.3, 1.4, 2.3, and 2.6 kcal/mol at 300 K, respectively, compared to experimental results for the thermodynamic cycles of CCGG → iCiCiGiG, GGCC → iGiGiCiC, GCGC → iGiCiGiC, and CGCG → iCiGiCiG. In contrast, unmodified AMBER99 predictions for GGCC → iGiGiCiC and GCGC → iGiCiGiC differ from experiment by 11.7 and 12.6 kcal/mol, respectively. In order to test the dynamic stability of the above duplexes with AMBER99TOR, four individual 50 ns molecular dynamics (MD) simulations in explicit solvent were run. All except r(CCGG)(2) retained A-form conformation for ≥82% of the time. This is consistent with NMR spectra of r(iGiGiCiC)(2), which reveal an A-form conformation. In MD simulations, r(CCGG)(2) retained A-form conformation 52% of the time, suggesting that its terminal base pairs may fray. The results indicate that revised backbone parameters improve predictions of RNA properties and that comparisons to measured sequence dependent thermodynamics provide useful benchmarks for testing force fields and computational methods.

Structure of an RNA/DNA dodecamer corresponding to the HIV-1 polypurine tract at 1.6 Å resolution

The crystal structure of an RNA/DNA hybrid dodecamer, r(5'-uaaaagaaaagg):d(5'-CCTTTTCTTTTA), which contains three-quarters of the polypurine tract (PPT) sequence of the HIV RNA genome is reported. The hybrid structure was determined at 1.6 Å resolution and was found to have the A-form conformation. However, the presence of alternate conformations along the RNA template strand indicated increased flexibility of the PPT sequence. Two segments (at nucleotides 1-2 and 6-8) of the RNA chain have two conformations exhibiting differences in torsion and pseudorotation angles. For conformation I((1-2), (6-8)), 25% of the RNA sugars have the C2'-exo pucker and the rest have the expected C3'-endo pucker. The II(1-2) and II(6-8) conformations of the RNA strand have one sugar with the C2'-exo pucker. None of the ribose rings exist in the C2'-endo form, in contrast to a previous report which postulated a C2'-endo ribose as a key structural element of the PPT. The widths of the minor groove for conformations I((1-2), (6-8)) and II((1-2), (6-8)) of the RNA strand are 9.2-10.5 and 9.4-10.7 Å, respectively. Both ranges are very close to the intervals accepted for A-form RNA duplexes. On the opposing DNA primer strand most of the sugars are C3'-endo, except for the 3'-terminal sugars, which are C2'-endo (T22) or O4'-endo (T23 and A24). The duplex includes a noncanonical u1(anti)·A24(syn) base interaction with only one hydrogen bond between the bases. This noncanonical base interaction at the 5'-end of the template distorts the values of the helical parameters of the adjacent base pair.

Nonparametric clustering for studying RNA conformations

The local conformation of RNA molecules is an important factor in determining their catalytic and binding properties. The analysis of such conformations is particularly difficult due to the large number of degrees of freedom, such as the measured torsion angles per residue and the interatomic distances among interacting residues. In this work, we use a nearest-neighbor search method based on the statistical mechanical Potts model to find clusters in the RNA conformational space. The proposed technique is mostly automatic and may be applied to problems, where there is no prior knowledge on the structure of the data space in contrast to many other clustering techniques. Results are reported for both single residue conformations, where the parameter set of the data space includes four to seven torsional angles, and base pair geometries, where the data space is reduced to two dimensions. Moreover, new results are reported for base stacking geometries. For the first two cases, i.e., single residue conformations and base pair geometries, we get a very good match between the results of the proposed clustering method and the known classifications with only few exceptions. For the case of base stacking geometries, we validate our classification with respect to geometrical constraints and describe the content, and the geometry of the new clusters.

Benchmarking AMBER force fields for RNA: comparisons to NMR spectra for single-stranded r(GACC) are improved by revised χ torsions

Accurately modeling unpaired regions of RNA is important for predicting structure, dynamics, and thermodynamics of folded RNA. Comparisons between NMR data and molecular dynamics simulations provide a test of force fields used for modeling. Here, NMR spectroscopy, including NOESY, (1)H-(31)P HETCOR, DQF-COSY, and TOCSY, was used to determine conformational preferences for single-stranded GACC RNA. The spectra are consistent with a conformational ensemble containing major and minor A-form-like structures. In a series of 50 ns molecular dynamics (MD) simulations with the AMBER99 force field in explicit solvent, initial A-form-like structures rapidly evolve to disordered conformations. A set of 50 ns simulations with revised χ torsions (AMBER99χ force field) gives two primary conformations, consistent with the NMR spectra. A single 1.9 μs MD simulation with the AMBER99χ force field showed that the major and minor conformations are retained for almost 68% of the time in the first 700 ns, with multiple transformations from A-form to non-A-form conformations. For the rest of the simulation, random-coil structures and a stable non-A-form conformation inconsistent with NMR spectra were seen. Evidently, the AMBER99χ force field improves structural predictions for single-stranded GACC RNA compared to the AMBER99 force field, but further force field improvements are needed.

A folding "framework structure" of Tetrahymena group I intron

We have published the dynamic extended folding (DEF) method, which is a RNA secondary structure prediction approach-to simulate the in vivo RNA co-transcriptional folding process. In order to verify the reliability of the method, we selected the X-ray-determined Tetrahymena group I intron as a sample to construct the framework of its folding secondary structure. Our prediction coincides well with the secondary structure predicted by T.R. Cech and the X-ray diffraction crystal structure determined by Lehnert V. Our results show that the DEF framework structure of Tetrahymena group I intron reflects its function sites in a concise and straightforward manner, and the scope of the simulation was expanded.

Quantitative 2D and 3D Gamma-HCP experiments for the determination of the angles alpha and zeta in the phosphodiester backbone of oligonucleotides

The quantitative Gamma-(HCP) experiment, a novel heteronuclear NMR pulse sequence for the determination of the RNA backbone angles alpha(O3'(i-1)-P(i)-O5'(i)-C5'(i)) and zeta(C3'(i)-O3'(i)-P(i+1)-O5'(i+1)) in (13)C-labeled RNA, is introduced. The experiment relies on the interaction between the CH bond vector dipole and the (31)P chemical shift anisotropy (CSA), which affects the relaxation of the (13)C,(31)P double- and zero-quantum coherence and thus the intensity of the detectable magnetization. With the new pulse sequence, five different cross-correlated relaxation rates along the phosphodiester backbone can be measured in a quantitative manner, allowing projection-angle and torsion-angle restraints for the two backbone angles alpha and zeta to be extracted. Two versions of the pulse sequence optimized for the CH and CH(2) groups are introduced and demonstrated for a 14-mer cUUCGg tetraloop RNA model system and for a 27-mer RNA with a previously unknown structure. The restraints were incorporated into the calculation of a very high resolution structure of the RNA model system (Nozinovic, S.; et al. Nucleic Acids Res. 2010, 38, 683). Comparison with the X-ray structure of the cUUCGg tetraloop confirmed the high quality of the data, suggesting that the method can significantly improve the quality of RNA structure determination.

Geometric properties of nucleic acids with potential for autobuilding

Algorithms and geometrical properties are described for the automated building of nucleic acids in experimental electron density.

Medium- to high-resolution X-ray structures of DNA and RNA molecules were investigated to find geometric properties useful for automated model building in crystallographic electron-density maps. We describe a simple method, starting from a list of electron-density ‘blobs’, for identifying backbone phosphates and nucleic acid bases based on properties of the local electron-density distribution. This knowledge should be useful for the automated building of nucleic acid models into electron-density maps. We show that the distances and angles involving C1′ and the P atoms, using the pseudo-torsion angles and that describe the …P—C1′—P—C1′… chain, provide a promising basis for building the nucleic acid polymer. These quantities show reasonably narrow distributions with asymmetry that should allow the direction of the phosphate backbone to be established.

Reparameterization of RNA chi Torsion Parameters for the AMBER Force Field and Comparison to NMR Spectra for Cytidine and Uridine

A reparameterization of the torsional parameters for the glycosidic dihedral angle, chi, for the AMBER99 force field in RNA nucleosides is used to provide a modified force field, AMBER99chi. Molecular dynamics simulations of cytidine, uridine, adenosine, and guanosine in aqueous solution using the AMBER99 and AMBER99chi force fields are compared with NMR results. For each nucleoside and force field, 10 individual molecular dynamics simulations of 30 ns each were run. For cytidine with AMBER99chi force field, each molecular dynamics simulation time was extended to 120 ns for convergence purposes. Nuclear magnetic resonance (NMR) spectroscopy, including one-dimensional (1D) (1)H, steady-state 1D (1)H nuclear Overhauser effect (NOE), and transient 1D (1)H NOE, was used to determine the sugar puckering and preferred base orientation with respect to the ribose of cytidine and uridine. The AMBER99 force field overestimates the population of syn conformations of the base orientation and of C2'-endo sugar puckering of the pyrimidines, while the AMBER99chi force field's predictions are more consistent with NMR results. Moreover, the AMBER99 force field prefers high anti conformations with glycosidic dihedral angles around 310 degrees for the base orientation of purines. The AMBER99chi force field prefers anti conformations around 185 degrees , which is more consistent with the quantum mechanical calculations and known 3D structures of folded ribonucleic acids (RNAs). Evidently, the AMBER99chi force field predicts the structural characteristics of ribonucleosides better than the AMBER99 force field and should improve structural and thermodynamic predictions of RNA structures.

Identification of dynamical hinge points of the L1 ligase molecular switch

The L1 ligase is an in vitro selected ribozyme that uses a noncanonically base-paired ligation site to catalyze regioselectively and regiospecifically the 5' to 3' phosphodiester bond ligation, a reaction relevant to origin of life hypotheses that invoke an RNA world scenario. The L1 ligase crystal structure revealed two different conformational states that were proposed to represent the active and inactive forms. It remains an open question as to what degree these two conformers persist as stable conformational intermediates in solution, and along what pathway are they able to interconvert. To explore these questions, we have performed a series of molecular dynamics simulations in explicit solvent of the inactive-active conformational switch in L1 ligase. Four simulations were performed departing from both conformers in both the reactant and product states, in addition to a simulation where local unfolding in the active state was induced. From these simulations, along with crystallographic data, a set of four virtual torsion angles that span two evolutionarily conserved and restricted regions were identified as dynamical hinge points in the conformational switch transition. The ligation site visits three distinct states characterized by hydrogen bond patterns that are correlated with the formation of specific contacts that may promote catalysis. The insights gained from these simulations contribute to a more detailed understanding of the coupled catalytic/conformational switch mechanism of L1 ligase that may facilitate the design and engineering of new catalytic riboswitches.

Comparison of intrinsic stacking energies of ten unique dinucleotide steps in A-RNA and B-DNA duplexes. Can we determine correct order of stability by quantum-chemical calculations?

High level ab initio methods have been used to study stacking interactions in ten unique base pair steps both in A-RNA and in B-DNA duplexes. The protocol for selection of geometries based on molecular dynamics (MD) simulations is proposed, and its suitability is demonstrated by comparison with stacking in steps at fiber diffraction geometries. It is shown that fiber diffraction geometries are not sufficiently accurate for interaction energy calculations. In addition, the protocol for selection of geometries based on MD simulations allows for the evaluation of the variability of the intrinsic stacking energies along the MD trajectories. The uncertainty in stacking energies (difference between the most and least stable geometry) due to the dynamical nature of systems can be, in some cases, as large as 3.0 kcal x mol(-1), which is almost 50% of the actual sequence dependence of base stacking energies (the energy difference between the most and least stable sequences). Thus, assessing the relative magnitude of the gas phase stacking energy using a single geometry for each sequence is insufficient to obtain an unambiguous order of gas phase stacking energies in canonical double helices. Though the ordering of ten unique dinucleotide steps cannot be definitive, some general conclusions were drawn. The stacking energies of base pair steps in A-RNA are more evenly separated compared to B-DNA, and their ordering is less sensitive to the dynamics of the system compared to be B-DNA. The most stable step both in B-DNA and A-RNA is the GC/GC [corrected] step that is well separated from the second most stable step CG/CG. [corrected] Also the least stable step (the CC/GG step) is well separated from the rest of the structures. The calculations further show that B-DNA stacking is favorable only marginally (on average by 1.14 kcal x mol(-1) per base pair step) over A-RNA stacking, and this difference vanishes after subtracting the stabilizing van der Waals effect of the thymine 5-methyl group that is absent in RNA. Basically, no correlation between the sequence dependence of gas phase stacking energies and the sequence dependence of DeltaG degrees(37) free energies used in nearest-neighbor models was found either for B-DNA or for A-RNA. This reflects the complexity of the balance of forces that are responsible for the sequence dependence of thermodynamics stability of nucleic acids, which masks the effect of the intrinsic interactions between the stacked base pairs.

RNA structural motifs that entail hydrogen bonds involving sugar-phosphate backbone atoms of RNA

The growing number of high-resolution crystal structures of large RNA molecules provides much information for understanding the principles of structural organization of these complex molecules. Several in-depth analyses of nucleobase-centered RNA structural motifs and backbone conformations have been published based on this information, including a systematic classification of base pairs by Leontis and Westhof. However, hydrogen bonds involving sugar-phosphate backbone atoms of RNA have not been analyzed systematically until recently, although such hydrogen bonds appear to be common both in local and tertiary interactions. Here we review some backbone structural motifs discussed in the literature and analyze a set of eight high-resolution multi-domain RNA structures. The analyzed RNAs are highly structured: among 5372 nucleotides in this set, 89% are involved in at least one "long-range" RNA-RNA hydrogen bond, i.e., hydrogen bonds between atoms in the same residue or sequential residues are ignored. These long-range hydrogen bonds frequently use backbone atoms as hydrogen bond acceptors, i.e., OP1, OP2, O2', O3', O4', or O5', or as a donor (2'OH). A surprisingly large number of such hydrogen bonds are found, considering that neither single-stranded nor double-stranded regions will contain such hydrogen bonds unless additional interactions with other residues exist. Among 8327 long-range hydrogen bonds found in this set of structures, 2811, or about one-third, are hydrogen bonds entailing RNA backbone atoms; they involve 39% of all nucleotides in the structures. The majority of them (2111) are hydrogen bonds entailing ribose hydroxyl groups, which can be used either as a donor or an acceptor; they constitute 25% of all hydrogen bonds and involve 31% of all nucleotides. The phosphate oxygens OP1 or OP2 are used as hydrogen bond acceptors in 12% of all nucleotides, and the ribose ring oxygen O4' and phosphodiester oxygens O3' and O5' are used in 4%, 4%, and 1% of all nucleotides, respectively. Distributions of geometric parameters and some examples of such hydrogen bonds are presented in this report. A novel motif involving backbone hydrogen bonds, the ribose-phosphate zipper, is also identified.

An RNA molecular switch: Intrinsic flexibility of 23S rRNA Helices 40 and 68 5'-UAA/5'-GAN internal loops studied by molecular dynamics methods

Functional RNA molecules such as ribosomal RNAs frequently contain highly conserved internal loops with a 5'-UAA/5'-GAN (UAA/GAN) consensus sequence. The UAA/GAN internal loops adopt distinctive structure inconsistent with secondary structure predictions. The structure has a narrow major groove and forms a trans Hoogsteen/Sugar edge (tHS) A/G base pair followed by an unpaired stacked adenine, a trans Watson-Crick/Hoogsteen (tWH) U/A base pair and finally by a bulged nucleotide (N). The structure is further stabilized by a three-adenine stack and base-phosphate interaction. In the ribosome, the UAA/GAN internal loops are involved in extensive tertiary contacts, mainly as donors of A-minor interactions. Further, this sequence can adopt an alternative 2D/3D pattern stabilized by a four-adenine stack involved in a smaller number of tertiary interactions. The solution structure of an isolated UAA/GAA internal loop shows substantially rearranged base pairing with three consecutive non-Watson-Crick base pairs. Its A/U base pair adopts an incomplete cis Watson-Crick/Sugar edge (cWS) A/U conformation instead of the expected Watson-Crick arrangement. We performed 3.1 µs of explicit solvent molecular dynamics (MD) simulations of the X-ray and NMR UAA/GAN structures, supplemented by MM-PBSA free energy calculations, locally enhanced sampling (LES) runs, targeted MD (TMD) and nudged elastic band (NEB) analysis. We compared parm99 and parmbsc0 force fields and net-neutralizing Na(+) vs. excess salt KCl ion environments. Both force fields provide a similar description of the simulated structures, with the parmbsc0 leading to modest narrowing of the major groove. The excess salt simulations also cause a similar effect. While the NMR structure is entirely stable in simulations, the simulated X-ray structure shows considerable widening of the major groove, loss of base-phosphate interaction and other instabilities. The alternative X-ray geometry even undergoes conformational transition towards the solution 2D structure. Free energy calculations confirm that the X-ray arrangement is less stable than the solution structure. LES, TMD and NEB provide a rather consistent pathway for interconversion between the X-ray and NMR structures. In simulations, the incomplete cWS A/U base pair of the NMR structure is water mediated and alternates with the canonical A-U base pair, which is not indicated by the NMR data. Completion of full cWS A/U base pair is prevented by the overall internal loop arrangement. In summary, the simulations confirm that the UAA/GAN internal loop is a molecular switch RNA module that adopts its functional geometry upon specific tertiary contexts.

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

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

MD simulations of ligand-bound and ligand-free aptamer: molecular level insights into the binding and switching mechanism of the add A-riboswitch

Riboswitches are structural cis-acting genetic regulatory elements in 5' UTRs of mRNAs, consisting of an aptamer domain that regulates the behavior of an expression platform in response to its recognition of, and binding to, specific ligands. While our understanding of the ligand-bound structure of the aptamer domain of the adenine riboswitches is based on crystal structure data and is well characterized, understanding of the structure and dynamics of the ligand-free aptamer is limited to indirect inferences from physicochemical probing experiments. Here we report the results of 15-nsec-long explicit-solvent molecular dynamics simulations of the add A-riboswitch crystal structure (1Y26), both in the adenine-bound (CLOSED) state and in the adenine-free (OPEN) state. Root-mean-square deviation, root-mean-square fluctuation, dynamic cross-correlation, and backbone torsion angle analyses are carried out on the two trajectories. These, along with solvent accessible surface area analysis of the two average structures, are benchmarked against available experimental data and are shown to constitute the basis for obtaining reliable insights into the molecular level details of the binding and switching mechanism. Our analysis reveals the interaction network responsible for, and conformational changes associated with, the communication between the binding pocket and the expression platform. It further highlights the significance of a, hitherto unreported, noncanonical W:H trans base pairing between A73 and A24, in the OPEN state, and also helps us to propose a possibly crucial role of U51 in the context of ligand binding and ligand discrimination.

Conformational evaluation of labeled C3'-O-P-(13)CH(2)-O-C4'' phosphonate internucleotide linkage, a phosphodiester isostere

Modified internucleotide linkage featuring the C3'-O-P-CH(2)-O-C4'' phosphonate grouping as an isosteric alternative to the phosphodiester C3'-O-P-O-CH(2)-C4'' bond was studied in order to learn more on its stereochemical arrangement, which we showed earlier to be of prime importance for the properties of the respective oligonucleotide analogues. Two approaches were pursued: First, the attempt to prepare the model dinucleoside phosphonate with (13)C-labeled CH(2) group present in the modified internucleotide linkage that would allow for a more detailed evaluation of the linkage conformation by NMR spectroscopy. Second, the use of ab initio calculations along with molecular dynamics (MD) simulations in order to observe the most populated conformations and specify main structural elements governing the conformational preferences. To deal with the former aim, a novel synthesis of key labeled reagent (CH(3)O)(2)P(O)(13)CH(2)OH for dimer preparation had to be elaborated using aqueous (13)C-formaldehyde. The results from both approaches were compared and found consistent.

A probabilistic model of RNA conformational space

The increasing importance of non-coding RNA in biology and medicine has led to a growing interest in the problem of RNA 3-D structure prediction. As is the case for proteins, RNA 3-D structure prediction methods require two key ingredients: an accurate energy function and a conformational sampling procedure. Both are only partly solved problems. Here, we focus on the problem of conformational sampling. The current state of the art solution is based on fragment assembly methods, which construct plausible conformations by stringing together short fragments obtained from experimental structures. However, the discrete nature of the fragments necessitates the use of carefully tuned, unphysical energy functions, and their non-probabilistic nature impairs unbiased sampling. We offer a solution to the sampling problem that removes these important limitations: a probabilistic model of RNA structure that allows efficient sampling of RNA conformations in continuous space, and with associated probabilities. We show that the model captures several key features of RNA structure, such as its rotameric nature and the distribution of the helix lengths. Furthermore, the model readily generates native-like 3-D conformations for 9 out of 10 test structures, solely using coarse-grained base-pairing information. In conclusion, the method provides a theoretical and practical solution for a major bottleneck on the way to routine prediction and simulation of RNA structure and dynamics in atomic detail.

Structure and dynamics of the ApA, ApC, CpA, and CpC RNA dinucleoside monophosphates resolved with NMR scalar spin-spin couplings

The measured NMR scalar coupling constants (J-couplings) in the XpY, (X,Y = adenine (A) or cytosine (C)) RNA dinucleoside monophosphates (DMPs) were assigned to the backbone (alpha, beta, gamma, delta, epsilon, zeta) and glycosidic (chi) torsion angles in order to resolve the global structure of the DMP molecules. The experimental J-couplings were correlated with the theoretical J-couplings obtained as the dynamical averages of the Karplus equations relevant to the torsion angles. The dynamical information was captured using the molecular dynamics (MD) calculation method. The individual conformational flexibility of the four DMP molecules was thus consistently probed with the NMR J-couplings. The calculated structure and flexibility of the DMP molecules depend on the sequence considered with respect to the 5' and 3' end of the DMP molecules (5'-XpY-3'). The dynamical characteristics of the two nucleosides are not equivalent even for the ApA and CpC homologues. An enhancement of the sampling in the MD calculations was achieved using five different starting structural motives classified previously for the RNA backbone in the solid phase (Richardson et al. RNA 2008, 14, 465-481). The initial structures were selected on the basis of a database search for RNA oligonucleotides. Frequent interconversions between the conformers during the MD calculations were actually observed. The structural interpretation of the NMR spectroscopic data based on the MD simulations combined with the Karplus equations indicates that the dominant conformation of the DMP molecules in solution corresponds to the A-RNA form. For 52% of the total simulation time (1000 ns), the zeta(g-)-alpha(g-)-gamma(g+) backbone topology corresponding to the canonical A-RNA form was observed, with roughly equally populated C2'- and C3'-endo sugar puckers interconverting on the nanosecond time scale. However, other noncanonical patterns were also found and thus indicate their relatively high potential to be populated in the dynamical regime. For approximately 72% of the time portion when the A-RNA of the zeta-alpha-gamma combination occurred, the nucleobases were classified as being mutually stacked. The geometries of the nucleobases classified in this work as stacked were significantly more populated for the DMP molecules with adenosine at the 3' end (ApA and CpA DMPs) than the ApC or CpC RNA molecules with C at the 3' end.

DNA conformations and their sequence preferences

The geometry of the phosphodiester backbone was analyzed for 7739 dinucleotides from 447 selected crystal structures of naked and complexed DNA. Ten torsion angles of a near-dinucleotide unit have been studied by combining Fourier averaging and clustering. Besides the known variants of the A-, B- and Z-DNA forms, we have also identified combined A + B backbone-deformed conformers, e.g. with alpha/gamma switches, and a few conformers with a syn orientation of bases occurring e.g. in G-quadruplex structures. A plethora of A- and B-like conformers show a close relationship between the A- and B-form double helices. A comparison of the populations of the conformers occurring in naked and complexed DNA has revealed a significant broadening of the DNA conformational space in the complexes, but the conformers still remain within the limits defined by the A- and B- forms. Possible sequence preferences, important for sequence-dependent recognition, have been assessed for the main A and B conformers by means of statistical goodness-of-fit tests. The structural properties of the backbone in quadruplexes, junctions and histone-core particles are discussed in further detail.

Molecular modeling of nucleic acid structure

This unit is the first in a series of four units covering the analysis of nucleic acid structure by molecular modeling. This unit provides an overview of computer simulation of nucleic acids. Topics include the static structure model, computational graphics and energy models, generation of an initial model, and characterization of the overall three-dimensional structure.

RNA backbone: consensus all-angle conformers and modular string nomenclature (an RNA Ontology Consortium contribution)

A consensus classification and nomenclature are defined for RNA backbone structure using all of the backbone torsion angles. By a consensus of several independent analysis methods, 46 discrete conformers are identified as suitably clustered in a quality-filtered, multidimensional dihedral angle distribution. Most of these conformers represent identifiable features or roles within RNA structures. The conformers are given two-character names that reflect the seven-angle delta epsilon zeta alpha beta gamma delta combinations empirically found favorable for the sugar-to-sugar "suite" unit within which the angle correlations are strongest (e.g., 1a for A-form, 5z for the start of S-motifs). Since the half-nucleotides are specified by a number for delta epsilon zeta and a lowercase letter for alpha beta gamma delta, this modular system can also be parsed to describe traditional nucleotide units (e.g., a1) or the dinucleotides (e.g., a1a1) that are especially useful at the level of crystallographic map fitting. This nomenclature can also be written as a string with two-character suite names between the uppercase letters of the base sequence (N1aG1gN1aR1aA1cN1a for a GNRA tetraloop), facilitating bioinformatic comparisons. Cluster means, standard deviations, coordinates, and examples are made available, as well as the Suitename software that assigns suite conformer names and conformer match quality (suiteness) from atomic coordinates. The RNA Ontology Consortium will combine this new backbone system with others that define base pairs, base-stacking, and hydrogen-bond relationships to provide a full description of RNA structural motifs.

Conformations of flanking bases in HIV-1 RNA DIS kissing complexes studied by molecular dynamics

Explicit solvent molecular dynamics simulations (in total almost 800 ns including locally enhanced sampling runs) were applied with different ion conditions and with two force fields (AMBER and CHARMM) to characterize typical geometries adopted by the flanking bases in the RNA kissing-loop complexes. We focus on flanking base positions in multiple x-ray and NMR structures of HIV-1 DIS kissing complexes and kissing complex from the large ribosomal subunit of Haloarcula marismortui. An initial x-ray open conformation of bulged-out bases in HIV-1 DIS complexes, affected by crystal packing, tends to convert to a closed conformation formed by consecutive stretch of four stacked purine bases. This is in agreement with those recent crystals where the packing is essentially avoided. We also observed variants of the closed conformation with three stacked bases, while nonnegligible populations of stacked geometries with bulged-in bases were detected, too. The simulation results reconcile differences in positions of the flanking bases observed in x-ray and NMR studies. Our results suggest that bulged-out geometries are somewhat more preferred, which is in accord with recent experiments showing that they may mediate tertiary contacts in biomolecular assemblies or allow binding of aminoglycoside antibiotics.

FR3D: finding local and composite recurrent structural motifs in RNA 3D structures

New methods are described for finding recurrent three-dimensional (3D) motifs in RNA atomic-resolution structures. Recurrent RNA 3D motifs are sets of RNA nucleotides with similar spatial arrangements. They can be local or composite. Local motifs comprise nucleotides that occur in the same hairpin or internal loop. Composite motifs comprise nucleotides belonging to three or more different RNA strand segments or molecules. We use a base-centered approach to construct efficient, yet exhaustive search procedures using geometric, symbolic, or mixed representations of RNA structure that we implement in a suite of MATLAB programs, "Find RNA 3D" (FR3D). The first modules of FR3D preprocess structure files to classify base-pair and -stacking interactions. Each base is represented geometrically by the position of its glycosidic nitrogen in 3D space and by the rotation matrix that describes its orientation with respect to a common frame. Base-pairing and base-stacking interactions are calculated from the base geometries and are represented symbolically according to the Leontis/Westhof basepairing classification, extended to include base-stacking. These data are stored and used to organize motif searches. For geometric searches, the user supplies the 3D structure of a query motif which FR3D uses to find and score geometrically similar candidate motifs, without regard to the sequential position of their nucleotides in the RNA chain or the identity of their bases. To score and rank candidate motifs, FR3D calculates a geometric discrepancy by rigidly rotating candidates to align optimally with the query motif and then comparing the relative orientations of the corresponding bases in the query and candidate motifs. Given the growing size of the RNA structure database, it is impossible to explicitly compute the discrepancy for all conceivable candidate motifs, even for motifs with less than ten nucleotides. The screening algorithm that we describe finds all candidate motifs whose geometric discrepancy with respect to the query motif falls below a user-specified cutoff discrepancy. This technique can be applied to RMSD searches. Candidate motifs identified geometrically may be further screened symbolically to identify those that contain particular basepair types or base-stacking arrangements or that conform to sequence continuity or nucleotide identity constraints. Purely symbolic searches for motifs containing user-defined sequence, continuity and interaction constraints have also been implemented. We demonstrate that FR3D finds all occurrences, both local and composite and with nucleotide substitutions, of sarcin/ricin and kink-turn motifs in the 23S and 5S ribosomal RNA 3D structures of the H. marismortui 50S ribosomal subunit and assigns the lowest discrepancy scores to bona fide examples of these motifs. The search algorithms have been optimized for speed to allow users to search the non-redundant RNA 3D structure database on a personal computer in a matter of minutes.

Indirect NMR spin-spin coupling constants 3J(P,C) and 2J(P,H) across the P-O...H-C link can be used for structure determination of nucleic acids

Calculated indirect NMR spin-spin coupling constants (3)J(P,C) and (2)J(P,H) were correlated with the local structure of the P-O...H-C linkage between the nucleic acid (NA) backbone phosphate and the H-C group(s) of a nucleic acid base. The calculations were carried out for selected nucleotides from the large ribosomal subunit (Ban et al. Science 2000, 289, 905) with the aim of identifying NMR parameters suitable for detection of certain noncanonical RNA structures. As calculations in the model system, dimethyl-phosphate-guanine, suggest, the calculated indirect spin-spin couplings across the linkage are sensitive to the mutual orientation and distance between the phosphate and nucleic acid base. A short distance between the nucleic acid base and phosphate group and the angles C...P-O and P...C-H smaller than 50 degrees are prerequisites for a measurable spin-spin interaction of either coupling (|J| > 1 Hz). A less favorable arrangement of the P-O...H-C motif, e.g., in nucleotides of the canonical A-RNA, results in an effective dumping of both spin-spin interactions and insignificant values of the NMR coupling constants. The present work indicates that quantum chemical calculations of the indirect spin-spin couplings across the P-O...H-C motif can help detect some rare but important backbone topologies, as seen for example in the reverse kink-turn. Measuring of (3)J(P,C) and (2)J(P,H) couplings can therefore provide critical constraints on the NA base and phosphate geometry and help to determine the structure of NAs.

Evaluating and learning from RNA pseudotorsional space: quantitative validation of a reduced representation for RNA structure

Quantitatively describing RNA structure and conformational elements remains a formidable problem. Seven standard torsion angles and the sugar pucker are necessary to characterize the conformation of an RNA nucleotide completely. Progress has been made toward understanding the discrete nature of RNA structure, but classifying simple and ubiquitous structural elements such as helices and motifs remains a difficult task. One approach for describing RNA structure in a simple, mathematically consistent, and computationally accessible manner involves the invocation of two pseudotorsions, eta (C4'(n-1), P(n), C4'(n), P(n+1)) and theta (P(n), C4'(n), P(n+1), C4'(n+1)), which can be used to describe RNA conformation in much the same way that varphi and psi are used to describe backbone configuration of proteins. Here, we conduct an exploration and statistical evaluation of pseudotorsional space and of the Ramachandran-like eta-theta plot. We show that, through the rigorous quantitative analysis of the eta-theta plot, the pseudotorsional descriptors eta and theta, together with sugar pucker, are sufficient to describe RNA backbone conformation fully in most cases. These descriptors are also shown to contain considerable information about nucleotide base conformation, revealing a previously uncharacterized interplay between backbone and base orientation. A window function analysis is used to discern statistically relevant regions of density in the eta-theta scatter plot and then nucleotides in colocalized clusters in the eta-theta plane are shown to have similar 3-D structures through RMSD analysis of the RNA structural constituents. We find that major clusters in the eta-theta plot are few, underscoring the discrete nature of RNA backbone conformation. Like the Ramachandran plot, the eta-theta plot is a valuable system for conceptualizing biomolecular conformation, it is a useful tool for analyzing RNA tertiary structures, and it is a vital component of new approaches for solving the 3-D structures of large RNA molecules and RNA assemblies.

Principles of RNA base pairing: structures and energies of the trans Watson-Crick/sugar edge base pairs

Due to the presence of the 2'-OH hydroxyl group of ribose, RNA molecules utilize an astonishing variability of base pairing patterns to build up their structures and perform the biological functions. Many of the key RNA base pairing families have no counterparts in DNA. In this study, the trans Watson-Crick/sugar edge (trans WC/SE) RNA base pair family has been characterized using quantum chemical and molecular mechanics calculations. Gas-phase optimized geometries from density functional theory (DFT) calculations and RIMP2 interaction energies are reported for the 10 crystallographically identified trans WC/SE base pairing patterns. Further, stable structures are predicted for all of the remaining six possible members of this family not seen in RNAs so far. Among these novel six base pairs, the computations substantially refine two structures suggested earlier based on simple isosteric considerations. For two additional trans WC/SE base pairs predicted in this study, no arrangement was suggested before. Thus, our study brings a complete set of trans WC/SE base pairing patterns. The present results are also contrasted with calculations reported recently for the cis WC/SE base pair family. The computed base pair sizes are in sound correlation with the X-ray data for all WC/SE pairing patterns including both their cis and trans isomers. This confirms that the isostericity of RNA base pairs, which is one of the key factors determining the RNA sequence conservation patterns, originates in the properties of the isolated base pairs. In contrast to the cis structures, however, the isosteric subgroups of the trans WC/SE family differ not only in their H-bonding patterns and steric dimensions but also in the intrinsic strength of the intermolecular interactions. The distribution of the total interaction energy over the sugar-base and base-base contributions is controlled by the cis-trans isomerism.

AUTOMATIC STRUCTURE INTERPRETATION OF SINGLE PARTICLE CRYO-ELECTRON MICROSCOPY: FROM IMAGES TO PSUEDO-ATOMIC MODELS

Three dimensional Electron Microscopy (EM) and in particular single particle reconstruction using cryo-EM, has rapidly advanced over recent years, such that increasingly several macromolecular complexes can be resolved at subnanometer resolution (6-10 Å). This paper reviews some of the main volumetric image and geometric post-processing steps once a three dimensional EM map (henceforth a 3D map) has been reconstructed from single particle Cryo-EM, as essential steps in an enhanced and automated computational structure interpretation pipeline. In particular the paper addresses automated filtering, critical point calculations, symmetric and non-symmetric molecular domain segmentation, molecular surface selection, curation, and protein secondary structure (α- helices and β-sheets) elucidation from 3D maps.

RNABC: forward kinematics to reduce all-atom steric clashes in RNA backbone

Although accurate details in RNA structure are of great importance for understanding RNA function, the backbone conformation is difficult to determine, and most existing RNA structures show serious steric clashes (>or= 0.4 A overlap) when hydrogen atoms are taken into account. We have developed a program called RNABC (RNA Backbone Correction) that performs local perturbations to search for alternative conformations that avoid those steric clashes or other local geometry problems. Its input is an all-atom coordinate file for an RNA crystal structure (usually from the MolProbity web service), with problem areas specified. RNABC rebuilds a suite (the unit from sugar to sugar) by anchoring the phosphorus and base positions, which are clearest in crystallographic electron density, and reconstructing the other atoms using forward kinematics. Geometric parameters are constrained within user-specified tolerance of canonical or original values, and torsion angles are constrained to ranges defined through empirical database analyses. Several optimizations reduce the time required to search the many possible conformations. The output results are clustered and presented to the user, who can choose whether to accept one of the alternative conformations. Two test evaluations show the effectiveness of RNABC, first on the S-motifs from 42 RNA structures, and second on the worst problem suites (clusters of bad clashes, or serious sugar pucker outliers) in 25 unrelated RNA structures. Among the 101 S-motifs, 88 had diagnosed problems, and RNABC produced clash-free conformations with acceptable geometry for 71 of those (about 80%). For the 154 worst problem suites, RNABC proposed alternative conformations for 72. All but 8 of those were judged acceptable after examining electron density (where available) and local conformation. Thus, even for these worst cases, nearly half the time RNABC suggested corrections suitable to initiate further crystallographic refinement. The program is available from http://kinemage.biochem.duke.edu .

Relationships between 31P chemical shift tensors and conformation of nucleic acid backbone: a DFT study

Density functional theory (DFT) has been applied to study the conformational dependence of 31P chemical shift tensors in B-DNA. The gg and gt conformations of backbone phosphate groups representing BI- and BII-DNA have been examined. Calculations have been carried out on static models of dimethyl phosphate (dmp) and dinucleoside-3',5'-monophosphate with bases replaced by hydrogen atoms in vacuo as well as in an explicit solvent. Trends in 31P chemical shift anisotropy (CSA) tensors with respect to the backbone torsion angles alpha, zeta, beta, and epsilon are presented. Although these trends do not change qualitatively upon solvation, quantitative changes result in the reduction of the chemical shift anisotropy. For alpha and zeta in the range from 270 degrees to 330 degrees and from 240 degrees to 300 degrees , respectively, the delta22 and delta33 principal components vary within as much as 30 ppm, showing a marked dependence on backbone conformation. The calculated 31P chemical shift tensor principal axes deviate from the axes of O-P-O bond angles by at most 5 degrees . For solvent models, our results are in a good agreement with experimental estimates of relative gg and gt isotropic chemical shifts. Solvation also brings the theoretical deltaiso of the gg conformation closer to the experimental gg data of barium diethyl phosphate.

Calculation of structural behavior of indirect NMR spin-spin couplings in the backbone of nucleic acids

Calculated indirect NMR spin-spin coupling constants (J-couplings) between (31)P, (13)C, and (1)H nuclei were related to the backbone torsion angles of nucleic acids (NAs), and it was shown that J-couplings can facilitate accurate and reliable structural interpretation of NMR measurements and help to discriminate between their distinct conformational classes. A proposed stepwise procedure suggests assignment of the J-couplings to torsion angles from the sugar part to the phosphodiester link. Some J-couplings show multidimensional dependence on torsion angles, the most prominent of which is the effect of the sugar pucker. J-couplings were calculated in 16 distinct nucleic acid conformations, two principal double-helical DNAs, B- and A-, the main RNA form, A-RNA, as well as in 13 other RNA conformations. High-level quantum mechanics calculations used a baseless dinucleoside phosphate as a molecular model, and the effect of solvent was included. The predicted J-couplings correlate reliably with available experimental data from the literature.

Solution structure of a purine rich hexaloop hairpin belonging to PGY/MDR1 mRNA and targeted by antisense oligonucleotides

A preferential target of antisense oligonucleotides directed against human PGY/MDR1 mRNA is a hairpin containing a stem with a G*U wobble pair, capped by the purine-rich 5'r(GGGAUG)3' hexaloop. This hairpin is studied by multidimensional NMR and restrained molecular dynamics, with special emphasis on the conformation of south sugars and non-standard phosphate linkages evidenced in both the stem and the loop. The hairpin is found to be highly structured. The G*U wobble pair, a strong counterion binding site, displays structural particularities that are characteristic of this type of mismatch. The upper part of the stem undergoes distortions that optimize its interactions with the beginning of the loop. The loop adopts a new fold in which the single-stranded GGGA purine tract is structured in A-like conformation stacked in continuity of the stem and displays an extensive hydrogen bonding surface for recognition. The remarkable hairpin stability results from classical inter- and intra-strand interactions reinforced by numerous hydrogen bonds involving unusual backbone conformations and ribose 2'-hydroxyl groups. Overall, this work emphasizes numerous features that account for the well-ordered structure of the whole hairpin and highlights the loop properties that facilitate interaction with antisense oligonucleotides.

Structure, dynamics, and elasticity of free 16s rRNA helix 44 studied by molecular dynamics simulations

Molecular dynamics (MD) simulations were employed to investigate the structure, dynamics, and local base-pair step deformability of the free 16S ribosomal helix 44 from Thermus thermophilus and of a canonical A-RNA double helix. While helix 44 is bent in the crystal structure of the small ribosomal subunit, the simulated helix 44 is intrinsically straight. It shows, however, substantial instantaneous bends that are isotropic. The spontaneous motions seen in simulations achieve large degrees of bending seen in the X-ray structure and would be entirely sufficient to allow the dynamics of the upper part of helix 44 evidenced by cryo-electron microscopic studies. Analysis of local base-pair step deformability reveals a patch of flexible steps in the upper part of helix 44 and in the area proximal to the bulge bases, suggesting that the upper part of helix 44 has enhanced flexibility. The simulations identify two conformational substates of the second bulge area (bottom part of the helix) with distinct base pairing. In agreement with nuclear magnetic resonance (NMR) and X-ray studies, a flipped out conformational substate of conserved 1492A is seen in the first bulge area. Molecular dynamics (MD) simulations reveal a number of reversible alpha-gamma backbone flips that correspond to transitions between two known A-RNA backbone families. The flipped substates do not cumulate along the trajectory and lead to a modest transient reduction of helical twist with no significant influence on the overall geometry of the duplexes. Despite their considerable flexibility, the simulated structures are very stable with no indication of substantial force field inaccuracies.

The building blocks and motifs of RNA architecture

RNA motifs can be defined broadly as recurrent structural elements containing multiple intramolecular RNA-RNA interactions, as observed in atomic-resolution RNA structures. They constitute the modular building blocks of RNA architecture, which is organized hierarchically. Recent work has focused on analyzing RNA backbone conformations to identify, define and search for new instances of recurrent motifs in X-ray structures. One current view asserts that recurrent RNA strand segments with characteristic backbone configurations qualify as independent motifs. Other considerations indicate that, to characterize modular motifs, one must take into account the larger structural context of such strand segments. This follows the biologically relevant motivation, which is to identify RNA structural characteristics that are subject to sequence constraints and that thus relate RNA architectures to sequences.

Automated extraction and classification of RNA tertiary structure cyclic motifs

A minimum cycle basis of the tertiary structure of a large ribosomal subunit (LSU) X-ray crystal structure was analyzed. Most cycles are small, as they are composed of 3- to 5 nt, and repeated across the LSU tertiary structure. We used hierarchical clustering to quantify and classify the 4 nt cycles. One class is defined by the GNRA tetraloop motif. The inspection of the GNRA class revealed peculiar instances in sequence. First is the presence of UA, CA, UC and CC base pairs that substitute the usual sheared GA base pair. Second is the revelation of GNR(X_n)A tetraloops, where X_n is bulged out of the classical GNRA structure, and of GN/RA formed by the two strands of interior-loops. We were able to unambiguously characterize the cycle classes using base stacking and base pairing annotations. The cycles identified correspond to small and cyclic motifs that compose most of the LSU RNA tertiary structure and contribute to its thermodynamic stability. Consequently, the RNA minimum cycles could well be used as the basic elements of RNA tertiary structure prediction methods.

Cations and hydration in catalytic RNA: molecular dynamics of the hepatitis delta virus ribozyme

The hepatitis delta virus (HDV) ribozyme is an RNA enzyme from the human pathogenic HDV. Cations play a crucial role in self-cleavage of the HDV ribozyme, by promoting both folding and chemistry. Experimental studies have revealed limited but intriguing details on the location and structural and catalytic functions of metal ions. Here, we analyze a total of approximately 200 ns of explicit-solvent molecular dynamics simulations to provide a complementary atomistic view of the binding of monovalent and divalent cations as well as water molecules to reaction precursor and product forms of the HDV ribozyme. Our simulations find that an Mg2+ cation binds stably, by both inner- and outer-sphere contacts, to the electronegative catalytic pocket of the reaction precursor, in a position to potentially support chemistry. In contrast, protonation of the catalytically involved C75 in the precursor or artificial placement of this Mg2+ into the product structure result in its swift expulsion from the active site. These findings are consistent with a concerted reaction mechanism in which C75 and hydrated Mg2+ act as general base and acid, respectively. Monovalent cations bind to the active site and elsewhere assisted by structurally bridging long-residency water molecules, but are generally delocalized.

The RNA Ontology Consortium: an open invitation to the RNA community

The aim of the RNA Ontology Consortium (ROC) is to create an integrated conceptual framework-an RNA Ontology (RO)-with a common, dynamic, controlled, and structured vocabulary to describe and characterize RNA sequences, secondary structures, three-dimensional structures, and dynamics pertaining to RNA function. The RO should produce tools for clear communication about RNA structure and function for multiple uses, including the integration of RNA electronic resources into the Semantic Web. These tools should allow the accurate description in computer-interpretable form of the coupling between RNA architecture, function, and evolution. The purposes for creating the RO are, therefore, (1) to integrate sequence and structural databases; (2) to allow different computational tools to interoperate; (3) to create powerful software tools that bring advanced computational methods to the bench scientist; and (4) to facilitate precise searches for all relevant information pertaining to RNA. For example, one initial objective of the ROC is to define, identify, and classify RNA structural motifs described in the literature or appearing in databases and to agree on a computer-interpretable definition for each of these motifs. To achieve these aims, the ROC will foster communication and promote collaboration among RNA scientists by coordinating frequent face-to-face workshops to discuss, debate, and resolve difficult conceptual issues. These meeting opportunities will create new directions at various levels of RNA research. The ROC will work closely with the PDB/NDB structural databases and the Gene, Sequence, and Open Biomedical Ontology Consortia to integrate the RO with existing biological ontologies to extend existing content while maintaining interoperability.

Molecular dynamics simulations of sarcin-ricin rRNA motif

Explicit solvent molecular dynamics (MD) simulations were carried out for sarcin-ricin domain (SRD) motifs from 23S (Escherichia coli) and 28S (rat) rRNAs. The SRD motif consists of GAGA tetraloop, G-bulged cross-strand A-stack, flexible region and duplex part. Detailed analysis of the overall dynamics, base pairing, hydration, cation binding and other SRD features is presented. The SRD is surprisingly static in multiple 25 ns long simulations and lacks any non-local motions, with root mean square deviation (r.m.s.d.) values between averaged MD and high-resolution X-ray structures of 1-1.4 A. Modest dynamics is observed in the tetraloop, namely, rotation of adenine in its apex and subtle reversible shift of the tetraloop with respect to the adjacent base pair. The deformed flexible region in low-resolution rat X-ray structure is repaired by simulations. The simulations reveal few backbone flips, which do not affect positions of bases and do not indicate a force field imbalance. Non-Watson-Crick base pairs are rigid and mediated by long-residency water molecules while there are several modest cation-binding sites around SRD. In summary, SRD is an unusually stiff rRNA building block. Its intrinsic structural and dynamical signatures seen in simulations are strikingly distinct from other rRNA motifs such as Loop E and Kink-turns.

Statistical analysis of RNA backbone

Local conformation is an important determinant of RNA catalysis and binding. The analysis of RNA conformation is particularly difficult due to the large number of degrees of freedom (torsion angles) per residue. Proteins, by comparison, have many fewer degrees of freedom per residue. In this work, we use and extend classical tools from statistics and signal processing to search for clusters in RNA conformational space. Results are reported both for scalar analysis, where each torsion angle is separately studied, and for vectorial analysis, where several angles are simultaneously clustered. Adapting techniques from vector quantization and clustering to the RNA structure, we find torsion angle clusters and RNA conformational motifs. We validate the technique using well-known conformational motifs, showing that the simultaneous study of the total torsion angle space leads to results consistent with known motifs reported in the literature and also to the finding of new ones.