Crystal structure of RNase T1 with 3'-guanylic acid and guanosine

RNA Cleavage Properties of Nucleobase-Specific RNase MC1 and Cusativin Are Determined by the Dinucleotide-Binding Interactions in the Enzyme-Active Site

Knowledge of the cleavage specificity of ribonucleases is critical for their application in RNA modification mapping or RNA-protein binding studies. Here, we detail the cleavage specificity and efficiency of ribonuclease MC1 and cusativin using a customized RNA sequence that contained all dinucleotide combinations and homopolymer sequences. The sequencing of the oligonucleotide digestion products by a semi-quantitative liquid chromatography coupled with mass spectrometry (LC-MS) analysis documented as little as 0.5-1% cleavage levels for a given dinucleotide sequence combination. While RNase MC1 efficiently cleaved the [A/U/C]pU dinucleotide bond, no cleavage was observed for the GpU bond. Similarly, cusativin efficiently cleaved Cp[U/A/G] dinucleotide combinations along with UpA and [A/U]pU, suggesting a broader specificity of dinucleotide preferences. The molecular interactions between the substrate and active site as determined by the dinucleotide docking studies of protein models offered additional evidence and support for the observed substrate specificity. Targeted alteration of the key amino acid residues in the nucleotide-binding site confirms the utility of this in silico approach for the identification of key interactions. Taken together, the use of bioanalytical and computational approaches, involving LC-MS and ligand docking of tertiary structural models, can form a powerful combination to help explain the RNA cleavage behavior of RNases.

Reactivity and Specificity of RNase T1, RNase A, and RNase H toward Oligonucleotides of RNA Containing 8-Oxo-7,8-dihydroguanosine

Understanding how oxidatively damaged RNA interacts with ribonucleases is important because of its proposed role in the development and progression of disease. Thus, understanding structural aspects of RNA containing lesions generated under oxidative stress, as well as its interactions with other biopolymers, is fundamental. We explored the reactivity of RNase A, RNase T₁, and RNase H toward oligonucleotides of RNA containing 8-oxo-7,8-dihydroguanosine (8oxoG). This is the first example that addresses this relationship and will be useful for understanding (1) how these RNases can be used to characterize the structural impact that this lesion has on RNA and (2) how oxidatively modified RNA may be handled intracellularly. 8-OxoG was incorporated into 10-16-mers of RNA, and its reactivity with each ribonuclease was assessed via electrophoretic analyses, circular dichroism, and the use of other C8-purine-modified analogues (8-bromoguanosine, 8-methoxyguanosine, and 8-oxoadenosine). RNase T₁ does not recognize sites containing 8-oxoG, while RNase A recognizes and cleaves RNA at positions containing this lesion while differentiating if it is involved in H-bonding. The selectivity of RNase A followed the order C > 8-oxoG ≈ U. In addition, isothermal titration calorimetry showed that an 8-oxoG-C3'-methylphosphate derivative can inhibit RNase A activity. Cleavage patterns obtained from RNase H displayed changes in reactivity in a sequence- and concentration-dependent manner and displayed recognition at sites containing the modification in some cases. These data will aid in understanding how this modification affects reactivity with ribonucleases and will enable the characterization of global and local structural changes in oxidatively damaged RNA.

Mechanism of endonuclease cleavage by the HigB toxin

Bacteria encode multiple type II toxin-antitoxin modules that cleave ribosome-bound mRNAs in response to stress. All ribosome-dependent toxin family members structurally characterized to date adopt similar microbial RNase architectures despite possessing low sequence identities. Therefore, determining which residues are catalytically important in this specialized RNase family has been a challenge in the field. Structural studies of RelE and YoeB toxins bound to the ribosome provided significant insights but biochemical experiments with RelE were required to clearly demonstrate which residues are critical for acid-base catalysis of mRNA cleavage. Here, we solved an X-ray crystal structure of the wild-type, ribosome-dependent toxin HigB bound to the ribosome revealing potential catalytic residues proximal to the mRNA substrate. Using cell-based and biochemical assays, we further determined that HigB residues His54, Asp90, Tyr91 and His92 are critical for activity in vivo, while HigB H54A and Y91A variants have the largest effect on mRNA cleavage in vitro Comparison of X-ray crystal structures of two catalytically inactive HigB variants with 70S-HigB bound structures reveal that HigB active site residues undergo conformational rearrangements likely required for recognition of its mRNA substrate. These data support the emerging concept that ribosome-dependent toxins have diverse modes of mRNA recognition.

