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

Reconstructing dynamics correlation network to simultaneously improve activity and stability of 2,3-butanediol dehydrogenase by design of distal interchain disulfide bonds

The complete enzyme catalytic cycle includes substrate binding, chemical reaction and product release, in which different dynamic conformations are adopted. Due to the complex relationship among enzyme activity, stability and dynamics, the directed evolution of enzymes for improved activity or stability commonly leads to a trade-off in stability or activity. It hence remains a challenge to engineer an enzyme to have both enhanced activity and stability. Here, we have attempted to reconstruct the dynamics correlation network involved with active center to improve both activity and stability of a 2,3-butanediol dehydrogenase (2,3-BDH) by introducing inter-chain disulfide bonds. A computational strategy was first applied to evaluate the effect of introducing inter-chain disulfide bond on activity and stability of three 2,3-BDHs, and the N258C mutation of 2,3-BDH from Corynebacterium glutamicum (CgBDH) was proved to be effective in improving both activity and stability. In the results, CgBDH-N258C showed a different unfolding curve from the wild type, with two melting temperatures (Tm) of 68.3 °C and 50.8 °C, 19.7 °C and 2 °C higher than 48.6 °C of the wild type. Its half-life was also improved by 14.8-fold compared to the wild type. Catalytic efficiency (kcat/Km) of the mutant was increased by 7.9-fold toward native substrate diacetyl and 8.8-fold toward non-native substrate 2,5-hexanedione compared to the wild type. Molecular dynamics simulations revealed that an interaction network formed by Cys258, Arg162, Ala144 and the catalytic residues was reconstructed in the mutant and the dynamics change caused by the disulfide bond could be propagated through the interactions network. This improved the enzyme stability and activity by decreasing the flexibility and locking more "reactive" pose, respectively. Further construction of mutations including A144G showing a 44-fold improvement in catalytic efficiency toward meso-2,3-BD confirmed the role of modifying dynamics correlation network in tunning enzyme activity and selectivity. This study provided important insights into the relationship among dynamics, enzyme catalysis and stability, and will be useful in the designing new enzymes with co-evolution of stability, activity and selectivity.

Multiscale Workflow for Modeling Ligand Complexes of Zinc Metalloproteins

Zinc metalloproteins are ubiquitous, with protein zinc centers of structural and functional importance, involved in interactions with ligands and substrates and often of pharmacological interest. Biomolecular simulations are increasingly prominent in investigations of protein structure, dynamics, ligand interactions, and catalysis, but zinc poses a particular challenge, in part because of its versatile, flexible coordination. A computational workflow generating reliable models of ligand complexes of biological zinc centers would find broad application. Here, we evaluate the ability of alternative treatments, using (nonbonded) molecular mechanics (MM) and quantum mechanics/molecular mechanics (QM/MM) at semiempirical (DFTB3) and density functional theory (DFT) levels of theory, to describe the zinc centers of ligand complexes of six metalloenzyme systems differing in coordination geometries, zinc stoichiometries (mono- and dinuclear), and the nature of interacting groups (specifically the presence of zinc-sulfur interactions). MM molecular dynamics (MD) simulations can overfavor octahedral geometries, introducing additional water molecules to the zinc coordination shell, but this can be rectified by subsequent semiempirical (DFTB3) QM/MM MD simulations. B3LYP/MM geometry optimization further improved the accuracy of the description of coordination distances, with the overall effectiveness of the approach depending upon factors, including the presence of zinc-sulfur interactions that are less well described by semiempirical methods. We describe a workflow comprising QM/MM MD using DFTB3 followed by QM/MM geometry optimization using DFT (e.g., B3LYP) that well describes our set of zinc metalloenzyme complexes and is likely to be suitable for creating accurate models of zinc protein complexes when structural information is more limited.

Characterisation of putative lactate synthetic pathways of Coxiella burnetii

The zoonotic pathogen Coxiella burnetii, the causative agent of the human disease Q fever, is an ever-present danger to global public health. Investigating novel metabolic pathways necessary for C. burnetii to replicate within its unusual intracellular niche may identify new therapeutic targets. Recent studies employing stable isotope labelling established the ability of C. burnetii to synthesize lactate, despite the absence of an annotated synthetic pathway on its genome. A noncanonical lactate synthesis pathway could provide a novel anti-Coxiella target if it is essential for C. burnetii pathogenesis. In this study, two C. burnetii proteins, CBU1241 and CBU0823, were chosen for analysis based on their similarities to known lactate synthesizing enzymes. Recombinant GST-CBU1241, a putative malate dehydrogenase (MDH), did not produce measurable lactate in in vitro lactate dehydrogenase (LDH) activity assays and was confirmed to function as an MDH. Recombinant 6xHis-CBU0823, a putative NAD+-dependent malic enzyme, was shown to have both malic enzyme activity and MDH activity, however, did not produce measurable lactate in either LDH or malolactic enzyme activity assays in vitro. To examine potential lactate production by CBU0823 more directly, [13C]glucose labelling experiments compared label enrichment within metabolic pathways of a cbu0823 transposon mutant and the parent strain. No difference in lactate production was observed, but the loss of CBU0823 significantly reduced 13C-incorporation into glycolytic and TCA cycle intermediates. This disruption to central carbon metabolism did not have any apparent impact on intracellular replication within THP-1 cells. This research provides new information about the mechanism of lactate biosynthesis within C. burnetii, demonstrating that CBU1241 is not multifunctional, at least in vitro, and that CBU0823 also does not synthesize lactate. Although critical for normal central carbon metabolism of C. burnetii, loss of CBU0823 did not significantly impair replication of the bacterium inside cells.

Potential biomarkers in septic shock besides lactate

IMPACT STATEMENT

Elevated lactate has been commonly considered as a biomarker and a useful prognostic tool for resuscitation in septic shock, facilitating physician more rapid intervention and treatment. However, it can be initiated by hypoxia, but persistent hyperlactatemia may not represent persistent hypoxia only. In the article, it is the first time to review potential biomarkers in septic shock from the point of view of energy metabolism including intermediates of TCA cycle, MAS, the NAD⁺/NADH ratio, NAD⁺, NADH, malate, and MDH. And the combination of lactate and MDH is also proposed in septic shock for the first time, as MDH in cytoplasm and mitochondria participates in both MAS and TCA cycle for ATP generation. Its feasibility in clinic has been analyzed at the end, although related research is still limited. It is reasonable the combination of lactate and MDH will be more comprehensive to reflex hypoxia in septic shock.