Structural and catalytic effects of an invariant purine substitution in the hammerhead ribozyme: implications for the mechanism of acid-base catalysis.

The hammerhead ribozyme catalyzes RNA cleavage via acid-base catalysis. Whether it does so by general acid-base catalysis, in which the RNA itself donates and abstracts protons in the transition state, as is typically assumed, or by specific acid-base catalysis, in which the RNA plays a structural role and proton transfer is mediated by active-site water molecules, is unknown. Previous biochemical and crystallographic experiments implicate an invariant purine in the active site, G12, as the general base. However, G12 may play a structural role consistent with specific base catalysis. To better understand the role of G12 in the mechanism of hammerhead catalysis, a 2.2 Šresolution crystal structure of a hammerhead ribozyme from Schistosoma mansoni with a purine substituted for G12 in the active site of the ribozyme was obtained. Comparison of this structure (PDB entry 3zd4), in which A12 is substituted for G, with three previously determined structures that now serve as important experimental controls, allows the identification of structural perturbations that are owing to the purine substitution itself. Kinetic measurements for G12 purine-substituted schistosomal hammerheads confirm a previously observed dependence of rate on the pK(a) of the substituted purine; in both cases inosine, which is similar to G in pK(a) and hydrogen-bonding properties, is unexpectedly inactive. Structural comparisons indicate that this may primarily be owing to the lack of the exocyclic 2-amino group in the G12A and G12I substitutions and its structural effect upon both the nucleotide base and phosphate of A9. The latter involves the perturbation of a previously identified and well characterized metal ion-binding site known to be catalytically important in both minimal and full-length hammerhead ribozyme sequences. The results permit it to be suggested that G12 plays an important role in stabilizing the active-site structure. This result, although not inconsistent with the potential role of G12 as a general base, indicates that an alternative hammerhead cleavage mechanism involving specific base catalysis may instead explain the observed rate dependence upon purine substitutions at G12. The crystallographic results, contrary to previous assumptions, therefore cannot be interpreted to favor the general base catalysis mecahnism over the specific base catalysis mechanism. Instead, both of these mutually exclusive mechanistic alternatives must be considered in light of the current structural and biochemical data.

A discontinuous hammerhead ribozyme embedded in a mammalian messenger RNA.

Structured RNAs embedded in the untranslated regions (UTRs) of messenger RNAs can regulate gene expression. In bacteria, control of a metabolite gene is mediated by the self-cleaving activity of a ribozyme embedded in its 5' UTR. This discovery has raised the question of whether gene-regulating ribozymes also exist in eukaryotic mRNAs. Here we show that highly active hammerhead ribozymes are present in the 3' UTRs of rodent C-type lectin type II (Clec2) genes. Using a hammerhead RNA motif search with relaxed delimitation of the non-conserved regions, we detected ribozyme sequences in which the invariant regions, in contrast to the previously identified continuous hammerheads, occur as two fragments separated by hundreds of nucleotides. Notably, a fragment pair can assemble to form an active hammerhead ribozyme structure between the translation termination and the polyadenylation signals within the 3' UTR. We demonstrate that this hammerhead structure can self-cleave both in vitro and in vivo, and is able to reduce protein expression in mouse cells. These results indicate that an unrecognized mechanism of post-transcriptional gene regulation involving association of discontinuous ribozyme sequences within an mRNA may be modulating the expression of several CLEC2 proteins that function in bone remodelling and the immune response of several mammals.

Challenges and surprises that arise with nucleic acids during model building and refinement.

The process of building and refining crystal structures of nucleic acids, although similar to that for proteins, has some peculiarities that give rise to both various complications and various benefits. Although conventional isomorphous replacement phasing techniques are typically used to generate an experimental electron-density map for the purposes of determining novel nucleic acid structures, it is also possible to couple the phasing and model-building steps to permit the solution of complex and novel RNA three-dimensional structures without the need for conventional heavy-atom phasing approaches.

Identification of dynamical hinge points of the L1 ligase molecular switch.

