Nuclear magnetic resonance methods for the detection and study of low-barrier hydrogen bonds on enzymes

Ensemble-function relationships to dissect mechanisms of enzyme catalysis

Decades of structure-function studies have established our current extensive understanding of enzymes. However, traditional structural models are snapshots of broader conformational ensembles of interchanging states. We demonstrate the need for conformational ensembles to understand function, using the enzyme ketosteroid isomerase (KSI) as an example. Comparison of prior KSI cryogenic x-ray structures suggested deleterious mutational effects from a misaligned oxyanion hole catalytic residue. However, ensemble information from room-temperature x-ray crystallography, combined with functional studies, excluded this model. Ensemble-function analyses can deconvolute effects from altering the probability of occupying a state (P-effects) and changing the reactivity of each state (k-effects); our ensemble-function analyses revealed functional effects arising from weakened oxyanion hole hydrogen bonding and substrate repositioning within the active site. Ensemble-function studies will have an integral role in understanding enzymes and in meeting the future goals of a predictive understanding of enzyme catalysis and engineering new enzymes.

Proton Bridging in Catalysis by and Inhibition of Serine Proteases of the Blood Cascade System

Inquiries into the participation of short hydrogen bonds in stabilizing transition states and intermediate states in the thrombin, factor Xa, plasmin and activated protein C-catalyzed reactions revealed that specific binding of effectors at Sn, n = 1-4 and S'n, n = 1-3 and at remote exosites elicit complex patterns of hydrogen bonding and involve water networks. The methods employed that yielded these discoveries include; (1) kinetics, especially partial or full kinetic deuterium solvent isotope effects with short cognate substrates and also with the natural substrates, (2) kinetic and structural probes, particularly low-field high-resolution nuclear magnetic resonance (1H NMR), of mechanism-based inhibitors and substrate-mimic peptide inhibitors. Short hydrogen bonds form at the transition states of the catalytic reactions at the active site of the enzymes as they do with mechanism-based covalent inhibitors of thrombin. The emergence of short hydrogen bonds at the binding interface of effectors and thrombin at remote exosites has recently gained recognition. Herein, I describe our contribution, a confirmation of this discovery, by low-field 1H NMR. The principal conclusion of this review is that proton sharing at distances below the sum of van der Waals radii of the hydrogen and both donor and acceptor atoms contribute to the remarkable catalytic prowess of serine proteases of the blood clotting system and other enzymes that employ acid-base catalysis. Proton bridges also play a role in tight binding in proteins and at exosites, i.e., allosteric sites, of enzymes.

Hydrogen Bonds: Simple after All?

Hydrogen bonds play integral roles in biological structure, function, and conformational dynamics and are fundamental to life as it has evolved on Earth. However, our understanding of these fundamental and ubiquitous interactions has seemed fractured and incomplete, and it has been difficult to extract generalities and principles about hydrogen bonds despite thousands of papers published on this topic, perhaps in part because of the expanse of this subject and the density of studies. Fortunately, recent hydrogen bond proposals, discussions, and debates have stimulated new tests and models and have led to a remarkably simple picture of the structure of hydrogen bonds. This knowledge also provides clarity concerning hydrogen bond energetics, limiting and simplifying the factors that need be considered. Herein we recount the advances that have led to this simpler view of hydrogen bond structure, dynamics, and energetics. A quantitative predictive model for hydrogen bond length can now be broadly and deeply applied to evaluate current proposals and to uncover structural features of proteins, their conformational restraints, and their correlated motions. In contrast, a quantitative energetic description of molecular recognition and catalysis by proteins remains an important ongoing challenge, although our improved understanding of hydrogen bonds may aid in testing predictions from current and future models. We close by codifying our current state of understanding into five "Rules for Hydrogen Bonding" that may provide a foundation for understanding and teaching about these vital interactions and for building toward a deeper understanding of hydrogen bond energetics.

Mapping the Hydrogen Bond Networks in the Catalytic Subunit of Protein Kinase A Using H/D Fractionation Factors

Protein kinase A is a prototypical phosphoryl transferase, sharing its catalytic core (PKA-C) with the entire kinase family. PKA-C substrate recognition, active site organization, and product release depend on the enzyme's conformational transitions from the open to the closed state, which regulate its allosteric cooperativity. Here, we used equilibrium nuclear magnetic resonance hydrogen/deuterium (H/D) fractionation factors (φ) to probe the changes in the strength of hydrogen bonds within the kinase upon binding the nucleotide and a pseudosubstrate peptide (PKI5-24). We found that the φ values decrease upon binding both ligands, suggesting that the overall hydrogen bond networks in both the small and large lobes of PKA-C become stronger. However, we observed several important exceptions, with residues displaying higher φ values upon ligand binding. Notably, the changes in φ values are not localized near the ligand binding pockets; rather, they are radiated throughout the entire enzyme. We conclude that, upon ligand and pseudosubstrate binding, the hydrogen bond networks undergo extensive reorganization, revealing that the open-to-closed transitions require global rearrangements of the internal forces that stabilize the enzyme's fold.

Determination of Hydrogen Bond Structure in Water versus Aprotic Environments To Test the Relationship Between Length and Stability

