Structural insights into group II intron catalysis and branch-site selection

Single-Molecule X-ray Scattering Used to Visualize the Conformation Distribution of Biological Molecules via Single-Object Scattering Sampling

Biological macromolecules, the fundamental building blocks of life, exhibit dynamic structures in their natural environment. Traditional structure determination techniques often oversimplify these multifarious conformational spectra by capturing only ensemble- and time-averaged molecular structures. Addressing this gap, in this work, we extend the application of the single-object scattering sampling (SOSS) method to diverse biological molecules, including RNAs and proteins. Our approach, referred to as "Bio-SOSS", leverages ultrashort X-ray pulses to capture instantaneous structures. In Bio-SOSS, we employ two gold nanoparticles (AuNPs) as labels, which provide strong contrast in the X-ray scattering signal, to ensure precise distance determinations between labeled sites. We generated hypothetical Bio-SOSS images for RNAs, proteins, and an RNA-protein complex, each labeled with two AuNPs at specified positions. Subsequently, to validate the accuracy of Bio-SOSS, we extracted distances between these nanoparticle labels from the images and compared them with the actual values used to generate the Bio-SOSS images. Specifically, for a representative RNA (1KXK), the standard deviation in distance discrepancies between molecular dynamics snapshots and Bio-SOSS retrievals was found to be optimally around 0.2 Å, typically within 1 Å under practical experimental conditions at state-of-the-art X-ray free-electron laser facilities. Furthermore, we conducted an in-depth analysis of how various experimental factors, such as AuNP size, X-ray properties, and detector geometry, influence the accuracy of Bio-SOSS. This comprehensive investigation highlights the practicality and potential of Bio-SOSS in accurately capturing the diverse conformation spectrum of biological macromolecules, paving the way for deeper insights into their dynamic natures.

Translational and rotational diffusion of short ribonucleic acids

Inevitable precondition for ribonucleic acids to regulate gene expression and to perform gene editing is diffusion. Free three-dimensional translational diffusion velocity of RNA of up to 200 nucleotides could be predicted with high accuracy by the empirical model D = 4.58 10^-10 N^-0.39 m²s^-1. Furthermore, the biological function of ribonucleic acids is determined by rotational diffusion. In the presented work, an empirical model is derived applying atom-level shell-modeling of electron density maps, D_r = 1.62 10⁹ N^-1.20 s^-1, to predict the rotational diffusion coefficient of short ribonucleic acids based on the polymer size.

Selenium-atom-modified thymidine enhances the specificity and sensitivity of DNA polymerization and detection

Nucleobase mismatches can jeopardize DNA polymerization specificity, causing mutations and errors in DNA replication and detection. Herein we report the first synthesis of novel 2-Se-thymidine triphosphate (SeTTP), describe the single-selenium atom-specific modification strategy (SAM) against T/G mismatches, and demonstrate SAM-assisted polymerization and detection with much higher specificity and sensitivity. SAM can effectively suppress the formation of non-specific products in DNA polymerization and detection. Thus, SAM enhances the specificity of DNA synthesis by approximately 10 000 fold, and in turn, it allows the detection of clinical COVID-19 viral RNA in low copy numbers (single-digit copies), while the conventional RT-qPCR does not.

Branch site bulge conformations in domain 6 determine functional sugar puckers in group II intron splicing

Although group II intron ribozymes are intensively studied the question how structural dynamics affects splicing catalysis has remained elusive. We report for the first time that the group II intron domain 6 exists in a secondary structure equilibrium between a single- and a two-nucleotide bulge conformation, which is directly linked to a switch between sugar puckers of the branch site adenosine. Our study determined a functional sugar pucker equilibrium between the transesterification active C2'-endo conformation of the branch site adenosine in the 1nt bulge and an inactive C3'-endo state in the 2nt bulge fold, allowing the group II intron to switch its activity from the branching to the exon ligation step. Our detailed NMR spectroscopic investigation identified magnesium (II) ions and the branching reaction as regulators of the equilibrium populations. The tuneable secondary structure/sugar pucker equilibrium supports a conformational selection mechanism to up- and downregulate catalytically active and inactive states of the branch site adenosine to orchestrate the multi-step splicing process. The conformational dynamics of group II intron domain 6 is also proposed to be a key aspect for the directionality selection in reversible splicing.

QRNAS: software tool for refinement of nucleic acid structures

BACKGROUND

Computational models of RNA 3D structure often present various inaccuracies caused by simplifications used in structure prediction methods, such as template-based modeling or coarse-grained simulations. To obtain a high-quality model, the preliminary RNA structural model needs to be refined, taking into account atomic interactions. The goal of the refinement is not only to improve the local quality of the model but to bring it globally closer to the true structure.

RESULTS

We present QRNAS, a software tool for fine-grained refinement of nucleic acid structures, which is an extension of the AMBER simulation method with additional restraints. QRNAS is capable of handling RNA, DNA, chimeras, and hybrids thereof, and enables modeling of nucleic acids containing modified residues.

CONCLUSIONS

We demonstrate the ability of QRNAS to improve the quality of models generated with different methods. QRNAS was able to improve MolProbity scores of NMR structures, as well as of computational models generated in the course of the RNA-Puzzles experiment. The overall geometry improvement may be associated with increased model accuracy, especially on the level of correctly modeled base-pairs, but the systematic improvement of root mean square deviation to the reference structure should not be expected. The method has been integrated into a computational modeling workflow, enabling improved RNA 3D structure prediction.

Influence of pH and Mg(ii) on the catalytic core domain 5 of a bacterial group II intron

