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.

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.

Computational Design of Non-natural Sugar Alcohols to Increase Thermal Storage Density: Beyond Existing Organic Phase Change Materials.

Thermal storage, a technology that enables us to control thermal energy, makes it possible to reuse a huge amount of waste heat, and materials with the ability to treat larger thermal energy are in high demand for energy-saving societies. Sugar alcohols are now one promising candidate for phase change materials (PCMs) because of their large thermal storage density. In this study, we computationally design experimentally unknown non-natural sugar alcohols and predict their thermal storage density as a basic step toward the development of new high performance PCMs. The non-natural sugar alcohol molecules are constructed in silico in accordance with the previously suggested molecular design guidelines: linear elongation of a carbon backbone, separated distribution of OH groups, and even numbers of carbon atoms. Their crystal structures are then predicted using the random search method and first-principles calculations. Our molecular simulation results clearly demonstrate that the non-natural sugar alcohols have potential ability to have thermal storage density up to ∼450-500 kJ/kg, which is significantly larger than the maximum thermal storage density of the present known organic PCMs (∼350 kJ/kg). This computational study suggests that, even in the case of H-bonded molecular crystals where the electrostatic energy contributes mainly to thermal storage density, the molecular distortion and van der Waals energies are also important factors to increase thermal storage density. In addition, the comparison between the three eight-carbon non-natural sugar alcohol isomers indicates that the selection of preferable isomers is also essential for large thermal storage density.

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.

Short hydrogen bond between redox-active tyrosine Y(Z) and D1-His190 in the photosystem II crystal structure.

The crystal structure of photosystem II (PSII) analyzed at a resolution of 1.9 Å revealed a remarkably short H-bond between redox-active tyrosine Y(Z) and D1-His190 (2.46 Å donor-acceptor distance). Using large-scale quantum mechanical/molecular mechanical (QM/MM) calculations with the explicit PSII protein environment, we were able to reproduce this remarkably short H-bond in the original geometry of the crystal structure in the neutral [Y(Z)O···H···N(ε)-His-N(δ)H···O═Asn] state, but not in the oxidized states, indicating that the neutral state was the one observed in the crystal structure. In addition to the appropriate redox/protonation state of Y(Z) and D1-His190, we found that the presence of a cluster of water molecules played a key role in shortening the distance between Y(Z) and D1-His190. The orientations of the water molecules in the cluster were energetically stabilized by the highly polarized PSII protein environment, where the Ca ion of the oxygen-evolving complex (OEC) and the OEC ligand D1-Glu189 were also involved.

Distribution of the cationic state over the chlorophyll pair of the photosystem II reaction center.

The reaction center chlorophylls a (Chla) of photosystem II (PSII) are composed of six Chla molecules including the special pair Chla P(D1)/P(D2) harbored by the D1/D2 heterodimer. They serve as the ultimate electron abstractors for water oxidation in the oxygen-evolving Mn(4)CaO(5) cluster. Using the PSII crystal structure analyzed at 1.9 Šresolution, the redox potentials of P(D1)/P(D2) for one-electron oxidation (E(m)) were calculated by considering all PSII subunits and the protonation pattern of all titratable residues. The E(m)(Chla) values were calculated to be 1015-1132 mV for P(D1) and 1141-1201 mV for P(D2), depending on the protonation state of the Mn(4)CaO(5) cluster. The results showed that E(m)(P(D1)) was lower than E(m)(P(D2)), favoring localization of the charge of the cationic state more on P(D1). The P(D1)(•+)/P(D2)(•+) charge ratio determined by the large-scale QM/MM calculations with the explicit PSII protein environment yielded a P(D1)(•+)/P(D2)(•+) ratio of ~80/~20, which was found to be due to the asymmetry in electrostatic characters of several conserved D1/D2 residue pairs that cause the E(m)(P(D1))/E(m)(P(D2)) difference, e.g., D1-Asn181/D2-Arg180, D1-Asn298/D2-Arg294, D1-Asp61/D2-His61, D1-Glu189/D2-Phe188, and D1-Asp170/D2-Phe169. The larger P(D1)(•+) population than P(D2)(•+) appears to be an inevitable fate of the intact PSII that possesses water oxidation activity.

