Classical QSAR and Docking Simulation of 4-Pyridone Derivatives for Their Antimalarial Activity

In this work, the minimum energy structures of 22 4-pyridone derivatives have been optimized at Density Functional Theory level, and several quantum molecular, including electronic and thermodynamic descriptors, were computed for these substrates in order to obtain a statistical and meaningful QSAR equation. In this sense, by using multiple linear regressions, five mathematical models have been obtained. The best model with only four descriptors (r2 = 0.86, Q2 = 0.92, S.E.P = 0.38) was validated by the leave-one-out cross-validation method. The antimalarial activity can be explained by the combination of the four mentioned descriptors e.g., electronic potential, dipolar momentum, partition coefficient and molar refractivity. The statistical parameters of this model suggest that it is robust enough to predict the antimalarial activity of new possible compounds; consequently, three small chemical modifications into the structural core of these compounds were performed specifically on the most active compound of the series (compound 13). These three new suggested compounds were leveled as 13A, 13B and 13C, and the predicted biological antimalarial activity is 0.02 µM, 0.03 µM, and 0.07 µM, respectively. In order to complement these results focused on the possible action mechanism of the substrates, a docking simulation was included for these new structures as well as for the compound 13 and the docking scores (binding affinity) obtained for the interaction of these substrates with the cytochrome bc1, were −7.5, −7.2, −6.9 and −7.5 kcal/mol for 13A, 13B, 13C and compound 13, respectively, which suggests that these compounds are good candidates for its biological application in this illness.

[1]  J. Vidal,et al.  Potent antimalarial 4-pyridones with improved physico-chemical properties. , 2011, Bioorganic & medicinal chemistry letters.

[2]  M. Mather,et al.  Mitochondrial drug targets in apicomplexan parasites. , 2007, Current drug targets.

[3]  James M. Woolven,et al.  Synthesis and structure-activity relationships of 4-pyridones as potential antimalarials. , 2008, Journal of medicinal chemistry.

[4]  Gene H. Golub,et al.  Generalized cross-validation as a method for choosing a good ridge parameter , 1979, Milestones in Matrix Computation.

[5]  Molar refractivity and interactions in solutions I. Molar refractivity of some monovalent ions in aqueous and dimethyl sulfoxide solutions , 2012 .

[6]  N. Meanwell Synopsis of some recent tactical application of bioisosteres in drug design. , 2011, Journal of medicinal chemistry.

[7]  P. Geerlings,et al.  Conceptual density functional theory. , 2003, Chemical reviews.

[8]  V. Barton,et al.  Inhibiting Plasmodium cytochrome bc1: a complex issue. , 2010, Current opinion in chemical biology.

[9]  T. Egan Interactions of quinoline antimalarials with hematin in solution. , 2006, Journal of inorganic biochemistry.

[10]  E. LaVoie,et al.  Bioisosterism: A Rational Approach in Drug Design. , 1996, Chemical reviews.

[11]  P. Pascutti,et al.  Insights into cytochrome bc 1 complex binding mode of antimalarial 2-hydroxy-1,4-naphthoquinones through molecular modelling , 2017, Memorias do Instituto Oswaldo Cruz.

[12]  Dennis E. Grawoig,et al.  Statistics, a foundation for analysis , 1972, The Mathematical Gazette.

[13]  David A. Winkler,et al.  The role of quantitative structure-activity relationships (QSAR) in biomolecular discovery , 2002, Briefings Bioinform..

[14]  Nancy Fullman,et al.  Global malaria mortality between 1980 and 2010: a systematic analysis , 2012, The Lancet.

[15]  C. Davis,et al.  Preclinical drug metabolism and pharmacokinetic evaluation of GW844520, a novel anti-malarial mitochondrial electron transport inhibitor. , 2006, Journal of pharmaceutical sciences.

[16]  R. C. Young,et al.  Dipole moment in relation to H2 receptor histamine antagonist activity for cimetidine analogues. , 1986, Journal of medicinal chemistry.

