Analysis of metal ion dependence in glmS ribozyme self-cleavage and coenzyme binding

Characterization of the glmS Ribozymes from Listeria Monocytogenes and Clostridium Difficile

The glmS ribozyme regulates the expression of the essential GlmS enzyme being involved in cell wall biosynthesis. While >450 variants of the glmS ribozyme were identified by in silico approaches and homology searches, only a few have yet been experimentally investigated. Herein, we validate and characterize the glmS ribozymes of the human pathogens Clostridium difficile and Listeria monocytogenes. Both ribozymes, as their previous characterized homologs rely on glucosamine-6-phosphate as co-factor and the presence of divalent cations for exerting the cleavage reaction. The observed EC₅₀ values in turn were found to be in the submicromolar range, at least an order of magnitude lower than observed for glmS ribozymes from other bacteria. The glmS ribozyme of L. monocytogenes was further shown to bear unique properties. It discriminates between co-factors very stringently and other than the glmS ribozyme of C. difficile retains activity at low temperatures. This finding illustrates that albeit being highly conserved, glmS ribozymes have unique characteristics.

Coordination Chemistry of Nucleotides and Antivirally Active Acyclic Nucleoside Phosphonates, including Mechanistic Considerations

Considering that practically all reactions that involve nucleotides also involve metal ions, it is evident that the coordination chemistry of nucleotides and their derivatives is an essential corner stone of biological inorganic chemistry. Nucleotides are either directly or indirectly involved in all processes occurring in Nature. It is therefore no surprise that the constituents of nucleotides have been chemically altered-that is, at the nucleobase residue, the sugar moiety, and also at the phosphate group, often with the aim of discovering medically useful compounds. Among such derivatives are acyclic nucleoside phosphonates (ANPs), where the sugar moiety has been replaced by an aliphatic chain (often also containing an ether oxygen atom) and the phosphate group has been replaced by a phosphonate carrying a carbon-phosphorus bond to make the compounds less hydrolysis-sensitive. Several of these ANPs show antiviral activity, and some of them are nowadays used as drugs. The antiviral activity results from the incorporation of the ANPs into the growing nucleic acid chain-i.e., polymerases accept the ANPs as substrates, leading to chain termination because of the missing 3'-hydroxyl group. We have tried in this review to describe the coordination chemistry (mainly) of the adenine nucleotides AMP and ATP and whenever possible to compare it with that of the dianion of 9-[2-(phosphonomethoxy)ethyl]adenine (PMEA2- = adenine(N9)-CH2-CH2-O-CH2-PO32) [or its diphosphate (PMEApp4-)] as a representative of the ANPs. Why is PMEApp4- a better substrate for polymerases than ATP4-? There are three reasons: (i) PMEA2- with its anti-like conformation (like AMP2-) fits well into the active site of the enzyme. (ii) The phosphonate group has an enhanced metal ion affinity because of its increased basicity. (iii) The ether oxygen forms a 5-membered chelate with the neighboring phosphonate and favors thus coordination at the Pα group. Research on ANPs containing a purine residue revealed that the kind and position of the substituent at C2 or C6 has a significant influence on the biological activity. For example, the shift of the (C6)NH2 group in PMEA to the C2 position leads to 9-[2-(phosphonomethoxy)ethyl]-2-aminopurine (PME2AP), an isomer with only a moderate antiviral activity. Removal of (C6)NH2 favors N7 coordination, e.g., of Cu2+, whereas the ether O atom binding of Cu2+ in PMEA facilitates N3 coordination via adjacent 5- and 7-membered chelates, giving rise to a Cu(PMEA)cl/O/N3 isomer. If the metal ions (M2+) are M(α,β)-M(γ)-coordinated at a triphosphate chain, transphosphorylation occurs (kinases, etc.), whereas metal ion binding in a M(α)-M(β,γ)-type fashion is relevant for polymerases. It may be noted that with diphosphorylated PMEA, (PMEApp4-), the M(α)-M(β,γ) binding is favored because of the formation of the 5-membered chelate involving the ether O atom (see above). The self-association tendency of purines leads to the formation of dimeric [M2(ATP)]2(OH)- stacks, which occur in low concentration and where one half of the molecule undergoes the dephosphorylation reaction and the other half stabilizes the structure-i.e., acts as the "enzyme" by bridging the two ATPs. In accord herewith, one may enhance the reaction rate by adding AMP2- to the [Cu2(ATP)]2(OH)- solution, as this leads to the formation of mixed stacked Cu3(ATP)(AMP)(OH)- species, in which AMP2- takes over the structuring role, while the other "half" of the molecule undergoes dephosphorylation. It may be added that Cu3(ATP)(PMEA) or better Cu3(ATP)(PMEA)(OH)- is even a more reactive species than Cu3(ATP)(AMP)(OH)-. - The matrix-assisted self-association and its significance for cell organelles with high ATP concentrations is summarized and discussed, as is, e.g., the effect of tryptophanate (Trp-), which leads to the formation of intramolecular stacks in M(ATP)(Trp)3- complexes (formation degree about 75%). Furthermore, it is well-known that in the active-site cavities of enzymes the dielectric constant, compared with bulk water, is reduced; therefore, we have summarized and discussed the effect of a change in solvent polarity on the stability and structure of binary and ternary complexes: Opposite effects on charged O sites and neutral N sites are observed, and this leads to interesting insights.

