Computational studies on nonenzymatic and enzymatic pyridoxal phosphate catalyzed decarboxylations of 2-aminoisobutyrate

A computational study of nonenzymatic and enzymatic pyridoxal phosphate-catalyzed decarboxylation of 2-aminoisobutyrate (AIB) is presented. Four prototropic isomers of a model aldimine between AIB and 5'-deoxypyridoxal, with acetate interacting with the pyridine nitrogen, were employed in calculations of both gas phase and water model (PM3 and PM3-SM3) decarboxylation reaction paths. Calculations employing the transition state structures obtained for the four isomers allow the demonstration of stereoelectronic effects in transition state stabilization as well as a separation of the contributions of the Schiff base and pyridine ring moieties to this stabilization. The unprotonated Schiff base contribution (approximately 16 kcal/mol) is larger than that of the pyridine ring even when it is protonated (approximately 10 kcal/mol), providing an explanation of the catalytic power of pyruvoyl-dependent amino acid decarboxylases. An active site model of dialkylglycine decarboxylase was constructed and validated, and enzymatic decarboxylation reaction paths were calculated. The reaction coordinate is shown to be complex, with proton transfer from Lys272 to the coenzyme C4' likely simultaneous with C alpha--CO(2)(-) bond cleavage. The proposed concerted decarboxylation/proton-transfer mechanism provides a simple explanation for the observed specificity of this enzyme toward oxidative decarboxylation.     

Carbinolamine formation and dehydration in a DNA repair enzyme active site