[17]  E. Márquez,et al.  Molecular modeling studies of bromopyrrole alkaloids as potential antimalarial compounds: a DFT approach , 2018, Medicinal Chemistry Research.

[18]  B. Trumpower,et al.  Specific roles of protein–phospholipid interactions in the yeast cytochrome bc1 complex structure , 2001, The EMBO journal.

[19]  Joanne M. Morrisey,et al.  Specific role of mitochondrial electron transport in blood-stage Plasmodium falciparum , 2007, Nature.

[20]  S. Meshnick,et al.  Modeling the molecular basis of atovaquone resistance in parasites and pathogenic fungi. , 2007, Trends in parasitology.

[21]  Martin Hofmann-Apitius,et al.  In silico drug discovery approaches on grid computing infrastructures. , 2010, Current clinical pharmacology.

[22]  D. Warhurst,et al.  Antimalarial drug discovery: development of inhibitors of dihydrofolate reductase active in drug resistance , 1998 .

[23]  S. Meshnick,et al.  Molecular Basis for Atovaquone Binding to the Cytochrome bc1 Complex* , 2003, Journal of Biological Chemistry.

[24]  A. Toro‐Labbé,et al.  Relations between Potential Energy, Electronic Chemical Potential, and Hardness Profiles , 1997 .

[25]  W J Dunn,et al.  Linear relationships between lipophilic character and biological activity of drugs. , 1972, Journal of pharmaceutical sciences.

[26]  K. Roy,et al.  Advances in quantitative structure–activity relationship models of antimalarials , 2010, Expert opinion on drug discovery.

[27]  BioChem Press,et al.  Predictive Comparative QSAR Modeling of 4-Pyridones as Potent Antimalarials , 2010 .

[28]  L. Mammino,et al.  Computational study of naphthylisoquinoline alkaloids with antimalarial activity from Dioncophyllaceae and Ancistrodaceae in vacuo , 2016, Theoretical Chemistry Accounts.

[29]  T. J. Villalobos,et al.  2D, 3D-QSAR and molecular docking of 4(1H)-quinolones analogues with antimalarial activities. , 2013 .

[30]  G. Moreau,et al.  Rational approach to drug design , 2008 .

[31]  Andrew Owen,et al.  Acridinediones: Selective and Potent Inhibitors of the Malaria Parasite Mitochondrial bc1 Complex , 2008, Molecular Pharmacology.

[32]  R. Mannhold,et al.  Calculation of molecular lipophilicity: state of the art and comparison of methods on more than 96000 compounds , 2009, Journal of pharmaceutical sciences.

[33]  Jonathan L. Vennerstrom,et al.  Structural Analysis of Chloroquine Resistance Reversal by Imipramine Analogs , 2001, Antimicrobial Agents and Chemotherapy.

[34]  R. Parr,et al.  Principle of maximum hardness , 1991 .

[35]  Biological Activity of Chemical Compounds and Their Molecular Structure-Information Approach , 2014 .

[36]  R. Pearson Chemical hardness and the electronic chemical potential , 1992 .

[37]  José R. Mora,et al.  Theoretical calculations on the mechanisms of the gas phase elimination kinetics of chlorocyclohexane, 3-chlorocyclohexene and 4-chlorocyclohexene , 2009 .

[38]  Michael J. Frisch,et al.  The performance of the Becke-Lee-Yang-Parr (B-LYP) density functional theory with various basis sets , 1992 .

[39]  J. A. Paixão,et al.  Quinolone-Hydroxyquinoline Tautomerism in Quinolone 3-Esters. Preserving the 4-Oxoquinoline Structure To Retain Antimalarial Activity. , 2015, The Journal of organic chemistry.

[40]  Wei-Chun Kao,et al.  Structural analysis of atovaquone-inhibited cytochrome bc1 complex reveals the molecular basis of antimalarial drug action , 2014, Nature Communications.

[41]  K. Kita,et al.  Parasite Mitochondria as a Target of Chemotherapy: Inhibitory Effect of Licochalcone A on the Plasmodium falciparum Respiratory Chain , 2005, Annals of the New York Academy of Sciences.