Evidence from 13C NMR for polarization of the carbonyl of oxaloacetate in the active site of citrate synthase

The carbon-13 NMR spectrum of oxaloacetate bound in the active site of citrate synthase has been obtained at 90.56 MHz. In the binary complex with enzyme, the positions of the resonances of oxaloacetate are shifted relative to those of the free ligand as follows: C-1 (carboxylate), -2.5 ppm; C-2 (carbonyl), +4.3 ppm; C-3 (methylene), -0.6 ppm; C-4 (carboxylate), +1.3 ppm. The change observed in the carbonyl chemical shift is successively increased in ternary complexes with the product [coenzyme A (CoA)], a substrate analogue (S-acetonyl-CoA), and an acetyl-CoA enolate analogue (carboxymethyl-CoA), reaching a value of +6.8 ppm from the free carbonyl resonance. Binary complexes are in intermediate to fast exchange on the NMR time scale with free oxaloacetate; ternary complexes are in slow exchange. Line widths of the methylene resonance in the ternary complexes suggest complete immobilization of oxaloacetate in the active site. Analysis of line widths in the binary complex suggests the existence of a dynamic equilibrium between two or more forms of bound oxaloacetate, primarily involving C-4. The changes in chemical shifts of the carbonyl carbon indicate strong polarization of the carbonyl bond or protonation of the carbonyl oxygen. Some of this carbonyl polarization occurs even in the binary complex. Development of positive charge on the carbonyl carbon enhances reactivity toward condensation with the carbanion/enolate of acetyl-CoA in the mechanism which has been postulated for this enzyme. The very large change in the chemical shift of the reacting carbonyl in the presence of an analogue of the enolate of acetyl-CoA supports this interpretation.     

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.     

Acetyl-CoA enolization in citrate synthase: A quantum mechanical/molecular mechanical (QM/MM) study

Citrate synthase forms citrate by deprotonation of acetyl-CoA followed by nucleophilic attack of this substrate on oxaloacetate, and subsequent hydrolysis. The rapid reaction rate is puzzling because of the instability of the postulated nucleophilic intermediate, the enolate of acetyl-CoA. As alternatives, the enol of acetyl-CoA, or an enolic intermediate sharing a proton with His-274 in a “low-barrier” hydrogen bond have been suggested. Similar problems of intermediate instability have been noted in other enzymic carbon acid deprotonation reactions. Quantum mechanical/molecular mechanical calculations of the pathway of acetyl-CoA enolization within citrate synthase support the identification of Asp-375 as the catalytic base. His-274, the proposed general acid, is found to be neutral. The acetyl-CoA enolate is more stable at the active site than the enol, and is stabilized by hydrogen bonds from His-274 and a water molecule. The conditions for formation of a low-barrier hydrogen bond do not appear to be met, and the calculated hydrogen bond stabilization in the reaction is less than the gas-phase energy, due to interactions with Asp-375 at the active site. The enolate character of the intermediate is apparently necessary for the condensation reaction to proceed efficiently. Proteins 27:9–25 © 1997 Wiley-Liss, Inc.     

Catalytic strategy of citrate synthase: effects of amino acid changes in the acetyl-CoA binding site on transition-state analog inhibitor complexes.

Acetyl-CoA enol has been proposed as an intermediate in the citrate synthase (CS) reaction with Asp375 acting as a base, removing a proton from the methyl carbon of acetyl-CoA, and His274 acting as an acid, donating a proton to the carbonyl [Karpusas, M., Branchaud, B., & Remington, S.J. (1990) Biochemistry 29, 2213]. CS-oxaloacetate (OAA) complexes with the transition-state analog inhibitor, carboxymethyl-CoA (CMCoA), mimic those with acetyl-CoA enol. Asp375 and His274 interact intimately with the carboxyl of the bound inhibitor. While enzymes in which these residues have been changed to other amino acids have very low catalytic activity, we find that they retain their ability to form complexes with substrates and the transition-state analog inhibitor. In comparison with the value of the chemical shift of the protonated CMCoA carboxyl in acidic aqueous solutions or its value in the wild-type ternary complex, the values in the Asp375 mutants are unusually low. Model studies suggest that these low values result from complete absence of one hydrogen bond partner for the Gly mutant and distortions in the active site hydrogen bond systems for the Glu mutant. The high affinity of Asp375Gly-OAA for CMCoA suggests that the unfavorable proton uptake required to stabilize the CMCoA-OAA ternary complex of the wild-type enzyme [Kurz, L.C., Shah, S., Crane, B.R., Donald, L.J., Duckworth, H.W., & Drysdale, G.R. (1992) Biochemistry (preceding paper in this issue)] is not required by this mutant; the needed proton is most likely provided by His274. This supports the proposed role of His274 as a general acid.(ABSTRACT TRUNCATED AT 250 WORDS)     