Transition State Charge Stabilization and Acid-Base Catalysis of mRNA Cleavage by the Endoribonuclease RelE

The bacterial toxin RelE is a ribosome-dependent endoribonuclease. It is part of a type II toxin-antitoxin system that contributes to antibiotic resistance and biofilm formation. During amino acid starvation, RelE cleaves mRNA in the ribosomal A-site, globally inhibiting protein translation. RelE is structurally similar to microbial RNases that employ general acid-base catalysis to facilitate RNA cleavage. The RelE active site is atypical for acid-base catalysis, in that it is enriched with positively charged residues and lacks the prototypical histidine-glutamate catalytic pair, making the mechanism of mRNA cleavage unclear. In this study, we use a single-turnover kinetic analysis to measure the effect of pH and phosphorothioate substitution on the rate constant for cleavage of mRNA by wild-type RelE and seven active-site mutants. Mutation and thio effects indicate a major role for stabilization of increased negative change in the transition state by arginine 61. The wild-type RelE cleavage rate constant is pH-independent, but the reaction catalyzed by many of the mutants is strongly dependent on pH, suggestive of general acid-base catalysis. pH-rate curves indicate that wild-type RelE operates with the pK(a) of at least one catalytic residue significantly downshifted by the local environment. Mutation of any single active-site residue is sufficient to disrupt this microenvironment and revert the shifted pK(a) back above neutrality. pH-rate curves are consistent with K54 functioning as a general base and R81 as a general acid. The capacity of RelE to effect a large pK(a) shift and facilitate a common catalytic mechanism by uncommon means furthers our understanding of other atypical enzymatic active sites.

Defining the mRNA recognition signature of a bacterial toxin protein

Bacteria contain multiple type II toxins that selectively degrade mRNAs bound to the ribosome to regulate translation and growth and facilitate survival during the stringent response. Ribosome-dependent toxins recognize a variety of three-nucleotide codons within the aminoacyl (A) site, but how these endonucleases achieve substrate specificity remains poorly understood. Here, we identify the critical features for how the host inhibition of growth B (HigB) toxin recognizes each of the three A-site nucleotides for cleavage. X-ray crystal structures of HigB bound to two different codons on the ribosome illustrate how HigB uses a microbial RNase-like nucleotide recognition loop to recognize either cytosine or adenosine at the second A-site position. Strikingly, a single HigB residue and 16S rRNA residue C1054 form an adenosine-specific pocket at the third A-site nucleotide, in contrast to how tRNAs decode mRNA. Our results demonstrate that the most important determinant for mRNA cleavage by ribosome-dependent toxins is interaction with the third A-site nucleotide.

The Lepidopteran endoribonuclease-U domain protein P102 displays dramatically reduced enzymatic activity and forms functional amyloids

Hemocytes of Heliothis virescens (F.) (Lepidoptera, Noctuidae) larvae produce a protein, P102, with a putative endoribonuclease-U domain. In previous works we have shown that P102 is involved in Lepidopteran immune response by forming amyloid fibrils, which catalyze and localize melanin deposition around non-self intruders during encapsulation, preventing harmful systemic spreading. Here we demonstrate that P102 belongs to a new class of proteins that, at least in Lepidoptera, has a diminished endoribonuclease-U activity probably due to the lack of two out of five catalytically essential residues. We show that the P102 homolog from Trichoplusia ni (Lepidoptera, Noctuidae) displays catalytic site residues identical to P102, a residual endoribonuclease-U activity and the ability to form functional amyloids. On the basis of these results as well as sequence and structural analyses, we hypothesize that all the Lepidoptera endoribonuclease-U orthologs with catalytic site residues identical to P102 form a subfamily with similar function.

Mechanisms of toxin inhibition and transcriptional repression by Escherichia coli DinJ-YafQ