RNA molecules fold into complex structures that allow them to perform specific functions. To compensate the relative lack of diversity of functional groups within nucleotides, metal ions work as crucial co-factors. In addition, shifted pK_as are observed in RNA, enabling acid-base reactions at ambient pH. The central catalytic domain 5 (D5) hairpin of the Azotobacter vinelandii group II intron undergoes both metal ion binding and pH dependence, presumably playing an important functional role in the ribozyme's reaction. By NMR spectroscopy we have here characterized the metal ion binding sites and affinities for the hairpin's internal G-A mismatch, bulge, and pentaloop. The influence of Mg(ii) and pH on the local conformation of the catalytically crucial region is also explored by fluorescence spectroscopy.

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.

Molecular Mechanism and Evolution of Nuclear Pre-mRNA and Group II Intron Splicing: Insights from Cryo-Electron Microscopy Structures

Nuclear pre-mRNA splicing and group II intron self-splicing both proceed by two-step transesterification reactions via a lariat intron intermediate. Recently determined cryo-electron microscopy (cryo-EM) structures of catalytically active spliceosomes revealed the RNA-based catalytic core and showed how pre-mRNA substrates and reaction products are positioned in the active site. These findings highlight a strong structural similarity to the group II intron active site, strengthening the notion that group II introns and spliceosomes evolved from a common ancestor. Prp8, the largest and most conserved protein in the spliceosome, cradles the active site RNA. Prp8 and group II intron maturase have a similar domain architecture, suggesting that they also share a common evolutionary origin. The interactions between maturase and key group II intron RNA elements, such as the exon-binding loop and domains V and VI, are recapitulated in the interactions between Prp8 and key elements in the spliceosome's catalytic RNA core. Structural comparisons suggest that the extensive RNA scaffold of the group II intron was gradually replaced by proteins as the spliceosome evolved. A plausible model of spliceosome evolution is discussed.

Structure determination of group II introns

Group II introns are self-splicing catalytic RNAs that are able to excise themselves from pre-mRNAs using a mechanism identical to that utilized by the spliceosome. Both structural and phylogenetic data support the hypothesis that group II introns and the spliceosome share a common ancestor. Structures of group II introns have given insight into the active site required for the catalysis of RNA splicing. This review outlines crucial aspects of the structure determination of group II introns such as sample preparation and data processing. Given that group II introns are large RNAs that must be synthesized through in vitro transcription, there are special considerations that must be taken into account in terms of purification and crystallization, as compared to the isolation of large intact ribonucleoprotein complexes such as the ribosome. We specifically focus on the methodology used to determine the structure of the eukaryotic group II intron lariat from the brown algae Pylaiella littoralis. The techniques described in this review can also be applied for the structure determination of other large RNAs.

Crystal structures of a group II intron lariat primed for reverse splicing

The 2'-5' branch of nuclear premessenger introns is believed to have been inherited from self-splicing group II introns, which are retrotransposons of bacterial origin. Our crystal structures at 3.4 and 3.5 angstrom of an excised group II intron in branched ("lariat") form show that the 2'-5' branch organizes a network of active-site tertiary interactions that position the intron terminal 3'-hydroxyl group into a configuration poised to initiate reverse splicing, the first step in retrotransposition. Moreover, the branchpoint and flanking helices must undergo a base-pairing switch after branch formation. A group II-based model of the active site of the nuclear splicing machinery (the spliceosome) is proposed. The crucial role of the lariat conformation in active-site assembly and catalysis explains its prevalence in modern splicing.

Using Group II Introns for Attenuating the In Vitro and In Vivo Expression of a Homing Endonuclease

In Chaetomium thermophilum (DSM 1495) within the mitochondrial DNA (mtDNA) small ribosomal subunit (rns) gene a group IIA1 intron interrupts an open reading frame (ORF) encoded within a group I intron (mS1247). This arrangement offers the opportunity to examine if the nested group II intron could be utilized as a regulatory element for the expression of the homing endonuclease (HEase). Constructs were generated where the codon-optimized ORF was interrupted with either the native group IIA1 intron or a group IIB type intron. This study showed that the expression of the HEase (in vivo) in Escherichia coli can be regulated by manipulating the splicing efficiency of the HEase ORF-embedded group II introns. Exogenous magnesium chloride (MgCl₂) stimulated the expression of a functional HEase but the addition of cobalt chloride (CoCl₂) to growth media antagonized the expression of HEase activity. Ultimately the ability to attenuate HEase activity might be useful in precision genome engineering, minimizing off target activities, or where pathways have to be altered during a specific growth phase.

Using Molecular Replacement Phasing to Study the Structure and Function of RNA

In recent years a wide variety of RNA molecules regulating fundamental cellular processes has been discovered. Therefore, RNA structure determination is experiencing a boost and many more RNA structures are likely to be determined in the years to come. The broader availability of experimentally determined RNA structures implies that molecular replacement (MR) will be used more and more frequently as a method for phasing future crystallographic structures. In this report we describe various aspects relative to RNA structure determination by MR. First, we describe how to select and create MR search models for nucleic acids. Second, we describe how to perform MR searches on RNA using available crystallographic software. Finally, we describe how to refine and interpret the successful MR solutions. These protocols are applicable to determine novel RNA structures as well as to establish structural-functional relationships on existing RNA structures.

Re-refinement of the spliceosomal U4 snRNP core-domain structure

The core domain of small nuclear ribonucleoprotein (snRNP), comprised of a ring of seven paralogous proteins bound around a single-stranded RNA sequence, functions as the assembly nucleus in the maturation of U1, U2, U4 and U5 spliceosomal snRNPs. The structure of the human U4 snRNP core domain was initially solved at 3.6 Å resolution by experimental phasing using data with tetartohedral twinning. Molecular replacement from this model followed by density modification using untwinned data recently led to a structure of the minimal U1 snRNP at 3.3 Å resolution. With the latter structure providing a search model for molecular replacement, the U4 core-domain structure has now been re-refined. The U4 Sm site-sequence AAUUUUU has been shown to bind to the seven Sm proteins SmF-SmE-SmG-SmD3-SmB-SmD1-SmD2 in an identical manner as the U1 Sm-site sequence AAUUUGU, except in SmD1 where the bound U replaces G. The progression from the initial to the re-refined structure exemplifies a tortuous route to accuracy: where well diffracting crystals of complex assemblies are initially unavailable, the early model errors are rectified by exploiting preliminary interpretations in further experiments involving homologous structures. New insights are obtained from the more accurate model.

