The empirical valence bond as an effective strategy for computer-aided enzyme design

The ability of the empirical valence bond (EVB) to be used in screening active site residues in enzyme design is explored in a preliminary way. This validation is done by comparing the ability of this approach to evaluate the catalytic contributions of various residues in chorismate mutase. It is demonstrated that the EVB model can serve as an accurate tool in the final stages of computer-aided enzyme design (CAED). The ability of the model to predict quantitatively the catalytic power of enzymes should augment the capacity of current approaches for enzyme design.     

The catalytic effect of dihydrofolate reductase and its mutants is determined by reorganization energies

The effect of distant mutations on the catalytic reaction of dihydrofolate reductase (DHFR) is reexamined by empirical valence bond simulations. The simulations reproduce for the first time the changes in the observed rate constants (without the use of adjustable parameters for this purpose) and show that the changes in activation barriers are strongly correlated with the corresponding changes in the reorganization energy. The preorganization of the polar groups of enzymes is the key catalytic factor, and anticatalytic mutations destroy this preorganization. Some anticatalytic mutations in DHFR also increase the distance between the donor and acceptor, but this effect is not directly related to catalysis since the native enzyme and the uncatalyzed reaction in water have similar average donor-acceptor distances. Insight into the effect of a mutation is provided by constructing the relevant free energy surfaces in terms of the generalized solute-solvent coordinates. It is shown how the mutations change the reaction coordinate and the activation barrier, and it is clarified that the corresponding changes do not reflect dynamical effects. It is also pointed out that all reactions in a condensed phase involve correlated motions (both in enzymes and in solution) and that the change of such motions upon mutations is a result of the change in the shape of the multidimensional reaction path on the solute-solvent surface, rather than the reason for the change in rate constant. Thus, as far as catalysis is concerned, the change in the activation barrier is due to the change in the electrostatic preorganization energy.     

Methyltransferases do not work by compression,cratic, or desolvation effects,but by electrostatic preorganization

The enzyme catechol O‐methyltransferase (COMT) catalyzes the transfer of a methyl group from S‐adenosylmethionine to dopamine and related catechols. The search for the origin of COMT catalysis has led to different proposals and hypothesis, including the entropic, the NAC, and the compression proposals as well as the more reasonable electrostatic idea. Thus, it is important to understand the catalytic power of this enzyme and to examine the validity of different proposals and in particular the repeated recent implication of the compression idea. The corresponding analysis should be done by well‐defined physically‐based considerations that involve computations rather than circular interpretations of experimental results. Thus, we explore here the origin of the catalytic efficiency of COMT by using the empirical valence bond and the linear response approximation approaches. The results demonstrate that the catalytic effect of COMT is mainly due to electrostatic preorganization effects. It is also shown that the compression, NAC and entropic proposals do not account for the catalytic effect. Proteins 2015; 83:318–330. © 2014 Wiley Periodicals, Inc.     

Misunderstanding the preorganization concept can lead to confusions about the origin of enzyme catalysis

Understanding the origin of the catalytic power of enzymes has both conceptual and practical importance. One of the most important finding from computational studies of enzyme catalysis is that a major part of the catalytic power is due to the preorganization of the enzyme active site. Unfortunately, misunderstanding of the nontrivial preorganization idea lead some to assume that it does not consider the effect of the protein residues. This major confusion reflects a misunderstanding of the statement that the interaction energy of the enzyme group and the transition state (TS) is similar to the corresponding interaction between the water molecules (in the reference system) and the TS, and that the catalysis is due to the reorganization free energy of the water molecules. Obviously, this finding does not mean that we do not consider the enzyme groups. Another problem is the idea that catalysis is due to substrate preorganization. This more traditional idea is based in some cases on inconsistent interpretation of the action of model compounds, which unfortunately, do not reflect the actual situation in the enzyme active site. The present article addresses the above problems, clarifying first the enzyme polar preorganization idea and the current misunderstandings. Next we take a specific model compound that was used to promote the substrate preorganization proposal and establish its irrelevance to enzyme catalysis. Overall, we show that the origin of the catalytic power of enzymes cannot be assessed uniquely without computer simulations, since at present this is the only way of relating structure and energetics.     

Advances in optimizing enzyme electrostatic preorganization

