Modeling evolution of hydrogen bonding and stabilization of transition states in the process of cocaine hydrolysis catalyzed by human butyrylcholinesterase

Molecular dynamics (MD) simulations and quantum mechanical/molecular mechanical (QM/MM) calculations were performed on the prereactive enzyme-substrate complex, transition states, intermediates, and product involved in the process of human butyrylcholinesterase (BChE)-catalyzed hydrolysis of (-)-cocaine. The computational results consistently reveal a unique role of the oxyanion hole (consisting of G116, G117, and A199) in BChE-catalyzed hydrolysis of cocaine, compared to acetylcholinesterase (AChE)-catalyzed hydrolysis of acetylcholine. During BChE-catalyzed hydrolysis of cocaine, only G117 has a hydrogen bond with the carbonyl oxygen (O31) of the cocaine benzoyl ester in the prereactive BChE-cocaine complex, and the NH groups of G117 and A199 are hydrogen-bonded with O31 of cocaine in all of the transition states and intermediates. Surprisingly, the NH hydrogen of G116 forms an unexpected hydrogen bond with the carboxyl group of E197 side chain and, therefore, is not available to form a hydrogen bond with O31 of cocaine in the acylation. The NH hydrogen of G116 is only partially available to form a weak hydrogen bond with O31 of cocaine in some structures involved in the deacylation. The change of the estimated hydrogen-bonding energy between the oxyanion hole and O31 of cocaine during the reaction process demonstrates how the protein environment can affect the energy barrier for each step of the BChE-catalyzed hydrolysis of cocaine. These insights concerning the effects of the oxyanion hole on the energy barriers provide valuable clues on how to rationally design BChE mutants with a higher catalytic activity for the hydrolysis of (-)-cocaine.     

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).     

Predicted Michaelis-Menten complexes of cocaine-butyrylcholinesterase. Engineering effective butyrylcholinesterase mutants for cocaine detoxication

Butyrylcholinesterase (BChE) is important in cocaine metabolism, but it hydrolyzes (-)-cocaine only one-two thousandth as fast as the unnatural (+)-stereoisomer. A starting point in engineering BChE mutants that rapidly clear cocaine from the bloodstream, for overdose treatment, is to elucidate structural factors underlying the stereochemical difference in catalysis. Here, we report two three-dimensional Michaelis-Menten complexes of BChE liganded with natural and unnatural cocaine molecules, respectively, that were derived from molecular modeling and supported by experimental studies. Such complexes revealed that the benzoic ester group of both cocaine stereoisomers must rotate toward the catalytic Ser(198) for hydrolysis. Rotation of (-)-cocaine appears to be hindered by interactions of its phenyl ring with Phe(329) and Trp(430). These interactions do not occur with (+)-cocaine. Because the rate of (-)-cocaine hydrolysis is predicted to be determined mainly by the re-orientation step, it should not be greatly influenced by pH. In fact, measured rates of this reaction were nearly constant over the pH range from 5.5 to 8.5, despite large rate changes in hydrolysis of (+)-cocaine. Our models can explain why BChE hydrolyzes (+)-cocaine faster than (-)-cocaine, and they suggest that mutations of certain residues in the catalytic site could greatly improve catalytic efficiency and the potential for detoxication.     

Design of high-activity mutants of human butyrylcholinesterase against (-)-cocaine: structural and energetic factors affecting the catalytic efficiency