De novo tertiary structure prediction using RNA123--benchmarking and application to Macugen

The present benchmarking study utilizes the RNA123 program for de novo prediction of tertiary structures of a set of 50 RNA molecules for which X-ray/NMR structures are available, based on the nucleic acid sequence only. All molecules contain a hairpin loop motif and a helical structure of canonical and non-canonical base pairs, interrupted by bulges and internal loops to various degrees. RNA molecules with double helices made up purely by canonical base pairing, and molecules containing symmetric internal loops of non-canonical base pairing are, overall, very well predicted. Structures containing bulges and asymmetric internal loops, and more complex structures containing multiple bulges and internal loops in the same molecule, result in larger deviations from their X-ray/NMR predicted structures due to higher degree of flexibility of the nucleotide bases in these regions. In a majority of the molecules included herein, the RNA123 program was, however, able to predict the tertiary structure with a heavy atom RMSD of less than 5 Å to the X-ray/NMR structure, and the models were in most cases structurally closer to the X-ray/NMR structures than models predicted by MC-Fold and MC-Sym. A set of RNA molecules containing pseudoknot tertiary structure motifs were included, but neither of the programs was able to predict the folding of the single-stranded stem onto the helix without additional structural input. The RNA123 program was then applied to predict the tertiary structure of the RNA segment of Macugen®, the first RNA aptamer approved for clinical use, and for which no tertiary structure has yet been solved. Four possible tertiary structures were predicted for this 27-nucleic-acid-long RNA molecule, which will be used in constructing a full model of the PEGylated aptamer and its interaction with the vascular endothelial growth factor target.

Solving nucleic acid structures by molecular replacement: examples from group II intron studies

Structured RNA molecules are key players in ensuring cellular viability. It is now emerging that, like proteins, the functions of many nucleic acids are dictated by their tertiary folds. At the same time, the number of known crystal structures of nucleic acids is also increasing rapidly. In this context, molecular replacement will become an increasingly useful technique for phasing nucleic acid crystallographic data in the near future. Here, strategies to select, create and refine molecular-replacement search models for nucleic acids are discussed. Using examples taken primarily from research on group II introns, it is shown that nucleic acids are amenable to different and potentially more flexible and sophisticated molecular-replacement searches than proteins. These observations specifically aim to encourage future crystallographic studies on the newly discovered repertoire of noncoding transcripts.

Improving experimental phases for strong reflections prior to density modification

Experimental phasing of diffraction data from macromolecular crystals involves deriving phase probability distributions. These distributions are often bimodal, making their weighted average, the centroid phase, improbable, so that electron-density maps computed using centroid phases are often non-interpretable. Density modification brings in information about the characteristics of electron density in protein crystals. In successful cases, this allows a choice between the modes in the phase probability distributions, and the maps can cross the borderline between non-interpretable and interpretable. Based on the suggestions by Vekhter [Vekhter (2005), Acta Cryst. D61, 899-902], the impact of identifying optimized phases for a small number of strong reflections prior to the density-modification process was investigated while using the centroid phase as a starting point for the remaining reflections. A genetic algorithm was developed that optimizes the quality of such phases using the skewness of the density map as a target function. Phases optimized in this way are then used in density modification. In most of the tests, the resulting maps were of higher quality than maps generated from the original centroid phases. In one of the test cases, the new method sufficiently improved a marginal set of experimental SAD phases to enable successful map interpretation. A computer program, SISA, has been developed to apply this method for phase improvement in macromolecular crystallography.

Structure and function of the T-loop structural motif in noncoding RNAs

The T-loop is a frequently occurring five-nucleotide motif found in the structure of noncoding RNAs where it is commonly assumed to play an important role in stabilizing the tertiary RNA structure by facilitating long-range interactions between different regions of the molecule. T-loops were first identified in tRNA(Phe) and a formal consensus sequence for this motif was formulated and later revised based on analyses of the crystal structures of prokaryotic ribosomal RNAs and RNase P and the corresponding primary sequence of their orthologues. In the past decade, several new structures of large RNA molecules have been added to the RCSB Protein Data Bank, including the eukaryotic ribosome, a self-splicing group II intron, numerous synthetases in complex with their cognate transfer RNAs (tRNAs), transfer-messenger RNA (tmRNA) in complex with SmpB, several riboswitches, and a complex of bacterial RNase P bound to its tRNA substrate. In this review, the search for T-loops is extended to these new RNA molecules based on the previously established structure-based criteria. The review highlights and discusses the function and additional roles of T-loops in four broad categories of RNA molecules, namely tRNAs, ribosomal RNAs (rRNAs), P RNAs, and RNA genetic elements. Additionally, the potential application for T-loops as interaction modules is also discussed.

The structural stabilization of the κ three-way junction by Mg(II) represents the first step in the folding of a group II intron