Bacteria encounter environmental stresses that regulate a gene expression program required for adaptation and survival. Here, we report the 1.8-Å crystal structure of the Escherichia coli toxin-antitoxin complex YafQ-(DinJ)2-YafQ, a key component of the stress response. The antitoxin DinJ dimer adopts a ribbon-helix-helix motif required for transcriptional autorepression, and toxin YafQ contains a microbial RNase fold whose proposed active site is concealed by DinJ binding. Contrary to previous reports, our studies indicate that equivalent levels of transcriptional repression occur by direct interaction of either YafQ-(DinJ)2-YafQ or a DinJ dimer at a single inverted repeat of its recognition sequence that overlaps with the -10 promoter region. Surprisingly, multiple YafQ-(DinJ)2-YafQ complexes binding to the operator region do not appear to amplify the extent of repression. Our results suggest an alternative model for transcriptional autorepression that may be novel to DinJ-YafQ.

Bacterial toxin RelE: a highly efficient ribonuclease with exquisite substrate specificity using atypical catalytic residues

The toxin RelE is a ribosome-dependent endoribonuclease implicated in diverse cellular processes, including persistence. During amino acid starvation, RelE inhibits translation by cleaving ribosomal A-site mRNA. Although RelE is structurally similar to other microbial endoribonucleases, the active-site amino acid composition differs substantially and lacks obvious candidates for general acid-base functionality. Highly conserved RelE residues (Lys52, Lys54, Arg61, Arg81, and Tyr87) surround the mRNA scissile phosphate, and specific 16S rRNA contacts further contribute to substrate positioning. We used a single-turnover kinetic assay to evaluate the catalytic importance of individual residues in the RelE active site. Within the context of the ribosome, RelE rapidly cleaves A-site mRNA at a rate similar to those of traditional ribonucleases. Single-turnover rate constants decreased between 10(2)- and 10(6)-fold for the RelE active-site mutants of Lys52, Lys54, Arg61, and Arg81. RelE may principally promote catalysis via transition-state charge stabilization and leaving-group protonation, in addition to achieving in-line substrate positioning in cooperation with the ribosome. This kinetic analysis complements structural information to provide a foundation for understanding the molecular mechanism of this atypical endoribonuclease.

The structural basis for mRNA recognition and cleavage by the ribosome-dependent endonuclease RelE

Translational control is widely used to adjust gene expression levels. During the stringent response in bacteria, mRNA is degraded on the ribosome by the ribosome-dependent endonuclease, RelE. The molecular basis for recognition of the ribosome and mRNA by RelE and the mechanism of cleavage are unknown. Here, we present crystal structures of E. coli RelE in isolation (2.5 Å) and bound to programmed Thermus thermophilus 70S ribosomes before (3.3 Å) and after (3.6 Å) cleavage. RelE occupies the A site and causes cleavage of mRNA after the second nucleotide of the codon by reorienting and activating the mRNA for 2′-OH-induced hydrolysis. Stacking of A site codon bases with conserved residues in RelE and 16S rRNA explains the requirement for the ribosome in catalysis and the subtle sequence specificity of the reaction. These structures provide detailed insight into the translational regulation on the bacterial ribosome by mRNA cleavage.

Characterization of a heat-stable enzyme possessing GTP-dependent RNA ligase activity from a hyperthermophilic archaeon, Pyrococcus furiosus

Using an expression protein library of a hyperthermophilic archaeon, Pyrococcus furiosus, we identified a gene (PF0027) that encodes a protein with heat-stable cyclic nucleotide phosphodiesterase (CPDase) activity. The PF0027 gene encoded a 21-kDa protein and an amino acid sequence that showed approximately 27% identity to that of the 2'-5' tRNA ligase protein, ligT (20 kDa), from Escherichia coli. We found that the purified PF0027 protein possessed GTP-dependent RNA ligase activity and that synthetic tRNA halves bearing 2',3'-cyclic phosphate and 5'-OH termini were substrates for the ligation reaction in vitro. GTP hydrolysis was not required for the reaction, and GTPgammaS enhanced the tRNA ligation activity of PF0027 protein, suggesting that the ligation step is regulated by a novel mechanism. In comparison to the strong CPDase activity of the PF0027 protein, the RNA ligase activity itself was quite weak, and the ligation product was unstable during in vitro reaction. Finally, we used NMR to determine the solution structure of the PF0027 protein and discuss the implications of our results in understanding the role of the PF0027 protein.

AMMOS: Automated Molecular Mechanics Optimization tool for in silico Screening

Background

