Specificity of RNA-RNA helix recognition

Predicting RNA structure and dynamics with deep learning and solution scattering

Advanced deep learning and statistical methods can predict structural models for RNA molecules. However, RNAs are flexible, and it remains difficult to describe their macromolecular conformations in solutions where varying conditions can induce conformational changes. Small-angle x-ray scattering (SAXS) in solution is an efficient technique to validate structural predictions by comparing the experimental SAXS profile with those calculated from predicted structures. There are two main challenges in comparing SAXS profiles to RNA structures: the absence of cations essential for stability and charge neutralization in predicted structures and the inadequacy of a single structure to represent RNA's conformational plasticity. We introduce a solution conformation predictor for RNA (SCOPER) to address these challenges. This pipeline integrates kinematics-based conformational sampling with the innovative deep learning model, IonNet, designed for predicting Mg2+ ion binding sites. Validated through benchmarking against 14 experimental data sets, SCOPER significantly improved the quality of SAXS profile fits by including Mg2+ ions and sampling of conformational plasticity. We observe that an increased content of monovalent and bivalent ions leads to decreased RNA plasticity. Therefore, carefully adjusting the plasticity and ion density is crucial to avoid overfitting experimental SAXS data. SCOPER is an efficient tool for accurately validating the solution state of RNAs given an initial, sufficiently accurate structure and provides the corrected atomistic model, including ions.

Predicting RNA Structure and Dynamics with Deep Learning and Solution Scattering

Advanced deep learning and statistical methods can predict structural models for RNA molecules. However, RNAs are flexible, and it remains difficult to describe their macromolecular conformations in solutions where varying conditions can induce conformational changes. Small-angle X-ray scattering (SAXS) in solution is an efficient technique to validate structural predictions by comparing the experimental SAXS profile with those calculated from predicted structures. There are two main challenges in comparing SAXS profiles to RNA structures: the absence of cations essential for stability and charge neutralization in predicted structures and the inadequacy of a single structure to represent RNA's conformational plasticity. We introduce Solution Conformation Predictor for RNA (SCOPER) to address these challenges. This pipeline integrates kinematics-based conformational sampling with the innovative deep-learning model, IonNet, designed for predicting Mg2+ ion binding sites. Validated through benchmarking against fourteen experimental datasets, SCOPER significantly improved the quality of SAXS profile fits by including Mg2+ ions and sampling of conformational plasticity. We observe that an increased content of monovalent and bivalent ions leads to decreased RNA plasticity. Therefore, carefully adjusting the plasticity and ion density is crucial to avoid overfitting experimental SAXS data. SCOPER is an efficient tool for accurately validating the solution state of RNAs given an initial, sufficiently accurate structure and provides the corrected atomistic model, including ions.

The central role of transfer RNAs in mistranslation

Transfer RNAs (tRNA) are essential small non-coding RNAs that enable the translation of genomic information into proteins in all life forms. The principal function of tRNAs is to bring amino acid building blocks to the ribosomes for protein synthesis. In the ribosome, tRNAs interact with messenger RNA (mRNA) to mediate the incorporation of amino acids into a growing polypeptide chain following the rules of the genetic code. Accurate interpretation of the genetic code requires tRNAs to carry amino acids matching their anticodon identity and decode the correct codon on mRNAs. Errors in these steps cause the translation of codons with the wrong amino acids (mistranslation), compromising the accurate flow of information from DNA to proteins. Accumulation of mutant proteins due to mistranslation jeopardizes proteostasis and cellular viability. However, the concept of mistranslation is evolving, with increasing evidence indicating that mistranslation can be used as a mechanism for survival and acclimatization to environmental conditions. In this review, we discuss the central role of tRNAs in modulating translational fidelity through their dynamic and complex interplay with translation factors. We summarize recent discoveries of mistranslating tRNAs and describe the underlying molecular mechanisms and the specific conditions and environments that enable and promote mistranslation.

Endogenous trans-translation structure visualizes the decoding of the first tmRNA alanine codon

Ribosomes stall on truncated or otherwise damaged mRNAs. Bacteria rely on ribosome rescue mechanisms to replenish the pool of ribosomes available for translation. Trans-translation, the main ribosome-rescue pathway, uses a circular hybrid transfer-messenger RNA (tmRNA) to restart translation and label the resulting peptide for degradation. Previous studies have visualized how tmRNA and its helper protein SmpB interact with the stalled ribosome to establish a new open reading frame. As tmRNA presents the first alanine codon via a non-canonical mRNA path in the ribosome, the incoming alanyl-tRNA must rearrange the tmRNA molecule to read the codon. Here, we describe cryo-EM analyses of an endogenous Escherichia coli ribosome-tmRNA complex with tRNAAla accommodated in the A site. The flexible adenosine-rich tmRNA linker, which connects the mRNA-like domain with the codon, is stabilized by the minor groove of the canonically positioned anticodon stem of tRNAAla. This ribosome complex can also accommodate a tRNA near the E (exit) site, bringing insights into the translocation and dissociation of the tRNA that decoded the defective mRNA prior to tmRNA binding. Together, these structures uncover a key step of ribosome rescue, in which the ribosome starts translating the tmRNA reading frame.

Geometric deep learning for the prediction of magnesium-binding sites in RNA structures

Magnesium ions (Mg2+) are essential for the folding, functional expression, and structural stability of RNA molecules. However, predicting Mg2+-binding sites in RNA molecules based solely on RNA structures is still challenging. The molecular surface, characterized by a continuous shape with geometric and chemical properties, is important for RNA modelling and carries essential information for understanding the interactions between RNAs and Mg2+ ions. Here, we propose an approach named RNA-magnesium ion surface interaction fingerprinting (RMSIF), a geometric deep learning-based conceptual framework to predict magnesium ion binding sites in RNA structures. To evaluate the performance of RMSIF, we systematically enumerated decoy Mg2+ ions across a full-space grid within the range of 2 to 10 Å from the RNA molecule and made predictions accordingly. Visualization techniques were used to validate the prediction results and calculate success rates. Comparative assessments against state-of-the-art methods like MetalionRNA, MgNet, and Metal3DRNA revealed that RMSIF achieved superior success rates and accuracy in predicting Mg2+-binding sites. Additionally, in terms of the spatial distribution of Mg2+ ions within the RNA structures, a majority were situated in the deep grooves, while a minority occupied the shallow grooves. Collectively, the conceptual framework developed in this study holds promise for advancing insights into drug design, RNA co-transcriptional folding, and structure prediction.

Context-dependence of T-loop Mediated Long-range RNA Tertiary Interactions

