Catalytic reaction mechanism of L-lactate dehydrogenase: an ab initio study

Studies on the catalytic reaction mechanism of L-lactate dehydrogenase have been carried out by using quantum chemical ab initio calculation at HF/6-31G* level. It is found that the interconversion reaction of pyruvate to L-lactate is dominated by the hydride ion HR- transfer, and the transfers of the hydride ion HR and proton HR are a quasi-coupled process, in which the energy barrier of the transition state is about 168.37 kJ/mol. It is shown that the reactant complex is 87.61 kJ/mol lower, in energy, than the product complex. The most striking features in our calculated results are that pyridine ring of the model cofactor is a quasi-boat-like configuration in the transited state, which differs from a planar conformation in some previous semiempirical quantum chemical studies. On the other hand, the similarity in the structure and charge between the HR transfer process and the hydrogen bonding with lower barrier indicates that the HR transfer process occurs by means of an unusual manner. In addition,     

On Transition Structures for Hydride Transfer Step in Enzyme Catalysis. A Comparative Study on Models of Glutathione Reductase Derived from Semiempirical, HF, and DFT Methods

As a model of the chemical reactions that take place in the active site of gluthatione reductase, the nature of the molecular mechanism for the hydride transfer step has been characterized by means of accurate quantum chemical characterizations of transition structures. The calculations have been carried out with analytical gradients at AM1 and PM3 semiempirical procedures, ab initio at HF level with 3-21G, 4-31G, 6-31G, and 6-31G basis sets and BP86 and BLYP as density functional methods. The results of this study suggest that the endo relative orientation on the substrate imposed by the active site is optimal in polarizing the C4-Ht bond and situating the system in the neighborhood of the quadratic region of the transition structure associated to the hydride transfer step on potential energy surface. The endo arrangement of the transition structure results in optimal frontier HOMO orbital interaction between NADH and FAD partners. The geometries of the transition structures and the corresponding transition vectors, that contain the fundamental information relating reactive fluctuation patterns, are model independent and weakly dependent on the level of theory used to determine them. A comparison between simple and complex molecular models shows that there is a minimal set of coordinates describing the essentials of hydride transfer step. The analysis of transition vector components suggests that the primary and secondary kinetic isotope effects can be strongly coupled, and this prompted the calculation of deuterium and tritium primary, secondary, and primary and secondary kinetic isotope effects. The results obtained agree well with experimental data and demonstrate this coupling.     

Catalytic reaction mechanism of L-lactate dehydrogenase: an ab initio study

Studies on the catalytic reaction mechanism of L-lactate dehydrogenase have been carried out by using quantum chemical ab initio calculation at HF/6-31G* level. It is found that the interconversion reaction of pyruvate to L-lactate is dominated by the hydride ion H_r transfer, and the transfers of the hydride ionH _r and protonH _r are a quasi-coupled process, in which the energy barrier of the transition state is about 168.37 kJ/mol. It is shown that the reactant complex is 87.61 kJ/mol lower, in energy, than the product complex. The most striking features in our calculated results are that pyridine ring of the model cofactor is a quasi-boat-like configuration in the transited state, which differs from a planar conformation in some previous semiempirical quantum chemical studies. On the other hand, the similarity in the structure and charge between theH _r transfer process and the hydrogen bonding with lower barrier indicates that the H_r transfer process occurs by means of an unusual manner. In addition, in the transition state the electrostatic interaction between the substrate and the active site of LDH is quite strong and the polarization of the carbonyl in the substrate is gradually enhanced accompanying the formation of the transition state. These calculated results are well in accord with the previous experimental studies, and indicate that the charge on the hydride ion H_r is only +0.13e in the transition state, which is in agreement with the reported semiempirical quantum chemical calculations.     

Catalytic reaction mechanism of L-lactate dehydrogenase: anab initio study

Studies on the catalytic reaction mechanism of L-lactate dehydrogenase have been carried out by using quantum chemical ab initio calculation at HF/6-31G* level. It is found that the interconversion reaction of pyruvate to L-lactate is dominated by the hydride ion H_r transfer, and the transfers of the hydride ionH _r and protonH _r are a quasi-coupled process, in which the energy barrier of the transition state is about 168.37 kJ/mol. It is shown that the reactant complex is 87.61 kJ/mol lower, in energy, than the product complex. The most striking features in our calculated results are that pyridine ring of the model cofactor is a quasi-boat-like configuration in the transited state, which differs from a planar conformation in some previous semiempirical quantum chemical studies. On the other hand, the similarity in the structure and charge between theH _r transfer process and the hydrogen bonding with lower barrier indicates that the H_r transfer process occurs by means of an unusual manner. In addition, in the transition state the electrostatic interaction between the substrate and the active site of LDH is quite strong and the polarization of the carbonyl in the substrate is gradually enhanced accompanying the formation of the transition state. These calculated results are well in accord with the previous experimental studies, and indicate that the charge on the hydride ion H_r is only +0.13e in the transition state, which is in agreement with the reported semiempirical quantum chemical calculations.     

