Macrophomate Synthase: QM/MM Simulations Address the Diels−Alder versus Michael−Aldol Reaction Mechanism

Recent advances in enzymatic carbon-carbon bond formation

Enzymatic carbon-carbon (C-C) bond formation reactions have become an effective and invaluable tool for designing new biological and medicinal molecules, often with asymmetric features. This review provides a systematic overview of key C-C bond formation reactions and enzymes, with the focus of reaction mechanisms and recent advances. These reactions include the aldol reaction, Henry reaction, Knoevenagel condensation, Michael addition, Friedel-Crafts alkylation and acylation, Mannich reaction, Morita-Baylis-Hillman (MBH) reaction, Diels-Alder reaction, acyloin condensations via Thiamine Diphosphate (ThDP)-dependent enzymes, oxidative and reductive C-C bond formation, C-C bond formation through C1 resource utilization, radical enzymes for C-C bond formation, and other C-C bond formation reactions.

Iron-sulphur protein catalysed [4+2] cycloadditions in natural product biosynthesis

To the best of our knowledge, enzymes that catalyse intramolecular Diels-Alder ([4+2] cycloaddition) reactions are frequently reported in natural product biosynthesis; however, no native enzymes utilising Lewis acid catalysis have been reported. Verticilactam is a representative member of polycyclic macrolactams, presumably produced by spontaneous cycloaddition. We report that the intramolecular [4+2] cycloadditions can be significantly accelerated by ferredoxins (Fds), a class of small iron-sulphur (Fe-S) proteins. Through iron atom substitution by Lewis acidic gallium (Ga) iron and computational calculations, we confirm that the ubiquitous Fe-S cluster efficiently functions as Lewis acid to accelerate the tandem [4+2] cycloaddition and Michael addition reactions by lowering free energy barriers. Our work highlights Nature's ingenious strategy to generate complex molecule structures using the ubiquitous Fe-S protein. Furthermore, our study sheds light on the future design of Fd as a versatile Lewis acid catalyst for [4+2] cycloaddition reactions.

Computational insights into the binding modes, keto-enol tautomerization and stereo-electronically controlled decarboxylation of oxaloacetate in the active site of macrophomate synthase

Oxaloacetic acid (OAA) is a β-ketocarboxylic acid, which plays an important role as an intermediate in some metabolic pathways, including the tricarboxylic acid cycle, gluconeogenesis and fatty acid biosynthesis. Animal studies have indicated that supplementing oxaloacetic acid shows an increase of lifespan and other substantial health benefits including mitochondrial DNA protection, and protection of retinal, neural and pancreatic tissues. Most of the chemical transformations of OAA in the metabolic pathways have been extensively studied; however, the understanding of decarboxylation of OAA at the atomic level is relatively lacking. Here, we carried out MD simulations and combined quantum mechanical/molecular mechanical (QM/MM) calculations as an example to systematically elucidate the binding modes, keto-enol tautomerization and decarboxylation of OAA in the active site of macrophomate synthase (MPS), which is a Mg(II)-dependent bifunctional enzyme that catalyzes both the decarboxylation of OAA and [4+2] cycloaddition of 2-pyrone with the decarboxylated intermediate of OAA (pyruvate enolate). On the basis of our calculations, it was found that the Mg2+-coordinated oxaloacetate may exist in enol forms and keto forms. The four keto forms can be transformed into each other by simply rotating the C2-C3 single bond, nevertheless, the keto-enol tautomerization strictly requires the assistance of pocket water molecules. In addition, the decarboxylation is stereo-electronically controlled, i.e., it is the relative orientation of the terminal carboxyl anion that determines the rate of decarboxylation. As such, the chemistry of oxaloacetate in the active site of MPS is complex. On one hand, the most stable binding mode (K-I) may undergo enol-keto tautomerization to isomerize to the enol form, which may further react with the second substrate; on the other hand, K-I may isomerize to another binding mode K-II to proceed decarboxylation to generate pyruvate enolate and CO2. Starting from K-I, the enol-keto tautomerization corresponds to a barrier of 16.2 kcal mol-1, whereas the decarboxylation is associated with an overall barrier of 19.7 kcal mol-1. These findings may provide useful information for understanding the chemistry of OAA and the catalysis of related enzymes, and they are basically in agreement with the available experimental kinetic data.

