Probing protein environment in an enzymatic process: All-electron quantum chemical analysis combined with ab initio quantum mechanical/molecular mechanical modeling of chorismate mutase

The influence of model building schemes and molecular dynamics sampling on QM-cluster models: the chorismate mutase case study

Most QM-cluster models of enzymes are constructed based on X-ray crystal structures, which limits comparison to in vivo structure and mechanism. The active site of chorismate mutase from Bacillus subtilis and the enzymatic transformation of chorismate to prephenate is used as a case study to guide construction of QM-cluster models built first from the X-ray crystal structure, then from molecular dynamics (MD) simulation snapshots. The Residue Interaction Network ResidUe Selector (RINRUS) software toolkit, developed by our group to simplify and automate the construction of QM-cluster models, is expanded to handle MD to QM-cluster model workflows. Several options, some employing novel topological clustering from residue interaction network (RIN) information, are evaluated for generating conformational clustering from MD simulation. RINRUS then generates a statistical thermodynamic framework for QM-cluster modeling of the chorismate mutase mechanism via refining 250 MD frames with density functional theory (DFT). The 250 QM-cluster models sampled provide a mean ΔG‡ of 10.3 ± 2.6 kcal mol-1 compared to the experimental value of 15.4 kcal mol-1 at 25 °C. While the difference between theory and experiment is consequential, the level of theory used is modest and therefore "chemical" accuracy is unexpected. More important are the comparisons made between QM-cluster models designed from the X-ray crystal structure versus those from MD frames. The large variations in kinetic and thermodynamic properties arise from geometric changes in the ensemble of QM-cluster models, rather from the composition of the QM-cluster models or from the active site-solvent interface. The findings open the way for further quantitative and reproducible calibration in the field of computational enzymology using the model construction framework afforded with the RINRUS software toolkit.

Free Energy Profile Decomposition Analysis for QM/MM Simulations of Enzymatic Reactions

In enzyme mechanistic studies and mutant design, it is highly desirable to know the individual residue contributions to the reaction free energy and barrier. In this work, we show that such free energy contributions from each residue can be readily obtained by postprocessing ab initio quantum mechanical molecular mechanical (ai-QM/MM) free energy simulation trajectories. Specifically, through a mean force integration along the minimum free energy pathway, one can obtain the electrostatic, polarization, and van der Waals contributions from each residue to the free energy barrier. Separately, a similar analysis procedure allows us to assess the contribution from different collective variables along the reaction coordinate. The chorismate mutase reaction is used to demonstrate the utilization of these two trajectory analysis tools.

Insight into the Catalytic Mechanism of GH11 Xylanase: Computational Analysis of Substrate Distortion Based on a Neutron Structure

The reaction mechanism of biomass decomposition by xylanases remains the subject of debate. To clarify the mechanism we investigated the glycosylation step of GH11 xylanase, an enzyme that catalyzes the hydrolysis of lignocellulosic hemicellulose (xylan). Making use of a recent neutron crystal structure, which revealed the protonation states of relevant residues, we used ab initio quantum mechanics/molecular mechanics (QM/MM) calculations to determine the detailed reaction mechanism of the glycosylation step. In particular, our focus is on the controversial question of whether or not an oxocarbenium ion intermediate is formed on the reaction pathway. The calculations support the validity of a basic retaining mechanism within a double-displacement scheme. The estimated free energy barrier of this reaction is ∼18 kcal/mol with QM/MM-CCSD(T)/6-31(+)G**//MP2/6-31+G**/AMBER calculations, and the rate-determining step of the glycosylation is scission of the glycosidic bond after proton transfer from the acidic Glu177. The estimated lifetime of the oxocarbenium ion intermediate (on the order of tens of ps) and the secondary kinetic isotope effect suggest that there is no accumulation of this intermediate on the reaction path, although the intermediate can be transiently formed. In the enzyme-substrate (ES) complex, the carbohydrate structure of the xylose residue at the -1 subsite has a rather distorted (skewed) geometry, and this xylose unit at the active site has an apparent half-chair conformation when the oxocarbenium ion intermediate is formed. The major catalytic role of the protein environment is to orient residues that take part in the initial proton transfer. Because of a fine alignment of catalytic residues, the enzyme can accelerate the glycosylation reaction without paying a reorganization energy penalty.

How the Local Environment of Functional Sites Regulates Protein Function

Proteins form complex biological machineries whose functions in the cell are highly regulated at both the cellular and molecular levels. Cellular regulation of protein functions involves differential gene expressions, post-translation modifications, and signaling cascades. Molecular regulation, on the other hand, involves tuning an optimal local protein environment for the functional site. Precisely how a protein achieves such an optimal environment around a given functional site is not well understood. Herein, by surveying the literature, we first summarize the various reported strategies used by certain proteins to ensure their correct functioning. We then formulate three key physicochemical factors for regulating a protein's functional site, namely, (i) its immediate interactions, (ii) its solvent accessibility, and (iii) its conformational flexibility. We illustrate how these factors are applied to regulate the functions of free/metal-bound Cys and Zn sites in proteins.

