Structural basis of the enhanced stability of a mutant ribozyme domain and a detailed view of RNA--solvent interactions

NucleoFind: a deep-learning network for interpreting nucleic acid electron density

Nucleic acid electron density interpretation after phasing by molecular replacement or other methods remains a difficult problem for computer programs to deal with. Programs tend to rely on time-consuming and computationally exhaustive searches to recognise characteristic features. We present NucleoFind, a deep-learning-based approach to interpreting and segmenting electron density. Using an electron density map from X-ray crystallography obtained after molecular replacement, the positions of the phosphate group, sugar ring and nitrogenous base group can be predicted with high accuracy. On average, 78% of phosphate atoms, 85% of sugar atoms and 83% of base atoms are positioned in predicted density after giving NucleoFind maps produced following successful molecular replacement. NucleoFind can use the wealth of context these predicted maps provide to build more accurate and complete nucleic acid models automatically.

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.

The structural features of the ligand-free moaA riboswitch and its ion-dependent folding

Riboswitches are structural elements of mRNA involved in the regulation of gene expression by responding to specific cellular metabolites. To fulfil their regulatory function, riboswitches prefold into an active state, the so-called binding competent form, that guarantees metabolite binding and allows a consecutive refolding of the RNA. Here, we describe the folding pathway to the binding competent form as well as the ligand free structure of the moaA riboswitch of E. coli. This RNA proposedly responds to the molybdenum cofactor (Moco), a highly oxygen-sensitive metabolite, essential in the carbon and sulfur cycles of eukaryotes. K⁺- and Mg²⁺-dependent footprinting assays and spectroscopic investigations show a high degree of structure formation of this RNA already at very low ion-concentrations. Mg²⁺ facilitates additionally a general compaction of the riboswitch towards its proposed active structure. We show that this fold agrees with the earlier suggested secondary structure which included also a long-range tetraloop/tetraloop-receptor like interaction. Metal ion cleavage assays revealed specific Mg²⁺-binding pockets within the moaA riboswitch. These Mg²⁺ binding pockets are good indicators for the potential Moco binding site, since in riboswitches, Mg²⁺ was shown to be necessary to bind phosphate-carrying metabolites. The importance of the phosphate and of other functional groups of Moco is highlighted by binding assays with tetrahydrobiopterin, the reduced and oxygen-sensitive core moiety of Moco. We demonstrate that the general molecular shape of pterin by its own is insufficient for the recognition by the riboswitch.

Structure of the OMEGA nickase IsrB in complex with ωRNA and target DNA

RNA-guided systems, such as CRISPR-Cas, combine programmable substrate recognition with enzymatic function, a combination that has been used advantageously to develop powerful molecular technologies1,2. Structural studies of these systems have illuminated how the RNA and protein jointly recognize and cleave their substrates, guiding rational engineering for further technology development3. Recent work identified a new class of RNA-guided systems, termed OMEGA, which include IscB, the likely ancestor of Cas9, and the nickase IsrB, a homologue of IscB lacking the HNH nuclease domain4. IsrB consists of only around 350 amino acids, but its small size is counterbalanced by a relatively large RNA guide (roughly 300-nt ωRNA). Here, we report the cryogenic-electron microscopy structure of Desulfovirgula thermocuniculi IsrB (DtIsrB) in complex with its cognate ωRNA and a target DNA. We find the overall structure of the IsrB protein shares a common scaffold with Cas9. In contrast to Cas9, however, which uses a recognition (REC) lobe to facilitate target selection, IsrB relies on its ωRNA, part of which forms an intricate ternary structure positioned analogously to REC. Structural analyses of IsrB and its ωRNA as well as comparisons to other RNA-guided systems highlight the functional interplay between protein and RNA, advancing our understanding of the biology and evolution of these diverse systems.

Sub-3-Å cryo-EM structure of RNA enabled by engineered homomeric self-assembly

High-resolution structural studies are essential for understanding the folding and function of diverse RNAs. Herein, we present a nanoarchitectural engineering strategy for efficient structural determination of RNA-only structures using single-particle cryogenic electron microscopy (cryo-EM). This strategy-ROCK (RNA oligomerization-enabled cryo-EM via installing kissing loops)-involves installing kissing-loop sequences onto the functionally nonessential stems of RNAs for homomeric self-assembly into closed rings with multiplied molecular weights and mitigated structural flexibility. ROCK enables cryo-EM reconstruction of the Tetrahymena group I intron at 2.98-Å resolution overall (2.85 Å for the core), allowing de novo model building of the complete RNA, including the previously unknown peripheral domains. ROCK is further applied to two smaller RNAs-the Azoarcus group I intron and the FMN riboswitch, revealing the conformational change of the former and the bound ligand in the latter. ROCK holds promise to greatly facilitate the use of cryo-EM in RNA structural studies.

Towards the automatic crystal structure solution of nucleic acids: automated model building using the new CAB program

CAB, a recently described automated model-building (AMB) program, has been modified to work effectively with nucleic acids. To this end, several new algorithms have been introduced and the libraries have been updated. To reduce the input average phase error, ligand heavy atoms are now located before starting the CAB interpretation of the electron-density maps. Furthermore, alternative approaches are used depending on whether the ligands belong to the target or to the model chain used in the molecular-replacement step. Robust criteria are then applied to decide whether the AMB model is acceptable or whether it must be modified to fit prior information on the target structure. In the latter case, the model chains are rearranged to fit prior information on the target chains. Here, the performance of the new AMB program CAB applied to various nucleic acid structures is discussed. Other well documented programs such as Nautilus, ARP/wARP and phenix.autobuild were also applied and the experimental results are described.

Cryo-EM structures of full-length Tetrahymena ribozyme at 3.1 Å resolution

