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

Ultra-high resolution X-ray structure of orthorhombic bovine pancreatic Ribonuclease A at 100K

The crystal structure of orthorhombic Bovine Pancreatic Ribonuclease A has been determined to 0.85 Å resolution using low temperature, 100 K, synchrotron X-ray data collected at 16000 keV (λ = 0.77 Å). This is the first ultra-high-resolution structure of a native form of Ribonuclease A to be reported. Refinement carried out with anisotropic displacement parameters, stereochemical restraints, inclusion of H atoms in calculated positions, five [Formula: see text] moieties, eleven ethanol molecules and 293 water molecules, converged with final R values of R1(Free) = 0.129 (4279 reflections) and R1 = 0.112 (85,346 reflections). The refined structure was deposited in the Protein Data Bank as structure 7p4r. Conserved waters, using four high resolution structures, have been investigated. Cluster analysis identified clusters of water molecules that are associated with the active site of Bovine Ribonuclease A. Particular attention has been paid to making detailed comparisons between the present structure and other high quality Bovine Pancreatic Ribonuclease A X-ray crystal structures with special reference to the deposited classic monoclinic structure 3RN3 Howlin et al. (Acta Crystallogr A 45:851-861, 1989). Detailed studies of various aspects of hydrogen bonding and conformation have been carried out with particular reference to active site residues Lys-1, Lys-7, Gln-11, His-12, Lys-41, Asn-44, Thr-45, Lys-66, His-119 and Ser-123. For the two histidine residues in the active site the initial electron density map gives a clear confirmation that the position of His-12 is very similar in the orthorhombic structure to that in 3RN3. In 3RN3 His-119 exhibited poor electron density which was modelled and refined as two distinct sites, A (65%) and B (35%) but with respect to His-119 in the present ultra-high resolution orthorhombic structure there is clear electron density which was modelled and refined as a single conformation distinct from either conformation A or B in 3RN3. Other points of interest include Serine-32 which is disordered at the end of the sidechain in the present orthorhombic form but has been modelled as a single form in 3RN3. Lysine-66: there is density indicating a possible conformation for this residue. However, the density is relatively weak, and the conformation is unclear. Three types of amino acid representation in the ultra-high resolution electron density are examined: (i) sharp with very clearly resolved features, for example Lys-37; (ii) well resolved but clearly divided into two conformations which are well behaved in the refinement, both having high quality geometry, for example Tyr-76; (iii) poor density and difficult or impossible to model, an example is Lys-31 for which density is missing except for Cβ. The side chains of Gln-11, His-12, Lys-41, Thr-45 and His-119 are generally recognised as being closely involved in the enzyme activity. It has also been suggested that Lys-7, Asp-44, Lys-66, Phe-120, Asp-121 and Ser-123 may also have possible roles in this mechanism. A molecular dynamics study on both structures has investigated the conformations of His-119 which was modelled as two conformations in 3RN3 but is observed to have a single clearly defined conformation in the present orthorhombic structure. MD has also been used to investigate Lys-31, Lys-41 and Ser32. The form of the Ribonuclease A enzyme used in both the present study and in 3RN3 (Howlin et al. in Acta Crystallogr A 45:851-861, 1989) includes a sulphate anion which occupies approximately the same location as the [Formula: see text] phosphate group in protein nucleotide complexes (Borkakoti et al. in J Mol Biol 169:743-755, 1983). The present structure contains 5 [Formula: see text] groups SO41151-SO41155 two of which, SO41152 and SO41153 are disordered, SO41152 being in the active site, and 11 EtOH molecules, EOH A 201-EOH A 211 all of which have good geometry. H atoms were built into the EtOH molecules geometrically. Illustrations of these features in the present structure are included here. The sulphates are presumably present in the material purchased for use in the present study. 293 water molecules are included in the present structure compared to 134 in 3RN3 (Howlin et al. in Acta Crystallogr A 45:851-861, 1989).

Exploring the RNase A scaffold to combine catalytic and antimicrobial activities. Structural characterization of RNase 3/1 chimeras

Design of novel antibiotics to fight antimicrobial resistance is one of the first global health priorities. Novel protein-based strategies come out as alternative therapies. Based on the structure-function knowledge of the RNase A superfamily we have engineered a chimera that combines RNase 1 highest catalytic activity with RNase 3 unique antipathogen properties. A first construct (RNase 3/1-v1) was successfully designed with a catalytic activity 40-fold higher than RNase 3, but alas in detriment of its anti-pathogenic activity. Next, two new versions of the original chimeric protein were created showing improvement in the antimicrobial activity. Both second generation versions (RNases 3/1-v2 and -v3) incorporated a loop characteristic of RNase 3 (L7), associated to antimicrobial activity. Last, removal of an RNase 1 flexible loop (L1) in the third version enhanced its antimicrobial properties and catalytic efficiency. Here we solved the 3D structures of the three chimeras at atomic resolution by X-ray crystallography. Structural analysis outlined the key functional regions. Prediction by molecular docking of the protein chimera in complex with dinucleotides highlighted the contribution of the C-terminal region to shape the substrate binding cavity and determine the base selectivity and catalytic efficiency. Nonetheless, the structures that incorporated the key features related to RNase 3 antimicrobial activity retained the overall RNase 1 active site conformation together with the essential structural elements for binding to the human ribonuclease inhibitor (RNHI), ensuring non-cytotoxicity. Results will guide us in the design of the best RNase pharmacophore for anti-infective therapies.

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.

Spectroscopic/Computational Characterization and the X-ray Structure of the Adduct of the VIVO-Picolinato Complex with RNase A

The structure, stability, and enzymatic activity of the adduct formed upon the reaction of the V-picolinato (pic) complex [VIVO(pic)2(H2O)], with an octahedral geometry and the water ligand in cis to the V═O group, with the bovine pancreatic ribonuclease (RNase A) were studied. While electrospray ionization-mass spectrometry, circular dichroism, and ultraviolet-visible absorption spectroscopy substantiate the interaction between the metal moiety and RNase A, electron paramagnetic resonance (EPR) allows us to determine that a carboxylate group, stemming from Asp or Glu residues, and imidazole nitrogen from His residues are involved in the V binding at acidic and physiological pH, respectively. Crystallographic data demonstrate that the VIVO(pic)2 moiety coordinates the side chain of Glu111 of RNase A, by substituting the equatorial water molecule at acidic pH. Computational methods confirm that Glu111 is the most affine residue and interacts favorably with the OC-6-23-Δ enantiomer establishing an extended network of hydrogen bonds and van der Waals stabilizations. By increasing the pH around neutrality, with the deprotonation of histidine side chains, the binding of the V complex to His105 and His119 could occur, with that to His105 which should be preferred when compared to that to the catalytically important His119. The binding of the V compound affects the enzymatic activity of RNase A, but it does not alter its overall structure and stability.

The structural features of an ancient ribonuclease from Salmo salar reveal an intriguing case of auto-inhibition