Folding of group II introns is characterized by a first slow compaction of domain 1 (D1) followed by the rapid docking of other domains to this scaffold. D1 compaction initiates in a small subregion encompassing the κ and ζ elements. These two tertiary elements are also the major interaction sites with domain 5 to form the catalytic core. Here, we provide the first characterization of the structure adopted at an early folding step and show that the folding control element can be narrowed down to the three-way junction with the κ motif. In our nuclear magnetic resonance studies of this substructure derived from the yeast mitochondrial group II intron Sc.ai5γ, we show that a high affinity Mg(II) ion stabilizes the κ element and enables coaxial stacking between helices d′ and d′′, favoring a rigid duplex across the three-way junction. The κ-element folds into a stable GAAA-tetraloop motif and engages in A-minor interactions with helix d′. The addition of cobalt(III)hexammine reveals three distinct binding sites. The Mg(II)-promoted structural rearrangement and rigidification of the D1 core can be identified as the first micro-step of D1 folding.

Evolutionarily divergent spliceosomal snRNAs and a conserved non-coding RNA processing motif in Giardia lamblia

Non-coding RNAs (ncRNAs) have diverse essential biological functions in all organisms, and in eukaryotes, two such classes of ncRNAs are the small nucleolar (sno) and small nuclear (sn) RNAs. In this study, we have identified and characterized a collection of sno and snRNAs in Giardia lamblia, by exploiting our discovery of a conserved 12 nt RNA processing sequence motif found in the 3' end regions of a large number of G. lamblia ncRNA genes. RNA end mapping and other experiments indicate the motif serves to mediate ncRNA 3' end formation from mono- and di-cistronic RNA precursor transcripts. Remarkably, we find the motif is also utilized in the processing pathway of all four previously identified trans-spliced G. lamblia introns, revealing a common RNA processing pathway for ncRNAs and trans-spliced introns in this organism. Motif sequence conservation then allowed for the bioinformatic and experimental identification of additional G. lamblia ncRNAs, including new U1 and U6 spliceosomal snRNA candidates. The U6 snRNA candidate was then used as a tool to identity novel U2 and U4 snRNAs, based on predicted phylogenetically conserved snRNA-snRNA base-pairing interactions, from a set of previously identified G. lamblia ncRNAs without assigned function. The Giardia snRNAs retain the core features of spliceosomal snRNAs but are sufficiently evolutionarily divergent to explain the difficulties in their identification. Most intriguingly, all of these snRNAs show structural features diagnostic of U2-dependent/major and U12-dependent/minor spliceosomal snRNAs.

Molecular dynamics re-refinement of two different small RNA loop structures using the original NMR data suggest a common structure

Restrained molecular dynamics simulations are a robust, though perhaps underused, tool for the end-stage refinement of biomolecular structures. We demonstrate their utility-using modern simulation protocols, optimized force fields, and inclusion of explicit solvent and mobile counterions-by re-investigating the solution structures of two RNA hairpins that had previously been refined using conventional techniques. The structures, both domain 5 group II intron ribozymes from yeast ai5γ and Pylaiella littoralis, share a nearly identical primary sequence yet the published 3D structures appear quite different. Relatively long restrained MD simulations using the original NMR restraint data identified the presence of a small set of violated distance restraints in one structure and a possibly incorrect trapped bulge nucleotide conformation in the other structure. The removal of problematic distance restraints and the addition of a heating step yielded representative ensembles with very similar 3D structures and much lower pairwise RMSD values. Analysis of ion density during the restrained simulations helped to explain chemical shift perturbation data published previously. These results suggest that restrained MD simulations, with proper caution, can be used to "update" older structures or aid in the refinement of new structures that lack sufficient experimental data to produce a high quality result. Notable cautions include the need for sufficient sampling, awareness of potential force field bias (such as small angle deviations with the current AMBER force fields), and a proper balance between the various restraint weights.

Biomass production of site selective 13C/15N nucleotides using wild type and a transketolase E. coli mutant for labeling RNA for high resolution NMR

Characterization of the structure and dynamics of nucleic acids by NMR benefits significantly from position specifically labeled nucleotides. Here an E. coli strain deficient in the transketolase gene (tktA) and grown on glucose that is labeled at different carbon sites is shown to facilitate cost-effective and large scale production of useful nucleotides. These nucleotides are site specifically labeled in C1' and C5' with minimal scrambling within the ribose ring. To demonstrate the utility of this labeling approach, the new site-specific labeled and the uniformly labeled nucleotides were used to synthesize a 36-nt RNA containing the catalytically essential domain 5 (D5) of the brown algae group II intron self-splicing ribozyme. The D5 RNA was used in binding and relaxation studies probed by NMR spectroscopy. Key nucleotides in the D5 RNA that are implicated in binding Mg(2+) ions are well resolved. As a result, spectra obtained using selectively labeled nucleotides have higher signal-to-noise ratio compared to those obtained using uniformly labeled nucleotides. Thus, compared to the uniformly (13)C/(15)N-labeled nucleotides, these specifically labeled nucleotides eliminate the extensive (13)C-(13)C coupling within the nitrogenous base and ribose ring, give rise to less crowded and more resolved NMR spectra, and accurate relaxation rates without the need for constant-time or band-selective decoupled NMR experiments. These position selective labeled nucleotides should, therefore, find wide use in NMR analysis of biologically interesting RNA molecules.

Correlating SHAPE signatures with three-dimensional RNA structures

Selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) is a facile technique for quantitative analysis of RNA secondary structure. In general, low SHAPE signal values indicate Watson-Crick base-pairing, and high values indicate positions that are single-stranded within the RNA structure. However, the relationship of SHAPE signals to structural properties such as non-Watson-Crick base-pairing or stacking has thus far not been thoroughly investigated. Here, we present results of SHAPE experiments performed on several RNAs with published three-dimensional structures. This strategy allows us to analyze the results in terms of correlations between chemical reactivities and structural properties of the respective nucleotide, such as different types of base-pairing, stacking, and phosphate-backbone interactions. We find that the RNA SHAPE signal is strongly correlated with cis-Watson-Crick/Watson-Crick base-pairing and is to a remarkable degree not dependent on other structural properties with the exception of stacking. We subsequently generated probabilistic models that estimate the likelihood that a residue with a given SHAPE score participates in base-pairing. We show that several models that take SHAPE scores of adjacent residues into account perform better in predicting base-pairing compared with individual SHAPE scores. This underscores the context sensitivity of SHAPE and provides a framework for an improved interpretation of the response of RNA to chemical modification.

