Crystal structure of a hairpin ribozyme-inhibitor complex with implications for catalysis

Facile synthesis of nucleoside 5'-(α-P-seleno)-triphosphates and phosphoroselenoate RNA transcription

Phosphoroselenoate RNA (PSe-RNA) is nuclease resistant and has great potentials in X-ray crystal structure and function studies of noncoding RNAs and protein-RNA interactions. In order to conveniently synthesize PSe-RNA via transcription, we have developed a one-pot synthetic method for the nucleoside 5'-(α-P-seleno)-triphosphates (NTPαSe) analogs without protecting any functionality of the ribonucleosides. The NTPαSe diastereomers have been purified, fully characterized, and incorporated into RNAs by T7 RNA polymerase. The transcribed RNAs are diastereomerically pure, and the Se-derivatized ribozymes are generally active. Furthermore, we have established an affinity purification strategy by using immobilized boronate to conveniently purify NTPαSe analogs. Though the affinity-purified NTPαSe analogs are diastereomeric mixtures, they can be directly used in transcription without a significant impact on the transcription efficiency. Moreover, we found that the PSe-nucleotide is stable during polyacrylamide gel purification, indicating that the PSe-RNAs can be purified straightforwardly for crystal structural and functional studies.

New insights into the fundamental role of topological constraints as a determinant of two-way junction conformation

Recent studies have shown that topological constraints encoded at the RNA secondary structure level involving basic steric and stereochemical forces can significantly restrict the orientations sampled by helices across two-way RNA junctions. Here, we formulate these topological constraints in greater quantitative detail and use this topological framework to rationalize long-standing but poorly understood observations regarding the basic behavior of RNA two-way junctions. Notably, we show that the asymmetric nature of the A-form helix and the finite length of a bulge provide a physical basis for the experimentally observed directionality and bulge-length amplitude dependence of bulge induced inter-helical bends. We also find that the topologically allowed space can be modulated by variations in sequence, particularly with the addition of non-canonical GU base pairs at the junction, and, surprisingly, by the length of the 5′ and 3′ helices. A survey of two-way RNA junctions in the protein data bank confirms that junction residues have a strong preference to adopt looped-in, non-canonically base-paired conformations, providing a route for extending our bulge-directed framework to internal loop motifs and implying a simplified link between secondary and tertiary structure. Finally, our results uncover a new simple mechanism for coupling junction-induced topological constraints with tertiary interactions.

Thermodynamics and kinetics of the hairpin ribozyme from atomistic folding/unfolding simulations

We report a set of atomistic folding/unfolding simulations for the hairpin ribozyme using a Monte Carlo algorithm. The hairpin ribozyme folds in solution and catalyzes self-cleavage or ligation via a specific two-domain structure. The minimal active ribozyme has been studied extensively, showing stabilization of the active structure by cations and dynamic motion of the active structure. Here, we introduce a simple model of tertiary-structure formation that leads to a phase diagram for the RNA as a function of temperature and tertiary-structure strength. We then employ this model to capture many folding/unfolding events and to examine the transition-state ensemble (TSE) of the RNA during folding to its active "docked" conformation. The TSE is compact but with few tertiary interactions formed, in agreement with single-molecule dynamics experiments. To compare with experimental kinetic parameters, we introduce a novel method to benchmark Monte Carlo kinetic parameters to docking/undocking rates collected over many single molecular trajectories. We find that topology alone, as encoded in a biased potential that discriminates between secondary and tertiary interactions, is sufficient to predict the thermodynamic behavior and kinetic folding pathway of the hairpin ribozyme. This method should be useful in predicting folding transition states for many natural or man-made RNA tertiary structures.

Modeling the RNA 2'OH activation: possible roles of metal ion and nucleobase as catalysts in self-cleaving ribozymes

The RNA 2'OH activation as taking place in the first chemical step of self-cleaving ribozymes is studied theoretically by DFT and MP2 methods using a continuum solvation model (CPCM). The reaction of proton transfer is studied in the presence of two kinds of catalysts: a fully hydrated metal ion (Mg(2+)) or partially hydrated nucleobase (guanine), taken separately or together leading to three different modes of activation. The metal ion is either directly bound (inner-sphere) or indirectly bound (outer-sphere) to the 2'OH group and a hydroxide ion acts as a general or specific base; the nucleobase is taken in anionic or in neutral enol-tautomeric forms playing itself the role of general base. The presence of a close metal ion (outer-sphere) lowers the pK(a) value of the 2'OH group by several log units in both metal-ion and nuleobase catalysis. The direct metal coordination to the 2'OH group (inner-sphere) further stabilizes the developing negative charge on the nucleophile. The switching from the inner-sphere to the outer-sphere coordination appears to be driven by the energy cost for reorganizing the first coordination shell rather than by the electrostatic repulsion between the ligands. The metal-ion catalysis is more effective with a specific base in the dianionic mechanism. On the other hand, the nucleobase catalysis is more effective in the monoanionic mechanism and in the presence of a metal ion acting as a cofactor through nonspecific electrostatic interactions. The results establish a baseline to study the possible roles of metal and nucleobase catalysts and their environment in more realistic models for self-cleaving ribozymes.

glmS Riboswitch binding to the glucosamine-6-phosphate α-anomer shifts the pKa toward neutrality

