Reaction Free Energy Profile Research Articles

Reaction mechanism and free energy profile for acylation of Candida Antarctica lipase B with methylcaprylate and acetylcholine: Density functional theory calculations

Candida Antarctica lipase B (CALB), a specific enzyme to catalyze the hydrolysis of esters, can be a good candidate for acetylcholine (ACh) hydrolysis instead of acetylcholinesterase. The catalytic mechanism of the CALB acylation, as the first stage in the hydrolysis reaction, with ACh and methylcaprylate (MEC) has been examined by using density functional theory technique. The significant emphasis of this article is on the free energy barriers for the acylation step of hydrolysis reactions. Computed free energy barriers of the first step are 9.2 and 15.9kcalmol−1, but for the second step are 7.9 and 11.6kcalmol−1 for MEC and ACh respectively. Activation free energies are in the comparable and acceptable range and imply both of two reactions are theoretically possible. The stability role of the adjacent amino acids was examined by using two applied tools. It is exposed that the oxyanion hole residues decrease energy barriers by stabilizing the transition state structures.

Role of substrate positioning in the catalytic reaction of 4-hydroxyphenylpyruvate dioxygenase-A QM/MM Study.

Ring hydroxylation and coupled rearrangement reactions catalyzed by 4-hydroxyphenylpyruvate dioxygenase were studied with the QM/MM method ONIOM(B3LYP:AMBER). For electrophilic attack of the ferryl species on the aromatic ring, five channels were considered: attacks on the three ring atoms closest to the oxo ligand (C1, C2, C6) and insertion of oxygen across two bonds formed by them (C1-C2, C1-C6). For the subsequent migration of the carboxymethyl substituent, two possible directions were tested (C1→C2, C1→C6), and two different mechanisms were sought (stepwise radical, single-step heterolytic). In addition, formation of an epoxide (side)product and benzylic hydroxylation, as catalyzed by the closely related hydroxymandelate synthase, were investigated. From the computed reaction free energy profiles it follows that the most likely mechanism of 4-hydroxyphenylpyruvate dioxygenase involves electrophilic attack on the C1 carbon of the ring and subsequent single-step heterolytic migration of the substituent. Computed values of the kinetic isotope effect for this step are inverse, consistent with available experimental data. Electronic structure arguments for the preferred mechanism of attack on the ring are also presented.

Comparative study of water reactivity with Mo₂O(y)⁻ and W₂O(y)⁻ clusters: a combined experimental and theoretical investigation.

A computational investigation of the Mo2O(y)(-) + H2O (y = 4, 5) reactions as well as a photoelectron spectroscopic probe of the deuterated Mo2O6D2(-) product have been carried out to understand a puzzling question from a previous study: Why is the rate constant determined for the Mo2O5(-) + H2O/D2O reaction, the terminal reaction in the sequential oxidation of Mo2O(y)(-) by water, higher than the W2O5(-) + H2O/D2O reaction? This disparity was intriguing because W3O(y)(-) clusters were found to be more reactive toward water than their Mo3O(y)(-) analogs. A comparison of molecular structures reveals that the lowest energy structure of Mo2O5(-) provides a less hindered water addition site than the W2O5(-) ground state structure. Several modes of water addition to the most stable molecular and electronic structures of Mo2O4(-) and Mo2O5(-) were explored computationally. The various modes are discussed and compared with previous computational studies on W2O(y)(-) + H2O reactions. Calculated free energy reaction profiles show lower barriers for the initial Mo2O(y)(-) + H2O addition, consistent with the higher observed rate constant. The terminal Mo2O(y)(-) sequential oxidation product predicted computationally was verified by the anion photoelectron spectrum of Mo2O6D2(-). Based on the computational results, this anion is a trapped dihydroxide intermediate in the Mo2O5(-) + H2O/D2O → Mo2O6(-) + H2/D2 reaction.

Reaction pathways and free energy profiles for cholinesterase-catalyzed hydrolysis of 6-monoacetylmorphine

