The glmS ribozyme tunes the catalytically critical pK(a) of its coenzyme glucosamine-6-phosphate

The glmS ribozyme riboswitch is the first known natural catalytic RNA that employs a small-molecule cofactor. Binding of glucosamine-6-phosphate (GlcN6P) uncovers the latent self-cleavage activity of the RNA, which adopts a catalytically competent conformation that is nonetheless inactive in the absence of GlcN6P. Structural and analogue studies suggest that the amine of GlcN6P functions as a general acid-base catalyst, while its phosphate is important for binding affinity. However, the solution pK(a) of the amine is 8.06 ± 0.05, which is not optimal for proton transfer. Here we used Raman crystallography directly to determine the pK(a)'s of GlcN6P bound to the glmS ribozyme. Binding to the RNA lowers the pK(a) of the amine of GlcN6P to 7.26 ± 0.09 and raises the pK(a) of its phosphate to 6.35 ± 0.09. Remarkably, the pK(a)'s of these two functional groups are unchanged from their values for free GlcN6P (8.06 ± 0.05 and 5.98 ± 0.05, respectively) when GlcN6P binds to the catalytically inactive but structurally unperturbed G40A mutant of the ribozyme, thus implicating the ribozyme active site guanine in pK(a) tuning. This is the first demonstration that a ribozyme can tune the pK(a) of a small-molecule ligand. Moreover, the anionic glmS ribozyme in effect stabilizes the neutral amine of GlcN6P by lowering its pK(a). This is unprecedented and illustrates the chemical sophistication of ribozyme active sites.     

Essential role of an active-site guanine in glmS ribozyme catalysis

The glmS ribozyme is a catalytic riboswitch that is activated for endonucleolytic cleavage by the coenzyme glucosamine-6-phosphate. Using kinetic assays and X-ray crystallography, we identify an active-site mutation of a conserved guanine that abolishes catalysis without perturbing coenzyme binding. Our results provide evidence that coenzyme function requires a specific nucleobase to interact with the nucleophile of the cleavage reaction.     

Structural investigation of the GlmS ribozyme bound to Its catalytic cofactor

The GlmS riboswitch is located in the 5'-untranslated region of the gene encoding glucosamine-6-phosphate (GlcN6P) synthetase. The GlmS riboswitch is a ribozyme with activity triggered by binding of the metabolite GlcN6P. Presented here is the structure of the GlmS ribozyme (2.5 A resolution) with GlcN6P bound in the active site. The GlmS ribozyme adopts a compact double pseudoknot tertiary structure, with two closely packed helical stacks. Recognition of GlcN6P is achieved through coordination of the phosphate moiety by two hydrated magnesium ions as well as specific nucleobase contacts to the GlcN6P sugar ring. Comparison of this activator bound and the previously published apoenzyme complex supports a model in which GlcN6P does not induce a conformational change in the RNA, as is typical of other riboswitches, but instead functions as a catalytic cofactor for the reaction. This demonstrates that RNA, like protein enzymes, can employ the chemical diversity of small molecules to promote catalytic activity.     

Synthesis and evaluation of glucosamine-6-phosphate analogues as activators of glmS riboswitch

The glmS riboswitch is a ribozyme found in numerous Gram-positive bacteria and responds to the cellular concentrations of glucosamine 6-phosphate (GlcN6P). Given the importance of GlcN6P for cell wall biosynthesis, the glmS riboswitch has become a new drug target for the development of antibiotics. Herein, we describe the efficient synthesis of three GlcN6P analogues and their evaluation on inducing self-cleavage of the glmS riboswitch from Bacillus subtilis. Our results provide valuable information for further elucidation of the structure-activity relationships and drug design for glmS riboswitch antibiotics.     

Theoretical examination of Mg(2+)-mediated hydrolysis of a phosphodiester linkage as proposed for the hammerhead ribozyme