The glmS riboswitch regulates gene expression through a self-cleavage activity. The reaction is catalyzed with the assistance of the metabolite cofactor glucosamine-6-phosphate (GlcN6P), whose amino group is proposed to serve as the general acid during the reaction. This reaction is pH-dependent with a pK(a) that is lower than the observed pK(a) for the amine of GlcN6P in solution. GlcN6P, like other pyranose sugars, undergoes spontaneous and rapid interconversion between the α and β anomers at the C1 position. Here we demonstrate by NMR that the Bacillus anthracis glmS riboswitch selectively binds the α-anomer of GlcN6P with a maximum binding affinity of 0.36 mM and that binding is pH-dependent. We also report that the anomeric ratio between α and β is pH-dependent and the pK(a)s of the two amines differ by 0.5 pH units, α being the higher of the two (pK(a)=8.3). The pH dependence of binding reveals a pK(a) of 6.7, suggesting that the glmS RNA reduces the pK(a) of the GlcN6P amine by 1.6 units in the ground state. We reevaluated previously obtained kinetic data and found the reaction pK(a) is 6.9, within error of the binding data. The data support a model where the reaction pK(a) corresponds to that of the GlcN6P amine. This observation has broader relevance for considering how the microenvironment of an RNA, despite its anionic character, can reduce the pK(a)s of functional groups for use in catalysis.

Catalysis by the nucleolytic ribozymes

The nucleolytic ribozymes use general acid-base catalysis to contribute significantly to their rate enhancement. The VS (Varkud satellite) ribozyme uses a guanine and an adenine nucleobase as general base and acid respectively in the cleavage reaction. The hairpin ribozyme is probably closely similar, while the remaining nucleolytic ribozymes provide some interesting contrasts.

Nucleic acid X-ray crystallography via direct selenium derivatization

X-ray crystallography has proven to be an essential tool for structural studies of bio-macromolecules at the atomic level. There are two major bottle-neck problems in the macromolecular crystal structure determination: phasing and crystallization. Although the selenium derivatization is routinely used for solving novel protein structures through the MAD phasing technique, the phase problem is still a critical issue in nucleic acid crystallography. The background and current progress of using direct selenium-derivatization of nucleic acids (SeNA) to solve the phase problem and to facilitate nucleic acid crystallization for X-ray crystallography are summarized in this tutorial review.

The shape-shifting quasispecies of RNA: one sequence, many functional folds

E Unus pluribum, or "Of One, Many", may be at the root of decoding the RNA sequence-structure-function relationship. RNAs embody the large majority of genes in higher eukaryotes and fold in a sequence-directed fashion into three-dimensional structures that perform functions conserved across all cellular life forms, ranging from regulating to executing gene expression. While it is the most important determinant of the RNA structure, the nucleotide sequence is generally not sufficient to specify a unique set of secondary and tertiary interactions due to the highly frustrated nature of RNA folding. This frustration results in folding heterogeneity, a common phenomenon wherein a chemically homogeneous population of RNA molecules folds into multiple stable structures. Often, these alternative conformations constitute misfolds, lacking the biological activity of the natively folded RNA. Intriguingly, a number of RNAs have recently been described as capable of adopting multiple distinct conformations that all perform, or contribute to, the same function. Characteristically, these conformations interconvert slowly on the experimental timescale, suggesting that they should be regarded as distinct native states. We discuss how rugged folding free energy landscapes give rise to multiple native states in the Tetrahymena Group I intron ribozyme, hairpin ribozyme, sarcin-ricin loop, ribosome, and an in vitro selected aptamer. We further describe the varying degrees to which folding heterogeneity impacts function in these RNAs, and compare and contrast this impact with that of heterogeneities found in protein folding. Embracing that one sequence can give rise to multiple native folds, we hypothesize that this phenomenon imparts adaptive advantages on any functionally evolving RNA quasispecies.

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.

Physics-based de novo prediction of RNA 3D structures

Current experiments on structural determination cannot keep up the pace with the steadily emerging RNA sequences and new functions. This underscores the request for an accurate model for RNA three-dimensional (3D) structural prediction. Although considerable progress has been made in mechanistic studies, accurate prediction for RNA tertiary folding from sequence remains an unsolved problem. The first and most important requirement for the prediction of RNA structure from physical principles is an accurate free energy model. A recently developed three-vector virtual bond-based RNA folding model ("Vfold") has allowed us to compute the chain entropy and predict folding free energies and structures for RNA secondary structures and simple pseudoknots. Here we develop a free energy-based method to predict larger more complex RNA tertiary folds. The approach is based on a multiscaling strategy: from the nucleotide sequence, we predict the two-dimensional (2D) structures (defined by the base pairs and tertiary contacts); based on the 2D structure, we construct a 3D scaffold; with the 3D scaffold as the initial state, we combine AMBER energy minimization and PDB-based fragment search to predict the all-atom structure. A key advantage of the approach is the statistical mechanical calculation for the conformational entropy of RNA structures, including those with cross-linked loops. Benchmark tests show that the model leads to significant improvements in RNA 3D structure prediction.

Emerging structural themes in large RNA molecules

Extensive networks of tertiary interactions give rise to unique, highly organized domain architectures that characterize the three-dimensional structure of large RNA molecules. Formed by stacked layers of a near-planar arrangement of contiguous coaxial helices, large RNA molecules are relatively flat in overall shape. The functional core of these molecules is stabilized by a diverse set of tertiary interaction motifs that often bring together distant regions of conserved nucleotides. Although homologous RNAs from different organisms can be structurally diverse, they adopt a structurally conserved functional core that includes preassembled active and/or substrate binding sites. These findings broaden our understanding of RNA folding and tertiary structure stabilization, illustrating how large, complex RNAs assemble into unique structures to perform recognition and catalysis.

The pH dependence of hairpin ribozyme catalysis reflects ionization of an active site adenine

Understanding how self-cleaving ribozymes mediate catalysis is crucial in light of compelling evidence that human and bacterial gene expression can be regulated through RNA self-cleavage. The hairpin ribozyme catalyzes reversible phosphodiester bond cleavage through a mechanism that does not require divalent metal cations. Previous structural and biochemical evidence implicated the amidine group of an active site adenosine, A38, in a pH-dependent step in catalysis. We developed a way to determine microscopic pK(a) values in active ribozymes based on the pH-dependent fluorescence of 8-azaadenosine (8azaA). We compared the microscopic pK(a) for ionization of 8azaA at position 38 with the apparent pK(a) for the self-cleavage reaction in a fully functional hairpin ribozyme with a unique 8azaA at position 38. Microscopic and apparent pK(a) values were virtually the same, evidence that A38 protonation accounts for the decrease in catalytic activity with decreasing pH. These results implicate the neutral unprotonated form of A38 in a transition state that involves formation of the 5'-oxygen-phosphorus bond.

