Minimal transition state charge stabilization of the oxyanion during peptide bond formation by the ribosome

Peptide bond formation during ribosomal protein synthesis involves an aminolysis reaction between the aminoacyl α-amino group and the carbonyl ester of the growing peptide via a transition state with a developing negative charge, the oxyanion. Structural and molecular dynamic studies have suggested that the ribosome may stabilize the oxyanion in the transition state of peptide bond formation via a highly ordered water molecule. To biochemically investigate this mechanistic hypothesis, we estimated the energetic contribution to catalytic charge stabilization of the oxyanion using a series of transition state mimics that contain different charge distributions and hydrogen bond potential on the functional group mimicking the oxyanion. Inhibitors containing an oxyanion mimic that carried a neutral charge and a mimic that preserved the negative charge but could not form hydrogen bonds had less than a 3-fold effect on inhibitor binding affinity. These observations argue that the ribosome provides minimal transition state charge stabilization to the oxyanion during peptide bond formation via the water molecule. This is in contrast to the substantial level of oxyanion stabilization provided by serine proteases. This suggests that the oxyanion may be neutralized via a proton shuttle, resulting in an uncharged transition state.     

Catalytic site forms and controls in ATP synthase catalysis

A suggested minimal scheme for substrate binding by and interconversion of three forms of the catalytic sites of the ATP synthase is presented. Each binding change, that drives simultaneous interchange of the three catalytic site forms, requires a 120 degrees rotation of the gamma with respect to the beta subunits. The binding of substrate(s) at two catalytic sites is regarded as sufficing for near maximal catalytic rates to be attained. Although three sites do not need to be filled for rapid catalysis, during rapid bisite catalysis some enzyme may be transiently present with three sites filled. Forms with preferential binding for ADP and P(i) or for ATP are considered to arise from the transition state and participate in other steps of the catalysis. Intermediate forms and steps that may be involved are evaluated. Experimental evidence for energy-dependent steps and for control of coupling to proton translocation and transition state forms are reviewed. Impact of relevant past data on present understanding of catalytic events is considered. In synthesis a key step is suggested in which proton translocation begins to deform an open site so as to increase the affinity for ADP and P(i), that then bind and pass through the transition state, and yield tightly bound ATP in one binding change. ADP binding appears to be a key parameter controlling rotation during synthesis. In hydrolysis ATP binding to a loose site likely precedes any proton translocation, with proton movement occurring as the tight site form develops. Aspects needing further study are noted. Characteristics of the related MgADP inhibition of the F(1) ATPases that have undermined many observations are summarized, and relations of three-site filling to catalysis are assessed.     

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

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

Short-strong hydrogen bonds and a low barrier transition state for the proton transfer reaction in RNase A catalysis: a quantum chemical study

There is growing evidence that some enzymes catalyze reactions through the formation of short-strong hydrogen bonds as first suggested by Gerlt and Gassman. Support comes from several experimental and quantum chemical studies that include correlation energies on model systems. In the present study, the process of proton transfer between hydroxyl and imidazole groups, a model of the crucial step in the hydrolysis of RNA by the enzymes of the RNase A family, is investigated at the quantum mechanical level of density functional theory and perturbation theory at the MP2 level. The model focuses on the nature of the formation of a complex between the important residues of the protein and the hydroxyl group of the substrate. We have also investigated different configurations of the ground state that are important in the proton transfer reaction. The nature of bonding between the catalytic unit of the enzyme and the substrate in the model is investigated by Bader's atoms in molecule theory. The contributions of solvation and vibrational energies corresponding to the reactant, the transition state and the product configurations are also evaluated. Furthermore, the effect of protein environment is investigated by considering the catalytic unit surrounded by complete proteins--RNase A and Angiogenin. The results, in general, indicate the formation of a short-strong hydrogen bond and the formation of a low barrier transition state for the proton transfer model of the enzyme.     

The catalytic aspartate is protonated in the Michaelis complex formed between trypsin and an in vitro evolved substrate-like inhibitor: a refined mechanism of serine protease action