Kinetics and mechanism of the citrate synthase from the thermophilic archaeon Thermoplasma acidophilum

The kinetics and mechanism of the citrate synthase from a moderate thermophile, Thermoplasma acidophilum (TpCS), are compared with those of the citrate synthase from a mesophile, pig heart (PCS). All discrete steps in the mechanistic sequence of PCS can be identified in TpCS. The catalytic strategies identified in PCS, destabilization of the oxaloacetate substrate carbonyl and stabilization of the reactive species, acetyl-CoA enolate, are present in TpCS. Conformational changes, which allow the enzyme to efficiently catalyze both condensation of acetyl-CoA thioester and subsequently hydrolysis of citryl-CoA thioester within the same active site, occur in both enzymes. However, significant differences exist between the two enzymes. PCS is a characteristically efficient enzyme: no internal step is clearly rate-limiting and the condensation step is readily reversible. TpCS is a less efficient catalyst. Over a broad temperature range, inadequate stabilization of the transition state for citryl-CoA hydrolysis renders this step nearly rate-limiting for the forward reaction of TpCS. Further, excessive stabilization of the citryl-CoA intermediate renders the condensation step nearly irreversible. Values of substrate and solvent deuterium isotope effects are consistent with the kinetic model. Near its temperature optimum (70 degrees C), there is a modest increase in the reversibility of the condensation step for TpCS, but reversibility still falls short of that shown by PCS at 37 degrees C. The root cause of the catalytic inefficiency of TpCS may lie in the lack of protein flexibility imposed by the requirement for thermal stability of the protein itself or its temperature-labile substrate, oxaloacetate.     

A quantum mechanics/molecular mechanics study of the reaction mechanism of the hepatitis C virus NS3 protease with the NS5A/5B substrate

Combined quantum mechanics and molecular mechanics (QM/MM) calculations were carried out to characterize the reaction mechanism of the NS3 protease with its preferred substrate (NS5A/5B). The main purpose of this study was to locate the barrier states and intermediates along the distinguished coordinate path (DCP) involved in this process. These structures, and in particular the one corresponding to the first barrier state and intermediate (B1 and I1), could be a starting point for the synthesis of inhibitors of this protease, which could be used to treat hepatitis C. The two first steps of the reaction mechanism were studied, i.e., the acylation step and the breaking of the peptide bond. The first step takes place through a tetracoordinated intermediate, as suggested from previous works on other Serine proteases. The importance of the different amino acid residues was also considered (perturbation study where the MM charges of each residue were set to zero independently). The residues of the oxyanion hole were confirmed as the most important for the electrostatic stabilization of the tetracoordinate intermediate. Moreover, the role of other residues, e.g., Arg-155 and Asp-79, was also explained.     

Energetics of the proposed rate-determining step of the glyoxalase I reaction.

The proposed rate-limiting step of the reaction catalyzed by glyoxalase I is the proton abstraction from the C1 carbon atom of the substrate by a glutamate residue, resulting in a high-energy enolate intermediate. This proton transfer reaction was modelled using molecular dynamics and free energy perturbation simulations, with the empirical valence bond method describing the potential energy surface of the system. The calculated rate constant for the reaction is approximately 300-1500 s(-1) with Zn2+, Mg2+ or Ca2+ bound to the active site, which agrees well with observed kinetics of the enzyme. Furthermore, the results imply that the origin of the catalytic rate enhancement is mainly associated with enolate stabilization by the metal ion.     

Reaction-path energetics and kinetics of the hydride transfer reaction catalyzed by dihydrofolate reductase