As the most active metabolite of heroin, 6-monoacetylmorphine (6-MAM) can penetrate into the brain for the rapid onset of heroin effects. The primary enzymes responsible for the metabolism of 6-MAM to the less potent morphine in humans are acetylcholinesterase (AChE) and butyrylcholinesterase (BChE). The detailed reaction pathways for AChE- and BChE-catalyzed hydrolysis of 6-MAM to morphine have been explored, for the first time, in the present study by performing first-principles quantum mechanical/molecular mechanical free energy calculations. It has been demonstrated that the two enzymatic reaction processes follow similar catalytic reaction mechanisms, and the whole catalytic reaction pathway for each enzyme consists of four reaction steps. According to the calculated results, the second reaction step associated with the transition state TS2(a)/TS2(b) should be rate-determining for the AChE/BChE-catalyzed hydrolysis, and the free energy barrier calculated for the AChE-catalyzed hydrolysis (18.3 kcal mol(-1)) is 2.5 kcal mol(-1) lower than that for the BChE-catalyzed hydrolysis (20.8 kcal mol(-1)). The free energy barriers calculated for the AChE- and BChE-catalyzed reactions are in good agreement with the experimentally derived activation free energies (17.5 and 20.7 kcal mol(-1) for the AChE- and BChE-catalyzed reactions, respectively). Further structural analysis reveals that the aromatic residues Phe295 and Phe297 in the acyl pocket of AChE (corresponding to Leu286 and Val288 in BChE) contribute to the lower energy of TS2(a) relative to TS2(b). The obtained structural and mechanistic insights could be valuable for use in future rational design of a novel therapeutic treatment of heroin abuse.

Tests of an Adaptive QM/MM Calculation on Free Energy Profiles of Chemical Reactions in Solution

We present reaction free energy calculations using the adaptive buffered force mixing quantum mechanics/molecular mechanics (bf-QM/MM) method. The bf-QM/MM method combines nonadaptive electrostatic embedding QM/MM calculations with extended and reduced QM regions to calculate accurate forces on all atoms, which can be used in free energy calculation methods that require only the forces and not the energy. We calculate the free energy profiles of two reactions in aqueous solution: the nucleophilic substitution reaction of methyl chloride with a chloride anion and the deprotonation reaction of the tyrosine side chain. We validate the bf-QM/MM method against a full QM simulation, and show that it correctly reproduces both geometrical properties and free energy profiles of the QM model, while the electrostatic embedding QM/MM method using a static QM region comprising only the solute is unable to do so. The bf-QM/MM method is not explicitly dependent on the details of the QM and MM methods, so long as it is possible to compute QM forces in a small region and MM forces in the rest of the system, as in a conventional QM/MM calculation. It is simple, with only a few parameters needed to control the QM calculation sizes, and allows (but does not require) a varying and adapting QM region which is necessary for simulating solutions.

Fundamental Reaction Pathway and Free Energy Profile for Butyrylcholinesterase-Catalyzed Hydrolysis of Heroin

The pharmacological function of heroin requires an activation process that transforms heroin into 6-monoacetylmorphine (6-MAM), which is the most active form. The primary enzyme responsible for this activation process in human plasma is butyrylcholinesterase (BChE). The detailed reaction pathway of the activation process via BChE-catalyzed hydrolysis has been explored computationally, for the first time, in this study via molecular dynamics simulation and first-principles quantum mechanical/molecular mechanical free energy calculations. It has been demonstrated that the whole reaction process includes acylation and deacylation stages. The acylation consists of two reaction steps, i.e., the nucleophilic attack on the carbonyl carbon of the 3-acetyl group of heroin by the hydroxyl oxygen of the Ser198 side chain and the dissociation of 6-MAM. The deacylation also consists of two reaction steps, i.e., the nucleophilic attack on the carbonyl carbon of the acyl-enzyme intermediate by a water molecule and the dissociation of the acetic acid from Ser198. The calculated free energy profile reveals that the second transition state (TS2) should be rate-determining. The structural analysis reveals that the oxyanion hole of BChE plays an important role in the stabilization of rate-determining TS2. The free energy barrier (15.9 ± 0.2 or 16.1 ± 0.2 kcal/mol) calculated for the rate-determining step is in good agreement with the experimentally derived activation free energy (~16.2 kcal/mol), suggesting that the mechanistic insights obtained from this computational study are reliable. The obtained structural and mechanistic insights could be valuable for use in the future rational design of a novel therapeutic treatment of heroin abuse.