The present study was aimed to explore the correlation between the protein structure and catalytic efficiency of butyrylcholinesterase (BChE) mutants against (-)-cocaine by modeling the rate-determining transition state (TS1), i.e., the transition state for the first step of chemical reaction process, of (-)-cocaine hydrolysis catalyzed by various mutants of human BChE in comparison with the wild type. Molecular modeling of the TS1 structures revealed that mutations on certain nonactive site residues can indirectly affect the catalytic efficiency of the enzyme against (-)-cocaine through enhancing or weakening the overall hydrogen bonding between the carbonyl oxygen of (-)-cocaine benzoyl ester and the oxyanion hole of the enzyme. Computational insights and predictions were supported by the catalytic activity data obtained from wet experimental tests on the mutants of human BChE, including five new mutants reported for the first time. The BChE mutants with at least ～1000-fold improved catalytic efficiency against (-)-cocaine compared to the wild-type BChE are all associated with the TS1 structures having stronger overall hydrogen bonding between the carbonyl oxygen of (-)-cocaine benzoyl ester and the oxyanion hole of the enzyme. The combined computational and experimental data demonstrate a reasonable correlation relationship between the hydrogen-bonding distances in the TS1 structure and the catalytic efficiency of the enzyme against (-)-cocaine.     

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.     

Comparison of a QM/MM force field and molecular mechanics force fields in simulations of alanine and glycine "dipeptides" (Ace-Ala-Nme and Ace-Gly-Nme) in water in relation to the problem of modeling the unfolded peptide backbone in solution

We compare the conformational distributions of Ace-Ala-Nme and Ace-Gly-Nme sampled in long simulations with several molecular mechanics (MM) force fields and with a fast combined quantum mechanics/molecular mechanics (QM/MM) force field, in which the solute's intramolecular energy and forces are calculated with the self-consistent charge density functional tight binding method (SCCDFTB), and the solvent is represented by either one of the well-known SPC and TIP3P models. All MM force fields give two main states for Ace-Ala-Nme, beta and alpha separated by free energy barriers, but the ratio in which these are sampled varies by a factor of 30, from a high in favor of beta of 6 to a low of 1/5. The frequency of transitions between states is particularly low with the amber and charmm force fields, for which the distributions are noticeably narrower, and the energy barriers between states higher. The lower of the two barriers lies between alpha and beta at values of psi near 0 for all MM simulations except for charmm22. The results of the QM/MM simulations vary less with the choice of MM force field; the ratio beta/alpha varies between 1.5 and 2.2, the easy pass lies at psi near 0, and transitions between states are more frequent than for amber and charmm, but less frequent than for cedar. For Ace-Gly-Nme, all force fields locate a diffuse stable region around phi = pi and psi = pi, whereas the amber force field gives two additional densely sampled states near phi = +/-100 degrees and psi = 0, which are also found with the QM/MM force field. For both solutes, the distribution from the QM/MM simulation shows greater similarity with the distribution in high-resolution protein structures than is the case for any of the MM simulations.     

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.     

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.     

Coherent singlet-triplet transitions at the initial step of tubulin assembly into microtubules

Quantum chemistry calculations [DFT-B3LYP QM/MM method, 6-31G** basis set, + ab initio molecular dynamics] were used to study the action of Mg2+ on tubulin properties. It was shown that the hydration of the guanosine triphosphate-tubulin forms a protein zone structure, which includes a electron-occupied zone and a conductivity zone. The binding of Mg2+ to guanosine triphosphate-tubulin results in the unpairing of electrons in the occupied zone (triplet state formation) followed by their transition to the conductivity zone in which the inversion of spin occurs (singlet state formation). The formation of triplet state is the initial step in the subsequent protein dynamics in the picosecond range of time. The dynamics shows up as a coherent oscillating transition of tubulin between the triplet and singlet states, which is evidence of a simultaneous adjustment between nuclear and electron configurations of the protein (ab initio molecular dynamics calculations). The barrier between the triplet and singlet states does not exceed 0.60 kcal x mol(-1). The barrier overcome is considered as electron tunneling through the Fermi surface, which separates the occupied and conductivity zones. Zone formation occurs in the presence of the shell of biological water surrounding the protein.     

Theoretical insights into the protonation states of active site cysteine and citrullination mechanism of Porphyromonas gingivalis peptidylarginine deiminase