The mechanism of serine proteases prominently illustrates how charged amino acid residues and proton transfer events facilitate enzyme catalysis. Here we present an ultrahigh resolution (0.93 Å) x-ray structure of a complex formed between trypsin and a canonical inhibitor acting through a substrate-like mechanism. The electron density indicates the protonation state of all catalytic residues where the catalytic histidine is, as expected, in its neutral state prior to the acylation step by the catalytic serine. The carboxyl group of the catalytic aspartate displays an asymmetric electron density so that the O_δ2–C_γ bond appears to be a double bond, with O_δ2 involved in a hydrogen bond to His-57 and Ser-214. Only when Asp-102 is protonated on O_δ1 atom could a density functional theory simulation reproduce the observed electron density. The presence of a putative hydrogen atom is also confirmed by a residual mF_obs − DF_calc density above 2.5 σ next to O_δ1. As a possible functional role for the neutral aspartate in the active site, we propose that in the substrate-bound form, the neutral aspartate residue helps to keep the pK_a of the histidine sufficiently low, in the active neutral form. When the histidine receives a proton during the catalytic cycle, the aspartate becomes simultaneously negatively charged, providing additional stabilization for the protonated histidine and indirectly to the tetrahedral intermediate. This novel proposal unifies the seemingly conflicting experimental observations, which were previously seen as either supporting the charge relay mechanism or the neutral pK_a histidine theory.     

Catalytic roles for two water bridged residues (Asp-98 and His-255) in the active site of copper-containing nitrite reductase

Two active site residues, Asp-98 and His-255, of copper-containing nitrite reductase (NIR) from Alcaligenes faecalis have been mutated to probe the catalytic mechanism. Three mutations at these two sites (D98N, H255D, and H255N) result in large reductions in activity relative to native NIR, suggesting that both residues are involved intimately in the reaction mechanism. Crystal structures of these mutants have been determined using data collected to better than 1. 9-A resolution. In the native structure, His-255 Nepsilon2 forms a hydrogen bond through a bridging water molecule to the side chain of Asp-98, which also forms a hydrogen bond to a water or nitrite oxygen ligated to the active site copper. In the D98N mutant, reorientation of the Asn-98 side chain results in the loss of the hydrogen bond to the copper ligand water, consistent with a negatively charged Asp-98 directing the binding and protonation of nitrite in the native enzyme. An additional solvent molecule is situated between residues 255 and the bridging water in the H255N and H255D mutants and likely inhibits nitrite binding. The interaction of His-255 with the bridging water appears to be necessary for catalysis and may donate a proton to reaction intermediates in addition to Asp-98.     

Reaction mechanism of the binuclear zinc enzyme glyoxalase II - A theoretical study

The glyoxalase system catalyzes the conversion of toxic methylglyoxal to nontoxic d-lactic acid using glutathione (GSH) as a coenzyme. Glyoxalase II (GlxII) is a binuclear Zn enzyme that catalyzes the second step of this conversion, namely the hydrolysis of S-d-lactoylglutathione, which is the product of the Glyoxalase I (GlxI) reaction. In this paper we use density functional theory method to investigate the reaction mechanism of GlxII. A model of the active site is constructed on the basis of the X-ray crystal structure of the native enzyme. Stationary points along the reaction pathway are optimized and the potential energy surface for the reaction is calculated. The calculations give strong support to the previously proposed mechanism. It is found that the bridging hydroxide is capable of performing nucleophilic attack at the substrate carbonyl to form a tetrahedral intermediate. This step is followed by a proton transfer from the bridging oxygen to Asp58 and finally C-S bond cleavage. The roles of the two zinc ions in the reaction mechanism are analyzed. Zn2 is found to stabilize the charge of tetrahedral intermediate thereby lowering the barrier for the nucleophilic attack, while Zn1 stabilizes the charge of the thiolate product, thereby facilitating the C-S bond cleavage. Finally, the energies involved in the product release and active-site regeneration are estimated and a new possible mechanism is suggested.     

A computational method for design of connected catalytic networks in proteins