Effects of point mutation on enzymatic activity: correlation between protein electronic structure and motion in chorismate mutase reaction.

Assignment of particular roles to catalytic residues is an important requirement in clearly understanding enzyme functions. Therefore, predicting the catalytic activities of mutant variants is a fundamental challenge in computational biochemistry. Although site-directed mutagenesis is widely used for studying enzymatic activities and other important classes of protein function, interpreting mutation experiments is usually difficult mainly due to side effects induced by point mutations. Because steric and, in many cases, electrostatic effects may affect the local, fine geometries conserved in wild-type proteins that are usually believed to be thermodynamically stable, simply reducing a loss in catalytic activity into clear elements is difficult. To address these important but difficult issues, we performed a systematic ab initio QM/MM computational analysis combined with MD-FEP simulations and all-electron QM calculations for the entire protein matrix. We selected chorismate mutase, one of the simplest and well-known enzymes, to discuss the details of mutational effects on the enzymatic reaction process. On the basis of the reliable free energy profiles of the wild-type enzyme and several mutant variants, we analyzed the effects of point mutations relative to electronic structure and protein dynamics. In general, changes in geometrical parameters introduced by a mutation were usually limited to the local mutational site. However, this local structural modification could affect the global protein dynamics through correlated motions of particular amino acid residues even far from the mutation site. Even for mutant reactions with low catalytic activity, transition state stabilization was observed as a result of conformational modifications and reorganization around the active site. As for the electrostatic effect created by the polar protein environment, the wild-type enzyme was most effectively designed to stabilize the transition state of the reactive substrate, and the effect of global polarization in the electronic structure was found to be a small catalytic element during the process. As electrostatic media for optimum catalysis, both wild-type and mutant variant proteins were generally robust against external electrostatic perturbations. Protein structures have a certain flexibility, which allows them to slightly modulate their conformations to maximize the transition state stabilization in response to the steric perturbations induced by mutations.

Fragment molecular orbital study of the electronic excitations in the photosynthetic reaction center of Blastochloris viridis.

All electron calculations were performed on the photosynthetic reaction center of Blastochloris viridis, using the fragment molecular orbital (FMO) method. The protein complex of 20,581 atoms and 77,754 electrons was divided into 1398 fragments, and the two-body expansion of FMO/6-31G* was applied to calculate the ground state. The excited electronic states of the embedded electron transfer system were separately calculated by the configuration interaction singles approach with the multilayer FMO method. Despite the structural symmetry of the system, asymmetric excitation energies were observed, especially on the bacteriopheophytin molecules. The asymmetry was attributed to electrostatic interaction with the surrounding proteins, in which the cytoplasmic side plays a major role.

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

In this study, we investigated the electronic character of protein environment in enzymatic processes by performing all-electron QM calculations based on the fragment molecular orbital (FMO) method. By introducing a new computational strategy combining all-electron QM analysis with ab initio QM/MM modeling, we investigated the details of molecular interaction energy between a reactive substrate and amino acid residues at a catalytic site. For a practical application, we selected the chorismate mutase catalyzed reaction as an example. Because the computational time required to perform all-electron QM reaction path searches was very large, we employed the ab initio QM/MM modeling technique to construct reliable reaction profiles and performed all-electron FMO calculations for the selected geometries. The main focus of the paper is to analyze the details of electrostatic stabilization, which is considered to be the major feature of enzymatic catalyses, and to clarify how the electronic structure of proteins is polarized in response to the change in electron distribution of the substrate. By performing interaction energy decomposition analysis from a quantum chemical viewpoint, we clarified the relationship between the location of amino acid residues on the protein domain and the degree of electronic polarization of each residue. In particular, in the enzymatic transition state, Arg7, Glu78, and Arg90 are highly polarized in response to the delocalized electronic character of the substrate, and as a result, a large amount of electrostatic stabilization energy is stored in the molecular interaction between the enzyme and the substrate and supplied for transition state stabilization.

Molecular Mechanism of the Reaction Specificity in Threonine Synthase: Importance of the Substrate Conformations.