Computational analysis of carbohydrate recognition based on hybrid QM/MM modeling: a case study of norovirus capsid protein in complex with Lewis antigen

Norovirus is a major pathogen of nonbacterial acute gastroenteritis in humans and animals. Carbohydrate recognition between norovirus capsid proteins and Lewis antigens is considered to play a critical role in initiating infection of eukaryotic cells. In this article, we first report a detailed atomistic simulation study of the norovirus capsid protein in complex with the Lewis antigen based on ab initio QM/MM combined with MD-FEP simulations. To understand the mechanistic details of ligand binding, we analyzed and compared the carbohydrate recognition mechanism of the wild-type P domain protein with a mutant protein. Small structural differences between two capsid proteins are observed on the weak interaction site of residue 389, which is located on the solvent exposed surface of the P domain. To further clarify affinity differences in ligand binding, we directly evaluated free energy changes of the ligand binding process. Although the mutant protein loses its interaction energy with the Lewis antigen, this small amount of energy penalty is compensated for by an increase in the solvation stability, which is induced by structural reorganization at the ligand binding site on the protein surface. As a sum of these opposite energy components, the mutant P domain obtains a slightly enhanced binding affinity for the Lewis antigen. The present computational study clearly demonstrated that a detailed free energy balance of the interaction energy between the capsid protein and the surrounding aqueous solvent is the mechanistic basis of carbohydrate recognition in the norovirus capsid protein.

Exploration of swapping enzymatic function between two proteins: a simulation study of chorismate mutase and isochorismate pyruvate lyase

The enzyme chorismate mutase EcCM from Escherichia coli catalyzes one of the few pericyclic reactions in biology, the transformation of chorismate to prephenate. The isochorismate pyruvate lyase PchB from Pseudomonas aeroginosa catalyzes another pericyclic reaction, the isochorismate to salicylate transformation. Interestingly, PchB possesses weak chorismate mutase activity as well thus being able to catalyze two distinct pericyclic reactions in a single active site. EcCM and PchB possess very similar folds, despite their low sequence identity. Using molecular dynamics simulations of four combinations of the two enzymes (EcCM and PchB) with the two substrates (chorismate and isochorismate) we show that the electrostatic field due to EcCM at atoms of chorismate favors the chorismate to prephenate transition and that, analogously, the electrostatic field due to PchB at atoms of isochorismate favors the isochorismate to salicylate transition. The largest differences between EcCM and PchB in electrostatic field strengths at atoms of the substrates are found to be due to residue side chains at distances between 0.6 and 0.8 nm from particular substrate atoms. Both enzymes tend to bring their non-native substrate in the same conformation as their native substrate. EcCM and to a lower extent PchB fail in influencing the forces on and conformations of the substrate such as to favor the other chemical reaction (isochorismate pyruvate lyase activity for EcCM and chorismate mutase activity for PchB). These observations might explain the difficulty of engineering isochorismate pyruvate lyase activity in EcCM by solely mutating active site residues.

Substrate distortion contributes to the catalysis of orotidine 5'-monophosphate decarboxylase

Orotidine 5'-monophosphate decarboxylase (ODCase) accelerates the decarboxylation of orotidine 5'-monophosphate (OMP) to uridine 5'-monophosphate (UMP) by 17 orders of magnitude. Eight new crystal structures with ligand analogues combined with computational analyses of the enzyme's short-lived intermediates and the intrinsic electronic energies to distort the substrate and other ligands improve our understanding of the still controversially discussed reaction mechanism. In their respective complexes, 6-methyl-UMP displays significant distortion of its methyl substituent bond, 6-amino-UMP shows the competition between the K72 and C6 substituents for a position close to D70, and the methyl and ethyl esters of OMP both induce rotation of the carboxylate group substituent out of the plane of the pyrimidine ring. Molecular dynamics and quantum mechanics/molecular mechanics computations of the enzyme-substrate complex also show the bond between the carboxylate group and the pyrimidine ring to be distorted, with the distortion contributing a 10-15% decrease of the ΔΔG(⧧) value. These results are consistent with ODCase using both substrate distortion and transition-state stabilization, primarily exerted by K72, in its catalysis of the OMP decarboxylation reaction.

Analysis of chorismate mutase catalysis by QM/MM modelling of enzyme-catalysed and uncatalysed reactions