Porphyromonas gingivalis peptidylarginine deiminase (PPAD) catalyzes the citrullination of peptidylarginine, which plays a critical role in the rheumatoid arthritis (RA) and gene regulation. For a better understanding of citrullination mechanism of PPAD, it is required to establish the protonation states of active site cysteine, which is still a controversial issue for the members of guanidino‐group‐modifying enzyme superfamily. In this work, we first explored the transformation between the two states: State N (both C351 and H236 are neutral) and State I (both residues exist as a thiolate–imidazolium ion pair), and then investigated the citrullination reaction of peptidylarginine, using a combined QM/MM approach. State N is calculated to be more stable than State I by 8.46 kcal/mol, and State N can transform to State I via two steps of substrate‐assisted proton transfer. Citrullination of the peptidylarginine contains deamination and hydrolysis. Starting from State N, the deamination reaction corresponds to an energy barrier of 18.82 kcal/mol. The deprotonated C351 initiates the nucleophilic attack to the substrate, which is the key step for deamination reaction. The hydrolysis reaction contains two chemical steps. Both the deprotonated D238 and H236 can act as the bases to activate the hydrolytic water, which correspond to similar energy barriers (～17 kcal/mol). On the basis of our calculations, C351, D238, and H236 constitute a catalytic triad, and their protonation states are critical for both the deamination and hydrolysis processes. In view of the sequence similarity, these findings may be shared with human PAD1–PAD4 and other guanidino‐group‐modifying enzymes. Proteins 2017; 85:1518–1528. © 2017 Wiley Periodicals, Inc.     

Quantum chemical modeling of the GTP hydrolysis by the RAS-GAP protein complex

We present results of the modeling for the hydrolysis reaction of guanosine triphosphate (GTP) in the RAS-GAP protein complex using essentially ab initio quantum chemistry methods. One of the approaches considers a supermolecular cluster composed of 150 atoms at a consistent quantum level. Another is a hybrid QM/MM method based on the effective fragment potential technique, which describes interactions between quantum and molecular mechanical subsystems at the ab initio level of the theory. Our results show that the GTP hydrolysis in the RAS-GAP protein complex can be modeled by a substrate-assisted catalytic mechanism. We can locate a configuration on the top of the barrier corresponding to the transition state of the hydrolysis reaction such that the straightforward descents from this point lead either to reactants GTP+H(2)O or to products guanosine diphosphate (GDP)+H(2)PO(4)(-). However, in all calculations such a single-step process is characterized by an activation barrier that is too high. Another possibility is a two-step reaction consistent with formation of an intermediate. Here the Pgamma-O(Pbeta) bond is already broken, but the lytic water molecule is still in the pre-reactive state. We present arguments favoring the assumption that the first step of the GTP hydrolysis reaction in the RAS-GAP protein complex may be assigned to the breaking of the Pgamma-O(Pbeta) bond prior to the creation of the inorganic phosphate.     

QM/MM modeling the Ras-GAP catalyzed hydrolysis of guanosine triphosphate

The mechanism of the hydrolysis reaction of guanosine triphosphate (GTP) by the protein complex Ras-GAP (p21(ras) - p120(GAP)) has been modeled by the quantum mechanical-molecular mechanical (QM/MM) and ab initio quantum calculations. Initial geometry configurations have been prompted by atomic coordinates of a structural analog (PDBID:1WQ1). It is shown that the minimum energy reaction path is consistent with an assumption of two-step chemical transformations. At the first stage, a unified motion of Arg789 of GAP, Gln61, Thr35 of Ras, and the lytic water molecule results in a substantial spatial separation of the gamma-phosphate group of GTP from the rest of the molecule (GDP). This phase of hydrolysis process proceeds through the low-barrier transition state TS1. At the second stage, Gln61 abstracts and releases protons within the subsystem including Gln61, the lytic water molecule and the gamma-phosphate group of GTP through the corresponding transition state TS2. Direct quantum calculations show that, in this particular environment, the reaction GTP + H(2)O --> GDP + H(2)PO(4) (-) can proceed with reasonable activation barriers of less than 15 kcal/mol at every stage. This conclusion leads to a better understanding of the anticatalytic effect of cancer-causing mutations of Ras, which has been debated in recent years.     

