Unveiling the Catalytic Mechanism of NADP+-Dependent Isocitrate Dehydrogenase with QM/MM Calculations

Gaining Insight into the Catalytic Mechanism of the R132H IDH1 Mutant: A Synergistic DFT Cluster and Experimental Investigation

Human isocitrate dehydrogenase 1 (IDH1) is an enzyme that is found in humans that plays a critical role in aerobic metabolism. As a part of the citric acid cycle, IDH1 becomes responsible for catalyzing the oxidative decarboxylation of isocitrate to form α-ketoglutarate (αKG), with nicotinamide adenine dinucleotide phosphate (NADP+) as a cofactor. Strikingly, mutations of the IDH1 enzyme have been discovered in several cancers including glioblastoma multiforme (GBM), a highly aggressive form of brain cancer. It has been experimentally determined that single-residue IDH1 mutations occur at a very high frequency in GBM. Specifically, the IDH1 R132H mutation is known to produce (D)2-hydroxyglutarate (2HG), a recognized oncometabolite. Using the previously determined catalytic mechanism of IDH1, a DFT QM model was developed to study the mechanistic properties of IDH1 R132H compared to wild type enzyme. Validating these insights, biochemical in vitro assays of metabolites produced by mutant vs wild type enzymes were measured and compared. From the results discussed herein, we discuss the mechanistic impact of mutations in IDH1 on its ability to catalyze the formation of αKG and 2HG.

DszA Catalyzes C-S Bond Cleavage through N5-Hydroperoxyl Formation

Due to its detrimental impact on human health and the environment, regulations demand ultralow sulfur levels on fossil fuels, in particular in diesel. However, current desulfurization techniques are expensive and cannot efficiently remove heteroaromatic sulfur compounds, which are abundant in crude oil and concentrate in the diesel fraction after distillation. Biodesulfurization via the four enzymes of the metabolic 4S pathway of the bacterium Rhodococcus erythropolis (DszA-D) is a possible solution. However, the 4S pathway needs to operate at least 500 times faster for industrial applicability, a goal currently pursued through enzyme engineering. In this work, we unveil the catalytic mechanism of the flavin monooxygenase DszA. Surprisingly, we found that this enzyme follows a recently proposed atypical mechanism that passes through the formation of an N5OOH intermediate at the re side of the cofactor, aided by a well-defined, predominantly hydrophobic O2 pocket. Besides clarifying the unusual chemical mechanism of the complex DszA enzyme, with obvious implications for understanding the puzzling chemistry of flavin-mediated catalysis, the result is crucial for the rational engineering of DszA, contributing to making biodesulfurization attractive for the oil refining industry.

Long-term conservation tillage increase cotton rhizosphere sequestration of soil organic carbon by changing specific microbial CO2 fixation pathways in coastal saline soil

Coastal saline soil is an important reserve resource for arable land globally. Data from 10 years of continuous stubble return and subsoiling experiments have revealed that these two conservation tillage measures significantly improve cotton rhizosphere soil organic carbon sequestration in coastal saline soil. However, the contribution of microbial fixation of atmospheric carbon dioxide (CO2) has remained unclear. Here, metagenomics and metabolomics analyses were used to deeply explore the microbial CO2 fixation process in rhizosphere soil of coastal saline cotton fields under long-term stubble return and subsoiling. Metagenomics analysis showed that stubble return and subsoiling mainly optimized CO2 fixing microorganism (CFM) communities by increasing the abundance of Acidobacteria, Gemmatimonadetes, and Chloroflexi, and improving composition diversity. Conjoint metagenomics and metabolomics analyses investigated the effects of stubble return and subsoiling on the reverse tricarboxylic acid (rTCA) cycle. The conversion of citrate to oxaloacetate was inhibited in the citrate cleavage reaction of the rTCA cycle. More citrate was converted to acetyl-CoA, which enhanced the subsequent CO2 fixation process of acetyl-CoA conversion to pyruvate. In the rTCA cycle reductive carboxylation reaction from 2-oxoglutarate to isocitrate, synthesis of the oxalosuccinate intermediate product was inhibited, with strengthened CO2 fixation involving the direct conversion of 2-oxoglutarate to isocitrate. The collective results demonstrate that stubble return and subsoiling optimizes rhizosphere CFM communities by increasing microbial diversity, in turn increasing CO2 fixation by enhancing the utilization of rTCA and 3-hydroxypropionate/4-hydroxybutyrate cycles by CFMs. These events increase the microbial CO2 fixation in the cotton rhizosphere, thereby promoting the accumulation of microbial biomass, and ultimately improving rhizosphere soil organic carbon. This study clarifies the impact of conservation tillage measures on microbial CO2 fixation in cotton rhizosphere of coastal saline soil, and provides fundamental data for the improvement of carbon sequestration in saline soil in agricultural ecosystems.