The hammerhead ribozyme is an RNA molecule capable of self-cleavage at a unique site within its sequence. Hydrolysis of this phosphodiester linkage has been proposed to occur via an in-line attack geometry for nucleophilic displacement by the 2'-hydroxyl on the adjoining phosphorus to generate a 2',3'-cyclic phosphate ester with elimination of the 5'-hydroxyl group, requiring a divalent metal ion under physiological conditions. The proposed S(N)2(P) reaction mechanism was investigated using density functional theory calculations incorporating the hybrid functional B3LYP to study this metal ion-dependent reaction with a tetraaquo magnesium (II)-bound hydroxide ion. For the Mg(2+)-catalyzed reaction, the gas-phase geometry optimized calculations predict two transition states with a kinetically insignificant, yet clearly defined, pentacoordinate intermediate. The first transition state located for the reaction is characterized by internal nucleophilic attack coupled to proton transfer. The second transition state, the rate-determining step, involves breaking of the exocyclic P-O bond where a metal-ligated water molecule assists in the departure of the leaving group. These calculations demonstrate that the reaction mechanism incorporating a single metal ion, serving as a Lewis acid, functions as a general base and can afford the necessary stabilization to the leaving group by orienting a water molecule for catalysis.     

An active-site guanine participates in glmS ribozyme catalysis in its protonated state

Active-site guanines that occupy similar positions have been proposed to serve as general base catalysts in hammerhead, hairpin, and glmS ribozymes, but no specific roles for these guanines have been demonstrated conclusively. Structural studies place G33(N1) of the glmS ribozyme of Bacillus anthracis within hydrogen-bonding distance of the 2'-OH nucleophile. Apparent pK(a) values determined from the pH dependence of cleavage kinetics for wild-type and mutant glmS ribozymes do not support a role for G33, or any other active-site guanine, in general base catalysis. Furthermore, discrepancies between apparent pK(a) values obtained from functional assays and microscopic pK(a) values obtained from pH-fluorescence profiles with ribozymes containing a fluorescent guanosine analogue, 8-azaguanosine, at position 33 suggest that the pH-dependent step in catalysis does not involve G33 deprotonation. These results point to an alternative model in which G33(N1) in its neutral, protonated form donates a hydrogen bond to stabilize the transition state.     

An unusual pH-independent and metal-ion-independent mechanism for hairpin ribozyme catalysis

Background: Hairpin ribozymes (RNA enzymes) catalyze the same chemical reaction as ribonuclease A and yet RNAs do not usually have functional groups analogous to the catalytically essential histidine and lysine sidechains of protein ribonucleases. Some RNA enzymes appear to recruit metal ions to act as Lewis acids in charge stabilization and metal-bound hydroxide for general base catalysis, but it has been reported that the hairpin ribozyme functions in the presence of metal ion chelators. This led us to investigate whether the hairpin ribozyme exploits a metal-ion-independent catalytic strategy.Results: Substitution of sulfur for nonbridging oxygens of the reactive phosphate of the hairpin ribozyme has small, stereospecific and metal-ionindependent effects on cleavage and ligation mediated by this ribozyme. Cobalt hexammine, an exchange-inert metal complex, supports full hairpin ribozyme activity, and the ribozyme's catalytic rate constants display only a shallow dependence on pH.Conclusions: Direct metal ion coordination to phosphate oxygens is not essential for hairpin ribozyme catalysis and metal-bound hydroxide does not serve as the general base in this catalysis. Several models might account for the unusual pH and metal ion independence: hairpin cleavage and ligation might be limited by a slow conformational change; a pH-independent or metalcation-independent chemical step, such as breaking the 5′ oxygen-phosphorus bond, might be rate determining; or finally, functional groups within the ribozyme might participate directly in catalytic chemistry. Whichever the case, the hairpin ribozyme appears to employ a unique strategy for RNA catalysis.     

A unique mechanism for RNA catalysis: the role of metal cofactors in hairpin ribozyme cleavage

