Classical molecular dynamics (MD) simulations and combined quantum mechanics/molecular mechanics (QM/MM) calculations were used to investigate the origin of the enantioselectivity of the Candida antarctica lipase B (CalB) catalyzed O-acetylation of (R,S)-propranolol. The reaction is a two-step process. The initial step is the formation of a reactive acyl enzyme (AcCalB) via a tetrahedral intermediate (TI-1). The stereoselectivity originates from the second step, when AcCalB reacts with the racemic substrate via a second tetrahedral intermediate (TI-2). Reaction barriers for the conversion of (R)- and (S)-propranolol to O-acetylpropranolol were computed for several distinct conformations of TI-2. In QM/MM geometry optimizations and reaction path calculations the QM region was described by density functional theory (B3LYP/TZVP) and the MM region by the CHARMM force field. The QM/MM calculations show that the formation of TI-2 is the rate-determining step. The energy barrier for transformation of (R)-propranolol to O-acetylpropranolol is 4.5 kcal/mol lower than that of the reaction of (S)-propranolol. Enzyme–substrate interactions were identified that play an important role in the enantioselectivity of the reaction. Our QM/MM calculations reproduce and rationalize the experimentally observed enantioselectivity in favor of (R)-propranolol. Furthermore, in contrast to what is commonly suggested for lipase-catalyzed reactions, our results indicate that the tetrahedral intermediate is not a good approximation of the corresponding transition states

Natural sesquiterpene synthases have evolved to make complex terpenoids by quenching reactive carbocations either by proton transfer or by hydroxylation (water capture), depending on their active site. Germacradien-11-ol synthase (Gd11olS) from Streptomyces coelicolor catalyzes the cyclization of farnesyl diphosphate (FDP) into the hydroxylated sesquiterpene germacradien-11-ol. Here, we combine experiment and simulation to guide the redesign of its active site pocket to avoid hydroxylation of the product. Molecular dynamics simulations indicate two regions between which water molecules can flow that are responsible for hydroxylation. Point mutations of selected residues result in variants that predominantly form a complex nonhydroxylated product, which we identify as isolepidozene. Our results indicate how these mutations subtly change the molecular choreography in the Gd11olS active site and thereby pave the way for the engineering of terpene synthases to make complex terpenoid products.

Molecular dynamics (MD) simulations have been performed to study the dynamic behavior of noncovalent enzyme carbocation complexes involved in the cyclization of geranylgeranyl diphosphate to taxadiene catalyzed by taxadiene synthase (TXS). Taxadiene and the observed four side products originate from the deprotonation of carbocation intermediates. The MD simulations of the TXS carbocation complexes provide insights into potential deprotonation mechanisms of such carbocations. The MD results do not support a previous hypothesis that carbocation tumbling is a key factor in the deprotonation of the carbocations by pyrophosphate. Instead water bridges are identified which may allow the formation of side products via multiple proton transfer reactions. A novel reaction path for taxadiene formation is proposed on the basis of the simulations. © 2018 Wiley Periodicals, Inc.

No abstract

Combined quantum mechanics/molecular mechanics (QM/MM) calculations were used to investigate the reaction mechanism of taxadiene synthase (TXS). TXS catalyzes the cyclization of geranylgeranyl diphosphate (GGPP) to taxadiene (T) and four minor cyclic products. All these products originate from the deprotonation of carbocation intermediates. The reaction profiles for the conversion of GGPP to T as well as to minor products were calculated for different configurations of relevant TXS carbocation complexes. The QM region was treated at the M06-2X/TZVP level, while the CHARMM27 force field was used to describe the MM region. The QM/MM calculations suggest a reaction pathway for the conversion of GGPP to T, which slightly differs from previous proposals regarding the number of reaction steps and the conformation of the carbocations. The QM/MM results also indicate that the formation of minor products via water-assisted deprotonation of the carbocations is highly exothermic, by about −7 to −23 kcal/mol. Curiously, however, the computed barriers and reaction energies indicate that the formation of some of the minor products is more facile than the formation of T. Thus, the present QM/MM calculations provide detailed insights into possible reaction pathways and into the origin of the promiscuity of TXS, but they do not reproduce the product distribution observed experimentally.

No abstract

[NiFe]-hydrogenases are highly efficient catalysts for the heterolytic splitting of molecular hydrogen (H2). The heterobimetallic cysteine-coordinated active site of these enzymes is covered by a highly conserved arginine residue, whose role in the reaction is not fully resolved yet. The structural and catalytic role of this arginine is investigated here using QM/MM calculations with various exchange-correlation functionals. All of them give a very consistent picture of the thermodynamics of H2 oxidation. The concept of the presence of a neutral arginine and its direct involvement as a Frustrated Lewis Pair (FLP) in the reaction is critically evaluated. The arginine, however, would exist in its standard protonation state and perform a critical role in positioning and slightly polarizing the substrate H2. It is not directly involved in the heterolytic processing of H2 but guides its approach and reduces its flexibility during binding. Upon substitution of the positively charged arginine by a charge-conserving lysine residue, the H2 binding position remains unaffected. However, critical hydrogen bonding interactions with nearby aspartate residues are lost. In addition, the H2 polarization is unfavorable and the reduced side-chain volume may negatively affect the kinetics of the catalytic process.

Monoterpene synthases are often promiscuous enzymes, yielding product mixtures rather than pure compounds due to the nature of the branched reaction mechanism involving reactive carbocations. Two previously identified bacterial monoterpene synthases, a linalool synthase (bLinS) and a cineole synthase (bCinS), produce nearly pure linalool and cineole from geranyl diphosphate, respectively. We used a combined experimental and computational approach to identify critical residues involved in bacterial monoterpenoid synthesis. Phe77 is essential for bCinS activity, guiding the linear carbocation intermediate towards the formation of the cyclic α‐terpinyl intermediate; removal of the aromatic ring results in variants that produce acyclic products only. Computational chemistry confirmed the importance of Phe77 in carbocation stabilisation. Phe74, Phe78 and Phe179 are involved in maintaining the active site shape in bCinS without a specific role for the aromatic ring. Phe295 in bLinS, and the equivalent Ala301 in bCinS, are essential for linalool and cineole formation, respectively. Where Phe295 places steric constraints on the carbocation intermediates, Ala301 is essential for bCinS initial cyclisation and activity. Our multidisciplinary approach gives unique insights into how carefully placed amino acid residues in the active site can direct carbocations down specific paths, by placing steric constraints or offering stabilisation via cation‐π interactions.