Computational design of new active sites has generally proceeded by geometrically defining interactions between the reaction transition state(s) and surrounding side‐chain functional groups which maximize transition‐state stabilization, and then searching for sites in protein scaffolds where the specified side‐chain–transition‐state interactions can be realized. A limitation of this approach is that the interactions between the side chains themselves are not constrained. An extensive connected hydrogen bond network involving the catalytic residues was observed in a designed retroaldolase following directed evolution. Such connected networks could increase catalytic activity by preorganizing active site residues in catalytically competent orientations, and enabling concerted interactions between side chains during catalysis, for example, proton shuffling. We developed a method for designing active sites in which the catalytic side chains, in addition to making interactions with the transition state, are also involved in extensive hydrogen bond networks. Because of the added constraint of hydrogen‐bond connectivity between the catalytic side chains, to find solutions, a wider range of interactions between these side chains and the transition state must be considered. Our new method starts from a ChemDraw‐like two‐dimensional representation of the transition state with hydrogen‐bond donors, acceptors, and covalent interaction sites indicated, and all placements of side‐chain functional groups that make the indicated interactions with the transition state, and are fully connected in a single hydrogen‐bond network are systematically enumerated. The RosettaMatch method can then be used to identify realizations of these fully‐connected active sites in protein scaffolds. The method generates many fully‐connected active site solutions for a set of model reactions that are promising starting points for the design of fully‐preorganized enzyme catalysts.     

Testing electrostatic complementarity in enzyme catalysis: hydrogen bonding in the ketosteroid isomerase oxyanion hole

A longstanding proposal in enzymology is that enzymes are electrostatically and geometrically complementary to the transition states of the reactions they catalyze and that this complementarity contributes to catalysis. Experimental evaluation of this contribution, however, has been difficult. We have systematically dissected the potential contribution to catalysis from electrostatic complementarity in ketosteroid isomerase. Phenolates, analogs of the transition state and reaction intermediate, bind and accept two hydrogen bonds in an active site oxyanion hole. The binding of substituted phenolates of constant molecular shape but increasing p K _a models the charge accumulation in the oxyanion hole during the enzymatic reaction. As charge localization increases, the NMR chemical shifts of protons involved in oxyanion hole hydrogen bonds increase by 0.50–0.76 ppm/p K _a unit, suggesting a bond shortening of ˜0.02 Å/p K _a unit. Nevertheless, there is little change in binding affinity across a series of substituted phenolates (ΔΔG = −0.2 kcal/mol/p K _a unit). The small effect of increased charge localization on affinity occurs despite the shortening of the hydrogen bonds and a large favorable change in binding enthalpy (ΔΔH = −2.0 kcal/mol/p K _a unit). This shallow dependence of binding affinity suggests that electrostatic complementarity in the oxyanion hole makes at most a modest contribution to catalysis of ˜300-fold. We propose that geometrical complementarity between the oxyanion hole hydrogen-bond donors and the transition state oxyanion provides a significant catalytic contribution, and suggest that KSI, like other enzymes, achieves its catalytic prowess through a combination of modest contributions from several mechanisms rather than from a single dominant contribution.     

Proton transport in carbonic anhydrase: Insights from molecular simulation

This article reviews the insights gained from molecular simulations of human carbonic anhydrase II (HCA II) utilizing non-reactive and reactive force fields. The simulations with a reactive force field explore protein transfer and transport via Grotthuss shuttling, while the non-reactive simulations probe the larger conformational dynamics that underpin the various contributions to the rate-limiting proton transfer event. Specific attention is given to the orientational stability of the His64 group and the characteristics of the active site water cluster, in an effort to determine both of their impact on the maximal catalytic rate. The explicit proton transfer and transport events are described by the multistate empirical valence bond (MS-EVB) method, as are alternative pathways for the excess proton charge defect to enter/leave the active site. The simulation results are interpreted in light of experimental results on the wild-type enzyme and various site-specific mutations of HCA II in order to better elucidate the key factors that contribute to its exceptional efficiency.     

Evaluation of catalytic free energies in genetically modified proteins

A combination of the empirical valence bond method and a free energy perturbation approach is used to simulate the activity of genetically modified enzymes. The simulations reproduce in a semiquantitative way the observed effects of mutations on the activity and binding free energies of trypsin and subtilisin. This suggests that we are approaching a stage of quantitative structure-function correlation of enzymes. The analysis of the calculations points towards the electrostatic energy of the reacting system as the key factor in enzyme catalysis. The changes in the charges of the reacting system and the corresponding changes in "solvation" free energy (generalized here as the interaction between the charges and the given microenvironment) are emphasized. It is argued that a reliable evaluation of these changes might be sufficient for correlating structure and catalysis. The use of free energy perturbation methods and thermodynamic cycles for evaluation of solvation energies and reactivity is discussed, pointing out our early contributions. The apparent elaborated nature of our treatment is clarified, explaining that such a treatment is essential for consistent calculations of chemical reactions in polar environments. The problems associated with seemingly more rigorous quantum mechanical methods are discussed, emphasizing the inconsistency associated with using gas phase charge distributions. The importance of dynamic aspects is examined by evaluating the autocorrelation of the protein "reaction field" on the reacting substrate. It is found that, at least in the present case, dynamic effects are not important. The nature of the catalytic free energy is considered, arguing that the protein provides preoriented dipoles (polarized to stabilize the transition state charge distribution) and small reorganization energy, thus reducing the activation free energy. The corresponding catalytic free energy is related to the folding free energy, which is being invested in aligning the active site dipoles.     