Solution structure of the K-turn and Specifier Loop domains from the Bacillus subtilis tyrS T-box leader RNA

In Gram-positive bacteria, the RNA transcripts of many amino acid biosynthetic and aminoacyl tRNA synthetase genes contain 5' untranslated regions, or leader RNAs, that function as riboswitches. These T-box riboswitches bind cognate tRNA molecules and regulate gene expression by a transcription attenuation mechanism. The Specifier Loop domain of the leader RNA contains nucleotides that pair with nucleotides in the tRNA anticodon loop and is flanked on one side by a kink-turn (K-turn), or GA, sequence motif. We have determined the solution NMR structure of the K-turn sequence element within the context of the Specifier Loop domain. The K-turn sequence motif has several noncanonical base pairs typical of K-turn structures but adopts an extended conformation. The Specifier Loop domain contains a loop E structural motif, and the single-strand Specifier nucleotides stack with their Watson-Crick edges displaced toward the minor groove. Mg(2+) leads to a significant bending of the helix axis at the base of the Specifier Loop domain, but does not alter the K-turn. Isothermal titration calorimetry indicates that the K-turn sequence causes a small enhancement of the interaction between the tRNA anticodon arm and the Specifier Loop domain. One possibility is that the K-turn structure is formed and stabilized when tRNA binds the T-box riboswitch and interacts with Stem I and the antiterminator helix. This motif in turn anchors the orientation of Stem I relative to the 3' half of the leader RNA, further stabilizing the tRNA-T box complex.

Predicting and modeling RNA architecture

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

Do the hairpin and VS ribozymes share a common catalytic mechanism based on general acid-base catalysis? A critical assessment of available experimental data

The active centers of the hairpin and VS ribozymes are both generated by the interaction of two internal loops, and both ribozymes use guanine and adenine nucleobases to accelerate cleavage and ligation reactions. The centers are topologically equivalent and the relative positioning of key elements the same. There is good evidence that the cleavage reaction of the VS ribozyme is catalyzed by the guanine (G638) acting as general base and the adenine (A756) as general acid. We now critically evaluate the experimental mechanistic evidence for the hairpin ribozyme. We conclude that all the available data are fully consistent with a major contribution to catalysis by general acid-base catalysis involving the adenine (A38) and guanine (G8). It appears that the two ribozymes are mechanistically equivalent.

NMR structure of the A730 loop of the Neurospora VS ribozyme: insights into the formation of the active site

The Neurospora VS ribozyme is a small nucleolytic ribozyme with unique primary, secondary and global tertiary structures, which displays mechanistic similarities to the hairpin ribozyme. Here, we determined the high-resolution NMR structure of a stem–loop VI fragment containing the A730 internal loop, which forms part of the active site. In the presence of magnesium ions, the A730 loop adopts a structure that is consistent with existing biochemical data and most likely reflects its conformation in the VS ribozyme prior to docking with the cleavage site internal loop. Interestingly, the A730 loop adopts an S-turn motif that is also present in loop B within the hairpin ribozyme active site. The S-turn appears necessary to expose the Watson–Crick edge of a catalytically important residue (A756) so that it can fulfill its role in catalysis. The A730 loop and the cleavage site loop of the VS ribozyme display structural similarities to internal loops found in the active site of the hairpin ribozyme. These similarities provided a rationale to build a model of the VS ribozyme active site based on the crystal structure of the hairpin ribozyme.

Protonation of base pairs in RNA: context analysis and quantum chemical investigations of their geometries and stabilities

Base pairs involving protonated nucleobases play important roles in mediating global macromolecular conformational changes and in facilitation of catalysis in a variety of functional RNA molecules. Here we present our attempts at understanding the role of such base pairs by detecting possible protonated base pairs in the available RNA crystal structures using BPFind software, in their specific structural contexts, and by the characterization of their geometries, interaction energies, and stabilities using advanced quantum chemical computations. We report occurrences of 18 distinct protonated base pair combinations from a representative data set of RNA crystal structures and propose a theoretical model for one putative base pair combination. Optimization of base pair geometries was carried out at the B3LYP/cc-pVTZ level, and the BSSE corrected interaction energies were calculated at the MP2/aug-cc-pVDZ level of theory. The geometries for each of the base pairs were characterized in terms of H-bonding patterns observed, rmsd values observed on optimization, and base pair geometrical parameters. In addition, the intermolecular interaction in these complexes was also analyzed using Morokuma energy decomposition. The gas phase interaction energies of the base pairs range from -24 to -49 kcal/mol and reveal the dominance of Hartree-Fock component of interaction energy constituting 73% to 98% of the total interaction energy values. On the basis of our combined bioinformatics and quantum chemical analysis of different protonated base pairs, we suggest resolution of structural ambiguities and correlate their geometric and energetic features with their structural and functional roles. In addition, we also examine the suitability of specific base pairs as key elements in molecular switches and as nucleators for higher order structures such as base triplets and quartets.

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.

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.

Prediction of a new class of RNA recognition motif