Single-particle cryogenic electron microscopy (cryo-EM) has become a standard technique for determining protein structures at atomic resolution^1-3. However, cryo-EM studies of protein-free RNA are in their early days. The Tetrahymena thermophila group I self-splicing intron was the first ribozyme to be discovered and has been a prominent model system for the study of RNA catalysis and structure-function relationships⁴, but its full structure remains unknown. Here we report cryo-EM structures of the full-length Tetrahymena ribozyme in substrate-free and bound states at a resolution of 3.1 Å. Newly resolved peripheral regions form two coaxially stacked helices; these are interconnected by two kissing loop pseudoknots that wrap around the catalytic core and include two previously unforeseen (to our knowledge) tertiary interactions. The global architecture is nearly identical in both states; only the internal guide sequence and guanosine binding site undergo a large conformational change and a localized shift, respectively, upon binding of RNA substrates. These results provide a long-sought structural view of a paradigmatic RNA enzyme and signal a new era for the cryo-EM-based study of structure-function relationships in ribozymes.

The Halastavi árva Virus Intergenic Region IRES Promotes Translation by the Simplest Possible Initiation Mechanism

Dicistrovirus intergenic region internal ribosomal entry sites (IGR IRESs) do not require initiator tRNA, an AUG codon, or initiation factors and jumpstart translation from the middle of the elongation cycle via formation of IRES/80S complexes resembling the pre-translocation state. eEF2 then translocates the [codon-anticodon]-mimicking pseudoknot I (PKI) from ribosomal A sites to P sites, bringing the first sense codon into the decoding center. Halastavi árva virus (HalV) contains an IGR that is related to previously described IGR IRESs but lacks domain 2, which enables these IRESs to bind to individual 40S ribosomal subunits. By using in vitro reconstitution and cryoelectron microscopy (cryo-EM), we now report that the HalV IGR IRES functions by the simplest initiation mechanism that involves binding to 80S ribosomes such that PKI is placed in the P site, so that the A site contains the first codon that is directly accessible for decoding without prior eEF2-mediated translocation of PKI.

Hexahydrated Mg2+ Binding and Outer-Shell Dehydration on RNA Surface

The interaction between metal ions, especially Mg²⁺ ions, and RNA plays a critical role in RNA folding. Upon binding to RNA, a metal ion that is fully hydrated in bulk solvent can become dehydrated. Here we use molecular dynamics simulation to investigate the dehydration of bound hexahydrated Mg²⁺ ions. We find that a hydrated Mg²⁺ ion in the RNA groove region can involve significant dehydration in the outer hydration shell. The first or innermost hydration shell of the Mg²⁺ ion, however, is retained during the simulation because of the strong ion-water electrostatic attraction. As a result, water-mediated hydrogen bonding remains an important form for Mg²⁺-RNA interaction. Analysis for ions at different binding sites shows that the most pronounced water deficiency relative to the fully hydrated state occurs at a radial distance of around 11 Å from the center of the ion. Based on the independent 200 ns molecular dynamics simulations for three different RNA structures (Protein Data Bank: 1TRA, 2TPK, and 437D), we find that Mg²⁺ ions overwhelmingly dominate over monovalent ions such as Na⁺ and K⁺ in ion-RNA binding. Furthermore, application of the free energy perturbation method leads to a quantitative relationship between the Mg²⁺ dehydration free energy and the local structural environment. We find that ΔΔG_hyd, the change of the Mg²⁺ hydration free energy upon binding to RNA, varies linearly with the inverse distance between the Mg²⁺ ion and the nearby nonbridging oxygen atoms of the phosphate groups, and ΔΔG_hyd can reach -2.0 kcal/mol and -3.0 kcal/mol for an Mg²⁺ ion bound to the surface and to the groove interior, respectively. In addition, the computation results in an analytical formula for the hydration ratio as a function of the average inverse Mg²⁺-O distance. The results here might be useful for further quantitative investigations of ion-RNA interactions in RNA folding.

Nucleobase carbonyl groups are poor Mg2+ inner-sphere binders but excellent monovalent ion binders-a critical PDB survey

Precise knowledge of Mg²⁺ inner-sphere binding site properties is vital for understanding the structure and function of nucleic acid systems. Unfortunately, the PDB, which represents the main source of Mg²⁺ binding sites, contains a substantial number of assignment issues that blur our understanding of the functions of these ions. Here, following a previous study devoted to Mg²⁺ binding to nucleobase nitrogens, we surveyed nucleic acid X-ray structures from the PDB with resolutions ≤2.9 Å to classify the Mg²⁺ inner-sphere binding patterns to nucleotide carbonyl, ribose hydroxyl, cyclic ether, and phosphodiester oxygen atoms. From this classification, we derived a set of "prior-knowledge" nucleobase Mg²⁺ binding sites. We report that crystallographic examples of trustworthy nucleobase Mg²⁺ binding sites are fewer than expected since many of those are associated with misidentified Na⁺ or K⁺ We also emphasize that binding of Na⁺ and K⁺ to nucleic acids is much more frequent than anticipated. Overall, we provide evidence derived from X-ray structures that nucleobases are poor inner-sphere binders for Mg²⁺ but good binders for monovalent ions. Based on strict stereochemical criteria, we propose an extended set of guidelines designed to help in the assignment and validation of ions directly contacting nucleobase and ribose atoms. These guidelines should help in the interpretation of X-ray and cryo-EM solvent density maps. When borderline Mg²⁺ stereochemistry is observed, alternative placement of Na⁺, K⁺, or Ca²⁺ must be considered. We also critically examine the use of lanthanides (Yb³⁺, Tb³⁺) as Mg²⁺ substitutes in crystallography experiments.

Revealing the distinct folding phases of an RNA three-helix junction

Remarkable new insight has emerged into the biological role of RNA in cells. RNA folding and dynamics enable many of these newly discovered functions, calling for an understanding of RNA self-assembly and conformational dynamics. Because RNAs pass through multiple structures as they fold, an ensemble perspective is required to visualize the flow through fleetingly populated sets of states. Here, we combine microfluidic mixing technology and small angle X-ray scattering (SAXS) to measure the Mg-induced folding of a small RNA domain, the tP5abc three helix junction. Our measurements are interpreted using ensemble optimization to select atomically detailed structures that recapitulate each experimental curve. Structural ensembles, derived at key stages in both time-resolved studies and equilibrium titrations, reproduce the features of known intermediates, and more importantly, offer a powerful new structural perspective on the time-progression of folding. Distinct collapse phases along the pathway appear to be orchestrated by specific interactions with Mg ions. These key interactions subsequently direct motions of the backbone that position the partners of tertiary contacts for later bonding, and demonstrate a remarkable synergy between Mg and RNA across numerous time-scales.

