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.     

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.     

On the relationship between thermal stability and catalytic power of enzymes

The possible relationship between the thermal stability and the catalytic power of enzymes is of great current interest. In particular, it has been suggested that thermophilic or hyperthermophilic (Tm) enzymes have lower catalytic power at a given temperature than the corresponding mesophilic (Ms) enzymes, because the thermophilic enzymes are less flexible (assuming that flexibility and catalysis are directly correlated). These suggestions presume that the reduced dynamics of the thermophilic enzymes is the reason for their reduced catalytic power. The present paper takes the specific case of dihydrofolate reductase (DHFR) and explores the validity of the above argument by simulation approaches. It is found that the Tm enzymes have restricted motions in the direction of the folding coordinate, but this is not relevant to the chemical process, since the motions along the reaction coordinate are perpendicular to the folding motions. Moreover, it is shown that the rate of the chemical reaction is determined by the activation barrier and the corresponding reorganization energy, rather than by dynamics or flexibility in the ground state. In fact, as far as flexibility is concerned, we conclude that the displacement along the reaction coordinate is larger in the Tm enzyme than in the Ms enzyme and that the general trend in enzyme catalysis is that the best catalyst involves less motion during the reaction than the less optimal catalyst. The relationship between thermal stability and catalysis appears to reflect the fact that to obtain small electrostatic reorganization energy it is necessary to invest some folding energy in the overall preorganization process. Thus, the optimized catalysts are less stable. This trend is clearly observed in the DHFR case.     

Electrostatic contributions to binding of transition state analogues can be very different from the corresponding contributions to catalysis: phenolates binding to the oxyanion hole of ketosteroid isomerase

The relationship between binding of transition state analogues (TSAs) and catalysis is an open problem. A recent study of the binding of phenolate TSAs to ketosteroid isomerase (KSI) found a small change in the binding energy with a change in charge delocalization of the TSAs. This has been taken as proof that electrostatic effects do not contribute in a major way to catalysis. Here we reanalyze the relationship between the binding of the TSAs and the chemical catalysis by KSI as well as the binding of the transition state (TS), by computer simulation approaches. Since the simulations reproduce the relevant experimental results, they can be used to quantify the different contributions to the observed effects. It is found that the binding of the TSA and the chemical catalysis represent different thermodynamic cycles with very different electrostatic contributions. While the binding of the TSA involves a small electrostatic contribution, the chemical catalysis involves a charge transfer process and a major electrostatic contribution due to the preorganization of the active site. Furthermore, it is found that the electrostatic preorganization contributions to the binding of the enolate intermediate of KSI and the TS are much larger than the corresponding effect for the binding of the TSAs. This reflects the dependence of the preorganization on the orientation of the nonpolar form of the TSAs relative to the oxyanion hole. It seems to us that this work provides an excellent example of the need for computational studies in analyzing key experimental findings about enzyme catalysis.     

Phospholipase A2 at the bilayer interface.

Interfacial catalysis is a necessary consequence for all enzymes that act on amphipathic substrates with a strong tendency to form aggregates in aqueous dispersions. In such cases the catalytic event occurs at the interface of the aggregated substrate, the overall turnover at the interface is processive, and it is influenced the molecular organization and dynamics of the interface. Such enzymes can access the substrate only at the interface because the concentration of solitary monomers of the substrate in the aqueous phase is very low. Moreover, the microinterface between the bound enzyme and the organized substrate not only facilitates formation of the enzyme-substrate complex, but a longer residence time of the enzyme at the substrate interface also promotes high catalytic processivity. Binding of the enzyme to the substrate interface as an additional step in the overall catalytic turnover permits adaptation of the Michaelis-Menten formalism as a basis to account for the kinetics of interfacial catalysis. As shown for the action of phospholipase A2 on bilayer vesicles, binding equilibrium has two extreme kinetic consequences. During catalysis in the scooting mode the enzyme does not leave the surface of the vesicle to which it is bound. On the other hand, in the hopping mode the absorption and desorption steps are a part of the catalytic turnover. In this minireview we elaborate on the factors that control binding of pig pancreatic phospholipase A2 to the bilayer interface. Binding of PLA2 to the interface occurs through ionic interactions and is further promoted by hydrophobic interactions which probably occur along a face of the enzyme, with a hydrophobic collar and a ring of cationic residues, through which the catalytic site is accessible to substrate molecules in the bilayer. An enzyme molecule binds to the surface occupied by about 35 lipid molecules with an apparent dissociation constant of less than 0.1 pM for the enzyme on anionic vesicles compared to 10 mM on zwitterionic vesicles. Results at hand also show that aggregation or acylation of the protein is not required for the high affinity binding or catalytic interaction at the interface.     

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.     

