Probing the B- & C-rings of the antimalarial tetrahydro-β-carboline MMV008138 for steric and conformational constraints

The antimalarial candidate MMV008138 (1a) is of particular interest because its target enzyme (IspD) is absent in human. To achieve higher potency, and to probe for steric demand, a series of analogs of 1a were prepared that featured methyl-substitution of the B- and C-rings, as well as ring-chain transformations. X-ray crystallography, NMR spectroscopy and calculation were used to study the effects of these modifications on the conformation of the C-ring and orientation of the D-ring. Unfortunately, all the B- and C-ring analogs explored lost in vitro antimalarial activity. The possible role of steric effects and conformational changes on target engagement are discussed.     

3D-QSAR pharmacophore-based virtual screening,molecular docking and molecular dynamics simulation toward identifying lead compounds for NS2B–NS3 protease inhibitors

NS2B–NS3 protease has been identified to serve as lead drug design target due to its significant role in West Nile viral (WNV) and dengue virus (DENV) reproduction and replication. There are currently no approved chemotherapeutic drugs and effective vaccines to inhibit DENV and WNV infections. In this work, 3D-QSAR pharmacophore model has been developed to discover potential inhibitory candidates. Validation through Fischer’s model and decoy test indicate that the developed 3D pharmacophore model is highly predictive for DENV inhibitors, which was then employed to screen ZINC chemical library to obtain reasonable hits. Following ADMET filtering, 15 hits were subjected to further filter through molecular docking and CoMFA modeling. Finally, top three hits were identified as lead compounds or potential inhibitory candidates with IC₅₀ values of ～0.4637?µM and fitness of ～57.73. It is implied from CoMFA modeling that substituents at the side site of benzotriazole such as a p-nitro group (e.g. biphenyl head) and a carbonyl (e.g. carboxylate function) at the side site of furan or amino group may improve bioactivity of ZINC85645245, respectively. Molecular dynamics simulations (MDS) were performed to discover new interactions and reinforce the binding modes from docking for the hits also. The QSAR and MDS results obtained from this work should be useful in determining structural requirements for inhibitor development as well as in designing more potential inhibitors for NS2B–NS3 protease.     

Moleculer dynamics simulaiton revealed reciever domain of Acinetobacter baumannii BfmR enzyme as the hot spot for future antibiotics designing

Acinetobacter baumannii is an alarming nosocomial pathogen that is resistant to multiple drugs. The pathogen is forefront of scientific attention because of high mortality and morbidity found for its complications in the past decade. As a consequence, identification of novel drug candidates and subsequent designing of novel chemical scaffolds is an imperative need of time. In the present study, we used a recently reported structure of BfmR enzyme and performed structure based virtual screening, MD simulation and binding free energies calculations. MD simulation revealed a profound movement of the best-characterized inhibitor towards the α4-β5-α5 face of the enzyme receiver domain, thus indicating its high affinity for this site compared to phosphorylation. Furthermore, it was observed that the enzyme and enzyme-inhibitor complex have high structure stability with mean RMSD of 1.2 and 1.1 Å, respectively. Binding free energy calculations for the complex unraveled high stability with MMGBSA score of ?26.21?kcal/mol and MMPBSA score of ?1.47?kcal/mol. Van der Waal energy was found highly favorable with value of ?30.25?kcal/mol and dominated significantly the overall binding energy. Furthermore, a novel WaterSwap assay was used to circumvent the limitations of MMGB/PBSA that complements the inhibitor affinity for enzyme active pocket as depicted by the low convergence of Bennett, TI and FEP algorithms. Results yielded from this study will not only give insight into the phenomena of inhibitor movement towards the enzyme receiver domain, but will also provide a useful baseline for designing derivatives with improved biological and pharmacokinetics profiles.

Communicated by Ramaswamy H. Sarma     

An update on microbial carotenoid production: application of recent metabolic engineering tools

Carotenoids are ubiquitous pigments synthesized by plants, fungi, algae, and bacteria. Industrially, carotenoids are used in pharmaceuticals, neutraceuticals, and animal feed additives, as well as colorants in cosmetics and foods. Scientific interest in dietary carotenoids has increased in recent years because of their beneficial effects on human health, such as lowering the risk of cancer and enhancement of immune system function, which are attributed to their antioxidant potential. The availability of carotenoid genes from carotenogenic microbes has made possible the synthesis of carotenoids in non-carotenogenic microbes. The increasing interest in microbial sources of carotenoid is related to consumer preferences for natural additives and the potential cost effectiveness of creating carotenoids via microbial biotechnology. In this review, we will describe the recent progress made in metabolic engineering of non-carotenogenic microorganisms with particular focus on the potential of Escherichia coli for improved carotenoid productivity. Amitabha Das and Sang-Hwal Yoon contributed equally to this work.     

Genetic evidence for the role of isopentenyl diphosphate isomerases in the mevalonate pathway and plant development in Arabidopsis

Isopentenyl/dimethylallyl diphosphate isomerase (IPI) catalyzes the interconversion of isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP), which are the universal C(5) units of isoprenoids. In plants, IPP and DMAPP are synthesized via the cytosolic mevalonate (MVA) and plastidic methylerythritol phosphate (MEP) pathways, respectively. However, the role of IPI in each pathway and in plant development is unknown due to a lack of genetic studies using IPI-defective mutants. Here, we show that the atipi1atipi2 double mutant, which is defective in two Arabidopsis IPI isozymes, exhibits dwarfism and male sterility under long-day conditions and decreased pigmentation under continuous light, whereas the atipi1 and atipi2 single mutants are phenotypically normal. We also show that the sterol and ubiquinone levels in the double mutant are <50% of those in wild-type plants, and that the male-sterile phenotype is chemically complemented by squalene, a sterol precursor. In vivo isotope labeling experiments using the atipi1atipi2 double mutant revealed a decrease in the incorporation of MVA (in its lactone form) into sterols, with no decrease in the incorporation of MEP pathway intermediates into tocopherol. These results demonstrate a critical role for IPI in isoprenoid biosynthesis via the MVA pathway, and they imply that IPI is essential for the maintenance of appropriate levels of IPP and DMAPP in different subcellular compartments in plants.     