In Crystallo Selection to Establish New RNA Crystal Contacts

Crystallography is a major technique for determining large RNA structures. Obtaining diffraction-quality crystals has been the bottleneck. Although several RNA crystallization methods have been developed, the field strongly needs additional approaches. Here we invented an in crystallo selection strategy for identifying mutations that enhance a target RNA's crystallizability. The strategy includes constructing an RNA pool containing random mutations, obtaining crystals, and amplifying the sequences enriched by crystallization. We demonstrated a proof-of-principle application to the P4-P6 domain from the Tetrahymena ribozyme. We further determined the structures of four selected mutants. All four establish new crystal lattice contacts while maintaining the native structure. Three mutants achieve this by relocating bulges and one by making a helix more flexible. In crystallo selection provides opportunities to improve crystals of RNAs or RNA-ligand complexes. Our results also suggest that mutants may be rationally designed for crystallization by "walking" a bulge along the RNA chain.

Mg2+ ions: do they bind to nucleobase nitrogens?

Given the many roles proposed for Mg2+ in nucleic acids, it is essential to accurately determine their binding modes. Here, we surveyed the PDB to classify Mg2+ inner-sphere binding patterns to nucleobase imine N1/N3/N7 atoms. Among those, purine N7 atoms are considered to be the best nucleobase binding sites for divalent metals. Further, Mg2+ coordination to N7 has been implied in several ribozyme catalytic mechanisms. We report that Mg2+ assigned near imine nitrogens derive mostly from poor interpretations of electron density patterns and are most often misidentified Na+, K+, NH4+ ions, water molecules or spurious density peaks. Consequently, apart from few documented exceptions, Mg2+ ions do not bind to N7 atoms. Without much of a surprise, Mn2+, Zn2+ and Cd2+, which have a higher affinity for nitrogens, may contact N7 atoms when present in crystallization buffers. In this respect, we describe for the first time a potential Zn2+ ribosomal binding site involving two purine N7 atoms. Further, we provide a set of guidelines to help in the assignment of Mg2+ in crystallographic, cryo-EM, NMR and model building practices and discuss implications of our findings related to ion substitution experiments.

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.

Assessing the Potential Effects of Active Site Mg2+ Ions in the glmS Ribozyme-Cofactor Complex

Ribozymes employ diverse catalytic strategies in their self-cleavage mechanisms, including the use of divalent metal ions. This work explores the effects of Mg2+ ions in the active site of the glmS ribozyme-GlcN6P cofactor complex using computational methods. Deleterious and potentially beneficial effects of an active site Mg2+ ion on the self-cleavage reaction were identified. The presence of a Mg2+ ion near the scissile phosphate oxygen atoms at the cleavage site was determined to be deleterious, and thereby anticatalytic, due to electrostatic repulsion of the cofactor, disruption of key hydrogen-bonding interactions, and obstruction of nucleophilic attack. On the other hand, the presence of a Mg2+ ion at another position in the active site, the Hoogsteen face of the putative base, was found to avoid these deleterious effects and to be potentially catalytically favorable owing to the stabilization of negative charge and pKa shifting of the guanine base.

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.

Divalent ions attenuate DNA synthesis by human DNA polymerase α by changing the structure of the template/primer or by perturbing the polymerase reaction

DNA polymerases (pols) are sophisticated protein machines operating in the replication, repair and recombination of genetic material in the complex environment of the cell. DNA pol reactions require at least two divalent metal ions for the phosphodiester bond formation. We explore two understudied roles of metals in pol transactions with emphasis on polα, a crucial enzyme in the initiation of DNA synthesis. We present evidence that the combination of many factors, including the structure of the template/primer, the identity of the metal, the metal turnover in the pol active site, and the influence of the concentration of nucleoside triphosphates, affect DNA pol synthesis. On the poly-dT70 template, the increase of Mg(2+) concentration within the range typically used for pol reactions led to the severe loss of the ability of pol to extend DNA primers and led to a decline in DNA product sizes when extending RNA primers, simulating the effect of "counting" of the number of nucleotides in nascent primers by polα. We suggest that a high Mg(2+) concentration promotes the dynamic formation of unconventional DNA structure(s), thus limiting the apparent processivity of the enzyme. Next, we found that Zn(2+) supported robust polα reactions when the concentration of nucleotides was above the concentration of ions; however, there was only one nucleotide incorporation by the Klenow fragment of DNA pol I. Zn(2+) drastically inhibited polα, but had no effect on Klenow, when Mg(2+) was also present. It is possible that Zn(2+) perturbs metal-mediated transactions in pol active site, for example affecting the step of pyrophosphate removal at the end of each pol cycle necessary for continuation of polymerization.

RNA model evaluation based on MD simulation of four tRNA analogs

Characteristic of 3D structures of tRNA analogs generated from primary sequence, subjected to all-atom MD simulation and assessed using selected RNA model evaluation methods.

Crystal structure of group II intron domain 1 reveals a template for RNA assembly

Although the importance of large noncoding RNAs is increasingly appreciated, our understanding of their structures and architectural dynamics remains limited. In particular, we know little about RNA folding intermediates and how they facilitate the productive assembly of RNA tertiary structures. Here, we report the crystal structure of an obligate intermediate that is required during the earliest stages of group II intron folding. Composed of domain 1 from the Oceanobacillus iheyensis group II intron (266 nucleotides), this intermediate retains native-like features but adopts a compact conformation in which the active site cleft is closed. Transition between this closed and the open (native) conformation is achieved through discrete rotations of hinge motifs in two regions of the molecule. The open state is then stabilized by sequential docking of downstream intron domains, suggesting a 'first come, first folded' strategy that may represent a generalizable pathway for assembly of large RNA and ribonucleoprotein structures.

