Correlation proton magnetic resonance studies at 250 MHz of bovine pancreatic ribonuclease. I. Reinvestigation of the histidine peak assignments

Significance of Histidine Hydrogen-Deuterium Exchange Mass Spectrometry in Protein Structural Biology

Histidine residues play crucial roles in shaping the function and structure of proteins due to their unique ability to act as both acids and bases. In other words, they can serve as proton donors and acceptors at physiological pH. This exceptional property is attributed to the side-chain imidazole ring of histidine residues. Consequently, determining the acid-base dissociation constant (Ka) of histidine imidazole rings in proteins often yields valuable insights into protein functions. Significant efforts have been dedicated to measuring the pKa values of histidine residues in various proteins, with nuclear magnetic resonance (NMR) spectroscopy being the most commonly used technique. However, NMR-based methods encounter challenges in assigning signals to individual imidazole rings and require a substantial amount of proteins. To address these issues associated with NMR-based approaches, a mass-spectrometry-based method known as histidine hydrogen-deuterium exchange mass spectrometry (His-HDX-MS) has been developed. This technique not only determines the pKa values of histidine imidazole groups but also quantifies their solvent accessibility. His-HDX-MS has proven effective across diverse proteins, showcasing its utility. This review aims to clarify the fundamental principles of His-HDX-MS, detail the experimental workflow, explain data analysis procedures and provide guidance for interpreting the obtained results.

Beyond the Plateau: pL Dependence of Proton Inventories as a Tool for Studying Ribozyme and Ribonuclease Catalysis

Acid/base catalysis is an important catalytic strategy used by ribonucleases and ribozymes; however, understanding the number and identity of functional groups involved in proton transfer remains challenging. The proton inventory (PI) technique analyzes the dependence of the enzyme reaction rate on the ratio of D₂O to H₂O and can provide information about the number of exchangeable sites that produce isotope effects and their magnitude. The Gross-Butler (GB) equation is used to evaluate H/D fractionation factors from PI data typically collected under conditions (i.e., a "plateau" in the pH-rate profile) assuming minimal change in active site residue ionization. However, restricting PI analysis to these conditions is problematic for many ribonucleases, ribozymes, and their variants due to ambiguity in the roles of active site residues, the lack of a plateau within the accessible pL range, or cooperative interactions between active site functional groups undergoing ionization. Here, we extend the integration of species distributions for alternative enzyme states in noncooperative models of acid/base catalysis into the GB equation, first used by Bevilacqua and colleagues for the HDV ribozyme, to develop a general population-weighted GB equation that allows simulation and global fitting of the three-dimensional relationship of the D₂O ratio (n) versus pL versus k_n/k₀. Simulations using the GPW-GB equation of PI results for RNase A, HDVrz, and VSrz illustrate that data obtained at multiple selected pL values across the pL-rate profile can assist in the planning and interpreting of solvent isotope effect experiments to distinguish alternative mechanistic models.

Pac13 is a Small, Monomeric Dehydratase that Mediates the Formation of the 3'-Deoxy Nucleoside of Pacidamycins

The uridyl peptide antibiotics (UPAs), of which pacidamycin is a member, have a clinically unexploited mode of action and an unusual assembly. Perhaps the most striking feature of these molecules is the biosynthetically unique 3'-deoxyuridine that they share. This moiety is generated by an unusual, small and monomeric dehydratase, Pac13, which catalyses the dehydration of uridine-5'-aldehyde. Here we report the structural characterisation of Pac13 with a series of ligands, and gain insight into the enzyme's mechanism demonstrating that H42 is critical to the enzyme's activity and that the reaction is likely to proceed via an E1cB mechanism. The resemblance of the 3'-deoxy pacidamycin moiety with the synthetic anti-retrovirals, presents a potential opportunity for the utilisation of Pac13 in the biocatalytic generation of antiviral compounds.

Biomolecular NMR: Past and future

The editors of this special volume suggested this topic, presumably because of the perspective lent by our combined >90-year association with biomolecular NMR. What follows is our personal experience with the evolution of the field, which we hope will illustrate the trajectory of change over the years. As for the future, one can confidently predict that it will involve unexpected advances. Our narrative is colored by our experience in using the NMR Facility for Biomedical Studies at Carnegie-Mellon University (Pittsburgh) and in developing similar facilities at Purdue (1977-1984) and the University of Wisconsin-Madison (1984-). We have enjoyed developing NMR technology and making it available to collaborators and users of these facilities. Our group's association with the Biological Magnetic Resonance data Bank (BMRB) and with the Worldwide Protein Data Bank (wwPDB) has also been rewarding. Of course, many groups contributed to the early growth and development of biomolecular NMR, and our brief personal account certainly omits many important milestones.

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.

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.

