Rational design of an enzyme mutant for anti-cocaine therapeutics

(-)-Cocaine is a widely abused drug and there is no available anti-cocaine therapeutic. The disastrous medical and social consequences of cocaine addiction have made the development of an effective pharmacological treatment a high priority. An ideal anti-cocaine medication would be to accelerate (-)-cocaine metabolism producing biologically inactive metabolites. The main metabolic pathway of cocaine in body is the hydrolysis at its benzoyl ester group. Reviewed in this article is the state-of-the-art computational design of high-activity mutants of human butyrylcholinesterase (BChE) against (-)-cocaine. The computational design of BChE mutants have been based on not only the structure of the enzyme, but also the detailed catalytic mechanisms for BChE-catalyzed hydrolysis of (-)-cocaine and (+)-cocaine. Computational studies of the detailed catalytic mechanisms and the structure-and-mechanism-based computational design have been carried out through the combined use of a variety of state-of-the-art techniques of molecular modeling. By using the computational insights into the catalytic mechanisms, a recently developed unique computational design strategy based on the simulation of the rate-determining transition state has been employed to design high-activity mutants of human BChE for hydrolysis of (-)-cocaine, leading to the exciting discovery of BChE mutants with a considerably improved catalytic efficiency against (-)-cocaine. One of the discovered BChE mutants (i.e., A199S/S287G/A328W/Y332G) has a approximately 456-fold improved catalytic efficiency against (-)-cocaine. The encouraging outcome of the computational design and discovery effort demonstrates that the unique computational design approach based on the transition-state simulation is promising for rational enzyme redesign and drug discovery.     

Fundamental reaction mechanism for cocaine hydrolysis in human butyrylcholinesterase

Butyrylcholinesterase (BChE)-cocaine binding and the fundamental pathway for BChE-catalyzed hydrolysis of cocaine have been studied by molecular modeling, molecular dynamics (MD) simulations, and ab initio calculations. Modeling and simulations indicate that the structures of the prereactive BChE/substrate complexes for (-)-cocaine and (+)-cocaine are all similar to that of the corresponding prereactive BChE/butyrylcholine (BCh) complex. The overall binding of BChE with (-)-cocaine and (+)-cocaine is also similar to that proposed with butyrylthiocholine and succinyldithiocholine, i.e., (-)- or (+)-cocaine first slides down the substrate-binding gorge to bind to Trp-82 and stands vertically in the gorge between Asp-70 and Trp-82 (nonprereactive complex) and then rotates to a position in the catalytic site within a favorable distance for nucleophilic attack and hydrolysis by Ser-198 (prereactive complex). In the prereactive complex, cocaine lies horizontally at the bottom of the gorge. The fundamental catalytic hydrolysis pathway, consisting of acylation and deacylation stages similar to those for ester hydrolysis by other serine hydrolases, was proposed on the basis of the simulated prereactive complex and confirmed theoretically by ab initio reaction coordinate calculations. Both the acylation and deacylation follow a double-proton-transfer mechanism. The calculated energetic results show that within the chemical reaction process the highest energy barrier and Gibbs free energy barrier are all associated with the first step of deacylation. The calculated ratio of the rate constant (k(cat)) for the catalytic hydrolysis to that (k(0)) for the spontaneous hydrolysis is approximately 9.0 x 10(7). The estimated k(cat)/k(0) value of approximately 9.0 x 10(7) is in excellent agreement with the experimentally derived k(cat)/k(0) value of approximately 7.2 x 10(7) for (+)-cocaine, whereas it is approximately 2000 times larger than the experimentally derived k(cat)/k(0) value of approximately 4.4 x 10(4) for (-)-cocaine. All of the results suggest that the rate-determining step of the BChE-catalyzed hydrolysis of (+)-cocaine is the first step of deacylation, whereas for (-)-cocaine the change from the nonprereactive complex to the prereactive complex is rate-determining and has a Gibbs free energy barrier higher than that for the first step of deacylation by approximately 4 kcal/mol. A further analysis of the structural changes from the nonprereactive complex to the prereactive complex reveals specific amino acid residues hindering the structural changes, providing initial clues for the rational design of BChE mutants with improved catalytic activity for (-)-cocaine.     

Molecular dynamics simulation of cocaine binding with human butyrylcholinesterase and its mutants