Inverse thio effects in the hepatitis delta virus ribozyme reveal that the reaction pathway is controlled by metal ion charge density

The hepatitis delta virus (HDV) ribozyme self-cleaves in the presence of a wide range of monovalent and divalent ions. Prior theoretical studies provided evidence that self-cleavage proceeds via a concerted or stepwise pathway, with the outcome dictated by the valency of the metal ion. In the present study, we measure stereospecific thio effects at the nonbridging oxygens of the scissile phosphate under a wide range of experimental conditions, including varying concentrations of diverse monovalent and divalent ions, and combine these with quantum mechanical/molecular mechanical (QM/MM) free energy simulations on the stereospecific thio substrates. The R_P substrate gives large normal thio effects in the presence of all monovalent ions. The S_P substrate also gives normal or no thio effects, but only for smaller monovalent and divalent cations, such as Li⁺, Mg²⁺, Ca²⁺, and Sr²⁺; in contrast, sizable inverse thio effects are found for larger monovalent and divalent cations, including Na⁺, K⁺, NH₄⁺, and Ba²⁺. Proton inventories are found to be unity in the presence of the larger monovalent and divalent ions, but two in the presence of Mg²⁺. Additionally, rate–pH profiles are inverted for the low charge density ions, and only imidazole plus ammonium ions rescue an inactive C75Δ variant in the absence of Mg²⁺. Results from the thio effect experiments, rate–pH profiles, proton inventories, and ammonium/imidazole rescue experiments, combined with QM/MM free energy simulations, support a change in the mechanism of HDV ribozyme self-cleavage from concerted and metal ion-stabilized to stepwise and proton transfer-stabilized as the charge density of the metal ion decreases.

Automated nucleic acid chain tracing in real time

The crystallographic structure solution of nucleotides and nucleotide complexes is now commonplace. The resulting electron-density maps are often poorer than for proteins, and as a result interpretation in terms of an atomic model can require significant effort, particularly in the case of large structures. While model building can be performed automatically, as with proteins, the process is time-consuming, taking minutes to days depending on the software and the size of the structure. A method is presented for the automatic building of nucleotide chains into electron density which is fast enough to be used in interactive model-building software, with extended chain fragments built around the current view position in a fraction of a second. The speed of the method arises from the determination of the 'fingerprint' of the sugar and phosphate groups in terms of conserved high-density and low-density features, coupled with a highly efficient scoring algorithm. Use cases include the rapid evaluation of an initial electron-density map, addition of nucleotide fragments to prebuilt protein structures, and in favourable cases the completion of the structure while automated model-building software is still running. The method has been incorporated into the Coot software package.

Specific RNA-binding antibodies with a four-amino-acid code

Numerous large non-coding RNAs are rapidly being discovered, and many of them have been shown to play vital roles in gene expression, gene regulation, and human diseases. Given their often structured nature, specific recognition with an antibody fragment becomes feasible and may help define the structure and function of these non-coding RNAs. As demonstrated for protein antigens, specific antibodies may aid in RNA crystal structure elucidation or the development of diagnostic tools and therapeutic drugs targeting disease-causing RNAs. Recent success and limitation of RNA antibody development has made it imperative to generate an effective antibody library specifically targeting RNA molecules. Adopting the reduced chemical diversity design and further restricting the interface diversity to tyrosines, serines, glycines, and arginines only, we have constructed a RNA-targeting Fab library. Phage display selection and downstream characterization showed that this library yielded high-affinity Fabs for all three RNA targets tested. Using a quantitative specificity assay, we found that these Fabs are highly specific, possibly due to the alternate codon design we used to avoid consecutive arginines in the Fab interface. In addition, the effectiveness of the minimal Fab library may challenge our view of the protein-RNA binding interface and provide a unique solution for future design of RNA-binding proteins.

Cobalt(III)hexaammine-dependent photocrosslinks in the hairpin ribozyme

We have utilized the hairpin ribozyme, an RNA enzyme whose structure has been solved by high-resolution methods, to develop a new tool for mapping nucleobase-stacking interactions and potential metal-binding sites in RNA molecules. This tool involves the photoactivation of a specifically bound cobalt(III)hexaammine molecule at wavelengths corresponding to excitation of the metal ion complex only; no base excitation is involved. The photoexcitation initiates a process which strongly promotes the formation of a novel covalent bond or crosslink between one base (termed the "first base"), which is close in space to the excited cobalt(III)hexaammine complex, and another base upon which the first base is closely stacked. These crosslinked species can be isolated and sequenced; their activities can be analyzed to ensure that the crosslinked structures represent an active conformation of the molecule. We have shown that, as in electron transfer in DNA, several criteria must be met to result in the successful formation of these crosslinks. These include the appropriate oxidation potential of the first donor base, the stacking and close interaction of the two donor bases involved in the crosslink, and the binding of a specific cobalt(III)hexaammine molecule to the first donor base. Additionally, we have determined that this crosslinking is pH-sensitive, although the cause of this sensitivity remains unknown. This tool has proven useful in the past for the analysis of the hairpin ribozyme folded structure, and has been applied to identify potential metal-binding sites on the hairpin and extended hammerhead ribozymes.

An RNA folding motif: GNRA tetraloop-receptor interactions

Nearly two decades after Westhof and Michel first proposed that RNA tetraloops may interact with distal helices, tetraloop–receptor interactions have been recognized as ubiquitous elements of RNA tertiary structure. The unique architecture of GNRA tetraloops (N=any nucleotide, R=purine) enables interaction with a variety of receptors, e.g., helical minor grooves and asymmetric internal loops. The most common example of the latter is the GAAA tetraloop–11 nt tetraloop receptor motif. Biophysical characterization of this motif provided evidence for the modularity of RNA structure, with applications spanning improved crystallization methods to RNA tectonics. In this review, we identify and compare types of GNRA tetraloop–receptor interactions. Then we explore the abundance of structural, kinetic, and thermodynamic information on the frequently occurring and most widely studied GAAA tetraloop–11 nt receptor motif. Studies of this interaction have revealed powerful paradigms for structural assembly of RNA, as well as providing new insights into the roles of cations, transition states and protein chaperones in RNA folding pathways. However, further research will clearly be necessary to characterize other tetraloop–receptor and long-range tertiary binding interactions in detail – an important milestone in the quantitative prediction of free energy landscapes for RNA folding.