Complexity of protein energy landscapes studied by solution NMR relaxation dispersion experiments

The millisecond time scale motions in ribonuclease A (RNase A) were studied by solution NMR CPMG and off-resonance R1ρ relaxation dispersion experiments over a wide pH and temperature range. These experiments identify three separate protein regions termed Cluster 1, Cluster 2, and R33, whose motions are governed by distinct thermodynamic parameters. Moreover, each of these regions has motions with different pH dependencies. Cluster 1 shows an increase in activation enthalpy and activation entropy as the pH is lowered, whereas Cluster 2 exhibits the opposite behavior. In contrast, the activation enthalpy and entropy of R33 show no pH dependence. Compounding the differences, Δω values for Cluster 2 are characteristic of two-site conformational exchange, yet similar analysis for Cluster 1 indicates that this region of the enzyme exhibits conformational fluctuations between a major conformer and a pH-dependent average of protonated and deprotonated minor conformers.

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.

Structural heterogeneity among four subunits in pyranose 2-oxidase: A molecular dynamics simulation study

The homotetramer pyranose 2-oxidase (P2O) from Tetrametes multicolor contains flavin adenine dinucleotide (FAD) as a cofactor, and displays two conformers with different transient fluorescence spectra and lifetimes (ca. 0.1 ps and 360 ps). The ultrashort lifetimes of isoalloxazine (Iso) are ascribed to the photoinduced electron transfer (ET) from Trp168 to the excited Iso. Here, the structural heterogeneity among the four subunits in solution was studied by means of molecular dynamics simulation (MDS). The ET donor–acceptor distances in crystal and solution were compared. The distribution of the H-bond distances between Iso and the surrounding amino acids revealed appreciable differences among the four subunits. The structural fluctuations in two distant places were examined for the Iso-P and Iso-Q distances (where P and Q are Trp or Tyr) with the correlation coefficients between Iso-P and Iso-Q distances, revealing cooperative motions even though P and Q were more than 1 nm apart and located in different subunits. Moreover, distributions of the distances between Iso and its closest ionic amino acids markedly differed among the four subunits. Electrostatic (ES) energies between the Iso anion and the ionic amino acids in the entire protein were obtained using a static dielectric constant of 1. The ES energy in each subunit was strongly influenced by the other subunits, whilst the distributions of the ES energies greatly differed among the four subunits. This heterogeneous distribution of the ES energy between subunits may contribute to the large differences in the experimentally detected ET rates.

Urea's effect on the ribonuclease A catalytic efficiency: a kinetic, 1H NMR and molecular orbital study

Understanding of protein-urea interactions is one of the greatest challenges to modern structural protein chemistry. Based in enzyme kinetics experiments and (1)H NMR spectroscopic analysis we proposed that urea, at low concentrations, directly interacts with the protonated histidines of the active center of RNase A, following a simple model of competitive inhibition. These results were supported by theoretical analysis based on the frontier molecular orbital theory and suggest that urea might establish a favorable interaction with the cationic amino acids. Our experimental evidence and theoretical analysis indicate that the initials steps of the molecular mechanism of Urea-RNase A interaction passes through the establishment of a three center four electron adduct. Also, our results would explain the observed disruption of the (1)H NMR signals corresponding to H12 and H119 (involved in catalysis) of the RNase A studied in the presence of urea. Our interaction model of urea-amino acids (cationic) can be extended to explain the inactivation of other enzymes with cationic amino acids at the active site.

Quantitative measurement of the solvent accessibility of histidine imidazole groups in proteins

We report a method for expressing the solvent accessibility of histidine imidazole groups in proteins. The method is based on measuring the rate of the hydrogen exchange (HX) reaction of the imidazole C(ε1)-hydrogen. The rate profile of the HX reaction as a function of pH gives a sigmoidal curve, which reaches the maximal rate constant (k(max)) on the alkaline side of the sigmoidal curve. To quantitatively describe the solvent accessibility of imidazole groups in proteins, it is necessary to compare the k(max) of the imidazole groups with their intrinsic k(max) ((i)k(max)), the maximal rate constants for the given imidazole groups when they are fully exposed to the bulk solvent. However, the mechanism of the HX reaction suggests that the (i)k(max) of an imidazole group differs depending on its pK(a), and no systematic study has been conducted to clarify how the (i)k(max) is affected by pK(a). We therefore investigated the relationship between (i)k(max) and pK(a) using four imidazole derivatives at three different temperatures. The experimentally determined pK(a)-specific (i)k(max) values allowed us to derive a general formula to estimate the (i)k(max) value of any given imidazole group exhibiting a specific pK(a) at a specific temperature. Using the formula, the protection factors (PF), the ratio of (i)k(max) to k(max), of five imidazole groups in dihydrofolate reductase were obtained and used to express the magnitude of their solvent accessibility. In this definition, the smaller the PF value, the higher the solvent accessibility, and a value of 1 indicates full exposure to the bulk solvent. The solvent accessibility expressed by the PF values agreed well with the solvent accessible surface areas obtained from the X-ray diffraction data.

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.