Scutellarin activates IDH1 to exert antitumor effects in hepatocellular carcinoma progression

Isochlorate dehydrogenase 1 (IDH1) is an important metabolic enzyme for the production of α-ketoglutarate (α-KG), which has antitumor effects and is considered to have potential antitumor effects. The activation of IDH1 as a pathway for the development of anticancer drugs has not been attempted. We demonstrated that IDH1 can limit glycolysis in hepatocellular carcinoma (HCC) cells to activate the tumor immune microenvironment. In addition, through proteomic microarray analysis, we identified a natural small molecule, scutellarin (Scu), which activates IDH1 and inhibits the growth of HCC cells. By selectively modifying Cys297, Scu promotes IDH1 active dimer formation and increases α-KG production, leading to ubiquitination and degradation of HIF1a. The loss of HIF1a further leads to the inhibition of glycolysis in HCC cells. The activation of IDH1 by Scu can significantly increase the level of α-KG in tumor tissue, downregulate the HIF1a signaling pathway, and activate the tumor immune microenvironment in vivo. This study demonstrated the inhibitory effect of IDH1-α-KG-HIF1a on the growth of HCC cells and evaluated the inhibitory effect of Scu, the first IDH1 small molecule agonist, which provides a reference for cancer immunotherapy involving activated IDH1.

Origin of the enantioselectivity of alcohol dehydrogenase

Alcohol dehydrogenases (ADH) are a family of enzymes that catalyse the interconversion between ketones/aldehydes and alcohols in the presence of NADPH cofactor. It is challenging to desymmetrise the substituted cyclopentane-1,3-dione by engineering an ADH, while the reaction mechanism of the metal independent ADH remains elusive. Here we measured the conversion of a model substrate 2-benzyl-2-methylcyclopentane-1,3-dione by LbADH and found it predominately gave the (2R,3R) product. Binding mode analysis of the substrate in LbADH from molecular dynamics simulations disclosed the origin of the enantioselectivity of the enzyme; the opening and closing of the loop 191-205 above the substrate are responsible for shaping the binding pocket to orientate the substrate, so as to give different stereoisomer products. Using QM/MM calculations, we elucidated the reaction mechanism of LbADH. Furthermore, we demonstrated the reaction profile corresponding to the production of different stereoisomers, which is in accordance with our experimental observations. This research here will shed a light on the rational engineering of ADH to achieve stereodivergent stereoisomer products.

Drug Design in the Exascale Era: A Perspective from Massively Parallel QM/MM Simulations

The initial phases of drug discovery - in silico drug design - could benefit from first principle Quantum Mechanics/Molecular Mechanics (QM/MM) molecular dynamics (MD) simulations in explicit solvent, yet many applications are currently limited by the short time scales that this approach can cover. Developing scalable first principle QM/MM MD interfaces fully exploiting current exascale machines - so far an unmet and crucial goal - will help overcome this problem, opening the way to the study of the thermodynamics and kinetics of ligand binding to protein with first principle accuracy. Here, taking two relevant case studies involving the interactions of ligands with rather large enzymes, we showcase the use of our recently developed massively scalable Multiscale Modeling in Computational Chemistry (MiMiC) QM/MM framework (currently using DFT to describe the QM region) to investigate reactions and ligand binding in enzymes of pharmacological relevance. We also demonstrate for the first time strong scaling of MiMiC-QM/MM MD simulations with parallel efficiency of ∼70% up to >80,000 cores. Thus, among many others, the MiMiC interface represents a promising candidate toward exascale applications by combining machine learning with statistical mechanics based algorithms tailored for exascale supercomputers.

NAD(P)H binding configurations revealed by time-resolved fluorescence and two-photon absorption