The superfamily of vertebrate ribonucleases, a large group of evolutionarily related proteins, continues to provide interesting structural and functional information. In particular, the crystal structure of SS-RNase-2 from Salmo salar (SS2), here presented, has revealed a novel auto-inhibition mechanism that enriches the number of inhibition strategies observed in some members of the family. Within an essentially unmodified RNase folding, the SS2 active site cleft is in part obstructed by the collapse of an extra pentapeptide inserted in the C-terminal region. This unexpected intrusion alters the organization of the catalytic triad by pushing one catalytic histidine off the pocket. Possible mechanisms to remove the active site obstruction have also been studied through the production of two mutants that provide useful information on the functionality of this intriguing version of the ribonuclease superfamily.

Bridging solvent molecules mediate RNase A - Ligand binding

Due to its high catalytic activity and readily available supply, ribonuclease A (RNase A) has become a pivotal enzyme in the history of protein science. Moreover, this great interest has carried over to computational chemistry and molecular dynamics, where RNase A has become a model system for various types of studies, all the while being an important drug design target in its own right. Here, we present a detailed molecular dynamics study of RNase–ligand binding involving 22 compounds, spanning nearly five orders of magnitude in affinity, and totaling 8.8 μs of sampling with the standard Amber parameters and an additional 8.8 μs of sampling with a modified potential. We show that short-lived, solvent-mediated bridging interactions are crucial to RNase–ligand binding. We characterize the behavior of bridging solvent molecules, uncovering a power-law dependence between the lifetime of a solvent bridge and the probability of its occurrence. We also demonstrate that from an energetic perspective, bridging solvent in RNase A–ligand binding behaves like part of the enzyme, rather than the ligands. Moreover, we describe an automated pipeline for the detection and processing of bridging interactions, and offer an independent assessment of the performance of the state-of-the-art fixed-charge force fields. Thus, our work has broad implications for drug design and computational chemistry in general.

Evolutionary Trends in RNA Base Selectivity Within the RNase A Superfamily

There is a growing interest in the pharmaceutical industry to design novel tailored drugs for RNA targeting. The vertebrate-specific RNase A superfamily is nowadays one of the best characterized family of enzymes and comprises proteins involved in host defense with specific cytotoxic and immune-modulatory properties. We observe within the family a structural variability at the substrate-binding site associated to a diversification of biological properties. In this work, we have analyzed the enzyme specificity at the secondary base binding site. Towards this end, we have performed a kinetic characterization of the canonical RNase types together with a molecular dynamic simulation of selected representative family members. The RNases' catalytic activity and binding interactions have been compared using UpA, UpG and UpI dinucleotides. Our results highlight an evolutionary trend from lower to higher order vertebrates towards an enhanced discrimination power of selectivity for adenine respect to guanine at the secondary base binding site (B2). Interestingly, the shift from guanine to adenine preference is achieved in all the studied family members by equivalent residues through distinct interaction modes. We can identify specific polar and charged side chains that selectively interact with donor or acceptor purine groups. Overall, we observe selective bidentate polar and electrostatic interactions: Asn to N1/N6 and N6/N7 adenine groups in mammals versus Glu/Asp and Arg to N1/N2, N1/O6 and O6/N7 guanine groups in non-mammals. In addition, kinetic and molecular dynamics comparative results on UpG versus UpI emphasize the main contribution of Glu/Asp interactions to N1/N2 group for guanine selectivity in lower order vertebrates. A close inspection at the B2 binding pocket also highlights the principal contribution of the protein ß6 and L4 loop regions. Significant differences in the orientation and extension of the L4 loop could explain how the same residues can participate in alternative binding modes. The analysis suggests that within the RNase A superfamily an evolution pressure has taken place at the B2 secondary binding site to provide novel substrate-recognition patterns. We are confident that a better knowledge of the enzymes' nucleotide recognition pattern would contribute to identify their physiological substrate and eventually design applied therapies to modulate their biological functions.

DRoP: Automated detection of conserved solvent-binding sites on proteins

Water and ligand binding play critical roles in the structure and function of proteins, yet their binding sites and significance are difficult to predict a priori. Multiple solvent crystal structures (MSCS) is a method where several X-ray crystal structures are solved, each in a unique solvent environment, with organic molecules that serve as probes of the protein surface for sites evolved to bind ligands, while the first hydration shell is essentially maintained. When superimposed, these structures contain a vast amount of information regarding hot spots of protein-protein or protein-ligand interactions, as well as conserved water-binding sites retained with the change in solvent properties. Optimized mining of this information requires reliable structural data and a consistent, objective analysis tool. Detection of related solvent positions (DRoP) was developed to automatically organize and rank the water or small organic molecule binding sites within a given set of structures. It is a flexible tool that can also be used in conserved water analysis given multiple structures of any protein independent of the MSCS method. The DRoP output is an HTML format list of the solvent sites ordered by conservation rank in its population within the set of structures, along with renumbered and recolored PDB files for visualization and facile analysis. Here, we present a previously unpublished set of MSCS structures of bovine pancreatic ribonuclease A (RNase A) and use it together with published structures to illustrate the capabilities of DRoP.

Ligand-protein target screening from cell matrices using reactive desorption electrospray ionization-mass spectrometry via a native-denatured exchange approach

Native mass spectrometry has been recognized as a powerful tool for probing interactions between small molecules, such as drugs and natural products, and target proteins. However, the presence of heterogeneous proteins and metabolites in real biological systems can alter the conformations of target proteins or compete with candidate ligands, thus necessitating a method for measuring binding stoichiometries in matrices aside from the extensively used pure/recombinant protein systems. Furthermore, some small molecule-protein interactions have a transient and low-affinity nature and thus can be mis-assigned as nonspecific binding complexes that are often formed during the native ESI process. A native-denatured exchange (NDX) approach was recently developed using a reactive desorption electrospray ionization-mass spectrometer (DESI-MS) setup to screen specific interacting partners. The method works by gradually increasing the composition of denaturing solvents contained in the DESI spray and thus conferring a switch from a native to denatured ionization environment. This change impairs three-dimensional structures of target proteins and disrupts specific ligand-protein interactions, leading to decreased holo/apo ratios. In contrast, ligand-protein complexes exhibiting different trends are assigned as nonspecific interactions. Herein, we applied the NDX approach to probe specific ligand-protein interactions in biological matrices. We first used mixtures of model ligands and proteins to examine the use of reactive DESI-MS in recognizing ligand-target binding in mixtures. Subsequently, we used the NDX approach to analyze binding affinity curves of ligands to target proteins spiked in cell lysates with the aid of size exclusion chromatography and demonstrated its use in probing specific ligand-protein interactions from cell matrices.

Characterization of an RNase with two catalytic centers. Human RNase6 catalytic and phosphate-binding site arrangement favors the endonuclease cleavage of polymeric substrates

BACKGROUND

Human RNase6 is a small cationic antimicrobial protein that belongs to the vertebrate RNaseA superfamily. All members share a common catalytic mechanism, which involves a conserved catalytic triad, constituted by two histidines and a lysine (His15/His122/Lys38 in RNase6 corresponding to His12/His119/Lys41 in RNaseA). Recently, our first crystal structure of human RNase6 identified an additional His pair (His36/His39) and suggested the presence of a secondary active site.

METHODS

In this work we have explored RNase6 and RNaseA subsite architecture by X-ray crystallography, site-directed mutagenesis and kinetic characterization.

RESULTS

