Ligand requirements for glmS ribozyme self-cleavage

Natural RNA catalysts (ribozymes) perform essential reactions in biological RNA processing and protein synthesis, whereby catalysis is intrinsic to RNA structure alone or in combination with metal ion cofactors. The recently discovered glmS ribozyme is unique in that it functions as a glucosamine-6-phosphate (GlcN6P)-dependent catalyst believed to enable "riboswitch" regulation of amino-sugar biosynthesis in certain prokaryotes. However, it is unclear whether GlcN6P functions as an effector or coenzyme to promote ribozyme self-cleavage. Herein, we demonstrate that ligand is absolutely requisite for glmS ribozyme self-cleavage activity. Furthermore, catalysis both requires and is dependent upon the acid dissociation constant (pKa) of the amine functionality of GlcN6P and related compounds. The data demonstrate that ligand is integral to catalysis, consistent with a coenzyme role for GlcN6P and illustrating an expanded capacity for biological RNA catalysis.     

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.     

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.     

Direct pK(a) measurement of the active-site cytosine in a genomic hepatitis delta virus ribozyme

Hepatitis delta virus ribozymes have been proposed to perform self-cleavage via a general acid/base mechanism involving an active-site cytosine, based on evidence from both a crystal structure of the cleavage product and kinetic measurements. To determine whether this cytosine (C75) in the genomic ribozyme has an altered pK(a) consistent with its role as a general acid or base, we used (13)C NMR to determine its microscopic pK(a) in the product form of the ribozyme. The measured pK(a) is moderately shifted from that of a free nucleoside or a base-paired cytosine and has the same divalent metal ion dependence as the apparent reaction pK(a)'s measured kinetically. However, under all conditions tested, the microscopic pK(a) is lower than the apparent reaction pK(a), supporting a model in which C75 is deprotonated in the product form of the ribozyme at physiological pH. While additional results suggest that the pK(a) is not shifted in the reactant state of the ribozyme, these data cannot rule out elevation of the C75 pK(a) in an intermediate state of the transesterification reaction.     

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.     

Restoration of Ribozyme Tertiary Contact and Function by Using a Molecular Glue for RNA

Some RNA classes require folding into the proper higher‐order structures to exert their functions. Hammerhead ribozyme (HHR) requires a folding conformation stabilized by tertiary interaction for full activity. A rationally engineered HHR was developed that was inactive, but could be activated by a synthetic RNA‐binding ligand, naphthyridine carbamate tetramer with Z‐stilbene linker (Z‐NCTS). Binding of Z‐NCTS could induce the formation of an active folding structure and thereby restore ribozyme activity, where Z‐NCTS acts as a molecular glue to bring two isolated RNA loops into contact with each other. Next, we designed a Z‐NCTS‐responsive genetic switch using the HHR sequence inserted into the 3′ untranslated region as a cis‐acting element. We demonstrated that the rationally designed ribozyme switch enabled regulation of gene expression by Z‐NCTS and was functional in mammalian cells.     

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.     

A fluorine scan of the phenylamidinium needle of tricyclic thrombin inhibitors: effects of fluorine substitution on pKa and binding affinity and evidence for intermolecular C-F...CN interactions

The H-atoms of the phenylamidinium needle of tricyclic thrombin inhibitors, which interacts with Asp189 at the bottom of the selectivity pocket S1 of the enzyme, were systematically exchanged with F-atoms in an attempt to improve the pharmacokinetic properties by lowering the pK(a) value. Both the pK(a) values and the inhibitory constants K(i) against thrombin and trypsin were decreased upon F-substitution. Interestingly, linear free energy relationships (LFERs) revealed that binding affinity against thrombin is much more affected by a decrease in pK(a) than the affinity against trypsin. Surprising effects of F-substitutions in the phenylamidinium needle on the pK(a) value of the tertiary amine centre in the tricyclic scaffold of the inhibitors were observed and subsequently rationalised by X-ray crystallographic analysis and ab initio calculations. Evidence for highly directional intermolecular C-F...CN interactions was obtained by analysis of small-molecule X-ray crystal structures and investigations in the Cambridge Structural Database (CSD).     

2'-Fluoro substituents can mimic native 2'-hydroxyls within structured RNA

The ability of fluorine in a C-F bond to act as?a hydrogen bond acceptor is controversial. To test such ability in complex RNA macromolecules, we have replaced native 2'-OH groups with 2'-F and 2'-H groups in two related systems, the Tetrahymena group I ribozyme and the ΔC209 P4-P6 RNA domain. In three cases the introduced 2'-F mimics the native 2'-OH group, suggesting that the fluorine atom can accept a hydrogen bond. In each of these cases the native hydroxyl group interacts with a purine exocyclic amine. Our results give insight about the properties of organofluorine and suggest a possible general biochemical signature for tertiary interactions between 2'-hydroxyl groups and exocyclic amino groups within RNA.     