Molecular dynamics (MD) simulations were carried out to study cocaine binding with wild-type human butyrylcholinesterase (BChE) and its mutants based on a recently reported X-ray crystal structure of human BChE. For each BChE-cocaine system, we simulated both the nonprereactive and prereactive complexes in water. Despite the significant difference found at the acyl binding pocket, the simulated structures confirm the fundamental structural and mechanistic insights obtained from earlier computational studies of wild-type BChE with cocaine based on a homology model, e.g. the rate-determining step for BChE-catalyzed hydrolysis of biologically active (-)-cocaine is the (-)-cocaine rotation in the active site from the nonprereactive BChE-(-)-cocaine complex to the prereactive complex. It has been demonstrated that the MD simulations on both the nonprereactive and prereactive BChE-cocaine complexes can clearly reveal whether specific mutations produce the desired BChE-(-)-cocaine binding structures in which the (-)-cocaine rotation is less hindered while the required prereactive BChE-(-)-cocaine binding is maintained. Based on the MD simulations, both A328W/Y332A and A328W/Y332G BChE's are expected to have catalytic activity for (-)-cocaine hydrolysis higher than that of wild-type BChE and the activity of A328W/Y332G BChE should be slightly higher than that of A328W/Y332A BChE due to the less-hindered (-)-cocaine rotation in the mutant BChE's. However, the less-hindered (-)-cocaine rotation is only a necessary condition for a higher activity mutant BChE. The (-)-cocaine rotation is also less hindered in A328W/Y332A/Y419S BChE, but (-)-cocaine binds with A328W/Y332A/Y419S BChE in a way that is not suitable for the catalysis. Thus, A328W/Y332A/Y419S BChE is expected to lose the catalytic activity. The computational predictions were confirmed by our experimental kinetic data, demonstrating that the MD simulation-based computational protocol used in this study is reliable in prediction of the catalytic activity of BChE mutants for (-)-cocaine hydrolysis.     

Free energy perturbation (FEP) simulation on the transition states of cocaine hydrolysis catalyzed by human butyrylcholinesterase and its mutants

A novel computational protocol based on free energy perturbation (FEP) simulations on both the free enzyme and transition state structures has been developed and tested to predict the mutation-caused shift of the free energy change from the free enzyme to the rate-determining transition state for human butyrylcholinesterase (BChE)-catalyzed hydrolysis of (-)-cocaine. The calculated shift, denoted by DeltaDeltaG(1 --> 2), of such kind of free energy change determines the catalytic efficiency (kcat/KM) change caused by the simulated mutation transforming enzyme 1 to enzyme 2. By using the FEP-based computational protocol, the DeltaDeltaG(1 --> 2) values for the mutations A328W/Y332A --> A328W/Y332G and A328W/Y332G --> A328W/Y332G/A199S were calculated to be -0.22 and -1.94 kcal/mol, respectively. The calculated DeltaDeltaG(1 --> 2) values predict that the change from the A328W/Y332A mutant to the A328W/Y332G mutant should slightly improve the catalytic efficiency and that the change from the A328W/Y332G mutant to the A328W/Y332G/A199S mutant should significantly improve the catalytic efficiency of the enzyme for the (-)-cocaine hydrolysis. The predicted catalytic efficiency increases are supported by the experimental data showing that kcat/KM = 8.5 x 10(6), 1.4 x 10(7), and 7.2 x 10(7) min(-1) M(-1) for the A328W/Y332A, A328W/Y332G, and A328W/Y332G/A199S mutants, respectively. The qualitative agreement between the computational and experimental data suggests that the FEP simulations may provide a promising protocol for rational design of high-activity mutants of an enzyme. The general computational strategy of the FEP simulation on a transition state can be used to study the effects of a mutation on the activation free energy for any enzymatic reaction.     

Ab initio model study on acetylcholinesterase catalysis: potential energy surfaces of the proton transfer reactions

