Huiying Chu

Modular Pathway Engineering of Diterpenoid Synthases and the Mevalonic Acid Pathway for Miltiradiene Production

Microbial production can be advantageous over the extraction of phytoterpenoids from natural plant sources, but it remains challenging to rationally and rapidly access efficient pathway variants. Previous engineering attempts mainly focused on the mevalonic acid (MVA) or methyl-d-erythritol phosphate (MEP) pathways responsible for the generation of precursors for terpenoids biosynthesis, and potential interactions between diterpenoids synthases were unexplored. Miltiradiene, the product of the stepwise conversion of (E,E,E)-geranylgeranyl diphosphate (GGPP) catalyzed by diterpene synthases SmCPS and SmKSL, has recently been identified as the precursor to tanshionones, a group of abietane-type norditerpenoids rich in the Chinese medicinal herb Salvia miltiorrhiza . Here, we present the modular pathway engineering (MOPE) strategy and its application for rapid assembling synthetic miltiradiene pathways in the yeast Saccharomyces cerevisiae . We predicted and analyzed the molecular interactions between SmCPS and SmKSL, and engineered their active sites into close proximity for enhanced metabolic flux channeling to miltiradiene biosynthesis by constructing protein fusions. We show that the fusion of SmCPS and SmKSL, as well as the fusion of BTS1 (GGPP synthase) and ERG20 (farnesyl diphosphate synthase), led to significantly improved miltiradiene production and reduced byproduct accumulation. The MOPE strategy facilitated a comprehensive evaluation of pathway variants involving multiple genes, and, as a result, our best pathway with the diploid strain YJ2X reached miltiradiene titer of 365 mg/L in a 15-L bioreactor culture. These results suggest that terpenoids synthases and the precursor supplying enzymes should be engineered systematically to enable an efficient microbial production of phytoterpenoids.

Insights into the molecular recognition of the granuphilin C2A domain with PI(4,5)P2.

Phosphatidylinositol 4,5-bisphosphate (PI(4,5)P2) is a key player in regulating the process of excytosis, including insulin secretion. Granuphilin, a tandem C2 domain containing protein, mediates the docking of insulin granules onto plasma membrane. The C2A domain plays key roles in this process through interaction with PI(4,5)P2. In this study, we have investigated the molecular recognition mechanism of granuphilin-C2A domain to PI(4,5)P2 head group, and further to PI(4,5)P2-nanodisc by NMR, ITC, MST and SEC methods. Our results demonstrate that PI(4,5)P2 binds to the concave surface of granuphilin-C2A domain. The key residues involved in the binding were validated by mutation analysis.

Advancement of Polarizable Force Field and Its Use for Molecular Modeling and Design

The most important requirement of biomolecular modeling is to deal with electrostatic energies. The electrostatic polarizability is an important part of electrostatic interaction for simulation systems. However, AMBER, CHARMM, OPLS, GROMOS, MMFF force fields etc. used in the past mostly apply fixed atomic center point charge to describe electrostatic energies, and are not sufficient for considering the influence of the electrostatic polarization. The emergence of polarizable force fields has solved this problem. In recent years, quickly developed polarizable force fields have involved a lot of fields. The chapter relating to polarizable force fields spread over several aspects. Firstly, we reviewed the history of the classical force fields and compared with polarizable force fields to elucidate the advancements of polarizable force fields. Secondly, it is introduced that the application of polarizable force fields to small molecules and biological macromolecules simulation, including molecular design. Finally, a brief development trend and perspective is given on rapidly growing polarizable force fields.

Free energy simulations with the AMOEBA polarizable force field and metadynamics on GPU platform

The free energy calculation library PLUMED has been incorporated into the OpenMM simulation toolkit, with the purpose to perform enhanced sampling MD simulations using the AMOEBA polarizable force field on GPU platform. Two examples, (I) the free energy profile of water pair separation (II) alanine dipeptide dihedral angle free energy surface in explicit solvent, are provided here to demonstrate the accuracy and efficiency of our implementation. The converged free energy profiles could be obtained within an affordable MD simulation time when the AMOEBA polarizable force field is employed. Moreover, the free energy surfaces estimated using the AMOEBA polarizable force field are in agreement with those calculated from experimental data and ab initio methods. Hence, the implementation in this work is reliable and would be utilized to study more complicated biological phenomena in both an accurate and efficient way.

Accurate Evaluation of Ion Conductivity of the Gramicidin A Channel Using a Polarizable Force Field without Any Corrections

Classical molecular dynamic (MD) simulation of membrane proteins faces significant challenges in accurately reproducing and predicting experimental observables such as ion conductance and permeability due to its incapability of precisely describing the electronic interactions in heterogeneous systems. In this work, the free energy profiles of K(+) and Na(+) permeating through the gramicidin A channel are characterized by using the AMOEBA polarizable force field with a total sampling time of 1 μs. Our results indicated that by explicitly introducing the multipole terms and polarization into the electrostatic potentials, the permeation free energy barrier of K(+) through the gA channel is considerably reduced compared to the overestimated results obtained from the fixed-charge model. Moreover, the estimated maximum conductance, without any corrections, for both K(+) and Na(+) passing through the gA channel are much closer to the experimental results than any classical MD simulations, demonstrating the power of AMOEBA in investigating the membrane proteins.

An anisotropic coarse‐grained model based on Gay–Berne and electric multipole potentials and its application to simulate a DMPC bilayer in an implicit solvent model