Background: Ribozymes are biological catalysts that promote the hydrolysis and transesterification of phosphate diesters of RNA. They typically require divalent magnesium ions for activation, although it has proven difficult to differentiate structural from catalytic roles for the magnesium ions and to identify the molecular mechanism of catalysis. Direct inner-sphere coordination is usually invoked in the catalytic step, although there is no evidence to support the generality of such a pathway for all ribozymes.Results: We studied the catalytic pathway for the hairpin class of ribozyme. The substitutionally inert transition metal complex cobalt hexaammine [Co(NH₃)₆³⁺) was shown to be as active as Mg²⁺(aq) in promoting hairpin ribozyme activity, demonstrating that inner-sphere pathways are not used by this class of ribozyme. These results were confirmed by studies with R_P- and S_P-phosphorothioate substrate analogs which show a similar reactivity to that of the native substrate towards the magnesium-activated ribozyme. Monovalent cations enhance the activity of Co(NH₃)₆³⁺-promoted reactions, but inhibit Mg²⁺-activated catalysis, demonstrating a requirement for hydrated cations at several key sites in the ribozyme.Conclusions: These results provide clear support for a model of RNA catalysis that does not involve direct coordination of magnesium to the phosphate ester, nor activation of a bound water molecule. A mechanism in which catalysis is carried out by functional groups on the RNA ribozyme itself is possible; such functional groups are likely to have pK_a values that are appropriate for carrying out this catalysis. The metal cofactor would then serve to define the architecture of the catalytic pocket and contribute to the stabilization of transient species, as has been described earlier. Hydrolytic pathways in nucleic acid reactions are apparently more diverse than was previously thought, and the hairpin ribozyme falls into a mechanistically distinct class from the Tetrahymena and the hammerhead ribozymes.     

General base catalysis for cleavage by the active-site cytosine of the hepatitis delta virus ribozyme: QM/MM calculations establish chemical feasibility

The hepatitis delta virus (HDV) ribozyme is an RNA motif embedded in human pathogenic HDV RNA. Previous experimental studies have established that the active-site nucleotide C75 is essential for self-cleavage of the ribozyme, although its exact catalytic role in the process remains debated. Structural data from X-ray crystallography generally indicate that C75 acts as the general base that initiates catalysis by deprotonating the 2'-OH nucleophile at the cleavage site, while a hydrated magnesium ion likely protonates the 5'-oxygen leaving group. In contrast, some mechanistic studies support the role of C75 acting as general acid and thus being protonated before the reaction. We report combined quantum chemical/molecular mechanical calculations for the C75 general base pathway, utilizing the available structural data for the wild type HDV genomic ribozyme as a starting point. Several starting configurations differing in magnesium ion placement were considered and both one-dimensional and two-dimensional potential energy surface scans were used to explore plausible reaction paths. Our calculations show that C75 is readily capable of acting as the general base, in concert with the hydrated magnesium ion as the general acid. We identify a most likely position for the magnesium ion, which also suggests it acts as a Lewis acid. The calculated energy barrier of the proposed mechanism, approximately 20 kcal/mol, would lower the reaction barrier by approximately 15 kcal/mol compared with the uncatalyzed reaction and is in good agreement with experimental data.     

Capreomycin and hygromycin B modulate the catalytic activity of the delta ribozyme in a manner that depends on the protonation and complexation with Cu2+ ions of these antibiotics

Catalytic RNA molecules (ribozymes) have often been used for the testing of interactions of antibiotics with ribonucleic acids. We showed that the impact of capreomycin and hygromycin B on delta ribozyme catalysis might change dramatically, from stimulation to inhibition, depending on conditions. In order to evaluate possible mechanisms of modulation of the ribozyme catalytic activity we used our earlier data on species distribution for protonated forms of capreomycin and hygromycin B and their complexes with Cu(2+) ions at different pH values. We proposed that, upon inhibition, the protonated amino group of capreomycin was located in the ribozyme catalytic cleft interfering with binding catalytic Mg(2+). Such a mechanism was also supported by the results of ribozyme inhibition with capreomycin complexed with Cu(2+). The effects of stimulation of the delta ribozyme activity by capreomycin and hygromycin B were less pronounced than inhibition. Possibly, the amino functions of these antibiotics might be involved in a general acid-base catalysis performed by the ribozyme, acting as proton acceptors/donors.     