The observation that activation domains (AD) of procarboxypeptidases are rather long compared to the pro-regions of other zymogens raises the possibility that they could play additional roles apart from precluding enzymatic activity within the proenzyme and helping in its folding process. In the present work, we compared the overall pro-domain tertiary structure with several proteins belonging to the same fold in the structural classification of proteins (SCOP) database by using structure and sequence comparisons. The best score obtained was between the activation domain of human procarboxypeptidase A4 (ADA4h) and the human U1A protein from the U1 snRNP. Structural alignment revealed the existence of RNP1- and RNP2-related sequences in ADA4h. After modeling ADA4h on U1A, the new structure was used to extract a new sequence pattern characteristic for important residues at key positions. The new sequence pattern allowed scanning protein sequences to predict the RNA-binding function for 32 sequences undetected by PFAM. Unspecific RNA electrophoretic mobility shift assays experimentally supported the prediction that ADA4h binds an RNA motif similar to the U1A binding-motif of stem-loop II of U1 small nuclear RNA. The experiments carried out with ADA4h in the present work suggest the sharing of a common ancestor with other RNA recognition motifs. However, the fact that key residues preventing activity within the proenzyme are also key residues for RNA binding might have induced the activation domains of procarboxypeptidases to evolve from the canonical RNP1 and RNP2 sequences.

The glmS ribozyme: use of a small molecule coenzyme by a gene-regulatory RNA

The glmS ribozyme is the first known example of a natural ribozyme that has evolved to require binding of an exogenous small molecule for activity. In Gram-positive bacteria, this RNA domain is part of the messenger RNA (mRNA) encoding the essential enzyme that synthesizes glucosamine-6-phosphate (GlcN6P). When present at physiologic concentration, this small molecule binds to the glmS ribozyme and uncovers a latent self-cleavage activity that ultimately leads to degradation of the mRNA. Biochemical and structural studies reveal that the RNA adopts a rigid fold stabilized by three pseudoknots and the packing of a peripheral domain against the ribozyme core. GlcN6P binding to this pre-organized RNA does not induce conformational changes; rather, the small molecule functions as a coenzyme, providing a catalytically essential amine group to the active site. The ribozyme is not a passive player, however. Active site functional groups are essential for catalysis, even in the presence of GlcN6P. In addition to being a superb experimental system with which to analyze how RNA catalysts can exploit small molecule coenzymes to broaden their chemical versatility, the presence of the glmS ribozyme in numerous pathogenic bacteria make this RNA an attractive target for the development of new antibiotics and antibacterial strategies.

Taming free energy landscapes with RNA chaperones

Many non-coding RNAs fold into complex three-dimensional structures, yet the self-assembly of RNA structure is hampered by mispairing, weak tertiary interactions, electrostatic barriers, and the frequent requirement that the 5' and 3' ends of the transcript interact. This rugged free energy landscape for RNA folding means that some RNA molecules in a population rapidly form their native structure, while many others become kinetically trapped in misfolded conformations. Transient binding of RNA chaperone proteins destabilize misfolded intermediates and lower the transition states between conformations, producing a smoother landscape that increases the rate of folding and the probability that a molecule will find the native structure. DEAD-box proteins couple the chemical potential of ATP hydrolysis with repetitive cycles of RNA binding and release, expanding the range of conditions under which they can refold RNA structures.

Structural basis for 16S ribosomal RNA cleavage by the cytotoxic domain of colicin E3

The toxin colicin E3 targets the 30S subunit of bacterial ribosomes and cleaves a phosphodiester bond in the decoding center. We present the crystal structure of the 70S ribosome in complex with the cytotoxic domain of colicin E3 (E3-rRNase). The structure reveals how the rRNase domain of colicin binds to the A site of the decoding center in the 70S ribosome and cleaves the 16S ribosomal RNA (rRNA) between A1493 and G1494. The cleavage mechanism involves the concerted action of conserved residues Glu62 and His58 of the cytotoxic domain of colicin E3. These residues activate the 16S rRNA for 2' OH-induced hydrolysis. Conformational changes observed for E3-rRNase, 16S rRNA and helix 69 of 23S rRNA suggest that a dynamic binding platform is required for colicin E3 binding and function.

Use of the spliceosomal protein U1A to facilitate crystallization and structure determination of complex RNAs

The structure determination of complex RNA molecules such as ribozymes, riboswitches and aptamers by X-ray crystallography hinges on the preparation of well-ordered crystals. Success usually results from molecular engineering to facilitate crystallization. An approach that has resulted in 10 new RNA structures in the past decade is the use of the U1A crystallization module. In this approach, the cognate site for the U1A spliceosomal protein is introduced into a functionally dispensable location in the RNA of interest, and the RNA is cocrystallized with the basic RNA-binding protein. In addition to facilitating crystallization, the presence of U1A can be useful for de novo phase determination. In this paper, some general considerations for the use of this approach to RNA crystallization are presented, and specifics of the application of the U1A module to the crystallization of the hairpin ribozyme and the tetracycline aptamer are reviewed.

Promoting RNA helical stacking via A-minor junctions

RNA molecules take advantage of prevalent structural motifs to fold and assemble into well-defined 3D architectures. The A-minor junction is a class of RNA motifs that specifically controls coaxial stacking of helices in natural RNAs. A sensitive self-assembling supra-molecular system was used as an assay to compare several natural and previously unidentified A-minor junctions by native polyacrylamide gel electrophoresis and atomic force microscopy. This class of modular motifs follows a topological rule that can accommodate a variety of interchangeable A-minor interactions with distinct local structural motifs. Overall, two different types of A-minor junctions can be distinguished based on their functional self-assembling behavior: one group makes use of triloops or GNRA and GNRA-like loops assembling with helices, while the other takes advantage of more complex tertiary receptors specific for the loop to gain higher stability. This study demonstrates how different structural motifs of RNA can contribute to the formation of topologically equivalent helical stacks. It also exemplifies the need of classifying RNA motifs based on their tertiary structural features rather than secondary structural features. The A-minor junction rule can be used to facilitate tertiary structure prediction of RNAs and rational design of RNA parts for nanobiotechnology and synthetic biology.

Small self-cleaving ribozymes