Interrogation of an Enzyme Library Reveals the Catalytic Plasticity of Naturally Evolved [4+2] Cyclases

Stereoselective carbon-carbon bond forming reactions are quintessential transformations in organic synthesis. One example is the Diels-Alder reaction, a [4+2] cycloaddition between a conjugated diene and a dienophile to form cyclohexenes. The development of biocatalysts for this reaction is paramount for unlocking sustainable routes to a plethora of important molecules. To obtain a comprehensive understanding of naturally evolved [4+2] cyclases, and to identify hitherto uncharacterised biocatalysts for this reaction, we constructed a library comprising forty-five enzymes with reported or predicted [4+2] cycloaddition activity. Thirty-one library members were successfully produced in recombinant form. In vitro assays employing a synthetic substrate incorporating a diene and a dienophile revealed broad-ranging cycloaddition activity amongst these polypeptides. The hypothetical protein Cyc15 was found to catalyse an intramolecular cycloaddition to generate a novel spirotetronate. The crystal structure of this enzyme, along with docking studies, establishes the basis for stereoselectivity in Cyc15, as compared to other spirotetronate cyclases.

Mechanism of the Stereoselective Catalysis of Diels-Alderase PyrE3 Involved in Pyrroindomycin Biosynthesis

The biosynthesis of pyrroindomycins A and B features a complexity-building [4 + 2] cycloaddition cascade, which generates the spirotetramate core under the catalytic effects of monofunctional Diels-Alderases PyrE3 and PyrI4. We recently showed that the main functions of PyrI4 include acid catalysis and induced-fit/conformational selection. We now present quantum mechanical and molecular dynamics studies implicating a different mode of action by PyrE3, which prearranges an anionic polyene substrate into a high-energy reactive conformation at which an inverse-electron-demand Diels-Alder reaction can occur with a low barrier. Stereoselection is realized by strong binding interactions at the endo stereochemical relationship and a local steric constraint on the endo-1,3-diene unit. These findings, illustrating distinct mechanisms for PyrE3 and PyrI4, highlight how nature has evolved multiple ways to catalyze Diels-Alder reactions.

Discovery and investigation of natural Diels-Alderases

It has been proposed that biosyntheses of many natural products involve pericyclic reactions, including Diels-Alder (DA) reaction. However, only a small set of enzymes have been proposed to catalyze pericyclic reactions. Most surprisingly, there has been no formal identification of natural enzymes that can be defined to catalyze DA reactions (DAases), despite the wide application of the reaction in chemical syntheses of complex organic compounds. However, recent studies began to accumulate a growing body of evidence that supports the notion that enzymes that formally catalyze DA reactions, in fact exist. In this review, I will begin by describing a short history behind the discovery and characterization of macrophomate synthase, one of the earliest enzymes that was proposed to catalyze an intermolecular DA reaction during the biosynthesis of a substituted benzoic acid in a phytopathogenic fungus Macrophoma commelinae. Then, I will discuss representative enzymes that have been chemically authenticated to catalyze DA reactions, with emphasis on more recent discoveries of DAases involved mainly in fungal secondary metabolite biosynthesis except for one example from a marine streptomycete. The current success in identification of a series of DAases and enzymes that catalyze other pericyclic reactions owes to the combined efforts from both the experimental and theoretical approaches in discovering natural products. Such efforts typically involve identifying the chemical features derived from cycloaddition reactions, isolating the biosynthetic genes that encode enzymes that generate such chemical features and deciphering the reaction mechanisms for the enzyme-catalyzed pericyclic reactions.

Diverse evolutionary origins of microbial [4 + 2]-cyclases in natural product biosynthesis

