Nuclear magnetic resonance structure of the III-IV-V three-way junction from the Varkud satellite ribozyme and identification of magnesium-binding sites using paramagnetic relaxation enhancement

Upstream Flanking Sequence Assists Folding of an RNA Thermometer

Many heat shock genes in bacteria are regulated through a class of temperature-sensitive stem-loop (SL) RNAs called RNA thermometers (RNATs). One of the most widely studied RNATs is the Repression Of heat Shock Expression (ROSE) element associated with expression of heat shock proteins. Located in the 5'UTR, the RNAT contains one to three auxiliary hairpins upstream of it. Herein, we address roles of these upstream SLs in the folding and function of an RNAT. Bradyrhizobium japonicum is a nitrogen-fixing bacterium that experiences a wide range of temperatures in the soil and contains ROSE elements, each having multiple upstream SLs. The 5'UTR of the messenger (mRNA) for heat shock protein A (hspA) in B. japonicum has an intricate secondary structure containing three SLs upstream of the RNAT SL. While structure-function studies of the hspA RNAT itself have been reported, it has been unclear if these auxiliary SLs contribute to the temperature-sensing function of the ROSE elements. Herein, we show that the full length (FL) sequence has several melting transitions indicating that the ROSE element unfolds in a non-two-state manner. The upstream SLs are more stable than the RNAT itself, and a variant with disrupted base pairing in the SL immediately upstream of the RNAT has little influence on the melting of the RNAT. On the basis of these results and modeling of the co-transcriptional folding of the ROSE element, we propose that the upstream SLs function to stabilize the transcript and aid proper folding and dynamics of the RNAT.

Acylation probing of "generic" RNA libraries reveals critical influence of loop constraints on reactivity

The reactivity of RNA 2'-OH acylation is broadly useful both in probing structure and in preparing conjugates. To date, this reactivity has been analyzed in limited sets of biological RNA sequences, leaving open questions of how reactivity varies inherently without regard to sequence in structured contexts. We constructed and probed "generic" structured RNA libraries using homogeneous loop sequences, employing deep sequencing to carry out a systematic survey of reactivity. We find a wide range of RNA reactivities among single-stranded sequences, with nearest neighbors playing substantial roles. Remarkably, certain small loops are found to be far more reactive on average (up to 4,000-fold) than single-stranded RNAs, due to conformational constraints that enhance reactivity. Among loops, we observe large variations in reactivity based on size, type, and position. The results lend insights into RNA designs for achieving high-efficiency local conjugation and provide new opportunities to refine structure analysis.

An integrative NMR-SAXS approach for structural determination of large RNAs defines the substrate-free state of a trans-cleaving Neurospora Varkud Satellite ribozyme

The divide-and-conquer strategy is commonly used for protein structure determination, but its applications to high-resolution structure determination of RNAs have been limited. Here, we introduce an integrative approach based on the divide-and-conquer strategy that was undertaken to determine the solution structure of an RNA model system, the Neurospora VS ribozyme. NMR and SAXS studies were conducted on a minimal trans VS ribozyme as well as several isolated subdomains. A multi-step procedure was used for structure determination that first involved pairing refined NMR structures with SAXS data to obtain structural subensembles of the various subdomains. These subdomain structures were then assembled to build a large set of structural models of the ribozyme, which was subsequently filtered using SAXS data. The resulting NMR-SAXS structural ensemble shares several similarities with the reported crystal structures of the VS ribozyme. However, a local structural difference is observed that affects the global fold by shifting the relative orientation of the two three-way junctions. Thus, this finding highlights a global conformational change associated with substrate binding in the VS ribozyme that is likely critical for its enzymatic activity. Structural studies of other large RNAs should benefit from similar integrative approaches that allow conformational sampling of assembled fragments.

Structure and Mechanical Stabilities of the Three-Way Junction Motifs in Prohead RNA

