The ionic environment determines ribozyme cleavage rate by modulation of nucleobase pK a

Several small ribozymes employ general acid–base catalysis as a mechanism to enhance site-specific RNA cleavage, even though the functional groups on the ribonucleoside building blocks of RNA have pK_a values far removed from physiological pH. The rate of the cleavage reaction is strongly affected by the identity of the metal cation present in the reaction solution; however, the mechanism(s) by which different cations contribute to rate enhancement has not been determined. Using the Neurospora VS ribozyme, we provide evidence that different cations confer particular shifts in the apparent pK_a values of the catalytic nucleobases, which in turn determines the fraction of RNA in the protonation state competent for general acid–base catalysis at a given pH, which determines the observed rate of the cleavage reaction. Despite large differences in observed rates of cleavage in different cations, mathematical models of general acid–base catalysis indicate that k₁, the intrinsic rate of the bond-breaking step, is essentially constant irrespective of the identity of the cation(s) in the reaction solution. Thus, in contrast to models that invoke unique roles for metal ions in ribozyme chemical mechanisms, we find that most, and possibly all, of the ion-specific rate enhancement in the VS ribozyme can be explained solely by the effect of the ions on nucleobase pK_a. The inference that k₁ is essentially constant suggests a resolution of the problem of kinetic ambiguity in favor of a model in which the lower pK_a is that of the general acid and the higher pK_a is that of the general base.     

Role of an active site guanine in hairpin ribozyme catalysis probed by exogenous nucleobase rescue

The hairpin ribozyme is a small catalytic RNA with reversible phosphodiester cleavage activity. Biochemical and structural studies exclude a requirement for divalent metal cation cofactors and implicate one active site nucleobase in particular, G8, in the catalytic mechanism. Our previous work demonstrated that the cleavage activity that is lost when G8 is replaced by an abasic residue is restored when certain nucleobases are provided in solution. The specificity and pH dependence of exogenous nucleobase rescue were consistent with several models of the rescue mechanism, including general acid base catalysis, electrostatic stabilization of negative charge in the transition state or a requirement for protonation to facilitate exogenous nucleobase binding. Detailed analyses of exogenous nucleobase rescue for both cleavage and ligation reactions now allow us to refine models of the rescue mechanism. Activity increased with increasing pH for both unmodified ribozyme reactions and unrescued reactions of abasic variants lacking G8. This similarity in pH dependence argues against a role for G8 as a general base catalyst, because G8 deprotonation could not be responsible for the pH-dependent transition in the abasic variant. Exogenous nucleobase rescue of both cleavage and ligation activity increased with decreasing pH, arguing against a role for rescuing nucleobases in general acid catalysis, because a nucleobase that contributes general acid catalysis in the cleavage pathway should provide general base catalysis in ligation. Analysis of the concentration dependence of cytosine rescue at high and low pH demonstrated that protonation promotes catalysis within the nucleobase-bound ribozyme complex but does not stabilize nucleobase binding in the ground state. These results support an electrostatic stabilization mechanism in which exogenous nucleobase binding counters negative charge that develops in the transition state.     

Role of SLV in SLI substrate recognition by the Neurospora VS ribozyme

Substrate recognition by the VS ribozyme involves a magnesium-dependent loop/loop interaction between the SLI substrate and the SLV hairpin from the catalytic domain. Recent NMR studies of SLV demonstrated that magnesium ions stabilize a U-turn loop structure and trigger a conformational change for the extruded loop residue U700, suggesting a role for U700 in SLI recognition. Here, we kinetically characterized VS ribozyme mutants to evaluate the contribution of U700 and other SLV loop residues to SLI recognition. To help interpret the kinetic data, we structurally characterized the SLV mutants by NMR spectroscopy and generated a three-dimensional model of the SLI/SLV complex by homology modeling with MC-Sym. We demonstrated that the mutation of U700 by A, C, or G does not significantly affect ribozyme activity, whereas deletion of U700 dramatically impairs this activity. The U700 backbone is likely important for SLI recognition, but does not appear to be required for either the structural integrity of the SLV loop or for direct interactions with SLI. Thus, deletion of U700 may affect other aspects of SLI recognition, such as magnesium ion binding and SLV loop dynamics. As part of our NMR studies, we developed a convenient assay based on detection of unusual (31)P and (15)N N7 chemical shifts to probe the formation of U-turn structures in RNAs. Our model of the SLI/SLV complex, which is compatible with biochemical data, leads us to propose novel interactions at the loop I/loop V interface.     