The architecture and folding of complex RNAs is governed by a limited set of highly recurrent structural motifs that form long-range tertiary interactions. One of these motifs is the T-loop, which was first identified in tRNA but is broadly distributed across biological RNAs. While the T-loop has been examined in detail in different biological contexts, the various receptors that it interacts with are not as well defined. In this study, we use a cell-based genetic screen in concert with bioinformatic analysis to examine three different, but related, T-loop receptor motifs found in the flavin mononucleotide (FMN) and cobalamin (Cbl) riboswitches. As a host for different T-loop receptors, we employed the env8 class-II Cbl riboswitch, an RNA that uses two T-loop motifs for both folding and supporting the ligand binding pocket. A set of libraries was created in which select nucleotides that participate in the T-loop/T-loop receptor (TL/TLR) interaction were fully randomized. Library members were screened for their ability to support Cbl-dependent expression of a reporter gene. While T-loops appear to be variable in sequence, we find that the functional sequence space is more restricted in the Cbl riboswitch, suggesting that TL/TLR interactions are context dependent. Our data reveal clear sequence signatures for the different types of receptor motifs that align with phylogenic analysis of these motifs in the FMN and Cbl riboswitches. Finally, our data suggest the functional contribution of various nucleobase-mediated long-range interactions within the riboswitch subclass of TL/TLR interactions that are distinct from those found in other RNAs.

Features and Functions of the A-Minor Motif, the Most Common Motif in RNA Structure

A-minor motifs are RNA tertiary structure motifs that generally involve a canonical base pair and an adenine base forming hydrogen bonds with the minor groove of the base pair. Such motifs are among the most common tertiary interactions in known RNA structures, comparable in number with the non-canonical base pairs. They are often found in functionally important regions of non-coding RNAs and, in particular, play a central role in protein synthesis. Here, we review local variations of the A-minor geometry and discuss difficulties associated with their annotation, as well as various structural contexts and common A-minor co-motifs, and diverse functions of A-minors in various processes in a living cell.

RNA stabilization by a poly(A) tail 3'-end binding pocket and other modes of poly(A)-RNA interaction

Polyadenylate [poly(A)] tail addition to the 3' end of a wide range of RNAs is a highly conserved modification that plays a central role in cellular RNA function. Elements for nuclear expression (ENEs) are cis-acting RNA elements that stabilize poly(A) tails by sequestering them in RNA triplex structures. A crystal structure of a double ENE from a rice hAT transposon messenger RNA complexed with poly(A)₂₈ at a resolution of 2.89 angstroms reveals multiple modes of interaction with poly(A), including major-groove triple helices, extended minor-groove interactions with RNA double helices, a quintuple-base motif that transitions poly(A) from minor-groove associations to major-groove triple helices, and a poly(A) 3'-end binding pocket. Our findings both expand the repertoire of motifs involved in long-range RNA interactions and provide insights into how polyadenylation can protect an RNA's extreme 3' end.

Incorporating a Thiophosphate Modification into a Common RNA Tetraloop Motif Causes an Unanticipated Stability Boost

GNRA (N = A, C, G, or U; R = A or G) tetraloops are common RNA secondary structural motifs and feature a phosphate stacked atop a nucleobase. The rRNA sarcin/ricin loop (SRL) is capped by GApGA, and the phosphate p stacks on G. We recently found that regiospecific incorporation of a single dithiophosphate (PS2) but not a monothiophosphate (PSO) instead of phosphate in the backbone of RNA aptamers dramatically increases the binding affinity for their targets. In the RNA:thrombin complex, the key contribution to the 1000-fold tighter binding stems from an edge-on contact between PS2 and a phenylalanine ring. Here we investigated the consequences of replacing the SRL phosphate engaged in a face-on interaction with guanine with either PS2 or PSO for stability. We found that PS2···G and Rp-PSO···G contacts stabilize modified SRLs compared to the parent loop to unexpected levels: up to 6.3 °C in melting temperature T_m and -4.7 kcal/mol in ΔΔG°. Crystal structures demonstrate that the vertical distance to guanine for the closest sulfur is just 0.05 Å longer on average compared to that of oxygen despite the larger van der Waals radius of the former (1.80 Å for S vs 1.52 Å for O). The higher stability is enthalpy-based, and the negative charge as assessed by a neutral methylphosphonate modification plays only a minor role. Quantum mechanical/molecular mechanical calculations are supportive of favorable dispersion attraction interactions by sulfur making the dominant contribution. A stacking interaction between phosphate and guanine (SRL) or uracil (U-turn) is also found in newly classified RNA tetraloop families besides GNRA.

Age-Dependent and Tissue-Specific Alterations in the rDNA Clusters of the Panax ginseng C. A. Meyer Cultivated Cell Lines

Long-term cultivation of Panax ginseng cell lines leads to a decreasing synthesis of the biologically active substances used in traditional medicine. To gain insight into the cellular mechanisms which may influence this process, we analyzed variations within the rDNA cluster of the Oriental ginseng cell lines. The cell lines were cultivated for 6 and 24 years; the number of nucleoli and chromosomes was analyzed. The complete 18S rDNA sequences were cloned and sequenced. The nucleotide polymorphism and phylogenetic relations of the sequences were analyzed, and the secondary structures for separate 18S rRNA regions were modeled. The 18S rDNA accumulated mutations during cell cultivation that correlate well with an increase in the number of chromosomes and nucleoli. The patterns of nucleotide diversity are culture-specific and the increasing polymorphism associates with cytosine methylation sites. The secondary structures of some 18S rRNA regions and their interaction can alter during cultivation. The phylogenetic tree topologies are particular for each cell line.The observed alterations in rDNA clusters are associated with a somaclonal variation, leading to changes in the pattern of intracellular synthesis during cell cultivation. The identified divergent rRNAs could provide additional gene expression regulation in P. ginseng cells by forming heterogeneous ribosomes.

The SSU rRNA secondary structures of the Plagiorchiida species (Digenea), its applications in systematics and evolutionary inferences

The small subunit ribosomal RNA (SSU rRNA) is widely used phylogenetic marker in broad groups of organisms and its secondary structure increasingly attracts the attention of researchers as supplementary tool in sequence alignment and advanced phylogenetic studies. Its comparative analysis provides a great contribution to evolutionary biology, allowing find out how the SSU rRNA secondary structure originated, developed and evolved. Herein, we provide the first data on the putative SSU rRNA secondary structures of the Plagiorchiida species. The structures were found to be quite conserved across broad range of species studied, well compatible with those of others eukaryotic SSU rRNA and possessed some peculiarities: cross-shaped structure of the ES6b, additional shortened ES6c2 helix, and elongated ES6a helix and h39 + ES9 region. The secondary structures of variable regions ES3 and ES7 appeared to be tissue-specific while ES6 and ES9 were specific at a family level allowing considering them as promising markers for digenean systematics. Their uniqueness more depends on the length than on the nucleotide diversity of primary sequences which evolutionary rates well differ. The findings have important implications for understanding rRNA evolution, developing molecular taxonomy and systematics of Plagiorchiida as well as for constructing new anthelmintic drugs.

Face-time with TAR: Portraits of an HIV-1 RNA with diverse modes of effector recognition relevant for drug discovery