Characterization of a catalytic ligand bridging metal ions in phosphodiesterases 4 and 5 by molecular dynamics simulations and hybrid quantum mechanical/molecular mechanical calculations

Cyclic nucleotide phosphodiesterases (PDEs) constitute a large superfamily of enzymes regulating concentrations of intracellular second messengers cAMP and cGMP through PDE-catalyzed hydrolysis. Although three-dimensional x-ray crystal structures of PDE4 and PDE5 have been reported, it is uncertain whether a critical, second bridging ligand (BL2) in the active site is H2O or HO- because hydrogen atoms cannot be determined by x-ray diffraction. The identity of BL2 is theoretically determined by performing molecular dynamics simulations and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations, for the first time, on the protein structures resolved by x-ray diffraction. The computational results confirm our previous suggestion, which was based on QM calculations on a simplified active site model, that BL2 in PDE4 should be HO-, rather than H2O, serving as the nucleophile to initialize the catalytic hydrolysis of cAMP. The molecular dynamics simulations and QM/MM calculations on PDE5 demonstrate for the first time that the BL2 in PDE5 should also be HO- rather than H2O as proposed in recently published reports on the x-ray crystal structures, which serves as the nucleophile to initialize the PDE5-catalyzed hydrolysis of cGMP. These fundamental structural insights provide a rational basis for future structure-based drug design targeting PDEs.     

Binding free energy calculation with QM/MM hybrid methods for Abl-Kinase inhibitor

We report a Quantum mechanics/Molecular Mechanics–Poisson-Boltzmann/ Surface Area (QM/MM-PB/SA) method to calculate the binding free energy of c-Abl human tyrosine kinase by combining the QM and MM principles where the ligand is treated quantum mechanically and the rest of the receptor by classical molecular mechanics. To study the role of entropy and the flexibility of the protein ligand complex in a solvated environment, molecular dynamics calculations are performed using a hybrid QM/MM approach. This work shows that the results of the QM/MM approach are strongly correlated with the binding affinity. The QM/MM interaction energy in our reported study confirms the importance of electronic and polarization contributions, which are often neglected in classical MM-PB/SA calculations. Moreover, a comparison of semi-empirical methods like DFTB-SCC, PM3, MNDO, MNDO-PDDG, and PDDG-PM3 is also performed. The results of the study show that the implementation of a DFTB-SCC semi-empirical Hamiltonian that is derived from DFT gives better results than other methods. We have performed such studies using the AMBER molecular dynamic package for the first time. The calculated binding free energy is also in agreement with the experimentally determined binding affinity for c-Abl tyrosine kinase complex with Imatinib.     

Substrate polarization in enzyme catalysis: QM/MM analysis of the effect of oxaloacetate polarization on acetyl-CoA enolization in citrate synthase

Citrate synthase is an archetypal carbon-carbon bond forming enzyme. It promotes the conversion of oxaloacetate (OAA) to citrate by catalyzing the deprotonation (enolization) of acetyl-CoA, followed by nucleophilic attack of the enolate form of this substrate on OAA to form a citryl-CoA intermediate and subsequent hydrolysis. OAA is strongly bound to the active site and its alpha-carbonyl group is polarized. This polarization has been demonstrated spectroscopically, [(Kurz et al., Biochemistry 1985;24:452-457; Kurz and Drysdale, Biochemistry 1987;26:2623-2627)] and has been suggested to be an important catalytic strategy. Substrate polarization is believed to be important in many enzymes. The first step, formation of the acetyl-CoA enolate intermediate, is thought to be rate-limiting in the mesophilic (pig/chicken) enzyme. We have examined the effects of substrate polarization on this key step using quantum mechanical/molecular mechanical (QM/MM) methods. Free energy profiles have been calculated by AM1/CHARMM27 umbrella sampling molecular dynamics (MD) simulations, together with potential energy profiles. To study the influence of OAA polarization, profiles were calculated with different polarization of the OAA alpha-carbonyl group. The results indicate that OAA polarization influences catalysis only marginally but has a larger effect on intermediate stabilization. Different levels of treatment of OAA are compared (MM or QM), and its polarization in the protein and in water analyzed at the B3LYP/6-31+G(d)/CHARMM27 level. Analysis of stabilization by individual residues shows that the enzyme mainly stabilizes the enolate intermediate (not the transition state) through electrostatic (including hydrogen bond) interactions: these contribute much more than polarization of OAA.     