Conformational dynamics of the human propeller telomeric DNA quadruplex on a microsecond time scale

The human telomeric DNA sequence with four repeats can fold into a parallel-stranded propeller-type topology. NMR structures solved under molecular crowding experiments correlate with the crystal structures found with crystal-packing interactions that are effectively equivalent to molecular crowding. This topology has been used for rationalization of ligand design and occurs experimentally in a number of complexes with a diversity of ligands, at least in the crystalline state. Although G-quartet stems have been well characterized, the interactions of the TTA loop with the G-quartets are much less defined. To better understand the conformational variability and structural dynamics of the propeller-type topology, we performed molecular dynamics simulations in explicit solvent up to 1.5 μs. The analysis provides a detailed atomistic account of the dynamic nature of the TTA loops highlighting their interactions with the G-quartets including formation of an A:A base pair, triad, pentad and hexad. The results present a threshold in quadruplex simulations, with regards to understanding the flexible nature of the sugar-phosphate backbone in formation of unusual architecture within the topology. Furthermore, this study stresses the importance of simulation time in sampling conformational space for this topology.

Cation-dependent folding of 3' cap-independent translation elements facilitates interaction of a 17-nucleotide conserved sequence with eIF4G

The 3′-untranslated regions of many plant viral RNAs contain cap-independent translation elements (CITEs) that drive translation initiation at the 5′-end of the mRNA. The barley yellow dwarf virus-like CITE (BTE) stimulates translation by binding the eIF4G subunit of translation initiation factor eIF4F with high affinity. To understand this interaction, we characterized the dynamic structural properties of the BTE, mapped the eIF4G-binding sites on the BTE and identified a region of eIF4G that is crucial for BTE binding. BTE folding involves cooperative uptake of magnesium ions and is driven primarily by charge neutralization. Footprinting experiments revealed that functional eIF4G fragments protect the highly conserved stem–loop I and a downstream bulge. The BTE forms a functional structure in the absence of protein, and the loop that base pairs the 5′-untranslated region (5′-UTR) remains solvent-accessible at high eIF4G concentrations. The region in eIF4G between the eIF4E-binding site and the MIF4G region is required for BTE binding and translation. The data support the model in which the eIF4F complex binds directly to the BTE which base pairs simultaneously to the 5′-UTR, allowing eIF4F to recruit the 40S ribosomal subunit to the 5′-end.

Utilizing the GAAA tetraloop/receptor to facilitate crystal packing and determination of the structure of a CUG RNA helix

Myotonic dystrophy type 1 (DM1) is a microsatellite expansion disorder caused by the aberrant expansion of CTG repeats in the 3'-untranslated region of the DMPK gene. When transcribed, the toxic RNA CUG repeats sequester RNA binding proteins, which leads to disease symptoms. The expanded CUG repeats can adopt a double-stranded structure, and targeting this helix is a therapeutic strategy for DM1. To improve our understanding of the 5'CUG/3'GUC motif and how it may interact with proteins and small molecules, we designed a short CUG helix attached to a GAAA tetraloop/receptor to facilitate crystal packing. Here we report the highest-resolution structure (1.95 Å) to date of a GAAA tetraloop/receptor and the CUG helix it was used to crystallize. Within the CUG helix, we identify two different forms of noncanonical U-U pairs and reconfirm that CUG repeats are essentially A-form. An analysis of all noncanonical U-U pairs in the context of CUG repeats revealed six different classes of conformations that the noncanonical U-U pairs are able to adopt.

Metal-ion rescue revisited: biochemical detection of site-bound metal ions important for RNA folding

Within the three-dimensional architectures of RNA molecules, divalent metal ions populate specific locations, shedding their water molecules to form chelates. These interactions help the RNA adopt and maintain specific conformations and frequently make essential contributions to function. Defining the locations of these site-bound metal ions remains challenging despite the growing database of RNA structures. Metal-ion rescue experiments have provided a powerful approach to identify and distinguish catalytic metal ions within RNA active sites, but the ability of such experiments to identify metal ions that contribute to tertiary structure acquisition and structural stability is less developed and has been challenged. Herein, we use the well-defined P4-P6 RNA domain of the Tetrahymena group I intron to reevaluate prior evidence against the discriminatory power of metal-ion rescue experiments and to advance thermodynamic descriptions necessary for interpreting these experiments. The approach successfully identifies ligands within the RNA that occupy the inner coordination sphere of divalent metal ions and distinguishes them from ligands that occupy the outer coordination sphere. Our results underscore the importance of obtaining complete folding isotherms and establishing and evaluating thermodynamic models in order to draw conclusions from metal-ion rescue experiments. These results establish metal-ion rescue as a rigorous tool for identifying and dissecting energetically important metal-ion interactions in RNAs that are noncatalytic but critical for RNA tertiary structure.

Folding path of P5abc RNA involves direct coupling of secondary and tertiary structures

Folding mechanisms in which secondary structures are stabilized through the formation of tertiary interactions are well documented in protein folding but challenge the folding hierarchy normally assumed for RNA. However, it is increasingly clear that RNA could fold by a similar mechanism. P5abc, a small independently folding tertiary domain of the Tetrahymena thermophila group I ribozyme, is known to fold by a secondary structure rearrangement involving helix P5c. However, the extent of this rearrangement and the precise stage of folding that triggers it are unknown. We use experiments and simulations to show that the P5c helix switches to the native secondary structure late in the folding pathway and is directly coupled to the formation of tertiary interactions in the A-rich bulge. P5c mutations show that the switch in P5c is not rate-determining and suggest that non-native interactions in P5c aid folding rather than impede it. Our study illustrates that despite significant differences in the building blocks of proteins and RNA, there may be common ways in which they self-assemble.