The core structure of phi29 prohead RNA (pRNA) is composed of three major helices organized into three-way junction pRNA (3WJ-pRNA) and has stout structural rigidity along the coaxial helices. Prohead RNAs of the other Bacillus subtilis bacteriophages such as GA1 and SF5 share similar secondary structure and function with phi29; whether these pRNAs have similar mechanical rigidity remains to be elucidated. In this study, we constructed the tertiary structures of GA1 and SF5 3WJ-pRNAs by comparative modeling. Both GA1 and SF5 3WJ-pRNAs adopt a similar structure, in which three helices are organized as the three-way junction and two of the three helices are stacked coaxially. Moreover, detailed structural features of GA1 and SF5 3WJ-pRNAs are also similar to those of phi29 3WJ-pRNA: all of the bases of the coaxial helices are paired, and all of the adenines in the junction region are paired, which eliminates the interference of A-minor tertiary interactions. Hence, the coaxial helices tightly join to each other, and the major groove between them is very narrow. Two Mg²⁺ ions can thus fit into this major groove and form double Mg clamps. A steered molecular dynamics simulation was used to study the mechanical properties of these 3WJ-pRNAs. Both GA1 and SF5 3WJ-pRNAs show strong resistance to applied force in the direction of their coaxial helices. Such mechanical stability can be attributed to the Mg clamps.

The Varkud Satellite Ribozyme: A Thirty-Year Journey through Biochemistry, Crystallography, and Computation

The discovery of catalytic RNAs or ribozymes introduced a new class of enzymes to biology. In addition to their increasingly important roles in modern life, ribozymes are key players in the RNA World hypothesis, which posits that life started or flourished with RNA supporting both genetic and enzymatic functions. Therefore, investigations into the mechanisms of ribozyme function provide an exciting opportunity to examine the foundational principles of biological catalysis. Ribozymes are also attractive model systems to investigate the relationship between structure and function in RNA. Endonucleolytic ribozymes represent the largest class of catalytic RNA, of which the Varkud satellite (VS) ribozyme is structurally the most complex. The last ribozyme to be discovered by accident, the VS ribozyme had eluded structural determination for over two decades. When we solved the first crystal structures of the VS ribozyme, an extensive body of biochemical and biophysical data had accumulated over the years with which we could evaluate the functional relevance of the structure. Conversely, the structures provided a new perspective from which to reexamine the functional data and test new hypotheses. The VS ribozyme is organized in a modular fashion where independently folding domains assemble into the active conformation of the ribozyme via three-way junctions. Structures of the VS ribozyme in complex with its substrate at different stages of activation enabled us to map the structural reorganization of the substrate that must precede catalysis. In addition to defining the global architecture of the RNA, the essential interactions between the substrate and catalytic domains, and the rearrangements in the substrate prior to catalysis, these structures provided detailed snapshots of the ribozyme active site, revealing potential catalytic interactions. High resolution structures of the active site bolstered the view that the catalytic mechanism involved nucleobase-mediated general acid-base catalysis and uncovered additional catalytic interactions between the cleavage site and catalytic residues. Informed by the crystal structures of the VS ribozyme, an integrated experimental and computational approach identified the key players and essential interactions that define the active site of the ribozyme. This confluence of biochemical, structural, and computational studies revealed the catalytic mechanism of the ribozyme at unprecedented detail. Additionally, comparative analyses of the active site structures of the VS ribozyme and other nucleic acid-based endoribonucleases revealed common architectural motifs and strikingly similar catalytic strategies. In this Account, we document the progress of VS ribozyme research starting from its discovery and extending to the elucidation of its detailed catalytic mechanism 30 years later.

Charge-driven condensation of RNA and proteins suggests broad role of phase separation in cytoplasmic environments

Phase separation processes are increasingly being recognized as important organizing mechanisms of biological macromolecules in cellular environments. Well-established drivers of phase separation are multi-valency and intrinsic disorder. Here, we show that globular macromolecules may condense simply based on electrostatic complementarity. More specifically, phase separation of mixtures between RNA and positively charged proteins is described from a combination of multiscale computer simulations with microscopy and spectroscopy experiments. Phase diagrams were mapped out as a function of molecular concentrations in experiment and as a function of molecular size and temperature via simulations. The resulting condensates were found to retain at least some degree of internal dynamics varying as a function of the molecular composition. The results suggest a more general principle for phase separation that is based primarily on electrostatic complementarity without invoking polymer properties as in most previous studies. Simulation results furthermore suggest that such phase separation may occur widely in heterogenous cellular environment between nucleic acid and protein components.

Confluence of theory and experiment reveals the catalytic mechanism of the Varkud satellite ribozyme