Noncanonical hydrogen bonding in nucleic acids. Benchmark evaluation of key base-phosphate interactions in folded RNA molecules using quantum-chemical calculations and molecular dynamics simulations

RNA molecules are stabilized by a wide range of noncanonical interactions that are not present in DNA. Among them, the recently classified base-phosphate (BPh) interactions belong to the most important ones. Twelve percent of nucleotides in the ribosomal crystal structures are involved in BPh interactions. BPh interactions are highly conserved and provide major constraints on RNA sequence evolution. Here we provide assessment of the energetics of BPh interactions using MP2 computations extrapolated to the complete basis set of atomic orbitals and corrected for higher-order electron correlation effects. The reference computations are compared with DFT-D and DFT-D3 approaches, the SAPT method, and the molecular mechanics force field. The computations, besides providing the basic benchmark for the BPh interactions, allow some refinements of the original classification, including identification of some potential doubly bonded BPh patterns. The reference computations are followed by analysis of some larger RNA fragments that consider the context of the BPh interactions. The computations demonstrate the complexity of interaction patterns utilizing the BPh interactions in real RNA structures. The BPh interactions are often involved in intricate interaction networks. We studied BPh interactions of protonated adenine that can contribute to catalysis of hairpin ribozyme, the key BPh interaction in the S-turn motif of the sarcin-ricin loop, which may predetermine the S-turn topology and complex BPh patterns from the glmS riboswitch. Finally, the structural stability of BPh interactions in explicit solvent molecular dynamics simulations is assessed. The simulations well preserve key BPh interactions and allow dissection of structurally/functionally important water-meditated BPh bridges, which could not be considered in earlier bioinformatics classification of BPh interactions.     

Two-metal-ion mechanism for hammerhead-ribozyme catalysis

The hammerhead ribozyme is one of the best studied ribozymes, but it still presents challenges for our understanding of RNA catalysis. It catalyzes a transesterification reaction that converts a 5',3' diester to a 2',3' cyclic phosphate diester via an S(N)2 mechanism. Thus, the overall reaction corresponds to that catalyzed by bovine pancreatic ribonuclease. However, an essential distinguishing aspect is that metal ions are not involved in RNase catalysis but appear to be important in ribozymes. Although various techniques have been used to assign specific functions to metals in the hammerhead ribozyme, their number and roles in catalysis is not clear. Two recent theoretical studies on RNA catalysis examined the reaction mechanism of a single-metal-ion model. A two-metal-ion model, which is supported by experiment and based on ab initio and density functional theory calculations, is described here. The proposed mechanism of the reaction has four chemical steps with three intermediates and four transition states along the reaction pathway. Reaction profiles are calculated in the gas phase and in solution. The early steps of the reaction are found to be fast (with low activation barriers), and the last step, corresponding to the departure of the leaving group, is rate limiting. This two-metal-ion model differs from the models proposed previously in that the two metal ions function not only as Lewis acids but also as general acids/bases. Comparison with experiment shows good agreement with thermodynamic and kinetic data. A detailed analysis based on natural bond orbitals (NBOs) and natural energy decomposition (NEDA) provides insights into the role of metal ions and other factors important for catalysis.     

Peptidyl-transferase ribozymes: trans reactions,structural characterization and ribosomal RNA-like features