Mechanisms of RNA catalysis

Ribozymes are RNA molecules that act as chemical catalysts. In contemporary cells, most known ribozymes carry out phosphoryl transfer reactions. The nucleolytic ribozymes comprise a class of five structurally-distinct species that bring about site-specific cleavage by nucleophilic attack of the 2'-O on the adjacent 3'-P to form a cyclic 2',3'-phosphate. In general, they will also catalyse the reverse reaction. As a class, all these ribozymes appear to use general acid-base catalysis to accelerate these reactions by about a million-fold. In the Varkud satellite ribozyme, we have shown that the cleavage reaction is catalysed by guanine and adenine nucleobases acting as general base and acid, respectively. The hairpin ribozyme most probably uses a closely similar mechanism. Guanine nucleobases appear to be a common choice of general base, but the general acid is more variable. By contrast, the larger ribozymes such as the self-splicing introns and RNase P act as metalloenzymes.     

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.     

An alternative route for the folding of large RNAs: apparent two-state folding by a group II intron ribozyme

Despite a growing literature on the folding of RNA, our understanding of tertiary folding in large RNAs derives from studies on a small set of molecular examples, with primary focus on group I introns and RNase P RNA. To broaden the scope of RNA folding models and to better understand group II intron function, we have examined the tertiary folding of a ribozyme (D135) that is derived from the self-splicing ai5gamma intron from yeast mitochondria. The D135 ribozyme folds homogeneously and cooperatively into a compact, well-defined tertiary structure that includes all regions critical for active-site organization and substrate recognition. When D135 was treated with increasing concentrations of Mg(2+) and then subjected to hydroxyl radical footprinting, similar Mg(2+) dependencies were seen for internalization of all regions of the molecule, suggesting a highly cooperative folding behavior. In this work, we show that global folding and compaction of the molecule have the same magnesium dependence as the local folding previously observed. Furthermore, urea denaturation studies indicate highly cooperative unfolding of the ribozyme that is governed by thermodynamic parameters similar to those for forward folding. In fact, D135 folds homogeneously and cooperatively from the unfolded state to its native, active structure, thereby demonstrating functional reversibility in RNA folding. Taken together, the data are consistent with two-state folding of the D135 ribozyme, which is surprising given the size and multi-domain structure of the RNA. The findings establish that the accumulation of stable intermediates prior to formation of the native state is not a universal feature of RNA folding and that there is an alternative paradigm in which the folding landscape is relatively smooth, lacking rugged features that obstruct folding to the native state.     

Hydrolysis of organophosphate compounds by mutant butyrylcholinesterase: a story of two histidines