Natural [4 + 2]-cyclases catalyze concerted cycloaddition during biosynthesis of over 400 natural products reported. Microbial [4 + 2]-cyclases are structurally diverse with a broad range of substrates. Thus far, about 52 putative microbial [4 + 2]-cyclases of 13 different types have been characterized, with over 20 crystal structures. However, how these cyclases have evolved during natural product biosynthesis remains elusive. Structural and phylogenetic analyses suggest that these different types of [4 + 2]-cyclases might have diverse evolutionary origins, such as reductases, dehydratases, methyltransferases, oxidases, etc. Divergent evolution of enzyme function might have occurred in these different families. Understanding the independent evolutionary history of these cyclases would provide new insights into their catalysis mechanisms and the biocatalyst design.

Enzymatic dimerization in the biosynthetic pathway of microbial natural products

Covering: up to August 2020The dramatic increase in the identification of dimeric natural products generated by microorganisms and plants has played a significant role in drug discovery. The biosynthetic pathways of these products feature inherent dimerization reactions, which are valuable for biosynthetic applications and chemical transformations. The extraordinary mechanisms of the dimerization of secondary metabolites should advance our understanding of the uncommon chemical rules for natural product biosynthesis, which will, in turn, accelerate the discovery of dimeric reactions and molecules in nature and provide promising strategies for the total synthesis of natural products through dimerization. This review focuses on the enzymes involved in the dimerization in the biosynthetic pathway of microbial natural products, with an emphasis on cytochrome P450s, laccases, and intermolecular [4 + 2] cyclases, along with other atypical enzymes. The identification, characterization, and catalytic landscapes of these enzymes are also introduced.

Mechanism of 3-Methylglutaconyl CoA Decarboxylase AibA/AibB: Pericyclic Reaction versus Direct Decarboxylation

The enzyme 3-methylglutaconyl coenzyme A (CoA) decarboxylase (called AibA/AibB) catalyzes the decarboxylation of 3-methylglutaconyl CoA to generate 3,3-dimethylacrylyl-CoA, representing an important step in the biosynthesis of isovaleryl-coenzyme A in Myxococcus xanthus when the regular pathway is blocked. A novel mechanism involving a pericyclic transition state has previously been proposed for this enzyme, making AibA/AibB unique among decarboxylases. Herein, density functional calculations are used to examine the energetic feasibility of this mechanism. It is shown that the intramolecular pericyclic reaction is associated with a very high energy barrier that is similar to the barrier of the same reaction in the absence of the enzyme. Instead, the calculations show that a direct decarboxylation mechanism has feasible energy barriers that are in line with the experimental observations.

Artificial Enzymes for Diels-Alder Reactions

The Diels-Alder (DA) reaction is a cycloaddition of a conjugated diene and an alkene (dienophile) leading to the formation of a cyclohexene derivative through a concerted mechanism. As DA reactions generally proceed with a high degree of regio- and stereoselectivity, they are widely used in synthetic organic chemistry. Considering eco-conscious public and governmental movements, efforts are now directed towards the development of synthetic processes that meet environmental concerns. Artificial enzymes, which can be developed to catalyze abiotic reactions, appear to be important synthetic tools in the synthetic biology field. This review describes the different strategies used to develop protein-based artificial enzymes for DA reactions, including for in cellulo approaches.

The expanding world of biosynthetic pericyclases: cooperation of experiment and theory for discovery

Covering: 2000 to 2018 Pericyclic reactions are a distinct class of reactions that have wide synthetic utility. Before the recent discoveries described in this review, enzyme-catalyzed pericyclic reactions were not widely known to be involved in biosynthesis. This situation is changing rapidly. We define the scope of pericyclic reactions, give a historical account of their discoveries as biosynthetic reactions, and provide evidence that there are many enzymes in nature that catalyze pericyclic reactions. These enzymes, the "pericyclases," are the subject of this review.

Crystal structure of the putative cyclase IdmH from the indanomycin nonribosomal peptide synthase/polyketide synthase