Utilizing electric fields to catalyze chemical reactions is not a new idea, but in enzymology it undergoes a renaissance, inspired by Warhsel's concept of electrostatic preorganization. According to this concept, the source of the immense catalytic efficiency of enzymes is the intramolecular electric field that permanently favors the reaction transition state over the reactants. Within enzyme design, computational efforts have fallen short in designing enzymes with natural-like efficacy. The outcome could improve if long-range electrostatics (often omitted in current protocols) would be optimized. Here, we highlight the major developments in methods for analyzing and designing electric fields generated by the protein scaffolds, in order to both better understand how natural enzymes function, and aid artificial enzyme design.     

Computer simulation of the initial proton transfer step in human carbonic anhydrase I.

The initial water proteolysis step in the proton transfer "half-reaction" of human carbonic anhydrase I is simulated using the empirical valence bond method in combination with free energy perturbation molecular dynamics calculations. A free energy profile for the enzyme catalysed reaction and the corresponding pKa associated with ionization of the zinc-bound water is calculated. The obtained pKa value of 7 to 8 appears to be in good agreement with experimental observations and the calculated rate constant for this step is also compatible with kinetic data. The simulations clearly emphasize the important electrostatic effect associated with the catalytic zinc ion.     

Catalysis and linear free energy relationships in aspartic proteases

Aspartic proteases are receiving considerable attention as potential drug targets in several serious diseases, such as AIDS, malaria, and Alzheimer's disease. These enzymes cleave polypeptide chains, often between specific amino acid residues, but despite the common reaction mechanism, they exhibit large structural differences. Here, the catalytic mechanism of aspartic proteases plasmepsin II, cathepsin D, and HIV-1 protease is examined by computer simulations utilizing the empirical valence bond approach in combination with molecular dynamics and free energy perturbation calculations. Free energy profiles are established for four different substrates, each six amino acids long and containing hydrophobic side chains in the P1 and P1' positions. Our simulations reproduce the catalytic effect of these enzymes, which accelerate the reaction rate by a factor of approximately 10(10) compared to that of the corresponding uncatalyzed reaction in water. The calculations elucidate the origin of the catalytic effect and allow a rationalization of the fact that, despite large structural differences between plasmepsin II/cathepsin D and HIV-1 protease, the magnitude of their rate enhancement is very similar. Amino acid residues surrounding the active site together with structurally conserved water molecules are found to play an important role in catalysis, mainly through dipolar (electrostatic) stabilization. A linear free energy relationship for the reactions in the different enzymes is established that also demonstrates the reduced reorganization energy in the enzymes compared to that in the uncatalyzed water reaction.     

Semiquantitative calculations of catalytic free energies in genetically modified enzymes

The catalytic free energy and binding free energies of the native and the Asn-155----Thr, Asn-155----Leu, and Asn-155----Ala mutants of subtilisin are calculated by the empirical valence bond method and a free energy perturbation method. Two simple procedures are used; one "mutates" the substrate, and the other "mutates" the enzyme. The calculated changes in free energies (delta delta G not equal to cat and delta delta Gbind) between the mutant and native enzymes are within 1 kcal/mol of the corresponding observed values. This indicates that we are approaching a quantitative structure-function correlation. The calculated changes in catalytic free energies are almost entirely due to the electrostatic interaction between the enzyme-water system and the charges of the reacting system. This supports the idea that the electrostatic free energy associated with the changes of charges of the reacting system is the key factor in enzyme catalysis.     

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.     

Enhancement of 7-methoxyresorufin O-demethylation activity of human cytochrome P450 1A2 by molecular breeding

Alkylresorufins are model substrates for cytochrome P450 (P450) 1A2. The ability of human P450 1A2 to catalyze 7-methoxyresorufin O-demethylation was improved by screening of random mutant libraries (expressed in Escherichia coli) on the basis of 7-methoxyresorufin O-demethylation. After three rounds of mutagenesis and screening, the triple mutant E163K/V193M/K170Q yielded a kcat > five times faster than wild type P450 1A2 in steady-state kinetic analysis using either isolated membrane fractions or purified, reconstituted enzymes. The enhanced catalytic activity was not attributed to changes in substrate affinity. The kinetic hydrogen isotope effect of the triple mutant did not change from wild type enzyme and suggests that C-H bond cleavage is rate-limiting in both enzymes. Homology modeling, based on an X-ray structure of rabbit P450 2C5, suggests that the locations of mutated residues are not close to the substrate binding site and therefore that structural elements outside of this site play roles in changing the catalytic activity. This approach has potential value in understanding P450 1A2 and generating engineered enzymes with enhanced catalytic activity.     

