Concerted Cyclization of Lanosterol C-Ring and D-Ring Under Human Oxidosqualene Cyclase Catalysis: An ab Initio QM/MM MD Study

Mining coral-derived terpene synthases and mechanistic studies of the coral biflorane synthase

Biflorane diterpenoids are unique natural products often seen in marine animals. Recent studies have reported a small number of biflorane synthases. However, the catalytic mechanism and structural basis for biflorane formation remain unclear. To address these issues, we conducted genome mining of terpene synthases from the sea whip coral Paramuricea clavata, resulting in the discovery of a biflorane synthase PcTS1. We performed a series of isotope labeling, crystallography, quantum mechanics/molecular mechanics calculations, and mutagenesis studies toward PcTS1 to investigate the mechanism. Isotopic labeling studies, together with calculations, elucidate a cascade of 1,10-cyclization, 1,3-hydride shift, 1,6-cyclization, 1,2-hydride shift, 2,6-cyclization, cyclopropane ring opening, and deprotonation by the generated pyrophosphate, forming the biflorane scaffold. Crystallography, quantum mechanics/molecular mechanics, and mutagenesis studies confirmed the cascade and produced different terpene scaffolds. Our work demonstrated the mechanism of marine biflorane formation, elucidated the second crystal structure of a coral terpene synthase, and realized the terpene skeleton expansion.

A theoretical comparison of different third component content in ternary organic solar cells

Ternary solar cells have been rapidly developed in the realm of organic solar cells (OSCs). The incorporation of a third component into a cell results in a complicated active layer morphology, and the relation of this morphology to power conversion efficiency remains elusive. In this work, two ternary active layers, B1:Y7 (10 wt%):BO-4Cl and B1:Y7 (50 wt%):BO-4Cl are constructed, and the reasons for the differences in PCE caused by varying the Y7 content are investigated using theoretical calculations. Firstly, four groups of binary complexes (B1:BO-4Cl-10 wt%, B1:BO-4Cl-50 wt%, B1:Y7-10 wt%, B1:Y7-50 wt%) were examined using molecular dynamics simulation and the stacking patterns of the complexes could mainly be categorized into three groups (IC-T, IC-BDT, IC-RHD). The results showed that with an increase of the Y7 content, the proportion of IC-T stacking decreased while IC-BDT stacking increased. Moreover, the properties of each stacking pattern were calculated and IC-T stacking was found to have a greater charge separation coupling and rate, and a smaller interaction energy. With more IC-T stacking, the number of charge transfer (CT) states and CT mechanisms in B1:BO-4Cl-10 wt% and B1:Y7-10 wt% improves the PCE of B1:Y7 (10 wt%):BO-4Cl. For the trimers, a greater number of CT states and CT pathways can also facilitate efficient charge separation in B1:Y7 (10 wt%):BO-4Cl. Additionally, this work provides basic knowledge of the influence that the third component content on cell performance, providing theoretical instruction for experimental work based on Y-series non-fullerene acceptor materials.

Enzymatic-related network of catalysis, polyamine, and tumors for acetylpolyamine oxidase: from calculation to experiment

The regulation of enzymes and development of polyamine analogs capable of controlling the dynamics of endogenous polyamines to achieve anti-tumor effects is one of the biggest challenges in polyamine research. However, the root of the problem remains unsolved. This study represents a significant milestone as it unveils, for the first time, the comprehensive catalytic map of acetylpolyamine oxidase that includes chemical transformation and product release kinetics, by utilizing multiscale simulations with over six million dynamical snapshots. The transportation of acetylspermine is strongly exothermic, and high binding affinity of enzyme and reactant is observed. The transfer of hydride from polyamine to FAD is the rate-limiting step, via an H-shift coupled electron transfer mechanism. The two products are released in a detour stepwise mechanism, which also impacts the enzymatic efficiency. Inspired by these mechanistic insights into enzymatic catalysis, we propose a novel strategy that regulates the polyamine level and catalytic progress through the action of His64. Directly suppressing APAO by mutating His64 further inhibited growth and migration of tumor cells and tumor tissue in vitro and in vivo. Therefore, the network connecting microcosmic and macroscopic scales opens up new avenues for designing polyamine compounds and conducting anti-tumor research in the future.