Reaction pathway and free energy profile for papain-catalyzed hydrolysis of N-acetyl-Phe-Gly 4-nitroanilide.

Possible reaction pathways for papain-catalyzed hydrolysis of N-acetyl-Phe-Gly 4-nitroanilide (APGNA) have been studied by performing pseudobond first-principles quantum mechanical/molecular mechanical-free energy (QM/MM-FE) calculations. The whole hydrolysis process includes two stages: acylation and deacylation. For the acylation stage of the catalytic reaction, we have explored three possible paths (A, B, and C) and the corresponding free energy profiles along the reaction coordinates. It has been demonstrated that the most favorable reaction path in this stage is path B consisting of two reaction steps: the first step is a proton transfer to form a zwitterionic form (i.e., Cys-S⁻/His-H⁺ ion-pair), and the second step is the nucleophilic attack on the carboxyl carbon of the substrate accompanied by the dissociation of 4-nitroanilide. The deacylation stage includes the nucleophilic attack of a water molecule on the carboxyl carbon of the substrate and dissociation between the carboxyl carbon of the substrate and the sulfhydryl sulfur of Cys25 side chain. The free energy barriers calculated for the acylation and deacylation stages are 20.0 and 10.7 kcal/mol, respectively. Thus, the acylation is rate-limiting. The overall free energy barrier calculated for papain-catalyzed hydrolysis of APGNA is 20.0 kcal/mol, which is reasonably close to the experimentally derived activation free energy of 17.9 kcal/mol.

Sirtuin Deacetylation Mechanism and Catalytic Role of the Dynamic Cofactor Binding Loop.

Sirtuins constitute a novel family of protein deacetylases and play critical roles in epigenetics, cell death, and metabolism. In spite of numerous experimental studies, the key and most complicated stage of its NAD+-dependent catalytic mechanism remains to be elusive. Herein by employing Born-Oppenheimer ab initio QM/MM molecular dynamics simulations, a state-of-the-art computational approach to study enzyme reactions, we have characterized the complete deacetylation mechanism for a sirtuin enzyme, determined its multistep free-energy reaction profile, and elucidated essential catalytic roles of the conserved dynamic cofactor binding loop. These new detailed mechanistic insights would facilitate the design of novel mechanism-based sirtuin modulators.

Reaction pathways and free energy profiles for spontaneous hydrolysis of urea and tetramethylurea: unexpected substituent effects

It has been difficult to directly measure the spontaneous hydrolysis rate of urea and, thus, 1,1,3,3-tetramethylurea (Me4U) was used as a model to determine the "experimental" rate constant for urea hydrolysis. The use of Me4U was based on an assumption that the rate of urea hydrolysis should be 2.8 times that of Me4U hydrolysis because the rate of acetamide hydrolysis is 2.8 times that of N,N-dimethyl-acetamide hydrolysis. The present first-principles electronic-structure calculations on the competing non-enzymatic hydrolysis pathways have demonstrated that the dominant pathway is the neutral hydrolysis via the CN addition for both urea (when pH < ~11.6) and Me4U (regardless of pH), unlike the non-enzymatic hydrolysis of amides where alkaline hydrolysis is dominant. Based on the computational data, the substituent shift of the free energy barrier calculated for the neutral hydrolysis is remarkably different from that for the alkaline hydrolysis, and the rate constant for the urea hydrolysis should be ~1.3 × 10(9)-fold lower than that (4.2 × 10(-12) s(-1)) measured for the Me4U hydrolysis. As a result, the rate enhancement and catalytic proficiency of urease should be 1.2 × 10(25) and 3 × 10(27) M(-1), respectively, suggesting that urease surpasses proteases and all other enzymes in its power to enhance the rate of reaction. All of the computational results are consistent with available experimental data for Me4U, suggesting that the computational prediction for urea is reliable.