Calculations of free energy profiles for the staphylococcal nuclease catalyzed reaction

Calculations of the free energy profile for the first two (rate-limiting) steps of the staphylococcal nuclease catalyzed reaction are reported. The calculations are based on the empirical valence bond method in combination with free energy perturbation molecular dynamics simulations. The calculated activation free energy is in good agreement with experimental kinetic data, and the catalytic effect of the enzyme is reproduced without any arbitrary adjustment of parameters. The enormous reduction of the activation barrier (relative to the reference reaction in water) appears to be largely associated with the strong electrostatic effect of the Ca2+ ion and the two arginine residues in the active site. This favorable electrostatic environment reduces the cost of the general-base catalysis step by almost 15 kcal/mol (by stabilizing the OH- nucleophile) and then stabilizes the developing negative charge on the 5'-phosphate group in the second step of the reaction by about 19 kcal/mol. The basic features of the originally postulated enzyme mechanism (Cotton et al., 1979) are found to be compatible with the observed activation free energy. However, the proposed modification of the mechanism (Sepersu et al., 1987), in which Arg 87 interacts only with the pentacoordinated transition state, is supported by the simulations. Further calculations on the D21E mutant also give results in good agreement with kinetic data.     

基于序列和结构分析的酶热稳定性改造策略*

功能酶被广泛应用于食品、化工、医药等领域,但却容易受高温环境限制,导致催化效率降低。以分子改造为目的的蛋白质工程技术是解决这一问题的关键环节,其能够对酶结构和功能进行改造,获得热稳定性好的工业酶。传统的定向进化方法只能依靠随机突变进行人工筛选,具有效率低、针对性差等缺点;理性设计作为酶热稳定性改造的主要方法,可借助各种计算机程序和软件预测潜在突变位点,但其要求对酶的催化机制、热稳定性机制有深入了解。对于大多数天然酶而言,酶的序列和晶体结构是最容易获取的信息,也是预测功能的重要基础。从酶的序列和晶体结构入手,重点介绍了共识突变、基于序列偏好性的突变、截短柔性区域、优化分子内相互作用力、刚化催化活性区域及计算机辅助筛选柔性位点等常用策略,这些策略具有筛选效率高、改造准确性高、实用性强等优点。结合多种酶的热稳定性改造案例进行分析,旨在为不同酶的改造策略选择提供有效参考,同时也为工业酶的耐热性研究提供理论支持。     

Global and local molecular dynamics of a bacterial carboxylesterase provide insight into its catalytic mechanism

Carboxylesterases (CEs) are ubiquitous enzymes responsible for the detoxification of xenobiotics. In humans, substrates for these enzymes are far-ranging, and include the street drug heroin and the anticancer agent irinotecan. Hence, their ability to bind and metabolize substrates is of broad interest to biomedical science. In this study, we focused our attention on dynamic motions of a CE from B. subtilis (pnbCE), with emphasis on the question of what individual domains of the enzyme might contribute to its catalytic activity. We used a 10 ns all-atom molecular dynamics simulation, normal mode calculations, and enzyme kinetics to understand catalytic consequences of structural changes within this enzyme. Our results shed light on how molecular motions are coupled with catalysis. During molecular dynamics, we observed a distinct C-C bond rotation between two conformations of Glu310. Such a bond rotation would alternately facilitate and impede protonation of the active site His399 and act as a mechanism by which the enzyme alternates between its active and inactive conformation. Our normal mode results demonstrate that the distinct low-frequency motions of two loops in pnbCE, coil_5 and coil_21, are important in substrate conversion and seal the active site. Mutant CEs lacking these external loops show significantly reduced rates of substrate conversion, suggesting this sealing motion prevents escape of substrate. Overall, the results of our studies give new insight into the structure-function relationship of CEs and have implications for the entire family of α/β fold family of hydrolases, of which this CE is a member.     