Acetyl-CoA enolization in citrate synthase: A quantum mechanical/molecular mechanical (QM/MM) study

Citrate synthase forms citrate by deprotonation of acetyl-CoA followed by nucleophilic attack of this substrate on oxaloacetate, and subsequent hydrolysis. The rapid reaction rate is puzzling because of the instability of the postulated nucleophilic intermediate, the enolate of acetyl-CoA. As alternatives, the enol of acetyl-CoA, or an enolic intermediate sharing a proton with His-274 in a “low-barrier” hydrogen bond have been suggested. Similar problems of intermediate instability have been noted in other enzymic carbon acid deprotonation reactions. Quantum mechanical/molecular mechanical calculations of the pathway of acetyl-CoA enolization within citrate synthase support the identification of Asp-375 as the catalytic base. His-274, the proposed general acid, is found to be neutral. The acetyl-CoA enolate is more stable at the active site than the enol, and is stabilized by hydrogen bonds from His-274 and a water molecule. The conditions for formation of a low-barrier hydrogen bond do not appear to be met, and the calculated hydrogen bond stabilization in the reaction is less than the gas-phase energy, due to interactions with Asp-375 at the active site. The enolate character of the intermediate is apparently necessary for the condensation reaction to proceed efficiently. Proteins 27:9–25 © 1997 Wiley-Liss, Inc.     

Molecular orbital evaluation of charge flow dynamics in natural pigments based photosensitizers

The relationship between structure and photo electrochemical property of ten natural pigments from plants, insects and microbes has been analyzed using density functional theory (DFT) at the B3LYP/6-31G(d) level. The essential parameters for their photoelectrochemical behaviour such as ground state geometries, electronic transition energies and oxidation potentials are computed. The attachment tendency of the anchoring groups, expressed as the deprotonation order, is determined by calculating the proton affinities at different sites of the molecules. A thorough analysis of the charge flow dynamics in the molecular orbitals (HOMO and LUMO) of these molecules has been carried out and presented to emphasize the role of these orbitals in effective charge separation, the important feature of photosensitizers for DSSC. This study highlights that the flexible spatial orientation provided by the bridging aliphatic unsaturation favours the oscillator strength and the hydroxyl anchor group attached to the ring of delocalized π electron cloud acts as the effective anchor.     

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.     

Electrophilic catalysis in triosephosphate isomerase: the role of histidine-95

Electrophilic catalysis by histidine-95 in triosephosphate isomerase has been probed by using Fourier transform infrared spectroscopy and X-ray crystallography. The carbonyl stretching frequency of dihydroxyacetone phosphate bound to the wild-type enzyme is known to be 19 cm-1 lower (at 1713 cm-1) than that of dihydroxyacetone phosphate free in solution (at 1732 cm-1), and this decrease in stretching frequency has been ascribed to an enzymic electrophile that polarizes the substrate carbonyl group toward the transition state for the enolization. Infrared spectra of substrate bound to two site-directed mutants of yeast triosephosphate isomerase in which histidine-95 has been changed to glutamine or to asparagine show unperturbed carbonyl stretching frequencies between 1732 and 1742 cm-1. The lack of carbonyl polarization when histidine-95 is removed suggests that histidine-95 is indeed the catalytic electrophile, at least for dihydroxyacetone phosphate. Kinetic studies of the glutamine mutant (H95Q) have shown that the enzyme follows a subtly different mechanism of proton transfers involving only a single acid-base catalytic group. These findings suggest an additional role for histidine-95 as a general acid-base catalyst in the wild-type enzyme. The X-ray crystal structure of the H95Q mutant with an intermediate analogue, phosphoglycolohydroxamate, bound at the active site has been solved to 2.8-A resolution, and this structure clearly implicates glutamate-165, the catalytic base in the wild-type isomerase, as the sole acid-base catalyst for the mutant enzyme.(ABSTRACT TRUNCATED AT 250 WORDS)     