The L1 ligase is an in vitro selected ribozyme that uses a noncanonically base-paired ligation site to catalyze regioselectively and regiospecifically the 5' to 3' phosphodiester bond ligation, a reaction relevant to origin of life hypotheses that invoke an RNA world scenario. The L1 ligase crystal structure revealed two different conformational states that were proposed to represent the active and inactive forms. It remains an open question as to what degree these two conformers persist as stable conformational intermediates in solution, and along what pathway are they able to interconvert. To explore these questions, we have performed a series of molecular dynamics simulations in explicit solvent of the inactive-active conformational switch in L1 ligase. Four simulations were performed departing from both conformers in both the reactant and product states, in addition to a simulation where local unfolding in the active state was induced. From these simulations, along with crystallographic data, a set of four virtual torsion angles that span two evolutionarily conserved and restricted regions were identified as dynamical hinge points in the conformational switch transition. The ligation site visits three distinct states characterized by hydrogen bond patterns that are correlated with the formation of specific contacts that may promote catalysis. The insights gained from these simulations contribute to a more detailed understanding of the coupled catalytic/conformational switch mechanism of L1 ligase that may facilitate the design and engineering of new catalytic riboswitches.

Capturing hammerhead ribozyme structures in action by modulating general base catalysis.

We have obtained precatalytic (enzyme-substrate complex) and postcatalytic (enzyme-product complex) crystal structures of an active full-length hammerhead RNA that cleaves in the crystal. Using the natural satellite tobacco ringspot virus hammerhead RNA sequence, the self-cleavage reaction was modulated by substituting the general base of the ribozyme, G12, with A12, a purine variant with a much lower pKa that does not significantly perturb the ribozyme's atomic structure. The active, but slowly cleaving, ribozyme thus permitted isolation of enzyme-substrate and enzyme-product complexes without modifying the nucleophile or leaving group of the cleavage reaction, nor any other aspect of the substrate. The predissociation enzyme-product complex structure reveals RNA and metal ion interactions potentially relevant to transition-state stabilization that are absent in precatalytic structures.

Solving novel RNA structures using only secondary structural fragments.

The crystallographic phase problem is the primary bottleneck encountered when attempting to solve macromolecular structures for which no close crystallographic structural homologues are known. Typically, isomorphous "heavy-atom" replacement and/or anomalous dispersion methods must be used in such cases to obtain experimentally-determined phases. Even three-dimensional NMR structures of the same macromolecule are often not sufficient to solve the crystallographic phase problem. RNA crystal structures present additional challenges due to greater difficulty in obtaining suitable heavy-atom derivatives. We present a unique approach to solve the phase problem for novel RNA crystal structures that has enjoyed a reasonable degree of success. This approach involves modeling only those portions of the RNA sequence whose structure can be predicted readily, i.e., the individual A-form helical regions and well-known stem-loop sub-structures. We have found that no prior knowledge of how the helices and other structural elements are arranged with respect to one another in three-dimensional space, or in some cases, even the sequence, is required to obtain a useable solution to the phase problem, using simultaneous molecular replacement of a set of generic helical RNA fragments.

Structure of Escherichia coli RNase E catalytic domain and implications for RNA turnover.

The coordinated regulation of gene expression is required for homeostasis, growth and development in all organisms. Such coordination may be partly achieved at the level of messenger RNA stability, in which the targeted destruction of subsets of transcripts generates the potential for cross-regulating metabolic pathways. In Escherichia coli, the balance and composition of the transcript population is affected by RNase E, an essential endoribonuclease that not only turns over RNA but also processes certain key RNA precursors. RNase E cleaves RNA internally, but its catalytic power is determined by the 5' terminus of the substrate, even if this lies at a distance from the cutting site. Here we report crystal structures of the catalytic domain of RNase E as trapped allosteric intermediates with RNA substrates. Four subunits of RNase E catalytic domain associate into an interwoven quaternary structure, explaining why the subunit organization is required for catalytic activity. The subdomain encompassing the active site is structurally congruent to a deoxyribonuclease, making an unexpected link in the evolutionary history of RNA and DNA nucleases. The structure explains how the recognition of the 5' terminus of the substrate may trigger catalysis and also sheds light on the question of how RNase E might selectively process, rather than destroy, specific RNA precursors.