Chorismate mutase is at the centre of current controversy about fundamental features of biological catalysts. Some recent studies have proposed that catalysis in this enzyme does not involve transition state (TS) stabilization but instead is due largely to the formation of a reactive conformation of the substrate. To understand the origins of catalysis, it is necessary to compare equivalent reactions in different environments. The pericyclic conversion of chorismate to prephenate catalysed by chorismate mutase also occurs (much more slowly) in aqueous solution. In this study we analyse the origins of catalysis by comparison of multiple quantum mechanics/molecular mechanics (QM/MM) reaction pathways at a reliable, well tested level of theory (B3LYP/6-31G(d)/CHARMM27) for the reaction (i) in Bacillus subtilis chorismate mutase (BsCM) and (ii) in aqueous solvent. The average calculated reaction (potential energy) barriers are 11.3 kcal mol(-1) in the enzyme and 17.4 kcal mol(-1) in water, both of which are in good agreement with experiment. Comparison of the two sets of reaction pathways shows that the reaction follows a slightly different reaction pathway in the enzyme than in it does in solution, because of a destabilization, or strain, of the substrate in the enzyme. The substrate strain energy within the enzyme remains constant throughout the reaction. There is no unique reactive conformation of the substrate common to both environments, and the transition state structures are also different in the enzyme and in water. Analysis of the barrier heights in each environment shows a clear correlation between TS stabilization and the barrier height. The average differential TS stabilization is 7.3 kcal mol(-1) in the enzyme. This is significantly higher than the small amount of TS stabilization in water (on average only 1.0 kcal mol(-1) relative to the substrate). The TS is stabilized mainly by electrostatic interactions with active site residues in the enzyme, with Arg90, Arg7 and Glu78 generally the most important. Conformational effects (e.g. strain of the substrate in the enzyme) do not contribute significantly to the lower barrier observed in the enzyme. The results show that catalysis is mainly due to better TS stabilization by the enzyme.

Structural evolution of protein-biofilms: Simulations and experiments

The control of biofilm formation is a challenging goal that has not been reached yet in many aspects. One unsolved question is the role of van der Waals forces and another is the importance of mutual interactions between the adsorbing and the adsorbed biomolecules ("critical crowding"). In this study, a combined experimental and theoretical approach is presented, which fundamentally probes both aspects. On three model proteins-lysozyme, α-amylase, and bovine serum albumin-the adsorption kinetics is studied experimentally. Composite substrates are used enabling a separation of the short- and the long-range forces. Although usually neglected, experimental evidence is given for the influence of van der Waals forces on the protein adsorption as revealed by in situ ellipsometry. The three proteins were chosen for their different conformational stabilities in order to investigate the influence of conformational changes on the adsorption kinetics. Monte Carlo simulations are used to develop a model for these experimental results by assuming an internal degree of freedom to represent conformational changes. The simulations also provide data on the distribution of adsorption sites. By in situ atomic force microscopy we can also test this distribution experimentally, which opens the possibility to, e.g., investigate the interactions between adsorbed proteins.

Computational modeling of carbohydrate-recognition process in E-selectin complex: structural mapping of sialyl Lewis X onto ab initio QM/MM free energy surface

To advance our knowledge of carbohydrate recognition by lectins, we propose a systematic computational modeling strategy to identify complex sugar-chain conformations on the reduced free energy surface (FES). We selected the complex of E-selectin with sialyl Lewis X (denoted E-selectin/SLe(x) complex) as a first target molecule. First, we introduced the reduced 2D-FES that characterizes conformational changes in carbohydrate structure as well as the degree of solvation stability of the carbohydrate ligand, and evaluated the overall free energy profile by classical molecular dynamics simulation combined with ab initio QM/MM energy corrections. Second, we mapped flexible carbohydrate structures onto the reduced QM/MM 2D-FES, and identified the details of molecular interactions between each monosaccharide component and the amino acid residues at the carbohydrate-recognition domain. Finally, we confirmed the validity of our modeling strategy by evaluating the chemical shielding tensor by ab initio QM/MM-GIAO computations for several QM/MM-refined geometries sampled from the minimum free energy region in the 2D-FES, and compared this theoretical averaging data with the experimental 1D-NMR profile. The model clearly shows that the binding geometries of the E-selectin/SLe(x) complex are determined not by one single, rigid carbohydrate structure but rather by the sum of averaged conformations fluctuating around the minimum free energy region. For the E-selectin/SLe(x) complex, the major molecular interactions are hydrogen bonds between Fuc and the Ca(2+) binding site in the carbohydrate-recognition domain, and Gal is important in determining the ligand specificity.

Finite-temperature effects in enzymatic reactions — Insights from QM/MM free-energy simulations

We report potential-energy and free-energy data for three enzymatic reactions: carbon–halogen bond formation in fluorinase, hydrogen abstraction from camphor in cytochrome P450cam, and chorismate-to-prephenate Claisen rearrangement in chorismate mutase. The results were obtained by combined quantum mechanics/molecular mechanics (QM/MM) optimizations and two types of QM/MM free-energy simulations (free-energy perturbation and umbrella sampling) using semi-empirical or density-functional QM methods. Based on these results and our previously published free-energy data on electrophilic substitution in para-hydroxybenzoate hydroxylase, we discuss the importance of finite-temperature effects in the chemical step of enzyme reactions. We find that the entropic contribution to the activation barrier is generally rather small, usually of the order of 5 kJ mol–1 or less, consistent with the notion that enzymes bind and pre-organize the reactants in the active site. A somewhat larger entropic contribution is encountered in the case of chorismate mutase where the pericyclic transition state is intrinsically more rigid than the chorismate reactant (also in the enzyme). The present results suggest that barriers from QM/MM geometry optimization may often be close to free-energy barriers for the chemical step in enzymatic reactions.