Threonine synthase (ThrS) catalyzes the final chemical reaction of l-threonine biosynthesis from its precursor, O-phospho-l-homoserine. As the phosphate ion generated in its former half reaction assists its latter reaction, ThrS is recognized as one of the best examples of product-assisted catalysis. In our previous QM/MM study, the chemical reactions for the latter half reactions, which are critical for the product-assisted catalysis, were revealed. However, accurate free energy changes caused by the conformational ensembles and entrance of water molecules into the active site are unknown. In the present study, by performing long-time scale MD simulations, the free energy changes by the divalent anions (phosphate or sulfate ions) and conformational states of the intermediate states were theoretically investigated. We found that the calculated free energy double differences are in good agreement with the experimental results. We also revealed that the phosphate ion contributes to forming hydrogen bonds that are suitable for the main reaction progress. This means that the conformation of the active site amino acid residues and the substrate, and hence, the tunable catalysis, are controlled by the product phosphate ion, and this clearly demonstrates a molecular mechanism of the product-assisted catalysis in ThrS.

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.

Molecular mechanisms of substrate specificities of uridinecytidine kinase.

A uridine-cytidine kinase (UCK) catalyzes the phosphorylation of uridine (Urd) and cytidine (Cyd) and plays a significant role in the pyrimidine-nucleotide salvage pathway. Unlike ordinary ones, UCK from Thermus thermophilus HB8 (ttCK) loses catalytic activity on Urd due to lack of a substrate binding ability and possesses an unusual amino acid, i.e. tyrosine 93 (Tyr93) at the binding site, whereas histidine (His) is located in the other UCKs. Mutagenesis experiments revealed that a replacement of Tyr93 by His or glutamine (Gln) recovered catalytic activity on Urd. However, the detailed molecular mechanism of the substrate specificity has remained unclear. In the present study, we performed molecular dynamics simulations on the wild-type ttCK, two mutant ttCKs, and a human UCK bound to Cyd and three protonation forms of Urd to elucidate their substrate specificity. We found three residues, Tyr88, Tyr/His/Gln93 and Arg152 in ttCKs, are important for recognizing the substrates. Arg152 contributes to induce a closed form of the binding site to retain the substrate, and the N3 atom of Urd needed to be deprotonated. Although Tyr88 tightly bound Cyd, it did not sufficiently bind Urd because of lack of the hydrogen bonding. His/Gln93 complemented the interaction of Tyr88 and raised the affinity of ttCK to Urd. The crucial distinction between Tyr and His or Gln was a role in the hydrogen-bonding network. Therefore, the ability to form both hydrogen-bonding donor and accepter is required to bind both Urd and Cyd.

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.

Change in a protein's electronic structure induced by an explicit solvent: an ab initio fragment molecular orbital study of ubiquitin.

The effect of solvation on the electronic structure of the ubiquitin protein was analyzed using the ab initio fragment molecular orbital (FMO) method. FMO calculations were performed for the protein in vacuo, and the protein was immersed in an explicit solvent shell as thick as 12 A at the HF or MP2 level by using the 6-31G* basis set. The protein's physical properties examined were the net charge, the dipole moment, the internal energy, and the solvent interaction energy. Comparison of the computational results revealed the following changes in the protein upon solvation. First, the positively charged amino acid residues on the protein surface drew electrons from the solvent, while the negatively charged ones transfer electrons to the solvent. Second, the dipole moment of the protein was enhanced as a result of the polarization. Third, the internal energy of the protein was destabilized, but the destabilization was more than compensated for by the generation of a favorable protein-solvent interaction. Finally, the energetic changes were elicited both by the electron correlation effect of the first solvent shell and by the electrostatic effect of more distant solvent molecules. These findings were consistent with the picture of the solvated protein being a polarizable molecule dissolved in a dielectric media.

The fragment molecular orbital method for geometry optimizations of polypeptides and proteins.