The hammerhead, hairpin, hepatitis delta virus (HDV), Varkud Satellite (VS), and glmS ribozymes catalyze sequence-specific intramolecular cleavage of RNA. They range between 50 and 150 nucleotides in length, and are known as the "small self-cleaving ribozymes." Except for the glmS ribozyme that functions as a riboswitch in Gram-positive bacteria, they were originally discovered as domains of satellite RNAs. However, recent studies show that several of them are broadly distributed in genomes of organisms from many phyla. Each of these ribozymes has a unique overall architecture and active site organization. Crystal structures have revealed how RNA active sites can bind preferentially to the transition state of a reaction, whereas mechanistic studies have shown that nucleobases can efficiently perform general acid-base and electrostatic catalysis. This versatility explains the abundance of ribozymes in contemporary organisms and also supports a role for catalytic RNAs early in evolution.

Compact intermediates in RNA folding

Large noncoding RNAs fold into their biologically functional structures via compact yet disordered intermediates, which couple the stable secondary structure of the RNA with the emerging tertiary fold. The specificity of the collapse transition, which coincides with the assembly of helical domains, depends on RNA sequence and counterions. It determines the specificity of the folding pathways and the magnitude of the free energy barriers to the ensuing search for the native conformation. By coupling helix assembly with nascent tertiary interactions, compact folding intermediates in RNA also play a crucial role in ligand binding and RNA-protein recognition.

Nucleobase-mediated general acid-base catalysis in the Varkud satellite ribozyme

Existing evidence suggests that the Varkud satellite (VS) ribozyme accelerates the cleavage of a specific phosphodiester bond using general acid-base catalysis. The key functionalities are the nucleobases of adenine 756 in helix VI of the ribozyme, and guanine 638 in the substrate stem loop. This results in a bell-shaped dependence of reaction rate on pH, corresponding to groups with pK(a) = 5.2 and 8.4. However, it is not possible from those data to determine which nucleobase is the acid, and which the base. We have therefore made substrates in which the 5' oxygen of the scissile phosphate is replaced by sulfur. This labilizes the leaving group, removing the requirement for general acid catalysis. This substitution restores full activity to the highly impaired A756G ribozyme, consistent with general acid catalysis by A756 in the unmodified ribozyme. The pH dependence of the cleavage of the phosphorothiolate-modified substrates is consistent with general base catalysis by nucleobase at position 638. We conclude that cleavage of the substrate by the VS ribozyme is catalyzed by deprotonation of the 2'-O nucleophile by G638 and protonation of the 5'-O leaving group by A756.

Catalytic importance of a protonated adenosine in the hairpin ribozyme active site

The hairpin ribozyme accelerates the rate of phosphodiester transfer reactions by at least 5 orders of magnitude. To achieve this rate enhancement, the active site forms via a substrate helix docking event that constrains the scissile phosphate linkage and positions G8 and A38 for catalysis, both of which have been implicated as sites of proton transfer in general acid-base catalysis. To investigate the functional groups required for hairpin activity, we previously reported a series of nucleotide analogue interference mapping experiments [Ryder, S. P., et al. (2001) RNA 7, 1454-1463]. The critical functional groups implicated in those studies were largely consistent with subsequent X-ray crystal structures, but the lack of A38 interference with 8-azaadenosine (n(8)A), a pK(a) perturbed nucleotide analogue, argued against functional base ionization at this site. This is inconsistent with a transition state crystal structure and other biochemical studies. To address this discrepancy, we investigated the hairpin ribozyme with an expanded set of pK(a) perturbed adenosine analogues containing fluorine. A38 was the only site that showed persistent and strong interference with low pK(a) analogues across a variety of construct/substrate pairs. This interference pattern suggests that A38 base ionization is required for catalytic activity. The lack of n(8)A interference at A38, in spite of its reduced pK(a), likely results from n(8)A stabilization of the docked state, which requires an unusual syn glycosidic base conformation at A38 for active site assembly. The fluorinated adenosine analogues are better suited to identify sites of functional ionization in systems where structural rearrangements are closely coupled to catalytic steps. All pK(a) reduced analogues, including those of the previous study, produce selective interference at A38 when substrates are stably bound and docked, consistent with the importance of base ionization at this site.

Extensive molecular dynamics simulations showing that canonical G8 and protonated A38H+ forms are most consistent with crystal structures of hairpin ribozyme

The hairpin ribozyme is a prominent member of the group of small catalytic RNAs (RNA enzymes or ribozymes) because it does not require metal ions to achieve catalysis. Biochemical and structural data have implicated guanine 8 (G8) and adenine 38 (A38) as catalytic participants in cleavage and ligation catalyzed by the hairpin ribozyme, yet their exact role in catalysis remains disputed. To gain insight into dynamics in the active site of a minimal self-cleaving hairpin ribozyme, we have performed extensive classical, explicit-solvent molecular dynamics (MD) simulations on time scales of 50-150 ns. Starting from the available X-ray crystal structures, we investigated the structural impact of the protonation states of G8 and A38, and the inactivating A-1(2'-methoxy) substitution employed in crystallography. Our simulations reveal that a canonical G8 agrees well with the crystal structures while a deprotonated G8 profoundly distorts the active site. Thus MD simulations do not support a straightforward participation of the deprotonated G8 in catalysis. By comparison, the G8 enol tautomer is structurally well tolerated, causing only local rearrangements in the active site. Furthermore, a protonated A38H(+) is more consistent with the crystallography data than a canonical A38. The simulations thus support the notion that A38H(+) is the dominant form in the crystals, grown at pH 6. In most simulations, the canonical A38 departs from the scissile phosphate and substantially perturbs the structures of the active site and S-turn. Yet, we occasionally also observe formation of a stable A-1(2'-OH)...A38(N1) hydrogen bond, which documents the ability of the ribozyme to form this hydrogen bond, consistent with a potential role of A38 as general base catalyst. The presence of this hydrogen bond is, however, incompatible with the expected in-line attack angle necessary for self-cleavage, requiring a rapid transition of the deprotonated 2'-oxyanion to a position more favorable for in-line attack after proton transfer from A-1(2'-OH) to A38(N1). The simulations revealed a potential force field artifact, occasional but irreversible formation of "ladder-like", underwound A-RNA structure in one of the external helices. Although it does not affect the catalytic center of the hairpin ribozyme, further studies are under way to better assess possible influence of such force field behavior on long RNA simulations.