Bovine pancreatic ribonuclease: fifty years of the first enzymatic reaction mechanism

Fifty years ago, the group of Tony Mathias and Bob Rabin at University College London deduced the first mechanism for catalysis by an enzyme, ribonuclease [Findlay, D., Herries, D. G., Mathias, A. P., Rabin, B. R., and Ross, C. A. (1961) Nature 190, 781-784]. Here, we celebrate this historic accomplishment by surveying knowledge of enzymology and protein science at that time, facts that led to the formulation of the mechanism, criticisms and alternative mechanisms, data that supported the proposed mechanism, and some of the refinements that have since provided a more precise picture of catalysis of RNA cleavage by ribonucleases. The Mathias and Rabin mechanism has appeared in numerous textbooks, monographs, and reviews and continues to have a profound impact on biochemistry.

Alteration of hydrogen bonding in the vicinity of histidine 48 disrupts millisecond motions in RNase A

The motion of amino acid residues on the millisecond (ms) time scale is involved in the tight regulation of catalytic function in numerous enzyme systems. Using a combination of mutational, enzymological, and relaxation-compensated (15)N Carr-Purcell-Meiboom-Gill (CPMG) methods, we have previously established the conformational significance of the distant His48 residue and the neighboring loop 1 in RNase A function. These studies suggested that RNase A relies on an intricate network of hydrogen bonding interactions involved in propagating functionally relevant, long-range ms motions to the catalytic site of the enzyme. To further investigate the dynamic importance of this H-bonding network, this study focuses on the individual replacement of Thr17 and Thr82 with alanine, effectively altering the key H-bonding interactions that connect loop 1 and His48 to the rest of the protein. (15)N CPMG dispersion studies, nuclear magnetic resonance (NMR) chemical shift analysis, and NMR line shape analysis of point mutants T17A and T82A demonstrate that the evolutionarily conserved single H-bond linking His48 to Thr82 is essential for propagating ms motions from His48 to the active site of RNase A on the time scale of catalytic turnover, whereas the T17A mutation increases the off rate and conformational exchange motions in loop 1. Accumulating evidence from our mutational studies indicates that residues experiencing conformational exchange in RNase A can be grouped into two separate clusters displaying distinct dynamical features, which appear to be independently affected by mutation. Overall, this study illuminates how tightly controlled and finely tuned ms motions are in RNase A, suggesting that designed modulation of protein motions may be possible.

An introduction to NMR-based approaches for measuring protein dynamics

Proteins are inherently flexible at ambient temperature. At equilibrium, they are characterized by a set of conformations that undergo continuous exchange within a hierarchy of spatial and temporal scales ranging from nanometers to micrometers and femtoseconds to hours. Dynamic properties of proteins are essential for describing the structural bases of their biological functions including catalysis, binding, regulation and cellular structure. Nuclear magnetic resonance (NMR) spectroscopy represents a powerful technique for measuring these essential features of proteins. Here we provide an introduction to NMR-based approaches for studying protein dynamics, highlighting eight distinct methods with recent examples, contextualized within a common experimental and analytical framework. The selected methods are (1) Real-time NMR, (2) Exchange spectroscopy, (3) Lineshape analysis, (4) CPMG relaxation dispersion, (5) Rotating frame relaxation dispersion, (6) Nuclear spin relaxation, (7) Residual dipolar coupling, (8) Paramagnetic relaxation enhancement. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Molecular mechanism for the denaturation of proteins by urea

Understanding protein-solute interactions is one of the sizable challenges of protein chemistry; therefore, numerous experimental studies have attempted to explain the mechanism by which proteins unfold in aqueous urea solutions. On the basis of kinetic evidence at low urea concentrations, (1)H NMR spectroscopic analysis, and molecular orbital calculations, we propose a mechanistic model for the denaturation of RNase A in urea. Our results support a direct interaction between urea and protonated histidine as the initial step for protein inactivation followed by hydrogen bond formation with polar residues, and the breaking of hydrophobic collapse as the final steps for protein denaturation. With the proposed model, we can rationalize apparently conflicting results in the literature about the mechanism of protein denaturation with urea.

Optimizing pKa computation in proteins with pH adapted conformations

