The global motion affecting electron transfer in Plasmodium falciparum type II NADH dehydrogenases: a novel non-competitive mechanism for quinoline ketone derivative inhibitors

Reactions of Plasmodium falciparum Type II NADH: Ubiquinone Oxidoreductase with Nonphysiological Quinoidal and Nitroaromatic Oxidants

In order to detail the antiplasmodial effects of quinones (Q) and nitroaromatic compounds (ArNO2), we investigated their reduction mechanism by Plasmodium falciparum flavoenzyme type II NADH:ubiquinone oxidoreductase (PfNDH2). The reactivity of Q and ArNO2 (n = 29) follows a common trend and exhibits a parabolic dependence on their single-electron reduction potential (E71), albeit with significantly scattered data. The reactivity of quinones with similar E71 values increases with their lipophilicity. Quinones are reduced by PfNDH2 in a two-electron way, but ArNO2 are reduced in a single-electron way. The inhibition studies using NAD+ and ADP-ribose showed that quinones oxidize the complexes of reduced enzyme with NADH and NAD+. This suggests that, as in the case of other NDH2s, quinones and the nicotinamide ring of NAD(H) bind at separate sites. A scheme of PfNDH2 catalysis is proposed, consistent with both the observed 'ping-pong' mechanism and the presence of two substrate binding sites. Molecular docking showed that Q and ArNO2 bind in a similar manner and that lipophilic quinones have a higher affinity for the binding site. One may expect that PfNDH2 can be partially responsible for the previously observed enhanced antiplasmodial activity of aziridinylbenzoquinones caused by their two-electron reduction, as well as for the redox cycling and oxidative stress-type action of ArNO2.

Feasibility of the inhibitor development for cancer: A systematic approach for drug design

The traditional Chinese medicine (TCM) bupleurum-ginger-licorice formula presents significant anti-cancer effects, but its active ingredients and inhibitory mechanism remain unclear. In this work, the core effective ingredient quercetin and its signal transducer and activator of transcription 3 (Stat3) receptor both were identified by network pharmacology. Quercetin is a low-toxicity, non-carcinogenic flavonoid with antioxidant, anti-inflammatory and anticancer activities, which is widely distributed in edible plants. Stat3 can bind to specific DNA response elements and serves as a transcription factor to promote the translation of some invasion/migration-related target genes, considered as a potential anticancer target. Here, molecular docking and molecular dynamics (MD) simulation both were used to explore molecular recognition of quercetin with Stat3. The results show that quercetin impairs DNA transcription efficiency by hindering Stat3 dimerization, partially destroying DNA conformation. Specifically, when the ligand occupies the SH2 cavity of the enzyme, spatial rejection is not conductive to phosphokinase binding. It indirectly prevents the phosphorylation of Y705 and the formation of Stat3 dimer. When the inhibitor binds to the DT1005 position, it obviously shortens the distance between DNA and DBD, enhances their binding capacity, and thereby reduces the degree of freedom required for transcription. This work not only provides the binding modes between Stat3 and quercetin, but also contributes to the optimization and design of such anti-cancer inhibitors.

S288T mutation altering MmpL3 periplasmic domain channel and H-bond network: a novel dual drug resistance mechanism

CONTEXT

Mycobacterial membrane proteins Large 3 (MmpL3) is responsible for the transport of mycobacterial acids out of cell membrane to form cell wall, which is essential for the survival of Mycobacterium tuberculosis (Mtb) and has become a potent anti-tuberculosis target. SQ109 is an ethambutol (EMB) analogue, as a novel anti-tuberculosis drug, can effectively inhibit MmpL3, and has completed phase 2b-3 clinical trials. Drug resistance has always been the bottleneck problem in clinical treatment of tuberculosis. The S288T mutant of MmpL3 shows significant resistance to the inhibitor SQ109, while the specific action mechanism remains unclear. The results show that MmpL3 S288T mutation causes local conformational change with little effect on the global structure. With MmpL3 bound by SQ109 inhibitor, the distance between D710 and R715 increases resulting in H-bond destruction, but their interactions and proton transfer function are still restored. In addition, the rotation of Y44 in the S288T mutant leads to an obvious bend in the periplasmic domain channel and an increased number of contact residues, reducing substrate transport efficiency. This work not only provides a possible dual drug resistance mechanism of MmpL3 S288T mutant but also aids the development of novel anti-tuberculosis inhibitors.

METHODS

In this work, molecular dynamics (MD) and quantum mechanics (QM) simulations both were performed to compare inhibitor (i.e., SQ109) recognition, motion characteristics, and H-bond energy change of MmpL3 after S288T mutation. In addition, the WT_SQ109 complex structure was obtained by molecular docking program (Autodock 4.2); Molecular Mechanics/ Poisson Boltzmann Surface Area (MM-PBSA) and Solvated Interaction Energy (SIE) methods were used to calculate the binding free energies (∆Gbind); Geometric criteria were used to analyze the changes of hydrogen bond networks.

Uncovering the interactions between PME and PMEI at the gene and protein levels: Implications for the design of specific PMEI

CONTEXT

Pectin methylesterase inhibitor (PMEI) can specifically bind and inhibit the activity of pectin methylesterase (PME), which has been widely used in fruit and vegetable juice processing. However, the limited three-dimensional structure, unclear action mechanism, low thermal stability and biological activity of PMEI severely limited its application. In this work, molecular recognition and conformational changes of PME and PMEI were analyzed by various molecular simulation methods. Then suggestions were proposed for improving thermal stability and affinity maturation of PMEI through semi-rational design.

METHODS