Decoding Catalysis by Terpene Synthases

The review by Christianson, published in 2017 on the twentieth anniversary of the emergence of the field, summarizes the foundational discoveries and key advances in terpene synthase/cyclase (TS) biocatalysis (Christianson, D. W. Chem Rev2017, 117 (17), 11570-11648. DOI: 10.1021/acs.chemrev.7b00287). Here, we review the TS literature published since then, bringing the field up to date and looking forward to what could be the near future of TS rational design. Many revealing discoveries have been made in recent years, building on the knowledge and fundamental principles uncovered during those initial two decades of study. We use these to explore TS reaction chemistry and see how a combined experimental and computational approach helps to decipher the complexities of TS catalysis. Revealed are a suite of catalytic motifs which control product outcome in TSs, some obvious, some more subtle. We examine each in detail, using the most recent papers and insights to illustrate how exactly this fascinating class of enzymes takes a single acyclic substrate and turns it into the many thousands of complex terpenoids found in Nature. We then explore some of the recent strategies for TS engineering, including machine learning and other data-driven approaches. From this, rational and predictive engineering of TSs, "designer terpene synthases", will begin to emerge as a realistic goal.

QM/MM Modeling Aided Enzyme Engineering in Natural Products Biosynthesis

Natural products and their derivatives are widely used across various industries, particularly pharmaceuticals. Modern engineered biosynthesis provides an alternative way of producing and meeting the growing need for diverse natural products. Natural enzymes, on the other hand, often exhibit unsatisfactory catalytic characteristics and necessitate further enzyme engineering modifications. QM/MM, as a powerful and extensively used computational tool in the field of enzyme catalysis, has been increasingly applied in rational enzyme engineering over the past decade. In this review, we summarize recent advances in QM/MM computational investigation on enzyme catalysis and enzyme engineering for natural product biosynthesis. The challenges and perspectives for future QM/MM applications aided enzyme engineering in natural product biosynthesis will also be discussed.

Molecular insights into the catalytic promiscuity of a bacterial diterpene synthase

Diterpene synthase VenA is responsible for assembling venezuelaene A with a unique 5-5-6-7 tetracyclic skeleton from geranylgeranyl pyrophosphate. VenA also demonstrates substrate promiscuity by accepting geranyl pyrophosphate and farnesyl pyrophosphate as alternative substrates. Herein, we report the crystal structures of VenA in both apo form and holo form in complex with a trinuclear magnesium cluster and pyrophosphate group. Functional and structural investigations on the atypical 115DSFVSD120 motif of VenA, versus the canonical Asp-rich motif of DDXX(X)D/E, reveal that the absent second Asp of canonical motif is functionally replaced by Ser116 and Gln83, together with bioinformatics analysis identifying a hidden subclass of type I microbial terpene synthases. Further structural analysis, multiscale computational simulations, and structure-directed mutagenesis provide significant mechanistic insights into the substrate selectivity and catalytic promiscuity of VenA. Finally, VenA is semi-rationally engineered into a sesterterpene synthase to recognize the larger substrate geranylfarnesyl pyrophosphate.

Discovering a uniform functional trade-off of the CBC-type 2,3-oxidosqualene cyclases and deciphering its chemical logic

Many functionally promiscuous plant 2,3-oxidosqualene cyclases (OSCs) have been found, but complete functional reshaping is rarely reported. In this study, we have identified two new plant OSCs: a unique protostadienol synthase (AoPDS) and a common cycloartenol synthase (AoCAS) from Alisma orientale (Sam.) Juzep. Multiscale simulations and mutagenesis experiments revealed that threonine-727 is an essential residue responsible for protosta-13 (17),24-dienol biosynthesis in AoPDS and that the F726T mutant completely reshapes the native function of AoCAS into a PDS function to yield almost exclusively protosta-13 (17),24-dienol. Unexpectedly, various native functions were uniformly reshaped into a PDS function by introducing the phenylalanine → threonine substitution at this conserved position in other plant and non-plant chair-boat-chair-type OSCs. Further computational modeling elaborated the trade-off mechanisms of the phenylalanine → threonine substitution that leads to the PDS activity. This study demonstrates a general strategy for functional reshaping by using a plastic residue based on the decipherment of the catalytic mechanism.