Background: One of the most significant questions in understanding the origin of life concerns the order of appearance of DNA, RNA and protein during early biological evolution. If an ‘RNA world’ was a precursor to extant life, RNA must be able not only to catalyze RNA replication but also to direct peptide synthesis. Iterative Iterative RNA selection previously identified catalytic RNAs (ribozymes) that form amide bonds between RNA and an amino acid or between two amino acids.Results: We characterized peptidyl-transferase reactions catalyzed by two different families of ribozymes that use substrates that mimic A site and P site tRNAs. The family II ribozyme secondary structure was modeled using chemical modification, enzymatic digestion and mutational analysis. Two regions resemble the peptidyl-transferase region of 23S ribosomal RNA in sequence and structural context; these regions are important for peptide-bond formation. A shortened form of this ribozyme was engineered to catalyze intermolecular (‘trans’) peptide-bond formation, with the two amino-acid substrates binding through an attached AMP or oligonucleotide moiety.Conclusions: An in vitro-selected ribozyme can catalyze the same type of peptide-bond formation as a ribosome; the ribozyme resembles the ribosome because a very specific RNA structure is required for substrate binding and catalysis, and both amino acids are attached to nucleotides. It is intriguing that, although there are many different possible peptidyl-transferase ribozymes, the sequence and secondary structure of one is strikingly similar to the ‘helical wheel’ portion of 23S rRNA implicated in ribosomal peptidyl-transferase activity.     

Probing general base catalysis in the hammerhead ribozyme

Recent structural and computational studies have shed new light on the catalytic mechanism and active site structure of the RNA cleaving hammerhead ribozyme. Consequently, specific ribozyme functional groups have been hypothesized to be directly involved in general/acid base catalysis. In order to test this hypothesis, we have developed an affinity label to identify the functional general base in the S. mansoni hammerhead ribozyme. The ribozyme was reacted with a substrate analogue bearing a 2'-bromoacetamide group in place of the nucleophilic 2'-hydroxyl group which would normally be deprotonated by a general base. The electrophilic 2'-bromoacetamide group is poised to alkylate the general base, which is subsequently identified by footprinting analysis. Herein, we demonstrate alkylation of N1 of G12 in the hammerhead ribozyme in a pH and [Mg(2+)] dependent manner that is consistent with the native cleavage reaction. These results provide substantial evidence that deprotonated N1 of G12 functions directly as a general base in the hammerhead ribozyme; moreover, our experiments provide evidence that the pKa of G12 is perturbed downward in the context of the active site structure. We also observed other pH-independent alkylations, which do not appear to reflect the catalytic mechanism, but offer further insight into ribozyme conformation and structure.     

Constant pH Molecular Dynamics Simulations of Nucleic Acids in Explicit Solvent

The nucleosides of adenine and cytosine have pKa values of 3.50 and 4.08, respectively, and are assumed to be unprotonated under physiological conditions. However, evidence from recent NMR and X-Ray crystallography studies has revealed the prevalence of protonated adenine and cytosine in RNA macromolecules. Such nucleotides with elevated pKa values may play a role in stabilizing RNA structure and participate in the mechanism of ribozyme catalysis. With the work presented here, we establish the framework and demonstrate the first constant pH MD simulations (CPHMD) for nucleic acids in explicit solvent in which the protonation state is coupled to the dynamical evolution of the RNA system via λ-dynamics. We adopt the new functional form λ(Nexp) for λ that was recently developed for Multi-Site λ-Dynamics (MSλD) and demonstrate good sampling characteristics in which rapid and frequent transitions between the protonated and unprotonated states at pH = pKa are achieved. Our calculated pKa values of simple nucleotides are in a good agreement with experimentally measured values, with a mean absolute error of 0.24 pKa units. This work demonstrates that CPHMD can be used as a powerful tool to investigate pH-dependent biological properties of RNA macromolecules.     

A bridge between the RNA and protein worlds? Accelerating delivery of chemical reactivity to RNA and DNA by a specific short peptide (AAKK)(4)

BACKGROUND: RNA can catalyze diverse chemical reactions, leading to the hypothesis that an RNA world existed early in evolution. Today, however, catalysis by naturally occurring RNAs is rare and most chemical transformations within cells require proteins. This has led to interest in the design of small peptides capable of catalyzing chemical transformations. RESULTS: We demonstrate that a short lysine-rich peptide (AAKK)(4) can deliver a nucleophile to DNA or RNA and amplify the rate of chemical modification by up to 3400-fold. We also tested similar peptides that contain ornithine or arginine in place of lysine, peptides with altered stereochemistry or orientation, and peptides containing eight lysines but with different spacing. Surprisingly, these similar peptides function much less well, suggesting that specific combinations of amino acids, charge distribution, and stereochemistry are necessary for the rate enhancement by (AAKK)(4). CONCLUSIONS: By appending other reactive groups to (AAKK)(4) it should be possible to greatly expand the potential for small peptides to directly catalyze modification of DNA or RNA or to act as cofactors to promote ribozyme catalysis.     