pK(A) in proteins are determined by electrostatic energy computations using a small number of optimized protein conformations derived from crystal structures. In these protein conformations hydrogen positions and geometries of salt bridges on the protein surface were determined self-consistently with the protonation pattern at three pHs (low, ambient, and high). Considering salt bridges at protein surfaces is most relevant, since they open at low and high pH. In the absence of these conformational changes, computed pK(A)(comp) of acidic (basic) groups in salt bridges underestimate (overestimate) experimental pK(A)(exp), dramatically. The pK(A)(comp) for 15 different proteins with 185 known pK(A)(exp) yield an RMSD of 1.12, comparable with two other methods. One of these methods is fully empirical with many adjustable parameters. The other is also based on electrostatic energy computations using many non-optimized side chain conformers but employs larger dielectric constants at short distances of charge pairs that diminish their electrostatic interactions. These empirical corrections that account implicitly for additional conformational flexibility were needed to describe the energetics of salt bridges appropriately. This is not needed in the present approach. The RMSD of the present approach improves if one considers only strongly shifted pK(A)(exp) in contrast to the other methods under these conditions. Our method allows interpreting pK(A)(comp) in terms of pH dependent hydrogen bonding pattern and salt bridge geometries. A web service is provided to perform pK(A) computations.

Structural basis for catalysis by onconase

Onconase (ONC) is a homolog of bovine pancreatic ribonuclease (RNase A) from the frog Rana pipiens. ONC displays antitumoral activity and is in advanced clinical trials for the treatment of cancer. Here, we report the first atomic structures of ONC-nucleic acid complexes: a T89N/E91A ONC-5'-AMP complex at 1.65 A resolution and a wild-type ONC-d(AUGA) complex at 1.90 A resolution. The latter structure and site-directed mutagenesis were used to reveal the atomic basis for substrate recognition and turnover by ONC. The residues in ONC that are proximal to the scissile phosphodiester bond (His10, Lys31, and His97) and uracil nucleobase (Thr35, Asp67, and Phe98) are conserved from RNase A and serve to generate a similar bell-shaped pH versus k(cat)/K(M) profile for RNA cleavage. Glu91 of ONC forms two hydrogen bonds with the guanine nucleobase in d(AUGA), and Thr89 is in close proximity to that nucleobase. Installing a neutral or cationic residue at position 91 or an asparagine residue at position 89 virtually eliminated the 10(2)-fold guanine:adenine preference of ONC. A variant that combined such substitutions, T89N/E91A ONC, actually preferred adenine over guanine. In contrast, installing an arginine residue at position 91 increased the guanine preference and afforded an ONC variant with the highest known k(cat)/K(M) value. These data indicate that ONC discriminates between guanine and adenine by using Coulombic interactions and a network of hydrogen bonds. The structure of the ONC-d(AUGA) complex was also used to probe other aspects of catalysis. For example, the T5R substitution, designed to create a favorable Coulombic interaction between ONC and a phosphoryl group in RNA, increased ribonucleolytic activity by twofold. No variant, however, was more toxic to human cancer cells than wild-type ONC. Together, these findings provide a cynosure for understanding catalysis of RNA cleavage in a system of high medicinal relevance.

Roles of N-terminal pyroglutamate in maintaining structural integrity and pKa values of catalytic histidine residues in bullfrog ribonuclease 3

Many proteins and bioactive peptides contain an N-terminal pyroglutamate residue (Pyr1). This residue reduces the susceptibility of the protein to aminopeptidases and often has important functional roles. The antitumor ribonuclease RC-RNase 3 (RNase 3) from oocytes of Rana catesbeiana (bullfrog) is one such protein. We have produced recombinant RNase 3 containing the N-terminal Pyr1 (pRNase 3) and found it to be indistinguishable from the native RNase 3 by mass spectrometry and a variety of other biochemical and immunological criteria. We demonstrated by NMR analysis that the Pyr1 of pRNase 3 forms hydrogen bonds with Lys9 and Ile96 and stabilizes the N-terminal alpha-helix in a rigid conformation. In contrast, the N-terminal alpha-helix becomes flexible and the pKa values of the catalytic residues His10 and His97 altered when Pyr1 formation is blocked by an extra methionine at the N terminus in the recombinant mqRNase 3. Thus, our results provide a mechanistic explanation on the essential role of Pyr1 in maintaining the structural integrity, especially at the N-terminal alpha-helix, and in providing the proper environment for the ionization of His10 and His97 residues for catalysis and cytotoxicity against HeLa cells.

Sub-Angstrom resolution enzyme X-ray structures: is seeing believing?