Indanomycin is biosynthesized by a hybrid nonribosomal peptide synthase/polyketide synthase (NRPS/PKS) followed by a number of 'tailoring' steps to form the two ring systems that are present in the mature product. It had previously been hypothesized that the indane ring of indanomycin was formed by the action of IdmH using a Diels-Alder reaction. Here, the crystal structure of a selenomethionine-labelled truncated form of IdmH (IdmH-Δ99-107) was solved using single-wavelength anomalous dispersion (SAD) phasing. This truncated variant allows consistent and easy crystallization, but importantly the structure was used as a search model in molecular replacement, allowing the full-length IdmH structure to be determined to 2.7 Å resolution. IdmH is a homodimer, with the individual protomers consisting of an α+β barrel. Each protomer contains a deep hydrophobic pocket which is proposed to constitute the active site of the enzyme. To investigate the reaction catalysed by IdmH, 88% of the backbone NMR resonances were assigned, and using chemical shift perturbation of [15N]-labelled IdmH it was demonstrated that indanomycin binds in the active-site pocket. Finally, combined quantum mechanical/molecular mechanical (QM/MM) modelling of the IdmH reaction shows that the active site of the enzyme provides an appropriate environment to promote indane-ring formation, supporting the assignment of IdmH as the key Diels-Alderase catalysing the final step in the biosynthesis of indanomycin through a similar mechanism to other recently characterized Diels-Alderases involved in polyketide-tailoring reactions. An animated Interactive 3D Complement (I3DC) is available in Proteopedia at https://proteopedia.org/w/Journal:IUCrJ:S2052252519012399.

Computational Insights into an Enzyme-Catalyzed [4+2] Cycloaddition

The enzyme SpnF, involved in the biosynthesis of spinosyn A, catalyzes a formal [4+2] cycloaddition of a 22-membered macrolactone, which may proceed as a concerted [4+2] Diels–Alder reaction or a stepwise [6+4] cycloaddition followed by a Cope rearrangement. Quantum mechanics/molecular mechanics (QM/MM) calculations combined with free energy simulations show that the Diels–Alder pathway is favored in the enzyme environment. OM2/CHARMM free energy simulations for the SpnF-catalyzed reaction predict a free energy barrier of 22 kcal/mol for the concerted Diels–Alder process and provide no evidence of a competitive stepwise pathway. Compared with the gas phase, the enzyme lowers the Diels–Alder barrier significantly, consistent with experimental observations. Inspection of the optimized geometries indicates that the enzyme may prearrange the substrate within the active site to accelerate the [4+2] cycloaddition and impede the [6+4] cycloaddition through interactions with active-site residues. Judging from partial charge analysis, we find that the hydrogen bond between the Thr196 residue of SpnF and the substrate C15 carbonyl group contributes to the enhancement of the rate of the Diels–Alder reaction. QM/MM simulations show that the substrate can easily adopt a reactive conformation in the active site of SpnF because interconversion between the C5–C6 s-trans and s-cis conformers is facile. Our QM/MM study suggests that the enzyme SpnF does behave as a Diels-Alderase.

Computational Design of Protein Function

The computational design of enzymes is a tremendous challenge for both chemistry and biochemistry. The ability to design stable and functional biocatalysts that could operate under different conditions to perform chemical reactions without precedent in nature, allowing the large-scale production of chemicals à la carte, would revolutionise both synthetic, pharmacologic and materials chemistry. Despite the great advances achieved, this highly multidisciplinary area of research is still in its infancy. This chapter describes the ‘inside-out’ protocol for computational enzyme design and both the achievements and limitations of the current technology are highlighted. Furthermore, molecular dynamics simulations have proved to be invaluable in the enzyme design process, constituting an important tool for discovering elusive catalytically relevant conformations of the engineered or designed enzyme. As a complement to the ‘inside-out’ design protocol, different examples where hybrid QM/MM approaches have been directly applied to discover beneficial mutations in rational computational enzyme design are described.

Natural [4 + 2]-Cyclases

[4 + 2]-Cycloadditions are increasingly being recognized in the biosynthetic pathways of many structurally complex natural products. A relatively small collection of enzymes from these pathways have been demonstrated to increase rates of cyclization and impose stereochemical constraints on the reactions. While mechanistic investigation of these enzymes is just beginning, recent studies have provided new insights with implications for understanding their biosynthetic roles, mechanisms of catalysis, and evolutionary origin.

Biogenesis of Triterpene Dimers from Orthoquinones Related to Quinonemethides: Theoretical Study on the Reaction Mechanism