Ab initio molecular orbital (MO) and hybrid density functional theory (DFT) calculations have been applied to the initial step of the acylation reaction catalyzed by acetylcholinesterase (AChE), which is the nucleophiric addition of Ser200 in catalytic triads to a neurotransmitter acetylcholine (ACh). We focus our attention mainly on the effects of oxyanion hole and Glu327 on the potential energy surfaces (PESs) for the proton transfer reactions in the catalytic triad Ser200-His440-Glu327. The activation barrier for the addition reaction of Ser200 to ACh was calculated to be 23.4 kcal/mol at the B3LYP/6-31G(d)//HF/3-21G(d) level of theory. The barrier height under the existence of oxyanion hole, namely, Ser200-His440-Glu327-ACh-(oxyanion hole) system, decreased significantly to 14.2 kcal/mol, which is in reasonable agreement with recent experimental value (12.0 kcal/mol). Removal of Glu327 from the catalytic triad caused destabilization of both energy of transition state for the reaction and tetrahedral intermediate (product). PESs calculated for the proton transfer reactions showed that the first proton transfer process is the most important in the stabilization of tetrahedral intermediate complex. The mechanism of addition reaction of ACh was discussed on the basis of theoretical results.     

Role of the catalytic triad and oxyanion hole in acetylcholinesterase catalysis: an ab initio QM/MM study

The initial step of the acylation reaction catalyzed by acetylcholinesterase (AChE) has been studied by a combined ab initio quantum mechanical/molecular mechanical (QM/MM) approach. The reaction proceeds through the nucleophilic addition of the Ser203 O to the carbonyl C of acetylcholine, and the reaction is facilitated by simultaneous proton transfer from Ser203 to His447. The calculated potential energy barrier at the MP2(6-31+G) QM/MM level is 10.5 kcal/mol, consistent with the experimental reaction rate. The third residue of the catalytic triad, Glu334, is found to be essential in stabilizing the transition state through electrostatic interactions. The oxyanion hole, formed by peptidic NH groups from Gly121, Gly122, and Ala204, is also found to play an important role in catalysis. Our calculations indicate that, in the AChE-ACh Michaelis complex, only two hydrogen bonds are formed between the carbonyl oxygen of ACh and the peptidic NH groups of Gly121 and Gly122. As the reaction proceeds, the distance between the carbonyl oxygen of ACh and NH group of Ala204 becomes smaller, and the third hydrogen bond is formed both in the transition state and in the tetrahedral intermediate.     

Ab initio and density functional theory studies of the catalytic mechanism for ester hydrolysis in serine hydrolases

We present results from ab initio and density functional theory studies of the mechanism for serine hydrolase catalyzed ester hydrolysis. A model system containing both the catalytic triad and the oxyanion hole was studied. The catalytic triad was represented by formate anion, imidazole, and methanol. The oxyanion hole was represented by two water molecules. Methyl formate was used as the substrate. In the acylation step, our computations show that the cooperation of the Asp group and oxyanion hydrogen bonds is capable of lowering the activation barrier by about 15 kcal/mol. The transition state leading to the first tetrahedral intermediate in the acylation step is rate limiting with an activation barrier (ΔE₀) of 13.4 kcal/mol. The activation barrier in the deacylation step is smaller. The double-proton-transfer mechanism is energetically unfavorable by about 2 kcal/mol. The bonds between the Asp group and the His group, and the hydrogen bonds in the oxyanion hole, increase in strength going from the Michaelis complex toward the transition state and the tetrahedral intermediate. In the acylation step, the tetrahedral intermediate is a very shallow minimum on the energy surface and is not viable when molecular vibrations are included. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 69: 89–103, 1998     

Structure-and-mechanism-based design and discovery of therapeutics for cocaine overdose and addiction

(-)-Cocaine is a widely abused drug and there is currently no available anti-cocaine therapeutic. Promising agents, such as anti-cocaine catalytic antibodies and high-activity mutants of human butyrylcholinesterase (BChE), for therapeutic treatment of cocaine overdose have been developed through structure-and-mechanism-based design and discovery. In particular, a unique computational design strategy based on the modeling and simulation of the rate-determining transition state has been developed and used to design and discover desirable high-activity mutants of BChE. One of the discovered high-activity mutants of BChE has a approximately 456-fold improved catalytic efficiency against (-)-cocaine. The encouraging outcome of the structure-and-mechanism-based design and discovery effort demonstrates that the unique computational design approach based on transition state modeling and simulation is promising for rational enzyme redesign and drug discovery. The general approach of the structure-and-mechanism-based design and discovery may be used to design high-activity mutants of any enzyme or catalytic antibody.     

Studying enzyme binding specificity in acetylcholinesterase using a combined molecular dynamics and multiple docking approach