Small molecules and short peptides that potently and selectively bind RNA are rare, making the molecular structures of these complexes highly exceptional. Accordingly, several recent investigations have provided unprecedented structural insights into how peptides and proteins recognize the HIV-1 transactivation response (TAR) element, a 59-nucleotide-long, noncoding RNA segment in the 5' long terminal repeat region of viral transcripts. Here, we offer an integrated perspective on these advances by describing earlier progress on TAR binding to small molecules, and by drawing parallels to recent successes in the identification of compounds that target the hepatitis C virus internal ribosome entry site (IRES) and the flavin-mononucleotide riboswitch. We relate this work to recent progress that pinpoints specific determinants of TAR recognition by: (i) viral Tat proteins, (ii) an innovative lab-evolved TAR-binding protein, and (iii) an ultrahigh-affinity cyclic peptide. New structural details are used to model the TAR-Tat-super-elongation complex (SEC) that is essential for efficient viral transcription and represents a focal point for antiviral drug design. A key prediction is that the Tat transactivation domain makes modest contacts with the TAR apical loop, whereas its arginine-rich motif spans the entire length of the TAR major groove. This expansive interface has significant implications for drug discovery and design, and it further suggests that future lab-evolved proteins could be deployed to discover steric restriction points that block Tat-mediated recruitment of the host SEC to HIV-1 TAR.

Gene regulation by a glycine riboswitch singlet uses a finely tuned energetic landscape for helical switching

Riboswitches contain structured aptamer domains that, upon ligand binding, facilitate helical switching in their downstream expression platforms to alter gene expression. To fully dissect how riboswitches function requires a better understanding of the energetic landscape for helical switching. Here, we report a sequencing-based high-throughput assay for monitoring in vitro transcription termination and use it to simultaneously characterize the functional effects of all 522 single point mutants of a glycine riboswitch type-1 singlet. Mutations throughout the riboswitch cause ligand-dependent defects, but only mutations within the terminator hairpin alter readthrough efficiencies in the absence of ligand. A comprehensive analysis of the expression platform reveals that ligand binding stabilizes the antiterminator by just 2-3 kcal/mol, indicating that the competing expression platform helices must be extremely close in energy to elicit a significant ligand-dependent response. These results demonstrate that gene regulation by this riboswitch is highly constrained by the energetics of ligand binding and conformational switching. These findings exemplify the energetic parameters of RNA conformational rearrangements driven by binding events.

Ribozyme-catalysed RNA synthesis using triplet building blocks

RNA-catalyzed RNA replication is widely believed to have supported a primordial biology. However, RNA catalysis is dependent upon RNA folding, and this yields structures that can block replication of such RNAs. To address this apparent paradox, we have re-examined the building blocks used for RNA replication. We report RNA-catalysed RNA synthesis on structured templates when using trinucleotide triphosphates (triplets) as substrates, catalysed by a general and accurate triplet polymerase ribozyme that emerged from in vitro evolution as a mutualistic RNA heterodimer. The triplets cooperatively invaded and unraveled even highly stable RNA secondary structures, and support non-canonical primer-free and bidirectional modes of RNA synthesis and replication. Triplet substrates thus resolve a central incongruity of RNA replication, and here allow the ribozyme to synthesise its own catalytic subunit ‘+’ and ‘–’ strands in segments and assemble them into a new active ribozyme.

Life as we know it relies on three types of molecules: DNA, which stores genetic information; proteins that carry out the chemical reactions necessary for life; and RNA, which relays information between the two. However, some scientists think that before life adopted DNA and proteins, it relied primarily on RNA.

Like DNA, strands of RNA contain genetic data. Yet, some RNA strands can also fold to form ribozymes, 3D structures that could have guided life’s chemical processes the way proteins do now. For early life to be built on RNA, though, this molecule must have had the ability to make copies of itself.

This duplication is a chemical reaction that could be driven by an ‘RNA replicase’ ribozyme. RNA strands are made of four different letters attached to each other in a specific order. When RNA is copied, one strand acts as a template, and a replicase ribozyme would accurately guide which letters are added to the strand under construction. However, no replicase ribozyme has been observed in existing life forms; this has led scientists to try to artificially create RNA replicase ribozymes that could copy themselves.

Until now, the best approaches have assumed that a replicase would add building blocks formed of a single letter one by one to grow a new strand. Yet, although ribozymes can be made to copy straight RNA templates this way, folded RNA templates – including the replicase ribozyme itself – impede copying. In this apparent paradox, a ribozyme needs to fold to copy RNA, but when folded, is itself copied poorly. Here, Attwater et al. wondered if choosing different building blocks might overcome this contradiction.

Biochemical techniques were used to engineer a ribozyme that copies RNA strands by adding letters not one-by-one, but three-by-three. Using three-letter ‘triplet’ building blocks, this new ribozyme can copy various folded RNA strands, including the active part of its own sequence.

This is because triplet building blocks have different, and sometimes unexpected, chemical properties compared to single-letter blocks. For example, these triplets work together to bind tightly to RNA strands and unravel structures that block RNA copying.

All life on Earth today uses a triplet RNA code to make proteins from DNA, and these experiments showed how RNA triplets might have helped RNA sustain early life forms. Further work is now needed to improve the ribozyme designed by Attwater et al. for efficient self-copying.

The structure of an E. coli tRNAfMet A1-U72 variant shows an unusual conformation of the A1-U72 base pair

Translation initiation in eukaryotes and archaea involves a methionylated initiator tRNA delivered to the ribosome in a ternary complex with e/aIF2 and GTP. Eukaryotic and archaeal initiator tRNAs contain a highly conserved A₁-U₇₂ base pair at the top of the acceptor stem. The importance of this base pair to discriminate initiator tRNAs from elongator tRNAs has been established previously using genetics and biochemistry. However, no structural data illustrating how the A₁-U₇₂ base pair participates in the accurate selection of the initiator tRNAs by the translation initiation systems are available. Here, we describe the crystal structure of a mutant E. coli initiator tRNA_f^MetA₁-U₇₂, aminoacylated with methionine, in which the C₁:A₇₂ mismatch at the end of the tRNA acceptor stem has been changed to an A₁-U₇₂ base pair. Sequence alignments show that the mutant E. coli tRNA is a good mimic of archaeal initiator tRNAs. The crystal structure, determined at 2.8 Å resolution, shows that the A₁-U₇₂ pair adopts an unusual arrangement. A₁ is in a syn conformation and forms a single H-bond interaction with U₇₂ This interaction requires protonation of the N1 atom of A₁ Moreover, the 5' phosphoryl group folds back into the major groove of the acceptor stem and interacts with the N7 atom of G₂ A possible role of this unusual geometry of the A₁-U₇₂ pair in the recognition of the initiator tRNA by its partners during eukaryotic and archaeal translation initiation is discussed.

Unusual Base-Pairing Interactions in Monomer-Template Complexes