We have studied the hydride transfer reaction catalyzed by the enzyme dihydrofolate reductase (DHFR) and the coenzyme nicotinamide adenine dinucleotide phosphate (NADPH); the substrate is 5-protonated 7,8-dihydrofolate, and the product is tetrahydrofolate. The potential energy surface is modeled by a combined quantum mechanical-molecular mechanical (QM/MM) method employing Austin model 1 (AM1) and a simple valence bond potential for 69 QM atoms and employing the CHARMM22 and TIP3P molecular mechanics force fields for the other 21 399 atoms; the QM and MM regions are joined by two boundary atoms treated by the generalized hybrid orbital (GHO) method. All simulations are carried out using periodic boundary conditions at neutral pH and 298 K. In stage 1, a reaction coordinate is defined as the difference between the breaking and forming bond distances to the hydride ion, and a quasithermodynamic free energy profile is calculated along this reaction coordinate. This calculation includes quantization effects on bound vibrations but not on the reaction coordinate, and it is used to locate the variational transition state that defines a transition state ensemble. Then, the key interactions at the reactant, variational transition state, and product are analyzed in terms of both bond distances and electrostatic energies. The results of both analyses support the conclusion derived from previous mutational studies that the M20 loop of DHFR makes an important contribution to the electrostatic stabilization of the hydride transfer transition state. Third, transmission coefficients (including recrossing factors and multidimensional tunneling) are calculated and averaged over the transition state ensemble. These averaged transmission coefficients, combined with the quasithermodynamic free energy profile determined in stage 1, allow us to calculate rate constants, phenomenological free energies of activation, and primary and secondary kinetic isotope effects. A primary kinetic isotope effect (KIE) of 2.8 has been obtained, in good agreement with the experimentally determined value of 3.0 and with the value 3.2 calculated previously. The primary KIE is mainly a consequence of the quantization of bound vibrations. In contrast, the secondary KIE, with a value of 1.13, is almost entirely due to dynamical effects on the reaction coordinate, especially tunneling.     

The reaction mechanism of paraoxon hydrolysis by phosphotriesterase from combined QM/MM simulations

Molecular dynamics simulations employing combined quantum mechanical and molecular mechanical (QM/MM) potentials have been carried out to investigate the reaction mechanism of the hydrolysis of paraoxon by phosphotriesterase (PTE). We used a dual-level QM/MM approach that synthesizes accurate results from high-level electronic structure calculations with computational efficiency of semiempirical QM/MM potentials for free energy simulations. In particular, the intrinsic (gas-phase) energies of the active site in the QM region are determined by using density functional theory (B3LYP) and second-order M?ller-Plesset perturbation theory (MP2) and the molecular dynamics free energy simulations are performed by using the mixed AM1:CHARMM potential. The simulation results suggest a revised mechanism for the phosphotriester hydrolysis mechanism by PTE. The reaction free energy profile is mirrored by structural motions of the binuclear metal center in the active site. The two zinc ions occupy a compact conformation with an average zinc-zinc distance of 3.5 +/- 0.1 A in the Michaelis complex, whereas it is elongated to 5.3 +/- 0.3 A at the transition state and product state. The substrate is loosely bound to the more exposed zinc ion (Znbeta2+) at an average distance of 3.8 A +/- 0.3 A. The P=O bond of the substrate paraoxon is activated by adopting a tight coordination to the Znbeta2+, releasing the coordinate to the bridging hydroxide ion and increasing its nucleophilicity. It was also found that a water molecule enters into the binding pocket of the loosely bound binuclear center, originally occupied by the nucleophilic hydroxide ion. We suggest that the proton of this water molecule is taken up by His254 at low pH or released to the solvent at high pH, resulting in a hydroxide ion that pulls the Znbeta2+ ion closer to form the compact configuration and restores the resting state of the enzyme.     

Modelling flavin and substrate substituent effects on the activation barrier and rate of oxygen transfer by p-hydroxybenzoate hydroxylase