Recent technical advances in crystallographic analysis, particularly highly focused and high brilliance synchrotron beam lines, have significantly improved the resolutions that are attainable for many macromolecular crystal structures. The Protein Data Bank (http://www.rcsb.org/pdb/) contains an increasing number of atomic resolution structures, which are providing a wealth of structural information that was not previously visible in lower resolution electron density maps. Here, we review the importance of visualizing hydrogen atoms and multiple sidechain conformations or anisotropy, as well as substrate strain, at sub-Angstrom resolution. The additional structural features that are visible in the electron density maps as a result of atomic resolution data provide a better understanding of the catalytic mechanisms of cholesterol oxidase, ribonuclease A, beta-lactamase, serine proteases, triosephosphate isomerase and endoglucanase.

The structural and functional studies of His119 and His12 in RNase A via chemical modification

The histidyl residues of bovine pancreatic ribonuclease A (RNase A) play a crucial role in enzymatic activity. Diethylpyrocarbonate (DEPC) is a potent inhibitor of RNase A, and its precise sites of action on the imidazole rings of the four histidyl residues of RNase A are not clearly defined. We have used a multidisciplinary approach including enzyme assay, calculation of accessible surface area (ASA), isoelectric pH gradient technique, fluorescence investigations, circular dichroism spectroscopy, differential scanning calorimetry, and 1H NMR analysis to study the sites of DEPC interaction with the imidazole rings of the four histidyl residues. Our results demonstrate that among the histidyl residues of RNase A, His48 is not accessible to react with DEPC. However, the sequential carbethoxylation of the imidazole rings of His119, His105, and His12 occurs on the nitrogen atoms of Ndelta, Nepsilon, and Nepsilon, respectively. Carbethoxylation of His119 was followed by conversion of the A conformation to the B conformation in the active site. However, the carbethoxylation of His12 was accompanied by a second spatial rotation of the corresponding imidazole ring in the active site to adopt a new conformation. These conformation changes are accompanied by subsequent decrements in the thermal stability of the protein. Therefore, these findings reinforce the important structural roles of the spatial positions for His119 and His12 in the active site of RNase A.

Variability in the pKa of histidine side-chains correlates with burial within proteins

Acidic pKas of histidines buried within the protein interior are frequently rationalized on the contradictory basis of either polar interactions within the protein or the effects of a hydrophobic environment. To examine these relationships, we surveyed the buried surface area, depth of burial, polar interactions, and crystallographic temperature factors of histidines of known pKa. It has been found that buried environments of histidines do not always result in acidic pKas. Instead, the variability of histidine pKas increases for residues where the majority of the side-chain is buried. Because buried histidines are always found in mixed polar/apolar environments, multiple environmental contributions to pKa values must be considered. However, the quantitative relationships between heterogeneous environments and pKa values are not immediately apparent from the available data.

His ... Asp catalytic dyad of ribonuclease A: histidine pKa values in the wild-type, D121N, and D121A enzymes

Bovine pancreatic ribonuclease A (RNase A) has a conserved His ... Asp catalytic dyad in its active site. Structural analyses had indicated that Asp121 forms a hydrogen bond with His119, which serves as an acid during catalysis of RNA cleavage. The enzyme contains three other histidine residues including His12, which is also in the active site. Here, 1H-NMR spectra of wild-type RNase A and the D121N and D121A variants were analyzed thoroughly as a function of pH. The effect of replacing Asp121 on the microscopic pKa values of the histidine residues is modest: none change by more than 0.2 units. There is no evidence for the formation of a low-barrier hydrogen bond between His119 and either an aspartate or an asparagine residue at position 121. In the presence of the reaction product, uridine 3'-phosphate (3'-UMP), protonation of one active-site histidine residue favors protonation of the other. This finding is consistent with the phosphoryl group of 3'-UMP interacting more strongly with the two active-site histidine residues when both are protonated. Comparison of the titration curves of the unliganded enzyme with that obtained in the presence of different concentrations of 3'-UMP shows that a second molecule of 3'-UMP can bind to the enzyme. Together, the data indicate that the aspartate residue in the His ... Asp catalytic dyad of RNase A has a measurable but modest effect on the ionization of the adjacent histidine residue.

Site-directed mutagenesis of the catalytic residues of bovine pancreatic deoxyribonuclease I