The Varkud satellite ribozyme catalyses site-specific RNA cleavage and ligation, and serves as an important model system to understand RNA catalysis. Here, we combine stereospecific phosphorothioate substitution, precision nucleobase mutation and linear free-energy relationship measurements with molecular dynamics, molecular solvation theory and ab initio quantum mechanical/molecular mechanical free-energy simulations to gain insight into the catalysis. Through this confluence of theory and experiment, we unify the existing body of structural and functional data to unveil the catalytic mechanism in unprecedented detail, including the degree of proton transfer in the transition state. Further, we provide evidence for a critical Mg²⁺ in the active site that interacts with the scissile phosphate and anchors the general base guanine in position for nucleophile activation. This novel role for Mg²⁺ adds to the diversity of known catalytic RNA strategies and unifies functional features observed in the Varkud satellite, hairpin and hammerhead ribozyme classes.

Limits in accuracy and a strategy of RNA structure prediction using experimental information

RNA structural complexity and flexibility present a challenge for computational modeling efforts. Experimental information and bioinformatics data can be used as restraints to improve the accuracy of RNA tertiary structure prediction. Regarding utilization of restraints, the fundamental questions are: (i) What is the limit in prediction accuracy that one can achieve with arbitrary number of restraints? (ii) Is there a strategy for selection of the minimal number of restraints that would result in the best structural model? We address the first question by testing the limits in prediction accuracy using native contacts as restraints. To address the second question, we develop an algorithm based on the distance variation allowed by secondary structure (DVASS), which ranks restraints according to their importance to RNA tertiary structure prediction. We find that due to kinetic traps, the greatest improvement in the structure prediction accuracy is achieved when we utilize only 40-60% of the total number of native contacts as restraints. When the restraints are sorted by DVASS algorithm, using only the first 20% ranked restraints can greatly improve the prediction accuracy. Our findings suggest that only a limited number of strategically selected distance restraints can significantly assist in RNA structure modeling.

Dynamics of metal ions around an RNA molecule

Noncoding RNA molecules take part in many biological processes, while metal ions play crucial roles in helping RNAs to perform their functions. However, the statics and dynamics of these metal ions around RNA molecules are still not well understood. In this work, we report a detailed molecular dynamics study of the type-I preQ_{1}-bound riboswitch aptamer domain (PRAD) at different ionic conditions (K^{+}, Na^{+}, and Mg^{2+}). The results show that the structural properties and flexibility of the PRAD molecule greatly influence the distributions and dynamics of metal ions around it. Simultaneously, Na^{+} ions show a stronger competitiveness with Mg^{2+} ions than K^{+} ions, and the three types of metal ions have different modes of interaction with the RNA molecule. Furthermore, we have also investigated specific binding sites of metal ions on the PRAD molecule and found that the dynamics and hydration structures of metal ions located at the ion-binding sites were obviously affected by the RNA structure near these ion-binding sites. These results may be useful to understand the role of the metal ions in noncoding RNA functions.

A Small Helical Bundle Prepares Primer Synthesis by Binding Two Nucleotides that Enhance Sequence-Specific Recognition of the DNA Template

Primases have a fundamental role in DNA replication. They synthesize a primer that is then extended by DNA polymerases. Archaeoeukaryotic primases require for synthesis a catalytic and an accessory domain, the exact contribution of the latter being unresolved. For the pRN1 archaeal primase, this domain is a 115-amino acid helix bundle domain (HBD). Our structural investigations of this small HBD by liquid- and solid-state nuclear magnetic resonance (NMR) revealed that only the HBD binds the DNA template. DNA binding becomes sequence-specific after a major allosteric change in the HBD, triggered by the binding of two nucleotide triphosphates. The spatial proximity of the two nucleotides and the DNA template in the quaternary structure of the HBD strongly suggests that this small domain brings together the substrates to prepare the first catalytic step of primer synthesis. This efficient mechanism is likely general for all archaeoeukaryotic primases.

Thermodynamic investigation of kissing-loop interactions