Reaction Pathway and Free Energy Profile for Conversion of π-Conjugation Modes in Porphyrin Isomer

Porphycene is a structural isomer of porphyrin with 18π-conjugated aromatic character. Porphycene modified with trifluoromethyl (CF(3)) groups in the periphery of the framework readily affords the isolable 20π-conjugated antiaromatic form through a reaction with a proton-donating reductant. The 20π-conjugated form can be characterized by not only a variety of spectroscopies in solutions but also X-ray crystallography. This paper focuses on the free energy profile in the conversion of the 18π-conjugated porphycene into the 20π-conjugated form. From the results of kinetics, electrochemical measurements, and acid/base titrations, the 20π-conjugated CF(3) porphycene is formed by a concerted proton-electron transfer (CPET) from a hydroquinone reagent to the 18π-conjugated form. The hydrogen-atom affinity of the 18π-conjugated CF(3) porphycene (for two hydrogen atoms) was calculated to be -490 kJ mol(-1), indicating that the N-H bonds in the 20π-conjugated form are rather easily cleaved. This reflects the antiaromatic characteristics of the 20π-conjugated porphycene. We propose that the kinetic and thermochemical analysis using redox potentials and pK(a) data is applicable for determining the reaction pathway in conversion of aromatic/antiaromatic mode of π-conjugated macrocycles as well as popular investigations for oxidations of organic molecules.

Revelation of a catalytic calcium-binding site elucidates unusual metal dependence of a human apyrase.

Human soluble calcium-activated nucleotidase 1 (hSCAN-1) represents a new family of apyrase enzymes that catalyze the hydrolysis of nucleotide di- and triphosphates, thereby modulating extracellular purinergic and pyrimidinergic signaling. Among well-characterized phosphoryl transfer enzymes, hSCAN-1 is unique not only in its unusual calcium-dependent activation, but also in its novel phosphate-binding motif. Its catalytic site does not utilize backbone amide groups to bind phosphate, as in the common P-loop, but contains a large cluster of acidic ionizable side chains. By employing a state-of-the-art computational approach, we have revealed a previously uncharacterized catalytic calcium-binding site in hSCAN-1, which elucidates the unusual calcium-dependence of its apyrase activity. In a high-order coordination shell, the newly identified calcium ion organizes the active site residues to mediate nucleotide binding, to orient the nucleophilic water, and to facilitate the phosphoryl transfer reaction. From ab initio QM/MM molecular dynamics simulations with umbrella sampling, we have characterized a reverse protonation catalytic mechanism for hSCAN-1 and determined its free energy reaction profile. Our results are consistent with available experimental studies and provide new detailed insight into the structure-function relationship of this novel calcium-activated phosphoryl transfer enzyme.

Fundamental Reaction Pathway and Free Energy Profile for Inhibition of Proteasome by Epoxomicin

First-principles quantum mechanical/molecular mechanical free energy calculations have been performed to provide the first detailed computational study on the possible mechanisms for reaction of proteasome with a representative peptide inhibitor, Epoxomicin (EPX). The calculated results reveal that the most favorable reaction pathway consists of five steps. The first is a proton transfer process, activating Thr1-O(γ) directly by Thr1-N(z) to form a zwitterionic intermediate. The next step is nucleophilic attack on the carbonyl carbon of EPX by the negatively charged Thr1-O(γ) atom, followed by a proton transfer from Thr1-N(z) to the carbonyl oxygen of EPX (third step). Then, Thr1-N(z) attacks on the carbon of the epoxide group of EPX, accompanied by the epoxide ring-opening (S(N)2 nucleophilic substitution) such that a zwitterionic morpholino ring is formed between residue Thr1 and EPX. Finally, the product of morpholino ring is generated via another proton transfer. Noteworthy, Thr1-O(γ) can be activated directly by Thr1-N(z) to form the zwitterionic intermediate (with a free energy barrier of only 9.9 kcal/mol), and water cannot assist the rate-determining step, which is remarkably different from the previous perception that a water molecule should mediate the activation process. The fourth reaction step has the highest free energy barrier (23.6 kcal/mol) which is reasonably close to the activation free energy (∼21-22 kcal/mol) derived from experimental kinetic data. The obtained novel mechanistic insights should be valuable for not only future rational design of more efficient proteasome inhibitors but also understanding the general reaction mechanism of proteasome with a peptide or protein.