Site-directed mutagenesis and substrate compatibility to reveal the structure-function relationships of plant oxidosqualene cyclases

Covering: up to May 2020Oxidosqualene cyclases (OSCs) catalyze one of the most complex polycyclization reactions in nature, using the linear 2,3-oxidosqualene to generate an array of triterpene skeletons in plants. Despite the structural diversity of the products, the protein sequences of plant OSCs are highly conserved, where a few key amino acids could govern the product selectivity. Due to the absence of crystal structures, site-directed mutagenesis and substrate structural modification become key approaches to understand the cyclization mechanism. In this review, 98 mutation sites in 25 plant OSCs have been summarized, and the conserved key residues have been identified by sequence alignment. Structure-function relationships are further discussed. Meanwhile, the substrate selectivity has been summarized to probe the active site cavity of plant OSCs. A total of 77 references are included.

A simpler method affords evaluation of π stabilization by phenylalanine of several biochemical carbocations

Carbocations are important intermediates in the biosynthesis of terpenes and steroids, and it is challenging to try to understand how these relatively unstable species survive even transiently during biochemical reactions. Carbocation-π interaction with aromatic amino acid residues is an important factor in helping to stabilize these positively charged species. However, the short lifetimes of these active site carbocations makes experimental evaluation of the stabilization afforded by such interaction impossible. Computational studies, however, have provided some insight into this phenomenon. Herein we report a simple, computationally efficient method to estimate such stabilization energies afforded by phenylalanine to biochemical carbocation intermediates. A model is constructed in which the biochemical carbocation is replaced by an appropriate carbocation mimic (t-butyl or dimethylallyl). This substitute carbocation is then aligned with an ethylbenzene serving as a surrogate for each proximate phenylalanine in a geometry that replicates as closely as possible the orientation of that phenylalanine using measurements made on an X-ray structure of an enzyme active site in which a carbocation surrogate is bound. Density functional theory computations on such models were then used to yield estimates of stabilization energies. Application of this method to the tertiary carbocation formed in the reaction catalyzed by geranyl diphosphate C-methyl transferase gave a stabilization energy (-12.3 kcal mol-1) that was essentially identical to that obtained previously by analysis of a much more computationally demanding model of the active site. As a check on the accuracy of the simpler method, it was applied with similar success to the farnesyl cation formed in the reaction catalyzed by aristolochene synthase that is stabilized by cation-π interaction with two phenylalanines. Application of this method is also described to estimate carbocation-π stabilization, by the same two phenylalanines, of the final carbocation intermediate leading to aristolochene through analysis of the X-ray structure of an inhibitor of that carbocation bound in the active site of aristolochene synthase. Finally, the stabilization, by either of two phenylalanines, of six carbocation intermediates in the oxidosqualene cyclase-catalyzed formation of lanosterol is estimated by comparable analysis of an X-ray structure of that reaction product bound in the enzyme active site.

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.

Theoretical Study of Sesterfisherol Biosynthesis: Computational Prediction of Key Amino Acid Residue in Terpene Synthase

The cyclization mechanisms involved in the biosynthesis of sesterterpenes are not fully understood. For example, there are two plausible reaction pathways for sesterfisherol biosynthesis, which differ in the order of ring cyclization: A-D-B/C (Path a) and A-B-C/D (Path b). It is difficult to capture intermediates of terpene cyclization, which is a complex, domino-type reaction, and so here we employed a combination of experimental and computational methods. Density functional theory calculations revealed unexpected intermediates and transition states, and implied that C-H···π interaction between a carbocation intermediate and an aromatic residue of sesterfisherol synthase (NfSS) plays a critical role, serving to accelerate the 1,2-H shift (thereby preventing triquinane carbocation formation) and to protect reactive carbocation intermediates from bases such as pyrophosphate or water in the active site. Site-directed mutagenesis of NfSS guided by docking simulations confirmed that phenylalanine F191 is a critical amino acid residue for sesterfisherol synthase, as the F191A mutant of NfSS produces novel sesterterpenes, but not sesterfisherol. Although both pathways are energetically viable, on the basis of our computational and experimental results, NfSS-mediated sesterfisherol biosynthesis appears to proceed via Path a. These findings may also provide new insight into the cyclization mechanisms in related sesterterpene synthases.