Enthalpy-driven RNA folding: single-molecule thermodynamics of tetraloop-receptor tertiary interaction

RNA folding thermodynamics are crucial for structure prediction, which requires characterization of both enthalpic and entropic contributions of tertiary motifs to conformational stability. We explore the temperature dependence of RNA folding due to the ubiquitous GAAA tetraloop-receptor docking interaction, exploiting immobilized and freely diffusing single-molecule fluorescence resonance energy transfer (smFRET) methods. The equilibrium constant for intramolecular docking is obtained as a function of temperature (T = 21-47 degrees C), from which a van't Hoff analysis yields the enthalpy (DeltaH degrees) and entropy (DeltaS degrees) of docking. Tetraloop-receptor docking is significantly exothermic and entropically unfavorable in 1 mM MgCl(2) and 100 mM NaCl, with excellent agreement between immobilized (DeltaH degrees = -17.4 +/- 1.6 kcal/mol, and DeltaS degrees = -56.2 +/- 5.4 cal mol(-1) K(-1)) and freely diffusing (DeltaH degrees = -17.2 +/- 1.6 kcal/mol, and DeltaS degrees = -55.9 +/- 5.2 cal mol(-1) K(-1)) species. Kinetic heterogeneity in the tetraloop-receptor construct is unaffected over the temperature range investigated, indicating a large energy barrier for interconversion between the actively docking and nondocking subpopulations. Formation of the tetraloop-receptor interaction can account for approximately 60% of the DeltaH degrees and DeltaS degrees of P4-P6 domain folding in the Tetrahymena ribozyme, suggesting that it may act as a thermodynamic clamp for the domain. Comparison of the isolated tetraloop-receptor and other tertiary folding thermodynamics supports a theme that enthalpy- versus entropy-driven folding is determined by the number of hydrogen bonding and base stacking interactions.

Effect of high hydrostatic pressure on hydration and activity of ribozymes

Formation and stabilization of RNA structure in the cell depends on its interaction with solvent and metal ions. High hydrostatic pressure (HHP) is a convenient tool in an analysis of the role of small molecules in the structure stabilization of biological macromolecules. Analysis of HHP effect and various concentrations of ions showed that water induce formation of the active ribozyme structure. So, it is clear that water is the driving force of conformational changes of nucleic acid.

Exploring RNA structure by integrative molecular modelling

RNA molecular modelling is adequate to rapidly tackle the structure of RNA molecules. With new structured RNAs constituting a central class of cellular regulators discovered every year, the need for swift and reliable modelling methods is more crucial than ever. The pragmatic method based on interactive all-atom molecular modelling relies on the observation that specific structural motifs are recurrently found in RNA sequences. Once identified by a combination of comparative sequence analysis and biochemical data, the motifs composing the secondary structure of a given RNA can be extruded in three dimensions (3D) and used as building blocks assembled manually during a bioinformatic interactive process. Comparing the models to the corresponding crystal structures has validated the method as being powerful to predict the RNA topology and architecture while being less accurate regarding the prediction of base-base interactions. These aspects as well as the necessary steps towards automation will be discussed.

NMR structure and dynamics of the Specifier Loop domain from the Bacillus subtilis tyrS T box leader RNA

Gram-positive bacteria utilize a tRNA-responsive transcription antitermination mechanism, designated the T box system, to regulate expression of many amino acid biosynthetic and aminoacyl-tRNA synthetase genes. The RNA transcripts of genes controlled by this mechanism contain 5′ untranslated regions, or leader RNAs, that specifically bind cognate tRNA molecules through pairing of nucleotides in the tRNA anticodon loop with nucleotides in the Specifier Loop domain of the leader RNA. We have determined the solution structure of the Specifier Loop domain of the tyrS leader RNA from Bacillus subtilis. Fifty percent of the nucleotides in the Specifier Loop domain adopt a loop E motif. The Specifier Sequence nucleotides, which pair with the tRNA anticodon, stack with their Watson–Crick edges rotated toward the minor groove and exhibit only modest flexibility. We also show that a Specifier Loop domain mutation that impairs the function of the B. subtilis glyQS T box RNA disrupts the tyrS loop E motif. Our results suggest a mechanism for tRNA–Specifier Loop binding in which the phosphate backbone kink created by the loop E motif causes the Specifier Sequence bases to rotate toward the minor groove, which increases accessibility for pairing with bases in the anticodon loop of tRNA.

Molecular dynamics and quantum mechanics of RNA: conformational and chemical change we can believe in

Structure and dynamics are both critical to RNA’s vital functions in biology. Numerous techniques can elucidate the structural dynamics of RNA, but computational approaches based on experimental data arguably hold the promise of providing the most detail. In this Account, we highlight areas wherein molecular dynamics (MD) and quantum mechanical (QM) techniques are applied to RNA, particularly in relation to complementary experimental studies.

We have expanded on atomic-resolution crystal structures of RNAs in functionally relevant states by applying explicit solvent MD simulations to explore their dynamics and conformational changes on the submicrosecond time scale. MD relies on simplified atomistic, pairwise additive interaction potentials (force fields). Because of limited sampling, due to the finite accessible simulation time scale and the approximated force field, high-quality starting structures are required.

Despite their imperfection, we find that currently available force fields empower MD to provide meaningful and predictive information on RNA dynamics around a crystallographically defined energy minimum. The performance of force fields can be estimated by precise QM calculations on small model systems. Such calculations agree reasonably well with the Cornell et al. AMBER force field, particularly for stacking and hydrogen-bonding interactions. A final verification of any force field is accomplished by simulations of complex nucleic acid structures.