In this work, we aim at optimizing the performance of the anisotropic GBEMP model, which adopts a framework by combining a Gay–Berne (GB) anisotropic potential with an electric multipole (EMP) potential, in simulating a DMPC lipid bilayer in an implicit solvent model. First, the Gay–Berne parameters were initially obtained by fitting to atomistic profiles of van der Waals interactions between homodimers of molecular fragments while EMP parameters was directly derived from the expansion of point multipoles at predefined EMP sites. Second, the GB and EMP parameters for DMPC molecule were carefully optimized to be comparable to AMBER atomistic model in the calculations of the dipole moments of DMPC monomers adopting different conformations as well as the nonbonded interactions between two DMPC molecules adopting different conformations and separated at various distances. Finally, the GB parameters for DMPC were slightly adjusted in simulating a 72 DMPC bilayer system so that our GBEMP model would be able to reproduce a few important structural properties, namely, thickness ( DHH ), area per lipid ( AL ) and volume per lipid ( VL ). Meanwhile, the atomistic and experimental results for electron density profiles and order parameters were reproduced reasonably well by the GBEMP model, demonstrating the promising feature of GBEMP model in modeling lipid systems. Finally, we have shown that current GBEMP model is more efficient by a factor of about 25 than AMBER atomistic point charge model.

Hypotaurine evokes a malignant phenotype in glioma through aberrant hypoxic signaling.

Metabolomics has shown significant potential in identifying small molecules specific to tumor phenotypes. In this study we analyzed resected tissue metabolites using capillary electrophoresis-mass spectrometry and found that tissue hypotaurine levels strongly and positively correlated with glioma grade. In vitro studies were conducted to show that hypotaurine activates hypoxia signaling through the competitive inhibition of prolyl hydroxylase domain-2. This leads to the activation of hypoxia signaling as well as to the enhancement of glioma cell proliferation and invasion. In contrast, taurine, the oxidation metabolite of hypotaurine, decreased intracellular hypotaurine and resulted in glioma cell growth arrest. Lastly, a glioblastoma xenograft mice model was supplemented with taurine feed and exhibited impaired tumor growth. Taken together, these findings suggest that hypotaurine is an aberrantly produced oncometabolite, mediating tumor molecular pathophysiology and progression. The hypotaurine metabolic pathway may provide a potentially new target for glioblastoma diagnosis and therapy.

Specific dephosphorylation of Janus Kinase 2 by protein tyrosine phosphatases

Many protein kinases are activated through phosphorylation of an activation loop thereby turning on downstream signaling pathways. Activation of JAK2, a nonreceptor tyrosine kinase with an important role in growth factor and cytokine signaling, requires phosphorylation of the 1007 and 1008 tyrosyl residues. Dephosphorylation of these two sites by phosphatases presumably inactivates the enzyme, but the underlying mechanism is not known. In this study, we employed MALDI‐TOF/TOF and triple quadrupole mass spectrometers to analyze qualitatively and quantitatively the dephosphorylation process by using synthetic peptides derived from the tandem autophosphorylation sites (Y1007 and Y1008) of human JAK2. We found that tyrosine phosphatases catalyzed the dephosphorylation reaction sequentially, but different enzymes exhibited different selectivity. Protein tyrosine phosphatase 1B caused rapid dephosphorylation of Y1008 followed by Y1007, while SHP1 and SHP2 selectively dephosphorylated Y1008 only, and yet HePTP randomly removed a single phosphate from either Y1007 or Y1008, leaving behind mono‐phosphorylated peptides. The specificity of dephosphorylation was further confirmed by molecular modeling. The data reveal multiple modes of JAK2 regulation by tyrosine phosphatases, reflecting a complex, and intricate interplay between protein phosphorylation and dephosphorylation.

Binding Modes and Interaction Mechanism Between Different Base Pairs and Methylene Blue Trihydrate: A Quantum Mechanics Study

Different quantum mechanic methods have been evaluated for the calculation of binding modes and interactions between intercalators with different DNA base pairs by comparing with the results of MP2, which is very expensive, indicating that WB97XD method under 6-311+G* basis set is able to efficiently reproduce MP2 results. We discovered that the methylene blue trihydrate intercalated into the DNA base pairs, and DNA intercalation increased the distance between DNA base pairs, depending on the types of DNA bases. According to the binding energy results, it was found that the intercalation of methylene blue trihydrate into AA-TT base pair was more favorable in the orientation of nitrogen than other directions and intercalation, and the electric charge was transferred from methylene blue trihydrate to the AA-TT base pair. The analysis of change in the charge density shows that changes often take place in the heavy atom in the middle of the system which the charge density changes most remarkable.

Mechanistic insight into the functional transition of the enzyme guanylate kinase induced by a single mutation

Dramatic functional changes of enzyme usually require scores of alterations in amino acid sequence. However, in the case of guanylate kinase (GK), the functional novelty is induced by a single (S→P) mutation, leading to the functional transition of the enzyme from a phosphoryl transfer kinase into a phosphorprotein interaction domain. Here, by using molecular dynamic (MD) and metadynamics simulations, we provide a comprehensive description of the conformational transitions of the enzyme after mutating serine to proline. Our results suggest that the serine plays a crucial role in maintaining the closed conformation of wild-type GK and the GMP recognition. On the contrary, the S→P mutant exhibits a stable open conformation and loses the ability of ligand binding, which explains its functional transition from the GK enzyme to the GK domain. Furthermore, the free energy profiles (FEPs) obtained by metadymanics clearly demonstrate that the open-closed conformational transition in WT GK is positive correlated with the process of GMP binding, indicating the GMP-induced closing motion of GK enzyme, which is not observed in the mutant. In addition, the FEPs show that the S→P mutation can also leads to the mis-recognition of GMP, explaining the vanishing of catalytic activity of the mutant.