Understanding substrate specificity in human and parasite phosphoribosyltransferases through calculation and experiment.

We present molecular dynamics (MD) simulations on two enzymes: a human hypoxanthine-guanine-phosphoribosyltransferase (HGPRTase) and its analogue in the protozoan parasite Tritrichomonas foetus. The parasite enzyme has an additional ability to process xanthine as a substrate, making it a hypoxanthine-guanine-xanthine phosphoribosyltransferase (HGXPRTase) [Chin, M. S., and Wang, C. C. (1994) Mol. Biochem. Parasitol. 63 (2), 221-229 (1)]. X-ray crystal structures of both enzymes complexed to guanine monoribosyl phosphate (GMP) have been solved, and show only subtle differences in the two active sites [Eads et al. (1994) Cell 78 (2), 325-334 (2); Somoza et al. (1996) Biochemistry 35 (22), 7032-7040 (3)]. Most of the direct contacts with the base region of the substrate are made by the protein backbone, complicating the identification of residues significantly associated with xanthine recognition. Our calculations suggest that the broader specificity of the parasite enzyme is due to a significantly more flexible base-binding region, and rationalize the effect of two mutations, R155E and D163N, that alter substrate specificity [Munagala, N. R., and Wang, C. C. (1998) Biochemistry 37 (47), 16612-16619 (4)]. In addition, our simulations suggested a double mutant (D106E/D163N) that might rescue the D163N mutant. This double mutant was expressed and assayed, and its catalytic activity was confirmed. Our molecular dynamics trajectories were also used with a structure-based design program, Pictorial Representation Of Free Energy Changes (PROFEC), to suggest parasite-selective derivatives of GMP. Our calculations here successfully rationalize the parasite-selectivity of two novel inhibitors derived from the computer-aided design of Somoza et al. (5) and demonstrate the utility of PROFEC in the design of species-selective inhibitors.     

Validating computer simulations of enantioselective catalysis; reproducing the large steric and entropic contributions in Candida Antarctica lipase B

The prospect for computer‐aided refinement of stereoselective enzymes is further validated by simulating the ester hydrolysis by the wild‐type and mutants of CalB, focusing on the challenge of dealing with strong steric effects and entropic contributions. This was done using the empirical valence bond (EVB) method in a quantitative screening of the enantioselectivity, considering both k_cat and k_cat/K_M of the R and S stereoisomers. Although the simulations require very extensive sampling for convergence they give encouraging results and major validation, indicating that our approach offers a powerful tool for computer‐aided design of enantioselective enzymes. This is particularly true in cases with large changes in steric effects where alternative approaches may have difficulties in capturing the interplay between steric clashes with the reacting substrate and protein flexibility. Proteins 2014; 82:1387–1399. © 2014 Wiley Periodicals, Inc.     

Computer simulation of primary kinetic isotope effects in the proposed rate-limiting step of the glyoxalase I catalyzed reaction

The proposed rate-limiting step of the glyoxalase I catalyzed reaction is the proton abstraction from the C1 carbon of the substrate by Glu(172). Here we examine primary kinetic isotope effects and the influence of quantum dynamics on this process by computer simulations. The calculations utilize the empirical valence bond method in combination with the molecular dynamics free energy perturbation technique and path integral simulations. For the enzyme-catalyzed reaction a H/D kinetic isotope effect of 5.0 +/- 1. 3 is predicted in reasonable agreement with the experimental result of about 3. Furthermore, the magnitude of quantum mechanical effects is found to be very similar for the enzyme reaction and the corresponding uncatalyzed process in solution, in agreement with other studies. The problems associated with attaining the required accuracy in order for the present approach to be useful as a diagnostic tool for the study of enzyme reactions are also discussed.     

Hydroxyl groups in the (beta)beta sandwich of metallo-beta-lactamases favor enzyme activity: a computational protein design study

