QM/MM study of the taxadiene synthase mechanism

Docking carbocations into terpene synthase active sites using chemically meaningful constraints-The TerDockin approach

Terpenes are a diverse class of natural products which have long been sought after for their chemical properties as medicine, perfumes, and for food flavoring. Computational docking studies of terpene mechanisms have been a challenge due to the lack of strong directing groups which many docking programs rely on. In this chapter, we dive into our computational method Terdockin (Terpene-Docking) as a successful methodology in modeling terpene synthase mechanisms. This method could also be used as inspiration for any multi-ligand docking project.

Catalytic Reaction Mechanism of Bacterial GH92 α-1,2-Mannosidase: A QM/MM Metadynamics Study

The catalytic mechanism of aC a + 2 ${C{a}^{+2}}$ -dependent family 92 α ${{\rm \alpha }}$ -mannosidase, which is abundantly present in human gut flora and malfunctions leading to the lysosomal storage disease α-mannosidosis, has been investigated using quantum mechanics/molecular mechanics and metadynamics methods. Computational efforts show that the enzyme follows a conformational itinerary of and theC a + 2 ${C{a}^{+2}}$ ion serves a dual purpose, as it not only distorts the sugar ring but also plays a crucial role in orchestrating the arrangement of catalytic residues. This orchestration, in turn, contributes to the facilitation of O S 2 ${{{\rm \ }}^{{\rm O}}{{\rm S}}_{2}}$ conformers for the ensuing reaction. This mechanistic insight is well-aligned with the experimental predictions of the catalytic pathway, and the computed energies are of the same order of magnitude as the experimental estimations. Hence, our results extend the mechanistic understanding of glycosidases.

Multiscale QM/MM modelling of catalytic systems with ChemShell

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods are a powerful computational tool for the investigation of all forms of catalysis, as they allow for an accurate description of reactions occurring at catalytic sites in the context of a complicated electrostatic environment. The scriptable computational chemistry environment ChemShell is a leading software package for QM/MM calculations, providing a flexible, high performance framework for modelling both biomolecular and materials catalysis. We present an overview of recent applications of ChemShell to problems in catalysis and review new functionality introduced into the redeveloped Python-based version of ChemShell to support catalytic modelling. These include a fully guided workflow for biomolecular QM/MM modelling, starting from an experimental structure, a periodic QM/MM embedding scheme to support modelling of metallic materials, and a comprehensive set of tutorials for biomolecular and materials modelling.

Improving solubility and copy number of taxadiene synthase to enhance the titer of taxadiene in Yarrowia lipolytica

Taxadiene is an important precursor for the biosynthesis of highly effective anticancer drug paclitaxel, but its microbial biosynthesis yield is very low. In this study, we employed Yarrowia lipolytica as a microbial host to produce taxadiene. First, a "push-pull" strategy was adopted to increase taxadiene production by 234%. Then taxadiene synthase was fused with five solubilizing tags respectively, leading a maximum increase of 62.3% in taxadiene production when fused with SUMO. Subsequently, a multi-copy iterative integration method was used to further increase taxadiene titer, achieving the maximum titer of 23.7 mg/L in shake flask culture after three rounds of integration. Finally, the taxadiene titer was increased to 101.4 mg/L by optimization of the fed-batch fermentation conditions. This is the first report of taxadiene biosynthesis accomplished in Y. lipolytica, serving as a good example for the sustainable production of taxadiene and other terpenoids in this oleaginous yeast.

Molecular Determinants of Carbocation Cyclisation in Bacterial Monoterpene Synthases

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.

A Benchmark Study of Quantum Mechanics and Quantum Mechanics-Molecular Mechanics Methods for Carbocation Chemistry