Fast cleavage kinetics of a natural hammerhead ribozyme

The hammerhead ribozyme is a small RNA motif that catalyzes the cleavage and ligation of RNA. The well-studied minimal hammerhead motif is inactive under physiological conditions and requires high Mg(2+) concentrations for efficient cleavage. In contrast, natural hammerheads are active under physiological conditions and contain motifs outside the catalytic core that lower the requirement for Mg(2+). Single-turnover kinetics were used here to characterize the Mg(2+) and pH dependence for cleavage of a trans-cleaving construct of the Schistosoma mansoni natural hammerhead ribozyme. Compared to the minimal hammerhead motif, the natural Schistosoma ribozyme requires 100-fold less Mg(2+) to achieve a cleavage rate of 1 min(-1). The improved catalysis results from tertiary interactions between loops in stems I and II and likely arises from increasing the population of the active conformation. Under optimum pH and Mg(2+) conditions this ribozyme cleaves at over 870 min(-1) at 25 degrees C, further demonstrating the impressive catalytic power of this ribozyme.     

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.     

Cooperative binding and self-assembling behavior of cationic low molecular-weight dendrons with RNA molecules

Tri(ethylene glycol) derived, low molecular-weight dendrons with various amine end groups were synthesized and characterized for their properties of binding and self-assembling with RNA using the Candida ribozyme as a model RNA molecule. These dendritic compounds form stable complexes and well-defined nanoscale particles with RNA molecules via electrostatic interactions and self-assembly process, while leaving the other terminal of the tri(ethylene glycol) chain accessible for targeting. This suggests that dendrimers of this type hold great promise for specific RNA targeting and RNA delivery.     

Transmission of electronic effects through the {closo-1-CB9} and {closo-1-CB11} cages: apparent dissociation constants for series of [closo-1-CB9H8-1-COOH-10-X] and [closo-1-CB11H10-1-COOH-12-X] acids

The apparent ionization constants pK(a)' for series of carboxylic acids [closo-1-CB(9)H(8)-1-COOH-10-X](-) (1) and [closo-1-CB(11)H(10)-1-COOH-12-X](-) (2), where X = H, I, n-C(6)H(13), (+)NMe(3), (+)N(2), (+)SMe(2), OC(5)H(11), were measured in EtOH/H(2)O (1/1, v/v) at 24 °C. Correlation analysis of the pK(a)' values using Hammett substituent constants σ(p)(X) gave the reaction constant ρ = 0.87 ± 0.04 for series 1 and ρ = 1.00 ± 0.09 for series 2. These values are higher than for derivatives of PhCH═CHCOOH (ρ = 0.70 ± 0.09 in 55% EtOH) and correspond to 56% and 65% efficiencies in transmission of electronic effects by [closo-1-CB(9)H(10)](-) (E) and [closo-1-CB(11)H(12)](-) (F), respectively, as compared to benzene (A). Experimental results were supported with DFT calculations of relative acidity for series of acids derived from A, E, and F in aqueous medium.     

Direct measurement of a pK(a) near neutrality for the catalytic cytosine in the genomic HDV ribozyme using Raman crystallography

The hepatitis delta virus (HDV) ribozyme uses a cytosine to facilitate general acid-base catalysis. Biochemical studies suggest that C75 has a pKa perturbed to near neutrality. To measure this pKa directly, Raman spectra were recorded on single ribozyme crystals using a Raman microscope. A spectral feature arising from a single neutral cytosine was identified at 1528 cm(-1). At low pH, this mode was replaced with a new spectral feature. Monitoring these features as a function of pH revealed pKa values for the cytosine that couple anticooperatively with Mg2+ binding, with values of 6.15 and 6.40 in the presence of 20 and 2 mM Mg2+, respectively. These pKa values agree well with those obtained from ribozyme activity experiments in solution. To correlate the observed pKa with a specific nucleotide, crystals of C75U, which is catalytically inactive, were examined. The Raman difference spectra show that this mutation does not affect the conformation of the ribozyme. However, crystals of C75U did not produce a signal from a protonatable cytosine, providing strong evidence that protonation of C75 is being monitored in the wild-type ribozyme. These studies provide the first direct physical measurement of a pKa near neutrality for a catalytic residue in a ribozyme and show that ribozymes, like their protein enzyme counterparts, can optimize the pKa of their side chains for proton transfer.     

Dissociation constants of thymolphthalexone

For the dissociation constants of thymolphthalexone the following values have been found: pK(3) = 7.03 +/- 0.02, pK(4) = 8.05 +/- 0.09 (by potentiometric titration), pK(5) = 10.83 +/- 0.10, pK(6) = 12.99 +/- 0.11 (by spectrophotometry). They were determined at I = 0.4 and at 25 degrees.     

Theoretical evaluation of pK(a) in phosphoranes: implications for phosphate ester hydrolysis