Challenges and advances in validating enzyme design proposals: the case of kemp eliminase catalysis

One of the fundamental challenges in biotechnology and biochemistry is the ability to design effective enzymes. Despite recent progress, most of the advances on this front have been made by placing the reacting fragments in the proper places, rather than by optimizing the preorganization of the environment, which is the key factor in enzyme catalysis. Thus, rational improvement of the preorganization would require approaches capable of evaluating reliably the actual catalytic effect. This work considers the catalytic effects in different Kemp eliminases as a benchmark for a computer-aided enzyme design. It is shown that the empirical valence bond provides a powerful screening tool, with significant advantages over current alternative strategies. The insights provided by the empirical valence bond calculations are discussed with an emphasis on the ability to analyze the difference between the linear free energy relationships obtained in solution and those found in the enzymes. We also point out the trade-off between the reliability and speed of the calculations and try to determine what it takes to realize reliable computer-aided screening.     

What are the roles of substrate-assisted catalysis and proximity effects in peptide bond formation by the ribosome?

The action of the peptidyl transferase center of the large ribosomal unit presents a fundamental step in the evolution from the RNA world to the protein world. Thus, it is important to understand the origin of the catalytic power of this ancient enzyme. Earlier studies suggested that the ribosome catalyzes peptide bond formation by using one of its groups as a general base, while more recent works have proposed that the catalysis is due to proximity effects or to substrate-assisted catalysis. However, the actual nature of the catalytic mechanism remains controversial. This work addresses the origin of the catalytic power of the ribosome by using computer simulation approaches and comparing the energetics of the peptide bond formation in the ribosome and in water. It is found that a significant part of the observed activation entropy of the reference solution reaction is due to solvation entropy, and that the proximity effect is smaller than previously thought. It is also found that the 2'-OH of the A76 ribose, which is associated with a large rate acceleration in the ribosome reaction, does not catalyze peptide bond formation in water. Thus, the catalytic effect cannot be attributed to substrate-assisted catalysis but rather to the effect of the ribosome on the reacting system. Overall, our calculations indicate that the reduction of the activation free energy is mainly due to electrostatic effects. The nature of these effects and their relationship to catalytic factors in modern enzymes is analyzed and discussed.     

Structural analyses of covalent enzyme-substrate analog complexes reveal strengths and limitations of de novo enzyme design

We report the cocrystal structures of a computationally designed and experimentally optimized retro-aldol enzyme with covalently bound substrate analogs. The structure with a covalently bound mechanism-based inhibitor is similar to, but not identical with, the design model, with an RMSD of 1.4 Å over active-site residues and equivalent substrate atoms. As in the design model, the binding pocket orients the substrate through hydrophobic interactions with the naphthyl moiety such that the oxygen atoms analogous to the carbinolamine and β-hydroxyl oxygens are positioned near a network of bound waters. However, there are differences between the design model and the structure: the orientation of the naphthyl group and the conformation of the catalytic lysine are slightly different; the bound water network appears to be more extensive; and the bound substrate analog exhibits more conformational heterogeneity than typical native enzyme–inhibitor complexes. Alanine scanning of the active-site residues shows that both the catalytic lysine and the residues around the binding pocket for the substrate naphthyl group make critical contributions to catalysis. Mutating the set of water-coordinating residues also significantly reduces catalytic activity. The crystal structure of the enzyme with a smaller substrate analog that lacks naphthyl ring shows the catalytic lysine to be more flexible than in the naphthyl–substrate complex; increased preorganization of the active site would likely improve catalysis. The covalently bound complex structures and mutagenesis data highlight the strengths and weaknesses of the de novo enzyme design strategy.     

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.     