Many high-resolution crystal structures have contributed to our understanding of the reaction pathway for catalysis by DNA and RNA polymerases, but the structural basis of nonenzymatic template-directed RNA replication has not been studied in comparable detail. Here we present crystallographic studies of the binding of ribonucleotide monomers to RNA primer-template complexes, with the goal of improving our understanding of the mechanism of nonenzymatic RNA copying, and of catalysis by polymerases. To explore how activated ribonucleotides recognize and bind to RNA templates, we synthesized an unreactive phosphonate-linked pyrazole analogue of guanosine 5'-phosphoro-2-methylimidazolide (2-MeImpG), a highly activated nucleotide that has been used extensively to study nonenzymatic primer extension. We cocrystallized this analogue with structurally rigidified RNA primer-template complexes carrying single or multiple monomer binding sites, and obtained high-resolution X-ray structures of these complexes. In addition to Watson-Crick base pairing, we repeatedly observed noncanonical guanine:cytidine base pairs in our crystal structures. In most structures, the phosphate and leaving group moieties of the monomers were highly disordered, while in others the distance from O3' of the primer to the phosphorus of the incoming monomer was too great to allow for reaction. We suggest that these effects significantly influence the rate and fidelity of nonenzymatic RNA replication, and that even primitive ribozyme polymerases could enhance RNA replication by enforcing Watson-Crick base pairing between monomers and primer-template complexes, and by bringing the reactive functional groups into closer proximity.

m(1)A and m(1)G disrupt A-RNA structure through the intrinsic instability of Hoogsteen base pairs

The B-DNA double helix can dynamically accommodate G-C and A-T base pairs in either Watson-Crick or Hoogsteen configurations. Here, we show that G-C(+) (in which + indicates protonation) and A-U Hoogsteen base pairs are strongly disfavored in A-RNA. As a result,N(1)-methyladenosine and N(1)-methylguanosine, which occur in DNA as a form of alkylation damage and in RNA as post-transcriptional modifications, have dramatically different consequences. Whereas they create G-C(+) and A-T Hoogsteen base pairs in duplex DNA, thereby maintaining the structural integrity of the double helix, they block base-pairing and induce local duplex melting in RNA. These observations provide a mechanism for disrupting RNA structure through post-transcriptional modifications. The different propensities to form Hoogsteen base pairs in B-DNA and A-RNA may help cells meet the opposing requirements of maintaining genome stability, on the one hand, and of dynamically modulating the structure of the epitranscriptome, on the other.

Kinetic and thermodynamic framework for P4-P6 RNA reveals tertiary motif modularity and modulation of the folding preferred pathway

The past decade has seen a wealth of 3D structural information about complex structured RNAs and identification of functional intermediates. Nevertheless, developing a complete and predictive understanding of the folding and function of these RNAs in biology will require connection of individual rate and equilibrium constants to structural changes that occur in individual folding steps and further relating these steps to the properties and behavior of isolated, simplified systems. To accomplish these goals we used the considerable structural knowledge of the folded, unfolded, and intermediate states of P4-P6 RNA. We enumerated structural states and possible folding transitions and determined rate and equilibrium constants for the transitions between these states using single-molecule FRET with a series of mutant P4-P6 variants. Comparisons with simplified constructs containing an isolated tertiary contact suggest that a given tertiary interaction has a stereotyped rate for breaking that may help identify structural transitions within complex RNAs and simplify the prediction of folding kinetics and thermodynamics for structured RNAs from their parts. The preferred folding pathway involves initial formation of the proximal tertiary contact. However, this preference was only ∼10 fold and could be reversed by a single point mutation, indicating that a model akin to a protein-folding contact order model will not suffice to describe RNA folding. Instead, our results suggest a strong analogy with a modified RNA diffusion-collision model in which tertiary elements within preformed secondary structures collide, with the success of these collisions dependent on whether the tertiary elements are in their rare binding-competent conformations.

Nuclease-Resistant c-di-AMP Derivatives That Differentially Recognize RNA and Protein Receptors

The ability of bacteria to sense environmental cues and adapt is essential for their survival. The use of second-messenger signaling molecules to translate these cues into a physiological response is a common mechanism employed by bacteria. The second messenger 3'-5'-cyclic diadenosine monophosphate (c-di-AMP) has been linked to a diverse set of biological processes involved in maintaining cell viability and homeostasis, as well as pathogenicity. A complex network of both protein and RNA receptors inside the cell activates specific pathways and mediates phenotypic outputs in response to c-di-AMP. Structural analysis of these RNA and protein receptors has revealed the different recognition elements employed by these effectors to bind the same small molecule. Herein, using a series of c-di-AMP analogues, we probed the interactions made with a riboswitch and a phosphodiesterase protein to identify the features important for c-di-AMP binding and recognition. We found that the ydaO riboswitch binds c-di-AMP in two discrete sites with near identical affinity and a Hill coefficient of 1.6. The ydaO riboswitch distinguishes between c-di-AMP and structurally related second messengers by discriminating against an amine at the C2 position more than a carbonyl at the C6 position. We also identified phosphate-modified analogues that bind both the ydaO RNA and GdpP protein with high affinity, whereas symmetrically modified ribose analogues exhibited a substantial decrease in ydaO affinity but retained high affinity for GdpP. These ligand modifications resulted in increased resistance to enzyme-catalyzed hydrolysis by the GdpP enzyme. Together, these data suggest that these c-di-AMP analogues could be useful as chemical tools to specifically target subsections of second-messenger signaling pathways.

RNA gymnastics in mammalian signal recognition particle assembly

More than one third of the cellular proteome is destined for incorporation into cell membranes or export from the cell. In all domains of life, the signal recognition particle (SRP) delivers these proteins to the membrane and protein traffic falls apart without SRP logistics. With the aid of a topogenic transport signal, SRP retrieves its cargo right at the ribosome, from where they are sorted to the translocation channel. Mammalian SRP is a ribonucleoprotein complex consisting of an SRP RNA of 300 nucleotides and 6 proteins bound to it. Assembly occurs in a hierarchical manner mainly in the nucleolus and only SRP54, which recognizes the signal sequence and regulates the targeting process, is added as the last component in the cytosol. Here we present an update on recent insights in the structure, function and dynamics of SRP RNA in SRP assembly with focus on the S domain, and present SRP as an example for the complex biogenesis of a rather small ribonucleoprotein particle.

Secondary structure encodes a cooperative tertiary folding funnel in the Azoarcus ribozyme

A requirement for specific RNA folding is that the free-energy landscape discriminate against non-native folds. While tertiary interactions are critical for stabilizing the native fold, they are relatively non-specific, suggesting additional mechanisms contribute to tertiary folding specificity. In this study, we use coarse-grained molecular dynamics simulations to explore how secondary structure shapes the tertiary free-energy landscape of the Azoarcus ribozyme. We show that steric and connectivity constraints posed by secondary structure strongly limit the accessible conformational space of the ribozyme, and that these so-called topological constraints in turn pose strong free-energy penalties on forming different tertiary contacts. Notably, native A-minor and base-triple interactions form with low conformational free energy, while non-native tetraloop/tetraloop–receptor interactions are penalized by high conformational free energies. Topological constraints also give rise to strong cooperativity between distal tertiary interactions, quantitatively matching prior experimental measurements. The specificity of the folding landscape is further enhanced as tertiary contacts place additional constraints on the conformational space, progressively funneling the molecule to the native state. These results indicate that secondary structure assists the ribozyme in navigating the otherwise rugged tertiary folding landscape, and further emphasize topological constraints as a key force in RNA folding.

