Ions in water: characterizing the forces that control chemical processes and biological structure.

The performance of different link atom based frontier treatments in QM/MM simulations was evaluated critically with SCC-DFTB as the QM method. In addition to the analysis of gas-phase molecules as in previous studies, an important element of the present work is that chemical reactions in realistic enzyme systems were also examined. The schemes tested include all options available in the program CHARMM for SCC-DFTB/MM simulation, which treat electrostatic interactions due to the MM atoms close to the QM/MM boundary in different ways. In addition, a new approach, the divided frontier charge (DIV), has been implemented in which the partial charge associated with the frontier MM atom ("link host") is evenly distributed to the other MM atoms in the same group. The performance of these schemes was evaluated based on properties including proton affinities, deprotonation energies, dipole moments, and energetics of proton transfer reactions. Similar to previous work, it was found that calculated proton affinities and deprotonation energies of alcohols, carbonic acids, amino acids, and model DNA bases are very sensitive to the link atom scheme; the commonly used single link atom approach often gives error on the order of 15 to 20 kcal/mol. Other schemes give better and, on average, mutually comparable results. For proton transfer reactions, encouragingly, both activation barriers and reaction energies are fairly insensitive (within a typical range of 2-4 kcal/mol) to the link atom scheme due to error cancellation, and this was observed for both gas-phase and enzyme systems. Therefore, the effect of using different link atom schemes in QM/MM simulations is rather small for chemical reactions that conserve the total charge. Although the current study used an approximate DFT method as the QM level, the observed trends are expected to be applicable to QM/MM methods with use of other QM approaches. This observation does not mean to encourage QM/MM simulations without careful benchmark in the study of specific systems, rather it emphasizes that other technical details, such as the treatment of long-range electrostatics, tend to play a more important role and need to be handled carefully.

A critical evaluation of different QM/MM frontier treatments with SCC-DFTB as the QM method.

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 a...

Hydrogen Bonds: Simple after All?

In enzyme catalysis, where exquisitely positioned functionality is the sine qua non, atomic coordinates for a Michaelis complex can provide powerful insights into activation of the substrate. We focus here on the initial proton transfer of the isomerization reaction catalyzed by triosephosphate isomerase and present the crystal structure of its Michaelis complex with the substrate dihydroxyacetone phosphate at near-atomic resolution. The active site is highly compact, with unusually short and bifurcated hydrogen bonds for both catalytic Glu-165 and His-95 residues. The carboxylate oxygen of the catalytic base Glu-165 is positioned in an unprecedented close interaction with the ketone and the alpha-hydroxy carbons of the substrate (C em leader O approximately 3.0 A), which is optimal for the proton transfer involving these centers. The electrophile that polarizes the substrate, His-95, has close contacts to the substrate's O1 and O2 (N em leader O < or = 3.0 and 2.6 A, respectively). The substrate is conformationally relaxed in the Michaelis complex: the phosphate group is out of the plane of the ketone group, and the hydroxy and ketone oxygen atoms are not in the cisoid configuration. The epsilon ammonium group of the electrophilic Lys-12 is within hydrogen-bonding distance of the substrate's ketone oxygen, the bridging oxygen, and a terminal phosphate's oxygen, suggesting a role for this residue in both catalysis and in controlling the flexibility of active-site loop.

/pdf/optimal-alignment-for-enzymatic-proton-transfer-structure-of-3j7ekeoen8.pdf

Optimal alignment for enzymatic proton transfer: structure of the Michaelis complex of triosephosphate isomerase at 1.2-A resolution.

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 A/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.

/pdf/testing-electrostatic-complementarity-in-enzyme-catalysis-1ito3ihlct.pdf

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

Short, strong hydrogen bonds on enzymes: NMR and mechanistic studies