A twin-track approach has optimized proton and hydride transfer by dynamically coupled tunneling during the evolution of protochlorophyllide oxidoreductase

Protein dynamics are crucial for realizing the catalytic power of enzymes, but how enzymes have evolved to achieve catalysis is unknown. The light-activated enzyme protochlorophyllide oxidoreductase (POR) catalyzes sequential hydride and proton transfers in the photoexcited and ground states, respectively, and is an excellent system for relating the effects of motions to catalysis. Here, we have used the temperature dependence of isotope effects and solvent viscosity measurements to analyze the dynamics coupled to the hydride and proton transfer steps in three cyanobacterial PORs and a related plant enzyme. We have related the dynamic profiles of each enzyme to their evolutionary origin. Motions coupled to light-driven hydride transfer are conserved across all POR enzymes, but those linked to thermally activated proton transfer are variable. Cyanobacterial PORs require complex and solvent-coupled dynamic networks to optimize the proton donor-acceptor distance, but evolutionary pressures appear to have minimized such networks in plant PORs. POR from Gloeobacter violaceus has features of both the cyanobacterial and plant enzymes, suggesting that the dynamic properties have been optimized during the evolution of POR. We infer that the differing trajectories in optimizing a catalytic structure are related to the stringency of the chemistry catalyzed and define a functional adaptation in which active site chemistry is protected from the dynamic effects of distal mutations that might otherwise impact negatively on enzyme catalysis.     

pH Dependence of catalysis by Pseudomonas aeruginosa isochorismate-pyruvate lyase: implications for transition state stabilization and the role of lysine 42

An isochorismate-pyruvate lyase with adventitious chorismate mutase activity from Pseudomonas aerugionsa (PchB) achieves catalysis of both pericyclic reactions in part by the stabilization of reactive conformations and in part by electrostatic transition-state stabilization. When the active site loop Lys42 is mutated to histidine, the enzyme develops a pH dependence corresponding to a loss of catalytic power upon deprotonation of the histidine. Structural data indicate that the change is not due to changes in active site architecture, but due to the difference in charge at this key site. With loss of the positive charge on the K42H side chain at high pH, the enzyme retains lyase activity at ～100-fold lowered catalytic efficiency but loses detectable mutase activity. We propose that both substrate organization and electrostatic transition state stabilization contribute to catalysis. However, the dominant reaction path for catalysis is dependent on reaction conditions, which influence the electrostatic properties of the enzyme active site amino acid side chains.     

Effect of oxidatively damaged DNA on the active site preorganization during nucleotide incorporation in a high fidelity polymerase from Bacillus stearothermophilus