Ribosome Structure Reveals Preservation of Active Sites in the Presence of a P-Site Wobble Mismatch

Translation initiation in the P site occasionally occurs at atypical (non-AUG) start codons, including those forming a mismatch in the third (wobble) position. During elongation, however, a pyrimidine-pyrimidine wobble mismatch may trigger a translation quality-control mechanism, whereby the P-site mismatch is thought to perturb the downstream A-site codon or the decoding center, thereby reducing translation fidelity and inducing termination of aberrant translation. We report a crystal structure of the 70S initiation complex containing an AUC codon in the ribosomal P site. Remarkably, the ribosome stabilizes the mismatched codon-anticodon helix, arranging a normally disruptive cytosine-cytosine pair into a Watson-Crick-like conformation. Translation-competent conformations of the tRNA, mRNA, and decoding center suggest that a P-site wobble-position mismatch in the 70S initiation complex does not pre-arrange the mRNA or decoding center to favor subsequent miscoding events.

Conserved residues in yeast initiator tRNA calibrate initiation accuracy by regulating preinitiation complex stability at the start codon

Eukaryotic initiator tRNA (tRNAi) contains several highly conserved, unique sequence features, yet their importance in accurate start codon selection is unknown. Using genetic and biochemical analyses, Dong et al. show that conserved bases throughout tRNAi, from the anticodon stem to the acceptor stem, play key roles in ensuring the fidelity of start codon recognition. This work delineates specific molecular functions for signature initiator tRNA residues and establishes their importance for initiation accuracy in living eukaryotic cells.

Eukaryotic initiator tRNA (tRNA_i) contains several highly conserved unique sequence features, but their importance in accurate start codon selection was unknown. Here we show that conserved bases throughout tRNA_i, from the anticodon stem to acceptor stem, play key roles in ensuring the fidelity of start codon recognition in yeast cells. Substituting the conserved G31:C39 base pair in the anticodon stem with different pairs reduces accuracy (the Sui⁻ [suppressor of initiation codon] phenotype), whereas eliminating base pairing increases accuracy (the Ssu⁻ [suppressor of Sui⁻] phenotype). The latter defect is fully suppressed by a Sui⁻ substitution of T-loop residue A54. These genetic data are paralleled by opposing effects of Sui⁻ and Ssu⁻ substitutions on the stability of methionylated tRNA_i (Met-tRNA_i) binding (in the ternary complex [TC] with eIF2-GTP) to reconstituted preinitiation complexes (PICs). Disrupting the C3:G70 base pair in the acceptor stem produces a Sui⁻ phenotype and also reduces the rate of TC binding to 40S subunits in vitro and in vivo. Both defects are suppressed by an Ssu⁻ substitution in eIF1A that stabilizes the open/P_OUT conformation of the PIC that exists prior to start codon recognition. Our data indicate that these signature sequences of tRNA_i regulate accuracy by distinct mechanisms, promoting the open/P_OUT conformation of the PIC (for C3:G70) or destabilizing the closed/P_IN state (for G31:C39 and A54) that is critical for start codon recognition.

MicroRNA miR-92a-1 biogenesis and mRNA targeting is modulated by a tertiary contact within the miR-17~92 microRNA cluster

While functional mature microRNAs (miRNAs) are small ∼22 base oligonucleotides that target specific mRNAs, miRNAs are initially expressed as long transcripts (pri-miRNAs) that undergo sequential processing to yield the mature miRNAs. We have previously reported that the pri-miR-17∼92 cluster adopts a compact globular folded structure that internalizes a 3' core domain resulting in reduced miRNA maturation and subsequent mRNA targeting. Using a site-specific photo-cross-linker we have identified a tertiary contact within the 3' core domain of the pri-miRNA between a non-miRNA stem-loop and the pre-miR-19b hairpin. This tertiary contact is involved in the formation of the compact globular fold of the cluster while its disruption enhances miR-92a expression and mRNA targeting. We propose that this tertiary contact serves as a molecular scaffold to restrict expression of the proposed antiangiogenic miR-92a, allowing for the overall pro-angiogenic effect of miR-17∼92 expression.

The scanning mechanism of eukaryotic translation initiation

In eukaryotes, the translation initiation codon is generally identified by the scanning mechanism, wherein every triplet in the messenger RNA leader is inspected for complementarity to the anticodon of methionyl initiator transfer RNA (Met-tRNAi). Binding of Met-tRNAi to the small (40S) ribosomal subunit, in a ternary complex (TC) with eIF2-GTP, is stimulated by eukaryotic initiation factor 1 (eIF1), eIF1A, eIF3, and eIF5, and the resulting preinitiation complex (PIC) joins the 5' end of mRNA preactivated by eIF4F and poly(A)-binding protein. RNA helicases remove secondary structures that impede ribosome attachment and subsequent scanning. Hydrolysis of eIF2-bound GTP is stimulated by eIF5 in the scanning PIC, but completion of the reaction is impeded at non-AUG triplets. Although eIF1 and eIF1A promote scanning, eIF1 and possibly the C-terminal tail of eIF1A must be displaced from the P decoding site to permit base-pairing between Met-tRNAi and the AUG codon, as well as to allow subsequent phosphate release from eIF2-GDP. A second GTPase, eIF5B, catalyzes the joining of the 60S subunit to produce an 80S initiation complex that is competent for elongation.

Steric complementarity in the decoding center is important for tRNA selection by the ribosome

Accurate tRNA selection by the ribosome is essential for the synthesis of functional proteins. Previous structural studies indicated that the ribosome distinguishes between cognate and near-cognate tRNAs by monitoring the geometry of the codon-anticodon helix in the decoding center using the universally conserved 16S ribosomal RNA bases G530, A1492 and A1493. These bases form hydrogen bonds with the 2'-hydroxyl groups of the codon-anticodon helix, which are expected to be disrupted with a near-cognate codon-anticodon helix. However, a recent structural study showed that G530, A1492 and A1493 form hydrogen bonds in a manner identical with that of both cognate and near-cognate codon-anticodon helices. To understand how the ribosome discriminates between cognate and near-cognate tRNAs, we made 2'-deoxynucleotide and 2'-fluoro substituted mRNAs, which disrupt the hydrogen bonds between the A site codon and G530, A1492 and A1493. Our results show that multiple 2'-deoxynucleotide substitutions in the mRNA substantially inhibit tRNA selection, whereas multiple 2'-fluoro substitutions in the mRNA have only modest effects on tRNA selection. Furthermore, the miscoding antibiotics paromomycin and streptomycin rescue the defects in tRNA selection with the multiple 2'-deoxynucleotide substituted mRNA. These results suggest that steric complementarity in the decoding center is more important than the hydrogen bonds between the A site codon and G530, A1492 and A1493 for tRNA selection.