Knowledge of the pK(a) of phosphoranes is important for the interpretation of phosphate ester hydrolysis. Calculated pK(a)'s of the model phosphorane, ethylene phosphorane, are reported. The method of calculation is based on the use of dimethyl phosphate as a reference state for evaluating relative pK(a) values, and on the optimization of the oxygen and acidic hydrogen van der Waals radii to give reasonable pK(1)(a), pK(2)(a), and pK(3)(a) for phosphoric acid in solution. Density functional theory is employed to calculate the gas-phase protonation energies, and continuum dielectric methods are used to determine the solvation corrections. The calculated pK(1)(a) and p(2)(a) for the model phosphorane are 7.9 and 14.3, respectively. These values are within the range of proposed experimental values, 6.5-11.0 for pK(1)(a), and 11.3-15.0 for pK(2)(a). The mechanistic implications of the calculated pK(a)'s are discussed.     

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.     

Highly non-planar dendritic porphyrin for pH sensing: observation of porphyrin monocation

Metal-free porphyrin-dendrimers provide a convenient platform for the construction of membrane-impermeable ratiometric probes for pH measurements in compartmentalized biological systems. In all previously reported molecules, electrostatic stabilization (shielding) of the core porphyrin by peripheral negative charges (carboxylates) was required to shift the intrinsically low porphyrin protonation pK(a)'s into the physiological pH range (pH 6-8). However, binding of metal cations (e.g., K(+), Na(+), Ca(2+), Mg(2+)) by the carboxylate groups on the dendrimer could affect the protonation behavior of such probes in biological environments. Here we present a dendritic pH nanoprobe based on a highly non-planar tetraaryltetracyclohexenoporphyrin (Ar(4)TCHP), whose intrinsic protonation pK(a)'s are significantly higher than those of regular tetraarylporphyrins, thereby eliminating the need for electrostatic core shielding. The porphyrin was modified with eight Newkome-type dendrons and PEGylated at the periphery, rendering a neutral water-soluble probe (TCHpH), suitable for measurements in the physiological pH range. The protonation of TCHpH could be followed by absorption (e.g., ε(Soret)(dication)～270,000 M(-1) cm(-1)) or by fluorescence. Unlike most tetraarylporphyrins, TCHpH is protonated in two distinct steps (pK(a)'s 7.8 and 6.0). In the region between the pK(a)'s, an intermediate species with a well-defined spectroscopic signature, presumably a TCHpH monocation, could be observed in the mixture. The performance of TCHpH was evaluated by pH gradient measurements in large unilamellar vesicles. The probe was retained inside the vesicles and did not pass through and/or interact with vesicle membranes, proving useful for quantification of proton transport across phospholipid bilayers. To interpret the protonation behavior of TCHpH we developed a model relating structural changes on the porphyrin macrocycle upon protonation to its basicity. The model was validated by density functional theory (DFT) calculations performed on a planar and non-planar porphyrin, making it possible to rationalize higher protonation pK(a)'s of non-planar porphyrins as well as the easier observation of their monocations.     

Nonenzymatic breakdown of the tetrahedral (alpha-carboxyketal phosphate) intermediates of MurA and AroA, two carboxyvinyl transferases. Protonation of different functional groups controls the rate and fate of breakdown

The mechanisms of nonenzymatic breakdown of the tetrahedral intermediates (THIs) of the carboxyvinyl transferases MurA and AroA were examined in order to illuminate the interplay between the inherent reactivities of the THIs and the enzymatic strategies used to promote catalysis. THI degradation was through phosphate departure, with C-O bond cleavage. It was acid catalyzed and dependent on the protonation state of the carboxyl of the alpha-carboxyketal phosphate functionality, with ionizations at pK(a) = 3.2 +/- 0.1 and 4.3 +/- 0.1 for MurA and AroA THIs, respectively. The solvent deuterium kinetic isotope effect for MurA THI at pL 2.0 was 1.3 +/- 0.4, consistent with general acid catalysis. The pK(a)'s suggested intramolecular general acid catalysis through protonation of the bridging oxygen of the phosphate, though H(3)O(+) catalysis was also possible. The product distribution varied with pH. The dominant breakdown products were pyruvate + phosphate + R-OH (R-OH = UDP-GlcNAc or shikimate 3-phosphate) at all pH's, particularly low pH. At higher pH's, increasing proportions of ketal, arising from intramolecular substitution of phosphate by the adjacent hydroxyl and the enolpyruvyl products of phosphate elimination were observed. With MurA THI, the product distribution fitted to pK(a)'s 1.6 and 6.2, corresponding to the expected pK(a)'s of a phosphate monoester. C-O bond cleavage was demonstrated by the lack of monomethyl [(33)P]phosphate formed upon degrading MurA [(33)P]THI in 50% methanol. General acid catalysis through the bridging oxygen is consistent with the location of the previously proposed general acid catalyst for THI breakdown in AroA, Lys22.