Reaction pathway and free energy profile determined for specific recognition of oligosaccharide moiety of carboxypeptidase Y

The interaction of mannose specific lectin (from Lens culinaris, LcL) with the carbohydrate moiety of carboxypeptidase Y (CaY) was studied using both atomic force microscope (AFM) and quartz crystal microbalance with dissipation monitoring (QCM-D). The AFM enables to determine the positions of energy barriers present in the energy landscape of the single complex undergoing dissociation. The QCM-D measurements allow the estimation of the quantitative parameters characterizing the kinetics of the studied molecular interaction (namely the association and dissociation rate constants and the association constant). The use of both methods not only delivers the complementary characterization of kinetic and thermodynamic parameters but also permits to investigate the mechanism of the binding and unbinding of the molecules. The results for LcL were compared with those obtained for concanavalin A i.e. lectin, which interacts with the carbohydrate moiety on a similar way.

A Collective Coordinate to Obtain Free Energy Profiles for Complex Reactions in Condensed Phases.

Exploration of chemical reactions in complex explicit environments has become an affordable task with the use of hybrid quantum mechanics/molecular mechanics potentials which allow calculating free energy profiles of chemical reactions under the influence of the surroundings. Tracing these free energy profiles requires the selection of a reaction coordinate, which can be cumbersome for those processes involving more than a single chemical event in a concerted step. We here propose a collective coordinate to be used in the calculation of free energy profiles for complex reactions in condensed phases. This coordinate is based in the definition of the advance along a path introduced by Branduardi et al. (J. Chem. Phys.2007, 126, 054103) but modified to use internal coordinates which are more adequate for the description of chemical reactions. The coordinate is tested with the analysis of the isochorismate transformation to pyruvate and salycilate in aqueous solution and in the active site of PchB, a reaction that involves a CO bond breaking simultaneously with a proton transfer between two carbon atoms. The coordinate introduced here allows obtaining smooth and meaningful free energy profiles of the reaction.

Reaction Pathway and Free Energy Profile for Cocaine Hydrolase-Catalyzed Hydrolysis of (-)-Cocaine.

Reaction pathway of (-)-cocaine hydrolysis catalyzed by our recently discovered most efficient cocaine hydrolase, which is the A199S/F227A/S287G/A328W/Y332G mutant of human butyrylcholinesterase (BChE), and the corresponding free energy profile have been studied by performing first-principles pseudobond quantum mechanical/molecular mechanical (QM/MM)-free energy (FE) calculations. Based on the QM/MM-FE results, the catalytic hydrolysis process consists of four major reaction steps, including the nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester by hydroxyl group of S198, dissociation of (-)-cocaine benzoyl ester, nucleophilic attack on carbonyl carbon of (-)-cocaine benzoyl ester by water, and finally the dissociation between (-)-cocaine benzoyl group and S198 of the enzyme. The second reaction step is rate-determining. The calculated free energy barrier associated with the transition state for the rate-determining step is ~15.0 kcal/mol, which is in excellent agreement with the experimentally-derived activation free energy of ~14.7 kcal/mol. The mechanistic insights obtained from the present study will be valuable for rational design of more active cocaine hydrolase against (-)-cocaine. In particular, future efforts aiming at further increasing the catalytic activity of the enzyme against (-)-cocaine should focus on stabilization of the transition state for the second reaction step in which the benzoyl ester of (-)-cocaine dissociates.