Cooperative tertiary interaction network guides RNA folding

Noncoding RNAs form unique 3D structures, which perform many regulatory functions. To understand how RNAs fold uniquely despite a small number of tertiary interaction motifs, we mutated the major tertiary interactions in a group I ribozyme by single-base substitutions. The resulting perturbations to the folding energy landscape were measured using SAXS, ribozyme activity, hydroxyl radical footprinting, and native PAGE. Double- and triple-mutant cycles show that most tertiary interactions have a small effect on the stability of the native state. Instead, the formation of core and peripheral structural motifs is cooperatively linked in near-native folding intermediates, and this cooperativity depends on the native helix orientation. The emergence of a cooperative interaction network at an early stage of folding suppresses nonnative structures and guides the search for the native state. We suggest that cooperativity in noncoding RNAs arose from natural selection of architectures conducive to forming a unique, stable fold.

Molecular mechanism of scanning and start codon selection in eukaryotes

The correct translation of mRNA depends critically on the ability to initiate at the right AUG codon. For most mRNAs in eukaryotic cells, this is accomplished by the scanning mechanism, wherein the small (40S) ribosomal subunit attaches to the 5' end of the mRNA and then inspects the leader base by base for an AUG in a suitable context, using complementarity with the anticodon of methionyl initiator tRNA (Met-tRNAiMet) as the key means of identifying AUG. Over the past decade, a combination of yeast genetics, biochemical analysis in reconstituted systems, and structural biology has enabled great progress in deciphering the mechanism of ribosomal scanning. A robust molecular model now exists, describing the roles of initiation factors, notably eukaryotic initiation factor 1 (eIF1) and eIF1A, in stabilizing an "open" conformation of the 40S subunit with Met-tRNAiMet bound in a low-affinity state conducive to scanning and in triggering rearrangement into a "closed" conformation incompatible with scanning, which features Met-tRNAiMet more tightly bound to the "P" site and base paired with AUG. It has also emerged that multiple DEAD-box RNA helicases participate in producing a single-stranded "landing pad" for the 40S subunit and in removing the secondary structure to enable the mRNA to traverse the 40S mRNA-binding channel in the single-stranded form for base-by-base inspection in the P site.

2'-Fluoro substituents can mimic native 2'-hydroxyls within structured RNA

The ability of fluorine in a C-F bond to act as a hydrogen bond acceptor is controversial. To test such ability in complex RNA macromolecules, we have replaced native 2'-OH groups with 2'-F and 2'-H groups in two related systems, the Tetrahymena group I ribozyme and the ΔC209 P4-P6 RNA domain. In three cases the introduced 2'-F mimics the native 2'-OH group, suggesting that the fluorine atom can accept a hydrogen bond. In each of these cases the native hydroxyl group interacts with a purine exocyclic amine. Our results give insight about the properties of organofluorine and suggest a possible general biochemical signature for tertiary interactions between 2'-hydroxyl groups and exocyclic amino groups within RNA.

Removal of covalent heterogeneity reveals simple folding behavior for P4-P6 RNA

RNA folding landscapes have been described alternately as simple and as complex. The limited diversity of RNA residues and the ability of RNA to form stable secondary structures prior to adoption of a tertiary structure would appear to simplify folding relative to proteins. Nevertheless, there is considerable evidence for long-lived misfolded RNA states, and these observations have suggested rugged energy landscapes. Recently, single molecule fluorescence resonance energy transfer (smFRET) studies have exposed heterogeneity in many RNAs, consistent with deeply furrowed rugged landscapes. We turned to an RNA of intermediate complexity, the P4-P6 domain from the Tetrahymena group I intron, to address basic questions in RNA folding. P4-P6 exhibited long-lived heterogeneity in smFRET experiments, but the inability to observe exchange in the behavior of individual molecules led us to probe whether there was a non-conformational origin to this heterogeneity. We determined that routine protocols in RNA preparation and purification, including UV shadowing and heat annealing, cause covalent modifications that alter folding behavior. By taking measures to avoid these treatments and by purifying away damaged P4-P6 molecules, we obtained a population of P4-P6 that gave near-uniform behavior in single molecule studies. Thus, the folding landscape of P4-P6 lacks multiple deep furrows that would trap different P4-P6 molecules in different conformations and contrasts with the molecular heterogeneity that has been seen in many smFRET studies of structured RNAs. The simplicity of P4-P6 allowed us to reliably determine the thermodynamic and kinetic effects of metal ions on folding and to now begin to build more detailed models for RNA folding behavior.

Predicting and modeling RNA architecture

A general approach for modeling the architecture of large and structured RNA molecules is described. The method exploits the modularity and the hierarchical folding of RNA architecture that is viewed as the assembly of preformed double-stranded helices defined by Watson-Crick base pairs and RNA modules maintained by non-Watson-Crick base pairs. Despite the extensive molecular neutrality observed in RNA structures, specificity in RNA folding is achieved through global constraints like lengths of helices, coaxiality of helical stacks, and structures adopted at the junctions of helices. The Assemble integrated suite of computer tools allows for sequence and structure analysis as well as interactive modeling by homology or ab initio assembly with possibilities for fitting within electronic density maps. The local key role of non-Watson-Crick pairs guides RNA architecture formation and offers metrics for assessing the accuracy of three-dimensional models in a more useful way than usual root mean square deviation (RMSD) values.

Use of RNA tertiary interaction modules for the crystallisation of the spliceosomal snRNP core domain

RNA is known to perform diverse roles in the cell, often as ribonucleoprotein (RNP) particles. While the crystal structure of these RNP particles could provide crucial insights into their functions, crystallographic work on RNP complexes is often hampered by difficulties in obtaining well-diffracting crystals. The small nuclear ribonucleoprotein (snRNP) core domain, acting as an assembly nucleus for the maturation of snRNPs, plays a crucial role in the biogenesis of four of the spliceosomal snRNPs. We have succeeded in crystallising the human U4 snRNP core domain containing seven Sm proteins and a truncated U4 snRNA variant. The most critical factor in our success in the crystallisation was the introduction of various tertiary interaction modules into the RNA that could promote crystal packing without altering the core structure. Here, we describe various strategies employed in our crystallisation effort that could be applied to crystallisation of other RNP particles.

Correlation of RNA secondary structure statistics with thermodynamic stability and applications to folding

