Abstract:A comparative theoretical study of a bimolecular reaction in aqueous solution and catalyzed by the enzyme catechol O-methyltransferase (COMT) has been carried out by a combination of two hybrid QM/MM techniques: statistical simulation methods and internal energy minimizations. In contrast to previous studies by other workers, we have located and characterized transition structures for the reaction in the enzyme active site, in water and in a vacuum, and our potential of mean force calculations are based upon reaction coordinates obtained from features of the potential energy surfaces in the condensed media, not from the gas phase. The AM1/CHARMM calculated free energy of activation for the reaction of S-adenosyl methionine (SAM) with catecholate catalyzed by COMT is 15 kcal mol(-1) lower the AM1/TIP3P free-energy barrier for the reaction of the trimethylsulfonium cation with the catecholate anion in water at 300 K, in agreement with previous estimates. The thermodynamically preferred form of the reactants in the uncatalyzed model reaction in water is a solvent-separated ion pair (SSIP). Conversion of the SSIP into a contact ion pair, with a structure resembling that of the Michaelis complex (MC) for the reaction in the COMT active site, is unfavorable by 7 kcal mol(-1), largely due to reorganization of the solvent. We have considered alternative ways to estimate the so-called "cratic" free energy for bringing the reactant species together in the correct orientation for reaction but conclude that direct evaluation of the free energy of association by means of molecular dynamics simulation with a simple standard-state correction is probably the best approach. The latter correction allows for the fact that the size of the unit cell employed with the periodic boundary simulations does not correspond to the standard state concentration of 1 M. Consideration of MC-like species allows a helpful decomposition of the catalytic effect into preorganization and reorganization phases. In the preorganization phase, the substrates are brought together into the MC-like species, either in water or in the enzyme active site. In the reorganization phase, the roles of the enzymic and aqueous environments may be compared directly because reorganization of the substrate is about the same in both cases. Analysis of the electric field along the reaction coordinate demonstrates that in water the TS is destabilized with respect to the MC-like species because the polarity of the solute diminishes and consequently the reaction field is also decreased. In the enzyme, the electric field is mainly a permanent field and consequently there is only a small reorganization of the environment. Therefore, destabilization of the TS is lower than in solution, and the activation barrier is smaller.

Modeling of Transition State by Molecular Dynamics. Prediction of Catalytic Efficiency of the Mutants of Mandelate Racemase

The Role Of Residue S139 Of Mandelate Racemase: Synergistic Effect Of S139 And E317 On Transition State Stabilization

Virtual screening of mandelate racemase mutants with enhanced activity based on binding energy in the transition state

The effect of Val22 on the catalytic activity of mandelate racemase from Pseudomonas putida

Directed Evolution of Mandelate Racemase by a Novel High-Throughput Screening Method

A Computational Method to Predict Effects of Residue Mutations on the Catalytic Efficiency of Hydrolases

Quantifying the error propagation in microkinetic modeling of catalytic reactions with model-predicted binding energies

A Regressor-Guided Graph Diffusion Model for Predicting Enzyme Mutations to Enhance Turnover Number

Capturing the free energy of transition state stabilization: insights from the inhibition of mandelate racemase

Improving the Catalytic Activity of Lipase LipK107 from Proteus Sp. by Site-Directed Mutagenesis in the Lid Domain Based on Computer Simulation

On the Importance of Preserving Active Induced-Fit Conformations during Enzyme-Catalyzed Reactions: a Case Study of CARNMT1

Molecular Dynamics Optimization of a Computational Model of TACE and its Substrate Peptide

Electrostatic transition state stabilization rather than reactant destabilization provides the chemical basis for efficient chorismate mutase catalysis

Unraveling the energetic significance of chemical events in enzyme catalysis via machine-learning based regression approach

Rapid screening of multi-point mutations for enzyme thermostability modification by utilizing computational tools

Steady State Kinetics for Enzymes with Multiple Binding Sites Upstream of the Catalytic Site

Theoretical modeling of enzyme catalytic power: analysis of "cratic" and electrostatic factors in catechol O-methyltransferase

The Influence of Model Building Schemes and Molecular Dynamics Sampling on QM-cluster Models: The Chorismate Mutase Case Study

The Infuence of Model Building Schemes and Molecular Dynamics Sampling on QM-cluster Models: The Chorismate Mutase Case Study

Precise redesign for improving enzyme robustness based on coevolutionary analysis and multidimensional virtual screening

Mechanism of Mutation-Induced Effects on the Catalytic Function of TEV Protease: A Molecular Dynamics Study