Group II introns: mobile ribozymes that invade DNA

Group II introns are mobile ribozymes that self-splice from precursor RNAs to yield excised intron lariat RNAs, which then invade new genomic DNA sites by reverse splicing. The introns encode a reverse transcriptase that stabilizes the catalytically active RNA structure for forward and reverse splicing, and afterwards converts the integrated intron RNA back into DNA. The characteristics of group II introns suggest that they or their close relatives were evolutionary ancestors of spliceosomal introns, the spliceosome, and retrotransposons in eukaryotes. Further, their ribozyme-based DNA integration mechanism enabled the development of group II introns into gene targeting vectors ("targetrons"), which have the unique feature of readily programmable DNA target specificity.

Linking the branchpoint helix to a newly found receptor allows lariat formation by a group II intron

Like spliceosomal introns, the ribozyme-containing group II introns are excised as branched, lariat structures: a 2'-5' bond is created between the first nucleotide of the intron and an adenosine in domain VI, a component which is missing from available crystal structures of the ribozyme. Comparative sequence analysis, modelling and nucleotide substitutions point to the existence, and probable location, of a specific RNA receptor for the section of domain VI that lies just distal to the branchpoint adenosine. By designing oligonucleotides that tether domain VI to this novel binding site, we have been able to specifically activate lariat formation in an engineered, defective group II ribozyme. The location of the newly identified receptor implies that prior to exon ligation, the distal part of domain VI undergoes a major translocation, which can now be brought under control by the system of anchoring oligonucleotides we have developed. Interestingly, these oligonucleotides, which link the branchpoint helix and the binding site for intron nucleotides 3-4, may be viewed as counterparts of U2-U6 helix III in the spliceosome.

Metal ions: supporting actors in the playbook of small ribozymes

Since the 1980s, several small RNA motifs capable of chemical catalysis have been discovered. These small ribozymes, composed of between approximately 40 and 200 nucleotides, have been found to play vital roles in the replication of subviral and viral pathogens, as well as in gene regulation in prokaryotes, and have recently been discovered in noncoding eukaryotic RNAs. All of the known natural small ribozymes - the hairpin, hammerhead, hepatitis delta virus, Varkud satellite, and glmS ribozymes--catalyze the same self-cleavage reaction as RNase A, resulting in two products, one bearing a 2'-3' cyclic phosphate and the other a 5'-hydroxyl group. Although originally thought to be obligate metalloenzymes like the group I and II self-splicing introns, the small ribozymes are now known to support catalysis in a wide variety of cations that appear to be only indirectly involved in catalysis. Nevertheless, under physiologic conditions, metal ions are essential for the proper folding and function of the small ribozymes, the most effective of these being magnesium. Metal ions contribute to catalysis in the small ribozymes primarily by stabilizing the catalytically active conformation, but in some cases also by activating RNA functional groups for catalysis, directly participating in catalytic acid-base chemistry, and perhaps by neutralizing the developing negative charge of the transition state. Although interactions between the small ribozymes and cations are relatively nonspecific, ribozyme activity is quite sensitive to the types and concentrations of metal ions present in solution, suggesting a close evolutionary relationship between cellular metal ion homeostasis and cation requirements of catalytic RNAs, and perhaps RNA in general.

Probing adenine rings and backbone linkages using base specific isotope-edited Raman spectroscopy: application to group II intron ribozyme domain V

Raman difference spectroscopy is used to probe the properties of a 36-nt RNA molecule, "D5", which lies at the heart of the catalytic apparatus in group II introns. For D5 that has all of its adenine residues labeled with (13)C and (15)N and utilizing Raman difference spectroscopy, we identify the conformationally sensitive -C-O-P-O-C- stretching modes of the unlabeled bonds adjacent to adenine bases, as well as the adenine ring modes themselves. The phosphodiester modes can be assigned to individual adenine residues based on earlier NMR data. The effect of Mg(2+) binding was explored by analyzing the Raman difference spectra for [D5 + Mg(2+)] minus [D5 no Mg(2+)], for D5 unlabeled, or D5 labeled with (13)C/(15)N-enriched adenine. In both sets of data we assign differential features to G ring modes perturbed by Mg(2+) binding at the N7 position. In the A-labeled spectra we attribute a Raman differential near 1450 cm(-1) and changes of intensity at 1296 cm(-1) to Mg binding at the N7 position of adenine bases. The A and G bases involved in Mg(2+) binding again can be identified using earlier NMR results. For the unlabeled D5, a change in the C-O-P-O-C stretch profile at 811 cm(-1) upon magnesium binding is due to a "tightening up" (in the sense of a more rigid molecule with less dynamic interchange among competing ribose conformers) of the D5 structure. For adenine-labeled D5, small changes in the adenine backbone bond signatures in the 810-830 cm(-1) region suggest that small conformational changes occur in the tetraloop and bulge regions upon binding of Mg(2+). The PO(2)(-) stretching vibration, near 1100 cm(-1), from the nonbridging phosphate groups, probes the effect of Mg(2+)-hydrate inner-sphere interactions that cause an upshift. In turn, the upshift is modulated by the presence of monovalent cations since in the presence of Na(+) and Li(+) the upshift is 23 +/- 2 cm(-1) while in the presence of K(+) and Cs(+) it is 13 +/- 3 cm(-1), a finding that correlates with the differences in hydration radii. These subtle differences in electrostatic interactions may be related to observed variations in catalytic activity. For a reconstructed ribozyme comprising domains 1-3 (D123) connected in cis plus domain 5 (D5) supplied in trans, cleavage of spliced exon substrates in the presence of magnesium and K(+) or Cs(+) is more efficient than that in the presence of magnesium with Na(+) or Li(+).