The analysis of two novel crystal structures of RNase6 in complex with phosphate anions at atomic resolution locates a total of nine binding sites and reveals the contribution of Lys87 to phosphate-binding at the secondary active center. Contribution of the second catalytic triad residues to the enzyme activity is confirmed by mutagenesis. RNase6 catalytic site architecture has been compared with an RNaseA engineered variant where a phosphate-binding subsite is converted into a secondary catalytic center (RNaseA-K7H/R10H).

CONCLUSIONS

We have identified the residues that participate in RNase6 second catalytic triad (His36/His39/Lys87) and secondary phosphate-binding sites. To note, residues His39 and Lys87 are unique within higher primates. The RNaseA/RNase6 side-by-side comparison correlates the presence of a dual active site in RNase6 with a favored endonuclease-type cleavage pattern.

GENERAL SIGNIFICANCE

An RNase dual catalytic and extended binding site arrangement facilitates the cleavage of polymeric substrates. This is the first report of the presence of two catalytic centers in a single monomer within the RNaseA superfamily.

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.

Chemical synthesis and enzymatic properties of RNase A analogues designed to enhance second-step catalytic activity

In this paper, we have used total chemical synthesis of RNase A analogues in order to probe the molecular basis of enzyme catalysis. Our goal was to obligately fill the adenine-binding pocket on the enzyme molecule, and to thus pre-orient the imidazole side chain of His¹¹⁹ in its catalytically productive orientation. Two designed analogues of the RNase A protein molecule that contained an adenine moiety covalently bound to distinct amino acid side chains adjacent to the adenine binding pocket were prepared. A crystal structure of one analogue was determined at 2.3 Å resolution. Kinetic data for RNA transphosporylation and 2',3' cyclic mononucleotide hydrolysis were acquired for the adenine-containing RNase A analogue proteins. As anticipated, the presence of a covalently attached adenine on the enzyme molecule decreased the rate of transphosphorylation and increased the rate of hydrolysis, although the magnitude of the effects was small. This work illustrates the use of total protein synthesis to investigate the chemistry of enzyme catalysis in ways not possible through traditional biochemistry or molecular biology.

A Multidimensional B-Spline Correction for Accurate Modeling Sugar Puckering in QM/MM Simulations

The computational efficiency of approximate quantum mechanical methods allows their use for the construction of multidimensional reaction free energy profiles. It has recently been demonstrated that quantum models based on the neglect of diatomic differential overlap (NNDO) approximation have difficulty modeling deoxyribose and ribose sugar ring puckers and thus limit their predictive value in the study of RNA and DNA systems. A method has been introduced in our previous work to improve the description of the sugar puckering conformational landscape that uses a multidimensional B-spline correction map (BMAP correction) for systems involving intrinsically coupled torsion angles. This method greatly improved the adiabatic potential energy surface profiles of DNA and RNA sugar rings relative to high-level ab initio methods even for highly problematic NDDO-based models. In the present work, a BMAP correction is developed, implemented, and tested in molecular dynamics simulations using the AM1/d-PhoT semiempirical Hamiltonian for biological phosphoryl transfer reactions. Results are presented for gas-phase adiabatic potential energy surfaces of RNA transesterification model reactions and condensed-phase QM/MM free energy surfaces for nonenzymatic and RNase A-catalyzed transesterification reactions. The results show that the BMAP correction is stable, efficient, and leads to improvement in both the potential energy and free energy profiles for the reactions studied, as compared with ab initio and experimental reference data. Exploration of the effect of the size of the quantum mechanical region indicates the best agreement with experimental reaction barriers occurs when the full CpA dinucleotide substrate is treated quantum mechanically with the sugar pucker correction.

Tautomeric stabilities of 4-fluorohistidine shed new light on mechanistic experiments with labeled ribonuclease A

Ribonuclease A is the oldest model for studying enzymatic mechanisms, yet questions remain about proton transfer within the active site. Seminal work by Jackson et al. (Science, 1994) labeled Ribonuclease A with 4-fluorohistidine, concluding that active-site histidines act as general acids and bases. Calculations of 4-fluorohistidine indicate that the π-tautomer is predominant in all simulated environments (by ~17 kJ/mol), strongly suggesting that fluoro-labeled ribonuclease A functions with His119 in π-tautomer. The tautomeric form of His119 during proton transfer and tautomerism as a putative mechanistic step in wild-type RNase A remain open questions and should be considered in future mechanistic studies.

The first crystal structure of human RNase 6 reveals a novel substrate-binding and cleavage site arrangement

We describe the first human RNase 6 crystal structure in complex with sulfate anions. Kinetic analysis, site-directed mutagenesis and molecular dynamics simulations identified novel substrate recognition and cleavage sites.

Human RNase 6 is a cationic secreted protein that belongs to the RNase A superfamily. Its expression is induced in neutrophils and monocytes upon bacterial infection, suggesting a role in host defence. We present here the crystal structure of RNase 6 obtained at 1.72 Å (1 Å=0.1 nm) resolution, which is the first report for the protein 3D structure and thereby setting the basis for functional studies. The structure shows an overall kidney-shaped globular fold shared with the other known family members. Three sulfate anions bound to RNase 6 were found, interacting with residues at the main active site (His¹⁵, His¹²² and Gln¹⁴) and cationic surface-exposed residues (His³⁶, His³⁹, Arg⁶⁶ and His⁶⁷). Kinetic characterization, together with prediction of protein–nucleotide complexes by molecular dynamics, was applied to analyse the RNase 6 substrate nitrogenous base and phosphate selectivity. Our results reveal that, although RNase 6 is a moderate catalyst in comparison with the pancreatic RNase type, its structure includes lineage-specific features that facilitate its activity towards polymeric nucleotide substrates. In particular, enzyme interactions at the substrate 5′ end can provide an endonuclease-type cleavage pattern. Interestingly, the RNase 6 crystal structure revealed a novel secondary active site conformed by the His³⁶–His³⁹ dyad that facilitates the polynucleotide substrate catalysis.

Using NMR spectroscopy to elucidate the role of molecular motions in enzyme function

Conformational motions play an essential role in enzyme function, often facilitating the formation of enzyme-substrate complexes and/or product release. Although considerable debate remains regarding the role of molecular motions in the conversion of enzymatic substrates to products, numerous examples have found motions to be crucial for optimization of enzyme scaffolds, effective substrate binding, and product dissociation. Conformational fluctuations are often rate-limiting to enzyme catalysis, primarily through product release, with the chemical reaction occurring much more quickly. As a result, the direct involvement of motions at various stages along the enzyme reaction coordinate remains largely unknown and untested. In the following review, we describe the use of solution NMR techniques designed to probe various timescales of molecular motions and detail examples in which motions play a role in propagating catalytic effects from the active site and directly participate in essential aspects of enzyme function.

fireball/amber: An Efficient Local-Orbital DFT QM/MM Method for Biomolecular Systems