The accurate prediction of the secondary and tertiary structure of an RNA with different folding algorithms is dependent on several factors, including the energy functions. However, an RNA higher-order structure cannot be predicted accurately from its sequence based on a limited set of energy parameters. The inter- and intramolecular forces between this RNA and other small molecules and macromolecules, in addition to other factors in the cell such as pH, ionic strength, and temperature, influence the complex dynamics associated with transition of a single stranded RNA to its secondary and tertiary structure. Since all of the factors that affect the formation of an RNAs 3D structure cannot be determined experimentally, statistically derived potential energy has been used in the prediction of protein structure. In the current work, we evaluate the statistical free energy of various secondary structure motifs, including base-pair stacks, hairpin loops, and internal loops, using their statistical frequency obtained from the comparative analysis of more than 50,000 RNA sequences stored in the RNA Comparative Analysis Database (rCAD) at the Comparative RNA Web (CRW) Site. Statistical energy was computed from the structural statistics for several datasets. While the statistical energy for a base-pair stack correlates with experimentally derived free energy values, suggesting a Boltzmann-like distribution, variation is observed between different molecules and their location on the phylogenetic tree of life. Our statistical energy values calculated for several structural elements were utilized in the Mfold RNA-folding algorithm. The combined statistical energy values for base-pair stacks, hairpins and internal loop flanks result in a significant improvement in the accuracy of secondary structure prediction; the hairpin flanks contribute the most.

The dimeric proto-ribosome: Structural details and possible implications on the origin of life

A symmetric pocket-like entity, composed of two L-shaped RNA units, encircles the peptide synthesis site within the contemporary ribosome. This entity was suggested to be the vestige of a dimeric proto-ribosome, which could have formed spontaneously in the prebiotic world, catalyzing non-coded peptide bond formation and elongation. This structural element, beyond offering the initial step in the evolution of translation, is hypothesized here to be linked to the origin of life. By catalyzing the production of random peptide chains, the proto-ribosome could have enabled the formation of primary enzymes, launching a process of co-evolution of the translation apparatus and the proteins, thus presenting an alternative to the RNA world hypothesis.

Comparing the three-dimensional structures of Dicistroviridae IGR IRES RNAs with other viral RNA structures

The intergenic region (IGR) internal ribosome entry site (IRES) RNAs do not require any of the canonical translation initiation factors to recruit the ribosome to the viral RNA, they eliminate the need for initiator tRNA, and they begin translation from the A-site. The function of these IRESs depends on a specific three-dimensional folded RNA structure. Thus, a complete understanding of the mechanisms of action of these IRESs requires that we understand their structure in detail. Recently, the structures of both domains of the IGR IRES RNAs were solved by X-ray crystallography, providing the first glimpse into an entire IRES RNA structure. Here, I present an analysis of these structures, emphasizing how the structures explain many aspects of IGR IRES function, discussing how these structures have similarities to motifs found in other viral RNAs, and illustrating how these structures give rise to new mechanistic hypotheses.

Annotation of tertiary interactions in RNA structures reveals variations and correlations

RNA tertiary motifs play an important role in RNA folding and biochemical functions. To help interpret the complex organization of RNA tertiary interactions, we comprehensively analyze a data set of 54 high-resolution RNA crystal structures for motif occurrence and correlations. Specifically, we search seven recognized categories of RNA tertiary motifs (coaxial helix, A-minor, ribose zipper, pseudoknot, kissing hairpin, tRNA D-loop/T-loop, and tetraloop-tetraloop receptor) by various computer programs. For the nonredundant RNA data set, we find 613 RNA tertiary interactions, most of which occur in the 16S and 23S rRNAs. An analysis of these motifs reveals the diversity and variety of A-minor motif interactions and the various possible loop-loop receptor interactions that expand upon the tetraloop-tetraloop receptor. Correlations between motifs, such as pseudoknot or coaxial helix with A-minor, reveal higher-order patterns. These findings may ultimately help define tertiary structure restraints for RNA tertiary structure prediction. A complete annotation of the RNA diagrams for our data set is available at http://www.biomath.nyu.edu/motifs/.

Navigating the RNA folding landscape

Analysis of individual RNA folding reactions reveals that, as in proteins, cooperative interactions selectively drive RNA toward its biologically active, native conformation. This new work establishes a platform for future investigations of the physical principles underlying the assembly of large RNA enzymes.

Direct measurement of tertiary contact cooperativity in RNA folding

All structured biological macromolecules must overcome the thermodynamic folding problem to populate a unique functional state among a vast ensemble of unfolded and alternate conformations. The exploration of cooperativity in protein folding has helped reveal and distinguish the underlying mechanistic solutions to this folding problem. Analogous dissections of RNA tertiary stability remain elusive, however, despite the central biological importance of folded RNA molecules and the potential to reveal fundamental properties of structured macromolecules via comparisons of protein and RNA folding. We report a direct quantitative measure of tertiary contact cooperativity in a folded RNA. We precisely measured the stability of an independently folding P4-P6 domain from the Tetrahymena thermophila group I intron by single molecule fluorescence resonance energy transfer (smFRET). Using wild-type and mutant RNAs, we found that cooperativity between the two tertiary contacts enhances P4-P6 stability by 3.2 +/- 0.2 kcal/mol.

Comprehensive features of natural and in vitro selected GNRA tetraloop-binding receptors

Specific recognitions of GNRA tetraloops by small helical receptors are among the most widespread long-range packing interactions in large ribozymes. However, in contrast to GYRA and GAAA tetraloops, very few GNRA/receptor interactions have yet been identified to involve GGAA tetraloops in nature. A novel in vitro selection scheme based on a rigid self-assembling tectoRNA scaffold designed for isolation of intermolecular interactions with A-minor motifs has yielded new GGAA tetraloop-binding receptors with affinity in the nanomolar range. One of the selected receptors is a novel 12 nt RNA motif, (CCUGUG ... AUCUGG), that recognizes GGAA tetraloop hairpin with a remarkable specificity and affinity. Its physical and chemical characteristics are comparable to those of the well-studied '11nt' GAAA tetraloop receptor motif. A second less specific motif (CCCAGCCC ... GAUAGGG) binds GGRA tetraloops and appears to be related to group IC3 tetraloop receptors. Mutational, thermodynamic and comparative structural analysis suggests that natural and in vitro selected GNRA receptors can essentially be grouped in two major classes of GNRA binders. New insights about the evolution, recognition and structural modularity of GNRA and A-minor RNA-RNA interactions are proposed.

Structure of the SAM-II riboswitch bound to S-adenosylmethionine

In bacteria, numerous genes harbor regulatory elements in the 5' untranslated regions of their mRNA, termed riboswitches, which control gene expression by binding small-molecule metabolites. These sequences influence the secondary and tertiary structure of the RNA in a ligand-dependent manner, thereby directing its transcription or translation. The crystal structure of an S-adenosylmethionine-responsive riboswitch found predominantly in proteobacteria, SAM-II, has been solved to reveal a second means by which RNA interacts with this important cellular metabolite. Notably, this is the first structure of a complete riboswitch containing all sequences associated with both the ligand binding aptamer domain and the regulatory expression platform. Chemical probing of this RNA in the absence and presence of ligand shows how the structure changes in response to S-adenosylmethionine to sequester the ribosomal binding site and affect translational gene regulation.