Hydrogen bonds profoundly influence the architecture and activity of biological macromolecules. Deep appreciation of hydrogen bond contributions to biomolecular function thus requires a detailed understanding of hydrogen bond structure and energetics and the relationship between these properties. Hydrogen bond formation energies (ΔGf) are enormously more favorable in aprotic solvents than in water, and two classes of contributing factors have been proposed to explain this energetic difference, focusing respectively on the isolated and hydrogen-bonded species: (I) water stabilizes the dissociated donor and acceptor groups much better than aprotic solvents, thereby reducing the driving force for hydrogen bond formation; and (II) water lengthens hydrogen bonds compared to aprotic environments, thereby decreasing the potential energy within the hydrogen bond. Each model has been proposed to provide a dominant contribution to ΔGf, but incisive tests that distinguish the importance of these contributions are lacking. Here we directly test the structural basis of model II. Neutron crystallography, NMR spectroscopy, and quantum mechanical calculations demonstrate that O-H···O hydrogen bonds in crystals, chloroform, acetone, and water have nearly identical lengths and very similar potential energy surfaces despite ΔGf differences >8 kcal/mol across these solvents. These results rule out a substantial contribution from solvent-dependent differences in hydrogen bond structure and potential energy after association (model II) and thus support the conclusion that differences in hydrogen bond ΔGf are predominantly determined by solvent interactions with the dissociated groups (model I). These findings advance our understanding of universal hydrogen-bonding interactions and have important implications for biology and engineering.

Proton bridging in the interactions of thrombin with hirudin and its mimics

Thrombin is the pivotal serine protease enzyme in the blood cascade system and thus a target of drug design for control of its activity. The most efficient nonphysiologic inhibitor of thrombin is hirudin, a naturally occurring small protein. Hirudin and its synthetic mimics employ a range of hydrogen bonding, salt bridging, and hydrophobic interactions with thrombin to achieve tight binding with K(i) values in the nano- to femtomolar range. The one-dimensional (1)H nuclear magnetic resonance spectrum recorded at 600 MHz reveals a resonance 15.33 ppm downfield from silanes in complexes between human α-thrombin and r-hirudin in pH 5.6-8.8 buffers and between 5 and 35 °C. There is also a resonance between 15.17 and 15.54 ppm seen in complexes of human α-thrombin with hirunorm IV, hirunorm V, an Nα(Me)Arg peptide, RGD-hirudin, and Nα-2-naphthylsulfonyl-glycyl-DL-4-amidinophenylalanyl-piperidide acetate salt (NAPAP), while there is no such low-field resonance observed in a complex of porcine trypsin and NAPAP. The chemical shifts suggest that these resonances represent H-bonded environments. H-Donor-acceptor distances in the corresponding H-bonds are estimated to be <2.7 Å. Addition of Phe-Pro-Arg-chloromethylketone (PPACK) to a complex of human α-thrombin with r-hirudin results in an additional signal at 18.03 ppm, which is 0.10 ppm upfield from the observed signal [Kovach, I. M., et al. (2009) Biochemistry 48, 7296-7304] for thrombin covalently modified with PPACK. In contrast, the peak at 15.33 ppm remains unchanged. The fractionation factors for the thrombin-hirudin complexes are near 1.0 within 20% error. The most likely site of the short H-bond in complexes of thrombin with the hirudin family of inhibitors is in the hydrophobic patch of the C-terminus of hirudin where Glu(57') and Glu(58') are embedded and interact with Arg(75) and Arg(77) and their solvate water (on thrombin). Glu(57') and Glu(58') present in the hirudin family of inhibitors make up a key binding epitope of fibrinogen, thrombin's prime substrate, which lends substantial interest to the short hydrogen bond as a binding element at the fibrinogen recognition site.

Enzyme active site interactions by Raman/FTIR, NMR, and ab initio calculations

Characterization of enzyme active site structure and interactions at high resolution is important for the understanding of the enzyme catalysis. Vibrational frequency and NMR chemical shift measurements of enzyme-bound ligands are often used for such purpose when X-ray structures are not available or when higher resolution active site structures are desired. This review is focused on how ab initio calculations may be integrated with vibrational and NMR chemical shift measurements to quantitatively determine high-resolution ligand structures (up to 0.001 Å for bond length and 0.01 Å for hydrogen bonding distance) and how interaction energies between bound ligand and its surroundings at the active site may be determined. Quantitative characterization of substrate ionic states, bond polarizations, tautomeric forms, conformational changes and its interactions with surroundings in enzyme complexes that mimic ground state or transition state can provide snapshots for visualizing the substrate structural evolution along enzyme-catalyzed reaction pathway. Our results have shown that the integration of spectroscopic studies with theoretical computation greatly enhances our ability to interpret experimental data and significantly increases the reliability of the theoretical analysis.

Role of glutamate 64 in the activation of the prodrug 5-fluorocytosine by yeast cytosine deaminase