Crystal structure of an RNA polymerase ribozyme in complex with an antibody fragment

All models of the RNA world era invoke the presence of ribozymes that can catalyse RNA polymerization. The class I ligase ribozyme selected in vitro 15 years ago from a pool of random RNA sequences catalyses formation of a 3',5'-phosphodiester linkage analogous to a single step of RNA polymerization. Recently, the three-dimensional structure of the ligase was solved in complex with U1A RNA-binding protein and independently in complex with an antibody fragment. The RNA adopts a tripod arrangement and appears to use a two-metal ion mechanism similar to protein polymerases. Here, we discuss structural implications for engineering a true polymerase ribozyme and describe the use of the antibody framework both as a portable chaperone for crystallization of other RNAs and as a platform for exploring steps in evolution from the RNA world to the RNA-protein world.

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

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

Characterization of metal ion-nucleic acid interactions in solution

Metal ions are inextricably involved with nucleic acids due to their polyanionic nature. In order to understand the structure and function of RNAs and DNAs, one needs to have detailed pictures on the structural, thermodynamic, and kinetic properties of metal ion interactions with these biomacromolecules. In this review we first compile the physicochemical properties of metal ions found and used in combination with nucleic acids in solution. The main part then describes the various methods developed over the past decades to investigate metal ion binding by nucleic acids in solution. This includes for example hydrolytic and radical cleavage experiments, mutational approaches, as well as kinetic isotope effects. In addition, spectroscopic techniques like EPR, lanthanide(III) luminescence, IR and Raman as well as various NMR methods are summarized. Aside from gaining knowledge about the thermodynamic properties on the metal ion-nucleic acid interactions, especially NMR can be used to extract information on the kinetics of ligand exchange rates of the metal ions applied. The final section deals with the influence of anions, buffers, and the solvent permittivity on the binding equilibria between metal ions and nucleic acids. Little is known on some of these aspects, but it is clear that these three factors have a large influence on the interaction between metal ions and nucleic acids.

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.

Two distinct catalytic strategies in the hepatitis δ virus ribozyme cleavage reaction

The hepatitis delta virus (HDV) ribozyme and related RNAs are widely dispersed in nature. This RNA is a small nucleolytic ribozyme that self-cleaves to generate products with a 2',3'-cyclic phosphate and a free 5'-hydroxyl. Although small ribozymes are dependent on divalent metal ions under biologically relevant buffer conditions, they function in the absence of divalent metal ions at high ionic strengths. This characteristic suggests that a functional group within the covalent structure of small ribozymes is facilitating catalysis. Structural and mechanistic analyses have demonstrated that the HDV ribozyme active site contains a cytosine with a perturbed pK(a) that serves as a general acid to protonate the leaving group. The reaction of the HDV ribozyme in monovalent cations alone never approaches the velocity of the Mg(2+)-dependent reaction, and there is significant biochemical evidence that a Mg(2+) ion participates directly in catalysis. A recent crystal structure of the HDV ribozyme revealed that there is a metal binding pocket in the HDV ribozyme active site. Modeling of the cleavage site into the structure suggested that this metal ion can interact directly with the scissile phosphate and the nucleophile. In this manner, the Mg(2+) ion can serve as a Lewis acid, facilitating deprotonation of the nucleophile and stabilizing the conformation of the cleavage site for in-line attack of the nucleophile at the scissile phosphate. This catalytic strategy had previously been observed only in much larger ribozymes. Thus, in contrast to most large and small ribozymes, the HDV ribozyme uses two distinct catalytic strategies in its cleavage reaction.

Improvement of the crystallizability and expression of an RNA crystallization chaperone

Crystallizing RNA has been an imperative and challenging task in the world of RNA research. Assistive methods such as chaperone-assisted RNA crystallography (CARC), employing monoclonal antibody fragments (Fabs) as crystallization chaperones have enabled us to obtain RNA crystal structures by forming crystal contacts and providing initial phasing information. Despite the early successes, the crystallization of large RNA-Fab complex remains a challenge in practice. The possible reason for this difficulty is that the Fab scaffold has not been optimized for crystallization in complex with RNA. Here, we have used the surface entropy reduction (SER) technique for the optimization of ΔC209 P4-P6/Fab2 model system. Protruding lysine and glutamate residues were mutated to a set of alanines or serines to construct Fab2SMA or Fab2SMS. Expression with the shake flask approach was optimized to allow large scale production for crystallization. Crystal screening shows that significantly higher crystal-forming ratio was observed for the mutant complexes. As the chosen SER residues are far away from the CDR regions of the Fab, the same set of mutations can now be directly applied to other Fabs binding to a variety of ribozymes and riboswitches to improve the crystallizability of Fab-RNA complex.

A new way to see RNA

Unlike proteins, the RNA backbone has numerous degrees of freedom (eight, if one counts the sugar pucker), making RNA modeling, structure building and prediction a multidimensional problem of exceptionally high complexity. And yet RNA tertiary structures are not infinite in their structural morphology; rather, they are built from a limited set of discrete units. In order to reduce the dimensionality of the RNA backbone in a physically reasonable way, a shorthand notation was created that reduced the RNA backbone torsion angles to two (η and θ, analogous to φ and ψ in proteins). When these torsion angles are calculated for nucleotides in a crystallographic database and plotted against one another, one obtains a plot analogous to a Ramachandran plot (the η/θ plot), with highly populated and unpopulated regions. Nucleotides that occupy proximal positions on the plot have identical structures and are found in the same units of tertiary structure. In this review, we describe the statistical validation of the η/θ formalism and the exploration of features within the η/θ plot. We also describe the application of the η/θ formalism in RNA motif discovery, structural comparison, RNA structure building and tertiary structure prediction. More than a tool, however, the η/θ formalism has provided new insights into RNA structure itself, revealing its fundamental components and the factors underlying RNA architectural form.

Insights into metalloregulation by M-box riboswitch RNAs via structural analysis of manganese-bound complexes