The performance of the Cornell et al. AMBER force field generally corresponds well with and augments experimental data, but one notable exception could be the capping loops of double-helical stems. In addition, the performance of pairwise additive force fields is obviously unsatisfactory for inclusion of divalent cations, because their interactions lead to major polarization and charge-transfer effects neglected by the force field. Neglect of polarization also limits, albeit to a lesser extent, the description accuracy of other contributions, such as interactions with monovalent ions, conformational flexibility of the anionic sugar−phosphate backbone, hydrogen bonding, and solute polarization by solvent. Still, despite limitations, MD simulations are a valid tool for analyzing the structural dynamics of existing experimental structures. Careful analysis of MD simulations can identify problematic aspects of an experimental RNA structure, unveil structural characteristics masked by experimental constraints, reveal functionally significant stochastic fluctuations, evaluate the structural role of base ionization, and predict structurally and potentially functionally important details of the solvent behavior, including the presence of tightly bound water molecules. Moreover, combining classical MD simulations with QM calculations in hybrid QM/MM approaches helps in the assessment of the plausibility of chemical mechanisms of catalytic RNAs (ribozymes).

In contrast, the reliable prediction of structure from sequence information is beyond the applicability of MD tools. The ultimate utility of computational studies in understanding RNA function thus requires that the results are neither blindly accepted nor flatly rejected, but rather considered in the context of all available experimental data, with great care given to assessing limitations through the available starting structures, force field approximations, and sampling limitations. The examples given in this Account showcase how the judicious use of basic MD simulations has already served as a powerful tool to help evaluate the role of structural dynamics in biological function of RNA.

Direct Raman measurement of an elevated base pKa in the active site of a small ribozyme in a precatalytic conformation

Catalytic RNA molecules can achieve rate acceleration by shifting base pK(a) values toward neutrality. Prior evidence has suggested that base A38 of the hairpin ribozyme plays an important role in phosphoryl transfer, possibly functioning as a general acid, or by orienting a specific water molecule for proton transfer. To address the role of A38, we used Raman spectroscopy to measure directly the pK(a) of the N1-imino moiety in the context of hairpin ribozyme crystals representative of a "precatalytic" conformation. The results revealed that the pK(a) of A38 is shifted to 5.46 +/- 0.05 relative to 3.68 +/- 0.06 derived from a reference solution of the nucleotide AMP. The elevated pK(a) correlates well with the first titration point of the macroscopic pH-rate profile of the hairpin ribozyme in solution and strongly supports A38 as a general acid catalyst in bond scission. The results confirm that A38 is protonated before the transition state, which would promote phosphorane development. Overall, the results establish a cogent structure-function paradigm that expands our understanding of how RNA structure can enhance nucleobase reactivity to catalyze biological reactions.

Thermodynamic analysis of ligand binding and ligand binding-induced tertiary structure formation by the thiamine pyrophosphate riboswitch

The thi-box riboswitch regulates gene expression in response to the intracellular concentration of thiamine pyrophosphate (TPP) in archaea, bacteria, and eukarya. To complement previous biochemical, genetic, and structural studies of this phylogenetically widespread RNA domain, we have characterized its interaction with TPP by isothermal titration calorimetry. This shows that TPP binding is highly dependent on Mg(2+) concentration. The dissociation constant decreases from approximately 200 nM at 0.5 mM Mg(2+) concentration to approximately 9 nM at 2.5 mM Mg(2+) concentration. Binding is enthalpically driven, but the unfavorable entropy of binding decreases as Mg(2+) concentration rises, suggesting that divalent cations serve to pre-organize the RNA. Mutagenesis, biochemical analysis, and a new crystal structure of the riboswitch suggest that a critical element that participates in organizing the riboswitch structure is the tertiary interaction formed between the P3 and L5 regions. This tertiary contact is distant from the TPP binding site, but calorimetric analysis reveals that even subtle mutations in L5 can have readily detectable effects on TPP binding. The thermodynamic signatures of these mutations, namely decreased favorable enthalpy of binding and small effects on entropy of binding, are consistent with the P3-L5 association contributing allosterically to TPP-induced compaction of the RNA.

Recognition of the bacterial second messenger cyclic diguanylate by its cognate riboswitch

The cyclic diguanylate (bis-(3'-5')-cyclic dimeric guanosine monophosphate, c-di-GMP) riboswitch is the first known example of a gene-regulatory RNA that binds a second messenger. c-di-GMP is widely used by bacteria to regulate processes ranging from biofilm formation to the expression of virulence genes. The cocrystal structure of the c-di-GMP responsive GEMM riboswitch upstream of the tfoX gene of Vibrio cholerae reveals the second messenger binding the RNA at a three-helix junction. The two-fold symmetric second messenger is recognized asymmetrically by the monomeric riboswitch using canonical and noncanonical base-pairing as well as intercalation. These interactions explain how the RNA discriminates against cyclic diadenylate (c-di-AMP), a putative bacterial second messenger. Small-angle X-ray scattering and biochemical analyses indicate that the RNA undergoes compaction and large-scale structural rearrangement in response to ligand binding, consistent with organization of the core three-helix junction of the riboswitch concomitant with binding of c-di-GMP.

Modulating RNA structure and catalysis: lessons from small cleaving ribozymes

RNA is a key molecule in life, and comprehending its structure/function relationships is a crucial step towards a more complete understanding of molecular biology. Even though most of the information required for their correct folding is contained in their primary sequences, we are as yet unable to accurately predict both the folding pathways and active tertiary structures of RNA species. Ribozymes are interesting molecules to study when addressing these questions because any modifications in their structures are often reflected in their catalytic properties. The recent progress in the study of the structures, the folding pathways and the modulation of the small ribozymes derived from natural, self-cleaving, RNA motifs have significantly contributed to today's knowledge in the field.

Comparative gel electrophoresis analysis of helical junctions in RNA