Catalytic promiscuity of the non-native FPP substrate in the TEAS enzyme: non-negligible flexibility of the carbocation intermediate

By QM(DFT)/MM MD simulations, it has been revealed that the non-native substrate catalytic promiscuity of TEAS (one of the sesquiterpene cyclases) is mostly attributable to its notable conformational flexibility of the branching intermediate bisabolyl cation.

Importance of Inherent Substrate Reactivity in Enzyme-Promoted Carbocation Cyclization/Rearrangements

The importance of inherent substrate reactivity for terpene synthase enzymes is discussed, with a focus on recent experimental tests of predictions derived from computations on gas-phase reactivity of carbocations.

Computational studies on the cyclization of squalene to the steroids and hopenes

A review of computational studies of the related biosyntheses of steroids and hopenes reported during the last two decades is presented.

Mechanistically informed predictions of binding modes for carbocation intermediates of a sesquiterpene synthase reaction

Sesquiterpenoids comprise a class of terpenoid natural products with thousands of compounds that are highly diverse in structure, generally containing a polycyclic carbon backbone that is constructed by a sesquiterpene synthase. Decades of experimental and computational studies have demonstrated that these enzymes generate a carbocation in the active site, which undergoes a series of structural rearrangements until a product is formed via deprotonation or nucleophile attack. However, for the vast majority of these enzymes the productive binding orientation of the intermediate carbocations has remained unclear. In this work, a method that combines quantum mechanics and computational docking is used to generate an all-atom model of every putative intermediate formed in the context of the enzyme active site for tobacco epi-aristolochene synthase (TEAS). This method identifies a single pathway that links the first intermediate to the last, enabling us to propose the first high-resolution model for the reaction intermediates in the active site of TEAS, and providing testable predictions.

Crucial Role of the Double Bond Isomerism in the Steroid B-Ring on the Membrane Properties of Sterols. Grazing Incidence X-Ray Diffraction and Brewster Angle Microscopy Studies

Three cholesterol precursors-desmosterol, zymosterol, and lanosterol-were comprehensively characterized in monolayers formed at the air/water interface. The studies were based on registration of the surface pressure (π)-area (A) isotherms complemented with in situ analysis performed with application of modern physicochemical techniques: grazing incidence X-ray diffraction (GIXD) and Brewster angle microscopy (BAM). In this approach we were interested in the correlation between molecular structures of the studied sterols found in the cholesterol biosynthetic pathway and their membrane properties. Our results revealed that only desmosterol behaves in Langmuir monolayers comparably to cholesterol, the molecules of which arrange in the monolayers into a hexagonal lattice, while the two remaining sterols possess extremely different properties. We found that molecules of both zymosterol and lanosterol are organized on the water surface in the two-dimensional oblique unit cells despite the fact that they are oriented perpendicular to the monolayer plane. The comparison of chemical structures of the investigated sterols leads to the conclusion that the only structural motive that can be responsible for such unusual behavior is the double bond in the B sterol ring, which is located in desmosterol in a different position from in the other two sterols. This issue, which was neglected in the scientific literature, seems to have crucial importance for sterol activity in biomembranes. We showed that this structural modification in sterol molecules is directly responsible for their adaptation to proper functioning in biomembranes.

A full picture of enzymatic catalysis by hydroxynitrile lyases from Hevea brasiliensis: protonation dependent reaction steps and residue-gated movement of the substrate and the product

Hydroxynitrile lyases (HNLs) defend plants from herbivores and microbial attack by releasing cyanide from hydroxynitriles.