In recent years, quantum mechanics/molecular mechanics (QM/MM) methods have become an important computational tool for the study of chemical reactions and other processes in biomolecular systems. In the QM/MM technique, the active region is described by means of QM calculations, while the remainder of the system is described using a MM approach. Because of the complexity of biomolecules and the desire to achieve converged sampling, it is important that the QM method presents a good balance between accuracy and computational efficiency. Here, we report on the implementation of a QM/MM technique that combines a DFT approach specially designed for the study of complex systems using first-principles molecular dynamics simulations (fireball) with the amber force fields and simulation programs. We also present examples of the application of this QM/MM approach to three representative biomolecular systems: the analysis of the effect of electrostatic embedding in the behavior of a salt bridge between an aspartic acid and a lysine residue, a study of the intermediate states for the triosephosphate isomerase catalyzed conversion of dihydroxyacetone phosphate into glyceraldehyde 3-phosphate, and the detailed description, using DFT QM/MM molecular dynamics, of the cleavage of a phosphodiester bond in RNA catalyzed by the enzyme RNase A.

Integration of kinetic isotope effect analyses to elucidate ribonuclease mechanism

The well-studied mechanism of ribonuclease A is believed to involve concerted general acid-base catalysis by two histidine residues, His12 and His119. The basic features of this mechanism are often cited to explain rate enhancement by both protein and RNA enzymes that catalyze RNA 2'-O-transphosphorylation. Recent kinetic isotope effect analyses and computational studies are providing a more chemically detailed description of the mechanism of RNase A and the rate limiting transition state. Overall, the results support an asynchronous mechanism for both solution and ribonuclease catalyzed reactions in which breakdown of a transient dianoinic phosphorane intermediate by 5'OP bond cleavage is rate limiting. Relative to non-enzymatic reactions catalyzed by specific base, a smaller KIE on the 5'O leaving group and a less negative βLG are observed for RNase A catalysis. Quantum mechanical calculations consistent with these data support a model in which electrostatic and H-bonding interactions with the non-bridging oxygens and proton transfer from His119 render departure of the 5'O less advanced and stabilize charge buildup in the transition state. Both experiment and computation indicate advanced 2'OP bond formation in the rate limiting transition state. However, this feature makes it difficult to resolve the chemical steps involved in 2'O activation. Thus, modeling the transition state for RNase A catalysis underscores those elements of its chemical mechanism that are well resolved, as well as highlighting those where ambiguity remains. This article is part of a Special Issue entitled: Enzyme Transition States from Theory and Experiment.

Sequence-specific backbone (1)H, (13)C, and (15)N resonance assignments of human ribonuclease 4

Human ribonuclease 4 (RNase 4) is the most evolutionarily conserved member of the 8 canonical human pancreatic-like RNases, showing more than 90% identity with bovine and porcine homologues. The enzyme displays ribonucleolytic activity with a strong preference for uracil-containing RNA substrates, a feature only shared with human eosinophil derived-neurotoxin (EDN, or RNase 2) and eosinophil cationic protein (ECP, or RNase 3). It is also the shortest member of the human family, with a significantly truncated C-terminal tail. Its unique active-site pocket and high degree of conservation among vertebrates suggest that the enzyme plays a crucial biological function. Here, we report on the (1)H, (13)C and (15)N backbone resonance assignments of RNase 4, providing means to characterize its molecular function at the atomic level by NMR.

Characterization of disulfide bonds by planned digestion and tandem mass spectrometry

The identification of disulfide bonds provides critical information regarding the structure and function of a protein and is a key aspect in understanding signaling cascades in biological systems. Recent proteomic approaches using digestion enzymes have facilitated the characterization of disulfide-bonds and/or oxidized products from cysteine residues, although these methods have limitations in the application of MS/MS. For example, protein digestion to obtain the native form of disulfide bonds results in short lengths of amino acids, which can cause ambiguous MS/MS analysis due to false positive identifications. In this study we propose a new approach, termed planned digestion, to obtain sufficient amino acid lengths after cleavage for proteomic approaches. Application of the DBond software to planned digestion of specific proteins accurately identified disulfide-linked peptides. RNase A was used as a model protein in this study because the disulfide bonds of this protein have been well characterized. Application of this approach to peptides digested with Asp-N/C (chemical digestion) and trypsin under acid hydrolysis conditions identified the four native disulfide bonds of RNase A. Missed cleavages introduced by trypsin treatment for only 3 hours generated sufficient lengths of amino acids for identification of the disulfide bonds. Analysis using MS/MS successfully showed additional fragmentation patterns that are cleavage products of S-S and C-S bonds of disulfide-linkage peptides. These fragmentation patterns generate thioaldehydes, persulfide, and dehydroalanine. This approach of planned digestion with missed cleavages using the DBond algorithm could be applied to other proteins to determine their disulfide linkage and the oxidation patterns of cysteine residues.

Interpretation of pH-activity profiles for acid-base catalysis from molecular simulations

The measurement of reaction rate as a function of pH provides essential information about mechanism. These rates are sensitive to the pK(a) values of amino acids directly involved in catalysis that are often shifted by the enzyme active site environment. Experimentally observed pH-rate profiles are usually interpreted using simple kinetic models that allow estimation of "apparent pK(a)" values of presumed general acid and base catalysts. One of the underlying assumptions in these models is that the protonation states are uncorrelated. In this work, we introduce the use of constant pH molecular dynamics simulations in explicit solvent (CpHMD) with replica exchange in the pH-dimension (pH-REMD) as a tool to aid in the interpretation of pH-activity data of enzymes and to test the validity of different kinetic models. We apply the methods to RNase A, a prototype acid-base catalyst, to predict the macroscopic and microscopic pK(a) values, as well as the shape of the pH-rate profile. Results for apo and cCMP-bound RNase A agree well with available experimental data and suggest that deprotonation of the general acid and protonation of the general base are not strongly coupled in transphosphorylation and hydrolysis steps. Stronger coupling, however, is predicted for the Lys41 and His119 protonation states in apo RNase A, leading to the requirement for a microscopic kinetic model. This type of analysis may be important for other catalytic systems where the active forms of the implicated general acid and base are oppositely charged and more highly correlated. These results suggest a new way for CpHMD/pH-REMD simulations to bridge the gap with experiments to provide a molecular-level interpretation of pH-activity data in studies of enzyme mechanisms.

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.

QM/MM simulation (B3LYP) of the RNase A cleavage-transesterification reaction supports a triester A(N) + D(N) associative mechanism with an O2' H internal proton transfer

The mechanism of the backbone cleavage-transesterification step of the RNase A enzyme remains controversial even after 60 years of study. We report quantum mechanics/molecule mechanics (QM/MM) free energy calculations for two optimized reaction paths based on an analysis of all structural data and identified by a search for reaction coordinates using a reliable quantum chemistry method (B3LYP), equilibrated structural optimizations, and free energy estimations. Both paths are initiated by nucleophilic attack of the ribose O2' oxygen on the neighboring diester phosphate bond, and both reach the same product state (PS) (a O3'-O2' cyclic phosphate and a O5' hydroxyl terminated fragment). Path 1, resembles the widely accepted dianionic transition-state (TS) general acid (His119)/base (His12) classical mechanism. However, this path has a barrier (25 kcal/mol) higher than that of the rate-limiting hydrolysis step and a very loose TS. In Path 2, the proton initially coordinating the O2' migrates to the nonbridging O1P in the initial reaction path rather than directly to the general base resulting in a triester (substrate as base) AN + DN mechanism with a monoanionic weakly stable intermediate. The structures in the transition region are associative with low barriers (TS1 10, TS2 7.5 kcal/mol). The Path 2 mechanism is consistent with the many results from enzyme and buffer catalyzed and uncatalyzed analog reactions and leads to a PS consistent with the reactive state for the following hydrolysis step. The differences between the consistently estimated barriers in Path 1 and 2 lead to a 10(11) difference in rate strongly supporting the less accepted triester mechanism.