The biogenetic origin of triterpene dimers from the Celastraceae family has been proposed as assisted hetero-Diels-Alder reaction (HDA). In this work, computational calculation of HDA between natural quinonemethides (tingenone and isopristimerol) and hypothetical orthoquinones has been performed at the M06-2X/6-31G(d) level of theory. We have located all the HDA transition states supporting the biogenetic route via HDA cycloadditions. We found that all reactions take place through a concerted inverse electron demand and asynchronous mechanism. The enzymatic assistance for dimer formation was analyzed in terms of the calculated transition state energy barrier.

A Counterion-Directed Approach to the Diels-Alder Paradigm: Cascade Synthesis of Tricyclic Fused Cyclopropanes

An approach to the intramolecular Diels–Alder reaction has led to a cascade synthesis of complex carbocycles composed of three fused rings and up to five stereocenters with complete stereocontrol. Computational analysis reveals that the reaction proceeds by a Michael/Michael/cyclopropanation/epimerization cascade in which size and coordination of the counterion is key.

Natural Diels-Alderases: Elusive and Irresistable

Eight examples of biosynthetic pathways wherein a natural enzyme has been identified and claimed to function as a catalyst for the [4 + 2] cycloaddition reaction, namely, Diels-Alderases, are briefly reviewed. These are discussed in the context of the mechanistic challenges associated with the technical difficulty of proving that the net formal [4 + 2] cycloaddition under study indeed proceeds through a synchronous mechanism and that the putative biosynthetic enzyme deploys the pericyclic transition state required for a Diels-Alder cycloaddition reaction.

Application of a BOSS-Gaussian interface for QM/MM simulations of Henry and methyl transfer reactions

Hybrid quantum mechanics and molecular mechanics (QM/MM) computer simulations have become an indispensable tool for studying chemical and biological phenomena for systems too large to treat with QM alone. For several decades, semiempirical QM methods have been used in QM/MM simulations. However, with increased computational resources, the introduction of ab initio and density function methods into on-the-fly QM/MM simulations is being increasingly preferred. This adaptation can be accomplished with a program interface that tethers independent QM and MM software packages. This report introduces such an interface for the BOSS and Gaussian programs, featuring modification of BOSS to request QM energies and partial atomic charges from Gaussian. A customizable C-shell linker script facilitates the interprogram communication. The BOSS-Gaussian interface also provides convenient access to Charge Model 5 (CM5) partial atomic charges for multiple purposes including QM/MM studies of reactions. In this report, the BOSS-Gaussian interface is applied to a nitroaldol (Henry) reaction and two methyl transfer reactions in aqueous solution. Improved agreement with experiment is found by determining free-energy surfaces with MP2/CM5 QM/MM simulations than previously reported investigations using semiempirical methods.

Effective use of heterologous hosts for characterization of biosynthetic enzymes allows production of natural products and promotes new natural product discovery

In the past few years, there has been impressive progress in elucidating the mechanism of biosynthesis of various natural products accomplished through the use of genetic, molecular biological and biochemical techniques. Here, we present a comprehensive overview of the current results from our studies on fungal natural product biosynthetic enzymes, including nonribosomal peptide synthetase and polyketide synthase-nonribosomal peptide synthetase hybrid synthetase, as well as auxiliary enzymes, such as methyltransferases and oxygenases. Specifically, biosynthesis of the following compounds is described in detail: (i) Sch210972, potentially involving a Diels-Alder reaction that may be catalyzed by CghA, a functionally unknown protein identified by targeted gene disruption in the wild type fungus; (ii) chaetoglobosin A, formed via multi-step oxidations catalyzed by three redox enzymes, one flavin-containing monooxygenase and two cytochrome P450 oxygenases as characterized by in vivo biotransformation of relevant intermediates in our engineered Saccharomyces cerevisiae; (iii) (-)-ditryptophenaline, formed by a cytochrome P450, revealing the dimerization mechanism for the biosynthesis of diketopiperazine alkaloids; (iv) pseurotins, whose variations in the C- and O-methylations and the degree of oxidation are introduced combinatorially by multiple redox enzymes; and (v) spirotryprostatins, whose spiro-carbon moiety is formed by a flavin-containing monooxygenase or a cytochrome P450 as determined by heterologous de novo production of the biosynthetic intermediates and final products in Aspergillus niger. We close our discussion by summarizing some of the key techniques that have facilitated the discovery of new natural products, production of their analogs and identification of biosynthetic mechanisms in our study.