Crystal structure of a translation termination complex formed with release factor RF2.

We report the crystal structure of a translation termination complex formed by the Thermus thermophilus 70S ribosome bound with release factor RF2, in response to a UAA stop codon, solved at 3 A resolution. The backbone of helix alpha5 and the side chain of serine of the conserved SPF motif of RF2 recognize U1 and A2 of the stop codon, respectively. A3 is unstacked from the first 2 bases, contacting Thr-216 and Val-203 of RF2 and stacking on G530 of 16S rRNA. The structure of the RF2 complex supports our previous proposal that conformational changes in the ribosome in response to recognition of the stop codon stabilize rearrangement of the switch loop of the release factor, resulting in docking of the universally conserved GGQ motif in the PTC of the 50S subunit. As seen for the RF1 complex, the main-chain amide nitrogen of glutamine in the GGQ motif is positioned to contribute directly to catalysis of peptidyl-tRNA hydrolysis, consistent with mutational studies, which show that most side-chain substitutions of the conserved glutamine have little effect. We show that when the H-bonding capability of the main-chain N-H of the conserved glutamine is eliminated by substitution with proline, peptidyl-tRNA esterase activity is abolished, consistent with its proposed role in catalysis.

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.

The RNA WikiProject: community annotation of RNA families.

The online encyclopedia Wikipedia has become one of the most important online references in the world and has a substantial and growing scientific content. A search of Google with many RNA-related keywords identifies a Wikipedia article as the top hit. We believe that the RNA community has an important and timely opportunity to maximize the content and quality of RNA information in Wikipedia. To this end, we have formed the RNA WikiProject (http://en.wikipedia.org/wiki/Wikipedia:WikiProject_RNA) as part of the larger Molecular and Cellular Biology WikiProject. We have created over 600 new Wikipedia articles describing families of noncoding RNAs based on the Rfam database, and invite the community to update, edit, and correct these articles. The Rfam database now redistributes this Wikipedia content as the primary textual annotation of its RNA families. Users can, therefore, for the first time, directly edit the content of one of the major RNA databases. We believe that this Wikipedia/Rfam link acts as a functioning model for incorporating community annotation into molecular biology databases.

Novel ring structure in the gp41 trimer of human immunodeficiency virus type 1 that modulates sensitivity and resistance to broadly neutralizing antibodies.

The identification of the determinants of sensitivity and resistance to broadly neutralizing antibodies is a high priority for human immunodeficiency virus (HIV) research. An analysis of the swarm of closely related envelope protein variants in an HIV-infected individual revealed a mutation that markedly affected sensitivity to neutralization by antibodies and antiviral entry inhibitors targeting both gp41 and gp120. This mutation mapped to the C34 helix of gp41 and disrupted an unexplored structural feature consisting of a ring of hydrogen bonds in the gp41 trimer. This mutation appeared to affect the assembly of the six-helix bundle required for virus fusion and to alter the conformational equilibria so as to favor the prehairpin intermediate conformation required for the binding of the membrane proximal external region-specific neutralizing antibodies 2F5 and 4E10 and the antiviral drug enfuvirtide (Fuzeon). The "swarm analysis" method we describe furthers our understanding of the relationships among the structure, function, and antigenicity of the HIV envelope protein and represents a new approach to the identification of vaccine antigens.

Ribozymes.

The structural molecular biology of ribozymes took another great leap forward during the past two years. Before ribozymes were discovered in the early 1980s, all enzymes were thought to be proteins. No detailed structural information on ribozymes became available until 1994. Now, within the past two years, near atomic resolution crystal structures are available for almost all of the known ribozymes. The latest additions include ribonuclease P, group I intron structures, the ribosome (the peptidyl transferase appears to be a ribozyme) and several smaller ribozymes, including a Diels-Alderase, the glmS ribozyme and a new hammerhead ribozyme structure that reconciles 12 years of discord. Although not all ribozymes are metalloenzymes, acid-base catalysis appears to be a universal property shared by all ribozymes as well as many of their protein cousins.

The crystal structure of the Escherichia coli RNase E apoprotein and a mechanism for RNA degradation.

RNase E is an essential bacterial endoribonuclease involved in the turnover of messenger RNA and the maturation of structured RNA precursors in Escherichia coli. Here, we present the crystal structure of the E. coli RNase E catalytic domain in the apo-state at 3.3 A. This structure indicates that, upon catalytic activation, RNase E undergoes a marked conformational change characterized by the coupled movement of two RNA-binding domains to organize the active site. The structural data suggest a mechanism of RNA recognition and cleavage that explains the enzyme's preference for substrates possessing a 5'-monophosphate and accounts for the protective effect of a triphosphate cap for most transcripts. Internal flexibility within the quaternary structure is also observed, a finding that has implications for recognition of structured RNA substrates and for the mechanism of internal entry for a subset of substrates that are cleaved without 5'-end requirements.

Solvent structure and hammerhead ribozyme catalysis.

Although the hammerhead ribozyme is regarded as a prototype for understanding RNA catalysis, the mechanistic roles of associated metal ions and water molecules in the cleavage reaction remain controversial. We have investigated the catalytic potential of observed divalent metal ions and water molecules bound to a 2 A structure of the full-length hammerhead ribozyme by using X-ray crystallography in combination with molecular dynamics simulations. A single Mn(2+) is observed to bind directly to the A9 phosphate in the active site, accompanying a hydrogen-bond network involving a well-ordered water molecule spanning N1 of G12 (the general base) and 2'-O of G8 (previously implicated in general acid catalysis) that we propose, based on molecular dynamics calculations, facilitates proton transfer in the cleavage reaction. Phosphate-bridging metal interactions and other mechanistic hypotheses are also tested with this approach.

Role of Mg2+ in hammerhead ribozyme catalysis from molecular simulation.

Molecular dynamics simulations have been performed to investigate the role of Mg2+ in the full-length hammerhead ribozyme cleavage reaction. In particular, the aim of this work is to characterize the binding mode and conformational events that give rise to catalytically active conformations and stabilization of the transition state. Toward this end, a series of eight 12 ns molecular dynamics simulations have been performed with different divalent metal binding occupations for the reactant, early and late transition state using recently developed force field parameters for metal ions and reactive intermediates in RNA catalysis. In addition, hybrid QM/MM calculations of the early and late transition state were performed to study the proton-transfer step in general acid catalysis that is facilitated by the catalytic Mg2+ ion. The simulations suggest that Mg2+ is profoundly involved in the hammerhead ribozyme mechanism both at structural and catalytic levels. Binding of Mg2+ in the active site plays a key structural role in the stabilization of stem I and II and to facilitate formation of near attack conformations and interactions between the nucleophile and G12, the implicated general base catalyst. In the transition state, Mg2+ binds in a bridging position where it stabilizes the accumulated charge of the leaving group while interacting with the 2'OH of G8, the implicated general acid catalyst. The QM/MM simulations provide support that, in the late transition state, the 2'OH of G8 can transfer a proton to the leaving group while directly coordinating the bridging Mg2+ ion. The present study provides evidence for the role of Mg2+ in hammerhead ribozyme catalysis. The proposed simulation model reconciles the interpretation of available experimental structural and biochemical data, and provides a starting point for more detailed investigation of the chemical reaction path with combined QM/MM methods.

The structure of a rigorously conserved RNA element within the SARS virus genome.

We have solved the three-dimensional crystal structure of the stem-loop II motif (s2m) RNA element of the SARS virus genome to 2.7-A resolution. SARS and related coronaviruses and astroviruses all possess a motif at the 3' end of their RNA genomes, called the s2m, whose pathogenic importance is inferred from its rigorous sequence conservation in an otherwise rapidly mutable RNA genome. We find that this extreme conservation is clearly explained by the requirement to form a highly structured RNA whose unique tertiary structure includes a sharp 90 degrees kink of the helix axis and several novel longer-range tertiary interactions. The tertiary base interactions create a tunnel that runs perpendicular to the main helical axis whose interior is negatively charged and binds two magnesium ions. These unusual features likely form interaction surfaces with conserved host cell components or other reactive sites required for virus function. Based on its conservation in viral pathogen genomes and its absence in the human genome, we suggest that these unusual structural features in the s2m RNA element are attractive targets for the design of anti-viral therapeutic agents. Structural genomics has sought to deduce protein function based on three-dimensional homology. Here we have extended this approach to RNA by proposing potential functions for a rigorously conserved set of RNA tertiary structural interactions that occur within the SARS RNA genome itself. Based on tertiary structural comparisons, we propose the s2m RNA binds one or more proteins possessing an oligomer-binding-like fold, and we suggest a possible mechanism for SARS viral RNA hijacking of host protein synthesis, both based upon observed s2m RNA macromolecular mimicry of a relevant ribosomal RNA fold.

Structural Simplicity and Mechanistic Complexity in the Hammerhead Ribozyme.

Natural or full-length hammerhead ribozymes are up to 1000-fold more active than their minimal counterparts that lack a complex tertiary interaction that pre-organizes and stabilizes the ribozyme active site, positioning RNA functional groups to facilitate acid-base catalysis. The recent discovery that a single tertiary contact (an AU Hoogsteen pair) between Stems I and II confers essentially all of the enhanced activity greatly simplifies our understanding of the structural requirements for hammerhead ribozyme activity. In contrast, the simplest mechanistic interpretations are challenged with the presentation of more complex alternatives. These alternatives are elucidated and critically analyzed in the context of several of the active hammerhead ribozyme structures now available.

A pH-dependent conformational change, rather than the chemical step, appears to be rate-limiting in the hammerhead ribozyme cleavage reaction.

We have investigated the chemical basis for a previously observed 7.8 A conformational change in the hammerhead ribozyme that positions the substrate for in-line attack. We have found that the conformational change can only be observed at or above pH 8.5 (in the presence of Co(2+)) and requires the presence of an ionizable 2'-OH at the cleavage site, and note that this observed apparent pK(a) of 8.5 for the conformational change is within experimental error (+/-0.5) of the previously reported apparent kinetic pK(a) of 8.5 for the hammerhead ribozyme in the presence of Co(2+). We have solved two crystal structures of hammerhead ribozymes having 2'-OCH(3) or 2'-F substitutions at the cleavage site and have found that these will not undergo a conformational change equivalent to that observed for the hammerhead ribozyme having an unmodified attacking nucleophile under otherwise identical conditions. We have also characterized the kinetics of cleavage in the crystal. In addition to verifying that the particular sequence of RNA that we crystallized cleaves faster in the crystal than in solution, we also find that the extent of cleavage in the crystal is complete, unlike in solution where this and most other hammerhead ribozyme substrates are cleaved only to about 70 % completion. The initial cleavage rate in the crystal obeys the expected log-linear relation between cleavage-rate and pH with a slope of 0.7, as has been observed for other hammerhead ribozyme sequences in solution, indicating that in both the crystal and in solution the pH-dependent step is rate-limiting. However, the cleavage rate in the crystal is biphasic, with the most dramatic distinction between initial (slower) and final (faster) phases appearing at pH 6.0. The initial phase corresponds to the pH-dependent cleavage rate observed in solution, but the second, faster phase is roughly pH-independent and closely parallels the cleavage rate observed at pH 8 (0.4/minute). This result is particularly remarkable because it entails that the rapidly cleaving phase at pH 6 is comparable to the cleavage rate for the fastest cleaving hammerhead ribozymes at pH 6. Based upon these observations, we conclude that the pH-dependent conformational change is the rate-determining step under standard conditions for the hammerhead ribozyme self-cleavage reaction, and that an ionizable 2'-proton at cleavage site is required for this conformational change. We further hypothesize that deprotonation of the cleavage-site 2'-oxygen drives this conformational change.