We study the effect of the oxidative lesion 8-oxoguanine (8oxoG) on the preorganization of the active site for DNA replication in the closed (active) state of the Bacillus fragment (BF), a Klenow analog from Bacillus stearothermophilus. Our molecular dynamics and free energy simulations of explicitly solvated model ternary complexes of BF bound to correct dCTP/incorrect dATP opposite guanine (G) and 8oxoG bases in DNA suggest that the lesion introduces structural and energetic changes at the catalytic site to favor dATP insertion. Despite the formation of a stable Watson-Crick pairing in the 8oxoG:dCTP system, the catalytic geometry is severely distorted to possibly slow down catalysis. Indeed, our calculated free energy landscapes associated with active site preorganization suggest additional barriers to assemble an efficient catalytic site, which need to be overcome during dCTP incorporation opposite 8oxoG relative to that opposite undamaged G. In contrast, the catalytic geometry for the Hoogsteen pairing in the 8oxoG:dATP system is highly organized and poised for efficient nucleotide incorporation via the "two-metal-ion" catalyzed phosphoryl transfer mechanism. However, the free energy calculations suggest that the catalytic geometry during dATP incorporation opposite 8oxoG is considerably less plastic than that during dCTP incorporation opposite G despite a very similar, well organized catalytic site for both systems. A correlation analysis of the dynamics trajectories suggests the presence of significant coupling between motions of the polymerase fingers and the primary distance for nucleophilic attack (i.e., between the terminal primer O3' and the dNTP P(alpha.) atoms) during correct dCTP incorporation opposite undamaged G. This coupling is shown to be disrupted during nucleotide incorporation by the polymerase with oxidatively damaged DNA/dNTP substrates. We also suggest that the lesion affects DNA interactions with key polymerase residues, thereby affecting the enzymes ability to discriminate against non-complementary DNA/dNTP substrates. Taken together, our results provide a unified structural, energetic, and dynamic platform to rationalize experimentally observed relative nucleotide incorporation rates for correct dCTP/incorrect dATP insertion opposite an undamaged/oxidatively damaged template G by BF.     

Transition state theory can be used in studies of enzyme catalysis: lessons from simulations of tunnelling and dynamical effects in lipoxygenase and other systems

The idea that enzyme catalysis involves special factors such as coherent fluctuations, quantum mechanical tunnelling and non-equilibrium solvation (NES) effects has gained popularity in recent years. It has also been suggested that transition state theory (TST) cannot be used in studies of enzyme catalysis. The present work uses reliable state of the art simulation approaches to examine the above ideas. We start by demonstrating that we are able to simulate any of the present catalytic proposals using the empirical valence bond (EVB) potential energy surfaces, the dispersed polaron model and the quantized classical path (QCP) approach, as well as the approximate vibronic method. These approaches do not treat the catalytic effects by phenomenological treatments and thus can be considered as first principles approaches (at least their ability to compare enzymatic reaction to the corresponding solution reactions). This work will consider the lipoxygenase reaction, and to lesser extent other enzymes, for specific demonstration. It will be pointed out that our study of the lipoxygenase reaction reproduces the very large observed isotope effect and the observed rate constant while obtaining no catalytic contribution from nuclear quantum mechanical (NQM) effects. Furthermore, it will be clarified that our studies established that the NQM effect decreases rather than increases when the donor-acceptor distance is compressed. The consequences of these findings in terms of the temperature dependence of the kinetic isotope effect and in terms of different catalytic proposals will be discussed. This paper will also consider briefly the dynamical effects and conclude that such effects do not contribute in a significant way to enzyme catalysis. Furthermore, it will be pointed out that, in contrast to recent suggestions, NES effects are not dynamical effects and should therefore be part of the activation free energy rather than the transmission factor. In view of findings of the present work and our earlier works, it seems that TST provides a quantitative tool for studies of enzyme catalysis and that the key open questions are related to the nature of the factors that lead to transition state stabilization.     

Pairwise specificity and sequential binding in enzyme catalysis: thymidylate synthase

The structures of thymidylate synthase (TS) from Escherichia coli, in ternary complexes with substrate and an analogue of the cofactor, are the basis of a stereochemical model for a key reaction intermediate in the catalyzed reaction. This model is used to compare the reaction chemistry and chirality of the transferred methyl group with structures of the components, to identify those residues that participate, and to propose a stereochemical mechanism for catalysis by TS. Effects of chemical modification of specific amino acid residues and site-directed mutations of residues are correlated with structure and effects on enzyme mechanism. The ordered binding sequence of substrate deoxyuridine monophosphate and methylenetetrahydrofolate can be understood from the structure, where each forms a large part of the binding site for the other. The catalytic site serves to orient the reactants, which are sequestered along with many water molecules within a cavernous active center. Conformational changes during the reaction could involve nearby residues in ways that are not obvious in this complex.     

Nonconserved residues Ala287 and Ser290 of the Cryptosporidium hominis thymidylate synthase domain facilitate its rapid rate of catalysis

Cryptosporidium hominis TS-DHFR exhibits an unusually high rate of catalysis at the TS domain, at least 10-fold greater than those of other TS enzymes. Using site-directed mutagenesis, we have mutated residues Ala287 and Ser290 in the folate-binding helix to phenylalanine and glycine, respectively, the corresponding residues in human and most other TS enzymes. Our results show that the mutant A287F, the mutant S290G, and the double mutant all have reduced affinities for methylene tetrahydrofolate and reduced rates of reaction at the TS domain. Interestingly, the S290G mutant enzyme had the lowest TS activity, with a catalytic efficiency approximately 200-fold lower than that of the wild type (WT). The rate of conformational change of the S290G mutant is approximately 80 times slower than that of WT, resulting in a change in the rate-limiting step from hydride transfer to covalent ternary complex formation. We have determined the crystal structure of ligand-bound S290G mutant enzyme, which shows that the primary effect of the mutation is an increase in the distance between the TS ligands. The kinetic and crystal structure data presented here provide the first evidence explaining the unusually fast TS rate in C. hominis.     

Catalytic cysteine of thymidylate synthase is activated upon substrate binding

The role of Ser 167 of Escherichia coli thymidylate synthase (TS) in catalysis has been characterized by kinetic and crystallographic studies. Position 167 variants including S167A, S167N, S167D, S167C, S167G, S167L, S167T, and S167V were generated by site-directed mutagenesis. Only S167A, S167G, S167T, and S167C complemented the growth of thymidine auxotrophs of E. coli in medium lacking thymidine. Steady-state kinetic analysis revealed that mutant enzymes exhibited k(cat) values 1.1-95-fold lower than that of the wild-type enzyme. Relative to wild-type TS, K(m) values of the mutant enzymes for 2'-deoxyuridylate (dUMP) were 5-90 times higher, while K(m) values for 5,10-methylenetetrahydrofolate (CH(2)H(4)folate) were 1.5-16-fold higher. The rate of dehalogenation of 5-bromo-2'-deoxyuridine 5'-monophosphate (BrdUMP), a reaction catalyzed by TS that does not require CH(2)H(4)folate as cosubstrate, by mutant TSs was analyzed and showed that only S167A and S167G catalyzed the dehalogenation reaction and values of k(cat)/K(m) for the mutant enzymes were decreased by 10- and 3000-fold, respectively. Analysis of pre-steady-state kinetics of ternary complex formation revealed that the productive binding of CH(2)H(4)folate is weaker to mutant TSs than to the wild-type enzyme. Chemical transformation constants (k(chem)) for the mutant enzymes were lower by 1.1-6.0-fold relative to the wild-type enzyme. S167A, S167T, and S167C crystallized in the I2(1)3 space group and scattered X-rays to either 1.7 A (S167A and S167T) or 2.6 A (S167C). The high-resolution data sets were refined to a R(crys) of 19.9%. In the crystals some cysteine residues were derivatized with 2-mercaptoethanol to form S,S-(2-hydroxyethyl)thiocysteine. The pattern of derivatization indicates that in the absence of bound substrate the catalytic cysteine is not more reactive than other cysteines. It is proposed that the catalytic cysteine is activated by substrate binding by a proton-transfer mechanism in which the phosphate group of the nucleotide neutralizes the charge of Arg 126', facilitating the transfer of a proton from the catalytic cysteine to a His 207-Asp 205 diad via a system of ordered water molecules.     

Sequence and structural differences between enzyme and nonenzyme homologs

To improve our understanding of the evolution of novel functions, we performed a sequence, structural, and functional analysis of homologous enzymes and nonenzymes of known three-dimensional structure. In most examples identified, the nonenzyme is derived from an ancestral catalytic precursor (as opposed to the reverse evolutionary scenario, nonenzyme to enzyme), and the active site pocket has been disrupted in some way, owing to the substitution of critical catalytic residues and/or steric interactions that impede substrate binding and catalysis. Pairwise sequence identity is typically insignificant, and almost one-half of the enzyme and nonenzyme pairs do not share any similarity in function. Heterooligomeric enzymes comprising homologous subunits in which one chain is catalytically inactive and enzyme polypeptides that contain internal catalytic and noncatalytic duplications of an ancient enzyme domain are also discussed.