All electron quantum chemical calculation of the entire enzyme system confirms a collective catalytic device in the chorismate mutase reaction

To elucidate the catalytic power of enzymes, we analyzed the reaction profile of Claisen rearrangement of Bacillus subtilis chorismate mutase (BsCM) by all electron quantum chemical calculations using the fragment molecular orbital (FMO) method. To the best of our knowledge, this is the first report of ab initio-based quantum chemical calculations of the entire enzyme system, where we provide a detailed analysis of the catalytic factors that accomplish transition-state stabilization (TSS). FMO calculations deliver an ab initio-level estimate of the intermolecular interaction between the substrate and the amino acid residues of the enzyme. To clarify the catalytic role of Arg90, we calculated the reaction profile of the wild-type BsCM as well as Lys90 and Cit90 mutant BsCMs. Structural refinement and the reaction path determination were performed at the ab initio QM/MM level, and FMO calculations were applied to the QM/MM refined structures. Comparison between three types of reactions established two collective catalytic factors in the BsCM reaction: (1) the hydrogen bonds connecting the Glu78-Arg90-substrate cooperatively control the stability of TS relative to the ES complex and (2) the positive charge on Arg90 polarizes the substrate in the TS region to gain more electrostatic stabilization.     

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.     

Transition state stabilization and substrate strain in enzyme catalysis: ab initio QM/MM modelling of the chorismate mutase reaction