Yeast cytosine deaminase (yCD) catalyzes the hydrolytic deamination of cytosine to uracil as well as the deamination of the prodrug 5-fluorocytosine (5FC) to the anticancer drug 5-fluorouracil. In this study, the role of Glu64 in the activation of the prodrug 5FC was investigated by site-directed mutagenesis, biochemical, nuclear magnetic resonance (NMR), and computational studies. Steady-state kinetics studies showed that the mutation of Glu64 causes a dramatic decrease in k(cat) and a dramatic increase in K(m), indicating Glu64 is important for both binding and catalysis in the activation of 5FC. (19)F NMR experiments showed that binding of the inhibitor 5-fluoro-1H-pyrimidin-2-one (5FPy) to the wild-type yCD causes an upfield shift, indicating that the bound inhibitor is in the hydrated form, mimicking the transition state or the tetrahedral intermediate in the activation of 5FC. However, binding of 5FPy to the E64A mutant enzyme causes a downfield shift, indicating that the bound 5FPy remains in an unhydrated form in the complex with the mutant enzyme. (1)H and (15)N NMR analysis revealed trans-hydrogen bond D/H isotope effects on the hydrogen of the amide of Glu64, indicating that the carboxylate of Glu64 forms two hydrogen bonds with the hydrated 5FPy. ONIOM calculations showed that the wild-type yCD complex with the hydrated form of the inhibitor 1H-pyrimidin-2-one is more stable than the initial binding complex, and in contrast, with the E64A mutant enzyme, the hydrated inhibitor is no longer favored and the conversion has a higher activation energy, as well. The hydrated inhibitor is stabilized in the wild-type yCD by two hydrogen bonds between it and the carboxylate of Glu64 as revealed by (1)H and (15)N NMR analysis. To explore the functional role of Glu64 in catalysis, we investigated the deamination of cytosine catalyzed by the E64A mutant by ONIOM calculations. The results showed that without the assistance of Glu64, both proton transfers before and after the formation of the tetrahedral reaction intermediate become partially rate-limiting steps. The results of the experimental and computational studies together indicate that Glu64 plays a critical role in both the binding and the chemical transformation in the conversion of the prodrug 5FC to the anticancer drug 5-fluorouracil.

Coupling of the distal hydrogen bond network to the exogenous ligand in substrate-bound, resting state human heme oxygenase