Bovine pancreatic deoxyribonuclease I (DNase I) is a well characterised endonuclease which cleaves double-stranded DNA to yield 5' phosphorylated polynucleotides. Co-crystal structures of DNase I with two different oligonucleotides have revealed the presence of several residues (R9, E78, H134, D168, D212 and H252) close to the scissile phosphate. The roles that these amino acids play in the catalytic mechanism have been investigated using site-directed mutagenesis. The following variants were used: R9A, E78T, H134Q, D168S, D212S and H252Q. The kinetics of all six mutants with both DNA and a small chromophoric substrate, thymidine-3',5'-di(p-nitrophenyl)-phosphate, were studied. Only R9A and E78T showed any significant turnover of the two substrates. D168S, H134Q, D212S and H252Q showed vanishingly low activities towards DNA and no detectable activity with thymidine-3',5'-di(p-nitrophenyl)-phosphate. These results demonstrate that H134, D168, D212 and H252 play a critical role in the catalytic mechanism. It is suggested that H134 and H252 (which are hydrogen-bonded to E78 and D212, respectively) provided general acid and general base catalysis. DNase I also requires Mg2+ and E39 has been identified as a ligand for this metal ion. We propose that D168 serves as a ligand for a second Mg2+, and thus DNase I, uses a two metal-ion hydrolytic mechanism. Both magnesium ions are used to supply electrophilic catalysis. Role assignment is based on the mutagenesis results, structural information, homologies between DNase I from different species and a comparison with exonuclease III. However, it is still not feasible to unequivocally assign a particular catalytic role to each amino acid/metal ion.

Catalytic activity of bovine seminal ribonuclease is essential for its immunosuppressive and other biological activities

Bovine seminal ribonuclease (BS-RNase) is a homologue of RNase A with special biological properties, including potent immunosuppressive activity. A mutant BS-RNase was created in which His-119, the active-site residue that acts as a general acid during catalysis, was changed to an aspartic acid. H119D BS-RNase formed a dimer with quaternary structure similar to that of the wild-type enzyme but with values of kcat. and kcat./Km for the cleavage of UpA [uridylyl(3'-->5')adenosine] that were 4 x 10(3)-fold lower. The mutant protein also demonstrated dramatically decreased immunosuppressive, anti-tumour, aspermatogenic, and embryotoxic activities. The catalytic activity of BS-RNase is therefore necessary for its special biological properties.

A misfolded but active dimer of bovine seminal ribonuclease

Bovine seminal ribonuclease (BS-RNase) is an unusual homolog of RNase A. Isolated from bulls as a dimer, BS-RNase has special biological properties including antispermatogenic, antitumor and immunosuppressive activities. The structural bases for these properties are unknown. Four forms of BS-RNase were isolated after folding and air oxidation of the denatured and reduced protein produced in Escherichia coli: two dimers (M = M and M x I, where x signifies an active site composed of residues from both subunits) and two monomers (M and I). Considerable ribonuclease activity was generated by air oxidation of an equimolar mixture of two inactive mutant proteins ([H12D]BS-RNase and [H119D]BS-RNase) prepared by site-directed mutagenesis. This activity came from a dimer (M x I) with a composite active site. 1H-NMR spectroscopy revealed that this dimer contained one correctly folded subunit (M), and one incorrectly folded subunit (I). Form I, which is a poor catalyst, was activated by ribonuclease S-protein, suggesting that the C-terminal portion of I was not folded properly. Electrospray-ionization mass spectrometry and sulfhydryl group titration indicated that I contains a single oxidized sulfhydryl group, which cannot participate in a disulfide bond. These results show that quaternary structure in BS-RNase is attained by the initial formation of two monomers, M and I, which then combine with another M to form M = M and M x I, respectively. Adventitious oxidation can thus lead to the formation of a misfolded but active enzyme (M x I).

NMR study of the positions of His-12 and His-119 in the ribonuclease A-uridine vanadate complex

The binding of uridine vanadate to ribonuclease A has been investigated by one- and two-dimensional 1H NMR. The homonuclear Nuclear Overhauser and exchange spectroscopy spectrum of the uridine vanadate/RNase A complex exhibits cross peaks between both the C5H and C6H protons of uridine vanadate and the H epsilon 1 proton of His-12 of ribonuclease A. These cross peaks suggest that the H epsilon 1 proton of His-12 is in the vicinity of the uracil base of uridine vanadate, as observed in the crystallographic structure of the uridine vanadate/RNase A complex. However, no cross peaks are observed between the C5H and C6H protons of uridine vanadate and the H epsilon 1 proton of His-119 of ribonuclease A, although they were predicted based upon the distances calculated from coordinates of the crystallographic structure of the complex. These results suggest that there is a significant difference between the positioning of the His-119 side chain in the solution and in the crystallographic structures.

Comparative studies of the interaction of human and bovine platelet factor 4 with heparin using histidine NMR resonances as spectroscopic probes

The pKa values of His-38 and His-50 of the heparin-binding protein, bovine platelet factor 4, are 5.6 and 6.5, respectively, as determined by 1H NMR spectroscopy. The 1H NMR resonance of His-38 of bovine platelet factor 4 which exhibits the lower pKa value is perturbed upon heparin binding to a greater degree than the resonance of His-50. Human platelet factor 4 contains the homologous residues His-23 and His-35. The pKa values of the two histidine residues of human platelet factor 4 are 5.3 and 6.4. The 1H NMR resonance of the histidine of human platelet factor 4 exhibiting the lower pKa value also is perturbed upon heparin binding to a greater degree than the histidine resonance exhibiting the higher pKa, thereby suggesting comparable heparin-protein interactions in bovine and human platelet factor 4.