The simulation of enzymatic reactions, using computer models, is becoming a powerful tool in the most fundamental challenge in biochemistry: to relate the catalytic activity of enzymes to their structure. In the present study, various computed parameters were correlated with the natural logarithm of experimental rate constants for the hydroxylation of various substrate derivatives catalysed by wild-type para-hydroxybenzoate hydroxylase (PHBH) as well as for the hydroxylation of the native substrate (p-hydroxybenzoate) by PHBH reconstituted with a series of 8-substituted flavins. The following relative parameters have been calculated and tested: (a) energy barriers from combined quantum mechanical/molecular mechanical (QM/MM) (AM1/CHARMM) reaction pathway calculations, (b) gas-phase reaction enthalpies (AM1) and (c) differences between the HOMO and LUMO energies of the isolated substrate and cofactor molecules (AM1 and B3LYP/6-31+G(d)). The gas-phase approaches yielded good correlations, as long as similarly charged species are involved. The QM/MM approach resulted in a good correlation, even including differently charged species. This indicates that the QM/MM model accounts quite well for the solvation effects of the active site surroundings, which vary for differently charged species. The correlations obtained demonstrate quantitative structure activity relationships for an enzyme-catalysed reaction including, for the first time, substitutions on both substrate and cofactor.     

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.     

Substrate distortion in the Michaelis complex of Bacillus 1,3-1,4-beta-glucanase. Insight from first principles molecular dynamics simulations

The structure and dynamics of the enzyme-substrate complex of Bacillus 1,3-1,4-beta-glucanase, one of the most active glycoside hydrolases, is investigated by means of Car-Parrinello molecular dynamics simulations (CPMD) combined with force field molecular dynamics (QM/MM CPMD). It is found that the substrate sugar ring located at the -1 subsite adopts a distorted 1S3 skew-boat conformation upon binding to the enzyme. With respect to the undistorted 4C1 chair conformation, the 1S3 skew-boat conformation is characterized by: (a) an increase of charge at the anomeric carbon (C1), (b) an increase of the distance between C1 and the leaving group, and (c) a decrease of the intraring O5-C1 distance. Therefore, our results clearly show that the distorted conformation resembles both structurally and electronically the transition state of the reaction in which the substrate acquires oxocarbenium ion character, and the glycosidic bond is partially broken. Together with analysis of the substrate conformational dynamics, it is concluded that the main determinants of substrate distortion have a structural origin. To fit into the binding pocket, it is necessary that the aglycon leaving group is oriented toward the beta region, and the skew-boat conformation naturally fulfills this premise. Only when the aglycon is removed from the calculation the substrate recovers the all-chair conformation, in agreement with the recent determination of the enzyme product structure. The QM/MM protocol developed here is able to predict the conformational distortion of substrate binding in glycoside hydrolases because it accounts for polarization and charge reorganization at the -1 sugar ring. It thus provides a powerful tool to model E.S complexes for which experimental information is not yet available.     

Enzyme polarization of substrates of dihydrofolate reductase by different theoretical methods

We have investigated the importance of polarization by the enzyme dihydrofolate reductase (DHFR) on its substrates, folate and dihydrofolate, using a series of quantum mechanical (QM) techniques (Hartree-Fock (HF), M?ller-Plesset second-order perturbation theory (MP2), local density approximation (LDA) and generalized gradient approximation (GGA) density functional theory (DFT) calculations) in which the bulk enzyme is included in the calculations as point charges. Polarization, in terms of both charges on components (residues) of the folate and dihydrofolate molecules and changes in the electron density, particularly of the pterin ring of the substrates, and the implications for the catalytic reduction are discussed. Significant differences in polarization behavior are observed for the different theoretical methods employed. The consequences of this, particularly for choosing an appropriate model for quantum mechanical/molecular mechanical (QM/MM) calculations, are pointed out. The HF and MP2 QM methods show small polarizations (approximately 0.04 electrons) of the pterin ring but quite large polarizations with both LDA and GGA DFT methods (0.3-0.5 electrons). This large difference in polarization for both folate and dihydrofolate arises as a result of substantial differences between the charge distributions for the gasphase DFT and HF calculations, specifically the charges on the dianionic glutamate side chain. Some recent literature reports of incorrect representation of anionic systems by DFT methods are noted. The DFT results are similar to the previously reported LDA DFT results of Bajorath et al. predicting a large polarization of the pterin ring of folate (Proteins 9:217-224, 1991) and dihydrofolate (PNAS 88:6423-6426, 1991) of approximately 0.5-0.6 electrons.     