The structure of SpnF, a standalone enzyme that catalyzes [4 + 2] cycloaddition

In the biosynthetic pathway of the spinosyn insecticides, the tailoring enzyme SpnF performs a [4 + 2] cycloaddition on a 22-membered macrolactone to forge an embedded cyclohexene ring. To learn more about this reaction, which could potentially proceed through a Diels-Alder mechanism, we determined the 1.50-Å-resolution crystal structure of SpnF bound to S-adenosylhomocysteine. This sets the stage for advanced experimental and computational studies to determine the precise mechanism of SpnF-mediated cyclization.

Diradical ring closing reactions displaying Woodward–Hoffmann behaviour and torquoselectivity

The applicability of the Woodward–Hoffmann rules seems to expand beyond classical closed shell structures.

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.

Fundamental reaction pathway for peptide metabolism by proteasome: insights from first-principles quantum mechanical/molecular mechanical free energy calculations

Proteasome is the major component of the crucial non-lysosomal protein degradation pathway in the cells, but the detailed reaction pathway is unclear. In this study, first-principles quantum mechanical/molecular mechanical free energy calculations have been performed to explore, for the first time, possible reaction pathways for proteasomal proteolysis/hydrolysis of a representative peptide, succinyl-leucyl-leucyl-valyl-tyrosyl-7-amino-4-methylcoumarin (Suc-LLVY-AMC). The computational results reveal that the most favorable reaction pathway consists of six steps. The first is a water-assisted proton transfer within proteasome, activating Thr1-O(γ). The second is a nucleophilic attack on the carbonyl carbon of a Tyr residue of substrate by the negatively charged Thr1-O(γ), followed by the dissociation of the amine AMC (third step). The fourth step is a nucleophilic attack on the carbonyl carbon of the Tyr residue of substrate by a water molecule, accompanied by a proton transfer from the water molecule to Thr1-N(z). Then, Suc-LLVY is dissociated (fifth step), and Thr1 is regenerated via a direct proton transfer from Thr1-N(z) to Thr1-O(γ). According to the calculated energetic results, the overall reaction energy barrier of the proteasomal hydrolysis is associated with the transition state (TS3(b)) for the third step involving a water-assisted proton transfer. The determined most favorable reaction pathway and the rate-determining step have provided a reasonable interpretation of the reported experimental observations concerning the substituent and isotopic effects on the kinetics. The calculated overall free energy barrier of 18.2 kcal/mol is close to the experimentally derived activation free energy of ∼18.3-19.4 kcal/mol, suggesting that the computational results are reasonable.

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.

Hybrid dynamics simulation engine for metalloproteins

Quality computational description of metalloproteins is a great challenge due to the vast span of time- and lengthscales characteristic of their existence. We present an efficient new method that allows for robust characterization of metalloproteins. It combines quantum mechanical (QM) description of the metal-containing active site, and extensive dynamics of the protein captured by discrete molecular dynamics (DMD) (QM/DMD). DMD samples the entire protein, including the backbone, and most of the active site, except for the immediate coordination region of the metal. QM operates on the part of the protein of electronic and chemical significance, which may include tens to hundreds of atoms. The breathing quantum-classical boundary provides a continuous mutual feedback between the two machineries. We test QM/DMD using the Fe-containing electron transporter protein, rubredoxin, and its three mutants as a model. QM/DMD can provide a reliable balanced description of metalloproteins' structure, dynamics, and electronic structure in a reasonable amount of time. As an illustration of QM/DMD capabilities, we then predict the structure of the Ca(2+) form of the enzyme catechol O-methyl transferase, which, unlike the native Mg(2+) form, is catalytically inactive. The Mg(2+) site is ochtahedral, but the Ca(2+) is 7-coordinate and features the misalignment of the reacting parts of the system. The change is facilitated by the backbone adjustment. QM/DMD is ideal and fast for providing this level of structural insight.