DRoP: a water analysis program identifies Ras-GTP-specific pathway of communication between membrane-interacting regions and the active site

Ras GTPase mediates several cellular signal transduction pathways and is found mutated in a large number of cancers. It is active in the GTP-bound state, where it interacts with effector proteins, and at rest in the GDP-bound state. The catalytic domain is tethered to the membrane, with which it interacts in a nucleotide-dependent manner. Here we present the program Detection of Related Solvent Positions (DRoP) for crystallographic water analysis on protein surfaces and use it to study Ras. DRoP reads and superimposes multiple Protein Data Bank coordinates, transfers symmetry-related water molecules to the position closest to the protein surface, and ranks the waters according to how well conserved and tightly clustered they are in the set of structures. Coloring according to this rank allows visualization of the results. The effector-binding region of Ras is hydrated with highly conserved water molecules at the interface between the P-loop, switch I, and switch II, as well as at the Raf-RBD binding pocket. Furthermore, we discovered a new conserved water-mediated H-bonding network present in Ras-GTP, but not in Ras-GDP, that links the nucleotide sensor residues R161 and R164 on helix 5 to the active site. The double mutant RasN85A/N86A, where the final link between helix 5 and the nucleotide is not possible, is a severely impaired enzyme, while the single mutant RasN86A, with partial connection to the active site, has a wild-type hydrolysis rate. DRoP was instrumental in determining the water-mediated connectivity networks that link two lobes of the catalytic domain in Ras.

Experimental and computational analysis of the transition state for ribonuclease A-catalyzed RNA 2'-O-transphosphorylation

Enzymes function by stabilizing reaction transition states; therefore, comparison of the transition states of enzymatic and nonenzymatic model reactions can provide insight into biological catalysis. Catalysis of RNA 2'-O-transphosphorylation by ribonuclease A is proposed to involve electrostatic stabilization and acid/base catalysis, although the structure of the rate-limiting transition state is uncertain. Here, we describe coordinated kinetic isotope effect (KIE) analyses, molecular dynamics simulations, and quantum mechanical calculations to model the transition state and mechanism of RNase A. Comparison of the (18)O KIEs on the 2'O nucleophile, 5'O leaving group, and nonbridging phosphoryl oxygens for RNase A to values observed for hydronium- or hydroxide-catalyzed reactions indicate a late anionic transition state. Molecular dynamics simulations using an anionic phosphorane transition state mimic suggest that H-bonding by protonated His12 and Lys41 stabilizes the transition state by neutralizing the negative charge on the nonbridging phosphoryl oxygens. Quantum mechanical calculations consistent with the experimental KIEs indicate that expulsion of the 5'O remains an integral feature of the rate-limiting step both on and off the enzyme. Electrostatic interactions with positively charged amino acid site chains (His12/Lys41), together with proton transfer from His119, render departure of the 5'O less advanced compared with the solution reaction and stabilize charge buildup in the transition state. The ability to obtain a chemically detailed description of 2'-O-transphosphorylation transition states provides an opportunity to advance our understanding of biological catalysis significantly by determining how the catalytic modes and active site environments of phosphoryl transferases influence transition state structure.

Structural and functional importance of local and global conformational fluctuations in the RNase A superfamily

Understanding the relationship between protein structure and flexibility is of utmost importance for deciphering the tremendous rates of reactions catalyzed by enzyme biocatalysts. It has been postulated that protein homologs have evolved similar dynamic fluctuations to promote catalytic function, a property that would presumably be encoded in their structural fold. Using one of the best-characterized enzyme systems of the past century, we explore this hypothesis by comparing the numerous and diverse flexibility reports available for a number of structural and functional homologs of the pancreatic-like RNase A superfamily. Using examples from the literature and from our own work, we cover recent and historical evidence pertaining to the highly dynamic nature of this important structural fold, as well as the presumed importance of local and global concerted motions on the ribonucleolytic function. This minireview does not pretend to cover the overwhelming RNase A literature in a comprehensive manner; rather, efforts have been made to focus on the characterization of multiple timescale motions observed in the free and/or ligand-bound structural homologs as they proceed along the reaction coordinates. Although each characterized enzyme of this architectural fold shows unique motional features on a local scale, accumulating evidence from X-ray crystallography, NMR spectroscopy and molecular dynamics simulations suggests that global dynamic fluctuations, such as the functionally relevant hinge-bending motion observed in the prototypical RNase A, are shared between homologs of the pancreatic-like RNase superfamily. These observations support the hypothesis that analogous dynamic residue clusters are evolutionarily conserved among structural and functional homologs catalyzing similar enzymatic reactions.

Nucleotide binding architecture for secreted cytotoxic endoribonucleases

Vertebrate secreted RNases are small cationic protein endowed with an endoribonuclease activity that belong to the RNase A superfamily and display diverse cytotoxic activities. In an effort to unravel their mechanism of action, we have analysed their nucleotide binding recognition patterns. General shared features with other nucleotide binding proteins were deduced from overall statistics on the available structure complexes at the Protein Data Bank and compared with the particularities of selected representative endoribonuclease families. Results were compared with other endoribonuclease representative families and with the overall protein-nucleotide interaction features. Preferred amino acids and atom types involved in pair bonding interactions were identified, defining the spatial motives for phosphate, base and ribose building blocks. Together with the conserved catalytic triad at the active site, variability was observed for secondary binding subsites that may contribute to the proper substrate alignment and could explain the distinct substrate preference patterns. Highly conserved binding patterns were identified for the pyrimidine and purine subsites at the main and secondary base subsites. Particular substitution could be ascribed to specific adenine or guanine specificities. Distribution of evolutionary conserved residues were compared to search for the structure determinants that underlie their diverse catalytic efficiency and those that may account for putative physiological substrate targets or other non-catalytic biological activities that contribute to the antipathogen role of the RNases involved in the host defence system. A side by side comparison with another endoribonuclease superfamily of secreted cytotoxic proteins, the microbial RNases, was carried on to analyse the common features and peculiarities that rule their substrate recognition. The data provides the structural basis for the development of applied therapies targeting cellular nucleotide polymers.

Crystal structures of murine angiogenin-2 and -3-probing 'structure--function' relationships amongst angiogenin homologues