A combined molecular dynamics simulation and multiple ligand docking approach is applied to study the binding specificity of acetylcholinesterase (AChE) with its natural substrate acetylcholine (ACh), a family of substrate analogues, and choline. Calculated docking energies are well correlated to experimental k(cat)/K(M) values, as well as to experimental binding affinities of a related series of TMTFA inhibitors. The "esteratic" and "anionic" subsites are found to act together to achieve substrate binding specificity. We find that the presence of ACh in the active site of AChE not only stabilizes the setup of the catalytic triad but also tightens both subsites to achieve better binding. The docking energy gained from this induced fit is 0.7 kcal/mol for ACh. For the binding of the substrate tailgroup to the anionic subsite, both the size and the positive charge of the tailgroup are important. The removal of the positive charge leads to a weaker binding of 1 kcal/mol loss in docking energy. Substituting each tail methyl group with hydrogen results in both an incremental loss in docking energy and also a decrease in the percentage of structures docked in the active site correctly set up for catalysis.     

Transition state stabilization and substrate strain in enzyme catalysis: ab initio QM/MM modelling of the chorismate mutase reaction

To investigate fundamental features of enzyme catalysis, there is a need for high-level calculations capable of modelling crucial, unstable species such as transition states as they are formed within enzymes. We have modelled an important model enzyme reaction, the Claisen rearrangement of chorismate to prephenate in chorismate mutase, by combined ab initio quantum mechanics/molecular mechanics (QM/MM) methods. The best estimates of the potential energy barrier in the enzyme are 7.4-11.0 kcal mol(-1)(MP2/6-31+G(d)//6-31G(d)/CHARMM22) and 12.7-16.1 kcal mol(-1)(B3LYP/6-311+G(2d,p)//6-31G(d)/CHARMM22), comparable to the experimental estimate of Delta H(++)= 12.7 +/- 0.4 kcal mol(-1). The results provide unequivocal evidence of transition state (TS) stabilization by the enzyme, with contributions from residues Arg90, Arg7, and Arg63. Glu78 stabilizes the prephenate product (relative to substrate), and can also stabilize the TS. Examination of the same pathway in solution (with a variety of continuum models), at the same ab initio levels, allows comparison of the catalyzed and uncatalyzed reactions. Calculated barriers in solution are 28.0 kcal mol(-1)(MP2/6-31+G(d)/PCM) and 24.6 kcal mol(-1)(B3LYP/6-311+G(2d,p)/PCM), comparable to the experimental finding of Delta G(++)= 25.4 kcal mol(-1) and consistent with the experimentally-deduced 10(6)-fold rate acceleration by the enzyme. The substrate is found to be significantly distorted in the enzyme, adopting a structure closer to the transition state, although the degree of compression is less than predicted by lower-level calculations. This apparent substrate strain, or compression, is potentially also catalytically relevant. Solution calculations, however, suggest that the catalytic contribution of this compression may be relatively small. Consideration of the same reaction pathway in solution and in the enzyme, involving reaction from a 'near-attack conformer' of the substrate, indicates that adoption of this conformation is not in itself a major contribution to catalysis. Transition state stabilization (by electrostatic interactions, including hydrogen bonds) is found to be central to catalysis by the enzyme. Several hydrogen bonds are observed to shorten at the TS. The active site is clearly complementary to the transition state for the reaction, stabilizing it more than the substrate, so reducing the barrier to reaction.     

Effect of electrostatic interaction on the mechanism of dehalogenation catalyzed by haloalkane dehalogenase

Theoretical calculation has been carried out for the nucleophilic displacement reaction of 1,2‐dichloroethane catalyzed by haloalkane dehalogenase. The results indicate that different hydrogen bond patterns of the oxyanion hole and the halide‐stabilizing residues play an important role in the dehalogenation reaction. They cause concertedly an earlier transition state (TS) with the activation barrier of 16.60 kcal/mol. The stabilization effect of Trp125 and Trp175 on chlorine atom in the TS is larger than that of the reactant complex by 15.67 kcal/mol so that, they make contribution to the stabilization of the TS. Moreover, the reaction shows the enzymatic action can be attributed to a combination of reactant‐state destabilization and transition‐state electrostatic stabilization. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Molecular dynamics simulations of the catalytic pathway of a cysteine protease: a combined QM/MM study of human cathepsin K

Molecular dynamics simulations using a combined QM/MM potential have been performed to study the catalytic mechanism of human cathepsin K, a member of the papain family of cysteine proteases. We have determined the two-dimensional free energy surfaces of both acylation and deacylation steps to characterize the reaction mechanism. These free energy profiles show that the acylation step is rate limiting with a barrier height of 19.8 kcal/mol in human cathepsin K and of 29.3 kcal/mol in aqueous solution. The free energy of activation for the deacylation step is 16.7 kcal/mol in cathepsin K and 17.8 kcal/mol in aqueous solution. The reduction of free energy barrier is achieved by stabilization of the oxyanion in the transition state. Interestingly, although the "oxyanion hole" has been formed in the Michaelis complex, the amide units do not donate hydrogen bonds directly to the carbonyl oxygen of the substrate, but they stabilize the thiolate anion nucleophile. Hydrogen-bonding interactions are induced as the substrate amide group approaches the nucleophile, moving more than 2 A and placing the oxyanion in contact with Gln19 and the backbone amide of Cys25. The hydrolysis of peptide substrate shares a common mechanism both for the catalyzed reaction in human cathepsin K and for the uncatalyzed reaction in water. Overall, the nucleophilic attack by Cys25 thiolate and the proton-transfer reaction from His162 to the amide nitrogen are highly coupled, whereas a tetrahedral intermediate is formed along the nucleophilic reaction pathway.     

Modeling the catalysis of anti-cocaine catalytic antibody: competing reaction pathways and free energy barriers

The competing reaction pathways and the corresponding free energy barriers for cocaine hydrolysis catalyzed by an anti-cocaine catalytic antibody, mAb15A10, were studied by using a novel computational strategy based on the binding free energy calculations on the antibody binding with cocaine and transition states. The calculated binding free energies were used to evaluate the free energy barrier shift from the cocaine hydrolysis in water to the antibody-catalyzed cocaine hydrolysis for each reaction pathway. The free energy barriers for the antibody-catalyzed cocaine hydrolysis were predicted to be the corresponding free energy barriers for the cocaine hydrolysis in water plus the calculated free energy barrier shifts. The calculated free energy barrier shift of -6.87 kcal/mol from the dominant reaction pathway of the cocaine benzoyl ester hydrolysis in water to the dominant reaction pathway of the antibody-catalyzed cocaine hydrolysis is in good agreement with the experimentally derived free energy barrier shift of -5.93 kcal/mol. The calculated mutation-caused shifts of the free energy barrier are also reasonably close to the available experimental activity data. The good agreement suggests that the protocol for calculating the free energy barrier shift from the cocaine hydrolysis in water to the antibody-catalyzed cocaine hydrolysis may be used in future rational design of possible high-activity mutants of the antibody as anti-cocaine therapeutics. The general strategy of the free energy barrier shift calculation may also be valuable in studying a variety of chemical reactions catalyzed by other antibodies or proteins through noncovalent bonding interactions with the substrates.     

Theoretical study on the reaction pathways of HFCO + H2O

For the reaction of methanoyl fluoride with water, both optimized structures and vibrational wavenumbers of reaction intermediates, transition structures and product complexes were calculated and characterized with theory at the MP2/6-311++G(d,p) level. Including a catalytic path and concerted and stepwise hydrolysis paths, possible reaction mechanisms were also investigated. The catalytic reaction of HFCO yielding HF and CO has the smallest activation barrier, 29.6 kcal/mol, whereas for the concerted hydrolysis 33.0 kcal/mol is required to overcome the barrier to form transoid HCOOH + HF, which is less than for the stepwise counterpart, 42.0 kcal/mol.     

First-principle studies of intermolecular and intramolecular catalysis of protonated cocaine

We have performed a series of first-principles electronic structure calculations to examine the reaction pathways and the corresponding free energy barriers for the ester hydrolysis of protonated cocaine in its chair and boat conformations. The calculated free energy barriers for the benzoyl ester hydrolysis of protonated chair cocaine are close to the corresponding barriers calculated for the benzoyl ester hydrolysis of neutral cocaine. However, the free energy barrier calculated for the methyl ester hydrolysis of protonated cocaine in its chair conformation is significantly lower than for the methyl ester hydrolysis of neutral cocaine and for the dominant pathway of the benzoyl ester hydrolysis of protonated cocaine. The significant decrease of the free energy barrier, approximately 4 kcal/mol, is attributed to the intramolecular acid catalysis of the methyl ester hydrolysis of protonated cocaine, because the transition state structure is stabilized by the strong hydrogen bond between the carbonyl oxygen of the methyl ester moiety and the protonated tropane N. The relative magnitudes of the free energy barriers calculated for different pathways of the ester hydrolysis of protonated chair cocaine are consistent with the experimental kinetic data for cocaine hydrolysis under physiologic conditions. Similar intramolecular acid catalysis also occurs for the benzoyl ester hydrolysis of (protonated) boat cocaine in the physiologic condition, although the contribution of the intramolecular hydrogen bonding to transition state stabilization is negligible. Nonetheless, the predictability of the intramolecular hydrogen bonding could be useful in generating antibody-based catalysts that recruit cocaine to the boat conformation and an analog that elicited antibodies to approximate the protonated tropane N and the benzoyl O more closely than the natural boat conformer might increase the contribution from hydrogen bonding. Such a stable analog of the transition state for intramolecular catalysis of cocaine benzoyl-ester hydrolysis was synthesized and used to successfully elicit a number of anticocaine catalytic antibodies.     

Synthesis,enzymes inhibitory properties and characterization of 2- (bis (4-aminophenyl) methyl) butan-1-ol compound: Quantum simulations,and in-silico molecular docking studies

In this study to be presented, the BAPMB molecule was synthesized and structurally characterized. All calculations were applied using the detailed DFT/B3LYP method with 6-311G (d, p) and SDD depended on the stable phase geometry of the molecule. Also, various HOMO-LUMO energy gaps, inter-orbital intramolecular interactions of the natural bond, and electro-static surface mapping actions were also realized. The molecule was characterized by NMR spectroscopic analysis. Besides, LC/MS data were acquired. In this analysis, fragment ions (m/z 168.9) of BAPMB were obtained as 137.8, 179.9, 198, 201.5, and 228.0 approximately. Also, molecular docking was performed while examining the exact binding site and binding mechanism of the ligand on the protein. In the study, glide scores in binding affinity with BAPMB – AChE, BAPMB – BChE, and BAPMB – GST, respectively; It was found to be −7.228 ​cal/mol, −7.205 ​cal/mol, −6.07 ​cal/mol and BAPMB - AChE was found to be more effective with receptor binding score. BAPMB was analyzed for its inhibition features counter AChE, BChE, and GST enzymes that exhibit effective inhibition. AChE, BChE, GST enzymes were powerfully inhibited by BAPMB. BChE showed excellent activity, particularly in comparison to standard tacrine. Eventually, AIM analysis was performed to search intermolecular interactions in the BAPMB compound.     

Theoretical study of oxidation of cyclohexane diol to adipic anhydride by [Ru(IV)(O)(tpa)(H2O)]2+ complex (tpa ═ tris(2-pyridylmethyl)amine)

The catalytic conversion of 1,2-cyclohexanediol to adipic anhydride by Ru(IV)O(tpa) (tpa ═ tris(2-pyridylmethyl)amine) is discussed using density functional theory calculations. The whole reaction is divided into three steps: (1) formation of α-hydroxy cyclohexanone by dehydrogenation of cyclohexanediol, (2) formation of 1,2-cyclohexanedione by dehydrogenation of α-hydroxy cyclohexanone, and (3) formation of adipic anhydride by oxygenation of cyclohexanedione. In each step the two-electron oxidation is performed by Ru(IV)O(tpa) active species, which is reduced to bis-aqua Ru(II)(tpa) complex. The Ru(II) complex is reactivated using Ce(IV) and water as an oxygen source. There are two different pathways of the first two steps of the conversion depending on whether the direct H-atom abstraction occurs on a C-H bond or on its adjacent oxygen O-H. In the first step, the C-H (O-H) bond dissociation occurs in TS1 (TS2-1) with an activation barrier of 21.4 (21.6) kcal/mol, which is followed by abstraction of another hydrogen with the spin transition in both pathways. The second process also bifurcates into two reaction pathways. TS3 (TS4-1) is leading to dissociation of the C-H (O-H) bond, and the activation barrier of TS3 (TS4-1) is 20.2 (20.7) kcal/mol. In the third step, oxo ligand attack on the carbonyl carbon and hydrogen migration from the water ligand occur via TS5 with an activation barrier of 17.4 kcal/mol leading to a stable tetrahedral intermediate in a triplet state. However, the slightly higher energy singlet state of this tetrahedral intermediate is unstable; therefore, a spin crossover spontaneously transforms the tetrahedral intermediate into a dione complex by a hydrogen rebound and a C-C bond cleavage. Kinetic isotope effects (k(H)/k(D)) for the electronic processes of the C-H bond dissociations calculated to be 4.9-7.4 at 300 K are in good agreement with experiment values of 2.8-9.0.     

Theoretical studies on the elimination of hydrogen fluoride from alkyl fluoride and its substituent effect

The reaction mechanism of the α, α and α, β elimination of hydrogen fluorides from alkyl fluorides has been studied theoretically. For fluoroethane as a reactant, the transition state (TS) optimized at the level of the 6-31G** basis set shows that the α, β elimination proceeds via a four membered-ring TS with a barrier height 64.6 kcal/mol, while the α, α elimination, via a three-membered ring TS with a 83.7 kcal/mol barrier. Four substituents, CH₃, CN, F, and NH₂, were used to investigate the substituent effect of elimination by using the 3-21G basis set. The calculated barriers show that NH₂-substituted alkyl fluorides favor both the α, α and α, β elimination and these two reactions would be expected to proceed simultaneously. © 1996 John Wiley & Sons, Inc.     

Synergism of catalysis and reaction center rehybridization in nucleophilic additions to cumulenes: the one-, two-, three-water hydrolyses of carbodiimide and methyleneimine

The results of a theoretical study of the one-, two- and three-water hydrolyses of carbodiimide and the one- and two-water hydrolyses of methyleneimine are presented. All structures were optimized and characterized at the MP2(full)/6-31G* level of theory. Energies for the one-water hydrolysis of carbodiimide were determined at numerous higher levels of theory, up to the QCISD(T)(fc)/6-311+G(3df,2p)//MP2(full)/6-31G* level. The Delta E0(Delta G298) activation barriers for the rate-determining steps of the one-, two- and three-water hydrolyses of carbodiimide, respectively, are 44.8 (46.3), 29.3 (32.3) and 22.9 (26.2) kcal mol(-1) at the MP2(full)/6-31G* level. The consideration of a second water molecule catalyzes the hydrolysis by 15.5 kcal mol(-1) on the E0 surface and by 14.0 kcal mol(-1) on the G298 surface with respect to the one-water hydrolysis. Placement of a third water molecule opposite the site of proton transfer catalyzes the reaction by an additional 6.4 kcal mol(-1) on the E0 surface and by 6.1 kcal mol(-1) on the G298 surface. The catalytic effect of the third water molecule results from the synergistic effects of rehybridization and charge relaxation in the transition state. The charge relaxation in the transition state is illustrated through natural population analysis calculations on the pre-coordination complexes and the transition state structures. We also consider the placement of the third water molecule in the proton transfer chain and we show this to be of little catalytic relevance. The activation barriers determined for the one- and two-water hydrolyses of methyleneimine are Delta G298=51.9 and Delta G298=35.5 kcal mol(-1), respectively, and they are larger than for carbodiimide. The results are compared with the hydrolyses of carbon dioxide and formaldehyde.     

Theoretical Study on the Dehydrogenation Reaction of H2S by VS+ (3∑-)

The dehydrogenation reaction of H2S by the 3∑- ground state of VS+: VS+ + H2S → VS2+ + H2 has been studied by using Density Functional Theory (DFT) at the B3LYP/DZVP level. It is found that the reaction proceeds along two possible pathways (A and B) yielding two isomer dehydrogenation products VS2+-1 (3B2) and VS2+-2 (3A1), respectively. For both pathways,the reaction has a two-step-reaction mechanism that involves the migration of two hydrogen atoms from S2 to V+, respectively. The migration of the second hydrogen via TS3 and that of the first via TS4 are the rate-determining steps for pathways A and B, respectively. The activation energy is 17.4 kcal/mol for pathway A and 22.8 kcal/mol for pathway B relative to the reactants. The calculated reaction heat of 9.9 kcal/mol indicates the endothermicity of pathway A and that of -11.9 kcal/mol suggests the exothermicity of pathway B.