The GlcN6P cofactor plays multiple catalytic roles in the glmS ribozyme

RNA enzymes (ribozymes) have remarkably diverse biological roles despite having limited chemical diversity. Protein enzymes enhance their reactivity through recruitment of cofactors; likewise, the naturally occurring glmS ribozyme uses the glucosamine-6-phosphate (GlcN6P) organic cofactor for phosphodiester bond cleavage. Prior structural and biochemical studies have implicated GlcN6P as the general acid. Here we describe new catalytic roles of GlcN6P through experiments and calculations. Large stereospecific normal thio effects and a lack of metal-ion rescue in the holoribozyme indicate that nucleobases and the cofactor play direct chemical roles and align the active site for self-cleavage. Large stereospecific inverse thio effects in the aporibozyme suggest that the GlcN6P cofactor disrupts an inhibitory interaction of the nucleophile. Strong metal-ion rescue in the aporibozyme reveals that this cofactor also provides electrostatic stabilization. Ribozyme organic cofactors thus perform myriad catalytic roles, thereby allowing RNA to compensate for its limited functional diversity.

A divalent cation-dependent variant of the glmS ribozyme with stringent Ca2+ selectivity co-opts a preexisting nonspecific metal ion-binding site

Ribozymes use divalent cations for structural stabilization, as catalytic cofactors, or both. Because of the prominent role of Ca²⁺ in intracellular signaling, engineered ribozymes with stringent Ca²⁺ selectivity would be important in biotechnology. The wild-type glmS ribozyme (glmS^WT) requires glucosamine-6-phosphate (GlcN6P) as a catalytic cofactor. Previously, a glmS ribozyme variant with three adenosine mutations (glmS^AAA) was identified, which dispenses with GlcN6P and instead uses, with little selectivity, divalent cations as cofactors for site-specific RNA cleavage. We now report a Ca²⁺-specific ribozyme (glmS^Ca) evolved from glmS^AAA that is >10,000 times more active in Ca²⁺ than Mg²⁺, is inactive in even 100 mM Mg²⁺, and is not responsive to GlcN6P. This stringent selectivity, reminiscent of the protein nuclease from Staphylococcus, allows rapid and selective ribozyme inactivation using a Ca²⁺ chelator such as EGTA. Because glmS^Ca functions in physiologically relevant Ca²⁺ concentrations, it can form the basis for intracellular sensors that couple Ca²⁺ levels to RNA cleavage. Biochemical analysis of glmS^Ca reveals that it has co-opted for selective Ca²⁺ binding a nonspecific cation-binding site responsible for structural stabilization in glmS^WT and glmS^AAA Fine-tuning of the selectivity of the cation site allows repurposing of this preexisting molecular feature.

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

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

Two Divalent Metal Ions and Conformational Changes Play Roles in the Hammerhead Ribozyme Cleavage Reaction

The hammerhead ribozyme is a self-cleaving RNA broadly dispersed across all kingdoms of life. Although it was the first of the small, nucleolytic ribozymes discovered, the mechanism by which it catalyzes its reaction remains elusive. The nucleobase of G12 is well positioned to be a general base, but it is unclear if or how this guanine base becomes activated for proton transfer. Metal ions have been implicated in the chemical mechanism, but no interactions between divalent metal ions and the cleavage site have been observed crystallographically. To better understand how this ribozyme functions, we have solved crystal structures of wild-type and G12A mutant ribozymes. We observe a pH-dependent conformational change centered around G12, consistent with this nucleotide becoming deprotonated. Crystallographic and kinetic analysis of the G12A mutant reveals a Zn(2+) specificity switch suggesting a direct interaction between a divalent metal ion and the purine at position 12. The metal ion specificity switch and the pH-rate profile of the G12A mutant suggest that the minor imino tautomer of A12 serves as the general base in the mutant ribozyme. We propose a model in which the hammerhead ribozyme rearranges prior to the cleavage reaction, positioning two divalent metal ions in the process. The first metal ion, positioned near G12, becomes directly coordinated to the O6 keto oxygen, to lower the pKa of the general base and organize the active site. The second metal ion, positioned near G10.1, bridges the N7 of G10.1 and the scissile phosphate and may participate directly in the cleavage reaction.