Phylogenetic trees of PME and PMEI were established using the Maximum likelihood (ML) method. The results show that PME and PMEI have good sequence and structure conservation in various plants, and the simulated data can be widely adopted. In this work, MD simulations were performed using AMBER20 package and ff14SB force field. Protein interaction analysis indicates that H-bonds, van der Waals forces, and the salt bridge formed of K224 with ID116 are the main driving forces for mutual molecular recognition of PME and PMEI. According to the analyses of free energy landscape (FEL), conformational cluster, and motion, the association with PMEI greatly disrupts PME's dispersed functional motion mode and biological function. By monitoring the changes of residue contact number and binding free energy, IG35M/ IG35R: IT93F and IT113W/ IT113W: ID116W mutations contribute to thermal stability and affinity maturation of the PME-PMEI complex system, respectively. This work reveals the interaction between PME and PMEI at the gene and protein levels and provides options for modifying specific PMEI.

Quantitative Estimation of the Hydrogen-Atom-Donating Ability of 4-Substituted Hantzsch Ester Radical Cations

The purpose of this study is to investigate thermodynamic and kinetic properties on the hydrogen-atom-donating ability of 4-substituted Hantzsch ester radical cations (XRH•+), which are excellent NADH coenzyme models. Gibbs free energy changes and activation free energies of 17 XRH•+ releasing H• [denoted as ΔG HD o(XRH•+) and ΔG HD ≠(XRH•+)] were calculated using density functional theory (DFT) and compared with that of Hantzsch ester (HEH2) and NADH. ΔG HD o(XRH•+) range from 19.35 to 31.25 kcal/mol, significantly lower than that of common antioxidants (such as ascorbic acid, BHT, the NADH coenzyme, and so forth). ΔG HD ≠(XRH•+) range from 29.81 to 39.00 kcal/mol, indicating that XRH•+ spontaneously releasing H• are extremely slow unless catalysts or active intermediate radicals exist. According to the computed data, it can be inferred that the Gibbs free energies and activation free energies of the core 1,4-dihydropyridine radical cation structure (DPH•+) releasing H• [ΔG HD o(DPH•+) and ΔG HD ≠(DPH•+)] should be 19-32 kcal/mol and 29-39 kcal/mol in acetonitrile, respectively. The correlations between the thermodynamic driving force [ΔG HD o(XRH•+)] and the activation free energy [ΔG HD ≠(XRH•+)] are also explored. Gibbs free energy is the important and decisive parameter, and ΔG HD ≠(XRH•+) increases in company with the increase of ΔG HD o(XRH•+), but no simple linear correlations are found. Even though all XRH•+ are judged as excellent antioxidants from the thermodynamic view, the computed data indicate that whether XRH•+ is an excellent antioxidant in reaction is decided by the R substituents in 4-position. XRH•+ with nonaromatic substituents tend to release R• instead of H• to quench radicals. XRH•+ with aromatic substituents tend to release H• and be used as antioxidants, but not all aromatic substituted Hantzsch esters are excellent antioxidants.

A High-Efficiency Artificial Synthetic Pathway for 5-Aminovalerate Production From Biobased L-Lysine in Escherichia coli

Bioproduction of 5-aminovalerate (5AVA) from renewable feedstock can support a sustainable biorefinery process to produce bioplastics, such as nylon 5 and nylon 56. In order to achieve the biobased production of 5AVA, a 2-keto-6-aminocaproate-mediated synthetic pathway was established. Combination of L-Lysine α-oxidase from Scomber japonicus, α-ketoacid decarboxylase from Lactococcus lactis and aldehyde dehydrogenase from Escherichia coli could achieve the biosynthesis of 5AVA from biobased L-Lysine in E. coli. The H₂O₂ produced by L-Lysine α-oxidase was decomposed by the expression of catalase KatE. Finally, 52.24 g/L of 5AVA were obtained through fed-batch biotransformation. Moreover, homology modeling, molecular docking and molecular dynamic simulation analyses were used to identify mutation sites and propose a possible trait-improvement strategy: the expanded catalytic channel of mutant and more hydrogen bonds formed might be beneficial for the substrates stretch. In summary, we have developed a promising artificial pathway for efficient 5AVA synthesis.

Inhibition Mechanism of Indoleamine 2, 3-Dioxygenase 1 (IDO1) by Amidoxime Derivatives and Its Revelation in Drug Design: Comparative Molecular Dynamics Simulations

For cancer treatment, in addition to the three standard therapies of surgery, chemotherapy, and radiotherapy, immunotherapy has become the fourth internationally-recognized alternative treatment. Indoleamine 2, 3-dioxygenase 1 (IDO1) catalyzes the conversion of tryptophan to kynurenine causing lysine depletion, which is an important target in the research and development of anticancer drugs. Epacadostat (INCB024360) is currently one of the most potent IDO1 inhibitors, nevertheless its inhibition mechanism still remains elusive. In this work, comparative molecular dynamics simulations were performed to reveal that the high inhibitory activity of INCB024360 mainly comes from two aspects: disturbing the ligand delivery tunnel and then preventing small molecules such as oxygen and water molecules from accessing the active site, as well as hindering the shuttle of substrate tryptophan with product kynurenine through the heme binding pocket. The scanning of key residues showed that L234 and R231 residues both were crucial to the catalytic activity of IDO1. With the association with INCB024360, L234 forms a stable hydrogen bond with G262, which significantly affects the spatial position of G262-A264 loop and then greatly disturbs the orderliness of ligand delivery tunnel. In addition, the cleavage of hydrogen bond between G380 and R231 increases the mobility of the GTGG conserved region, leading to the closure of the substrate tryptophan channel. This work provides new ideas for understanding action mechanism of amidoxime derivatives, improving its inhibitor activity and developing novel inhibitors of IDO1.