High fidelity of base pairing by 2-selenothymidine in DNA

The base pairs are the contributors to the sequence-dependent recognition of nucleic acids, genetic information storage, and high fidelity of DNA polymerase replication. However, the wobble base pairing, where T pairs with G instead of A, reduces specific base-pairing recognition and compromises the high fidelity of the enzymatic polymerization. Via the selenium atomic probing at the 2-position of thymidine, we have investigated the wobble discrimination by manipulating the steric and electronic effects at the 2-exo position, providing a unique chemical strategy to enhance the base pair specificity. We report here the first synthesis of the novel 2-Se-thymidine ((Se)T) derivative, its phosphoramidite, and the Se-DNAs. Our biophysical and structural studies of the 2-Se-T DNAs reveal that the bulky 2-Se atom with a weak hydrogen-bonding ability can largely increase mismatch discriminations (including T/G wobble and T/C mismatched base pairs) while maintaining the (Se)T/A virtually identical to the native T/A base pair. The 2-Se atom bulkiness and the electronic effect are probably the main factors responsible for the discrimination against the formation of the wobble (Se)T/G base pair. Our investigations provide a potential novel tool to investigate the specific recognition of base pairs, which is the basis of high fidelity during replication, transcription, and translation. Furthermore, this Se-atom-specific substitution and probing are useful for X-ray crystal structure and function studies of nucleic acids.

A structural analysis of the group II intron active site and implications for the spliceosome

Group II introns are self-splicing, mobile genetic elements that have fundamentally influenced the organization of terrestrial genomes. These large ribozymes remain important for gene expression in almost all forms of bacteria and eukaryotes and they are believed to share a common ancestry with the eukaryotic spliceosome that is required for processing all nuclear pre-mRNAs. The three-dimensional structure of a group IIC intron was recently determined by X-ray crystallography, making it possible to visualize the active site and the elaborate network of tertiary interactions that stabilize the molecule. Here we describe the molecular features of the active site in detail and evaluate their correspondence with prior biochemical, genetic, and phylogenetic analyses on group II introns. In addition, we evaluate the structural significance of RNA motifs within the intron core, such as the major-groove triple helix and the domain 5 bulge. Having combined what is known about the group II intron core, we then compare it with known structural features of U6 snRNA in the eukaryotic spliceosome. This analysis leads to a set of predictions for the molecular structure of the spliceosomal active site.

Knowledge-based instantiation of full atomic detail into coarse-grain RNA 3D structural models

Motivation: The recent development of methods for modeling RNA 3D structures using coarse-grain approaches creates a need to bridge low- and high-resolution modeling methods. Although they contain topological information, coarse-grain models lack atomic detail, which limits their utility for some applications.

Results: We have developed a method for adding full atomic detail to coarse-grain models of RNA 3D structures. Our method [Coarse to Atomic (C2A)] uses geometries observed in known RNA crystal structures. Our method rebuilds full atomic detail from ideal coarse-grain backbones taken from crystal structures to within 1.87–3.31 Å RMSD of the full atomic crystal structure. When starting from coarse-grain models generated by the modeling tool NAST, our method builds full atomic structures that are within 1.00 Å RMSD of the starting structure. The resulting full atomic structures can be used as starting points for higher resolution modeling, thus bridging high- and low-resolution approaches to modeling RNA 3D structure.

Availability: Code for the C2A method, as well as the examples discussed in this article, are freely available at www.simtk.org/home/c2a.

Contact:russ.altman@stanford.edu

Spliceosome structure: piece by piece

Processing of pre-mRNAs by RNA splicing is an essential step in the maturation of protein coding RNAs in eukaryotes. Structural studies of the cellular splicing machinery, the spliceosome, are a major challenge in structural biology due to the size and complexity of the splicing ensemble. Specifically, the structural details of splice site recognition and the architecture of the spliceosome active site are poorly understood. X-ray and NMR techniques have been successfully used to address these questions defining the structure of individual domains, isolated splicing proteins, spliceosomal RNA fragments and recently the U1 snRNP multiprotein.RNA complex. These results combined with extant biochemical and genetic data have yielded important insights as well as posing fresh questions with respect to the regulation and mechanism of this critical gene regulatory process.

Decision-making in structure solution using Bayesian estimates of map quality: the PHENIX AutoSol wizard

Ten measures of experimental electron-density-map quality are examined and the skewness of electron density is found to be the best indicator of actual map quality. A Bayesian approach to estimating map quality is developed and used in the PHENIX AutoSol wizard to make decisions during automated structure solution.

Estimates of the quality of experimental maps are important in many stages of structure determination of macromolecules. Map quality is defined here as the correlation between a map and the corresponding map obtained using phases from the final refined model. Here, ten different measures of experimental map quality were examined using a set of 1359 maps calculated by re-analysis of 246 solved MAD, SAD and MIR data sets. A simple Bayesian approach to estimation of map quality from one or more measures is presented. It was found that a Bayesian estimator based on the skewness of the density values in an electron-density map is the most accurate of the ten individual Bayesian estimators of map quality examined, with a correlation between estimated and actual map quality of 0.90. A combination of the skewness of electron density with the local correlation of r.m.s. density gives a further improvement in estimating map quality, with an overall correlation coefficient of 0.92. The PHENIX AutoSol wizard carries out automated structure solution based on any combination of SAD, MAD, SIR or MIR data sets. The wizard is based on tools from the PHENIX package and uses the Bayesian estimates of map quality described here to choose the highest quality solutions after experimental phasing.