The M-box riboswitch couples intracellular magnesium levels to expression of bacterial metal transport genes. Structural analyses on other riboswitch RNA classes, which typically respond to a small organic metabolite, have revealed that ligand recognition occurs through a combination of base-stacking, electrostatic, and hydrogen-bonding interactions. In contrast, the M-box RNA triggers a change in gene expression upon association with an undefined population of metals, rather than responding to only a single ligand. Prior biophysical experimentation suggested that divalent ions associate with the M-box RNA to promote a compacted tertiary conformation, resulting in sequestration of a short sequence tract otherwise required for downstream gene expression. Electrostatic shielding from loosely associated metals is undoubtedly an important influence during this metal-mediated compaction pathway. However, it is also likely that a subset of divalent ions specifically occupies cation binding sites and promotes proper positioning of functional groups for tertiary structure stabilization. To better elucidate the role of these metal binding sites, we resolved a manganese-chelated M-box RNA complex to 1.86 Å by X-ray crystallography. These data support the presence of at least eight well-ordered cation binding pockets, including several sites that had been predicted by biochemical studies but were not observed in prior structural analysis. Overall, these data support the presence of three metal-binding cores within the M-box RNA that facilitate a network of long-range interactions within the metal-bound, compacted conformation.

Bidentate RNA-magnesium clamps: on the origin of the special role of magnesium in RNA folding

Magnesium plays a special role in RNA function and folding. Although water is magnesium's most common first-shell ligand, the oxyanions of RNA have significant affinity for magnesium. Here we provide a quantum mechanical description of first-shell RNA-magnesium and DNA-magnesium interactions, demonstrating the unique features that characterize the energetics and geometry of magnesium complexes within large folded RNAs. Our work focuses on bidentate chelation of magnesium by RNA or DNA, where multiple phosphate oxyanions enter the first coordination shell of magnesium. These bidentate RNA clamps of magnesium occur frequently in large RNAs. The results here suggest that magnesium, compared to calcium and sodium, has an enhanced ability to form bidentate clamps with RNA. Bidentate RNA-sodium clamps, in particular, are unstable and spontaneously open. Due to magnesium's size and charge density it binds more intimately than other cations to the oxyanions of RNA, so that magnesium clamps are stabilized not only by electrostatic interactions, but also by charge transfer, polarization, and exchange interactions. These nonelectrostatic components of the binding are quite substantial with the high charge and small interatomic distances within the magnesium complexes, but are less pronounced for calcium due to its larger size, and for sodium due to its smaller charge. Additionally, bidentate RNA clamps of magnesium are more stable than those with DNA. The source of the additional stability of RNA complexes is twofold: there is a slightly attenuated energetic penalty for ring closure in the formation of RNA bidentate chelation complexes and elevated electrostatic interactions between the RNA and cations. In sum, it can be seen why sodium and calcium cannot replicate the structures or energetics of RNA-magnesium complexes.

Engineering cooperative tecto-RNA complexes having programmable stoichiometries

High affinity and specificity RNA-RNA binding interfaces can be constructed by combining pairs of GNRA loop/loop-receptor interaction motifs. These interactions can be fused using flexible four-way junction motifs to create divalent, self-assembling scaffolding units ('tecto-RNA') that have favorable properties for nanomedicine and other applications. We describe the design and directed assembly of tecto-RNA units ranging from closed, cooperatively assembling ring-shaped complexes of programmable stoichiometries (dimers, trimers and tetramers) to open multimeric structures. The novelty of this work is that tuning of the stoichiometries of self-assembled complexes is achieved by precise positioning of the interaction motifs in the monomer units rather than changing their binding specificities. Structure-probing and transmission electron microscopy studies as well as thermodynamic analysis support formation of closed cooperative complexes that are highly resistant to nuclease digestion. The present designs provide two helical arms per RNA monomer for further functionalization aims.

Structural basis for recognition of S-adenosylhomocysteine by riboswitches

S-adenosyl-(L)-homocysteine (SAH) riboswitches are regulatory elements found in bacterial mRNAs that up-regulate genes involved in the S-adenosyl-(L)-methionine (SAM) regeneration cycle. To understand the structural basis of SAH-dependent regulation by RNA, we have solved the structure of its metabolite-binding domain in complex with SAH. This structure reveals an unusual pseudoknot topology that creates a shallow groove on the surface of the RNA that binds SAH primarily through interactions with the adenine ring and methionine main chain atoms and discriminates against SAM through a steric mechanism. Chemical probing and calorimetric analysis indicate that the unliganded RNA can access bound-like conformations that are significantly stabilized by SAH to direct folding of the downstream regulatory switch. Strikingly, we find that metabolites bearing an adenine ring, including ATP, bind this aptamer with sufficiently high affinity such that normal intracellular concentrations of these compounds may influence regulation of the riboswitch.

Identification and characterization of anion binding sites in RNA

Although RNA molecules are highly negatively charged, anions have been observed bound to RNA in crystal structures. It has been proposed that anion binding sites found within isolated RNAs represent regions of the molecule that could be involved in intermolecular interactions, indicating potential contact points for negatively charged amino acids from proteins or phosphate groups from an RNA. Several types of anion binding sites have been cataloged based on available structures. However, currently there is no method for unambiguously assigning anions to crystallographic electron density, and this has precluded more detailed analysis of RNA-anion interaction motifs and their significance. We therefore soaked selenate into two different types of RNA crystals and used the anomalous signal from these anions to identify binding sites in these RNA molecules unambiguously. Examination of these sites and comparison with other suspected anion binding sites reveals features of anion binding motifs, and shows that selenate may be a useful tool for studying RNA-anion interactions.

A quantitative Raman spectroscopic signal for metal-phosphodiester interactions in solution