Structural and mechanistic studies of enolase

The high-resolutionstructure of yeast enolase cocrystallized with its equilibrium mixture of substrate and product reveals the stereochemistry of substrate/product binding and therefore the groups responsible for acid/base catalysis and stabilization of the enolate intermediate. Expression and characterization of site-specific mutant forms of the enzyme have confirmed the roles of amino acid side chains in the catalysis of the first and second steps of the reaction. Coordination of both required magnesium ions to the carboxylate of the substrate/product indicates a role for these cations in stabilization of the intermediate.     

Exploring the RING-Catalyzed Ubiquitin Transfer Mechanism by MD and QM/MM Calculations

Ubiquitylation is a universal mechanism for controlling cellular functions. A large family of ubiquitin E3 ligases (E3) mediates Ubiquitin (Ub) modification. To facilitate Ub transfer, RING E3 ligases bind both the substrate and ubiquitin E2 conjugating enzyme (E2) linked to Ub via a thioester bond to form a catalytic complex. The mechanism of Ub transfer catalyzed by RING E3 remains elusive. By employing a combined computational approach including molecular modeling, molecular dynamics (MD) simulations, and quantum mechanics/molecular mechanics (QM/MM) calculations, we characterized this catalytic mechanism in detail. The three-dimensional model of dimeric RING E3 ligase RNF4 RING, E2 ligase UbcH5A, Ub and the substrate SUMO2 shows close contact between the substrate and Ub transfer catalytic center. Deprotonation of the substrate lysine by D117 on UbcH5A occurs with almost no energy barrier as calculated by MD and QM/MM calculations. Then, the side chain of the activated lysine gets close to the thioester bond via a conformation change. The Ub transfer pathway begins with a nucleophilic addition that forms an oxyanion intermediate of a 4.23 kcal/mol energy barrier followed by nucleophilic elimination, resulting in a Ub modified substrate by a 5.65 kcal/mol energy barrier. These results provide insight into the mechanism of RING-catalyzed Ub transfer guiding the discovery of Ub system inhibitors.     

QM/MM study of the mechanism of enzymatic limonene 1,2-epoxide hydrolysis

Limonene 1,2-epoxide hydrolase (LEH) is completely different from those of classic epoxide hydrolases (EHs) which catalyze the hydrolysis of epoxides to vicinal diols. A novel concerted general acid catalysis step involving the Asp101-Arg99-Asp132 triad is proposed to play an important role in the mechanism. Combined quantum-mechanical/molecular-mechanical (QM/MM) calculations gave activation barriers of 16.9 and 25.1 kcal/mol at the B3LYP/6-31G(d,p)//CHARMM level for nucleophilic attack on the more and less substituted epoxide carbons, respectively. Furthermore, the important roles of residues Arg99, Tyr53 and Asn55 on mutated LEH were evaluated by QM/MM-scanned energy mapping. These results may provide an explanation for site-directed mutagenesis.     

Substrate diffusion and oxidation in GMC oxidoreductases: an experimental and computational study on fungal aryl-alcohol oxidase

AAO (aryl-alcohol oxidase) provides H?O? in fungal degradation of lignin, a process of high biotechnological interest. The crystal structure of AAO does not show open access to the active site, where different aromatic alcohols are oxidized. In the present study we investigated substrate diffusion and oxidation in AAO compared with the structurally related CHO (choline oxidase). Cavity finder and ligand diffusion simulations indicate the substrate-entrance channel, requiring side-chain displacements and involving a stacking interaction with Tyr?2. Mixed QM (quantum mechanics)/MM (molecular mechanics) studies combined with site-directed mutagenesis showed two active-site catalytic histidine residues, whose substitution strongly decreased both catalytic and transient-state reduction constants for p-anisyl alcohol in the H502A (over 1800-fold) and H546A (over 35-fold) variants. Combination of QM/MM energy profiles, protonation predictors, molecular dynamics, mutagenesis and pH profiles provide a robust answer regarding the nature of the catalytic base. The histidine residue in front of the FAD ring, AAO His??2 (and CHO His???), acts as a base. For the two substrates assayed, it was shown that proton transfer preceded hydride transfer, although both processes are highly coupled. No stable intermediate was observed in the energy profiles, in contrast with that observed for CHO. QM/MM, together with solvent KIE (kinetic isotope effect) results, suggest a non-synchronous concerted mechanism for alcohol oxidation by AAO.     