Virtual or in silico ligand screening combined with other computational methods is one of the most promising methods to search for new lead compounds, thereby greatly assisting the drug discovery process. Despite considerable progresses made in virtual screening methodologies, available computer programs do not easily address problems such as: structural optimization of compounds in a screening library, receptor flexibility/induced-fit, and accurate prediction of protein-ligand interactions. It has been shown that structural optimization of chemical compounds and that post-docking optimization in multi-step structure-based virtual screening approaches help to further improve the overall efficiency of the methods. To address some of these points, we developed the program AMMOS for refining both, the 3D structures of the small molecules present in chemical libraries and the predicted receptor-ligand complexes through allowing partial to full atom flexibility through molecular mechanics optimization.

Results

The program AMMOS carries out an automatic procedure that allows for the structural refinement of compound collections and energy minimization of protein-ligand complexes using the open source program AMMP. The performance of our package was evaluated by comparing the structures of small chemical entities minimized by AMMOS with those minimized with the Tripos and MMFF94s force fields. Next, AMMOS was used for full flexible minimization of protein-ligands complexes obtained from a mutli-step virtual screening. Enrichment studies of the selected pre-docked complexes containing 60% of the initially added inhibitors were carried out with or without final AMMOS minimization on two protein targets having different binding pocket properties. AMMOS was able to improve the enrichment after the pre-docking stage with 40 to 60% of the initially added active compounds found in the top 3% to 5% of the entire compound collection.

Conclusion

The open source AMMOS program can be helpful in a broad range of in silico drug design studies such as optimization of small molecules or energy minimization of pre-docked protein-ligand complexes. Our enrichment study suggests that AMMOS, designed to minimize a large number of ligands pre-docked in a protein target, can successfully be applied in a final post-processing step and that it can take into account some receptor flexibility within the binding site area.

Analysis of internal motions of RNase T1 complexed with a productive substrate involving 15N NMR relaxation measurements

The backbone dynamics of RNase T1 in the presence of exo-guanosine 2',3'-cyclophosphorothioate (exo-cGPS isomer), which is a productive substrate, and in the presence of 3'-guanylic acid (3'GMP), which is an nonproductive substrate, were examined using (15)N nuclear magnetic resonance. Although the X-ray crystal structure suggests that the modes of binding of these substrates to the active-site cleft are very similar, the order parameters in a number of regions in RNase T1 complexed with exo-cGPS isomer were different from those with 3'GMP. Moreover, the chemical exchange in line width observed for RNase T1 complexed with exo-cGPS isomer was also different from that observed for RNase T1 complexed with 3'GMP. From these results, we concluded that the internal motions in RNase T1 complexed with a productive substrate were not always identical to those in RNase T1 complexed with a nonproductive substrate.

The structure of the endoribonuclease XendoU: From small nucleolar RNA processing to severe acute respiratory syndrome coronavirus replication

Small nucleolar RNAs (snoRNAs) play a key role in eukaryotic ribosome biogenesis. In most cases, snoRNAs are encoded in introns and are released through the splicing reaction. Some snoRNAs are, instead, produced by an alternative pathway consisting of endonucleolytic processing of pre-mRNA. XendoU, the endoribonuclease responsible for this activity, is a U-specific, metal-dependent enzyme that releases products with 2'-3' cyclic phosphate termini. XendoU is broadly conserved among eukaryotes, and it is a genetic marker of nidoviruses, including the severe acute respiratory syndrome coronavirus, where it is essential for replication and transcription. We have determined by crystallography the structure of XendoU that, by refined search methodologies, appears to display a unique fold. Based on sequence conservation, mutagenesis, and docking simulations, we have identified the active site. The conserved structural determinants of this site may provide a framework for attempting to design antiviral drugs to interfere with the infectious nidovirus life cycle.

A nucleophile activation dyad in ribonucleases. A combined X-ray crystallographic/ab initio quantum chemical study

Ribonucleases (RNases) catalyze the cleavage of the phosphodiester bond in RNA up to 10(15)-fold, as compared with the uncatalyzed reaction. High resolution crystal structures of these enzymes in complex with 3'-mononucleotide substrates demonstrate the accommodation of the nucleophilic 2'-OH group in a binding pocket comprising the catalytic base (glutamate or histidine) and a charged hydrogen bond donor (lysine or histidine). Ab initio quantum chemical calculations performed on such Michaelis complexes of the mammalian RNase A (EC ) and the microbial RNase T(1) (EC ) show negative charge build up on the 2'-oxygen upon substrate binding. The increased nucleophilicity results from stronger hydrogen bonding to the catalytic base, which is mediated by a hydrogen bond from the charged donor. This hitherto unrecognized catalytic dyad in ribonucleases constitutes a general mechanism for nucleophile activation in both enzymic and RNA-catalyzed phosphoryl transfer reactions.