To investigate fundamental features of enzyme catalysis, there is a need for high-level calculations capable of modelling crucial, unstable species such as transition states as they are formed within enzymes. We have modelled an important model enzyme reaction, the Claisen rearrangement of chorismate to prephenate in chorismate mutase, by combined ab initio quantum mechanics/molecular mechanics (QM/MM) methods. The best estimates of the potential energy barrier in the enzyme are 7.4-11.0 kcal mol(-1)(MP2/6-31+G(d)//6-31G(d)/CHARMM22) and 12.7-16.1 kcal mol(-1)(B3LYP/6-311+G(2d,p)//6-31G(d)/CHARMM22), comparable to the experimental estimate of Delta H(++)= 12.7 +/- 0.4 kcal mol(-1). The results provide unequivocal evidence of transition state (TS) stabilization by the enzyme, with contributions from residues Arg90, Arg7, and Arg63. Glu78 stabilizes the prephenate product (relative to substrate), and can also stabilize the TS. Examination of the same pathway in solution (with a variety of continuum models), at the same ab initio levels, allows comparison of the catalyzed and uncatalyzed reactions. Calculated barriers in solution are 28.0 kcal mol(-1)(MP2/6-31+G(d)/PCM) and 24.6 kcal mol(-1)(B3LYP/6-311+G(2d,p)/PCM), comparable to the experimental finding of Delta G(++)= 25.4 kcal mol(-1) and consistent with the experimentally-deduced 10(6)-fold rate acceleration by the enzyme. The substrate is found to be significantly distorted in the enzyme, adopting a structure closer to the transition state, although the degree of compression is less than predicted by lower-level calculations. This apparent substrate strain, or compression, is potentially also catalytically relevant. Solution calculations, however, suggest that the catalytic contribution of this compression may be relatively small. Consideration of the same reaction pathway in solution and in the enzyme, involving reaction from a 'near-attack conformer' of the substrate, indicates that adoption of this conformation is not in itself a major contribution to catalysis. Transition state stabilization (by electrostatic interactions, including hydrogen bonds) is found to be central to catalysis by the enzyme. Several hydrogen bonds are observed to shorten at the TS. The active site is clearly complementary to the transition state for the reaction, stabilizing it more than the substrate, so reducing the barrier to reaction.     

What is so special about Arg 55 in the catalysis of cyclophilin A? insights from hybrid QM/MM simulations

Potential of mean force (PMF) simulations with a hybrid QM/MM potential function were used to analyze the catalytic mechanism of human cyclophilin A (CypA). PMF calculations were performed for proline isomerization of peptides in solution, the wild-type CypA, and several CypA mutants. With an approximate density functional theory, the self-consistent-charge density functional tight binding (SCC-DFTB) as the QM level, and CHARMM 22 force field as MM, satisfactory energetics compared to available experiments were obtained. Calculations for the Arg55Ala and zero-charge-Arg55 mutants clearly indicated that Arg 55 significantly stabilizes the isomerization transition state through electrostatic interactions. However, the decrease in the average distance (thus the increase in interaction) between Arg 55 and the substrate amide N in going from the stable states to the transition state is mainly due to the pyramidalization of the amide N rather than motions associated with Arg 55. Although the nanosecond simulations cannot exclude the existence of sub-millisecond collective motions proposed on the basis of recent elegant NMR relaxation and line-shape analyses, the energetics obtained for the various enzyme systems here indicate that the contribution from motions of active site residues to catalysis is expected to be small. Instead, the present simulations support that the structural stability rather than mobility of the preorganized active site is more important. Through hydrogen-bonding interactions among the substrate, Arg 55, Gln 63, and Asn 102, the active site of the wild-type enzyme is structurally very stable and puts Arg 55 in a favorable position to perform its catalytic role in the transition state. This is further illustrated with the somewhat unexpected prediction that Arg55Lys is largely catalytically inactive, because Lys does not have the unique bifurcating construct of the guanidino group in Arg and thus the active site of Arg55Lys cannot accommodate Lys in a position capable of providing electrostatic stabilization of the isomerization transition state. Among all the enzyme systems studied, the wild-type CypA is the only one that selects the syn/exo transition state, while the syn/endo conformation is also present in the mutants, which is another reason for their higher barriers. Finally, the present analysis indicated that the population of near-attack-conformations (NAC) is not relevant to catalysis in CypA.     

Selective stabilization of the chorismate mutase transition state by a positively charged hydrogen bond donor

Citrulline was incorporated via chemical semisynthesis at position 90 in the active site of the AroH chorismate mutase from Bacillus subtilis. The wild-type arginine at this position makes hydrogen-bonding interactions with the ether oxygen of chorismate. Replacement of the positively charged guanidinium group with the isosteric but neutral urea has a dramatic effect on the ability of the enzyme to convert chorismate into prephenate. The Arg90Cit variant exhibits a >104-fold decrease in the catalytic rate constant kcat with a 2.7-fold increase in the Michaelis constant Km. In contrast, its affinity for a conformationally constrained inhibitor molecule that effectively mimics the geometry but not the dissociative character of the transition state is only reduced by a factor of approximately 6. These results show that an active site merely complementary to the reactive conformation of chorismate is insufficient for catalysis of the mutase reaction. Instead, electrostatic stabilization of the polarized transition state by provision of a cationic hydrogen bond donor proximal to the oxygen in the breaking C-O bond is essential for high catalytic efficiency.     

Physical nature of intermolecular interactions within cAMP-dependent protein kinase active site: differential transition state stabilization in phosphoryl transfer reaction

The origin of enzyme catalytic activity may be effectively explored within the nonempirical theory of intermolecular interactions. The knowledge of electrostatic, exchange, delocalization, and correlation components of the transition state and substrates stabilization energy arising from each enzyme active site residue allows to examine the most essential physical effects involved in enzymatic catalysis. Consequently, one can build approximate models of the catalytic activity in a systematic and legitimate manner. Whenever the dominant role of electrostatic interactions is recognized or assumed, the properties of an optimal catalytic environment could be simply generalized and visualized by means of catalytic fields that, in turn, aids the design of new catalysts. Differential transition state stabilization (DTSS) methodology has been applied herein to the phosphoryl transfer reaction catalyzed by cAMP-dependent protein kinase (PKA). The MP2 results correlate well with the available experimental data and theoretical findings indicating that Lys72, Asp166, and the two magnesium ions contribute -22.7, -13.3, -32.4, and -15.2 kcal/mol to differential transition state stabilization, respectively. Although all interaction energy components except that of electron correlation contribution are meaningful, the first-order electrostatic term correlates perfectly with MP2 catalytic activity. Catalytic field technique was also employed to visualize crucial electrostatic features of an ideal catalyst and to compare the latter with the environment provided by PKA active site. The map of regional electronic chemical potential was used to analyze the unfavorable catalytic effect of Lys168. It was found that locally induced polarization of TS atoms thermodynamically destabilizes electrons, pulling them to regions displaying higher electronic chemical potential.     

Contributions of conformational compression and preferential transition state stabilization to the rate enhancement by chorismate mutase

The rate enhancement provided by the chorismate mutase (CM) enzyme for the Claisen rearrangement of chorismate to prephenate has been investigated by application of the concept of near attack conformations (NACs). Using a combined QM/MM Monte Carlo/free-energy perturbation (MC/FEP) method, 82% and 100% of chorismate conformers were found to be NAC structures in water and in the CM active site, respectively. Consequently, the conversion of non-NACs to NACs does not contribute to the free energy of activation from preorganization of the substrate into NACs. The FEP calculations yielded differences in free energies of activation that well reproduce the experimental data. Additional calculations indicate that the rate enhancement by CM over the aqueous phase results primarily from conformational compression of NACs by the enzyme and that this process is enthalpically controlled. This suggests that preferential stabilization of the transition state in the enzyme environment relative to water plays a secondary role in the catalysis by CM.     

A hybrid potential reaction path and free energy study of the chorismate mutase reaction.

We present a combination of two techniques--QM/MM statistical simulation methods and QM/MM internal energy minimizations--to get a deeper insight into the reaction catalyzed by the enzyme chorismate mutase. Structures, internal energies and free energies, taken from the paths of the reaction in solution and in the enzyme have been analyzed in order to estimate the relative importance of the reorganization and preorganization effects. The results we obtain for this reaction are in good agreement with experiment and show that chorismate mutase achieves its catalytic efficiency in two ways; first, it preferentially binds the active conformer of the substrate and, second, it reduces the free energy of activation for the reaction relative to that in solution by providing an environment which stabilizes the transition state.     

A comparative study of claisen and cope rearrangements catalyzed by chorismate mutase. An insight into enzymatic efficiency: transition state stabilization or substrate preorganization?

In this work we present a detailed analysis of the activation free energies and averaged interactions for the Claisen and Cope rearrangements of chorismate and carbachorismate catalyzed by Bacillus subtilischorismate mutase (BsCM) using quantum mechanics/molecular mechanics (QM/MM) simulation methods. In gas phase, both reactions are described as concerted processes, with the activation free energy for carbachorismate being about 10-15 kcal mol(-)(1) larger than for chorismate, at the AM1 and B3LYP/6-31G levels. Aqueous solution and BsCM active site environments reduce the free energy barriers for both reactions, due to the fact that in these media the two carboxylate groups can be approached more easily than in the gas phase. The enzyme specifically reduces the activation free energy of the Claisen rearrangement about 3 kcal mol(-)(1) more than that for the Cope reaction. This result is due to a larger transition state stabilization associated to the formation of a hydrogen bond between Arg90 and the ether oxygen. When this oxygen atom is changed by a methylene group, the interaction is lost and Arg90 moves inside the active site establishing stronger interactions with one of the carboxylate groups. This fact yields a more intense rearrangement of the substrate structure. Comparing two reactions in the same enzyme, we have been able to obtain conclusions about the relative magnitude of the substrate preorganization and transition state stabilization effects. Transition state stabilization seems to be the dominant effect in this case.     

How does the cAMP-dependent protein kinase catalyze the phosphorylation reaction: an ab initio QM/MM study

We have carried out density functional theory QM/MM calculations on the catalytic subunit of cAMP-dependent protein kinase (PKA). The QM/MM calculations indicate that the phosphorylation reaction catalyzed by PKA is mainly dissociative, and Asp166 serves as the catalytic base to accept the proton delivered by the substrate peptide. Among the key interactions in the active site, the Mg(2+) ions, glycine rich loop, and Lys72 are found to stabilize the transition state through electrostatic interactions. On the other hand, Lys168, Asn171, Asp184, and the conserved waters bound to Mg(2+) ions do not directly contribute to lower the energy barrier of the phosphorylation reaction, and possible roles for these residues are proposed. The QM/MM calculations with different QM/MM partition schemes or different initial structures yield consistent results. In addition, we have carried out 12 ns molecular dynamics simulations on both wild type and K168A mutated PKA, respectively, to demonstrate that the catalytic role of Lys168 is to keep ATP and substrate peptide in the near-attack reactive conformation.     

Catalytic mechanism of glycosyltransferases: hybrid quantum mechanical/molecular mechanical study of the inverting N-acetylglucosaminyltransferase I

The Golgi glycosyltransferase, N-acetylglucosaminyltransferase I (GnT-I), catalyzes the transfer of a GlcNAc residue from the donor UDP-GlcNAc to the C2-hydroxyl group of a mannose residue in the trimannosyl core of the Man5GlcNAc2-Asn-X oligosaccharide. The catalytic mechanism of GnT-I was investigated using a hybrid quantum mechanical/molecular mechanical (QM/MM) method with a QM part containing 88 atoms treated with density functional theory (DFT) at the BP/TZP level. The remaining parts of a GnT-I complex, altogether 5633 atoms, were modeled using the AMBER molecular force field. A theoretical model of a Michaelis complex was built using the X-ray structure of GnT-I in complex with the donor having geometrical features consistent with kinetic studies. The QM(DFT)/MM model identified a concerted SN2-type of transition state with D291 as the catalytic base for the reaction in the enzyme active site. The TS model features nearly simultaneous nucleophilic addition and dissociation steps accompanied by the transfer of the nucleophile proton Hb2 to the catalytic base D291. The structure of the TS model is characterized by the Ob2-C1 and C1-O1 bond distances of 1.912 and 2.542 A, respectively. The activation energy for the proposed reaction mechanism was estimated to be approximately 19 kcal mol-1. The calculated alpha-deuterium kinetic isotope effect of 1.060 is consistent with the proposed reaction mechanism. Theoretical results also identified interactions between the Hb6 and beta-phosphate oxygen of the UDP and a low-barrier hydrogen bond between the nucleophile and the catalytic base D291. It is proposed that these interactions contribute to a stabilization of TS. This modeling study provided detailed insight into the mechanism of the GlcNAc transfer catalyzed by GnT-I, which is the first step in the conversion of high mannose oligosaccharides to complex and hybrid N-glycan structures.     

A designed synthetic analogue of 4-OT is specific for a non-natural substrate

The substrate specificity of 4-oxalocrotonate tautomerase (4-OT) is characterized by electrostatic interactions between positively charged arginine (Arg) side chains on the enzyme and the dianionic substrate, 4-oxalocrotonate. To generate specific hydrogen-bonding interactions with a monoanionic substrate analogue, we have introduced a urea functional group into the active site by replacing arginine side chains with isosteric citrulline (Cit) residues. This design was based on the complementarity between the urea functionality of citrulline and the uncharged amide function of the substrate, as opposed to the guanidinium-carboxylate electrostatic interaction between the wild-type enzyme and the natural substrate. Indeed, the synthetic (Arg39Cit)4-OT analogue catalyzed the tautomerization of the non-natural monoamide-monoacid substrate while it was a poor catalyst for the natural diacid substrate. The specificity of (Arg39Cit)4-OT for the monoamide-monoacid substrate relative to that of the diacid substrate was found to be 740-fold greater than that of the wild-type enzyme for tautomerization of the non-natural substrate as compared with the natural one. The role of electrostatic interactions in the tautomerization of the monoamide-monoacid substrate was probed in detail with several other Arg to Cit analogues of this enzyme. This study has demonstrated that chemical manipulation of the functional groups within the active site of an enzyme can modify its catalytic activity and substrate specificity in a predictable way, suggesting that the incorporation of noncoded amino acids into proteins has great promise for the development of new enzymatic mechanisms and new binding interactions.     

Role of Mg2+ in hammerhead ribozyme catalysis from molecular simulation

Molecular dynamics simulations have been performed to investigate the role of Mg2+ in the full-length hammerhead ribozyme cleavage reaction. In particular, the aim of this work is to characterize the binding mode and conformational events that give rise to catalytically active conformations and stabilization of the transition state. Toward this end, a series of eight 12 ns molecular dynamics simulations have been performed with different divalent metal binding occupations for the reactant, early and late transition state using recently developed force field parameters for metal ions and reactive intermediates in RNA catalysis. In addition, hybrid QM/MM calculations of the early and late transition state were performed to study the proton-transfer step in general acid catalysis that is facilitated by the catalytic Mg2+ ion. The simulations suggest that Mg2+ is profoundly involved in the hammerhead ribozyme mechanism both at structural and catalytic levels. Binding of Mg2+ in the active site plays a key structural role in the stabilization of stem I and II and to facilitate formation of near attack conformations and interactions between the nucleophile and G12, the implicated general base catalyst. In the transition state, Mg2+ binds in a bridging position where it stabilizes the accumulated charge of the leaving group while interacting with the 2'OH of G8, the implicated general acid catalyst. The QM/MM simulations provide support that, in the late transition state, the 2'OH of G8 can transfer a proton to the leaving group while directly coordinating the bridging Mg2+ ion. The present study provides evidence for the role of Mg2+ in hammerhead ribozyme catalysis. The proposed simulation model reconciles the interpretation of available experimental structural and biochemical data, and provides a starting point for more detailed investigation of the chemical reaction path with combined QM/MM methods.     

Hybrid QM/MM and DFT investigations of the catalytic mechanism and inhibition of the dinuclear zinc metallo-beta-lactamase CcrA from Bacteroides fragilis

Based on hybrid QM/MM molecular dynamics simulation and density functional theoretical (DFT) calculations, we investigate the mechanistic and energetic features of the catalytic action of dizinc metallo-beta-lactamase CcrA from Bacteroides fragilis. The 200 ps QM/MM simulation of the CcrA enzyme in complex with nitrocefin shows that the substrate beta-lactam moiety is directed toward the active site dizinc center through the interactions of aminocarbonyl and carboxylate groups with the two active site zinc ions and the two conserved residues, Lys167 and Asn176. From the determination of the potential energy profile of a relevant enzymatic reaction model, it is found that the nucleophilic displacement reaction step proceeds with a low-barrier height, leading to the formation of an energetically favored reaction intermediate. The results also show that the high catalytic activity of the CcrA enzyme stems from a simultaneous operation of three catalytic components: activation of the bridging hydroxide nucleophile by zinc-coordinated Asp86; polarization of the substrate aminocarbonyl group by the first zinc ion; stabilization of the negative charge developed on the departing amide nitrogen by the second zinc ion. Consistent with the previous experimental finding that the proton-transfer reaction step is rate-limiting, the activation energy of the second step is found to be 1.6 kcal/mol higher than that of the first step. Finally, through an examination of the structural and energetic features of binding of a thiazolidinecarboxylic acid inhibitor to the active site dizinc center, a two-step inhibition mechanism involving a protonation-induced ligand exchange reaction is proposed for the inhibitory action of a tight-binding inhibitor possessing a thiol group.     

Quantum chemical study on the coordination environment of the catalytic zinc ion in matrix metalloproteinases

X-ray analyses of matrix metalloproteinases (MMPs) have shown that the catalytic zinc ion (Zn1) can bind to one to three water molecules in addition to three conserved histidine residues. To estimate the relative stability of the possible Zn1 coordination structures in the active site of the MMPs, we carry out computational analyses on the coordination environment of the Zn1 ion in the gelatinase A enzyme (or matrix metalloproteinase 2; MMP-2). Four-, five-, and six-coordinated complexes representative of the Zn1 site are fully characterized by means of quantum mechanical (QM) methodologies. On one hand, B3LYP/LACVP* minimizations of various cluster models of the MMP-2 active site show that the trigonal bipyramidal geometry is energetically favored in the gas phase and that continuum solvent effects stabilize preferentially the tetrahedral complexes. On the other hand, B3LYP/OPLS-AA hybrid QM/molecular mechanical calculations in the solvated catalytic domain of the MMP-2 enzyme complemented with electrostatic Poisson-Boltzmann calculations show that the mature enzyme presents most likely a Zn1 ion coordinated by three histidine residues and two water molecules, while the active site glutamic acid is negatively charged. In consonance with X-ray diffraction data, other possible Zn1 configurations, a six-coordinated structure with Zn1-water as well as four- and five-coordinated complexes with a Zn1-bound hydroxide, are predicted to be very close in energy.     

Long‐Range Electrostatics‐Induced Two‐Proton Transfer Captured by Neutron Crystallography in an Enzyme Catalytic Site

Neutron crystallography was used to directly locate two protons before and after a pH‐induced two‐proton transfer between catalytic aspartic acid residues and the hydroxy group of the bound clinical drug darunavir, located in the catalytic site of enzyme HIV‐1 protease. The two‐proton transfer is triggered by electrostatic effects arising from protonation state changes of surface residues far from the active site. The mechanism and pH effect are supported by quantum mechanics/molecular mechanics (QM/MM) calculations. The low‐pH proton configuration in the catalytic site is deemed critical for the catalytic action of this enzyme and may apply more generally to other aspartic proteases. Neutrons therefore represent a superb probe to obtain structural details for proton transfer reactions in biological systems at a truly atomic level.     

Aspects of the mechanism of catalysis of glucose oxidase: A docking, molecular mechanics and quantum chemical study

The complex structure of glucose oxidase (GOX) with the substrate glucose was determined using a docking algorithm and subsequent molecular dynamics simulations. Semiempirical quantum chemical calculations were used to investigate the role of the enzyme and FAD co-enzyme in the catalytic oxidation of glucose. On the basis of a small active site model, substrate binding residues were determined and heats of formation were computed for the enzyme substrate complex and different potential products of the reductive half reaction. The influence of the protein environment on the active site model was estimated with a point charge model using a mixed QM/MM method. Solvent effects were estimated with a continuum model. Possible modes of action are presented in relation to experimental data and discussed with respect to related enzymes. The calculations indicate that the redox reaction of GOX differs from the corresponding reaction of free flavins as a consequence of the protein environment. One of the active site histidines is involved in substrate binding and stabilization of potential intermediates, whereas the second histidine is a proton acceptor. The former one, being conserved in a series of oxidoreductases, is also involved in the stabilization of a C4a-hydroperoxy dihydroflavin in the course of the oxidative half reaction.     

Physical nature of interactions within the active site of cytosine-5-methyltransferase

The physical nature of interactions within the active site of cytosine-5-methyltransferase (CMT) was studied using a variation-perturbation energy decomposition scheme defining a sequence of approximate intermolecular interaction energy models. These models have been used to analyze the catalytic activity of residues constituting cytosine-5-methyltransferase active site as well their role in the binding group of de novo designed inhibitors. Our results indicate that Glu119, Arg163, and Arg165 appear to play the dominant role in stabilizing the protonated transition state structure and their influence can be qualitatively approximated by electrostatic interactions alone. The stabilization of neutral structures of the alternative reaction pathway is small, which might suggest the protonated pathway as preferred by the enzyme. Exchange and delocalization terms are negligible in most cases, or they cancel each other to some extent. Interactions of inhibitors with the CMT active site are dominated by electrostatic multipole contributions in analogy with previously studied transition state analogue inhibitors of leucyl aminopeptidase.     

On the origin of the stabilization of the zwitterionic resting state of cysteine proteases: a theoretical study

Papain-like cysteine proteases are ubiquitous proteolytic enzymes. The protonated His199/deprotonated Cys29 ion pair (cathepsin B numbering) in the active site is essential for their proper functioning. The presence of this ion pair stands in contrast to the corresponding intrinsic residue p K a values, indicating a strong influence of the enzyme environment. In the present work we show by molecular dynamics simulations on quantum mechanical/molecular mechanical (QM/MM) potentials that the ion pair is stabilized by a complex hydrogen bond network which comprises several amino acids situated in the active site of the enzyme and 2-4 water molecules. QM/MM reaction path computations for the proton transfer from His199 to the thiolate of the Cys29 moiety indicate that the ion pair is about 32-36 kJ mol (-1) more stable than the neutral form if the whole hydrogen bonding network is active. Without any hydrogen bonding network the ion pair is predicted to be significantly less stable than the neutral form. QM/MM charge deletion analysis and QM model calculations are used to quantify the stabilizing effect of the active-site residues and the L1 helix in favor of the zwitterionic form. The active-site water molecules contribute about 30 kJ mol (-1) to the overall stabilization. Disruption of the hydrogen bonding network upon substrate binding is expected to enhance the nucleophilic reactivity of the thiolate.     

Substrate orientation in 4-oxalocrotonate tautomerase and its effect on QM/MM energy profiles

The tautomerization of 2-oxo-4E-hexendioate by 4-oxalocrotonate tautomerase has been studied by quantum mechanical/molecular mechanical (QM/MM) methods using three models, A-C, with different substrate orientations. The computed QM/MM energy profiles are rather different. Various energy partitioning analyses indicate the origin of these differences and the role of the active site residues for different substrate orientations. The proposed new model C is preferred over the previously used models A and B because it combines favorable substrate binding geometries with reasonable barriers and is consistent with the experimental evidence from mutation studies concerning the catalytic ability of specific residues in the binding site, especially R11'.