Comparative gel electrophoresis provides information on the relative angles subtended between helical arms at a branchpoint in RNA. It is based upon the comparison of electrophoretic mobility in polyacrylamide gels of species containing two long arms, with the remaining one(s) being significantly shorter. Although the method currently lacks a really well-established basis of physical theory, it is very powerful, yet simple to apply. It has had a number of significant successes in RNA, DNA and DNA-protein complexes, and in all cases to date the results have stood the test of time and eventual comparison with crystallographic analysis.

Raman crystallography of RNA

Raman crystallography is the application of Raman spectroscopy to single crystals. This technique has been applied to a variety of protein molecules where it has provided unique information about biopolymer folding, substrate binding, and catalysis. Here, we describe the application of Raman crystallography to functional RNA molecules. RNA represents unique opportunities and challenges for Raman crystallography. One issue that confounds studies of RNA is its tendency to adopt multiple non-functional folds. Raman crystallography has the advantage that it isolates a single state of the RNA within the crystal and can evaluate its fold, metal ion binding properties (ligand identity, stoichiometry, and affinity), proton binding properties (identity, stoichiometry, and affinity), and catalytic potential. In particular, base-specific stretches can be identified and then associated with the binding of metal ions and protons. Because measurements are carried out in the hanging drop at ambient, rather than cryo, conditions and because RNA crystals tend to be approximately 70% solvent, RNA dynamics and conformational changes become experimentally accessible. This review focuses on experimental setup and procedures, acquisition and interpretation of Raman data, and determination of physicochemical properties of the RNA. Raman crystallographic and solution biochemical experiments on the HDV RNA enzyme are summarized and found to be in excellent agreement. Remarkably, characterization of the crystalline state has proven to help rather than hinder functional characterization of functional RNA, most likely because the tendency of RNA to fold heterogeneously is limited in a crystalline environment. Future applications of Raman crystallography to RNA are briefly discussed.

Theoretical studies of RNA catalysis: hybrid QM/MM methods and their comparison with MD and QM

Hybrid QM/MM methods combine the rigor of quantum mechanical (QM) calculations with the low computational cost of empirical molecular mechanical (MM) treatment allowing to capture dynamic properties to probe critical atomistic details of enzyme reactions. Catalysis by RNA enzymes (ribozymes) has only recently begun to be addressed with QM/MM approaches and is thus still a field under development. This review surveys methodology as well as recent advances in QM/MM applications to RNA mechanisms, including those of the HDV, hairpin, and hammerhead ribozymes, as well as the ribosome. We compare and correlate QM/MM results with those from QM and/or molecular dynamics (MD) simulations, and discuss scope and limitations with a critical eye on current shortcomings in available methodologies and computer resources. We thus hope to foster mutual appreciation and facilitate collaboration between experimentalists and theorists to jointly advance our understanding of RNA catalysis at an atomistic level.

Exploring ribozyme conformational changes with X-ray crystallography

Relating three-dimensional fold to function is a central challenge in RNA structural biology. Toward this goal, X-ray crystallography has long been considered the "gold standard" for structure determinations at atomic resolution, although NMR spectroscopy has become a powerhouse in this arena as well. In the area of dynamics, NMR remains the dominant technique to probe the magnitude and timescales of molecular motion. Although the latter area remains largely unassailable by conventional crystallographic methods, inroads have been made on proteins using Laue radiation on timescales of ms to ns. Proposed 'fourth generation' radiation sources, such as free-electron X-ray lasers, promise ps- to fs-timescale resolution, and credible evidence is emerging that supports the feasibility of single molecule imaging. At present however, the preponderance of RNA structural information has been derived from timescale and motion insensitive crystallographic techniques. Importantly, developments in computing, automation and high-flux synchrotron sources have propelled the rapidity of 'conventional' RNA crystal structure determinations to timeframes of hours once a suitable set of phases is obtained. With a sufficient number of crystal structures, it is possible to create a structural ensemble that can provide insight into global and local molecular motion characteristics that are relevant to biological function. Here we describe techniques to explore conformational changes in the hairpin ribozyme, a representative non-protein-coding RNA catalyst. The approaches discussed include: (i) construct choice and design using prior knowledge to improve X-ray diffraction; (ii) recognition of long-range conformational changes and (iii) use of single-base or single-atom changes to create ensembles. The methods are broadly applicable to other RNA systems.

Single-atom imino substitutions at A9 and A10 reveal distinct effects on the fold and function of the hairpin ribozyme catalytic core

The hairpin ribozyme cleaves a phosphodiester bond within a cognate substrate. Structural and biochemical data indicate the conserved A9 and A10 bases reside close to the scissile bond but make distinct contributions to catalysis. To investigate these residues, we replaced the imino moiety of each base with N1-deazaadenosine. This single-atom change resulted in an 8-fold loss in k(obs) for A9 and displacement of the base from the active site; no effects were observed for A10. We propose that the imino moiety of A9 promotes a key water-mediated contact that favors transition-state formation, which suggests an enhanced chemical repertoire for RNA.

Structure and function of regulatory RNA elements: ribozymes that regulate gene expression

Since their discovery in the 1980s, it has gradually become apparent that there are several functional classes of naturally occurring ribozymes. These include ribozymes that mediate RNA splicing (the Group I and Group II introns, and possibly the RNA components of the spliceosome), RNA processing ribozymes (RNase P, which cleaves precursor tRNAs and other structural RNA precursors), the peptidyl transferase center of the ribosome, and small, self-cleaving genomic ribozymes (including the hammerhead, hairpin, HDV and VS ribozymes). The most recently discovered functional class of ribozymes include those that are embedded in the untranslated regions of mature mRNAs that regulate the gene's translational expression. These include the prokaryotic glmS ribozyme, a bacterial riboswitch, and a variant of the hammerhead ribozyme, which has been found embedded in mammalian mRNAs. With the discovery of a mammalian riboswitch ribozyme, the question of how an embedded hammerhead ribozyme's switching mechanism works becomes a compelling question. Recent structural results suggest several possibilities.

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.