Kissing loop interactions (KLIs) are a common motif that is critical in retroviral dimerization, viroid replication, mRNA, and riboswitches. In addition, KLIs are currently used in a variety of biotechnology applications, such as in aptamer sensors, RNA scaffolds and to stabilize vaccines for therapeutics. Here we describe the thermodynamics of a basic intramolecular DNA capable of engaging in a KLI, consisting of two hairpins connected by a flexible linker. Each hairpin loop has a five-nucleotide complementary sequence theoretically capable of engaging in a KLI. On either side of each loop is two thymines which will not engage in kissing but are present to provide more flexibility and optimal KLI positioning. Our results suggest that the KLI occurs even at physiological salt levels, and that the KLI does not alter the thermodynamics and stability of the two stem structures. The KLI does not involve all five nucleotides, or at least each base-pair stack is not making full contact. Adding a second strand complementary to the bottom of the kissing complex removes flexibility and causes destabilization of the stems. The KLI of this less flexible complex is maintained but the T_M is reduced, indicating an entopic penalty to its formation.

Principles for targeting RNA with drug-like small molecules

Recent studies have indicated the potential to develop small-molecule drugs that act on RNA targets, leading to burgeoning interest in the field. This article discusses general principles for discovering small-molecule drugs that target RNA and argues that the overarching challenge is to identify appropriate target structures in disease-causing RNAs that have high information content and, consequently, appropriate ligand-binding pockets.

RNA molecules are essential for cellular information transfer and gene regulation, and RNAs have been implicated in many human diseases. Messenger and non-coding RNAs contain highly structured elements, and evidence suggests that many of these structures are important for function. Targeting these RNAs with small molecules offers opportunities to therapeutically modulate numerous cellular processes, including those linked to 'undruggable' protein targets. Despite this promise, there is currently only a single class of human-designed small molecules that target RNA used clinically — the linezolid antibiotics. However, a growing number of small-molecule RNA ligands are being identified, leading to burgeoning interest in the field. Here, we discuss principles for discovering small-molecule drugs that target RNA and argue that the overarching challenge is to identify appropriate target structures — namely, in disease-causing RNAs that have high information content and, consequently, appropriate ligand-binding pockets. If focus is placed on such druggable binding sites in RNA, extensive knowledge of the typical physicochemical properties of drug-like small molecules could then enable small-molecule drug discovery for RNA targets to become (only) roughly as difficult as for protein targets.

Increased Flexibility between Stems of Intramolecular Three-Way Junctions by the Insertion of Bulges

Intramolecular junctions are a ubiquitous structure within DNA and RNA; three-way junctions in particular have high strain around the junction because of the lack of flexibility, preventing the junctions from adopting conformations that would allow for optimal folding. In this work, we used a combination of calorimetric and spectroscopic techniques to study the unfolding of four intramolecular three-way junctions. The control three-way junction, 3H, has the sequence d(GAAATTGCGCT5GCGCGTGCT5GCACAATTTC), which has three arms of different sequences. We studied three other three-way junctions in which one (2HS1H), two (HS12HS1), and three (HS1HS1HS1) cytosine bulges were placed at the junction to allow the arms to adopt a wider range of conformations that may potentially relieve strain. Through calorimetric studies, it was concluded that bulges produce only minor effects on the enthalpic and thermal stability at physiological salt concentrations for 2HS1H and HS1HS1HS1. HS12HS1 displays the strongest effect, with the GTGC stem lacking a defined transition. In addition to unfolding thermodynamics, the differential binding of counterions, water, and protons was determined. It was found that with each bulge, there was a large increase in the binding of counterions; this correlated with a decrease in the immobilization of structural water molecules. The increase in counterion uptake upon folding likely displaces binding of structural water, which is measured by the osmotic stress method, in favor of electrostricted waters. The cytosine bulges do not affect the binding of protons; this finding indicates that the bulges are not forming base-triplet stacks. These results indicate that bulges in junctions do not affect the unfolding profile or the enthalpy of oligonucleotides but do affect the number and amount of molecules immobilized by the junction.

Influence of Na+ and Mg2+ ions on RNA structures studied with molecular dynamics simulations