NADH and NADPH play key roles in the regulation of metabolism. Their endogenous fluorescence is sensitive to enzyme binding, allowing changes in cellular metabolic state to be determined using fluorescence lifetime imaging microscopy (FLIM). However, to fully uncover the underlying biochemistry, the relationships between their fluorescence and binding dynamics require greater understanding. Here we accomplish this through time- and polarization-resolved fluorescence and polarized two-photon absorption measurements. Two lifetimes result from binding of both NADH to lactate dehydrogenase and NADPH to isocitrate dehydrogenase. The composite fluorescence anisotropy indicates the shorter (1.3-1.6 ns) decay component to be accompanied by local motion of the nicotinamide ring, pointing to attachment solely via the adenine moiety. For the longer lifetime (3.2-4.4 ns), the nicotinamide conformational freedom is found to be fully restricted. As full and partial nicotinamide binding are recognized steps in dehydrogenase catalysis, our results unify photophysical, structural, and functional aspects of NADH and NADPH binding and clarify the biochemical processes that underlie their contrasting intracellular lifetimes.

Structures of a constitutively active mutant of human IDH3 reveal new insights into the mechanisms of allosteric activation and the catalytic reaction

Human NAD-dependent isocitrate dehydrogenase or IDH3 (HsIDH3) catalyzes the decarboxylation of isocitrate into α-ketoglutarate in the tricarboxylic acid cycle. It consists of three types of subunits (α, β, and γ) and exists and functions as the (αβαγ)₂ heterooctamer. HsIDH3 is regulated allosterically and/or competitively by numerous metabolites including CIT, ADP, ATP, and NADH. Our previous studies have revealed the molecular basis for the activity and regulation of the αβ and αγ heterodimers. However, the molecular mechanism for the allosteric activation of the HsIDH3 holoenzyme remains elusive. In this work, we report the crystal structures of the αβ and αγ heterodimers and the (αβαγ)₂ heterooctamer containing an α-Q139A mutation in the clasp domain, which renders all the heterodimers and the heterooctamer constitutively active in the absence of activators. Our structural analysis shows that the α-Q139A mutation alters the hydrogen-bonding network at the heterodimer-heterodimer interface in a manner similar to that in the activator-bound αγ heterodimer. This alteration not only stabilizes the active sites of both α_Q139Aβ and α_Q139Aγ heterodimers in active conformations but also induces conformational changes of the pseudo-allosteric site of the α_Q139Aβ heterodimer enabling it to bind activators. In addition, the α_Q139A^ICT+Ca+NADβ^NAD structure presents the first pseudo-Michaelis complex of HsIDH3, which allows us to identify the key residues involved in the binding of cofactor, substrate, and metal ion. Our structural and biochemical data together reveal new insights into the molecular mechanisms for allosteric regulation and the catalytic reaction of HsIDH3.

Characteristics of Escherichia coli Isolated from Intestinal Microbiota Children of 0–5 Years Old in the Commune of Abomey-Calavi

Escherichia coli is a commensal bacterium and one of the first bacteria to colonize the digestive tract of newborns after birth. It is characterized by great versatility and metabolic flexibility that allows its survival in different niches. The present study aims at analyzing the diversity of E. coli strains isolated from the intestinal microbiota of children aged from 0 to 5 years in the commune of Abomey-Calavi in Benin. For this purpose, a descriptive and analytical cross-sectional study was conducted. A total of 135 stool samples were collected from the pediatric clinic of Abomey-Calavi. Microbiological analyses were performed according to standard microbiology analytical techniques. The molecular characterization of E. coli was performed by investigating eight genes (dinB, icdA, pabB, polB, putP, trpA, trpB, and uidA) using the PCR technique. The results showed that the average loading rate on stool samples was 3.74 × 107 CFU/g for TAMF. A total of 7 species of bacteria were identified at different proportions: Staphylococcus spp (55.36%), E. coli (14.29%), Klebsiella ornithinolytica (12.5%), Serratia odorifera (5.36%), and Enterobacter aerogenes (5.36%). Interestingly, isolated E. coli presented a resistance of 100% to cefotaxime and aztreonam. In addition, resistances of 95.24% and 50% were observed against erythromycin and nalidixic acid, respectively. The molecular characterization of the isolated E. coli strains allowed us to discover another molecular variation within the isolated strains. Genes encoding the enzymes isocitrate dehydrogenase (icd) and DNA polymerase II (polB) were detected at 96.30% in the isolated E. coli strains. Moreover, the genes encoding the enzymes beta-D-glucuronidase (uidA) and DNA polymerase (dinB) were detected at 88.89% in the isolated E. coli strains. Interestingly, 81.48%, 85.19, 92.59%, and 100% of isolated E. coli strains expressed the genes encoding the enzymes tryptophan synthase subunit A (trpA), proline permease (putP), p-aminobenzoate synthase, and tryptophan synthase subunit B (trpB), respectively. The diversity of E. coli strains reflects the importance of regulatory mechanisms in the adaptation of bacteria to the gut microbiota.