This study is aimed at understanding the hydrolysis mechanism of organophosphate (OP) compounds by G117H-BChE. It is a theoretical study that focuses on the role of the G117H mutation in the dephosphorylation step. Various proposed mechanisms are examined. We show that His117 acts as a general base by activating a water molecule, and thus assisting its nucleophilic attack on the phosphate. The calculated reaction energy profile agrees well with the experimental data. Moreover, analysis of the reaction via its two hypothetical elementary steps, proton transfer and hydroxide attack, supports the role of His117 as a general base. Further support to the proposed mechanism is gained by structural comparison of the active site to RNAse A, which has similar composition of substrate and functional groups. The similarity between these enzymes extends beyond the structure and also becomes evident when comparing functionality of various active sites residues as well as rate-pH dependence obtained in the two cases. Moreover, it is demonstrated that an extended form of Bevilacqua's model (Biochemistry 2003;42:2259-2265) may resolve the apparent contradictions between the proposed mechanism and various experimental observations regarding rate-pH dependence. Finally, that same model is shown to rationalize the hydrolase activity of G117D BChE, an observation which is considered puzzling. It is concluded that G117H-BChE hydrolyzes echothiophate and possibly other OP compounds via a general acid-base mechanism. On the basis of this mechanism, one can now proceed with rational design aimed at improving the enzyme by exploiting both the structural and mechanistic knowledge.     

High-resolution structures of the SAMHD1 dGTPase homolog from Leeuwenhoekiella blandensis reveal a novel mechanism of allosteric activation by dATP

Deoxynucleoside triphosphate (dNTP) triphosphohydrolases (dNTPases) are important enzymes that may perform multiple functions in the cell, including regulating the dNTP pools and contributing to innate immunity against viruses. Among the homologs that are best studied are human sterile alpha motif and HD domain–containing protein 1 (SAMHD1), a tetrameric dNTPase, and the hexameric Escherichia coli dGTPase; however, it is unclear whether these are representative of all dNTPases given their wide distribution throughout life. Here, we investigated a hexameric homolog from the marine bacterium Leeuwenhoekiella blandensis, revealing that it is a dGTPase that is subject to allosteric activation by dATP, specifically. Allosteric regulation mediated solely by dATP represents a novel regulatory feature among dNTPases that may facilitate maintenance of cellular dNTP pools in L. blandensis. We present high-resolution X-ray crystallographic structures (1.80–2.26 Å) in catalytically important conformations as well as cryo-EM structures (2.1–2.7 Å) of the enzyme bound to dGTP and dATP ligands. The structures, the highest resolution cryo-EM structures of any SAMHD1-like dNTPase to date, reveal an intact metal-binding site with the dGTP substrate coordinated to three metal ions. These structural and biochemical data yield insights into the catalytic mechanism and support a conserved catalytic mechanism for the tetrameric and hexameric dNTPase homologs. We conclude that the allosteric activation by dATP appears to rely on structural connectivity between the allosteric and active sites, as opposed to the changes in oligomeric state upon ligand binding used by SAMHD1.     

Mitogenome polymorphism in a single branch sample revealed by SOLiD deep sequencing of the Lophelia pertusa coral genome

We present an initial genomic analysis of the non-symbiotic scleractinian coral Lophelia pertusa, the dominant cold-water reef-building coral species in the North Atlantic Ocean. A significant fraction of the deep sequencing reads was of mitochondrial and microbial origins. SOLiD deep sequencing reads from fragment library experiments of total DNA and PCR amplified mitogenome generated about 21,000 times and 136,000 times coverage, respectively, of the 16,150bp mitogenome. Five polymorphic sites that include two non-synonymous sites in the NADH dehydrogenase subunit 5 genes were detected in both experiments. This observation is surprising since anthozoans in general exhibit very low mtDNA sequence variation at intraspecific level compared to nuclear sequences. More than fifty bacterial species associated with the coral isolate were also sequence detected, representing at least ten complete genomes. Most reads, however, were predicted to originate from the Lophelia nuclear genome.     

Noncanonical G(syn)-G(anti) base pairs stabilized by sulphate anions in two X-ray structures of the (GUGGUCUGAUGAGGCC) RNA duplex