The fragment molecular orbital method (FMO) has been used with a large number of wave functions for single-point calculations, and its high accuracy in comparison to ab initio methods has been well established. We have developed the analytic derivative of the electrostatic interaction between far separated fragments and performed a number of restricted Hartree-Fock (RHF) geometry optimizations using FMO and ab initio methods. In particular, the alpha-helix, beta-turn, and extended conformers of a 10-residue polyalanine were studied and the good FMO accuracy was established (the rms deviations for the former two forms were about 0.2 A and for the latter structure about 0.001 A). Met-enkephalin dimer was used as a model for the polypeptide binding and computed at the 3-21G and 6-31G* levels with a similar accuracy achieved; the error in the binding energy predictions (FMO vs ab initio) was 1-3 kcal/mol. Chignolin (PDB: 1uao) and an agonist polypeptide of the erythropoietin receptor protein (emp1) were optimized at the 3-21(+)G level, with the rms deviation from ab initio of about 0.2 A, or 0.5 degrees in terms of bond angles. The effect of solvation on the structure optimization was studied in chignolin and the Trp-cage miniprotein construct (PDB:1l2y), by describing water with TIP3P. The computed structures in gas phase and solution are compared to each other and experiment.

Accuracy of the three-body fragment molecular orbital method applied to Møller-Plesset perturbation theory.

The three-body energy expansion in the fragment molecular orbital method (FMO) was applied to the 2nd order Møller-Plesset theory (MP2). The accuracy of both the two and three-body expansions was determined for water clusters, alanine n-mers (alpha-helices and beta-strands) and one synthetic protein, using the 6-31G* and 6-311G* basis sets. At the best level of theory (three-body, two molecules/residues per fragment), the absolute errors in energy relative to ab initio MP2 were at most 1.2 and 5.0 mhartree, for the 6-31G* and 6-311G* basis sets, respectively. The relative accuracy was at worst 99.996% and 99.96%, for 6-31G* and 6-311G*, respectively. A three-body approximation was introduced and the optimum threshold value was determined. The protein calculation (6-31G*) at the production level (FMO2/2) took 3 h on 36 3.2-GHz Pentium 4 nodes and had the absolute error in the MP2 correlation energy of only 2 kcal/mol.

Low-barrier hydrogen bond hypothesis in the catalytic triad residue of serine proteases: correlation between structural rearrangement and chemical shifts in the acylation process.

To elucidate the catalytic advantage of the low-barrier hydrogen bond (LBHB), we analyze the hydrogen bonding network of the catalytic triad (His57-Asp102-Ser195) of serine protease trypsin, one of the best examples of the LBHB reaction mechanism. Especially, we focus on the correlation between the change of the chemical shifts and the structural rearrangement of the active site in the acylation process. To clarify LBHB, we evaluate the two complementary properties. First, we calculate the NMR chemical shifts of the imidazole ring of His57 by the gauge-including atomic orbital (GIAO) approach within the ab initio QM/MM framework. Second, the free energy profile of the proton transfer from His57 to Asp102 in the tetrahedral intermediate is obtained by ab initio QM/MM calculations combined with molecular dynamics free energy perturbation (MD-FEP) simulations. The present analyses reveal that the calculated shifts reasonably reproduce the observed values for (1)H chemical shift of H(epsilon)(1) and H(delta)(1) in His57. The (15)N and (13)C chemical shifts are also consistent with the experiments. It is also shown that the proton between His57 and Asp102 is localized at the His57 side. This largely downfield chemical shift is originated from the strong electrostatic interaction, not a covalent-like bonding character between His57 and Asp102. Also, it is proved that a slight downfield character of H(epsilon)(1) is originated from a electrostatic interaction between His57 and the backbone carbonyl group of Val213 and Ser214. These downfield chemical shifts are observed only when the tetrahedral intermediate is formed in the acylation process.

Multilayer formulation of the fragment molecular orbital method (FMO).

The fragment molecular orbital method (FMO) has been generalized to allow for multilayer structure. Fragments are assigned to layers, and each layer can be described with a different basis set and/or level of electron correlation. Interlayer boundaries are treated in the general spirit of the FMO method since they also coincide with some interfragment boundaries. The question of the one- and two-layer FMO accuracy dependence upon the fragmentation scheme is also addressed. The new method has been applied to predict the reaction barrier and the reaction heat for the Diels-Alder reaction with a representative set of reactants based on dividing fragments in two layers. The 6-31G* basis set has been used for the active site and the 6-31G*, 6-31G, 3-21G, and STO-3G basis sets have been used for the substituents. Different levels of electron correlation (RHF, B3LYP, and MP2) have been applied to layers in systematic fashion. The one-layer FMO errors in the reaction barrier and the reaction heat were 2.0 kcal/mol or less for all levels applied (RHF, B3LYP, and MP2), relative to full ab initio methods. For the two-layer method the error was found to be several kcal/mol. Benchmark calculations of the activation barrier for the decarboxylation of phenylcyanoacetate by beta-cyclodextrin demonstrated that the two-layer calculations are efficient, being 36 times faster than the regular DFT, as well as accurate, with the error being 1.0 kcal/mol.