Computational Enzyme Design: Advances, hurdles and possible ways forward

This mini review addresses recent developments in computational enzyme design. Successful protocols as well as known issues and limitations are discussed from an energetic perspective. It will be argued that improved results can be obtained by including a dynamic treatment in the design protocol. Finally, a molecular dynamics-based approach for evaluating and refining computational designs is presented.

Computational design of a Diels-Alderase from a thermophilic esterase: the importance of dynamics

A novel computational Diels-Alderase design, based on a relatively rare form of carboxylesterase from Geobacillus stearothermophilus, is presented and theoretically evaluated. The structure was found by mining the PDB for a suitable oxyanion hole-containing structure, followed by a combinatorial approach to find suitable substrates and rational mutations. Four lead designs were selected and thoroughly modeled to obtain realistic estimates of substrate binding and prearrangement. Molecular dynamics simulations and DFT calculations were used to optimize and estimate binding affinity and activation energies. A large quantum chemical model was used to capture the salient interactions in the crucial transition state (TS). Our quantitative estimation of kinetic parameters was validated against four experimentally characterized Diels-Alderases with good results. The final designs in this work are predicted to have rate enhancements of ≈ 10(3)-10(6) and high predicted proficiencies. This work emphasizes the importance of considering protein dynamics in the design approach, and provides a quantitative estimate of the how the TS stabilization observed in most de novo and redesigned enzymes is decreased compared to a minimal, 'ideal' model. The presented design is highly interesting for further optimization and applications since it is based on a thermophilic enzyme (T (opt) = 70 °C).

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.

The Intramolecular Diels-Alder Reaction of Tryptamine-Derived Zincke Aldehydes Is a Stepwise Process

Computational studies show that the base-mediated intramolecular Diels-Alder of tryptamine-derived Zincke aldehydes, used as a key step in the synthesis of the Strychnos alkaloids norfluorocurarine and strychnine, proceeds via a stepwise pathway. The experimentally determined importance of a potassium counterion in the base is explained by its ability to preorganize the Zincke aldehyde diene in an s-cis conformation suitable to bicyclization. Computation also supports the thermodynamic importance of the generation of a stable enolate in the final reaction step. The thermal cycloreversion reaction of the Diels-Alder products is also found to proceed in a stepwise manner.

Current developments and challenges in the search for a naturally selected Diels-Alderase

Only a very few examples of enzymes known to catalyze pericyclic reactions have been reported, and presently no enzyme has been demonstrated unequivocally to catalyze a Diels-Alder reaction. Nevertheless, research into secondary metabolism has led to the discovery of numerous natural products exhibiting the structural hallmarks of [4+2] cycloadditions, prompting efforts to characterize the responsible enzymatic processes. These efforts have resulted in a growing collection of enzymes believed to catalyze pericyclic [4+2] cycloaddition reactions; however, in each case the complexity of the substrates and catalytic properties of these enzymes poses significant challenges in substantiating these hypotheses. Herein we consider the principles motivating these efforts and the enzymological systems currently under investigation.

Asymmetric total syntheses of (+)- and (-)-versicolamide B and biosynthetic implications

The Diels-Alder reaction is one of the most well-studied, synthetically useful organic transformations. While a significant number of naturally occurring substances are postulated to arise by biosynthetic Diels-Alder reactions, rigorous confirmation of a mechanistically distinct natural Diels-Alderase enzyme remains elusive. Within this context, several related fungi within the Aspergillus genus produce a number of metabolites of opposite absolute configuration including (+)- or (-)-versicolamide B. These alkaloids are hypothesized to arise via biosynthetic Diels-Alder reactions implying that each Aspergillus species possesses enantiomerically distinct Diels-Alderases. Herein, experimental validation of these biosynthetic proposals via deployment of the IMDA reaction as a key step in the asymmetric total syntheses of (+)- and (-)-versicolamide B is described. Laboratory validation of the proposed biosynthetic Diels-Alder construction, coupled with the secondary metabolite profile of the producing fungi, reveals that each Aspergillus species has evolved enantiomerically distinct indole oxidases, as well as enantiomerically distinct Diels-Alderases.