Quantum mechanical study on the activation mechanism of human carbonic anhydrase VII cluster model with bis-histamine schiff bases and bis-spinaceamine derivatives

The activation mechanism of human carbonic anhydrase (hCA) isoform VII, hCA VII, with histamine, histamine bis-Schiff bases and bis-spinaceamine derivatives has been investigated using quantum mechanical calculations. The DFT-D3 method has been employed to calculate in detail the electronic structure and electronic energy of different compounds and complexes throughout the reaction pathway. The model system of hCA VII included the core catalytic center, the Zn²⁺ ion, its three histidine ligands and a hydroxide ion or water molecule coordinated to it. Furthermore, Thr199, Glu106 and the deep water molecule were considered in the model. Five activators of this enzyme, including histamine as standard, in complex with the cluster model of hCA VII were investigated. Thermodynamic functions for the overall reaction and for the complexation between activators and hCA VII were evaluated. Our results demonstrate that the protonatable moiety of these activators participates in proton transfer reactions from the zinc-bound water molecule to the reaction medium, promoting the formation of the catalytically active zinc hydroxide species of the enzyme. The QM analysis revealed that the electrostatic interactions between activators and hCA VII are the driving force of the enzyme-activator complex formation.

THE CATALYTIC MECHANISM OF THE RETAINING GLYCOSYLTRANSFERASE MANNOSYLGLYCERATE SYNTHASE

To protect their intracellular proteins, extremophile microorganisms synthesize molecules called compatible solutes. These molecules are the result of the attachment of a small negatively charged molecule to a sugar molecule. It has been found that these molecules, not only protect the microorganism against osmotic stress, as initially thought, but also against other extreme conditions. The observation that these molecules can confer protection against extreme conditions to isolated enzymes from different organisms made them an exciting prospect for potential biotechnological applications. One of the most widespread compatible solute in hyperthermophile organisms is the molecule 2-O-α-D-mannosyl-D-glycerate (MG). In addition to confer protection to proteins against extreme conditions, MG was found to prevent Alzheimer's β-amyloid aggregation and reduce α-synuclein fibril formation in Parkinson's disease.  In this work we studied, using computational methods, the catalytic mechanism of the synthesis of MG by the enzyme mannosylglycerate synthase (MGS) from the thermophilic bacteria  Rhodothermus marinus . MGS is a promiscuous enzyme, accepting a variety of sugar donors and acceptors. This feature can be used to synthesize other molecules with potential biotechnological applications beyond MG. The unravelling of the catalytic mechanism with atomistic resolution and the associated free energies and electrostatic profiles of the stationary states obtained in the present work will help future investigations to full explore the potential of MGS.

Characterization of cancer-associated IDH2 mutations that differ in tumorigenicity, chemosensitivity and 2-hydroxyglutarate production

The family of isocitrate dehydrogenase (IDH) enzymes is vital for cellular metabolism, as IDH1 and IDH2 are required for the decarboxylation of isocitrate to α-ketoglutarate. Heterozygous somatic mutations in IDH1 or IDH2 genes have been detected in many cancers. They share the neomorphic production of the oncometabolite (R)-2-hydroxyglutarate [(R)-2-HG]. With respect to IDH2, it is unclear whether all IDH2 mutations display the same or differ in tumorigenic properties and degrees of chemosensitivity. Here, we evaluated the three most frequent IDH2 mutations occurring in cancer. The predicted changes to the enzyme structure introduced by these individual mutations are supported by the observed production of (R)-2-HG. However, their tumorigenic properties, response to chemotherapeutic agents, and baseline activation of STAT3 differed. Paradoxically, the varying levels of endogenous (R)-2-HG produced by each IDH2 mutant inversely correlated with their respective growth rates. Interestingly, while we found that (R)-2-HG stimulated the growth of non-transformed cells, (R)-2-HG also displayed antitumor activity by suppressing the growth of tumors harboring wild type IDH2. The mitogenic effect of (R)-2-HG in immortalized cells could be switched to antiproliferative by transformation with oncogenic RAS. Thus, our findings show that despite their shared (R)-2-HG production, IDH2 mutations are not alike and differ in shaping tumor cell behavior and response to chemotherapeutic agents. Our study also reveals that under certain conditions, (R)-2-HG has antitumor properties.