The structures of two crystal forms of the RNA 16-mer with the sequence GUGGUCUGAUGAGGCC, grown in the presence of a high concentration of sulphate ions, have been determined using synchrotron radiation at 1.4- and 2.0-Å resolution. RNA with this sequence is known as one of the two strands of the noncleavable form of the hammerhead ribozyme. In both crystal structures, two G(syn)–G(anti) noncanonical base pairs are observed in the middle of a 14 base-pair (bp) duplex having 5′-dangling GU residues. Both structures contain sulphate anions interacting with the G–G bp stabilizing G in its syn conformation and bridging the two RNA strands. In both cases the interactions take place in the major groove, although the anions are accommodated within different helix geometries, most pronounced in the changing width of the major groove. In one structure, where a single sulphate spans both G–G pairs, the major groove is closed around the anion, while in the other structure, where each of the two G–G pairs is associated with a separate sulphate, the groove is open. This work provides the first examples of a G–G pair in syn-anti conformation, which minimizes the purine–purine clash in the center of the duplex, while utilizing its residual hydrogen bonding potential in specific interactions with sulphate anions.     

Critical roles of the CuB site in efficient proton pumping as revealed by crystal structures of mammalian cytochrome c oxidase catalytic intermediates

Mammalian cytochrome c oxidase (CcO) reduces O₂ to water in a bimetallic site including Fe_a3 and Cu_B giving intermediate molecules, termed A-, P-, F-, O-, E-, and R-forms. From the P-form on, each reaction step is driven by single-electron donations from cytochrome c coupled with the pumping of a single proton through the H-pathway, a proton-conducting pathway composed of a hydrogen-bond network and a water channel. The proton-gradient formed is utilized for ATP production by F-ATPase. For elucidation of the proton pumping mechanism, crystal structural determination of these intermediate forms is necessary. Here we report X-ray crystallographic analysis at ∼1.8 Å resolution of fully reduced CcO crystals treated with O₂ for three different time periods. Our disentanglement of intermediate forms from crystals that were composed of multiple forms determined that these three crystallographic data sets contained ∼45% of the O-form structure, ∼45% of the E-form structure, and ∼20% of an oxymyoglobin-type structure consistent with the A-form, respectively. The O- and E-forms exhibit an unusually long Cu_B²⁺-OH⁻ distance and Cu_B¹⁺-H₂O structure keeping Fe_a3³⁺-OH⁻ state, respectively, suggesting that the O- and E-forms have high electron affinities that cause the O→E and E→R transitions to be essentially irreversible and thus enable tightly coupled proton pumping. The water channel of the H-pathway is closed in the O- and E-forms and partially open in the R-form. These structures, together with those of the recently reported P- and F-forms, indicate that closure of the H-pathway water channel avoids back-leaking of protons for facilitating the effective proton pumping.     

Accumulation of piRNAs in the chromatoid bodies purified by a novel isolation protocol

Haploid male germ cells are featured by an intriguing cytoplasmic cloud-like structure that has been named as chromatoid body (CB) on the basis of its staining properties and appearance under a microscope. Notwithstanding its early discovery in the late 19th century, the function of the CB is still largely obscure. Emerging evidence suggests a role for the CB and other similar RNA-containing granules, such as germ plasm in lower organism and processing bodies in somatic cells, in the control and organization of RNA processing and/or storage. Despite the increasing scientific demand, the lack of CB purification protocols has still been the main obstacle in the functional characterization of this structure. We have successfully isolated CBs from mouse testis by a novel immunoaffinity purification procedure and validated by several different methods that pure CB fractions are obtained. Analysis of the CB RNA content reveals enrichment of PIWI-interacting RNAs (piRNAs), further emphasizing the role of CB as the RNA processing body.     

A novel mutation in the myeloperoxidase gene in a Chinese female with complete myeloperoxidase deficiency: The role of nonsense-mediated mRNA decay