Angiogenin (Ang) is a potent inducer of neovascularization. Point mutations in human Ang have been linked to cancer progression and two neurodegenerative diseases: amyotrophic lateral sclerosis and Parkinson's disease. Intensive structural and functional analyses of Ang have been paramount in assigning functions to this novel homologue of bovine pancreatic RNase A. However, inhibitor-binding studies with crystalline Ang (for designing potential anti-cancer drugs) have been hampered as a result of the inaccessibility of the active site. Experiments with the murine homologues of Ang have not only overcome the obvious practical limitations encountered when studying the role of a human protein in healthy individuals, but also the crystal structures of murine angiogenins (mAng and mAng-4) have revealed themselves to have greater potential for the visualization of small-molecule inhibitor binding at the active site. In the present study, we report the crystal structures of two more murine Ang paralogues, mAng-2 and mAng-3, at 1.6 and 1.8 Å resolution, respectively. These constitute the first crystal structures of an Ang with a zinc ion bound at the active site and provide some insight into the possible mode of inhibition of the ribonucleolytic activity of the enzyme by these divalent cations. Both structures show that the residues forming the putative P₁, B₁ and B₂ subsites occupy positions similar to their counterparts in human Ang and are likely to have conserved roles. However, a less obtrusive conformation of the C-terminal segment in mAng-3 and the presence of a sulfate ion in the B₁ subsite of mAng-2 suggest that these proteins have the potential to be used for inhibitor-binding studies. We also discuss the biological relevance of the structural similarities and differences between the different Ang homologues.

Understanding the molecular mechanism of enzyme dynamics of ribonuclease A through protonation/deprotonation of HIS48

Molecular dynamics simulation is carried out to investigate the enzyme dynamics of RNase A with the HIS48 in three different states (HIP48 (protonated), HID48 (deprotonated), and H48A mutant). Insights derived from the current theoretical study, combined with the available experimental observations, enabled us to provide a microscopic picture for the efficient enzyme dynamics. Specifically, in the "closed" state or HIP48, the N-terminal hinge loop is intact and the enzyme remains in a relatively stable conformation which is preferred for catalytic reaction. Deprotonation of HIS48 induces the denaturing of this hinge-loop into a 3(10)-helix, causing it to break the original interaction network around the loop-1 and drive the partial unfolding of the N-terminal. The enhanced dynamic motion of the N-terminal helix facilitates the release of the catalytic product (the rate limiting step) and speeds up the overall catalytic process. The current study established that HIS49 acts as a modulator for the transformation of conformational states through the perturbing of hydrogen bond networks across loop-1, the N-terminal helix, and other residues nearby. Our study suggests that HIS48 may also serve to transport loop-1's kinetic energy to the reaction center.

Crystal structure of Onconase at 1.1 Å resolution--insights into substrate binding and collective motion

Onconase® (ONC) is an amphibian member of the pancreatic ribonuclease superfamily that is selectively toxic to tumor cells. It is a much less efficient enzyme than the archetypal ribonuclease A and, in an attempt to gain further insight, we report the first atomic resolution crystal structure of ONC, determined in complex with sulfate ions at 100 K. The electron density map is of a quality sufficient to reveal significant nonplanarity in several peptide bonds. The majority of active site residues are very well defined, with the exceptions being Lys31 from the catalytic triad and Lys33 from the B₁ subsite, which are relatively mobile but rigidify upon nucleotide binding. Cryocooling causes a compaction of the unit cell and the protein contained within. This is principally the result of an inward movement of one of the lobes of the enzyme (lobe 2), which also narrows the active site cleft. Binding a nucleotide in place of sulfate is associated with an approximately perpendicular movement of lobe 2 and has little further effect on the cleft width. Aspects of this deformation are present in the principal axes of anisotropy extracted from C^α atomic displacement parameters, indicating its intrinsic nature. The three lowest-frequency modes of ONC motion predicted by an anisotropic network model are compaction/expansion variations in which lobe 2 is the prime mover. Two of these have high similarity to the cryocooling response and imply that the essential ‘breathing’ motion of ribonuclease A is conserved in ONC. Instead, shifts in conformational equilibria may contribute to the reduced ribonucleolytic activity of ONC.

Site-selective artificial ribonucleases: oligonucleotide conjugates containing multiple imidazole residues in the catalytic domain

Design of site-selective artificial ribonucleases (aRNases) is one of the most challenging tasks in RNA targeting. Here, we designed and studied oligonucleotide-based aRNases containing multiple imidazole residues in the catalytic part and systematically varied structure of cleaving constructs. We demonstrated that the ribonuclease activity of the conjugates is strongly affected by the number of imidazole residues in the catalytic part, the length of a linker between the catalytic imidazole groups of the construct and the oligonucleotide, and the type of anchor group, connecting linker structure and the oligonucleotide. Molecular modeling of the most active aRNases showed that preferable orientation(s) of cleaving constructs strongly depend on the structure of the anchor group and length of the linker. The inclusion of deoxyribothymidine anchor group significantly reduced the probability of cleaving groups to locate near the cleavage site, presumably due to a stacking interaction with the neighbouring nucleotide residue. Altogether the obtained results show that dynamics factors play an important role in site-specific RNA cleavage. Remarkably high cleavage activity was displayed by the conjugates with the most flexible and extended cleaving construct, which presumably provides a better opportunity for imidazole residues to be correctly positioned in the vicinity of scissile phosphodiester bond.

Functional and structural analyses of N-acylsulfonamide-linked dinucleoside inhibitors of RNase A

Molecular probes are useful for both studying and controlling the functions of enzymes and other proteins. The most useful probes have high affinity for their target, along with small size and resistance to degradation. Here, we report on new surrogates for nucleic acids that fulfill these criteria. Isosteres in which phosphoryl [R–O–P(O₂⁻)–O–R′] groups are replaced with N-acylsulfonamidyl [R–C(O)–N⁻–S(O₂)–R′] or sulfonimidyl [R–S(O₂)–N⁻–S(O₂)–R′] groups increase the number of nonbridging oxygens from two (phosphoryl) to three (N-acylsulfonamidyl) or four (sulfonimidyl). Six such isosteres were found to be more potent inhibitors of catalysis by bovine pancreatic RNase A than are parent compounds containing phosphoryl groups. The atomic structures of two RNase A·N-acylsulfonamide complexes were determined at high resolution by X-ray crystallography. The N-acylsulfonamidyl groups were observed to form more hydrogen bonds with active site residues than did the phosphoryl groups in analogous complexes. These data encourage the further development and use of N-acylsulfonamides and sulfonimides as antagonists of nucleic acid-binding proteins.

Nucleotide docking: prediction of reactant state complexes for ribonuclease enzymes

Ribonuclease enzymes (RNases) play key roles in the maturation and metabolism of all RNA molecules. Computational simulations of the processes involved can help to elucidate the underlying enzymatic mechanism and is often employed in a synergistic approach together with biochemical experiments. Theoretical calculations require atomistic details regarding the starting geometries of the molecules involved, which, in the absence of crystallographic data, can only be achieved from computational docking studies. Fortunately, docking algorithms have improved tremendously in recent years, so that reliable structures of enzyme-ligand complexes can now be successfully obtained from computation. However, most docking programs are not particularly optimized for nucleotide docking. In order to assist our studies on the cleavage of RNA by the two most important ribonuclease enzymes, RNase A and RNase H, we evaluated four docking tools-MOE2009, Glide 5.5, QXP-Flo+0802, and Autodock 4.0-for their ability to simulate complexes between these enzymes and RNA oligomers. To validate our results, we analyzed the docking results with respect to the known key interactions between the protein and the nucleotide. In addition, we compared the predicted complexes with X-ray structures of the mutated enzyme as well as with structures obtained from previous calculations. In this manner, we were able to prepare the desired reaction state complex so that it could be used as the starting structure for further DFT/B3LYP QM/MM reaction mechanism studies.

