Molecular dynamics study of taxadiene synthase catalysis

Insights into taxadiene synthase catalysis and promiscuity facilitated by mutability landscape and molecular dynamics

MAIN CONCLUSION

Protein modeling, carbocation docking, and molecular dynamics along with structure-based mutability landscapes provided insight into taxadiene synthase catalysis (first step of the anticancer Taxol biosynthesis), protein structure-function correlations, and promiscuity. Plant terpenes belong to one of the largest and most diverse classes of natural products. This diversity is driven by the terpene synthase enzyme family which comprises numerous different synthases, several of which are promiscuous. Taxadiene synthase (TXS) is a class I diterpene synthase that catalyzes the first step in the biosynthesis pathway of the diterpene Taxol, an anticancer natural product produced by the Taxus plant. Exploring the molecular basis of TXS catalysis and its promiscuous potential garnered interest as a necessary means for understanding enzyme evolution and engineering possibilities to improve Taxol biosynthesis. A catalytically active closed conformation TXS model was designed using the artificial intelligence system, AlphaFold, accompanied by docking and molecular dynamics simulations. In addition, a mutability landscape of TXS including 14 residues was created to probe for structure-function relations. The mutability landscape revealed no mutants with improved catalytic activity compared to wild-type TXS. However, mutations of residues V584, Q609, V610, and Y688 showed high degree of promiscuity producing cembranoid-type and/or verticillene-type major products instead of taxanes. Mechanistic insights into V610F, V584M, Q609A, and Y688C mutants compared to the wild type revealed the trigger(s) for product profile change. Several mutants spanning residues V584, Q609, Y688, Y762, Q770, and F834 increased production of taxa-4(20),11(12)-diene which is a more favorable substrate for Taxol production compared to taxa-4(5),11(12)-diene. Finally, molecular dynamics simulations of the TXS reaction cascade revealed residues involved in ionization, carbocation stabilization, and cyclization ushering deeper understanding of the enzyme catalysis mechanism.

Unveiling the Methyl Transfer Mechanisms in the Epigenetic Machinery DNMT3A-3L: A Comprehensive Study Integrating Assembly Dynamics with Catalytic Reactions

In epigenetic mechanisms, DNA methyltransferase 3 alpha (DNMT3A) acts as an initiator for DNA methylation and prevents the downstream genes from expressing. Perturbations of DNMT3A functions may cause uncontrolled gene expression, resulting in pathogenic consequences such as cancers. It is, therefore, vitally important to understand the catalytic process of DNMT3A in its biological macromolecule assembly, viz., heterotetramer: (DNMT3A-3 L)^dimer. In this study, we utilized molecular dynamics (MD) simulations, Markov State Models (MSM), and quantum mechanics/molecular mechanics simulations (QM/MM) to investigate the de novo methyl transfer process. We identified the dynamics of the key residues relevant to the insertion of the target cytosine (dC) into the catalytic domain of DNMT3A, and the detailed potential energy surface of the seven-step reaction referring to methyl transfer. Our calculated potential energy barrier (22.51 kcal/mol) approximates the former experimental data (23.12 kcal/mol). The conformational change of the 5-methyl-cytosine (5mC) intermediate was found necessary in forming a four-water chain for the elimination step, which is unique to the other DNMTs. The biological assembly facilitates the creation of such a water chain, and the elimination occurs in an asynchronized mechanism in the two catalytic pockets. We anticipate the findings can enable a better understanding of the general mechanisms of the de novo methyl transfer for fulfilling the key enzymatic functions in epigenetics. And the unique elimination of DNMT3A might ignite novel methods for designing anti-cancer and tumor inhibitors of DNMTs.

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.

DFT, molecular docking and molecular dynamics simulation studies on some newly introduced natural products for their potential use against SARS-CoV-2

Throughout the history, natural products always give new paths to develop new drugs. As with many other diseases, natural compounds can be helpful in the treatment of COVID-19. SARS-CoV-2 main protease enzyme has an important role in viral replication and transcription. Therefore, inhibiting this enzyme may be helpful in the treatment of COVID-19. In this study, it is aimed to investigate eight natural compounds which have recently entered the literature, computationally for their potential use against SARS-CoV-2. For this purpose, first, density functional theory (DFT) calculations were performed on the investigated compounds, and energy minimizations, geometry optimizations, vibrational analyses, molecular electrostatic potential map calculations were carried out. After DFT calculations, geometry optimized structures were subjected to molecular docking calculations with the use of SARS-CoV-2 main protease (pdb id: 5r80) and top-scoring ligand-receptor complexes were obtained. In the next part of the study, molecular dynamics (MD) simulations were performed on the top-scoring ligand-receptor complexes to investigate the stability of the ligand-receptor complexes and the interactions between ligands and receptor in more detail. Additionally, in this part of the study, binding free energies are calculated with the use of molecular mechanics with Poisson-Boltzmann surface area (MM-PBSA) method. Results showed that, all ligand-receptor complexes remain stable during the MD simulations and most of the investigated compounds but especially two of them showed considerably high binding affinity to SARS-CoV-2 main protease. Finally, in the study, ADME (adsorption, desorption, metabolism, excretion) predictions and drug-likeness analyses were performed on the investigated compounds.