Role of the active site guanine in the glmS ribozyme self-cleavage mechanism: quantum mechanical/molecular mechanical free energy simulations

The glmS ribozyme catalyzes a self-cleavage reaction at the phosphodiester bond between residues A-1 and G1. This reaction is thought to occur by an acid-base mechanism involving the glucosamine-6-phosphate cofactor and G40 residue. Herein quantum mechanical/molecular mechanical free energy simulations and pKa calculations, as well as experimental measurements of the rate constant for self-cleavage, are utilized to elucidate the mechanism, particularly the role of G40. Our calculations suggest that an external base deprotonates either G40(N1) or possibly A-1(O2'), which would be followed by proton transfer from G40(N1) to A-1(O2'). After this initial deprotonation, A-1(O2') starts attacking the phosphate as a hydroxyl group, which is hydrogen-bonded to deprotonated G40, concurrent with G40(N1) moving closer to the hydroxyl group and directing the in-line attack. Proton transfer from A-1(O2') to G40 is concomitant with attack of the scissile phosphate, followed by the remainder of the cleavage reaction. A mechanism in which an external base does not participate, but rather the proton transfers from A-1(O2') to a nonbridging oxygen during nucleophilic attack, was also considered but deemed to be less likely due to its higher effective free energy barrier. The calculated rate constant for the favored mechanism is in agreement with the experimental rate constant measured at biological Mg(2+) ion concentration. According to these calculations, catalysis is optimal when G40 has an elevated pKa rather than a pKa shifted toward neutrality, although a balance among the pKa's of A-1, G40, and the nonbridging oxygen is essential. These results have general implications, as the hammerhead, hairpin, and twister ribozymes have guanines at a similar position as G40.

Metabolite-binding ribozymes

Catalysis in the biological context was largely thought to be a protein-based phenomenon until the discovery of RNA catalysts called ribozymes. These discoveries demonstrated that many RNA molecules exhibit remarkable structural and functional versatility. By virtue of these features, naturally occurring ribozymes have been found to be involved in catalyzing reactions for fundamentally important cellular processes such as translation and RNA processing. Another class of RNAs called riboswitches directly binds ligands to control downstream gene expression. Most riboswitches regulate downstream gene expression by controlling premature transcription termination or by affecting the efficiency of translation initiation. However, one riboswitch class couples ligand-sensing to ribozyme activity. Specifically, the glmS riboswitch is a nucleolytic ribozyme, whose self-cleavage activity is triggered by the binding of GlcN6P. The products of this self-cleavage reaction are then targeted by cellular RNases for rapid degradation, thereby reducing glmS expression under conditions of sufficient GlcN6P. Since the discovery of the glmS ribozyme, other metabolite-binding ribozymes have been identified. Together, these discoveries have expanded the general understanding of noncoding RNAs and provided insights that will assist future development of synthetic riboswitch-ribozymes. A very broad overview of natural and synthetic ribozymes is presented herein with an emphasis on the structure and function of the glmS ribozyme as a paradigm for metabolite-binding ribozymes that control gene expression. This article is part of a Special Issue entitled: Riboswitches.

An in vitro evolved glmS ribozyme has the wild-type fold but loses coenzyme dependence

Uniquely among known ribozymes, the glmS ribozyme-riboswitch requires a small-molecule coenzyme, glucosamine-6-phosphate (GlcN6P). Although consistent with its gene-regulatory function, the use of GlcN6P is unexpected because all of the other characterized self-cleaving ribozymes use RNA functional groups or divalent cations for catalysis. To determine what active site features make this ribozyme reliant on GlcN6P and to evaluate whether it might have evolved from a coenzyme-independent ancestor, we isolated a GlcN6P-independent variant through in vitro selection. Three active site mutations suffice to generate a highly reactive RNA that adopts the wild-type fold but uses divalent cations for catalysis and is insensitive to GlcN6P. Biochemical and crystallographic comparisons of wild-type and mutant ribozymes show that a handful of functional groups fine-tune the RNA to be either coenzyme or cation dependent. These results indicate that a few mutations can confer new biochemical activities on structured RNAs. Thus, families of structurally related ribozymes with divergent function may exist.

Metabolite recognition principles and molecular mechanisms underlying riboswitch function

Riboswitches are mRNA elements capable of modulating gene expression in response to specific binding by cellular metabolites. Riboswitches exert their function through the interplay of alternative ligand-free and ligand-bound conformations of the metabolite-sensing domain, which in turn modulate the formation of adjacent gene expression controlling elements. X-ray crystallography and NMR spectroscopy have determined three-dimensional structures of virtually all the major riboswitch classes in the ligand-bound state and, for several riboswitches, in the ligand-free state. The resulting spatial topologies have demonstrated the wide diversity of riboswitch folds and revealed structural principles for specific recognition by cognate metabolites. The available three-dimensional information, supplemented by structure-guided biophysical and biochemical experimentation, has led to an improved understanding of how riboswitches fold, what RNA conformations are required for ligand recognition, and how ligand binding can be transduced into gene expression modulation. These studies have greatly facilitated the dissection of molecular mechanisms underlying riboswitch action and should in turn guide the anticipated development of tools for manipulating gene regulatory circuits.

The structural and functional uniqueness of the glmS ribozyme

The glmS bacterial ribozyme/riboswitch is found in a number of Gram-positive bacteria, many of which are human pathogens. Investigation of the structure and function of the glmS catalyst will aid in the development of artificial agonists/antagonists that might function as novel antibiotics. The glmS ribozyme is mechanistically unique in that it is the first RNA catalyst identified to require a coenzyme, glucosamine-6-phosphate, for RNA self-cleavage. In addition, it is the first riboswitch identified to utilize self-cleavage as a mode of genetic regulation in metabolism. Significant biochemical and biophysical data exist for the glmS ribozyme and aid in mechanistically understanding the importance of RNA and coenzyme structure to function in acid-base catalysis.

Characterization of metal ion-nucleic acid interactions in solution

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

Mechanism and distribution of glmS ribozymes

Among the nine classes of ribozymes that have been experimentally validated to date is the metabolite-responsive self-cleaving ribozyme called glmS. This RNA is almost exclusively located in the 5'-untranslated region of bacterial mRNAs that code for the production of GlmS proteins, which catalyze the synthesis of the aminosugar glucosamine-6-phosphate (GlcN6P). Each glmS ribozyme forms a conserved catalytic core that selectively binds GlcN6P and uses this metabolite as a cofactor to promote ribozyme self-cleavage. Metabolite-induced self-cleavage results in down-regulation of glmS gene expression, and thus the ribozyme functions as a key riboswitch component to permit feedback regulation of GlcN6P levels. Representatives of glmS ribozymes also serve as excellent experimental models to elucidate how RNAs fold to recognize small molecule ligands and promote chemical transformations.

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.

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

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

An expanded collection and refined consensus model of glmS ribozymes

Self-cleaving glmS ribozymes selectively bind glucosamine-6-phosphate (GlcN6P) and use this metabolite as a cofactor to promote self-cleavage by internal phosphoester transfer. Representatives of the glmS ribozyme class are found in Gram-positive bacteria where they reside in the 5' untranslated regions (UTRs) of glmS messenger RNAs that code for the essential enzyme L-glutamine:D-fructose-6-phosphate aminotransferase. By using comparative sequence analyses, we have expanded the number of glmS ribozyme representatives from 160 to 463. All but two glmS ribozymes are present in glmS mRNAs and most exhibit striking uniformity in sequence and structure, which are features that make representatives attractive targets for antibacterial drug development. However, our discovery of rare variants broadens the consensus sequence and structure model. For example, in the Deinococcus-Thermus phylum, several structural variants exist that carry additional stems within the catalytic core and changes to the architecture of core-supporting substructures. These findings reveal that glmS ribozymes have a broader phylogenetic distribution than previously known and suggest that additional rare structural variants may remain to be discovered.

Rapid steps in the glmS ribozyme catalytic pathway: cation and ligand requirements

The glmS ribozyme is a conserved riboswitch found in numerous Gram-positive bacteria and responds to the cellular concentrations of glucosamine 6-phosphate (GlcN6P). GlcN6P binding promotes site-specific self-cleavage in the 5' UTR of the glmS mRNA, resulting in downregulation of gene expression. The glmS ribozyme has previously been shown to lack strong cation specificity when the rate-limiting folding step of the cleavage reaction pathway is measured. This does not provide data regarding cation and ligand specificities of the glmS ribozyme during the rapid ligand binding chemical catalysis events. Prefolding of the ribozyme in Mg(2+)-containing buffers effectively isolates the rapid ligand binding and catalytic events (k(obs) > 60 min(-1)) from rate-limiting folding (k(obs) < 4 min(-1)). Here we employ this experimental design to assay the cations and ligand requirements for rapid ligand binding and catalysis. We show that molar concentrations of monovalent cations are also capable of inducing the formation of the native GlcN6P binding structure but are unable to promote ligand binding and catalysis rates of >4 min(-1). Our data show that the sole obligatory role for divalent cations, for which there is crystallographic evidence, is coordination of the phosphate moiety of GlcN6P in the ligand-binding pocket. In further support of this hypothesis, our data show that a nonphosphorylated analogue of GlcN6P, glucosamine, is unable to promote rapid ligand binding and catalysis in the presence of divalent cations. Folding of the ribozyme is, therefore, relatively independent of cation identity, but the rapid initiation of catalysis upon the addition of ligand is stricter.