Carbocations play key roles in classical organic reactions and have also been implicated in several enzyme families. A hallmark of carbocation chemistry is multitudes of competing reaction pathways, and to be able to distinguish between pathways with quantum chemical calculations, it is necessary to approach chemical accuracy for relative energies between carbocations. Here, we present an extensive study of the performance of selected density functional theory (DFT) methods in describing the thermochemistry and kinetics of carbocations and their corresponding neutral alkenes both in the gas-phase and within a hybrid quantum mechanics-molecular mechanics (QM/MM) framework. The density functionals are benchmarked against accurate ab initio methods such as CBS-QB3 and DLPNO-CCSD(T). Based on the findings in the gas-phase calculations of carbocations and alkenes, the best functionals are chosen and tested further for non-covalent interactions in model systems using QM and QM/MM methods. We compute the interaction energies between a model carbocation/alkane and model π, dipole, and hydrophobic systems using DFT and QM(DFT)/MM and compare with DLPNO-CCSD(T). These latter model systems are representative of side chains of amino acids such as phenylalanine/tyrosine, tryptophan, asparagine/glutamine, serine/threonine, methionine, and other hydrophobic groups. The Lennard-Jones parameters of the QM atoms in QM(DFT)/MM calculations are modified to obtain an optimal fit with the QM energies. Finally, a selected carbocation reaction is studied in the gas phase and in implicit chloroform solvent using QM and in explicit chloroform solvent using QM/MM and umbrella sampling simulations. This study highlights the highest accuracy possible with selected density functionals and QM/MM methods but also some limitations in using QM/MM methods for carbocation systems.

Mining methods and typical structural mechanisms of terpene cyclases

Terpenoids, formed by cyclization and/or permutation of isoprenes, are the most diverse and abundant class of natural products with a broad range of significant functions. One family of the critical enzymes involved in terpenoid biosynthesis is terpene cyclases (TCs), also known as terpene synthases (TSs), which are responsible for forming the ring structure as a backbone of functionally diverse terpenoids. With the recent advances in biotechnology, the researches on terpene cyclases have gradually shifted from the genomic mining of novel enzyme resources to the analysis of their structures and mechanisms. In this review, we summarize both the new methods for genomic mining and the structural mechanisms of some typical terpene cyclases, which are helpful for the discovery, engineering and application of more and new TCs.

Chromosome-level genome of Himalayan yew provides insights into the origin and evolution of the paclitaxel biosynthetic pathway

Taxus, commonly known as yew, is a well-known gymnosperm with great ornamental and medicinal value. In this study, by assembling a chromosome-level genome of the Himalayan yew (Taxus wallichiana) with 10.9 Gb in 12 chromosomes, we revealed that tandem duplication acts as the driving force of gene family evolution in the yew genome, resulting in the main genes for paclitaxel biosynthesis, i.e. those encoding the taxadiene synthase, P450s, and transferases, being clustered on the same chromosome. The tandem duplication may also provide genetic resources for the nature to sculpt the core structure of taxoids at different positions and subsequently establish the complex pathway of paclitaxel by neofunctionalization. Furthermore, we confirmed that there are two genes in the cluster encoding isoenzymes of a known enzyme in the paclitaxel biosynthetic pathway. The reference genome of the Himalayan yew will serve as a platform for decoding the complete biosynthetic pathway of paclitaxel and understanding the chemodiversity of taxoids in gymnosperms.

Recent Advances in Metabolic Engineering, Protein Engineering, and Transcriptome-Guided Insights Toward Synthetic Production of Taxol

The diterpenoid paclitaxel (Taxol®) is a blockbuster anticancer agent that was originally isolated from the Pacific yew (Taxus brevifolia) five decades ago. Despite the wealth of information gained over the years on Taxol research, there still remains supply issues to meet increasing clinical demand. Although alternative Taxol production methods have been developed, they still face several drawbacks that cause supply shortages and high production costs. It is highly desired to develop biotechnological production platforms for Taxol, however, there are still gaps in our understanding of the biosynthetic pathway, catalytic enzymes, regulatory and control mechanisms that hamper production of this critical drug by synthetic biology approaches. Over the past 5 years, significant advances were made in metabolic engineering and optimization of the Taxol pathway in different hosts, leading to accumulation of taxane intermediates. Computational and experimental approaches were leveraged to gain mechanistic insights into the catalytic cycle of pathway enzymes and guide rational protein engineering efforts to improve catalytic fitness and substrate/product specificity, especially of the cytochrome P450s (CYP450s). Notable breakthroughs were also realized in engineering the pathway in plant hosts that are more promising in addressing the challenging CYP450 chemistry. Here, we review these recent advances and in addition, we summarize recent transcriptomic data sets of Taxus species and elicited culture cells, and give a bird's-eye view of the information that can be gleaned from these publicly available resources. Recent mining of transcriptome data sets led to discovery of two putative pathway enzymes, provided many lead candidates for the missing steps and provided new insights on the regulatory mechanisms governing Taxol biosynthesis. All these inferences are relevant to future biotechnological production of Taxol.