Omic analysis of the endangered Taxaceae species Pseudotaxus chienii revealed the differences in taxol biosynthesis pathway between Pseudotaxus and Taxus yunnanensis trees

BACKGROUND

Taxol is an efficient anticancer drug accumulated in Taxus species. Pseudotaxus chienii is an important member of Taxaceae, however, the level of six taxoids in P. chienii is largely unknown.

RESULTS

High accumulation of 10-DAB, taxol, and 7-E-PTX suggested that P. chienii is a good taxol-yielding species for large-scale cultivation. By the omics approaches, a total of 3,387 metabolites and 61,146 unigenes were detected and annotated. Compared with a representative Taxus tree (Taxus yunnanensis), most of the differentially accumulated metabolites and differential expressed genes were assigned into 10 primary and secondary metabolism pathways. Comparative analyses revealed the variations in the precursors and intermediate products of taxol biosynthesis between P. chienii and T. yunnanensis. Taxusin-like metabolites highly accumulated in P. chienii, suggesting a wider value of P. chienii in pharmaceutical industry.

CONCLUSIONS

In our study, the occurrence of taxoids in P. chienii was determined. The differential expression of key genes involved in the taxol biosynthesis pathway is the major cause of the differential accumulation of taxoids. Moreover, identification of a number of differentially expressed transcription factors provided more candidate regulators of taxol biosynthesis. Our study may help to reveal the differences between Pseudotaxus and Taxus trees, and promote resource utilization of the endangered and rarely studied P. chienii.

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.

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

The verticillenes. Pivotal intermediates in the biosynthesis of the taxanes and the phomactins

Recent enzymatic studies, quantum chemical calculations and biomimetic conversions consolidate the role of verticillenes in the biosynthesis of taxanes and phomactins.

Insights Into the Bifunctional Aphidicolan-16-ß-ol Synthase Through Rapid Biomolecular Modeling Approaches

Diterpene synthases catalyze complex, multi-step C-C coupling reactions thereby converting the universal, aliphatic precursor geranylgeranyl diphosphate into diverse olefinic macrocylces that form the basis for the structural diversity of the diterpene natural product family. Since catalytically relevant crystal structures of diterpene synthases are scarce, homology based biomolecular modeling techniques offer an alternative route to study the enzyme's reaction mechanism. However, precise identification of catalytically relevant amino acids is challenging since these models require careful preparation and refinement techniques prior to substrate docking studies. Targeted amino acid substitutions in this protein class can initiate premature quenching of the carbocation centered reaction cascade. The structural characterization of those alternative cyclization products allows for elucidation of the cyclization reaction cascade and provides a new source for complex macrocyclic synthons. In this study, new insights into structure and function of the fungal, bifunctional Aphidicolan-16-ß-ol synthase were achieved using a simplified biomolecular modeling strategy. The applied refinement methodologies could rapidly generate a reliable protein-ligand complex, which provides for an accurate in silico identification of catalytically relevant amino acids. Guided by our modeling data, ACS mutations lead to the identification of the catalytically relevant ACS amino acid network I626, T657, Y658, A786, F789, and Y923. Moreover, the ACS amino acid substitutions Y658L and D661A resulted in a premature termination of the cyclization reaction cascade en-route from syn-copalyl diphosphate to Aphidicolan-16-ß-ol. Both ACS mutants generated the diterpene macrocycle syn-copalol and a minor, non-hydroxylated labdane related diterpene, respectively. Our biomolecular modeling and mutational studies suggest that the ACS substrate cyclization occurs in a spatially restricted location of the enzyme's active site and that the geranylgeranyl diphosphate derived pyrophosphate moiety remains in the ACS active site thereby directing the cyclization process. Our cumulative data confirm that amino acids constituting the G-loop of diterpene synthases are involved in the open to the closed, catalytically active enzyme conformation. This study demonstrates that a simple and rapid biomolecular modeling procedure can predict catalytically relevant amino acids. The approach reduces computational and experimental screening efforts for diterpene synthase structure-function analyses.