Identity hydride-ion transfer from C-H donors to C acceptor sites. Enthalpies of hydride addition and enthalpies of activation. Comparison with C...H...C proton transfer. An ab initio study

Enthalpies of addition of hydride ion to eleven carbonyl acceptors (X-CHO), two conjugate addition sites (X-CH=CH2; X = CHO, NO2), eight carbenium ion acceptors, fulvene, borane, and SiH3(+) were calculated at the MP2/6-311+G level. Correlation between calculated and experimental enthalpies of addition of hydride ion is excellent. Transition states (ts) for the identity hydride transfers between the acceptors and their corresponding hydride adducts (hydride donors) were also calculated. The carbonyl and fulvene reactions have transition states with one imaginary frequency: the hydrogen transfer coordinate. The carbenium ions, borane, and SiH3(+) gave not transition states but stable compounds upon addition of the hydride donor. Computational differences between these hydride transfers and previously reported proton transfers include shorter partial C...H bonds and a tendency toward bent C...H...C angles for the hydride transfer ts and addition compound structures, particularly when a bent geometry improves interactions elsewhere in the structure. These and other differences are explained by modeling the hydride transfer ts and addition compounds as two-electron, three-center systems involving the transfer termini and the shared hydrogen but the proton transfer ts structures as four-electron, three-center systems. Charge and geometry measures suggest transition states in which these features change synchronously, again in contrast to many proton transfer reactions. For the X-CHO set, polar effects dominate enthalpies of hydride addition, with resonance effects also important for resonance donors; these preferentially stabilize the acceptor, reducing its hydride ion affinity. Activation enthalpies are dominated by resonance stabilization of the acceptors, greatly attenuated in the transition states.     

Comparison of coupled motions in Escherichia coli and Bacillus subtilis dihydrofolate reductase

Hybrid quantum/classical molecular dynamics simulations are used to compare the role of protein motion in the hydride transfer reaction catalyzed by Escherichia coli and Bacillus subtilis dihydrofolate reductase (DHFR). These two enzymes have 44% sequence identity, and the experimentally determined structures and hydride transfer rates are similar. The simulations indicate that the tertiary structures of both enzymes evolve in a similar manner during the hydride transfer reaction. In both enzymes, the donor-acceptor distance decreases to approximately 2.7 Angstroms at the transition state configurations to enable hydride transfer. Zero point energy and hydrogen tunneling effects are found to be significant for both enzymes. Covariance and rank correlation analyses of motions throughout the protein and ligands illustrate that E. coli and B. subtilis DHFR exhibit both similarities and differences in the equilibrium fluctuations and the conformational changes along the collective reaction coordinate for hydride transfer. A common set of residues that play a significant role in the network of coupled motions leading to configurations conducive to hydride transfer for both E. coli and B. subtilis DHFR was identified. These results suggest a balance between conservation and flexibility in the thermal motions and conformational changes during hydride transfer. Homologous protein structures, in conjunction with conformational sampling, enable enzymes with different sequences to catalyze the same hydride transfer reaction with similar efficiency.     

Hydride transfer reaction catalyzed by hyperthermophilic dihydrofolate reductase is dominated by quantum mechanical tunneling and is promoted by both inter- and intramonomeric correlated motions