Accurate identification and quantification of metal ion-phosphodiester interactions are essential for understanding the role of metal ions as determinants of three-dimensional folding of large RNAs and as cofactors in the active sites of both RNA and protein phosphodiesterases. Accomplishing this goal is difficult due to the dynamic and complex mixture of direct and indirect interactions formed with nucleic acids and other phosphodiesters in solution. To address this issue, Raman spectroscopy has been used to measure changes in bond vibrational energies due to metal interactions. However, the contributions of inner-sphere, H-bonding, and electrostatic interactions to the Raman spectrum of phosphoryl oxygens have not been analyzed quantitatively. Here, we report that all three forms of metal ion interaction result in attenuation of the Raman signal for the symmetric vibration of the nonbridging phosphate oxygens (nu(s)PO(2)(-)), while only inner-sphere coordination gives rise to an apparent shift of nu(s)PO(2)(-) to higher wavenumbers (nu(s)PO(2)(-)M) in solution. Formation of nu(s)PO(2)(-)M is shown to be both dependent on metal ion identity and an accurate measure of site-specific metal ion binding. In addition, the spectroscopic parameter reflecting the energetic difference between nu(s)PO(2)(-) and nu(s)PO(2)(-)M (DeltanuM) is largely insensitive to changes in phosphodiester structure but strongly dependent on the absolute electronegativity and hardness of the interacting metal ion. Together, these studies provide strong experimental support for the use of nu(s)PO(2)(-)M and DeltanuM as general spectroscopic features for the quantitative analysis of metal binding affinity and the identification of metal ions associated with phosphodiesters in solution.

Modulation of RNA metal binding by flanking bases: 15N NMR evaluation of GC, tandem GU, and tandem GA sites

(15)N NMR chemical shift changes in the presence of Mg(H(2)O)(6)(2+), Zn(2+), Cd(2+), and Co(NH(3))(6)(3+) were used to probe the effect of flanking bases on metal binding sites in three different RNA motifs. We found that: for GC pairs, the presence of a flanking purine creates a site for the soft metals Zn(2+) and Cd(2+) only; a GG.UU motif selectively binds only Co(NH(3))(6)(3+), while a UG.GU motif binds none of these metals; a 3' guanosine flanking the adenosine of a sheared GA.AG pair creates an unusually strong binding site that precludes binding to the cross-strand stacked guanosines within the tandem pair.

Discrimination between closely related cellular metabolites by the SAM-I riboswitch

The SAM-I riboswitch is a cis-acting element of genetic control found in bacterial mRNAs that specifically binds S-adenosylmethionine (SAM). We previously determined the 2.9-A X-ray crystal structure of the effector-binding domain of this RNA element, revealing details of RNA-ligand recognition. To improve this structure, variations were made to the RNA sequence to alter lattice contacts, resulting in a 0.5-A improvement in crystallographic resolution and allowing for a more accurate refinement of the crystallographic model. The basis for SAM specificity was addressed by a structural analysis of the RNA complexed to S-adenosylhomocysteine (SAH) and sinefungin and by measuring the affinity of SAM and SAH for a series of mutants using isothermal titration calorimetry. These data illustrate the importance of two universally conserved base pairs in the RNA that form electrostatic interactions with the positively charged sulfonium group of SAM, thereby providing a basis for discrimination between SAM and SAH.

EPR methods to study specific metal-ion binding sites in RNA

The properties of metal-ion interactions with RNA can be explored by spectroscopic methods. In this chapter, we describe the use of paramagnetic Mn(2+) ions and electron paramagnetic resonance (EPR)-based techniques to monitor the association of Mn(2+) with RNA and related nucleotides. Solution EPR methods are used to determine the numbers of Mn(2+) ions associating with RNA. For RNA poised with a single-bound Mn(2+), low-temperature EPR characteristics provide information about the asymmetry of the Mn(2+) coordination site. To identify the RNA groups coordinating to the Mn(2+) ion, ENDOR (electron nuclear double resonance) and ESEEM (electron spin echo envelope modulation) methods are applied. Both continuous-wave (CW) and electron spin echo (ESE)-detected ENDOR methods are described. This chapter includes practical details for RNA sample preparation, including isotope substitution and cryoprotection, and an overview of data acquisition and analysis methods used in these techniques, as well as examples from the current literature.

2'-amino-modified ribonucleotides as probes for local interactions within RNA

The 2'-hydroxyl group plays an integral role in RNA structure and catalysis. This ubiquitous component of the RNA backbone can participate in multiple interactions essential for RNA function, such as hydrogen bonding and metal ion coordination, but the multifunctional nature of the 2'-hydroxyl renders identification of these interactions a significant challenge. By virtue of their versatile physicochemical properties, such as distinct metal coordination preferences, hydrogen bonding properties, and ability to be protonated, 2'-amino-2'-deoxyribonucleotides can serve as tools for probing local interactions involving 2'-hydroxyl groups within RNA. The 2'-amino group can also serve as a chemoselective site for covalent modification, permitting the introduction of probes for investigation of RNA structure and dynamics. In this chapter, we describe the use of 2'-aminonucleotides for investigation of local interactions within RNA, focusing on interactions involving 2'-hydroxyl groups required for RNA structure, function, and catalysis.

Preparation of group I introns for biochemical studies and crystallization assays by native affinity purification

The study of functional RNAs of various sizes and structures requires efficient methods for their synthesis and purification. Here, 23 group I intron variants ranging in length from 246 to 341 nucleotides—some containing exons—were subjected to a native purification technique previously applied only to shorter RNAs (<160 nucleotides). For the RNAs containing both exons, we adjusted the original purification protocol to allow for purification of radiolabeled molecules. The resulting RNAs were used in folding assays on native gel electrophoresis and in self-splicing assays. The intron-only RNAs were subjected to the regular native purification scheme, assayed for folding and employed in crystallization screens. All RNAs that contained a 3′ overhang of one nucleotide were efficiently cleaved off from the support and were at least 90% pure after the non-denaturing purification. A representative subset of these RNAs was shown to be folded and self-splicing after purification. Additionally, crystals were grown for a 286 nucleotide long variant of the Clostridium botulinum intron. These results demonstrate the suitability of the native affinity purification method for the preparation of group I introns. We hope these findings will stimulate a broader application of this strategy to the preparation of other large RNA molecules.