Myeloperoxidase (MPO) is an important enzyme in innate immunity. Here, we describe the first identified Chinese individual with complete MPO deficiency. The proband was ascertained through routine automated complete blood analysis. Analysis of MPO function and immunogenicity revealed that MPO levels in neutrophils were significantly decreased. Mutational analysis revealed a novel premature termination codon p.(Trp602*) in exon 11 of the MPO gene, which was inherited in an autosomal recessive manner. We demonstrated that nonsense-mediated mRNA decay is involved in the molecular pathology of MPO deficiency in this case. The study of MPO deficiency can be helpful in understanding the function and biosynthesis mechanisms of MPO.     

Catalytic role of a conserved cysteine residue in the desulfonation reaction by the alkanesulfonate monooxygenase enzyme

Detailed kinetic studies were performed in order to determine the role of the single cysteine residue in the desulfonation reaction catalyzed by SsuD. Mutation of the conserved cysteine at position 54 in SsuD to either serine or alanine had little effect on FMNH₂ binding. The k_cat/K_m value for the C54S SsuD variant increased 3-fold, whereas the k_cat/K_m value for C54A SsuD decreased 6-fold relative to wild-type SsuD. An initial fast phase was observed in kinetic traces obtained for the oxidation of flavin at 370 nm when FMNH₂ was mixed against C54S SsuD (k_obs, 111 s^− 1) in oxygenated buffer that was 10-fold faster than wild-type SsuD (k_obs, 12.9 s^− 1). However, there was no initial fast phase observed in similar kinetic traces obtained for C54A SsuD. This initial fast phase was previously assigned to the formation of the C4a-(hydro)peroxyflavin in studies with wild-type SsuD. There was no evidence for the formation of the C4a-(hydro)peroxyflavin with either SsuD variant when octanesulfonate was included in rapid reaction kinetic studies, even at low octanesulfonate concentrations. The absence of any C4a-(hydro)peroxyflavin accumulation correlates with the increased catalytic activity of C54S SsuD. These results suggest that the conservative serine substitution is able to effectively take the place of cysteine in catalysis. Conversely, decreased accumulation of the C4a-(hydro)peroxyflavin intermediate with the C54A SsuD variant may be due to decreased activity. The data described suggest that Cys54 in SsuD may be either directly or indirectly involved in stabilizing the C4a-(hydro)peroxyflavin intermediate formed during catalysis through hydrogen bonding interactions.     

An RNA electrophoretic mobility shift and mutational analysis of rnp-4f 5′-UTR intron splicing regulatory proteins in Drosophila reveals a novel new role for a dADAR protein isoform

Alternative splicing greatly enhances the diversity of proteins encoded by eukaryotic genomes, and is also important in gene expression control. In contrast to the great depth of knowledge as to molecular mechanisms in the splicing pathway itself, relatively little is known about the regulatory events behind this process. The 5′-UTR and 3′-UTR in pre-mRNAs play a variety of roles in controlling eukaryotic gene expression, including translational modulation, and nearly 4000 of the roughly 14,000 protein coding genes in Drosophila contain introns of unknown functional significance in their 5′-UTR. Here we report the results of an RNA electrophoretic mobility shift analysis of Drosophila rnp-4f 5′-UTR intron 0 splicing regulatory proteins. The pre-mRNA potential regulatory element consists of an evolutionarily-conserved 177-nt stem-loop arising from pairing of intron 0 with part of adjacent exon 2. Incubation of in vitro transcribed probe with embryo protein extract is shown to result in two shifted RNA–protein bands, and protein extract from a dADAR null mutant fly line results in only one shifted band. A mutated stem-loop in which the conserved exon 2 primary sequence is changed but secondary structure maintained by introducing compensatory base changes results in diminished band shifts. To test the hypothesis that dADAR plays a role in intron splicing regulation in vivo, levels of unspliced rnp-4f mRNA in dADAR mutant were compared to wild-type via real-time qRT-PCR. The results show that during embryogenesis unspliced rnp-4f mRNA levels fall by up to 85% in the mutant, in support of the hypothesis. Taken together, these results demonstrate a novel role for dADAR protein in rnp-4f 5′-UTR alternative intron splicing regulation which is consistent with a previously proposed model.     