Simulations of hydride and deuteride transfer catalyzed by dihydrofolate reductase from the hyperthermophile Thermotoga maritima (TmDHFR) are presented. TmDHFR was modeled with its active homodimeric quaternary structure, where each monomer has three subdomains. The potential energy function was a combined quantum mechanical and molecular mechanical potential (69 atoms were treated quantum mechanically, and 35 287, by molecular mechanics). The calculations of the rate constants by ensemble-averaged variational transition state theory with multidimensional tunneling predicted that hydride and deuteride transfer at 278 K proceeded with 81 and 80% by tunneling. These percentages decreased to 50 and 49% at 338 K. The kinetic isotope effect was dominated by contributions of bound vibrations and decreased from 3.0 to 2.2 over the temperature range. The calculated rates for hydride and deuteride transfer catalyzed by the hypothetical monomer were smaller by approximately 2 orders of magnitude. At 298 K tunneling contributed 73 and 66% to hydride and deuteride transfer in the monomer. The decreased catalytic efficiency of the monomer was therefore not the result of a decrease of the tunneling contribution but an increase in the quasi-classical activation free energy. The catalytic effect was associated in the dimer with correlated motions between domains as well as within and between subunits. The intrasubunit correlated motions were decreased in the monomer when compared to both native dimeric TmDHFR and monomeric E. coli enzyme. TmDHFR and its E. coli homologue involve similar patterns of correlated interactions that affect the free energy barrier of hydride transfer despite only 27% sequence identity and different quaternary structures.     

Development of predictive models of π-facial selectivity; a critical study of nucleophilic addition to sterically unbiased ketones

Quantum chemical calculations at B3LYP/6-31G* and semiempirical levels have been performed on a series of sterically unbiased ketones, where facial differentiation during nucleophilic additions is electronically induced through distal functional groups. The face selectivity data for fifty-four substrates representing nine different skeleta were computed and compared with the available experimental data on thirty-eight of them. The predictive abilities of various computational methods such as, charge model, hydride model, LiH transition state model, Cieplak hyperconjugation effect estimated by NBO analysis and the cation complexation model have been evaluated. A comparison of the computed and experimental face-selectivity data indicates that the hydride model and the LiH transition state model at the semiempirical levels are the best choices to predict diastereoselectivity. Unexpectedly, the performance of charge, hydride and LiH transition state models are inferior at the B3LYP level compared to the semiempirical methods in predicting the facial selectivities. On the other hand, the Cieplak type hyperconjugation evaluated using the NBO analysis, and cation complexation model are less reliable despite the fact that these two involve higher (B3LYP/6-31G*) level calculations. The inadequate performance of the charge model, NBO and the cation complexation models were traced to their emphasis on only one or two factors which are responsible for stereodifferentiation and undermining of the other subtle aspects involving a combination of orbital and electrostatic effects. On the other hand, the hydride and LiH transition state models, at semiempirical levels, provide reliable results to model the face-selectivities.     

Freezing a single distal motion in dihydrofolate reductase

Constraining a single motion between distal residues separated by approximately 28 A in hybrid quantum/classical molecular dynamics simulations is found to increase the free energy barrier for hydride transfer in dihydrofolate reductase by approximately 3 kcal/mol. Our analysis indicates that a single distal constraint alters equilibrium motions throughout the enzyme on a wide range of time scales. This alteration of the conformational sampling of the entire system is sufficient to significantly increase the free energy barrier and decrease the rate of hydride transfer. Despite the changes in conformational sampling introduced by the constraint, the system assumes a similar transition state conformation with a donor-acceptor distance of approximately 2.72 A to enable the hydride transfer reaction. The modified thermal sampling leads to a substantial increase in the average donor-acceptor distance for the reactant state, however, thereby decreasing the probability of sampling the transition state conformations with the shorter distances required for hydride transfer. These simulations indicate that fast thermal fluctuations of the enzyme, substrate, and cofactor lead to conformational sampling of configurations that facilitate hydride transfer. The fast thermal motions are in equilibrium as the reaction progresses along the collective reaction coordinate, and the overall average equilibrium conformational changes occur on the slower time scale measured experimentally. Recent single molecule experiments suggest that at least some of these thermally averaged equilibrium conformational changes occur on the millisecond time scale of the hydride transfer reaction. Thus, introducing a constraint that modifies the conformational sampling of an enzyme could significantly impact its catalytic activity.     

Mechanism of gas-phase aldose-ketose isomerization: A study using tandem mass spectrometry and theoretical calculations

Fast-atom bombardment tandem mass spectrometry and the semiempirical molecular orbital method were used to investigate the mechanism of gas-phase aldose-ketose isomerization process in lithiated 1,3 linked disaccharide isomers. Both the 1,3 hydrogen shift and hydride transfer mechanisms were investigated. Our experimental and theoretical calculations support the latter. The hydride transfer mechanism in these lithium-coordinated systems is similar to the xylose isomerase catalyzed aldose-ketose isomerization.     

Quantum chemical study on enantioselective reduction of aromatic ketones catalyzed by chiral cyclic sulfur‐containing oxazaborolidines. Part 2. Structures of catalyst–borane–ketone adducts