Interaction of 2′-Deoxyadenosine with cis-2-Butene-1,4-dial: Computational Approach to Analysis of Multistep Chemical Reactions

The computational analysis of multistep chemical interactions between 2'-deoxyadenosine and cis-2-butene-1,4-dial has been performed. The applied protocol includes generation of a multistep Gibbs free-energy reaction profile (PCM/M05-2X/6-311+G(d) level) for the transformations of the reagents to products, followed by evaluation of the rate constants, construction of the corresponding kinetic equations, and solving them. Such a procedure allows one to significantly extend the number of experimentally determined steps by addition of the ones computationally predicted. The primary products of the reaction are found to be four diastereomeric adducts characterized by virtually the same stability. The acid-catalyzed dehydration of these adducts leads to a more stable secondary product. Computational verification of UV and NMR spectra has also been performed. It has been revealed that simulated UV and NMR spectra of primary and secondary 2'-deoxyadenosine adducts of cis-2-butene-1,4-dial are in agreement with the experimental observations.

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

Comprehensive DFT Study of the Mechanism of Vanadium-Catalyzed Amination of Benzene with Hydroxylamine

Reaction pathways and free-energy profiles for the conversion of benzene and hydroxylamine to aniline, catalyzed by NaVO3 and VOSO4, in acetic acid/water, were discussed using density functional th...

Fundamental reaction pathway and free energy profile for hydrolysis of intracellular second messenger adenosine 3',5'-cyclic monophosphate (cAMP) catalyzed by phosphodiesterase-4.

As important drug targets for a variety of human diseases, cyclic nucleotide phosphodiesterases (PDEs) are a superfamily of enzymes sharing a similar catalytic site. We have performed pseudobond first-principles quantum mechanical/molecular mechanical-free energy perturbation (QM/MM-FE) and QM/MM-Poisson-Boltzmann surface area (PBSA) calculations to uncover the detailed reaction mechanism for PDE4-catalyzed hydrolysis of adenosine 3',5'-cyclic monophosphate (cAMP). This is the first report on QM/MM reaction-coordinate calculations including the protein environment of any PDE-catalyzed reaction system, demonstrating a unique catalytic reaction mechanism. The QM/MM-FE and QM/MM-PBSA calculations revealed that the PDE4-catalyzed hydrolysis of cAMP consists of two reaction stages: cAMP hydrolysis (stage 1) and bridging hydroxide ion regeneration (stage 2). The stage 1 includes the binding of cAMP in the active site, nucleophilic attack of the bridging hydroxide ion on the phosphorus atom of cAMP, cleavage of O3'-P phosphoesteric bond of cAMP, protonation of the departing O3' atom, and dissociation of hydrolysis product (AMP). The stage 2 includes the binding of solvent water molecules with the metal ions in the active site and regeneration of the bridging hydroxide ion. The dissociation of the hydrolysis product is found to be rate-determining for the enzymatic reaction process. The calculated activation Gibbs free energy of ≥16.0 and reaction free energy of -11.1 kcal/mol are in good agreement with the experimentally derived activation free energy of 16.6 kcal/mol and reaction free energy of -11.5 kcal/mol, suggesting that the catalytic mechanism obtained from this study is reliable and provides a solid base for future rational drug design.

A model SN2 reaction ‘on water’ does not show rate enhancement

Molecular dynamics calculations of the benchmark nucleophilic substitution reaction (SN2) Cl−+CH3Cl are carried out at the water liquid/vapor interface. The reaction free energy profile and the activation free energy are determined as a function of the reactants’ location normal to the surface. The activation free energy remains almost constant relative to that in bulk water, despite the fact that the barrier is expected to significantly decrease as the reaction is carried out near the vapor phase. We show that this is due to the combined effects of a clustering of water molecules around the nucleophile and a relatively weak hydration of the transition state.