The structure of ribonucleic acid (RNA) polymers is strongly dependent on the presence of, in particular Mg2+ cations to stabilize structural features. Only in high-resolution X-ray crystallography structures can ions be identified reliably. Here, we perform molecular dynamics simulations of 24 RNA structures with varying ion concentrations. Twelve of the structures were helical and the others complex folded. The aim of the study is to predict ion positions but also to evaluate the impact of different types of ions (Na+ or Mg2+) and the ionic strength on structural stability and variations of RNA. As a general conclusion Mg2+ is found to conserve the experimental structure better than Na+ and, where experimental ion positions are available, they can be reproduced with reasonable accuracy. If a large surplus of ions is present the added electrostatic screening makes prediction of binding-sites less reproducible. Distinct differences in ion-binding between helical and complex folded structures are found. The strength of binding (ΔG‡ for breaking RNA atom-ion interactions) is found to differ between roughly 10 and 26 kJ/mol for the different RNA atoms. Differences in stability between helical and complex folded structures and of the influence of metal ions on either are discussed.

Mapping the Universe of RNA Tetraloop Folds

We report a map of RNA tetraloop conformations constructed by calculating pairwise distances among all experimentally determined four-nucleotide hairpin loops. Tetraloops with similar structures are clustered together and, as expected, the two largest clusters are the canonical GNRA and UNCG folds. We identify clusters corresponding to known tetraloop folds such as GGUG, RNYA, AGNN, and CUUG. These clusters are represented in a simple two-dimensional projection that recapitulates the relationship among the different folds. The cluster analysis also identifies 20 novel tetraloop folds that are peculiar to specific positions in ribosomal RNAs and that are stabilized by tertiary interactions. In our RNA tetraloop database we find a significant number of non-GNRA and non-UNCG sequences adopting the canonical GNRA and UNCG folds. Conversely, we find a significant number of GNRA and UNCG sequences adopting non-GNRA and non-UNCG folds. Our analysis demonstrates that there is not a simple one-to-one, but rather a many-to-many mapping between tetraloop sequence and tetraloop fold.

Insights into RNA structure and dynamics from recent NMR and X-ray studies of the Neurospora Varkud satellite ribozyme

Despite the large number of noncoding RNAs and their importance in several biological processes, our understanding of RNA structure and dynamics at atomic resolution is still very limited. Like many other RNAs, the Neurospora Varkud satellite (VS) ribozyme performs its functions through dynamic exchange of multiple conformational states. More specifically, the VS ribozyme recognizes and cleaves its stem-loop substrate via a mechanism that involves several structural transitions within its stem-loop substrate. The recent publications of high-resolution structures of the VS ribozyme, obtained by NMR spectroscopy and X-ray crystallography, offer an opportunity to integrate the data and closely examine the structural and dynamic properties of this model RNA system. Notably, these investigations provide a valuable example of the divide-and-conquer strategy for structural and dynamic characterization of a large RNA, based on NMR structures of several individual subdomains. The success of this divide-and-conquer approach reflects the modularity of RNA architecture and the great care taken in identifying the independently-folding modules. Together with previous biochemical and biophysical characterizations, the recent NMR and X-ray studies provide a coherent picture into how the VS ribozyme recognizes its stem-loop substrate. Such in-depth characterization of this RNA enzyme will serve as a model for future structural and engineering studies of dynamic RNAs and may be particularly useful in planning divide-and-conquer investigations. WIREs RNA 2017, 8:e1421. doi: 10.1002/wrna.1421 For further resources related to this article, please visit the WIREs website.

Rational engineering of the Neurospora VS ribozyme to allow substrate recognition via different kissing-loop interactions

The Neurospora VS ribozyme is a catalytic RNA that has the unique ability to specifically recognize and cleave a stem-loop substrate through formation of a highly stable kissing-loop interaction (KLI). In order to explore the engineering potential of the VS ribozyme to cleave alternate substrates, we substituted the wild-type KLI by other known KLIs using an innovative engineering method that combines rational and combinatorial approaches. A bioinformatic search of the protein data bank was initially performed to identify KLIs that are structurally similar to the one found in the VS ribozyme. Next, substrate/ribozyme (S/R) pairs that incorporate these alternative KLIs were kinetically and structurally characterized. Interestingly, several of the resulting S/R pairs allowed substrate cleavage with substantial catalytic efficiency, although with reduced activity compared to the reference S/R pair. Overall, this study describes an innovative approach for RNA engineering and establishes that the KLI of the trans VS ribozyme can be adapted to cleave other folded RNA substrates.

Crystal structure of the Varkud satellite ribozyme