The contribution of metal ions to the conformational stability of ribonuclease T1: crystal versus solution

In the crystalline state, ribonuclease T1 binds calcium ions at different lattice-dependent positions. In solution, its conformational stability is also remarkably increased in the presence of divalent metal ions. Combining urea unfolding studies and X-ray crystallography, we compared the presence of several metal ions at specific sites in the protein to their contribution to the overall stabilizing effect in solution. We constructed thermodynamic cycles involving particular metal ions and specific carboxylate functions. The resulting coupling energies indicate that some (but not all) metal ions found at lattice contacts in crystal structures may indeed significantly contribute to stability enhancement in the presence of metal ions in solution.

Analysis of a water mediated protein-protein interactions within RNase T1

Buried and well-ordered solvent molecules are an integral part of each folded protein. For a few individual water molecules, the exchange kinetics with solvent have been described in great detail. So far, little is known about the energetics of this exchange process. Here, we present an experimental approach to investigate water-mediated intramolecular protein-protein interactions by use of double mutant cycles. As a first example, we analyzed the interdependence of the contribution of two side chains (Asn9 and Thr93) to the conformational stability of RNase T1. In the folded state, both side chains are involved in the "solvation of the same water molecule WAT1. The coupling of the contributions of Asn9 and Thr93 to the conformational stability of RNase T1 was measured by urea unfolding and differential scanning calorimetry. The structural integrity of each mutant was analyzed by X-ray crystallography. We find that the effects of the Asn9Ala and the Thr93Ala mutations on the conformational stability are additive in the corresponding double mutant. We conclude that the free energy of the WAT1 mediated intramolecular protein-protein interaction in the folded state is very similar to solvent mediated protein-protein interaction in the unfolded state.

Hydrogen-exchange stabilities of RNase T1 and variants with buried and solvent-exposed Ala --> Gly mutations in the helix

Hydrogen-exchange rates were measured for RNase T1 and three variants with Ala --> Gly substitutions at a solvent-exposed (residue 21) and a buried (residue 23) position in the helix: A21G, G23A, and A21G + G23A. These results were used to measure the stabilities of the proteins. The hydrogen-exchange stabilities (DeltaG(HX)) for the most stable residues in each variant agree with the equilibrium conformational stability measured by urea denaturation (DeltaG(U)), if the effects of D(2)O and proline isomerization are included [Huyghues-Despointes, B. M. P., Scholtz, J. M., and Pace, C. N. (1999) Nat. Struct. Biol. 6, 210-212]. These residues also show similar changes in DeltaG(HX) upon Ala --> Gly mutations (DeltaDeltaG(HX)) as compared to equilibrium measurements (DeltaDeltaG(U)), indicating that the most stable residues are exchanging from the globally unfolded ensemble. Alanine is stabilizing compared to glycine by 1 kcal/mol at a solvent-exposed site 21 as seen by other methods for the RNase T1 protein and peptide helix [Myers, J. K., Pace, C. N., and Scholtz, J. M. (1997) Proc. Natl. Acad. Sci. U.S.A. 94, 3833-2837], while it is destabilizing at the buried site 23 by the same amount. For the A21G variant, only local NMR chemical shift perturbations are observed compared to RNase T1. For the G23A variant, large chemical shift changes are seen throughout the sequence, although X-ray crystal structures of the variant and RNase T1 are nearly superimposable. Ala --> Gly mutations in the helix of RNase T1 at both helical positions alter the native-state hydrogen-exchange stabilities of residues throughout the sequence.

Ribonuclease activity of rat liver perchloric acid-soluble protein, a potent inhibitor of protein synthesis