Characterization of a novel growth hormone receptor-encoding cDNA in rainbow trout and regulation of its expression by nutritional state

To clarify the divergence of the growth hormone receptor (GHR) family, we characterized a novel GHR from a teleost fish (rainbow trout). A 2357-nt cDNA was isolated and found to contain a single initiation site 71 nt from the most 5′ end, an open reading frame of 1971 nt encoding a 657-amino acid protein, and a single polyadenylation site 229 nt from the poly-A tail. Based on structural analysis, the protein was identified as a type 1 GHR (GHR1). The new GHR1 shares 42% and 43% amino acid identity, respectively, with GHR2a and GHR2b, the two type 2 GHRs isolated from trout previously. GHR1 mRNA was found in a wide array of tissues with the highest expression in the liver, red muscle, and white muscle. Fasting animals for 4 weeks reduced steady state levels of GHR1 in the liver, adipose, and red muscle. These findings help clarify the divergence and nomenclature of GHRs and provide insight into the function of duplicated GHR types.     

Crystal structure of the streptococcal superantigen SpeI and functional role of a novel loop domain in T cell activation by group V superantigens

Superantigens (SAgs) are potent microbial toxins that bind simultaneously to T cell receptors (TCRs) and class II major histocompatibility complex molecules, resulting in the activation and expansion of large T cell subsets and the onset of numerous human diseases. Within the bacterial SAg family, streptococcal pyrogenic exotoxin I (SpeI) has been classified as belonging to the group V SAg subclass, which are characterized by a unique, relatively conserved ∼15 amino acid extension (amino acid residues 154 to 170 in SpeI; herein referred to as the α3-β8 loop), absent in SAg groups I through IV. Here, we report the crystal structure of SpeI at 1.56 Å resolution. Although the α3-β8 loop in SpeI is several residues shorter than that of another group V SAg, staphylococcal enterotoxin serotype I, the C-terminal portions of these loops, which are located adjacent to the putative TCR binding site, are structurally similar. Mutagenesis and subsequent functional analysis of SpeI indicates that TCR β-chains are likely engaged in a similar general orientation as other characterized SAgs. We show, however, that the α3-β8 loop length, and the presence of key glycine residues, are necessary for optimal activation of T cells. Based on Vβ-skewing analysis of human T cells activated with SpeI and structural models, we propose that the α3-β8 loop is positioned to form productive intermolecular contacts with the TCR β-chain, likely in framework region 3, and that these contacts are required for optimal TCR recognition by SpeI, and likely all other group V SAgs.     

Development of a novel chemiluminescence method for the determination of cefazolin sodium in injectable powder and human urine based on a luminol‐Cu(III) complex reaction in alkaline medium

A novel chemiluminescence (CL) method was developed for the determination of cefazolin sodium based on the CL reaction between the [Cu(HIO₆)₂]^5‐Cu(III) complex and luminol in alkaline solution. Results showed that CL emission of Cu(III) complex–luminol in alkaline medium was significantly different from that in acidic medium. A possible mechanism of the enhanced effect of cefazolin on CL emission of the [Cu(HIO₆)₂]^5‐‐ luminol system was proposed. The effect of the reaction conditions on CL emissions was examined. Under optimized conditions, a good linear relationship was obtained between CL intensity and concentrations of cefazolin sodium in the range of 2.0 x 10^‐8 to 2.0 x 10^‐6 g/mL with a correlation coefficient of R² = 0.9978. The limit of detection was 4.58 x 10^‐9 g/mL. The proposed method was applied for the determination of cefazolin sodium in real samples with recoveries of 82.0‐109% with an RSD of 0.7‐2.1%. The proposed method was successfully used for the determination of cefazolin sodium in injectable powder preparations and human urine with satisfactory results. Copyright © 2012 John Wiley & Sons, Ltd.