Varkud Satellite (VS) ribozyme mediates rolling circle replication of a plasmid found in the Neurospora mitochondria. We report crystal structures of this ribozyme at 3.1Å resolution, revealing an intertwined dimer formed by an exchange of substrate helices. Within each protomer, an arrangement of three-way helical junctions organizes seven helices into a global fold that creates a docking site for the substrate helix of the other protomer, resulting in the formation of two active sites in trans. This mode of RNA-RNA association resembles the process of domain swapping in proteins and has implications for RNA regulation and evolution. Within each active site, adenine and guanine nucleobases abut the scissile phosphate, poised to serve direct roles in catalysis. Similarities to the active sites of the hairpin and hammerhead ribozymes highlight the functional significance of active site features, underscore the ability of RNA to access functional architectures from distant regions of sequence space, and suggest convergent evolution.

The NMR structure of the II-III-VI three-way junction from the Neurospora VS ribozyme reveals a critical tertiary interaction and provides new insights into the global ribozyme structure

As part of an effort to structurally characterize the complete Neurospora VS ribozyme, NMR solution structures of several subdomains have been previously determined, including the internal loops of domains I and VI, the I/V kissing-loop interaction and the III-IV-V junction. Here, we expand this work by determining the NMR structure of a 62-nucleotide RNA (J236) that encompasses the VS ribozyme II-III-VI three-way junction and its adjoining stems. In addition, we localize Mg(2+)-binding sites within this structure using Mn(2+)-induced paramagnetic relaxation enhancement. The NMR structure of the J236 RNA displays a family C topology with a compact core stabilized by continuous stacking of stems II and III, a cis WC/WC G•A base pair, two base triples and two Mg(2+) ions. Moreover, it reveals a remote tertiary interaction between the adenine bulges of stems II and VI. Additional NMR studies demonstrate that both this bulge-bulge interaction and Mg(2+) ions are critical for the stable folding of the II-III-VI junction. The NMR structure of the J236 RNA is consistent with biochemical studies on the complete VS ribozyme, but not with biophysical studies performed with a minimal II-III-VI junction that does not contain the II-VI bulge-bulge interaction. Together with previous NMR studies, our findings provide important new insights into the three-dimensional architecture of this unique ribozyme.

Fitness Landscapes of Functional RNAs

The notion of fitness landscapes, a map between genotype and fitness, was proposed more than 80 years ago. For most of this time data was only available for a few alleles, and thus we had only a restricted view of the whole fitness landscape. Recently, advances in genetics and molecular biology allow a more detailed view of them. Here we review experimental and theoretical studies of fitness landscapes of functional RNAs, especially aptamers and ribozymes. We find that RNA structures can be divided into critical structures, connecting structures, neutral structures and forbidden structures. Such characterisation, coupled with theoretical sequence-to-structure predictions, allows us to construct the whole fitness landscape. Fitness landscapes then can be used to study evolution, and in our case the development of the RNA world.

Structural Analysis of Protein-RNA Complexes in Solution Using NMR Paramagnetic Relaxation Enhancements

Biological activity in the cell is predominantly mediated by large multiprotein and protein-nucleic acid complexes that act together to ensure functional fidelity. Nuclear magnetic resonance (NMR) spectroscopy is the only method that can provide information for high-resolution three-dimensional structures and the conformational dynamics of these complexes in solution. Mapping of binding interfaces and molecular interactions along with the characterization of conformational dynamics is possible for very large protein complexes. In contrast, de novo structure determination by NMR becomes very time consuming and difficult for protein complexes larger than 30 kDa as data are noisy and sparse. Fortunately, high-resolution structures are often available for individual domains or subunits of a protein complex and thus sparse data can be used to define their arrangement and dynamics within the assembled complex. In these cases, NMR can therefore be efficiently combined with complementary solution techniques, such as small-angle X-ray or neutron scattering, to provide a comprehensive description of the structure and dynamics of protein complexes in solution. Particularly useful are NMR-derived paramagnetic relaxation enhancements (PREs), which provide long-range distance restraints (ca. 20Å) for structural analysis of large complexes and also report on conformational dynamics in solution. Here, we describe the use of PREs from sample production to structure calculation, focusing on protein-RNA complexes. On the basis of recent examples from our own research, we demonstrate the utility, present protocols, and discuss potential pitfalls when using PREs for studying the structure and dynamic features of protein-RNA complexes.