Relationship between catalytic activity and secondary structure of a hammerhead ribozyme: A study using thermodynamic parameters for RNA structure prediction

Abstract

The stabilization energy for the secondary structures of wild-type hammerhead and mutant ribozymes has been calculated at different salt conditions and temperatures by using the thermodynamic parameters for RNA structure prediction. The most stable structure at each condition has been searched and the obtained secondary structure is compared with the structure suggested phylogenetically or experimentally. The results indicate that the hammerhead-type secondary structure of the ribozyme and its reactivity correlate with each other. The multibranched loop containing the self-cleavage site of the ribozyme particularly should be a key structure in the hammerhead ribozyme reaction. The predicted secondary structures also suggest that the reactivity of the hammerhead ribozyme should be very much lower at 10°C than that at 37°C.     

Rational design of allosteric ribozymes

Background: Efficient operation of cellular processes relies on the strict control that each cell exerts over its metabolic pathways. Some protein enzymes are subject to allosteric regulation, in which binding sites located apart from the enzyme's active site can specifically recognize effector molecules and alter the catalytic rate of the enzyme via conformational changes. Although RNA also performs chemical reactions, no ribozymes are known to operate as true allosteric enzymes in biological systems. It has recently been established that small-molecule receptors can readily be made of RNA, as demonstrated by the in vitro selection of various RNA aptamers that can specifically bind corresponding ligand molecules. We set out to examine whether the catalytic activity of an existing ribozyme could be brought under the control of an effector molecule by designing conjoined aptamer-ribozyme complexes.Results: By joining an ATP-binding RNA to a self-cleaving ribozyme, we have created the first example of an allosteric ribozyme that has a catalytic rate that can be controlled by ATP. A 180-fold reduction in rate is observed upon addition of either adenosine or ATP, but no inhibition is detected in the presence of dATP or other nucleoside triphosphates. Mutations in the aptamer domain that are expected to eliminate ATP binding or that increase the distance between aptamer and ribozyme domains result in a loss of ATP-specific allosteric control. Using a similar design approach, allosteric hammerhead ribozymes that are activated in the presence of ATP were created and another ribozyme that can be controlled by theophylline was created.Conclusions: The catalytic features of these conjoined aptamer-ribozyme constructs demonstrate that catalytic RNAs can also be subject to allosteric regulation — a key feature of certain protein enzymes. Moreover, by using simple rational design strategies, it is now possible to engineer new catalytic polynucleotides which have rates that can be tightly and specifically controlled by small effector molecules.     

Active site labeling of G8 in the hairpin ribozyme: implications for structure and mechanism

There is mounting evidence that suggests that general acid/base catalysis is operative in the hairpin ribozyme, with analogy to the protein enzyme RNaseA. Nevertheless, the extent of general base catalysis as well as the identity of the specific chemical groups responsible remains the subject of some controversy. An affinity label has previously been used to alkylate histidine 12 (His12), the active general base in RNaseA. To date, no such experiment has been applied to a ribozyme. We have synthesized the analogous affinity label for the hairpin ribozyme with an electrophilic 2'-bromoacetamide group in lieu of the 2'-hydroxyl (2'OH) at the substrate cleavage site and show that guanosine 8 (G8) of the hairpin ribozyme is specifically alkylated, most likely at the N1 position. This evidence strongly implicates N1 of G8 in active site chemistry. By direct analogy to RNase A, these findings could be consistent with the hypothesis that deprotonated G8 residue functions as a general base in the hairpin ribozyme. Other mechanistic possibilities for N1 of G8 such as indirect general base catalysis mediated by a water molecule or transition state stabilization could also be consistent with our findings.