Redesigning the Molecular Choreography to Prevent Hydroxylation in Germacradien-11-ol Synthase Catalysis

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.

Enhanced production of taxadiene in Saccharomyces cerevisiae

Background

Cost-effective production of the highly effective anti-cancer drug, paclitaxel (Taxol^®), remains limited despite growing global demands. Low yields of the critical taxadiene precursor remains a key bottleneck in microbial production. In this study, the key challenge of poor taxadiene synthase (TASY) solubility in S. cerevisiae was revealed, and the strains were strategically engineered to relieve this bottleneck.

Results

Multi-copy chromosomal integration of TASY harbouring a selection of fusion solubility tags improved taxadiene titres 22-fold, up to 57 ± 3 mg/L at 30 °C at microscale, compared to expressing a single episomal copy of TASY. The scalability of the process was highlighted through achieving similar titres during scale up to 25 mL and 250 mL in shake flask and bioreactor cultivations, respectively at 20 and 30 °C. Maximum taxadiene titres of 129 ± 15 mg/L and 127 mg/L were achieved through shake flask and bioreactor cultivations, respectively, of the optimal strain at a reduced temperature of 20 °C.

Conclusions

The results of this study highlight the benefit of employing a combination of molecular biology and bioprocess tools during synthetic pathway development, with which TASY activity was successfully improved by 6.5-fold compared to the highest literature titre in S. cerevisiae cell factories.

Enzymatic control of product distribution in terpene synthases: insights from multiscale simulations

In this opinion, we review some recent work on terpene biosynthesis using multiscale simulation approaches, with special focus on contributions from our group. Terpene synthases generate terpenes employing rich carbocation chemistry, including highly specific ring formations, proton, hydride, methyl, and methylene migrations, followed by reaction quenching. In these enzymes, the main catalytic challenge is not rate enhancement, but rather control of intrinsically reactive carbocations and the resulting product distribution. Herein, we review multiscale simulations of selected mono-, sesqui-, and diterpene synthases. We point to the many tools adopted by terpene synthases to achieve correct substrate fold, carbocation formation, carbocation reaction environment, and reaction quenching. A better understanding of the toolbox employed by terpene synthases is expected to aid in the search for new enzymatic and biomimetic synthetic routes to natural and unnatural terpenes.

Exploring the catalytic cascade of cembranoid biosynthesis by combination of genetic engineering and molecular simulations

While chemical steps involved in bioactive cembranoid biosynthesis have been examined, the corresponding enzymatic mechanisms leading to their formation remain elusive. In the tobacco plant, Nicotiana tabacum, a putative cembratriene-ol synthase (CBTS) initiates the catalytic cascade that lead to the biosynthesis of cembratriene-4,6-diols, which displays antibacterial- and anti-proliferative activities. We report here on structural homology models, functional studies, and mechanistic explorations of this enzyme using a combination of biosynthetic and computational methods. This approach guided us to develop an efficient de novo production of five bioactive non- and monohydroxylated cembranoids. Our homology models in combination with quantum and classical simulations suggested putative principles of the CBTS catalytic cycle, and provided a possible rationale for the formation of premature olefinic side products. Moreover, the functional reconstruction of a N. tabacum-derived class II P450 with a cognate CPR, obtained by transcriptome mining provided for production of bioactive cembratriene-4,6-diols. Our combined findings provide mechanistic insights into cembranoid biosynthesis, and a basis for the sustainable industrial production of highly valuable bioactive cembranoids.

The Impact of Macrocycle Conformation on the Taxadiene-Forming Carbocation Cascade: Insight Gained from Sobralene, a Recently Discovered Verticillene Isomer

DFT calculations on the carbocation intermediates that connect the biosynthetic pathways leading to the sand fly pheromone sobralene and taxadiene have been made. Establishment of the conformation of the macrocyclic carbocation intermediate required to produce the (Z)-C8,C9 alkene bond in sobralene has identified new conformations of the verticillyl carbocation intermediates on the taxadiene biosynthetic pathway. These "sobralene-like" carbocation conformations provide an exothermic pathway to taxadiene and are validated by comparison to closely related structures (X-ray and NMR).