Rat liver perchloric acid-soluble protein (L-PSP) is a potent inhibitor of cell-free protein synthesis; however, its mechanism of action is not known. Here we show that the protein is a unique ribonuclease and that this activity is responsible for the inhibition of translation. The addition of perchloric acid-soluble protein to a rabbit reticulocyte cell-free system at a concentration of 6.2 microM led to an almost complete inhibition of protein synthesis. The kinetics are unlike those of hemin-controlled inhibitor, a protein that acts at the initiation step. The inhibition appears to be due to an endoribonucleolytic activity of perchloric acid-soluble protein because L-PSP directly affects mRNA template activity and induces disaggregation of the reticulocyte polysomes into 80 S ribosomes, even in the presence of cycloheximide. These effects were observed with authentic as well as recombinant L-PSP. Analysis by thin-layer chromatography of [alpha-32P]UTP-labeled mRNA incubated with the protein showed production of the ribonucleoside 3'-monophosphates Ap, Gp, Up, and Cp, providing direct evidence that the protein is an endoribonuclease. When either 5'- or 3'-32P-labeled 5 S rRNA was the substrate, L-PSP cleaved phosphodiester bonds only in the single-stranded regions of the molecule.

Conserved water molecules in a large family of microbial ribonucleases

We systematically analyzed the crystallographically determined water molecules of all known structures of RNase T1 and compared them to the ordered solvent in a large number of related microbial nucleases. To assess the crystallographers' impact on the interpretation of the solvent structure, we independently refined five validation structures from diffraction data derived from five isomorphous crystals of RNase T1. We also compared the positions of water molecules found in 11 published isomorphous RNase T1 inhibitor complexes. These data suggest that the positions of most of the waters located on the surface of a protein and that are well-determined in the experimental electron density maps are determined primarily by crystal packing forces. Water molecules with less well-defined electron density are in general unique to one or a small number of crystal structures. Only a small number of the well-defined waters are found to be independent of the crystal environment. These waters have a low accessible surface area and B-factor, and tend to be conserved in the crystal structures of a number of evolutionary related ribonucleases as well. A single water molecule is found conserved in all known microbial ribonucleases.

Dissection of the structural and functional role of a conserved hydration site in RNase T1

The reoccurrence of water molecules in crystal structures of RNase T1 was investigated. Five waters were found to be invariant in RNase T1 as well as in six other related fungal RNases. The structural, dynamical, and functional characteristics of one of these conserved hydration sites (WAT1) were analyzed by protein engineering, X-ray crystallography, and (17)O and 2H nuclear magnetic relaxation dispersion (NMRD). The position of WAT1 and its surrounding hydrogen bond network are unaffected by deletions of two neighboring side chains. In the mutant Thr93Gln, the Gln93N epsilon2 nitrogen replaces WAT1 and participates in a similar hydrogen bond network involving Cys6, Asn9, Asp76, and Thr91. The ability of WAT1 to form four hydrogen bonds may explain why evolution has preserved a water molecule, rather than a side-chain atom, at the center of this intricate hydrogen bond network. Comparison of the (17)O NMRD profiles from wild-type and Thr93Gln RNase T1 yield a mean residence time of 7 ns at 27 degrees C and an orientational order parameter of 0.45. The effects of mutations around WAT1 on the kinetic parameters of RNase T1 are small but significant and probably relate to the dynamics of the active site.

Three-dimensional structure of ribonuclease T1 complexed with an isosteric phosphonate substrate analogue of GpU: alternate substrate binding modes and catalysis

The X-ray crystal structure of a complex between ribonuclease T1 and guanylyl(3'-6')-6'-deoxyhomouridine (GpcU) has been determined at 2. 0 A resolution. This ligand is an isosteric analogue of the minimal RNA substrate, guanylyl(3'-5')uridine (GpU), where a methylene is substituted for the uridine 5'-oxygen atom. Two protein molecules are part of the asymmetric unit and both have a GpcU bound at the active site in the same manner. The protein-protein interface reveals an extended aromatic stack involving both guanines and three enzyme phenolic groups. A third GpcU has its guanine moiety stacked on His92 at the active site on enzyme molecule A and interacts with GpcU on molecule B in a neighboring unit via hydrogen bonding between uridine ribose 2'- and 3'-OH groups. None of the uridine moieties of the three GpcU molecules in the asymmetric unit interacts directly with the protein. GpcU-active-site interactions involve extensive hydrogen bonding of the guanine moiety at the primary recognition site and of the guanosine 2'-hydroxyl group with His40 and Glu58. On the other hand, the phosphonate group is weakly bound only by a single hydrogen bond with Tyr38, unlike ligand phosphate groups of other substrate analogues and 3'-GMP, which hydrogen-bonded with three additional active-site residues. Hydrogen bonding of the guanylyl 2'-OH group and the phosphonate moiety is essentially the same as that recently observed for a novel structure of a RNase T1-3'-GMP complex obtained immediately after in situ hydrolysis of exo-(Sp)-guanosine 2',3'-cyclophosphorothioate [Zegers et al. (1998) Nature Struct. Biol. 5, 280-283]. It is likely that GpcU at the active site represents a nonproductive binding mode for GpU [Steyaert, J., and Engleborghs (1995) Eur. J. Biochem. 233, 140-144]. The results suggest that the active site of ribonuclease T1 is adapted for optimal tight binding of both the guanylyl 2'-OH and phosphate groups (of GpU) only in the transition state for catalytic transesterification, which is stabilized by adjacent binding of the leaving nucleoside (U) group.