QM/MM studies on the catalytic mechanism of Phenylethanolamine N-methyltransferase

Epinephrine is a naturally occurring adrenomedullary hormone that transduces environmental stressors into cardiovascular actions. As the only route in the catecholamine biosynthetic pathway, Phenylethanolamine N-methyltransferase (PNMT) catalyzes the synthesis of epinephrine. To elucidate the detailed mechanism of enzymatic catalysis of PNMT, combined quantum-mechanical/molecular-mechanical (QM/MM) calculations were performed. The calculation results reveal that this catalysis contains three elementary steps: the deprotonation of protonated norepinphrine, the methyl transferring step and deprotonation of the methylated norepinphrine. The methyl transferring step was proved to be the rate-determining step undergoing a SN2 mechanism with an energy barrier of 16.4kcal/mol. During the whole catalysis, two glutamic acids Glu185 and Glu219 were proved to be loaded with different effects according to the calculations results of the mutants. These calculation results can be used to explain the experimental observations and make a good complementarity for the previous QM study.     

QM/MM studies on the catalytic mechanism of Phenylethanolamine N-methyltransferase

Epinephrine is a naturally occurring adrenomedullary hormone that transduces environmental stressors into cardiovascular actions. As the only route in the catecholamine biosynthetic pathway, Phenylethanolamine N-methyltransferase (PNMT) catalyzes the synthesis of epinephrine. To elucidate the detailed mechanism of enzymatic catalysis of PNMT, combined quantum-mechanical/molecular-mechanical (QM/MM) calculations were performed. The calculation results reveal that this catalysis contains three elementary steps: the deprotonation of protonated norepinphrine, the methyl transferring step and deprotonation of the methylated norepinphrine. The methyl transferring step was proved to be the rate-determining step undergoing a SN2 mechanism with an energy barrier of 16.4 kcal/mol. During the whole catalysis, two glutamic acids Glu185 and Glu219 were proved to be loaded with different effects according to the calculations results of the mutants. These calculation results can be used to explain the experimental observations and make a good complementarity for the previous QM study.     

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.     

Computational and experimental studies on the catalytic mechanism of biliverdin-IXbeta reductase

BVR-B (biliverdin-IXbeta reductase) also known as FR (flavin reductase) is a promiscuous enzyme catalysing the pyridine-nucleotide-dependent reduction of a variety of flavins, biliverdins, PQQ (pyrroloquinoline quinone) and ferric ion. Mechanistically it is a good model for BVR-A (biliverdin-IXalpha reductase), a potential pharmacological target for neonatal jaundice and also a potential target for adjunct therapy to maintain protective levels of biliverdin-IXalpha during organ transplantation. In a commentary on the structure of BVR-B it was noted that one outstanding issue remained: whether the mechanism was a concerted hydride transfer followed by protonation of a pyrrolic anion or protonation of the pyrrole followed by hydride transfer. In the present study we have attempted to address this question using QM/MM (quantum mechanics/molecular mechanics) calculations. QM/MM potential energy surfaces show that the lowest energy pathway proceeds with a positively charged pyrrole intermediate via two transition states. These initial calculations were performed with His(153) as the source of the proton. However site-directed mutagenesis studies with both the H153A and the H153N mutant reveal that His(153) is not required for catalytic activity. We have repeated the calculation with a solvent hydroxonium donor and obtain a similar energy landscape indicating that protonation of the pyrrole is the most likely first step followed by hydride transfer and that the required proton may come from bulk solvent. The implications of the present study for the design of inhibitors of BVR-A are discussed.     

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.