Mammalian heme oxygenase (HO) possesses catalytically implicated distal ordered water molecules within an extended H-bond network, with one of the ordered water molecules (#1) providing a bridge between the iron-coordinated ligand and the catalytically critical Asp140, that, in turn, serves as an acceptor for the Tyr58 OH H-bond. The degree of H-bonding by the ligated water molecule and the coupling of this water molecule to the H-bond network are of current interest and are herein investigated by (1)H NMR. Two-dimensional NMR allowed sufficient assignments to provide both the H-bond strength and hyperfine shifts, the latter of which were used to quantify the magnetic anisotropy in both the ferric high-spin aquo and low-spin hydroxo complexes. The anisotropy in the aquo complex indicates that the H-bond donation to water #1 is marginally stronger than in a bacterial HO, while the anisotropy for the hydroxo complex reveals a conventional (d(xz), d(yz))(1) ground state indicative of only moderate to weak H-bond acceptance by the ligated hydroxide. Mapping out the changes of the H-bond strengths in the network during the ligated water --> hydroxide conversion by correcting for the effects of magnetic anisotropy reveals a very substantial change in H-bond strength for Tyr58 OH and lesser effects on nearby H-bonds. The effect of pH on the H-bonding network in human HO is much larger and transmitted much further from the iron than in a pathogenic bacterial HO. The implications for the HO mechanism of the H-bond of Tyr58 to Asp140 are discussed.

The noncatalytic triad of alpha-amylases: a novel structural motif involved in conformational stability

Chloride-activated alpha-amylases contain a noncatalytic triad, independent of the glycosidic active site, perfectly mimicking the catalytic triad of serine-proteases and of other active serine hydrolytic enzymes. Mutagenesis of Glu, His, and Ser residues in various alpha-amylases shows that this pattern is a structural determinant of the enzyme conformation that cannot be altered without losing the intrinsic stability of the protein. (1)H-(15)N NMR spectra of a bacterial alpha-amylase reveal proton signals that are identical with the NMR signature of catalytic triads and especially a deshielded proton involving a protonated histidine and displaying properties similar to that of a low barrier hydrogen bond. It is proposed that the H-bond between His and Glu of the noncatalytic triad is an unusually strong interaction, responsible for the observed NMR signal and for the weak stability of the triad mutants. Furthermore, a stringent template-based search of the Protein Data Bank demonstrated that this motif is not restricted to alpha-amylases, but is also found in 80 structures from 33 different proteins, amongst which SH2 domain-containing proteins are the best representatives.

Solid-state NMR studies of weak interactions in supramolecular systems

The field of application of solid-state NMR to the study of supramolecular systems is growing rapidly, with many research groups involved in the development of techniques for the study of crystalline and amorphous phases. This Feature Article aims to provide an overview of the recent contributions of our research group to this field, paying particular attention to the study of the weak interactions such as hydrogen bonds in supramolecular systems through solid-state NMR investigations. The structure and dynamic behaviour of selected host-guest systems will be also discussed.

Ultrafast excited-state dynamics in the green fluorescent protein variant S65T/H148D. 1. Mutagenesis and structural studies

Wild type green fluorescent protein (wt-GFP) and the variant S65T/H148D each exhibit two absorption bands, A and B, which are associated with the protonated and deprotonated chromophores, respectively. Excitation of either band leads to green emission. In wt-GFP, excitation of band A ( approximately 395 nm) leads to green emission with a rise time of 10-15 ps, due to excited-state proton transfer (ESPT) from the chromophore hydroxyl group to an acceptor. This process produces an anionic excited-state intermediate I* that subsequently emits a green photon. In the variant S65T/H148D, the A band absorbance maximum is red-shifted to approximately 415 nm, and as detailed in the accompanying papers, when the A band is excited, green fluorescence appears with a rise time shorter than the instrument time resolution ( approximately 170 fs). On the basis of the steady-state spectroscopy and high-resolution crystal structures of several variants described herein, it is proposed that in S65T/H148D, the red shift of absorption band A and the ultrafast appearance of green fluorescence upon excitation of band A are due to a very short ( Collapse

Key Words Collapse

MESH Headings

• Green Fluorescent Proteins/chemistry
• Green Fluorescent Proteins/genetics
• Models, Molecular
• Mutagenesis
• Protein Conformation
• Spectrophotometry, Ultraviolet

Collapse

Grants

• R37 GM027738-22S1 NIGMS NIH HHS
• R01 GM042618-12S1 NIGMS NIH HHS
• R37 GM027738-20 NIGMS NIH HHS
• R01 GM042618-12 NIGMS NIH HHS
• R01 GM027738 NIGMS NIH HHS
• R01 GM042618-11 NIGMS NIH HHS
• R01 GM42618 NIGMS NIH HHS
• R01 GM042618 NIGMS NIH HHS
• R37 GM027738 NIGMS NIH HHS
• R37 GM027738-22 NIGMS NIH HHS
• R01 GM042618-08 NIGMS NIH HHS
• R37 GM027738-23 NIGMS NIH HHS
• R01 GM042618-14 NIGMS NIH HHS
• GM 27738 NIGMS NIH HHS
• R01 GM042618-13A1 NIGMS NIH HHS
• R01 GM042618-10 NIGMS NIH HHS
• R37 GM027738-21 NIGMS NIH HHS
• R01 GM042618-09 NIGMS NIH HHS

Collapse

Affiliation(s)

Xiaokun Shu Institute of Molecular Biology and Department of Physics, University of Oregon, Eugene, OR 97403-1229
Karen Kallio Institute of Molecular Biology and Department of Physics, University of Oregon, Eugene, OR 97403-1229
Xinghua Shi Department of Chemistry, Stanford University, Stanford, California 94305-5080
Paul Abbyad Department of Chemistry, Stanford University, Stanford, California 94305-5080
Pakorn Kanchanawong Department of Chemistry, Stanford University, Stanford, California 94305-5080
William Childs Department of Chemistry, Stanford University, Stanford, California 94305-5080
Steven G. Boxer Department of Chemistry, Stanford University, Stanford, California 94305-5080
S. James Remington Institute of Molecular Biology and Department of Physics, University of Oregon, Eugene, OR 97403-1229

Collapse

Crystal structures of the substrate free-enzyme, and reaction intermediate of the HAD superfamily member, haloacid dehalogenase DehIVa from Burkholderia cepacia MBA4

DehIVa is a haloacid dehalogenase (EC 3.8.1.2) from the soil and water borne bacterium Burkholderia cepacia MBA4, which belongs to the functionally variable haloacid dehalogenase (HAD) superfamily of enzymes. The haloacid dehalogenases catalyse the removal of halides from haloacids resulting in a hydroxlated product. These enzymes are of interest for their potential to degrade recalcitrant halogenated environmental pollutants and their use in the synthesis of industrial chemicals. The haloacid dehalogenases utilise a nucleophilic attack on the substrate by an aspartic acid residue to form an enzyme-substrate ester bond and concomitantly cleaving of the carbon-halide bond and release of a hydroxylated product following ester hydrolysis. We present the crystal structures of both the substrate-free DehIVa refined to 1.93 A resolution and DehIVa covalently bound to l-2-monochloropropanoate trapped as a reaction intermediate, refined to 2.7 A resolution. Electron density consistent with a previously unidentified yet anticipated water molecule in the active site poised to donate its hydroxyl group to the product and its proton to the catalytic Asp11 is evident. It has been unclear how substrate enters the active site of this and related enzymes. The results of normal mode analysis (NMA) are presented and suggest a means whereby the predicted global dynamics of the enzyme allow for entry of the substrate into the active site. In the context of these results, the possible role of Arg42 and Asn178 in a "lock down" mechanism affecting active site access is discussed. In silico substrate docking of enantiomeric substrates has been examined in order to evaluate the enzymes enantioselectivity.

Hydrogen bonding and dynamic behaviour in crystals and polymorphs of dicarboxylic-diamine adducts: a comparison between NMR parameters and X-ray diffraction studies

Fumaric, malonic, maleic, and hydromuconic (HOOCCH2(CH)2CH2COOH) acids were used to prepare a series of hydrogen-bonded adducts or salts, depending on whether acid-base proton transfer takes place, with the dibase [N(muCH2CH2)3N] in various stoichiometric ratios. The resulting compounds have been investigated by using the 1H MAS, 15N, and 13C cross polarisation magic-angle spinning (CPMAS) methods and discussed in relation to X-ray diffraction studies to ascertain the nature of the O-HO, NH-O, and N+-HO- hydrogen bonds between the various species. In addition, two polymorphic forms of the malonic compound and a hydrate in the maleic case were examined. We also present the correlations between the chemical shifts of the hydrogen-bonded protons and those from the proton transfer reaction (acid-to-base) with the heavy atom distances. The dynamic behaviour in the solid-state of the [N(muCH2CH2)3N] adducts with fumaric 2:1, maleic 1:1 hydrate, and hydromuconic acids, and a malonate 2:1 polymorph adduct have been investigated by using variable-temperature 1H spin-lattice relaxation times. A substantial agreement between the activation energies obtained from fitting the T1 data and the results of potential energy barrier calculations demonstrates that the facile reorientation of the [N(muCH2CH2)3N] molecule occurs in several of the adducts.

Artificial di-iron proteins: solution characterization of four helix bundles containing two distinct types of inter-helical loops

Peptide-based models have an enormous impact for the development of metalloprotein models, as they seem appropriate candidates to mimic both the structural characteristics and reactivity of the natural systems. Through the de novo design of four-helix bundles, we developed the DF (Due Ferri) family of artificial proteins, as models of di-iron and di-manganese metalloproteins. The goal of our research is to elucidate how the electrostatic environment, polarity and solvent accessibility of the metal-binding site, influence the functional properties of di-iron proteins. The first two subsets of the DF protein family, DF1 and DF2, consist of two non-covalently associated helix-loop-helix motifs, which bind the di-metal cofactor near the center of the structure. The DF2 subset was designed to improve the properties of DF1: DF2 and DF2t have several changes in their sequences to improve solubility and metal ion access, as well as a change in the loop connecting the two helices. In order to evaluate how these changes affect the overall structure of the model proteins, we solved the NMR structures of the di-Zn(II) complexes of DF2 and DF2t, and compared these structures with those recently obtained from X-ray crystallography. Further, we examined the thermodynamic consequences associated with the mutations, by measuring the stability of DF2t in the presence of different metal ions, and comparing the results with the data already obtained for DF2. Taken together, analysis of all the data showed the importance of the turn conformation in the design and stability of four-helix bundle.

Bisubstrate analog probes for the insulin receptor protein tyrosine kinase: molecular yardsticks for analyzing catalytic mechanism and inhibitor design

Bisubstrate analogs have the potential to provide enhanced specificity for protein kinase inhibition and tools to understand catalytic mechanism. Previous efforts led to the design of a peptide-ATP conjugate bisubstrate analog utilizing aminophenylalanine in place of tyrosine and a thioacetyl linker to the gamma-phosphate of ATP which was a potent inhibitor of the insulin receptor kinase (IRK). In this study, we have examined the contributions of various electrostatic and structural elements in the bisubstrate analog to IRK binding affinity. Three types of changes (seven specific analogs in all) were introduced: a Tyr isostere of the previous aminophenylalanine moiety, modifications of the spacer between the adenine and the peptide, and deletions and substitutions within the peptide moiety. These studies allowed a direct evaluation of the hydrogen bond strength between the anilino nitrogen of the bisubstrate analog and the enzyme catalytic base Asp and showed that it contributes 2.5 kcal/mol of binding energy, in good agreement with previous predictions. Modifications of the linker length resulted in weakened inhibitory affinity, consistent with the geometric requirements of an enzyme-catalyzed dissociative transition state. Alterations in the peptide motif generally led to diminished inhibitory potency, and only some of these effects could be rationalized based on prior kinetic and structural studies. Taken together, these results suggest that a combination of mechanism-based design and empirical synthetic manipulation will be necessary in producing optimized protein kinase bisubstrate analog inhibitors.

Observation of a Short, Strong Hydrogen Bond in the Active Site of Hydroxynitrile Lyase from Hevea brasiliensis Explains a Large pK Shift of the Catalytic Base Induced by the Reaction Intermediate

The hydroxynitrile lyase from Hevea brasiliensis (HbHNL) uses a catalytic triad consisting of Ser(80)-His(235)-Asp(207) to enhance the basicity of Ser(80)-O gamma for abstracting a proton from the OH group of the substrate cyanohydrin. Following the observation of a relatively short distance between a carboxyl oxygen of Asp(207) and the N delta(1)(His(235)) in a 1.1 A crystal structure of HbHNL, we here show by (1)H and (15)N-NMR spectroscopy that a short, strong hydrogen bond (SSHB) is formed between the two residues upon binding of the competitive inhibitor thiocyanate to HbHNL: the proton resonance of H-N delta 1(His(235)) moves from 15.41 ppm in the free enzyme to 19.35 ppm in the complex, the largest downfield shift observed so far upon inhibitor binding. Simultaneously, the D/H fractionation factor decreases from 0.98 to 0.35. In the observable pH range, i.e. between pH 4 and 10, no significant changes in chemical shifts (and therefore hydrogen bond strength) were observed for free HbHNL. For the complex with thiocyanate, the 19.35 ppm signal returned to 15.41 ppm at approximately pH 8, which indicates a pK(a) near this value for the H-N epsilon(2)(His(235)). These NMR results were analyzed on the basis of finite difference Poisson-Boltzmann calculations, which yielded the relative free energies of four protonation states of the His(235)-Asp(207) pair in solution as well as in the protein environment with and without bound inhibitor. The calculations explain all the NMR features, i.e. they suggest why a short, strong hydrogen bond is formed upon inhibitor binding and why this short, strong hydrogen bond reverts back to a normal one at approximately pH 8. Importantly, the computations also yield a shift of the free energy of the anionic state relative to the zwitterionic reference state by about 10.6 kcal/mol, equivalent to a shift in the apparent pK(a) of His(235) from 2.5 to 10. This huge inhibitor-induced increase in basicity is a prerequisite for His(235) to act as general base in the HbHNL-catalyzed cyanohydrin reaction.

1H NMR investigation of the solution structure of substrate-free human heme oxygenase: comparison to the cyanide-inhibited, substrate-bound complex

1H NMR was used to investigate the molecular structure, and dynamic properties of soluble, recombinant, substrate-free human heme oxygenase (apohHO) on a comparative basis with similar studies on the substrate complex. Limited but crucial sequence-specific assignments identify five conserved secondary structural elements, and the detection of highly characteristic dipolar or H-bond interactions among these elements together with insignificant chemical shift differences confirm a strongly conserved folding topology of helices C-H relative to that of substrate complexes in either solution or the crystal. The correction of the chemical shifts for paramagnetic and porphyrin ring current influences in the paramagnetic substrate complex reveals that the strength of all but one of the numerous relatively robust H-bonds are conserved in apohHO, and similar ordered water molecules are located near these H-bond donors as observed in the substrate complexes. The unique and significant weakening of the Tyr(58) OH hydrogen bond to the catalytically critical Asp(140) carboxylate in apohHO is suggested to arise from the removal of the axial H-bond acceptor ligand rather than the loss of substrate. The interhelical positions of the conserved strong H-bonds argue for a structural role in maintaining a conserved structure for helices C-H upon loss of substrate. While the structure and H-bond network are largely conserved upon loss of substrate, the variably increased rate of NH lability dictates a significant loss of dynamic stability in the conserved structure, particularly near the distal helix F.

A low-barrier hydrogen bond between histidine of secreted phospholipase A2 and a transition state analog inhibitor

This work describes in-depth NMR characterization of a unique low-barrier hydrogen bond (LBHB) between an active site residue from the enzyme and a bound inhibitor: the complex between secreted phospholipase A(2) (sPLA(2), from bee venom and bovine pancreas) and a transition-state analog inhibitor HK32. A downfield proton NMR resonance, at 17-18 ppm, was observed in the complex but not in the free enzyme. On the basis of site-specific mutagenesis and specific 15N-decoupling, this downfield resonance was assigned to the active site H48, which is part of the catalytic dyad D99-H48. These results led to a hypothesis that the downfield resonance represents the proton (H(epsilon 2) of H48) involved in the H-bonding between D99 and H48, in analogy with serine proteases. However, this was shown not to be the case by use of the bovine enzyme labeled with specific [15N(epsilon 2)]His. Instead, the downfield resonance arises from H(delta1) of H48, which forms a hydrogen bond with a non-bridging phosphonate oxygen of the inhibitor. Further studies showed that this proton displays a fractionation factor of 0.62(+/-0.06), and an exchange rate protection factor of >100 at 285 K and >40 at 298 K, which are characteristic of a LBHB. The pK(a) of the imidazole ring of H48 was shown to be shifted from 5.7 for the free enzyme to an apparent value of 9.0 in the presence of the inhibitor. These properties are very similar to those of the Asp em leader His LBHBs in serine proteases. Possible structural bases and functional consequences for the different locations of the LBHB between these two types of enzymes are discussed. The results also underscore the importance of using specific isotope labeling, rather than extrapolation of NMR results from other enzyme systems, to assign the downfield proton resonance to a specific hydrogen bond. Although our studies did not permit the strength of the LBHB to be accurately measured, the data do not provide support for an unusually strong hydrogen bond strength (i.e. >10 kcal/mol).

Investigation of the role of the histidine-aspartate pair in the human exonuclease III-like abasic endonuclease, Ape1

Hydrogen bonded histidine-aspartate (His-Asp) pairs are critical constituents in several key enzymatic reactions. To date, the role that these pairs play in catalysis is best understood in serine and trypsin-like proteases, where structural and biochemical NMR studies have revealed important pK(a) values and hydrogen bonding patterns within the catalytic pocket. However, the role of the His-Asp pair in metal-assisted catalysis is less clear. Here, we apply liquid-state NMR to investigate the role of a critical histidine residue of apurinic endonuclease 1 (Ape1), a human DNA repair enzyme that cleaves adjacent to abasic sites in DNA using one or more divalent cations and an active-site His-Asp pair. The results of these studies suggest that the Ape1 His-Asp pair does not function as either a general base catalyst or a metal ligand. Rather, the pair likely stabilizes the pentavalent transition state necessary for phospho-transfer.

NMR studies of the hydrogen bonds involving the catalytic triad ofEscherichia colithioesterase/protease I

Escherichia coli thioesterase/protease I (TEP-I) is a lipolytic enzyme of the serine protease superfamily with Ser(10), Asp(154) and His(157) as the catalytic triad residues. Based on comparison of the low-field (1)H nuclear magnetic resonance spectra of two mutants (S10G and S12G) and two transition state analogue complexes we have assigned the exchangeable proton resonances at 16.3 ppm, 14.3 ppm, and 12.8 ppm at pH 3.5 to His(157)-N(delta1)H, Ser(10)-O(gamma)H and His(157)-N(epsilon2)H, respectively. Thus, the presence of a strong Asp(154)-His(157) hydrogen bond in free TEP-I was observed. However, Ser(10)-O(gamma)H was shown to form a H-bond with a residue other than His(157)-N(epsilon2).

Five atomic resolution structures of endothiapepsin inhibitor complexes: implications for the aspartic proteinase mechanism

Endothiapepsin is derived from the fungus Endothia parasitica and is a member of the aspartic proteinase class of enzymes. This class of enzyme is comprised of two structurally similar lobes, each lobe contributing an aspartic acid residue to form a catalytic dyad that acts to cleave the substrate peptide bond. The three-dimensional structures of endothiapepsin bound to five transition state analogue inhibitors (H189, H256, CP-80,794, PD-129,541 and PD-130,328) have been solved at atomic resolution allowing full anisotropic modelling of each complex. The active sites of the five structures have been studied with a view to studying the catalytic mechanism of the aspartic proteinases by locating the active site protons by carboxyl bond length differences and electron density analysis. In the CP-80,794 structure there is excellent electron density for the hydrogen on the inhibitory statine hydroxyl group which forms a hydrogen bond with the inner oxygen of Asp32. The location of this proton has implications for the catalytic mechanism of the aspartic proteinases as it is consistent with the proposed mechanism in which Asp32 is the negatively charged aspartate. A number of short hydrogen bonds (approximately 2.6 A) with ESD values of around 0.01 A that may have a role in catalysis have been identified within the active site of each structure; the lengths of these bonds have been confirmed using NMR techniques. The possibility and implications of low barrier hydrogen bonds in the active site are considered.

The role of Tyr248 probed by mutant bovine carboxypeptidase A: insight into the catalytic mechanism of carboxypeptidase A

We have investigated the function of Tyr248 using bovine wild-type CPA and its Y248F and Y248A mutants to find that the K(M) values were increased by 4.5-11-fold and the k(cat) values were reduced by 4.5-10.7-fold by the replacement of Tyr248 with Phe for the hydrolysis of hippuryl-L-Phe (HPA) and N-[3-(2-furyl)acryloyl]-Phe-Phe (FAPP), respectively. In the case of O-(trans-p-chlorocinnamoyl)-L-beta-phenyllactate (ClCPL), an ester substrate, the K(M) value was increased by 2.5-fold, and the k(cat) was reduced by 20-fold. The replacement of Tyr248 with Ala decreased the k(cat) values by about 18- and 237-fold for HPA and ClCPL, respectively, demonstrating that the aromatic ring of Tyr248 plays a critical role in the enzymic reaction. The increases of the K(M) values were only 6- and 5-fold for HPA and ClCPL, respectively. Thus, the present study indicates clearly that Tyr248 plays an important role not only in the binding of substrate but also in the enzymic hydrolysis. The kinetic results may be rationalized by the proposition that the phenolic hydroxyl of Tyr248 forms a hydrogen bond with the zinc-bound water molecule, causing further activation of the water molecule by reducing its pK(a) value. The pH dependency study of k(cat) values and the solvent isotope effects also support the proposition. A unified catalytic mechanism is proposed that can account for the different kinetic behavior observed in the CPA-catalyzed hydrolysis of peptide and ester substrates.

Solution structure and dynamics of ribonuclease Sa

We have used NMR methods to characterize the structure and dynamics of ribonuclease Sa in solution. The solution structure of RNase Sa was obtained using the distance constraints provided by 2,276 NOEs and the C6-C96 disulfide bond. The 40 resulting structures are well determined; their mean pairwise RMSD is 0.76 A (backbone) and 1.26 A (heavy atoms). The solution structures are similar to previously determined crystal structures, especially in the secondary structure, but exhibit new features: the loop composed of Pro 45 to Ser 48 adopts distinct conformations and the rings of tyrosines 51, 52, and 55 have reduced flipping rates. Amide protons with greatly reduced exchange rates are found predominantly in interior beta-strands and the alpha-helix, but also in the external 3/10 helix and edge beta-strand linked by the disulfide bond. Analysis of (15)N relaxation experiments (R1, R2, and NOE) at 600 MHz revealed five segments, consisting of residues 1-5, 28-31, 46-50, 60-65, 74-77, retaining flexibility in solution. The change in conformation entropy for RNase SA folding is smaller than previously believed, since the native protein is more flexible in solution than in a crystal.

Substrate deformation in a hypoxanthine-guanine phosphoribosyltransferase ternary complex: the structural basis for catalysis

BACKGROUND

Hypoxanthine-guanine phosphoribosyltransferases (HGPRTs) are well-recognized antiparasitic drug targets. HGPRT is also a paradigmatic representative of the phosphoribosyltransferase family of enzymes, which includes other important biosynthetic and salvage enzymes and drug targets. To better understand the reaction mechanism of this enzyme, we have crystallized HGPRT from the apicomplexan protozoan Toxoplasma gondii as a ternary complex with a substrate and a substrate analog.

RESULTS

The crystal structure of T. gondii HGPRT with the substrate Mg2+-PRPP and a nonreactive substrate analog, 9-deazaguanine, bound in the active site has been determined at 1.05 A resolution and refined to a free R factor of 15.4%. This structure constitutes the first atomic-resolution structure of both a phosphoribosyltransferase and the central metabolic substrate PRPP. This pre-transition state complex provides a clearer understanding of the structural basis for catalysis by HGPRT.

CONCLUSIONS

Three types of substrate deformation, chief among them an unexpected C2'-endo pucker adopted by the PRPP ribose ring, raise the energy of the ground state. A cation-pi interaction between Tyr-118 and the developing oxocarbenium ion in the ribose ring helps to stabilize the transition state. Enforced substrate propinquity coupled with optimal reactive geometry for both the substrates and the active site residues with which they interact contributes to catalysis as well.

Low-barrier hydrogen bonds and enzymatic catalysis

Short, strong (low barrier) hydrogen bonds occur when the pK values of the atoms sharing the proton are similar. The overall distance is 2.5 A or less, the deuterium fractionation factor is less than 0.5, the proton NMR chemical shift can approach 20 ppm, and deuterium or tritium substitution causes an up-field change in the chemical shift. Such bonds can have deltaH values of 25 kcal/mol in the gas phase, and at least half that in water or other high-dielectric medium. The strength of the hydrogen bond in an active site drops by approximately 1 kcal/mol for each pH unit mismatch in pKs. When a weak hydrogen bond in the initial enzyme-substrate complex is converted into a low-barrier one by alteration of the pK of the substrate or catalytic group so that the pKs match, the increase in hydrogen bond strength can be used to help catalyze the reaction. A well-established example of this is the reaction catalyzed by serine proteases. The pK of neutral histidine is 14, while that of aspartate is approximately 6. Proton transfer from serine to permit attack on bound substrate produces protonated histidine, with a pK now matching that of aspartate. Studies with trifluoromethyl ketone inhibitors that form tetrahedral adducts show up to five orders of magnitude in binding strength as the result of formation of a low-barrier hydrogen bond between aspartate and histidine. Other enzymes whose mechanisms appear to involve low-barrier hydrogen bonds include liver alcohol dehydrogenase, steroid isomerase, triose-P isomerase, aconitase, citrate synthase, and zinc proteases. It is likely that low-barrier hydrogen bonds form at the transition state of any reaction involving general-acid or general-base catalysis, as at that point the pKs of the catalytic group and reactant will be equal.

The deuterium isotope effect on the NMR signal of the low-barrier hydrogen bond in a transition-state analog complex of chymotrypsin

The participation of a low-barrier hydrogen bond (LBHB) in the mechanism of action of chymotrypsin introduces a new role for Asp 102 and His 57 in catalysis [C. S. Cassidy, J. Lin, and P. A. Frey (1997) Biochemistry 36, 4576-4584]. It is postulated that the LBHB increases the basicity of His 57-N(epsilon2) in the transition state, thereby facilitating the abstraction of a proton from Ser 195, and stabilizes the tetrahedral intermediate in the acylation step. Evidence for this mechanism includes the downfield chemical shift of the proton bridging His 57 and Asp 102 in transition-state analog complexes and the low deuterium fractionation factors for this proton in the same complexes. We present additional spectroscopic evidence supporting the assignment of an LBHB between His 57 and Asp 102. The tetrahedral addition complex between Ser 195 of chymotrypsin and N-acetyl-l-leucyl-l-phenylalanyl trifluoromethylketone is regarded as a close structural analog of a tetrahedral intermediate. The deuterium NMR signal for the downfield deuteron bridging His 57 and Asp 102 in D(2)O has now been observed as a broad band centered at 17.8 +/- 0.5 ppm. The proton NMR signal in H(2)O is centered at 18.9 +/- 0.05 ppm. The two signals are clearly separated corresponding to a deuterium isotope effect of Delta[delta(H) - delta(D)] = 1.1 +/- 0.5 ppm. Deuterium isotope effects in this range are characteristic of LBHBs, and this observation provides further support for the assignment of the proton bridging His 57 and Asp 102 in transition-state analog complexes as an LBHB.