The ionic environment determines ribozyme cleavage rate by modulation of nucleobase pK a

Several small ribozymes employ general acid-base catalysis as a mechanism to enhance site-specific RNA cleavage, even though the functional groups on the ribonucleoside building blocks of RNA have pK (a) values far removed from physiological pH. The rate of the cleavage reaction is strongly affected by the identity of the metal cation present in the reaction solution; however, the mechanism(s) by which different cations contribute to rate enhancement has not been determined. Using the Neurospora VS ribozyme, we provide evidence that different cations confer particular shifts in the apparent pK (a) values of the catalytic nucleobases, which in turn determines the fraction of RNA in the protonation state competent for general acid-base catalysis at a given pH, which determines the observed rate of the cleavage reaction. Despite large differences in observed rates of cleavage in different cations, mathematical models of general acid-base catalysis indicate that k (1), the intrinsic rate of the bond-breaking step, is essentially constant irrespective of the identity of the cation(s) in the reaction solution. Thus, in contrast to models that invoke unique roles for metal ions in ribozyme chemical mechanisms, we find that most, and possibly all, of the ion-specific rate enhancement in the VS ribozyme can be explained solely by the effect of the ions on nucleobase pK (a). The inference that k (1) is essentially constant suggests a resolution of the problem of kinetic ambiguity in favor of a model in which the lower pK (a) is that of the general acid and the higher pK (a) is that of the general base.

A glimpse into the active site of a group II intron and maybe the spliceosome, too

The X-ray crystal structure of an excised group II self-splicing intron was recently solved by the Pyle group. Here we review some of the notable features of this structure and what they may tell us about the catalytic active site of the group II ribozyme and potentially the spliceosome. The new structure validates the central role of domain V in both the structure and catalytic function of the ribozyme and resolves several outstanding puzzles raised by previous biochemical, genetic and structural studies. While lacking both exons as well as the cleavage sites and nucleophiles, the structure reveals how a network of tertiary interactions can position two divalent metal ions in a configuration that is ideal for catalysis.

Specificity of Mg2+ binding at the Group II intron branch site

Metal ions play a crucial role in the conformation and splicing activity of Group II introns. Results from 2-aminopurine fluorescence and solution NMR studies suggest that metal ion binding within the branch site region of native D6 of the Group II intron is specific for alkaline earth metal ions and involves inner sphere coordination. Although Mg(2+) and Ca(2+) still bind to a mutant stem loop sequence from which the internal loop had been deleted, ion binding to the mutant RNA results in decreased, rather than increased, exposure of the branch site residue to solvent. These data further support the role of the internal loop in defining branch site conformation of the Group II intron. The specific bound Mg(2+) may play a bivalent role: facilitates the extrahelical conformation of the branch site and has the potential to act as a Lewis acid during splicing.

X-ray structures of U2 snRNA-branchpoint duplexes containing conserved pseudouridines

A pseudouridine-modified region of the U2 small nuclear (sn)RNA anneals with the intronic branchpoint sequence and positions a bulged adenosine to serve as the nucleophile in the first chemical step of pre-mRNA splicing. We have determined three X-ray structures of RNA oligonucleotides containing the pseudouridylated U2 snRNA and the branchpoint consensus sequences. The expected adenosine branchpoint is extrahelical in a 1.65 A resolution structure containing the mammalian consensus sequence variant and in a 2.10 A resolution structure containing a shortened Saccharomyces cerevisiae consensus sequence. The adenosine adjacent to the expected branchpoint is extrahelical in a third structure, which contains the intact yeast consensus sequence at 1.57 A resolution. The hydration and base stacking interactions mediated by the U2 snRNA pseudouridines correlate with the identity of the unpaired adenosine. The expected adenosine bulge is associated with a well-stacked pseudouridine, which is linked via an ordered water molecule to a neighboring nucleotide. In contrast, the bulge of the adjacent adenosine shifts the base stacking and disrupts the water-mediated interactions of the pseudouridine. These structural differences may contribute to the ability of the pseudouridine modification to promote the bulged conformation of the branch site adenosine and to enhance catalysis by snRNAs. Furthermore, iodide binding sites are identified adjacent to the unconventional bulged adenosine, and the structure of the mammalian consensus sequence variant provides a high-resolution view of a hydrated magnesium ion bound in a similar manner to a divalent cation binding site of the group II intron.

Crystal structure of a self-spliced group II intron

Group II introns are self-splicing ribozymes that catalyze their own excision from precursor transcripts and insertion into new genetic locations. Here we report the crystal structure of an intact, self-spliced group II intron from Oceanobacillus iheyensis at 3.1 angstrom resolution. An extensive network of tertiary interactions facilitates the ordered packing of intron subdomains around a ribozyme core that includes catalytic domain V. The bulge of domain V adopts an unusual helical structure that is located adjacent to a major groove triple helix (catalytic triplex). The bulge and catalytic triplex jointly coordinate two divalent metal ions in a configuration that is consistent with a two-metal ion mechanism for catalysis. Structural and functional analogies support the hypothesis that group II introns and the spliceosome share a common ancestor.

2'-Methylseleno-modified oligoribonucleotides for X-ray crystallography synthesized by the ACE RNA solid-phase approach

Site-specifically modified 2'-methylseleno RNA represents a valuable derivative for phasing of X-ray crystallographic data. Several successful applications in three-dimensional structure determination of nucleic acids, such as the Diels-Alder ribozyme, have relied on this modification. Here, we introduce synthetic routes to 2'-methylseleno phosphoramidite building blocks of all four standard nucleosides, adenosine, cytidine, guanosine and uridine, that are tailored for 2'-O-bis(acetoxyethoxy)methyl (ACE) RNA solid-phase synthesis. We additionally report on their incorporation into oligoribonucleotides including deprotection and purification. The methodological expansion of 2'-methylseleno labeling via ACE RNA chemistry is a major step to make Se-RNA generally accessible and to receive broad dissemination of the Se-approach for crystallographic studies on RNA. Thus far, preparation of 2'-methylseleno-modified oligoribonucleotides has been restricted to the 2'-O-[(triisopropylsilyl)oxy]methyl (TOM) and 2'-O-tert-butyldimethylsilyl (TBDMS) RNA synthesis methods.