In the present paper, the ab initio molecular orbital method is employed to study the structures of the adducts of borane and aromatic ketone to chiral cyclic sulfur‐containing oxazaborolidine used as a catalyst in the enantioselective reduction of aromatic ketone. The catalyst–borane–ketone adducts have four different structures. All the structures are optimized completely by means of the Hartree–Fock method at 6‐31g* basis sets. The structure which is of the greatest advantage to a hydride transfer from the borane moiety to the carbonyl carbon of aromatic ketone is the one with the next lowest formation energy, and the plausible transition state for the hydride transfer is predicted to be of a twisted boat structure. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 78: 252–260, 2000     

Impact of protein flexibility on hydride-transfer parameters in thermophilic and psychrophilic alcohol dehydrogenases

The hydride transfer catalyzed by thermophilic alcohol dehydrogenase (htADH) exhibits sharply different kinetic and activation parameters from that catalyzed by the more flexible psychrophilic alcohol dehydrogenase (psADH). In addition, the hydride transfer in htADH is affected by mutating two distal residues that are suggested to be responsible for the decreased local protein flexibility in htADH. These observations provide support for the view that protein dynamics is tightly coupled to the hydrogen-transfer processes in these enzymes.     

A comparative computational investigation on the proton and hydride transfer mechanisms of monoamine oxidase using model molecules

Monoamine oxidase (MAO) enzymes regulate the level of neurotransmitters by catalyzing the oxidation of various amine neurotransmitters, such as serotonin, dopamine and norepinephrine. Therefore, they are the important targets for drugs used in the treatment of depression, Parkinson, Alzeimer and other neurodegenerative disorders. Elucidation of MAO-catalyzed amine oxidation will provide new insights into the design of more effective drugs. Various amine oxidation mechanisms have been proposed for MAO so far, such as single electron transfer mechanism, polar nucleophilic mechanism and hydride mechanism. Since amine oxidation reaction of MAO takes place between cofactor flavin and the amine substrate, we focus on the small model structures mimicking flavin and amine substrates so that three model structures were employed. Reactants, transition states and products of the polar nucleophilic (proton transfer), the water-assisted proton transfer and the hydride transfer mechanisms were fully optimized employing various semi-empirical, ab initio and new generation density functional theory (DFT) methods. Activation energy barriers related to these mechanisms revealed that hydride transfer mechanism is more feasible.     

Effect of mutation on enzyme motion in dihydrofolate reductase

Hybrid quantum-classical molecular dynamics simulations of a mutant Escherichia coli dihydrofolate reductase enzyme are presented. Although residue 121 is on the exterior of the enzyme, experimental studies have shown that the mutation of Gly-121 to valine reduces the rate of hydride transfer by a factor of 163. The simulations indicate that the decrease in the hydride transfer rate for the G121V mutant is due to an increase in the free energy barrier. The calculated free energy barrier is higher for the mutant than for the wild-type enzyme by an amount that is consistent with the experimentally observed rate reduction. The calculated transmission coefficients are comparable for the wild-type and mutant enzymes. The simulations suggest that this mutation may interrupt a network of coupled promoting motions proposed to play an important role in DHFR catalysis. This phenomenon has broad implications for protein engineering and drug design.     

Nonperfect synchronization of reaction center rehybridization in the transition state of the hydride transfer catalyzed by dihydrofolate reductase

It has been suggested that the magnitudes of secondary kinetic isotope effects (2 degrees KIEs) of enzyme-catalyzed reactions are an indicator of the extent of reaction-center rehybridization at the transition state. A 2 degrees KIE value close to the corresponding secondary equilibrium isotope effects (2 degrees EIE) is conventionally interpreted as indicating a late transition state that resembles the final product. The reliability of using this criterion to infer the structure of the transition state is examined by carrying out a theoretical investigation of the hybridization states of the hydride donor and acceptor in the Escherichia coli dihydrofolate reductase (ecDHFR)-catalyzed reaction for which a 2 degrees KIE close to the 2 degrees EIE was reported. Our results show that the donor carbon at the hydride transfer transition state resembles the reactant state more than the product state, whereas the acceptor carbon is more productlike, which is a symptom of transition state imbalance. The conclusion that the isotopically substituted carbon is reactant-like disagrees with the conclusion that would have been derived from the criterion of 2 degrees KIEs and 2 degrees EIEs, but the breakdown of the correlation with the equilibrium isotope effect can be explained by considering the effect of tunneling.     