The crystal structure of ribonuclease A in complex with thymidine-3'-monophosphate provides further insight into ligand binding

Thymidine-3'-monophosphate (3'-TMP) is a competitive inhibitor analogue of the 3'-CMP and 3'-UMP natural product inhibitors of bovine pancreatic ribonuclease A (RNase A). Isothermal titration calorimetry experiments show that 3'-TMP binds the enzyme with a dissociation constant (K(d)) of 15 microM making it one of the strongest binding members of the five natural bases found in nucleic acids (A, C, G, T, and U). To further investigate the molecular properties of this potent natural affinity, we have determined the crystal structure of bovine pancreatic RNase A in complex with 3'-TMP at 1.55 A resolution and we have performed NMR binding experiments with 3'-CMP and 3'-TMP. Our results show that binding of 3'-TMP is very similar to other natural and non-natural pyrimidine ligands, demonstrating that single nucleotide affinity is independent of the presence or absence of a 2'-hydroxyl on the ribose moiety of pyrimidines and suggesting that the pyrimidine binding subsite of RNase A is not a significant contributor of inhibitor discrimination. Accumulating evidence suggests that very subtle structural, chemical, and potentially motional variations contribute to ligand discrimination in this enzyme.

Molecular dynamics simulation of bovine pancreatic ribonuclease A-CpA and transition state-like complexes

The mechanisms of enzymes are intimately connected with their overall structure and dynamics in solution. Experimentally, it is considerably challenging to provide detailed atomic level information about the conformational events that occur at different stages along the chemical reaction path. Here, theoretical tools may offer new potential insights that complement those obtained from experiments that may not yield an unambiguous mechanistic interpretation. In this study, we apply molecular dynamics simulations of bovine pancreatic ribonuclease A, an archetype ribonuclease, to study the conformational dynamics, structural relaxation, and differential solvation that occur at discrete stages of the transesterification and cleavage reaction. Simulations were performed with explicit solvation with rigorous electrostatics and utilize recently developed molecular mechanical force field parameters for transphosphorylation and hydrolysis transition state analogues. Herein, we present results for the enzyme complexed with the dinucleotide substrate cytidilyl-3',5'-adenosine (CpA) in the reactant, and transphosphorylation and hydrolysis transition states. A detailed analysis of active site structures and hydrogen-bond patterns is presented and compared. The integrity of the overall backbone structure is preserved in the simulations and supports a mechanism whereby His12 stabilizes accumulating negative charge at the transition states through hydrogen-bond donation to the nonbridge oxygens. Lys41 is shown to be highly versatile along the reaction coordinate and can aid in the stabilization of the dianionic transition state, while being poised to act as a general acid catalyst in the hydrolysis step.

Influence of naturally-occurring 5'-pyrophosphate-linked substituents on the binding of adenylic inhibitors to ribonuclease a: an X-ray crystallographic study

Ribonuclease A is the archetype of a functionally diverse superfamily of vertebrate-specific ribonucleases. Inhibitors of its action have potential use in the elucidation of the in vivo roles of these enzymes and in the treatment of pathologies associated therewith. Derivatives of adenosine 5'-pyrophosphate are the most potent nucleotide-based inhibitors known. Here, we use X-ray crystallography to visualize the binding of four naturally-occurring derivatives that contain 5'-pyrophosphate-linked extensions. 5'-ATP binds with the adenine occupying the B(2) subsite in the manner of an RNA substrate but with the gamma-phosphate at the P(1) subsite. Diadenosine triphosphate (Ap(3)A) binds with the adenine in syn conformation, the beta-phosphate as the principal P(1) subsite ligand and without order beyond the gamma-phosphate. NADPH and NADP(+) bind with the adenine stacked against an alternative rotamer of His119, the 2'-phosphate at the P(1) subsite, and without order beyond the 5'-alpha-phosphate. We also present the structure of the complex formed with pyrophosphate ion. The structural data enable existing kinetic data on the binding of these compounds to a variety of ribonucleases to be rationalized and suggest that as the complexity of the 5'-linked extension increases, the need to avoid unfavorable contacts places limitations on the number of possible binding modes.

Multiple solvent crystal structures of ribonuclease A: an assessment of the method

The multiple solvent crystal structures (MSCS) method uses organic solvents to map the surfaces of proteins. It identifies binding sites and allows for a more thorough examination of protein plasticity and hydration than could be achieved by a single structure. The crystal structures of bovine pancreatic ribonuclease A (RNAse A) soaked in the following organic solvents are presented: 50% dioxane, 50% dimethylformamide, 70% dimethylsulfoxide, 70% 1,6-hexanediol, 70% isopropanol, 50% R,S,R-bisfuran alcohol, 70% t-butanol, 50% trifluoroethanol, or 1.0M trimethylamine-N-oxide. This set of structures is compared with four sets of crystal structures of RNAse A from the protein data bank (PDB) and with the solution NMR structure to assess the validity of previously untested assumptions associated with MSCS analysis. Plasticity from MSCS is the same as from PDB structures obtained in the same crystal form and deviates only at crystal contacts when compared to structures from a diverse set of crystal environments. Furthermore, there is a good correlation between plasticity as observed by MSCS and the dynamic regions seen by NMR. Conserved water binding sites are identified by MSCS to be those that are conserved in the sets of structures taken from the PDB. Comparison of the MSCS structures with inhibitor-bound crystal structures of RNAse A reveals that the organic solvent molecules identify key interactions made by inhibitor molecules, highlighting ligand binding hot-spots in the active site. The present work firmly establishes the relevance of information obtained by MSCS.

The flexibility of a distant loop modulates active site motion and product release in ribonuclease A

The role of the flexible loop 1 in protein conformational motion and in the dissociation of enzymatic product from ribonuclease A (RNase A) was investigated by creation of a chimeric enzyme in which a 6-residue loop 1 from the RNase A homologue, eosinophil cationic protein (ECP), replaced the 12-residue loop 1 in RNase A. The chimera (RNase A(ECP)) experiences only local perturbations in NMR backbone chemical shifts compared to WT RNase A. Many of the flexible residues that were previously identified in WT as involved in an important conformational change now experience no NMR-detected millisecond motions in the chimera. Likewise, binding of the product analogue, 3'-CMP, to RNase A(ECP) results in only minor chemical shift changes in the enzyme similar to what is observed for the H48A mutant of RNase A and in contrast to WT enzyme. For both RNase A(ECP) and H48A there is a 10-fold decrease in the product release rate constant, k(off), compared to WT, in agreement with previous studies indicating the importance of flexibility in RNase A in the overall rate-limiting product release step. Together, these NMR and biochemical experiments provide additional insight into the mechanism of millisecond motions in the RNase A catalytic cycle.

Inhibitor design for ribonuclease A: the binding of two 5'-phosphate uridine analogues