Chemical basis of glycine riboswitch cooperativity

The glycine binding riboswitch forms a unique tandem aptamer structure that binds glycine cooperatively. We employed nucleotide analog interference mapping (NAIM) and mutagenesis to explore the chemical basis of glycine riboswitch cooperativity. Based on the interference pattern, at least two sites appear to facilitate cooperative tertiary interactions, namely, the minor groove of the P1 helix from aptamer 1 and the major groove of the P3a helix from both aptamers. Mutation of these residues altered both the cooperativity and binding affinity of the riboswitch. The data support a model in which the P1 helix of the first aptamer participates in a tertiary interaction important for cooperativity, while nucleotides in the P1 helix of the second aptamer interface with the expression platform. These data have direct analogy to well-characterized mutations in hemoglobin, which provides a framework for considering cooperativity in this RNA-based system.

Kinetic redistribution of native and misfolded RNAs by a DEAD-box chaperone

DExD/H-box proteins are ubiquitously involved in RNA-mediated processes and use ATP to accelerate conformational changes in RNA. However, their mechanisms of action, and what determines which RNA species are targeted, are not well understood. Here we show that the DExD/H-box protein CYT-19, a general RNA chaperone, mediates ATP-dependent unfolding of both the native conformation and a long-lived misfolded conformation of a group I catalytic RNA with efficiencies that depend on the stabilities of the RNA species but not on specific structural features. CYT-19 then allows the RNA to refold, changing the distribution from equilibrium to kinetic control. Because misfolding is favoured kinetically, conditions that allow unfolding of the native RNA yield large increases in the population of misfolded species. Our results suggest that DExD/H-box proteins act with sufficient breadth and efficiency to allow structured RNAs to populate a wider range of conformations than would be present at equilibrium. Thus, RNAs may face selective pressure to stabilize their active conformations relative to inactive ones to avoid significant redistribution by DExD/H-box proteins. Conversely, RNAs whose functions depend on forming multiple conformations may rely on DExD/H-box proteins to increase the populations of less stable conformations, thereby increasing their overall efficiencies.

Low specificity of metal ion binding in the metal ion core of a folded RNA

The structure and activity of nucleic acids depend on their interactions with metal ions. Fundamental to these interactions is the degree of specificity observed between the metal ions and nucleic acids, and a complete description of nucleic acid folding requires that we understand the nature of the interactions with metal ions, including specificity. The prior demonstration that high concentrations of monovalent cations prevent nonspecific association of divalent ions with nucleic acids provides a novel and powerful means to examine site-specific metal ion binding isolated from complicating effects of the ion atmosphere. Using these high monovalent cation solution conditions we have monitored the affinity of a series of divalent metal ions for two site-specific metal ion binding sites in the P4-P6 domain of the Tetrahymena group I intron ribozyme. The metal ion core of this highly structured RNA binds two divalent metal ions under these conditions. Despite multiple metal ion-RNA interactions observed in the X-ray crystallographic structure of P4-P6 RNA at the metal ion binding sites, these sites exhibit low specificity among Mn(2+), Mg(2+), Ca(2+), Ni(2+), and Zn(2+). Nevertheless, the largest divalent metal ions tested, Sr(2+) and Ba(2+), were excluded from binding, exhibiting affinities at least two orders of magnitude weaker than observed for the other metal ions. Thus, a picture emerges of two metal ion binding sites, each with a high tolerance for metal ions with different properties but also with limits to accommodation. These limits presumably arise from steric or electrostatic features of the metal ion binding sites.

Leading RNA tertiary interactions: structures, energies, and water insertion of A-minor and P-interactions. A quantum chemical view

Complex molecular shapes of ribosomal RNA molecules are stabilized by recurrent types of tertiary interactions involving highly specific and conserved non-Watson-Crick base pairs, triplets, and quartets. We analyzed the intrinsic structure and stability of the P-motif and the four basic types of A-minor interactions (types I, II, III, and 0), which represent the most prominent RNA tertiary interaction patterns refined in the course of evolution. In the studied interactions, the electron correlation component of the stabilization usually exceeds the Hartree-Fock (HF) term, leading to a strikingly different balance of forces as compared to standard base pairing stabilized primarily by the HF term. In other words, the A-minor and P-interactions are considerably more influenced by the dispersion energy as compared to canonical base pairs, which makes them particularly suitable to zip the folded RNA structures that are substantially hydrated even in their interior. Continuum solvent COSMO calculations confirm that the stability of the canonical GC base pair is affected (reduced) by the continuum solvent screening considerably more than the stability of the A-minor interaction. Among the studied systems, the strong A-minor II and weak A-minor III interactions require water molecules to stabilize the experimental geometry. Gas-phase optimization of the canonical A-minor II A/CG triplet without water results in a geometry that is clearly inconsistent with the RNA structure. The gas-phase structure of the P-interaction and the most stable A-minor I interaction nicely agrees with the geometries occurring in the ribosome. A-minor I can also adopt an alternative water-mediated substate rather often observed in X-ray and molecular dynamics studies. The A-minor I water bridge, however, does not appear to stabilize the tertiary contact, and its role is to provide structural flexibility to this binding pattern within the context of the RNA structure. Interestingly, the insertion of a polar water molecule in the A-minor I A/CG tertiary contact occurring in the A/C tertiary pair is stabilized primarily by the HF (electrostatic) interaction energy, while the dispersion-controlled A/G contact remains firmly bound. Thus, the intrinsic balance of forces as revealed by quantum mechanics (QM) calculations nicely correlates with many behavioral aspects of the studied interactions inside RNA. The comparison of interaction energies computed using quantum chemistry and an AMBER force field reveals that common molecular mechanics calculations perform rather well, except that the strength of the P-interaction is modestly overestimated. We also briefly discuss the non-negligible methodological differences when evaluating simple base-base nucleic acids base pairs and the complex RNA tertiary base pairing patterns using QM procedures.

Exploratory studies to investigate a linked prebiotic origin of RNA and coded peptides. 3rd communication. Behaviour of 5-amino-1H-imidazole-4-carbonitrile derivatives

The potentially prebiotic synthesis of pyrimidine ribonucleotides by stepwise nucleobase assembly on arabinose-3-phosphate derivatives has been demonstrated in previous work. The generation of xylose-2-phosphate derivatives by aldolisation, and the behaviour of these compounds under the conditions of pyrimidine nucleobase assembly have also been described. In this paper, the scope for generation of purine nucleotides via 3,3'-anhydro-xylo-nucleotides is investigated. In neutral D2O solution, the potential intermediate 47 (Schemes 6 and 8) undergoes H-C2 --> D-C2 exchange, but no appreciable reaction with cyanide or cyanamide occurs. The exchange chemistry expands options for purine nucleobase assembly on sugar phosphate scaffolds.