Solution structure of domain 6 from a self-splicing group II intron ribozyme: a Mg(2+) binding site is located close to the stacked branch adenosine

Group II intron self-splicing is essential for the correct expression of organellar genes in plants, fungi, and yeast, as well as of bacterial genes. Self-excision of these autocatalytic introns from the primary RNA transcript is achieved in a two-step mechanism that is apparently analogous to that of the eukaryotic spliceosome. The 2'-OH of a conserved adenosine (the branch point) located within domain 6 (D6) acts as the nucleophile in the first step of splicing. Despite the biological importance of group II introns, little is known about their structural organization and usage of metal ions in catalysis. Here we report the first solution structure of a catalytically active D6 construct encompassing the branch point and the neighboring helical regions from the mitochondrial yeast intron ai5gamma. The branch adenosine is the single unpaired nucleotide, and, in contrast to the spliceosomal branch site, resides within the helix, being partially stacked between two flanking GU wobble pairs. We identified a novel prominent Mg(2+) binding site in the major groove of the branch site. Importantly, Mg(2+) addition does not impair the stacking of the branch adenosine, rather it strengthens the interaction with the flanking uridines, as shown by NMR and fluorescence studies. This means that domain 6 presents the branch adenosine in a stacked fashion to the core of group II introns upon folding to the active conformation.

Three essential and conserved regions of the group II intron are proximal to the 5'-splice site

Despite the central role of group II introns in eukaryotic gene expression and their importance as biophysical and evolutionary model systems, group II intron tertiary structure is not well understood. In order to characterize the architectural organization of intron ai5gamma, we incorporated the photoreactive nucleotides s(4)U and s(6)dG at specific locations within the intron core and monitored the formation of cross-links in folded complexes. The resulting data reveal the locations for many of the most conserved, catalytically important regions of the intron (i.e., the J2/3 linker region, the IC1(i-ii) bulge in domain 1, the bulge of D5, and the 5'-splice site), showing that all of these elements are closely colocalized. In addition, we show by nucleotide analog interference mapping (NAIM) that a specific functional group in J2/3 plays a role in first-step catalysis, which is consistent with its apparent proximity to other first-step components. These results extend our understanding of active-site architecture during the first step of group II intron self-splicing and they provide a structural basis for spliceosomal comparison.

Automated de novo prediction of native-like RNA tertiary structures

RNA tertiary structure prediction has been based almost entirely on base-pairing constraints derived from phylogenetic covariation analysis. We describe here a complementary approach, inspired by the Rosetta low-resolution protein structure prediction method, that seeks the lowest energy tertiary structure for a given RNA sequence without using evolutionary information. In a benchmark test of 20 RNA sequences with known structure and lengths of approximately 30 nt, the new method reproduces better than 90% of Watson-Crick base pairs, comparable with the accuracy of secondary structure prediction methods. In more than half the cases, at least one of the top five models agrees with the native structure to better than 4 A rmsd over the backbone. Most importantly, the method recapitulates more than one-third of non-Watson-Crick base pairs seen in the native structures. Tandem stacks of "sheared" base pairs, base triplets, and pseudoknots are among the noncanonical features reproduced in the models. In the cases in which none of the top five models were native-like, higher energy conformations similar to the native structures are still sampled frequently but not assigned low energies. These results suggest that modest improvements in the energy function, together with the incorporation of information from phylogenetic covariance, may allow confident and accurate structure prediction for larger and more complex RNA chains.

Structural studies of the purine and SAM binding riboswitches

Riboswitches are recently discovered genetic regulatory elements found in the 5'-untranslated regions of bacterial mRNAs that act through their ability to specifically bind small-molecule metabolites. Binding of the ligand to the aptamer domain of the riboswitch is communicated to a second domain, the expression platform, which directs transcription or translation of the mRNA. To understand this process on a molecular level, structures of three of these riboswitches bound to their cognate ligands have been solved by X-ray crystallography: the purine, thiamine pyrophosphate (TPP), and S-adenosylmethionine (SAM-I) binding aptamer domains. These studies have uncovered three common themes between the otherwise different molecules. First, the natural RNA aptamers recognize directly or indirectly almost every feature of their ligand to achieve extraordinary specificity. Second, all of these RNAs use a complex tertiary architecture to establish the binding pocket. Finally, in each case, ligand binding serves to stabilize a helix that communicates the binding event to the expression platform. Here, we discuss these properties of riboswitches in the context of the purine and SAM-I riboswitches.

A Second Divalent Metal Ion in the Group II Intron Reaction Center

Group II introns are mobile genetic elements that have been implicated as agents of genetic diversity, and serve as important model systems for investigating RNA catalysis and pre-mRNA splicing. In the absence of an atomic-resolution structure of the intron, detailed understanding of its catalytic mechanism has remained elusive. Previous identification of a divalent metal ion stabilizing the leaving group in both splicing steps suggested that the group II intron may employ a "two-metal ion" mechanism, a catalytic strategy used by a number of protein phosphoester transfer enzymes. Using metal rescue experiments, we now reveal the presence of a second metal ion required for nucleophile activation in the exon-ligation step of group II intron splicing. Coupled with biochemical and structural evidence of at least two metal ions at the group I intron reaction center, these results suggest a mechanistic paradigm for describing catalysis by large ribozymes.