Screening for the selective inhibitors of MMP-9 from natural products based on pharmacophore modeling and molecular docking in combination with bioassay experiment,hybrid QM/MM calculation,and MD simulation

Matrix metalloproteinase-9 (MMP-9) has been considered as an attractive target involving cancer therapy. In this study, the 3D QSAR pharmacophore model of MMP-9 inhibitors is built, and its reliability is subsequently validated based on different methods. The built pharmacophore model consists of the four chemical features, including two hydrogen bond acceptors (HBA), one hydrophobic (HY), and one ring aromatic (RA). Among them, both HY and RA are found to be especially important features because they involve the interactions of inhibitors with the S1′ pocket of MMP-9, which determines the selectivity of MMP-9 inhibitors. By combining pharmacophore model with molecular docking, the virtual screening is carried out to identify the selective MMP-9 inhibitors from natural products. The four potential selective MMP-9 inhibitors of natural products are found. One of them was used to carry out the bioassay experiment inhibiting MMP-9, and the estimated IC₅₀ value of only 26.94 µM clearly shows its strongly inhibitory activity; besides, both the hybrid quantum mechanics/molecular mechanics (QM/MM) calculation and the molecular dynamics simulation are performed to examine the reliability regarding the binding mode of this inhibitor with MMP-9 active sites predicted by molecular docking. All the screened four natural products are found to well bind with the MMP-9 active sites by different kinds of interactions. Finally, the ADMET properties of screened four natural products are assessed. These screened MMP-9 inhibitors of natural products could be used as the lead compounds to perform structural modifications and optimizations in the future work.

Communicated by Ramaswamy H. Sarma     

88 Comparison of monovalent and divalent ion distributions around a DNA duplex with molecular dynamic simulation and Poisson–Boltzmann approach

Ion interactions with nucleic acids (both DNA and RNA) are an important and evolving field of investigation. Positively charged cations may interact with highly negatively charged nucleic acids via simple electrostatic interactions to help screen the electrostatic repulsion along the nucleic acids and assist their folding and/or compaction. Cations may also bind at specific sites and become integral parts of the structures, possibly playing important enzymatic roles. Two popular methods for computationally exploring a nucleic acid’s ion atmosphere are atomistic molecular dynamics (MD) simulations and the Poisson–Boltzmann (PB) equation. In general, monovalent ion results obtained from MD simulations and the PB equation agree well with experiment. However, Bai et al. (2007) observed discrepancies between experiment and the PB equation while examining the competitive binding of monovalent and divalent ions, with more significant discrepancies for divalent ions. The goal of this project was to thoroughly investigate monovalent (Na⁺) and divalent (Mg²⁺) ion distributions formed around a DNA duplex with MD simulations and the PB equation. We simulated three different cation concentrations, and matched the equilibrated bulk ion concentration for our theoretical calculations with the PB equation. Based on previous work, our Mg²⁺ ions were fully solvated, the expected state of Mg²⁺ ions when interacting with a duplex, when the production simulations began and remained throughout the simulations (Kirmizialtin, 2010; Robbins, 2012). Na⁺ ion distributions and number of Na⁺ ions within 10?Å of the DNA obtained from our two methods agreed well. However, results differed for Mg²⁺ ions, with a lower number of ions within the cut-off distance obtained from the PB equation when compared to MD simulations. The Mg²⁺ ion distributions around the DNA obtained via the two methods also differed. Based on our results, we conclude that the PB equation will systematically underestimate Mg²⁺ ions bound to DNA, and much of this deviation is attributed to dielectric saturation associated with high valency ions.     

In-Silico molecular docking and simulation studies on novel chalcone and flavone hybrid derivatives with 1, 2, 3-triazole linkage as vital inhibitors of Plasmodium falciparum dihydroorotate dehydrogenase

The structural motifs of chalcones, flavones, and triazoles with varied substitutions have been studied for the antimalarial activity. In this study, 25 novel derivatives of chalcone and flavone hybrid derivatives with 1, 2, 3-triazole linkage are docked with Plasmodium falciparum dihydroorotate dehydrogenase to establish their inhibitory activity against Plasmodium falciparum. The best binding conformation of the ligands at the catalytic site of dihydroorotate dehydrogenase are selected to characterize the best bound ligand using the best consensus score and the number of hydrogen bond interactions. The ligand namely (2E)-3-(4-{[1-(3-chloro-4-fluorophenyl)-1H-1, 2, 3-triazol-4-yl]methoxy}-3-methoxyphenyl-1-(2-hydroxy-4,6-dimethoxyphenyl)prop-2-en-1-one, is one the among the five best docked ligands, which interacts with the protein through nine hydrogen bonds and with a consensus score of five. To refine and confirm the docking study results, the stability of complexes is verified using Molecular Dynamics Simulations, Molecular Mechanics /Poisson–Boltzmann Surface Area free binding energy analysis, and per residue contribution for the binding energy. The study implies that the best docked Plasmodium falciparum dihydroorotate dehydrogenase–ligand complex is having high negative binding energy, most stable, compact, and rigid with nine hydrogen bonds. The study provides insight for the optimization of chalcone and flavone hybrids with 1, 2, 3-triazole linkage as potent inhibitors.