Trimeric domain-swapped barnase

The structure of a trimeric domain-swapped form of barnase (EC 3.1. 27.3) was determined by x-ray crystallography at a resolution of 2.2 A from crystals of space group R32. Residues 1-36 of one molecule associate with residues 41-110 from another molecule related through threefold symmetry. The resulting cyclic trimer contains three protein folds that are very similar to those in monomeric barnase. Both swapped domains contain a nucleation site for folding. The formation of a domain-swapped trimer is consistent with the description of the folding process of monomeric barnase as the formation and subsequent association of two foldons.

Hydrolysis of a slow cyclic thiophosphate substrate of RNase T1 analyzed by time-resolved crystallography

Here we present a time-resolved crystallographic analysis of the hydrolysis of exo (Sp) guanosine 2',3'-cyclophosphorothioate by RNase T1. The use of a slow substrate and fast crystallization methods made it possible to perform the study with conventional data-collection techniques. The results support the idea that the hydrolysis reaction proceeds through a mechanism that is the inverse of the transesterification reaction. In addition, the structures provide an explanation for the differential behavior of RNase T1 towards exo- and endo-cyclic thiophosphates.

A decade of protein engineering on ribonuclease T1--atomic dissection of the enzyme-substrate interactions

During the last decade, protein engineering has been used to identify the residues that contribute to the ribonuclease-T1-catalyzed transesterification. His40, Glu58 and His92 accelerate the associative nucleophilic displacement at the phosphate atom by the entering 2'-oxygen downstream guanosines in a highly cooperative manner. Glu58, assisted by the protonated His40 imidazole, abstracts a proton from the 2'-oxygen, while His92 protonates the leaving group. Tyr38, Arg77 and Phe100 further stabilize the transition state of the reaction. A functionally independent subsite, including Asn36 and Asn98, contributes to chemical turnover by aligning the substrate relative to the catalytic side chains upon binding of the leaving group. An invariant structural motive, involving residues 42-46, renders ribonuclease T1 guanine specific through a series of intermolar hydrogen bonds. Tyr42 contributes significantly to guanine binding through a parallel face-to-face stacking interaction. Tyr45, often referred to as the lid of the guanine-binding site, does not contribute to the binding of the base.

Additivity of protein-guanine interactions in ribonuclease T1

It has been established that Tyr-42, Tyr-45, and Glu-46 take part in a structural motif that renders guanine specificity to ribonuclease T1. We report on the impact of Tyr-42, Tyr-45, and Glu-46 substitutions on the guanine specificity of RNase T1. The Y42A and E46A mutations profoundly affect substrate binding. No such effect is observed for Y45A RNase T1. From the kinetics of the Y42A/Y45A and Y42A/E46A double mutants, we conclude that these pairs of residues contribute to guanine specificity in a mutually independent way. From our results, it appears that the energetic contribution of aromatic face-to-face stacking interactions may be significant if polycyclic molecules, such as guanine, are involved.

Limits of NMR structure determination using variable target function calculations: ribonuclease T1, a case study

Limits of NMR structure determination using multidimensional NMR spectroscopy, variable target function calculations and relaxation matrix analysis were explored using the model protein ribonuclease T1 (RNase T1). The enzyme consists of 104 amino acid residues and has a molecular mass of approximately 11 kDa. Primary experimental data comprise 1856 assigned NOE intensities, 493 3J coupling constants and 62 values of amid proton exchange rates. From these data, 2580 distance bounds, 168 allowed ranges for torsional angles and stereospecific assignments for 75% of beta-methylene protons as well as for 80% of diastereotopic methyl groups were derived. Whenever possible, the distance restraints were refined in a relaxation matrix analysis including amid proton exchange data for improvement of lower distance limits. Description of side-chain conformations were based on various models of motional averaging of 3J coupling constants. The final structure ensemble was selected from the starting ensemble comparing the global precision of structures with order parameters derived from 15N relaxation time measurements. Significant differences between the structure of RNase T1 in solution and in the crystal became apparent from a comparison of the two highly resolved structures.