Dynamic Description of the Catalytic Cycle of Malate Enzyme: Stereoselective Recognition of Substrate, Chemical Reaction, and Ligand Release

In protein engineering, investigations of catalytic cycle facilitate rational design of enzymes. In the present work, deeper analysis on the catalytic cycle of malate enzyme (EC 1.1.1.40), an enzyme involved in cancer metabolic and fatty acid synthesis, was performed. In substrate binding, stereoselective recognition of a substrate originates from distance and angle difference between two chiral substrates and Mn²⁺ as well as monodentate or coplanar ion reaction with Arg165. In catalytic transformation, the activation barrier for the hydride transfer of d-malate is 20.28 kcal/mol higher than that for l-malate. The activation barrier for β-decarboxylation of oxaloacetate is about 4.59 kcal/mol higher than the activation barrier for the hydride transfer of l-malate. The effective activation barrier is 16.44 kcal/mol, which is in close agreement with the value derived from the application of transition-state theory and the Eyring equation to k_cat. In ligand release, l/d-malate needs to overcome a higher barrier than pyruvate to break all bonds in parallel and then to escape from the binding pocket. Leu167 and Asn421 comprise a swinging gate to control the product release. The more open gate is possibly required in the direction of pyruvate to l-malate. Our studies are focused on extending structural knowledge regarding the malate enzyme and provided a powerful strategy for future experimental investigations.

A QM/MM study of the reaction mechanism of human β-ketoacyl reductase

Human fatty acid synthase (hFAS) is a multifunctional enzyme involved in a wide diversity of biological functions. For instance, it is a precursor of phospholipids and other complex processes such as the de novo synthesis of long chain fatty acid. Human FAS is also a component of biological membranes and it is implicated in the overexpression of several types of cancers. In this work, we describe the catalytic mechanism of β-ketoreductase (KR), which is a catalytic domain of the hFAS enzyme that catalyzes the reduction of β-ketoacyl to β-hydroxyacyl with the concomitant oxidation of the NADPH cofactor. The catalysis by KR is an intermediate step in the cycle of reactions that elongate the substrate's carbon chain until the final product is obtained. We study and propose the catalytic mechanism of the KR domain determined using the hybrid QM/MM methodology, at the ONIOM(B3LYP/6-311+G(2d,2p):AMBER) level of theory. The results indicate that the reaction mechanism occurs in two stages: (i) nucleophilic attack by a NADPH hydride to the β-carbon of the substrate, together with an asynchronous deprotonation of the Tyr2034 by the oxygen of the β-alkoxide to hold the final alcohol product; and (ii) an asynchronous deprotonation of the hydroxyl in the NADP⁺'s ribose by Tyr2034, and of the Lys1995 by the resulting alkoxide in the former ribose to restore the protonation state of Tyr2034. The reduction step occurs with a Gibbs energy barrier of 11.7 kcal mol^-1 and a Gibbs reaction energy of -10.6 kcal mol^-1. These results have provided an understanding of the catalytic mechanism of the KR hFAS domain, a piece of the heavy hFAS biosynthetic machinery.

Computational tools for the evaluation of laboratory-engineered biocatalysts

Biocatalysis is based on the application of natural catalysts for new purposes, for which enzymes were not designed. Although the first examples of biocatalysis were reported more than a century ago, biocatalysis was revolutionized after the discovery of an in vitro version of Darwinian evolution called Directed Evolution (DE). Despite the recent advances in the field, major challenges remain to be addressed. Currently, the best experimental approach consists of creating multiple mutations simultaneously while limiting the choices using statistical methods. Still, tens of thousands of variants need to be tested experimentally, and little information is available on how these mutations lead to enhanced enzyme proficiency. This review aims to provide a brief description of the available computational techniques to unveil the molecular basis of improved catalysis achieved by DE. An overview of the strengths and weaknesses of current computational strategies is explored with some recent representative examples. The understanding of how this powerful technique is able to obtain highly active variants is important for the future development of more robust computational methods to predict amino-acid changes needed for activity.