Calculation of isotope effects from first principles

Various means of calculating the effect of changing the mass of a given atom upon a chemical process are reviewed. Of particular interest is the deuterium isotope effect comparing the normal protium nucleus with its heavier deuterium congener. The replacement of the bridging protium in a neutral hydrogen bond such as the water dimer by a deuterium strengthens the interaction by a small amount via effects upon the vibrational energy. In an ionic H-bond such as the protonated water dimer, on the other hand, the reverse trend is observed in that replacement of the bridging protium by dimer weakens the interaction. In addition to the stability of a given complex, the rate at which a proton transfers from one group to another is likewise affected by deuterium substitution, viz. kinetic isotope effects (KIEs). The KIE is enlarged as the temperature drops, particularly so if the calculation of KIE includes proton tunneling. The KIE is also sensitive to any angular distortions or stretches present in the H-bond of interest. KIEs can be computed either by the standard transition state theory which is derived via only two points on the potential energy surface, or by more complete formalisms which take account of larger swaths of the surface. While more time intensive, the latter can also be applied to provide insights important in interpretation of experimental data.     

Near-ultraviolet magnetic circular dichroism spectroscopy of protein conformational states: correlation of tryptophan band position and intensity with hemoglobin allostery

The near-UV magnetic circular dichroism spectroscopy of the aromatic amino acid bands of hemoglobin was investigated as a potential probe of structural changes at the alpha(1)beta(2) interface during the allosteric transition. Allosteric effectors were used to direct carp and chemically modified human hemoglobins into the R (relaxed) or T (tense) state in order to determine the heme-ligation-independent spectral characteristics of the quaternary states. The tryptophan magnetic circular dichroism (MCD) peak observed at 293 nm in the R state of N-ethylsuccinimide- (NES-) des-Arg-modified human hemoglobin (Hb) was shifted to a slightly longer wavelength in the T state, consistent with the shift expected for tryptophan acting as a proton donor in a T-state hydrogen bond. Moreover, the increase observed in the T-state MCD intensity of this band relative to the R-state intensity was consistent with the effect expected for proton donation by tryptophan on the basis of the Michl perimeter model of aromatic MCD. The peak-to-trough magnitude of the R - T MCD difference spectrum is equal to 30% of the total R-state peak intensity contributed by all six tryptophans present in the human tetramer; the relative magnitude specific to the two beta37 tryptophans undergoing conformational change is estimated accordingly to be 3 times larger. The Trp-beta37 spectral shift, about 200 cm(-)(1), is in good agreement with the shifts observed in other H-bonded proton donors and provides corroborating spectral evidence for the formation in solution of a T-state Trp beta37-Asp alpha94 hydrogen bond observed in X-ray diffraction studies of deoxyHb crystals.     

Secondary 18O isotope effects as a tool for studying reactions of phosphate mono-, di-, and triesters

Secondary 18O isotope effects have been developed as a tool for determining transition state structures in enzymatic and nonenzymatic phosphoryl transfer reactions. 18O substitution in the nonbridge oxygens of a phosphoryl group makes the reaction go faster when the bond order is higher to these oxygens in the transition state than in the reactant, whereas the reaction goes slower if the bond order is less. The isotope effects are measured by the remote label method, using an isotope ratio mass spectrometer for analysis. The bond order to p-nitrophenolate ion when it is the leaving group is indicated by the secondary 15N isotope effect in the nitro group, with a value of 1.0028 representing nearly complete bond cleavage. It appears that the transition states for phosphoryl transfer have no more than one negative charge on the nonbridge oxygens, so that reactions of monoesters are dissociative, reactions of triesters are associative, and reactions of diesters are SN2 with half bond order to entering and leaving groups.     

A proposal for the mechanism of carbonic anhydrase action

A new catalytic mechanism is proposed for the hydration of CO₂ by the zinc metalloenzyme carbonic anhydrase. This mechanism identifies the group controlling catalytic activity as an active site histidine, in which the protonated imidazole ring coordinates to zinc, losing a proton. Geometric constraints on the histidine unit make the metal-ligand bond a strained and, therefore, labile one. In the hydration reaction, the metal-bound neutral histidine moiety serves as a proton acceptor for the transient ionization of metal-bound water. Zinc-bound hydroxide attacks the carbon of the substrate to generate metal-bound bicarbonate, and the system regenerates itself by losing the elements of carbonic acid.