In the quest for the rational design of selective and potent inhibitors for members of the pancreatic ribonuclease A (RNase A) family of biomedical interest, the binding of uridine 5'-phosphate (U5P) and uridine 5'-diphosphate (UDP) to RNase A have been investigated using kinetic studies and X-ray crystallography. Both nucleotides are competitive inhibitors of the enzyme, with K(i) values of 4.0 and 0.65 mM, respectively. They bind to the active site of the enzyme by anchoring two molecules connected to each other by hydrogen bonds and van der Waals interactions. While the first of the inhibitor molecules binds with its nucleobase in the pyrimidinyl-binding subsite, the second is bound at the purine-preferring subsite. The unexpected binding of a pyrimidine at the purine-binding subsite has added new important elements to the rational design approach for the discovery of new potent inhibitors of the RNase A superfamily.

Ribonuclease A homologues of the zebrafish: polymorphism, crystal structures of two representatives and their evolutionary implications

The widespread and functionally varied members of the ribonuclease A (RNase A) superfamily provide an excellent opportunity to study evolutionary forces at work on a conserved protein scaffold. Representatives from the zebrafish are of particular interest as the evolutionary distance from non-ichthyic homologues is large. We conducted an exhaustive survey of available zebrafish DNA sequences and found significant polymorphism among its four known homologues. In an extension of previous nomenclature, the variants have been named RNases ZF-1a–c,-2a–d,-3a–e and-4. We present the first X-ray crystal structures of zebrafish ribonucleases, RNases ZF-1a and-3e at 1.35-and 1.85 Å resolution, respectively. Structure-based clustering with ten other ribonuclease structures indicates greatest similarity to mammalian angiogenins and amphibian ribonucleases, and supports the view that all present-day ribonucleases evolved from a progenitor with three disulphide bonds. In their details, the two structures are intriguing melting-pots of features present in ribonucleases from other vertebrate classes. Whereas in RNase ZF-1a the active site is obstructed by the C-terminal segment (as observed in angiogenin), in RNase ZF-3e the same region is open (as observed in more catalytically efficient homologues). The progenitor of present-day ribonucleases is more likely to have had an obstructive C terminus, and the relatively high similarity (late divergence) of RNases ZF-1 and-3 infers that the active site unblocking event has happened independently in different vertebrate lineages.

Ribozyme catalysis revisited: is water involved?

Enzymatic catalysis by RNA was discovered 25 years ago, yet mechanistic insights are emerging only slowly. Thought to be metalloenzymes at first, some ribozymes proved more versatile than anticipated when shown to utilize their own functional groups for catalysis. Recent evidence suggests that some may also judiciously place structural water molecules to shuttle protons in acid-base catalyzed reactions.

Characterization of enzyme motions by solution NMR relaxation dispersion

In many enzymes, conformational changes that occur along the reaction coordinate can pose a bottleneck to the rate of conversion of substrates to products. Characterization of these rate-limiting protein motions is essential for obtaining a full understanding of enzyme-catalyzed reactions. Solution NMR experiments such as the Carr-Purcell-Meiboom-Gill (CPMG) spin-echo or off-resonance R 1rho pulse sequences enable quantitation of protein motions in the time range of microseconds to milliseconds. These experiments allow characterization of the conformational exchange rate constant, k ex, the equilibrium populations of the relevant conformations, and the chemical shift differences (Deltaomega) between the conformations. The CPMG experiments were applied to the backbone N-H positions of ribonuclease A (RNase A). To probe the role of dynamic processes in the catalytic cycle of RNase A, stable mimics of the apo enzyme (E), enzyme-substrate (ES) complex, and enzyme-product (EP) complex were formed. The results indicate that the ligand has relatively little influence on the kinetics of motion, which occurs at 1700 s (-1) and is the same as both k cat, and the product dissociation rate constant. Instead, the effect of ligand is to stabilize one of the pre-existing conformations. Thus, these NMR experiments indicate that the conformational change in RNase A is ligand-stabilized and does not appear to be ligand-induced. Further evidence for the coupling of motion and enzyme function comes from the similar solvent deuterium kinetic isotope effect on k ex derived from the NMR measurements and k cat from enzyme kinetic studies. This isotope effect of approximately 2 depends linearly on solvent deuterium content suggesting the involvement of a single proton in RNase A motion and function. Moreover, mutation of His48 to alanine eliminates motion in RNase A and decreases the catalytic turnover rate indicating the involvement of His48, which is far from the active site, in coupling motion and function. For the enzyme triosephosphate isomerase (TIM), the opening and closing motion of a highly conserved active site loop (loop 6) has been implicated in many studies to play an important role in the catalytic cycle of the enzyme. Off-resonance R 1rho experiments were performed on TIM, and results were obtained for amino acid residues in the N-terminal (Val167), and C-terminal (Lys174, Thr177) portions of loop 6. The results indicate that all three loop residues move between the open and closed conformation at about 10,000 s (-1), which is the same as the catalytic rate constant. The O (eta) atom of Tyr208 provides a hydrogen bond to stabilize the closed form of loop 6 by interacting with the amide nitrogen of Ala176; these atoms are outside of hydrogen bonding distance in the open form of the enzyme. Mutation of Tyr208 to phenylalanine results in significant loss of catalytic activity but does not appear to alter the kex value of the N-terminal part of loop 6. Instead, removal of this hydrogen bond appears to result in an increase in the equilibrium population of the open conformer of loop 6, thereby resulting in a loss of activity through a shift in the conformational equilibrium of loop 6. Solution NMR relaxation dispersion experiments are powerful experimental tools that can elucidate protein motions with atomic resolution and can provide insight into the role of these motions in biological function.

Ternary borate-nucleoside complex stabilization by ribonuclease A demonstrates phosphate mimicry

Phosphate esters play a central role in cellular energetics, biochemical activation, signal transduction and conformational switching. The structural homology of the borate anion with phosphate, combined with its ability to spontaneously esterify hydroxyl groups, suggested that phosphate ester recognition sites on proteins might exhibit significant affinity for nonenzymatically formed borate esters. (11)B NMR studies and activity measurements on ribonuclease A (RNase A) in the presence of borate and several cytidine analogs demonstrate the formation of a stable ternary RNase A.3'-deoxycytidine-2'-borate ternary complex that mimics the complex formed between RNase A and a 2'-cytidine monophosphate (2'-CMP) inhibitor. Alternatively, no slowly exchanging borate resonance is observed for a ternary RNase A, borate, 2'-deoxycytidine mixture, demonstrating the critical importance of the 2'-hydroxyl group for complex formation. Titration of the ternary complex with 2'-CMP shows that it can displace the bound borate ester with a binding constant that is close to the reported inhibition constant of RNase A by 2'-CMP. RNase A binding of a cyclic cytidine-2',3'-borate ester, which is a structural homolog of the cytidine-2',3'-cyclic phosphate substrate, could also be demonstrated. The apparent dissociation constant for the cytidine-2',3'-borate.RNase A complex is 0.8 mM, which compares with a Michaelis constant of 11 mM for cytidine-2',3'-cyclic phosphate at pH 7, indicating considerably stronger binding. However, the value is 1,000-fold larger than the reported dissociation constant of the RNase A complex with uridine-vanadate. These results are consistent with recent reports suggesting that in situ formation of borate esters that mimic the corresponding phosphate esters supports enzyme catalysis.