In order to suggest detailed mechanistic hypotheses for the formation and dehydration of a key carbinolamine intermediate in the T4 pyrimidine dimer glycosylase (T4PDG) reaction, we have investigated these reactions using steered molecular dynamics with a coupled quantum mechanics-molecular mechanics potential (QM/MM). We carried out simulations of DNA abasic site carbinolamine formation with and without a water molecule restrained to remain within the active site quantum region. We recovered potentials of mean force (PMF) from thirty replicate reaction trajectories using Jarzynski averaging. We demonstrated feasible pathways involving water, as well as those independent of water participation. The water-independent enzyme-catalyzed reaction had a bias-corrected Jarzynski-average barrier height of approximately (6.5 kcal mol(-1) (27.2 kJ mol(-1)) for the carbinolamine formation reaction and 44.5 kcal mol(-1) (186 kJ mol(-1)) for the reverse reaction at this level of representation. When the proton transfer was facilitated with an intrinsic quantum water, the barrier height was approximately 15 kcal mol(-1) (62.8 kJ mol(-1)) in the forward (formation) reaction and 19 kcal mol(-1) (79.5 kJ mol(-1)) for the reverse. In addition, two modes of unsteered (free dynamics) carbinolamine dehydration were observed: in one, the quantum water participated as an intermediate proton transfer species, and in the other, the active site protonated glutamate hydrogen was directly transferred to the carbinolamine oxygen. Water-independent unforced proton transfer from the protonated active site glutamate carboxyl to the unprotonated N-terminal amine was also observed. In summary, complex proton transfer events, some involving water intermediates, were studied in QM/MM simulations of T4PDG bound to a DNA abasic site. Imine carbinolamine formation was characterized using steered QM/MM molecular dynamics. Dehydration of the carbinolamine intermediate to form the final imine product was observed in free, unsteered, QM/MM dynamics simulations, as was unforced acid-base transfer between the active site carboxylate and the N-terminal amine.     

QM/MM study of the active site of free papain and of the NMA-papain complex.

Hybrid quantum mechanical/molecular mechanical (QM/MM) calculations using restricted and unrestricted Hartree-Fock and B3LYP ab initio (QM) and Amber force field (MM), respectively, have been applied to study the catalytic site of papain in both free and substrate bonded forms. Ab initio geometry optimizations have been performed for the active site of papain and the N-methyl-acetamide (NMA)-papain complex within the molecular mechanical treatment of the protein environment. A covalent tetrahedral intermediate structure could be obtained only when the amide N atom of the substrate molecule was protonated through a proton transfer from the His-159 in the catalytic site. Our results support the previous assumption that a proton transfer from His-159 to the amide N atom of the substrate occurs prior to or concerted with the nucleophilic attack of the Cys-25 sulfur atom to the carbonyl group of the substrate. The electron correlation effect will reduce the proton transfer barrier. Therefore, this proton transfer can be easily observed in the B3LYP/6-31G* calculations. The HF/6-31G* method overestimates the reaction barrier against this proton transfer. The sulfur atom of Cys-25 and the imidazole ring of His-159 are found to be coplanar in the free form of the enzyme. However, the rotation of the imidazole ring of His-159 was observed during the formation of the tetrahedral intermediate. Without the papain environment, the coplanar thiolate-imidazolium ion pair RS-...ImH+ is much less stable than the neutral form of RSH....Im. Within the protein environment, however, the thiolate-imidazolium ion pair becomes more stable than its neutral form by 4.1 and 0.4 kcal/mol in HF/6-31G* and B3LYP/6-31G* calculations, respectively. The barrier of proton transfer from S-H group of Cys-25 to the imidazole ring of His-159 was reduced from 22.0 kcal/mol to 15.2 kcal/mol by the protein environment in HF/6-31G* calculations. This barrier is found to be much smaller (2.5 kcal/mol) in B3LYP/6-31G* calculations.     

Catalytic mechanism of cellulose degradation by a cellobiohydrolase, CelS

The hydrolysis of cellulose is the bottleneck in cellulosic ethanol production. The cellobiohydrolase CelS from Clostridium thermocellum catalyzes the hydrolysis of cello-oligosaccharides via inversion of the anomeric carbon. Here, to examine key features of the CelS-catalyzed reaction, QM/MM (SCCDFTB/MM) simulations are performed. The calculated free energy profile for the reaction possesses a 19 kcal/mol barrier. The results confirm the role of active site residue Glu87 as the general acid catalyst in the cleavage reaction and show that Asp255 may act as the general base. A feasible position in the reactant state of the water molecule responsible for nucleophilic attack is identified. Sugar ring distortion as the reaction progresses is quantified. The results provide a computational approach that may complement the experimental design of more efficient enzymes for biofuel production.     

Protonation states of intermediates in the reaction mechanism of [NiFe] hydrogenase studied by computational methods

The [NiFe] hydrogenases catalyse the reversible conversion of H₂ to protons and electrons. The active site consists of a Fe ion with one carbon monoxide, two cyanide, and two cysteine (Cys) ligands. The latter two bridge to a Ni ion, which has two additional terminal Cys ligands. It has been suggested that one of the Cys residues is protonated during the reaction mechanism. We have used combined quantum mechanical and molecular mechanics (QM/MM) geometry optimisations, large QM calculations with 817 atoms, and QM/MM free energy simulations, using the TPSS and B3LYP methods with basis sets extrapolated to the quadruple zeta level to determine which of the four Cys residues is more favourable to protonate for four putative states in the reaction mechanism, Ni-SI_a, Ni-R, Ni-C, and Ni-L. The calculations show that for all states, the terminal Cys-546 residue is most easily protonated by 14–51 kJ/mol, owing to a more favourable hydrogen-bond pattern around this residue in the protein.     

Structure of a Michaelis complex analogue: propionate binds in the substrate carboxylate site of alanine racemase

The structure of alanine racemase from Bacillus stearothermophilus with the inhibitor propionate bound in the active site was determined by X-ray crystallography to a resolution of 1.9 A. The enzyme is a homodimer in solution and crystallizes with a dimer in the asymmetric unit. Both active sites contain a pyridoxal 5'-phosphate (PLP) molecule in aldimine linkage to Lys39 as a protonated Schiff base, and the pH-independence of UV-visible absorption spectra suggests that the protonated PLP-Lys39 Schiff base is the reactive form of the enzyme. The carboxylate group of propionate bound in the active site makes numerous interactions with active-site residues, defining the substrate binding site of the enzyme. The propionate-bound structure therefore approximates features of the Michaelis complex formed between alanine racemase and its amino acid substrate. The structure also provides evidence for the existence of a carbamate formed on the side-chain amino group of Lys129, stabilized by interactions with one of the residues interacting with the carboxylate group of propionate, Arg136. We propose that this novel interaction influences both substrate binding and catalysis by precisely positioning Arg136 and modulating its charge.     

Determinants of reactivity and selectivity in soluble epoxide hydrolase from quantum mechanics/molecular mechanics modeling

Soluble epoxide hydrolase (sEH) is an enzyme involved in drug metabolism that catalyzes the hydrolysis of epoxides to form their corresponding diols. sEH has a broad substrate range and shows high regio- and enantioselectivity for nucleophilic ring opening by Asp333. Epoxide hydrolases therefore have potential synthetic applications. We have used combined quantum mechanics/molecular mechanics (QM/MM) umbrella sampling molecular dynamics (MD) simulations (at the AM1/CHARMM22 level) and high-level ab initio (SCS-MP2) QM/MM calculations to analyze the reactions, and determinants of selectivity, for two substrates: trans-stilbene oxide (t-SO) and trans-diphenylpropene oxide (t-DPPO). The calculated free energy barriers from the QM/MM (AM1/CHARMM22) umbrella sampling MD simulations show a lower barrier for phenyl attack in t-DPPO, compared with that for benzylic attack, in agreement with experiment. Activation barriers in agreement with experimental rate constants are obtained only with the highest level of QM theory (SCS-MP2) used. Our results show that the selectivity of the ring-opening reaction is influenced by several factors, including proximity to the nucleophile, electronic stabilization of the transition state, and hydrogen bonding to two active site tyrosine residues. The protonation state of His523 during nucleophilic attack has also been investigated, and our results show that the protonated form is most consistent with experimental findings. The work presented here illustrates how determinants of selectivity can be identified from QM/MM simulations. These insights may also provide useful information for the design of novel catalysts for use in the synthesis of enantiopure compounds.     

Multiple hydrogen kinetic isotope effects for enzymes catalyzing exchange with solvent: application to alanine racemase

Alanine racemase catalyzes the pyridoxal phosphate-dependent interconversion of the D- and L-isomers of alanine. Previous studies have shown that the enzyme employs a two-base mechanism in which Lys39 and Tyr265 are the acid/base catalysts. It is thus possible that stereoisomerization of the external aldimine intermediates occurs through a concerted double proton transfer without the existence of a distinct carbanionic intermediate. This possibility was tested by the application of multiple kinetic isotope effect (KIE) methodology to alanine racemase. The mutual dependence of primary substrate and solvent deuterium KIEs has been measured using equilibrium perturbation-type experiments. The conceptually straightforward measurement of the substrate KIE in H(2)O is complemented with a less intuitive protium washout perturbation-type measurement in D(2)O. The primary substrate KIE in the D --> L direction at 25 degrees C is reduced from 1.297 in H(2)O to 1.176 in D(2)O, while in the L --> D direction it is reduced from 1.877 in H(2)O to 1.824 in D(2)O. Similar reductions are also observed at 65 degrees C, the temperature to which the Bacillus stearothermophilus enzyme is adapted. These data strongly support a stepwise racemization of stereoisomeric aldimine intermediates in which a substrate-based carbanion is an obligatory intermediate. The ionizations observed in k(cat)/K(M) pH profiles have been definitively assigned based on the DeltaH(ion) values of the observed pK(a)'s with alanine and on the pH dependence of k(cat)/K(M) for the alternative substrate serine. The acidic pK(a) in the bell-shaped curve is due to the phenolic hydroxyl of Tyr265, which must be unprotonated for reaction with either isomer of alanine. The basic pK(a) is due to the substrate amino group, which must be protonated to react with Tyr265-unprotonated enzyme. A detailed reaction mechanism incorporating these results is proposed.     

A quantum mechanics/molecular mechanics study of the catalytic mechanism and product specificity of viral histone lysine methyltransferase

There are three reaction steps in the S-adenosylmethionine (AdoMet) methylation of lysine-NH2 catalyzed by a methyltransferase. They are (i) combination of enzyme.Lys-NH3+ with AdoMet, (ii) substrate ionization to provide enzyme.AdoMet.Lys-NH2, and (iii) methyl transfer providing enzyme.AdoHcy.Lys-N(Me)H2+ and the dissociation of AdoHcy. In this study of the viral histone methyltransferase (vSET), we find that substrate ionization of vSET.Lys27-NH3+, vSET.Lys27-N(Me)H2+, and vSET.Lys27-N(Me)2H+ takes place upon combination with AdoMet. The presence of a water channel allows dissociation of a proton to the solvent. There is no water channel in the absence of AdoMet. That the formation of a water channel is combined with AdoMet binding was first discovered in our investigation of Rubisco large subunit methyltransferase. Via a quantum mechanics/molecular mechanics (QM/MM) approach, the calculated free energy barrier (DeltaG++) of the first methyl transfer reaction catalyzed by vSET [Lys27-NH2 + AdoMet --> Lys27-N(Me)H2+ + AdoHcy] equals 22.5 +/- 4.3 kcal/mol, which is in excellent agreement with the free energy barrier (21.7 kcal/mol) calculated from the experimental rate constant (0.047 min-1). The calculated DeltaG++ of the second methyl transfer reaction [AdoMet + Lys27-N(Me)H --> AdoHcy + Lys27-N(Me)2H+] at the QM/MM level is 22.6 +/- 3.6 kcal/mol, which is in agreement with the value of 22.4 kcal/mol determined from the experimental rate constant (0.015 min-1). The third methylation [Lys27-N(Me)2 + AdoMet --> Lys27-N(Me)3+ + AdoHcy] is associated with a DeltaG++ of 23.1 +/- 4.0 kcal/mol, which is in agreement with the value of 23.0 kcal/mol determined from the experimental rate constant (0.005 min-1). Our computations establish that the first, second, and third methyl transfer steps catalyzed by vSET are linear SN2 reactions with the bond making being approximately 50% associative.     

Second step of hydrolytic dehalogenation in haloalkane dehalogenase investigated by QM/MM methods

Mechanistic studies on the hydrolytic dehalogenation catalyzed by haloalkane dehalogenases are of importance for environmental and industrial applications. Here, Car-Parrinello (CP) and ONIOM hybrid quantum-mechanical/molecular mechanics (QM/MM) are used investigate the second reaction step of the catalytic cycle, which comprises a general base-catalyzed hydrolysis of an ester intermediate (EI) to alcohol and free enzyme. We focus on the enzyme LinB from Sphingomonas paucimobilis UT26, for which the X-ray structure at atomic resolution is available. In agreement with previous proposals, our calculations suggest that a histidine residue (His272), polarized by glutamate (Glu132), acts as a base, accepting a proton from the catalytic water molecule and transferring it to an alcoholate ion. The reaction proceeds through a metastable tetrahedral intermediate, which shows an easily reversed reaction to the EI. In the formation of the products, the protonated aspartic acid (Asp108) can easily adopt conformation of the relaxed state found in the free enzyme. The overall free energy barrier of the reaction calculated by potential of the mean force integration using CP-QM/MM calculations is equal to 19.5 +/- 2 kcal . mol(-1). The lowering of the energy barrier of catalyzed reaction with respect to the water reaction is caused by strong stabilization of the reaction intermediate and transition state and their preorganization by electrostatic field of the enzyme.     

A definitive mechanism for chorismate mutase

In previous research presentations, we have described the important features of the chorismate --> prephenate reaction using molecular dynamics (MD) and thermodynamic integration studies. This investigation of the reaction in Escherichia coli and water involves QM/MM procedures (SCCDFTB/MM two-dimensional reaction coordinates to identify transition state structures in the water, enzyme, and gas phase followed by B3LYP/6-31+G* single-point computations which allow the determination of activation energies in water and in the E. coli enzyme). Computed activation energies of 11.3 kcal/mol in enzyme and 20.3 kcal/mol in water may be compared to the experimental values of 12.7 and 20.7 kcal/mol, respectively. The transition state structures in the gas phase, water, and enzyme are much the same. The transition states are characteristic of a concerted pericyclic rearrangement. The very small differences in the partial charges of O13 in NAC and TS support only a small preferential (10%) electrostatic stabilization of TS. The free energy of NAC formation in water exceeds that in enzyme by 8.5 kcal/mol, and it is this favored formation of NAC that provides the major kinetic advantage to the enzymatic reaction. These findings compare most favorably with those previous observations of this laboratory employing molecular dynamics and thermodynamic integrations. A definitive mechanism for the chorismate mutase enzymes is provided.     

The 1.9 A crystal structure of alanine racemase from Mycobacterium tuberculosis contains a conserved entryway into the active site

We report the crystal structure of alanine racemase from Mycobacterium tuberculosis (Alr(Mtb)) at 1.9 A resolution. In our structure, Alr(Mtb) is found to be a dimer formed by two crystallographically different monomers, each comprising 384 residues. The domain makeup of each monomer is similar to that of Bacillus and Pseudomonas alanine racemases and includes both an alpha/beta-barrel at the N-terminus and a C-terminus primarily made of beta-strands. The hinge angle between these two domains is unique for Alr(Mtb), but the active site geometry is conserved. In Alr(Mtb), the PLP cofactor is covalently bound to the protein via an internal aldimine bond with Lys42. No guest substrate is noted in its active site, although some residual electron density is observed in the enzyme's active site pocket. Analysis of the active site pocket, in the context of other known alanine racemases, allows us to propose the inclusion of conserved residues found at the entrance to the binding pocket as additional targets in ongoing structure-aided drug design efforts. Also, as observed in other alanine racemase structures, PLP adopts a conformation that significantly distorts the planarity of the extended conjugated system between the PLP ring and the internal aldimine bond.     

The protein backbone makes important contributions to 4-oxalocrotonate tautomerase enzyme catalysis: understanding from theory and experiment

The role of polypeptide backbone interactions in 4-oxalocrotonate tautomerase (4OT) catalysis has been investigated using a combination of site-directed mutagenesis experiments with unnatural amino acids and quantum mechanical/molecular mechanical (QM/MM) calculations of the 4OT reaction mechanism. Energy barriers for the wild-type enzyme (wt-4OT) and for a 4OT analogue containing a backbone amide to ester bond mutation between Ile-7 and Leu-8 [(OL8)4OT] were determined by both theory and experiment. The amide to ester bond mutation in (OL8)4OT effectively deleted a putative hydrogen bonding interaction between the enzyme's polypeptide backbone and its substrate. Recent theoretical calculations for the 4OT reaction mechanism suggested that this hydrogen bonding interaction helps properly position the substrate in the active site [Cisneros, G. A., et al. (2003) J. Am. Chem. Soc. 125, 10384-10393]. Our experimental results for (OL8)4OT reveal that the energy barrier for the (OL8)4OT-catalyzed reaction was increased 1.8 kcal/mol over that of the wild-type enzyme. This increase was in good agreement with the 1.0 kcal/mol increase obtained from QM/MM calculations for this analogue. Our theoretical calculations further suggest the hydrogen bond deletion in (OL8)4OT results in a rearrangement of the substrate in the active site. In this rearrangement, an ordered water molecule loses its ability to stabilize the transition state (TS), and Arg-61 gains the ability to stabilize the TS. The predicted role of Arg-61 in (OL8)4OT catalysis was confirmed in kinetic experiments with an analogue of (OL8)4OT containing an Arg to Ala mutation at position 61.     

Energy of binding of Aspergillus oryzae beta-glucosidase with the substrate, and the mechanism of its enzymic action

Several beta-D-glucopyranosides (p-nitrophenyl, phenyl, and ethyl), 1-thio-beta-D-glucopyranosides, and phenyl 2-deoxy, 3-deoxy, 4-deoxy, and 6-deoxy beta-D-glucopyranosides were synthesized and used to study the mechanism of the enzymatic action of Taka-beta-glucosidase [EC 3.2.1.21 Aspergillus oryzae]. Kinetic constants of the enzyme for these glycosides were determined from S/V-S or 1/V-1/S plots, and the hydrolysis rates of these compounds with the enzyme, acid (3 N HCl) and alkali (3 N NaOH) were compared. Inhibition of the enzyme by 1,5-anhydroglucitol, glucal, dihydroglucal, and 1,6-anhydroglucopyranose was also examined. Glucal and 1,5-anhydroglucitol showed strong competitive inhibition. Free energy of binding of each hydroxyl group of glucosidic glucose with the enzyme was estimated from Kms of phenyl beta-glucoside and its deoxy analogues, and also Ki values of some inhibitors. The free energies of binding of 2-OH, 3-OH, 4-OH, and 6-OH were calculated to be 1.1, 2.4, 0.7, and 1.8 kcal/mol, respectively. The free energy of binding of phenoxide at C-1 (0.3 kcal/mol) was calculated from the Km of Ph-beta-Glc and Ki of 1,5-anhydroglucitol. The energy of binding of 5-CH2OH (2.3 kcal/mol) was obtained from the Km of Ph-beta-Glc and that of Ph-beta-Xyl. The sum (6.8 kcal/mol) of each partial binding free energy was close to the value of binding free energy of Ph-beta-Glc (7.0 kcal/mol) calculated by the equation; -delta Gbind = -RT ln Km-T delta Smix, showing that the methods of estimation of each binding energy used in the present study seemed reasonable. Glucal, having a pyranose form distorted slightly, showed strong competitive inhibition and the Ki of this inhibitor was smaller than the Km of Ph-beta-Glc, suggesting that the sugar ring bound to the active site was distored to a half chair form which is labile to acid hydrolysis.     

Product specificity and mechanism of protein lysine methyltransferases: insights from the histone lysine methyltransferase SET8

Molecular dynamics simulations employing a molecular mechanics (MM) force field and hybrid quantum mechanics (QM) and MM (QM/MM) have been carried out to investigate the product specificity and mechanism of the histone H4 lysine 20 (H4-K20) methylation by human histone lysine methyltransferase SET8. At neutral pH, the target lysine is available to only the enzyme in the protonated state. The first step in the methylation reaction must be deprotonation of the lysine target which is followed by the (+)AdoMet methylation of the neutral lysine [Enz.Lys-CH(2)-NH(3)(+).(+)AdoMet --> H(+) + Enz.Lys-CH(2)-NH(2).(+)AdoMet -->--> Enz.Lys-CH(2)-N(Me)H(2)(+).AdoHcy]. The electrostatic interactions between two positive charges on (+)AdoMet and Lys20-NH(3)(+) decrease the pK(a) of Lys20-NH(3)(+). Upon formation of Enz.Lys-NH(3)(+).(+)AdoMet, a water channel by which the proton escapes to the outer solvent phase is formed. The formation of a water channel for the escape of a proton from Lys20-N(Me)H(2)(+) in Enz.Lys20-N(Me)H(2)(+).(+)AdoMet is not formed because the methyl substituent blocks the starting of the water channel. Thus, a second methylation does not take place. The dependence of the occurrence of methyl transfer on the formation of a water channel in SET8 is in accord with our previous reports on product specificity by histone lysine monomethyltransferase SET7/9, large subunit lysine dimethyltransferase (LSMT), and viral histone lysine trimethyltransferase (vSET). The average value of the experimental DeltaG(E)() for the six lysine methyl transfer reactions catalyzed by vSET, LSMT, and SET7/9 with p53 as a substrate is 22.1 +/- 1.0 kcal/mol, and the computed average (DeltaG(C)()) is 22.2 +/- 0.8 kcal/mol. In this study, the computed free energy barrier of the methyl transfer reaction [Lys20-NH(2) + (+)AdoMet --> Lys20-N(Me)H(2)(+) + AdoHcy] catalyzed by SET8 is 20.8 kcal/mol. This is in agreement with the value of 20.6 kcal/mol calculated from the experimental rate constant (0.43 +/- 0.02 min(-1)). Our bond-order computations establish that the H4-K20 monomethylation in SET8 is a concerted linear S(N)2 displacement reaction.     

Reaction pathway and free energy profiles for butyrylcholinesterase-catalyzed hydrolysis of acetylthiocholine

The catalytic mechanism for butyrylcholineserase (BChE)-catalyzed hydrolysis of acetylthiocholine (ATCh) has been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical-free energy (QM/MM-FE) calculations on both acylation and deacylation of BChE. Additional quantum mechanical (QM) calculations have been carried out, along with the QM/MM-FE calculations, to understand the known substrate activation effect on the enzymatic hydrolysis of ATCh. It has been shown that the acylation of BChE with ATCh consists of two reaction steps including the nucleophilic attack on the carbonyl carbon of ATCh and the dissociation of thiocholine ester. The deacylation stage includes nucleophilic attack of a water molecule on the carboxyl carbon of substrate and dissociation between the carboxyl carbon of substrate and hydroxyl oxygen of Ser198 side chain. QM/MM-FE calculation results reveal that the acylation of BChE is rate-determining. It has also been demonstrated that an additional substrate molecule binding to the peripheral anionic site (PAS) of BChE is responsible for the substrate activation effect. In the presence of this additional substrate molecule at PAS, the calculated free energy barrier for the acylation stage (rate-determining step) is decreased by ~1.7 kcal/mol. All of our computational predictions are consistent with available experimental kinetic data. The overall free energy barriers calculated for BChE-catalyzed hydrolysis of ATCh at regular hydrolysis phase and substrate activation phase are ~13.6 and ~11.9 kcal/mol, respectively, which are in reasonable agreement with the corresponding experimentally derived activation free energies of 14.0 kcal/mol (for regular hydrolysis phase) and 13.5 kcal/mol (for substrate activation phase).     

Analysis of the decarboxylation step in mammalian histidine decarboxylase. A computational study

We report a hybrid quantum mechanics/molecular mechanics theoretical study on the reaction mechanism of mammalian histidine decarboxylase that allows us to obtain valuable insights on the structure of the cofactor-substrate adduct (external aldimine) in the active site of rat histidine decarboxylase. By means of molecular dynamics simulations, we traced the potential of mean force corresponding to the decarboxylation reaction of the adduct both in the active site of the enzyme and in aqueous solution. By comparing this process in both media, we have identified the key electrostatic interactions that explain the lowering of the free energy barrier in the enzyme. Our analysis also offers a validation of Dunathan's hypothesis (Dunathan, H. (1966) Proc. Natl. Acad. Sci. U. S. A. 55, 712-716) regarding the role of stereoelectronic effects in the enzymatic decarboxylation process.     

Computational characterization of the chemical step in the GTP hydrolysis by Ras‐GAP for the wild‐type and G13V mutated Ras

The free energy profiles for the chemical reaction of the guanosine triphosphate hydrolysis GTP + H₂O → GDP + Pi by Ras‐GAP for the wild‐type and G13V mutated Ras were computed by using molecular dynamics protocols with the QM(ab initio)/MM potentials. The results are consistent with the recent measurements of reaction kinetics in Ras‐GAP showing about two‐order reduction of the rate constant upon G13V mutation in Ras: the computed activation barrier on the free energy profile is increased by 3 kcal/mol upon the G13V replacement. The major reason for a higher energy barrier is a shift of the “arginine finger” (R789 from GAP) from the favorable position in the active site. The results of simulations provide support for the mechanism of the reference reaction according to which the Q61 side chain directly participates in chemical transformations at the proton transfer stage. Proteins 2015; 83:1046–1053. © 2015 Wiley Periodicals, Inc.     

Catalytic mechanism and product specificity of rubisco large subunit methyltransferase: QM/MM and MD investigations

Molecular dynamics (MD) simulations and hybrid quantum mechanics/molecular mechanics (QM/MM) calculations have been carried out in an investigation of Rubisco large subunit methyltransferase (LSMT). It was found that the appearance of a water channel is required for the stepwise methylation by S-adenosylmethionine (AdoMet). The water channel appears in the presence of AdoMet (LSMT.Lys-NH3+.AdoMet), but is not present immediately after methyl transfer (LSMT.Lys-N(Me)H2+.AdoHcy). The water channel allows proton dissociation from both LSMT.AdoMet.Lys-NH3+ and LSMT.AdoMet.Lys-N(Me)H2+. The water channel does not appear for proton dissociation from LSMT.AdoMet.Lys-N(Me)2H+, and a third methyl transfer does not occur. By QM/MM, the calculated free energy barrier of the first methyl transfer reaction catalyzed by LSMT (Lys-NH2 + AdoMet --> Lys-N(Me)H2+ + AdoHcy) is DeltaG++ = 22.8 +/- 3.3 kcal/mol. This DeltaG++ is in remarkable agreement with the value 23.0 kcal/mol calculated from the experimental rate constant (6.2 x 10-5 s-1). The calculated DeltaG++ of the second methyl transfer reaction (AdoMet + Lys-N(Me)H --> AdoHcy + Lys-N(Me)2H+) at the QM/MM level is 20.5 +/- 3.6 kcal/mol, which is in agreement with the value 22.0 kcal/mol calculated from the experimental rate constant (2.5 x 10-4 s-1). The third methyl transfer (Lys-N(Me)2 + AdoMet --> Lys-N(Me)3+ + AdoHcy) is associated with an allowed DeltaG++ of 25.9 +/- 3.2 kcal/mol. However, this reaction does not occur because a water channel does not form to allow the proton dissociation of Lys-N(Me)2H+. Future studies will determine whether the product specificity of lysine (mono, di, and tri) methyltransferases is determined by the formation of water channels.     

Role of conserved active site tryptophan-101 in functional activity and stability of phosphoserine aminotransferase from an enteric human parasite

Site-directed mutagenesis study was performed to elucidate the role of conserved tryptophan-101 present at the active site of phosphoserine aminotransferase from an enteric human parasite Entamoeba histolytica. Fluorescence resonance energy transfer and molecular dynamic simulation show that the indole ring of Trp101 stacks with the cofactor PLP. Loss of enzymatic activity and PLP polarization values suggest that Trp101 plays a major role in maintaining a defined PLP microenvironment essentially required for optimal enzymatic activity. Studies on W101F, W101H and W101A mutants show that only the indole ring of the conserved Trp101 forms most favorable stacking interaction with the pyridine ring of the cofactor PLP. Protein stability was compromised on substitution of Trp101 with Phe/His/Ala amino acids. A difference in conformational free energy of 1.65?kcal?mol(-1) was observed between WT-protein and W101A mutant.