Design of peptide-acridine mimics of ribonuclease activity

A series of peptide-acridine conjugates was designed and synthesized, based on three features of the proposed catalytic mechanism of RNase A: 2'-proton abstraction by His-12, proton donation to the leaving 5'-oxygen by His-119, and stabilization of the pentacoordinated phosphorous transition state by Lys-41. The substrate binding capability of RNase A was mimicked by the intercalator, acridine. Lysine served as a linker between acridine and the catalytic tripeptide. Cleavage of target RNA was monitored by agarose gel electrophoresis and by gel-permeation chromatography. The carboxyl-amidated conjugates HGHK(Acr)-NH2, HPHK(Acr)-NH2, and GGHK(Acr)-NH2 (where Acr indicates 2-methyl-9-acridinemethylene) all had similar hydrolytic activity. The catalytic mechanism most likely involved only the abstraction of the 2'-proton and stabilization of the transition state. These RNase mimics utilized rRNA and single-stranded RNA but not double-stranded RNA and tRNA as substrates.

Histidine pKa shifts accompanying the inactivating Asp121----Asn substitution in a semisynthetic bovine pancreatic ribonuclease

A semisynthetic RNase, RNase-(1-118).(111-124), consisting of a noncovalent complex between residues 1-118 of RNase (obtained from the proteolytic digestion of RNase A), and a synthetic 14-residue peptide containing residues 111-124 of RNase, exhibits 98% of the enzymatic activity of bovine pancreatic ribonuclease A (EC 3.1.27.5). The replacement of aspartic acid-121 by asparagine in this semisynthetic RNase to form the "D121N" analog reduces kcat/Km to 2.7% of the value for RNase A. In the present work, 1H NMR spectroscopy has been used to probe the ionization states of His12, His105, and His119 in this catalytically defective semisynthetic RNase. A comparison of the observed resonances of D121N with those previously determined by others for RNase A enabled us to assign the C2 proton NMR resonances to individual residues; the assignment of His119 was confirmed by titrating D121N with the fully deuterated peptide, [Asn121]-RNase-(111-124). The observed pKa values of His12, His105, and His119 decrease 0.18, 0.16, and 0.02 pH unit, respectively, as a result of the D121N replacement. Values calculated by using a finite difference algorithm to solve the Poisson-Boltzmann equation (the DELPHI program, version 3.0) and a refined 2.0-A coordinate set for the crystal structure of D121N differ significantly for active site residues His12 (delta pKa = -0.58) and His119 (delta pKa = -0.55) but not for His105 (delta pKa = -0.10). The elmination of bound water from the calculations reduced, but did not reconcile, these discrepancies (His12, delta pKa = -0.36; His119, delta pKa = -0.41).

Protein interactions with surface-immobilized metal ions: structure-dependent variations in affinity and binding capacity with temperature and urea concentration

We have used equilibrium binding analyses to evaluate the influence of temperature and urea on the affinity of hen egg white lysozyme and bovine pancreatic ribonuclease A for surface-immobilized Cu(II) ions. Linear Scatchard plots suggested that these model proteins were interacting with immobilized metal ions via a single class of intermediate-affinity (Kd = 10-40 microM) binding sites. Alterations in temperature had little or no effect on the immobilized Cu(II) binding capacity of either protein. Temperature effects on the interaction affinity, however, were protein-dependent and varied considerably. The affinity of lysozyme for immobilized Cu(II) ions was significantly decreased with increased temperature (0 degree C-37 degrees C), yet the affinity of ribonuclease did not vary measurably over the same temperature range. The van 't Hoff plot (1n K vs 1/T) for lysozyme suggests a straight line relationship (single mechanism) with a delta H of approximately -5.5 kcal/mol. Urea effects also varied in a protein-dependent manner. A 10-fold reduction in the affinity of lysozyme for the immobilized Cu(II) was observed with the urea concentrations up to 3 M; yet urea had no effect on the affinity of ribonuclease for the immobilized metal ions. Although the interaction capacity of lysozyme with the immobilized Cu(II) ions was decreased by 50% in 3 M urea, ribonuclease interaction capacity was not diminished in urea. Thus, temperature- and urea-dependent alterations in protein-metal ion interactions were observed for lysozyme but not ribonuclease A. The complete, yet reversible, inhibition of lysozyme- and ribonuclease-metal ion interactions by carboxyethylation with low concentrations of diethylpyrocarbonate provided direct evidence of histidyl involvement. The differential response of these proteins to the effects of temperature and urea was, therefore, interpreted based on calculated solvent-accessibilities and surface distributions of His residues, individual His residue pKa values, and specific features of the protein surface structure in the immediate environment of the surface-exposed histidyl residues. Possible interaction mechanisms involved in protein recognition of macromolecular surface-immobilized metal ions are presented.