Theoretical perspectives on the reaction mechanism of serine proteases: the reaction free energy profiles of the acylation process.

The reaction mechanism of serine proteases (trypsin), which catalyze peptide hydrolysis, is studied theoretically by ab initio QM/MM electronic structure calculations combined with Molecular Dynamics-Free Energy Perturbation calculations. We have calculated the entire reaction free energy profiles of the first reaction step of this enzyme (acylation process). The present calculations show that the rate-determining step of the acylation is the formation of the tetrahedral intermediate, and the breakdown of this intermediate has a small energy barrier. The calculated activation free energy for the acylation is approximately 17.8 kcal/mol at QM/MM MP2/(aug)-cc-pVDZ//HF/6-31(+)G/AMBER level, and this reaction is an exothermic process. MD simulations of the enzyme-substrate (ES) complex and the free enzyme in aqueous phase show that the substrate binding induces slight conformational changes around the active site, which favor the alignment of the reactive fragments (His57, Asp102, and Ser195) together in a reactive orientation. It is also shown that the proton transfer from Ser195 to His57 and the nucleophilic attack of Ser195 to the carbonyl carbon of the scissile bond of the substrate occur in a concerted manner. In this reaction, protein environment plays a crucial role to lowering the activation free energy by stabilizing the tetrahedral intermediate compared to the ES complex. The polarization energy calculations show that the enzyme active site is in a very polar environment because of the polar main chain contributions of protein. Also, the ground-state destabilization effect (steric strain) is not a major catalytic factor. The most important catalytic factor of stabilizing the tetrahedral intermediate is the electrostatic interaction between the active site and particular regions of protein: the main chain NH groups in Gly193 and Ser195 (so-called oxyanion hole region) stabilize negative charge generated on the carbonyl oxygen of the scissile bond, and the main chain carbonyl groups in Ile212 approximately Ser214 stabilize a positive charge generated on the imidazole ring of His57.

Role of Asp102 in the catalytic relay system of serine proteases: a theoretical study.

The role of Asp102 in the catalytic relay system of serine proteases is studied theoretically by calculating the free energy profiles of the single proton-transfer reaction by the Asn102 mutant trypsin and the concerted double proton-transfer reaction (so-called the charge-relay mechanism) of the wild-type trypsin. For each reaction, the reaction free energy profile of the rate-determining step (the tetrahedral intermediate formation step) is calculated by using ab initio QM/MM electronic structure calculations combined with molecular dynamics-free energy perturbation method. In the mutant reaction, the free energy monotonically increases along the reaction path. The rate-determining step of the mutant reaction is the formation of tetrahedral intermediate complex, not the base (His57) abstraction of the proton from Ser195. In contrast to the single proton-transfer reaction of the wild-type, MD simulations of the enzyme-substrate complex show that the catalytically favorable alignment of the relay system (the hydrogen bonding network between the mutant triad, His57, Asn102, and Ser195) is rarely observed even in the presence of a substrate at the active site. In the double proton-transfer reaction, the energy barrier is observed at the proton abstraction step, which corresponds to the rate-determining step of the single proton-transfer reaction of the wild-type. Although both reaction profiles show an increase of the activation barrier by several kcals/mol, these increases have different energetic origins: a large energetic loss of the electrostatic stabilization between His57 and Asn102 in the mutant reaction, while the lack of stabilization by the protein environment in the double proton-transfer reaction. Comparing the present results with the single proton transfer of the wild-type, Asp102 is proven to play two important roles in the catalytic process. One is to stabilize the protonated His57, or ionic intermediate, formed during the acylation, and the other is to fix the configuration around the active site, which is favorable to promote the catalytic process. These two factors are closely related to each other and are indispensable for the efficient catalysis. Also the present calculations suggest the importance of the remote site interaction between His57 and Val213-Ser214 at the catalytic transition state.