DFT-based prediction of reactivity of short-chain alcohol dehydrogenase

The reaction mechanism of ketone reduction by short chain dehydrogenase/reductase, (S)-1-phenylethanol dehydrogenase from Aromatoleum aromaticum, was studied with DFT methods using cluster model approach. The characteristics of the hydride transfer process were investigated based on reaction of acetophenone and its eight structural analogues. The results confirmed previously suggested concomitant transfer of hydride from NADH to carbonyl C atom of the substrate with proton transfer from Tyr to carbonyl O atom. However, additional coupled motion of the next proton in the proton-relay system, between O2′ ribose hydroxyl and Tyr154 was observed. The protonation of Lys158 seems not to affect the pKa of Tyr154, as the stable tyrosyl anion was observed only for a neutral Lys158 in the high pH model. The calculated reaction energies and reaction barriers were calibrated by calorimetric and kinetic methods. This allowed an excellent prediction of the reaction enthalpies (R²?=?0.93) and a good prediction of the reaction kinetics (R²?=?0.89). The observed relations were validated in prediction of log K _eq obtained for real whole-cell reactor systems that modelled industrial synthesis of S-alcohols.     

Comparison of AM1 and PM3 semiempirical to ab initio methods in the study of Diels-Alder reactions of butadiene and cyclopentadiene with cyanoethylenes

Assuming a concerted synchronous mechanism with one transition state of the Diels-Alder reactions, the structures of the transition states and the activation energies for the reactions of butadiene and cyclopentadiene with cyanoethylenes were calculated by AM1 and PM3 semiempirical methods. The structural parameters were compared with those obtained by high level Gaussian calculations, whereas the activation energies were compared both with the ab initio calculations and those obtained experimentally. The structural properties calculated with PM3 methods are in general in better agreement with the ab initio calculations. The low level ab initio calculations are in many cases worse than the semiempirical methods. All predicted activation energies with both semiempirical methods are up to 300% higher than the experimental values. The predicted reactivity is also opposite to the experimental data. Only the very high level Gaussian calculations are in good correlation with experimental results. The predicted selectivity of the reaction is also opposite to the experimental facts. Two explanations are offered for this discrepancy: AM1 and PM3 methods cannot handle the calculation of the concerted Diels-Alder transition states and are not recommended to be used for that purpose, or this Diels-Alder reaction is not concerted but is stepwise.     

Quantum chemical study of the intramolecular transfer of hydrogen in o- methylacetophenone and 1- alkylanthraquinones

The geometries and formation enthalpies of reagents, products, and triplet intermediates in the reaction of photoenolization of o- methylacetophenone and 1- methylanthraquinone are calculated by semiempirical quantum chemical methods. It is shown that the triplet intermediate is a triplet Σπ- biradical. The geometries and formation enthalpies of the transition state of the thermal transfer of the hydrogen atom in 9- hydroxy- 1,10- anthraquinone- 1- methide are calculated. The MNDO semiempirical methods give significantly exaggerated activation enthalpies in the thermal transfer of the hydrogen atom. Translated fromZhumal Strukturnoi Khimii, Vol. 38, No. 4, pp. 645–654, July–August, 1997.     

醇脱氢酶结构和作用机理研究进展

介绍了醇脱氢酶的种类, 酵母醇脱氢酶和肝醇脱氢酶等两类常用的醇脱氢酶的物理化学性质和活性位点结构. 归纳了对肝醇脱氢酶和酵母醇脱氢酶作用机理的研究, 重点评述了醇脱氢酶催化反应中的两个关键步骤质子转移和氢化物转移过程机理的研究进展.     

Direct hydride transfer in the reaction mechanism of quinoprotein alcohol dehydrogenases: a quantum mechanical investigation

Oxidation of alcohols by direct hydride transfer to the pyrroloquinoline quinone (PQQ) cofactor of quinoprotein alcohol dehydrogenases has been studied using ab initio quantum mechanical methods. Energies and geometries were calculated at the 6-31G(d,p) level of theory. Comparison of the results obtained for PQQ and several derivatives with available structural and spectroscopic data served to judge the feasibility of the calculations. The role of calcium in the enzymatic reaction mechanism has been investigated. Transition state searches have been conducted at the semiempirical and STO-3G(d) level of theory. It is concluded that hydride transfer from the Calpha-position of the substrate alcohol (or aldehyde) directly to the C(5) carbon of PQQ is energetically feasible. Copyright 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1732-1749, 2001