Metallo-beta-lactamases challenge antimicrobial therapies by their ability to hydrolyze and inactivate a broad spectrum of beta-lactam antibiotics. The potential of these enzymes to acquire enhanced catalytic efficiency through mutation is of great concern. Here, we explore the potential of computational protein design to predict mutants of the imipenemase IMP-1 that modulate the catalytic efficiency of the enzyme against a range of substrates. Focusing on the four amino acid positions 69, 121, 218, and 262, we carried out a number of design calculations. Two mutant enzymes were predicted: the single mutant S262A and the double mutant F218Y-S262A. Compared to IMP-1, the single mutant (S262A) results in the loss of a hydroxyl group and the double mutant (F218Y-S262A) results in a hydroxyl transfer from position 262 to position 218. The presence of both hydroxyl groups at positions 218 and 262 was tested by examining the mutant F218Y. Kinetic constants of IMP-1, the two computationally designed mutants (S262A and F218Y-S262A), and the hydroxyl addition mutant (F218Y) were determined with seven substrates. Catalytic efficiencies are highest for the enzyme with both hydroxyl groups (F218Y) and lowest for the enzyme lacking both hydroxyl groups (S262A). The catalytic efficiencies of the two enzymes with one hydroxyl group each are intermediate, with the F218Y-S262A double mutant exhibiting enhanced hydrolysis of nitrocefin, cephalothin, and cefotaxime relative to IMP-1.     

Molecular dynamics simulation of proton transport near the surface of a phospholipid membrane

The structural and dynamical properties of a hydrated proton near the surface of DMPC membrane were studied using a molecular dynamics simulation. The proton transport between water molecules was modeled using the second generation multistate empirical valence bond model. The proton diffusion was found to be inhibited at the membrane surface. The potential of mean force for the proton adsorption to the membrane surface and its release back into the bulk water was also determined, yielding a small barrier in each direction. An efficient algorithm for Ewald summation calculations for the multistate empirical valence bond model is also introduced.     

Computational method for the design of enzymes with altered substrate specificity.

A combination of enzyme kinetics and X-ray crystallographic analysis of site-specific mutants has been used to probe the determinants of substrate specificity for the enzyme alpha-lytic protease. We now present a generalized model for understanding the effects of mutagenesis on enzyme substrate specificity. This algorithm uses a library of side-chain rotamers to sample conformation space within the binding site for the enzyme-substrate complex. The free energy of each conformation is evaluated with a standard molecular mechanics force field, modified to include a solvation energy term. This rapid energy calculation based on coarse conformation sampling quite accurately predicts the relative catalytic efficiency of over 40 different alpha-lytic protease-substrate combinations. Unlike other computational approaches, with this method it is feasible to evaluate all possible mutations within the binding site. Using this algorithm, we have successfully designed a protease that is both highly active and selective for a non-natural substrate. These encouraging results indicate that it is possible to design altered enzymes solely on the basis of empirical energy calculations.     

Glutamine 53 is a gatekeeper residue in the FK506 binding protein

The structures of enzymes reflect two tendencies that appear opposed. On one hand, they fold into compact, stable structures; on the other hand, they bind a ligand and catalyze a reaction. To be stable, enzymes fold to maximize favorable interactions, forming a tightly packed hydrophobic core, exposing hydrophilic groups, and optimizing intramolecular hydrogen-bonding. To be functional, enzymes carve out an active site for ligand binding, exposing hydrophobic surface area, clustering like charges, and providing unfulfilled hydrogen bond donors and acceptors. Using AmpC beta-lactamase, an enzyme that is well-characterized structurally and mechanistically, the relationship between enzyme stability and function was investigated by substituting key active-site residues and measuring the changes in stability and activity. Substitutions of catalytic residues Ser64, Lys67, Tyr150, Asn152, and Lys315 decrease the activity of the enzyme by 10(3)-10(5)-fold compared to wild-type. Concomitantly, many of these substitutions increase the stability of the enzyme significantly, by up to 4.7kcal/mol. To determine the structural origins of stabilization, the crystal structures of four mutant enzymes were determined to between 1.90A and 1.50A resolution. These structures revealed several mechanisms by which stability was increased, including mimicry of the substrate by the substituted residue (S64D), relief of steric strain (S64G), relief of electrostatic strain (K67Q), and improved polar complementarity (N152H). These results suggest that the preorganization of functionality characteristic of active sites has come at a considerable cost to enzyme stability. In proteins of unknown function, the presence of such destabilized regions may indicate the presence of a binding site.