All in one: a highly detailed rotamer library improves both accuracy and speed in the modelling of sidechains by dead-end elimination

BACKGROUND

About a decade ago, the concept of rotamer libraries was introduced to model sidechains given known mainchain coordinates. Since then, several groups have developed methods to handle the challenging combinatorial problem that is faced when searching rotamer libraries. To avoid a combinatorial explosion, the dead-end elimination method detects and eliminates rotamers that cannot be members of the global minimum energy conformation (GMEC). Several groups have applied and further developed this method in the fields of homology modelling and protein design.

RESULTS

This work addresses at the same time increased prediction accuracy and calculation speed improvements. The proposed enhancements allow the elimination of more than one-third of the possible rotameric states before applying the dead-end elimination method. This is achieved by using a highly detailed rotamer library allowing the safe application of an energy-based rejection criterion without risking the elimination of a GMEC rotamer. As a result, we gain both in modelling accuracy and in computational speed. Being completely automated, the current implementation of the dead-end elimination prediction of protein sidechains can be applied to the modelling of sidechains of proteins of any size on the high-end computer systems currently used in molecular modelling. The improved accuracy is highlighted in a comparative study on a collection of proteins of varying size for which score results have previously been published by multiple groups. Furthermore, we propose a new validation method for the scoring of the modelled structure versus the experimental data based upon the volume overlap of the predicted and observed sidechains. This overlap criterion is discussed in relation to the classic RMSD and the frequently used +/- 40 degrees window in comparing chi 1 and chi 2 angles.

CONCLUSIONS

We have shown that a very detailed library allows the introduction of a safe energy threshold rejection criterion, thereby increasing both the execution speed and the accuracy of the modelling program. We speculate that the current method will allow the sidechain prediction of medium-sized proteins and complex protein interfaces involving up to 150 residues on low-end desktop computers.

A two-binding-site kinetic model for the ribonuclease-T1-catalysed transesterification of dinucleoside phosphate substrates

Ribonucleases have been found to have subsites that confer large rate enhancements but do not contribute to substrate binding. In this study, we present a kinetic model that formally explains how subsite binding energy is converted into chemical activation energy. The proposed mechanism takes into account a primary specificity site and a subsite, both of which must be occupied for chemical turnover. An unstable reaction intermediate is formed upon binding of the polymeric substrate monomers at the corresponding subsites. The structure of this reaction intermediate resembles the transition state of the catalysed transphosphorylation reaction. Similar mechanisms may be used by other depolymerizing enzymes including nucleases, glycosidases, and proteases.

The structures of RNase A complexed with 3'-CMP and d(CpA): active site conformation and conserved water molecules

The interactions of RNase A with cytidine 3'-monophosphate (3'-CMP) and deoxycytidyl-3',5'-deoxyadenosine (d(CpA)) were analyzed by X-ray crystallography. The 3'-CMP complex and the native structure were determined from trigonal crystals, and the d(CpA) complex from monoclinic crystals. The differences between the overall structures are concentrated in loop regions and are relatively small. The protein-inhibitor contacts are interpreted in terms of the catalytic mechanism. The general base His 12 interacts with the 2' oxygen, as does the electrostatic catalyst Lys 41. The general acid His 119 has 2 conformations (A and B) in the native structure and is found in, respectively, the A and the B conformation in the d(CpA) and the 3'-CMP complex. From the present structures and from a comparison with RNase T1, we propose that His 119 is active in the A conformation. The structure of the d(CpA) complex permits a detailed analysis of the downstream binding site, which includes His 119 and Asn 71. The comparison of the present RNase A structures with an inhibitor complex of RNase T1 shows that there are important similarities in the active sites of these 2 enzymes, despite the absence of any sequence homology. The water molecules were analyzed in order to identify conserved water sites. Seventeen water sites were found to be conserved in RNase A structures from 5 different space groups. It is proposed that 7 of those water molecules play a role in the binding of the N-terminal helix to the rest of the protein and in the stabilization of the active site.