Detailed reaction mechanism of macrophomate synthase. Extraordinary enzyme catalyzing five-step transformation from 2-pyrones to benzoates

Macrophomate synthase from the fungus Macrophoma commelinae IFO 9570 is a Mg(II)-dependent dimeric enzyme that catalyzes an extraordinary, complex five-step chemical transformation from 2-pyrone and oxalacetate to benzoate involving decarboxylation, C-C bond formation, and dehydration. The catalytic mechanism of the whole pathway was investigated in three separate chemical steps. In the first decarboxylation step, the enzyme loses oxalacetate decarboxylation activity upon incubation with EDTA. Activity is fully restored by addition of Mg(II) and is not restored with other divalent metal cations. The dissociation constant of 0.93 x 10(-)(7) for Mg(II) and atomic absorption analysis established a 1:1 stoichiometric complex. Inhibition of pyruvate formation with 2-pyrone revealed that the actual product in the first step is a pyruvate enolate, which undergoes C-C bond formation in the presence of 2-pyrone. Incubation of substrate analogs provided aberrant adducts that were produced via C-C bond formation and rearrangement. This strongly indicates that the second step is two C-C bond formations, affording a bicyclic intermediate. Based on the stereospecificity, involvement of a Diels-Alder reaction at the second step is proposed. Incubation of the stereospecifically deuterium-labeled malate with 2-pyrones in the presence of malate dehydrogenase provided information for the stereochemical course of the reaction catalyzed by macrophomate synthase, indicating that the first decarboxylation provides pyruvate (Z)-[3-(2)H]enolate and that dehydration at the final step occurs with anti-elimination accompanied by concomitant decarboxylation. Examination of kinetic parameters in the individual steps suggests that the third step is the rate-determining step of the overall transformation.     

Molecular determinants of xenobiotic metabolism: QM/MM simulation of the conversion of 1-chloro-2,4-dinitrobenzene catalyzed by M1-1 glutathione S-transferase

Modeling methods allow the identification and analysis of determinants of reactivity and specificity in enzymes. The reaction between glutathione and 1-chloro-2,4-dinitrobenzene (CDNB) is widely used as a standard activity assay for glutathione S-transferases (GSTs). It is important to understand the causes of differences between catalytic GST isoenzymes and the effects of mutations and genetic polymorphisms. Quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulations have been performed here to investigate the addition of the glutathione anion to CDNB in the wild-type M1-1 GST isoenzyme from rat and in three single point mutant (Tyr6Phe, Tyr115Phe, and Met108Ala) M1-1 GST enzymes. We have developed a specifically parameterized QM/MM method (AM1-SRP/CHARMM22) to model this reaction by fitting to experimental heats of formation and ionization potentials. Free energy profiles were obtained from molecular dynamics simulations of the reaction using umbrella sampling and weighted histogram analysis techniques. The reaction in solution has also been simulated and is compared to the enzymatic reaction. The free energies are in excellent agreement with experimental results. Overall the results of the present study show that QM/MM reaction pathway analysis provides detailed insight into the chemistry of GST and can be used to obtain mechanistic insight into the effects of specific mutations on this catalytic process.     

Application of a SCC-DFTB QM/MM approach to the investigation of the catalytic mechanism of fatty acid amide hydrolase

Self-consistent charge density functional tight binding (SCC-DFTB) is a promising method for hybrid quantum mechanics/molecular mechanics (QM/MM) simulations of enzyme-catalyzed reactions. The acylation reaction of fatty acid amide hydrolase (FAAH), a promising drug target, was investigated by applying a SCC-DFTB/CHARMM27 scheme. Calculated potential energy barriers resulted in reasonable agreement with experiments for oleamide (OA) and oleoylmethyl ester (OME) substrates, outperforming previous calculations performed at the PM3/CHARMM22 level. Furthermore, the experimental preference of FAAH in hydrolyzing OA faster than OME was adequately reproduced by calculations. All these findings indicate that the SCC-DFTB/CHARMM27 approach can be successfully applied to mechanistic investigations of FAAH-catalyzed reactions.