Identification of a new site of conformational heterogeneity in unfolded ribonuclease A

The results presented here indicate that there are two slowly exchanging conformational isomers in unfolded bovine pancreatic ribonuclease A (RNase A) in the vicinity of Lys-41. The conformational heterogeneity is not observed in the fully folded protein. Therefore, one of the isomers may correspond to one of the slow-folding forms of the protein observed when refolding is initiated. These results were obtained from a chemically modified form of the protein, CL(7-41) RNase A, that has a dinitrophenyl cross-link between the epsilon-amino groups of Lys-7 and Lys-41. Extensive physical studies have shown that the cross-link does not significantly perturb the structure or the folding pathways of the protein. Therefore, the results obtained from this modified form of the protein are relevant to intact RNase A. The one-dimensional (1D) NMR spectrum of heat-unfolded CL(7-41) RNase A reveals that the singlet resonance for the C3H ring proton of the dinitrophenyl cross-link has been split into two unequal peaks in a 3:1 ratio, indicating that there are two distinct environments for the dinitrophenyl group. Variations in temperature, and the addition of urea, do not affect the relative peak intensities. The two peaks collapse into one after the protein is refolded. The observed splitting must originate from a slow reversible isomerization (greater than 100 msec) in a neighboring bond. The two most likely candidates are either the cis/trans isomerization of the Lys-41-Pro-42 peptide bond or hindered rotation about the disulfide bond between Cys-40 and Cys-95.

Identification of two essential histidine residues of ribonuclease T2 from Aspergillus oryzae

Ribonuclease (RNase) T2 from Aspergillus oryzae was modified by diethyl pyrocarbonate and iodoacetic acid. RNase T2 was rapidly inactivated by diethyl pyrocarbonate above pH 6.0 and by incorporation of a carboxymethyl group. No inactivation occurred in the presence of 3'AMP. 1H-NMR titration and photo-chemically induced dynamic nuclear polarization experiments demonstrated that two histidine residues were involved in the active site of RNase T2. Furthermore, analysis of inactive carboxymethylated RNase T2 showed that both His53 and His115 were partially modified to yield a total of one mole of N tau-carboxymethylhistidine/mole enzyme. The results indicate that the two histidine residues in the active site of RNase T2 are essential for catalysis and that modification of either His53 or His115 inactivates the enzyme.

15N- and 1H-NMR investigations of the active-site amino acids in semisynthetic RNase S' and RNase A

Extensive 15N-NMR investigations of active-site amino acids were made possible by the solid-phase synthesis of the N-terminal pentadecapeptide of RNase A with selectively 15N-enriched amino acids. On complexation with S-protein a fully active RNase S' complex was obtained. The 15N resonances of the side chains of lysine-7 (N epsilon), glutamine-11 (N gamma), and histidine-12 (N pi, tau) were studied in the free synthetic peptide, in the RNase S' complex and in the nucleotide complexes RNase S' with 2'CMP, 3'CMP, and 5'AMP. The analysis of the 15N-1H couplings, the 15N line broadenings due to proton exchange, and the chemical shift values showed that, while the imidazole ring is directly involved in the peptide-protein interaction, the side chains of Lys-7 and Gln-11 do not contribute to this interaction. In the nucleotide complexes the resonances of His-12 and Gln-11 are shifted downfield. In the 2'CMP complex a doublet for the N tau signal of His-12 indicates a stable H bond between this nitrogen and the phosphate group of nucleotide. The other nucleotide influence the resonances of the imidazole group much less, possibly due to a slightly different orientation of the phosphate group. The downfield shift of the Gln-11 resonance indicates an interaction between the carbonyl oxygen of the amide group and the phosphate moiety of the nucleotide. The only observable effect of nucleotide complexation on the Lys-7 signal is line broadening due to reduced proton exchange. For comparison with the 15N-NMR titration curves of His-12 in RNase S' the 1H-NMR titration curves of RNase A were also recorded. Both shape and pK values were very similar for the 15N and the 1H titration curves. An extensive analysis of the protonation equilibria with several fitting models showed that a mutual interaction of the imidazole groups of the active-site histidines results in flat titration curves. The Hill plots of all resonances of the imidazole rings, including the 15N resonances, show a small inflection in the pH range 5.8-6.4. Since the existence of a diimidazole system is most likely in this pH range, the inflection could be interpreted as a disturbance of the mutual electrostatic interaction of the active-site histidines by a partial H-bond formation between the imidazole groups.