Methylglyoxal synthase (MGS) and triosephosphate isomerase (TIM) share neither sequence nor structural similarities, yet the reactions catalyzed by both enzymes are similar, in that both initially convert dihydroxyacetone phosphate to a cis-enediolic intermediate. This enediolic intermediate is formed from the abstraction of the pro-S C3 proton of DHAP by Asp-71 of MGS or the pro-R C3 proton of DHAP by Glu-165 of TIM. MGS then catalyzes the elimination of phosphate from this enediolic intermediate to form the enol of methylglyoxal, while TIM catalyzes proton donation to C2 to form d-glyceraldehyde phosphate. A competitive inhibitor of TIM, phosphoglycolohydroxamic acid (PGH) is found to be a tight binding competitive inhibitor of MGS with a Ki of 39 nM. PGH's high affinity for MGS may be due in part to a short, strong hydrogen bond (SSHB) from the NOH of PGH to the carboxylate of Asp-71. Evidence for this SSHB is found in X-ray, 1H NMR, and fractionation factor data. The X-ray structure of the MGS homohex...

Mechanistic implications of methylglyoxal synthase complexed with phosphoglycolohydroxamic acid as observed by X-ray crystallography and NMR spectroscopy.

Two detailed mechanisms [Marks et al. (2001) Biochemistry 40, 6805] have been proposed to explain the activity of methylglyoxal synthase (MGS), a homohexameric allosterically regulated enzyme that catalyzes the elimination of phosphate from DHAP to form enol pyruvaldehyde. This enol then tautomerizes to methylglyoxal in solution. In one of these mechanisms His 98 plays an obligate role in the transfer of a proton from the O(3) oxygen of DHAP to the O2 oxygen. To test this hypothesized mechanism, the variants H98N and H98Q were expressed and purified. Relative to the wild-type enzyme, the H98N variant shows a 50-fold decrease in k(cat) with all other kinetic parameters (i.e., K(m), K(PGA), etc.) being nearly the same. By contrast, the apparent catalytic rate for the H98Q variant is 250-fold lower than that of the wild-type enzyme. Inorganic phosphate acts as a competitive inhibitor (with a K(i) of 15 microM) rather than as an allosteric-type inhibitor as it does in the wild-type enzyme, and the competitive inhibitor phosphoglyolate (PGA) acts as an activator of this variant. These facts are explained by a shift in the conformational equilibrium toward an "inactive" state. When activation by PGA is accounted for, the catalytic rate for the "active" state of H98Q is estimated to be only 6-fold less than that of the wild-type enzyme, and thus His 98 is not essential for catalysis. Primary deuterium isotope effect data were measured for the wild-type enzyme and the two variants. At pH 7.0, the (D)V isotope effect (1.5) and the absence of a (D)(V/K) isotope effect for the wild-type enzyme suggest that the rate for the isotopically sensitive step is partially rate limiting but greater than the rate of substrate dissociation. Both the (D)V (2.0) and (D)(V/K) (3.4) isotope effects were more pronounced in the H98N variant, while the H98Q variant displayed an unusual inverse (D)V (0.8) isotope effect and a normal (D)(V/K) (1.5) isotope effect. The X-ray crystal structures of PGA bound to the H98Q and H98N variants were both determined to a resolution of 2.2 A. These mutations had little effect on the overall structure. However, the pattern of hydrogen bonding in the active site suggests an explanation as to how in the wild-type and H98N mutated enzymes the "inactive to active" equilibrium lies toward the active state, while with the H98Q mutated enzyme the equilibrium lies toward the inactive state. Thus, the role of His 98 appears to be more as a regulator of the enzyme's conformation rather than as a critical contributor to the catalytic mechanism.

Gregory T. Marks

Papers

Short, strong hydrogen bonds on enzymes: NMR and mechanistic studies

Mechanistic implications of methylglyoxal synthase complexed with phosphoglycolohydroxamic acid as observed by X-ray crystallography and NMR spectroscopy.

Mutagenic studies on histidine 98 of methylglyoxal synthase: effects on mechanism and conformational change.