1
|
Ren R, Gao L, Li G, Wang S, Zhao Y, Wang H, Liu J. 2D, 3D-QSAR study and docking of vascular endothelial growth factor receptor 3 (VEGFR3) inhibitors for potential treatment of retinoblastoma. Front Pharmacol 2023; 14:1177282. [PMID: 37089961 PMCID: PMC10119426 DOI: 10.3389/fphar.2023.1177282] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2023] [Accepted: 03/24/2023] [Indexed: 04/25/2023] Open
Abstract
Background: Retinoblastoma is currently the most common malignant tumor seen in newborns and children's eyes worldwide, posing a life-threatening hazard. Chemotherapy is an integral part of retinoblastoma treatment. However, the chemotherapeutic agents used in clinics often lead to drug resistance. Thus there is a need to investigate new chemotherapy-targeted agents. VEGFR3 inhibitors are anti-tumour-growth and could be used to develop novel retinoblastoma-targeted agents. Objective: To predict drug activity, discover influencing factors and design new drugs by building 2D, 3D-QSAR models. Method: First, linear and non-linear QSAR models were built using heuristic methods and gene expression programming (GEP). The comparative molecular similarity indices analysis (COMISA) was then used to construct 3D-QSAR models through the SYBYL software. New drugs were designed by changing drug activity factors in both models, and molecular docking experiments were performed. Result: The best linear model created using HM had an R2, S2, and R2cv of 0.82, 0.02, and 0.77, respectively. For the training and test sets, the best non-linear model created using GEP had correlation coefficients of 0.83 and 0.72 with mean errors of 0.02 and 0.04. The 3D model designed using SYBYL passed external validation due to its high Q2 (0.503), R2 (0.805), and F-value (76.52), as well as its low standard error of SEE value (0.172). This demonstrates the model's reliability and excellent predictive ability. Based on the molecular descriptors of the 2D model and the contour plots of the 3D model, we designed 100 new compounds using the best active compound 14 as a template. We performed activity prediction and molecular docking experiments on them, in which compound 14.d performed best regarding combined drug activity and docking ability. Conclusion: The non-linear model created using GEP was more stable and had a more substantial predictive power than the linear model built using the heuristic technique (HM). The compound 14.d designed in this experiment has the potential for anti-retinoblastoma treatment, which provides new design ideas and directions for retinoblastoma-targeted drugs.
Collapse
Affiliation(s)
- Rui Ren
- Affiliated Hospital of Weifang Medical University, School of Clinical Medicine, Weifang Medical University, Weifang, China
| | - Liyu Gao
- Affiliated Hospital of Weifang Medical University, School of Clinical Medicine, Weifang Medical University, Weifang, China
| | - Guoqi Li
- Affiliated Hospital of Weifang Medical University, School of Clinical Medicine, Weifang Medical University, Weifang, China
| | | | | | | | - Jianwei Liu
- Affiliated Hospital of Weifang Medical University, Weifang, China
- *Correspondence: Jianwei Liu,
| |
Collapse
|
2
|
Iwaloye O, Elekofehinti OO, Kikiowo B, Fadipe TM, Akinjiyan MO, Ariyo EO, Aiyeku OO, Adewumi NA. Discovery of Traditional Chinese Medicine Derived Compounds as Wild Type and Mutant Plasmodium falciparum Dihydrofolate Reductase Inhibitors: Induced Fit Docking and ADME Studies. Curr Drug Discov Technol 2021; 18:554-569. [PMID: 32729419 DOI: 10.2174/1570163817999200729122753] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2020] [Revised: 05/28/2020] [Accepted: 06/04/2020] [Indexed: 11/22/2022]
Abstract
BACKGROUND In a bid to come up with effective compounds as inhibitors for antimalarial treatment, we built a library of 2,000 traditional Chinese medicine(TCM)-derived compounds retrieved from TCM Database@Taiwan. METHODS The active sites of both the wild type and mutant Plasmodium falciparum dihydrofolatereductase (pfDHFR) were explored using computational tools. pfDHFR, one of the prime drug targets in the prevention of malaria infection induced by the female anopheles mosquito has continued to offer resistance to drugs (antifolates) due to mutation in some of the key amino acid residues crucial for its inhibition. RESULTS We utilized virtual throughput screening and glide XP docking to screen the compounds, and 8 compounds were found to have promising docking scores with both the wild type and mutant pfDHFR. They were further subjected to Induce Fit Docking (IFD) to affirm their inhibitory potency. The ADME properties and biological activity spectrum of the compounds were also considered. The inhibition profile of the compounds revealed that a number of compounds formed intermolecular interactions with ASP54, ILE14, LEU164, SER108/ASN108, ARG122 and ASP58. Most of the compounds can be considered as drug candidates due to their antiprotozoal activities and accordance with the Lipinski's Rule of Five (ROF). CONCLUSION The outcome of the present study should further be investigated to attest the efficacy of these compounds as better drug candidates than the antifolates.
Collapse
Affiliation(s)
- Opeyemi Iwaloye
- Department of Biochemistry, Bioinformatics and Molecular Biology Unit, Federal University of Technology Akure, Ondo State, Nigeria
| | - Olusola Olalekan Elekofehinti
- Department of Biochemistry, Bioinformatics and Molecular Biology Unit, Federal University of Technology Akure, Ondo State, Nigeria
| | - Babatomiwa Kikiowo
- Department of Biochemistry, Adekunle Ajasin University, Akungba Akoko, Ondo State, Nigeria
| | - Toyin Mary Fadipe
- Department of Bioscience, International Institute of Agriculture, Ibadan, Oyo State, Nigeria
| | - Moses Orimoloye Akinjiyan
- Department of Biochemistry, Bioinformatics and Molecular Biology Unit, Federal University of Technology Akure, Ondo State, Nigeria
| | - Esther Opeyemi Ariyo
- Department of Biochemistry, Bioinformatics and Molecular Biology Unit, Federal University of Technology Akure, Ondo State, Nigeria
| | - Olabisi Olapade Aiyeku
- Department of Biochemistry, Bioinformatics and Molecular Biology Unit, Federal University of Technology Akure, Ondo State, Nigeria
| | - Nicholas Adeyemi Adewumi
- Department of Biochemistry, Bioinformatics and Molecular Biology Unit, Federal University of Technology Akure, Ondo State, Nigeria
| |
Collapse
|
3
|
Two Decades of 4D-QSAR: A Dying Art or Staging a Comeback? Int J Mol Sci 2021; 22:ijms22105212. [PMID: 34069090 PMCID: PMC8156896 DOI: 10.3390/ijms22105212] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2021] [Revised: 05/11/2021] [Accepted: 05/12/2021] [Indexed: 01/01/2023] Open
Abstract
A key question confronting computational chemists concerns the preferable ligand geometry that fits complementarily into the receptor pocket. Typically, the postulated ‘bioactive’ 3D ligand conformation is constructed as a ‘sophisticated guess’ (unnecessarily geometry-optimized) mirroring the pharmacophore hypothesis—sometimes based on an erroneous prerequisite. Hence, 4D-QSAR scheme and its ‘dialects’ have been practically implemented as higher level of model abstraction that allows the examination of the multiple molecular conformation, orientation and protonation representation, respectively. Nearly a quarter of a century has passed since the eminent work of Hopfinger appeared on the stage; therefore the natural question occurs whether 4D-QSAR approach is still appealing to the scientific community? With no intention to be comprehensive, a review of the current state of art in the field of receptor-independent (RI) and receptor-dependent (RD) 4D-QSAR methodology is provided with a brief examination of the ‘mainstream’ algorithms. In fact, a myriad of 4D-QSAR methods have been implemented and applied practically for a diverse range of molecules. It seems that, 4D-QSAR approach has been experiencing a promising renaissance of interests that might be fuelled by the rising power of the graphics processing unit (GPU) clusters applied to full-atom MD-based simulations of the protein-ligand complexes.
Collapse
|
4
|
Chemoinformatics and QSAR. Adv Bioinformatics 2021. [DOI: 10.1007/978-981-33-6191-1_10] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022] Open
|
5
|
Sakkiah S, Guo W, Pan B, Kusko R, Tong W, Hong H. Computational prediction models for assessing endocrine disrupting potential of chemicals. JOURNAL OF ENVIRONMENTAL SCIENCE AND HEALTH. PART C, ENVIRONMENTAL CARCINOGENESIS & ECOTOXICOLOGY REVIEWS 2019; 36:192-218. [PMID: 30633647 DOI: 10.1080/10590501.2018.1537132] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Endocrine disrupting chemicals (EDCs) mimic natural hormones and disrupt endocrine function. Humans and wildlife are exposed to EDCs might alter endocrine functions through various mechanisms and lead to an adverse effects. Hence, EDCs identification is important to protect the ecosystem and to promote the public health. Leveraging in-vitro and in-vivo experiments to identify potential EDCs is time consuming and expensive. Hence, quantitative structure-activity relationship is applied to screen the potential EDCs. Here, we summarize the predictive models developed using various algorithms to forecast the binding activity of chemicals to the estrogen and androgen receptors, alpha-fetoprotein, and sex hormone binding globulin.
Collapse
Affiliation(s)
- Sugunadevi Sakkiah
- a Division of Bioinformatics and Biostatistics , National Center for Toxicological Research, U.S. Food and Drug Administration , Jefferson , Arkansas , USA
| | - Wenjing Guo
- a Division of Bioinformatics and Biostatistics , National Center for Toxicological Research, U.S. Food and Drug Administration , Jefferson , Arkansas , USA
| | - Bohu Pan
- a Division of Bioinformatics and Biostatistics , National Center for Toxicological Research, U.S. Food and Drug Administration , Jefferson , Arkansas , USA
| | - Rebecca Kusko
- b Immuneering Corporation , Cambridge , Massachusetts , USA
| | - Weida Tong
- a Division of Bioinformatics and Biostatistics , National Center for Toxicological Research, U.S. Food and Drug Administration , Jefferson , Arkansas , USA
| | - Huixiao Hong
- a Division of Bioinformatics and Biostatistics , National Center for Toxicological Research, U.S. Food and Drug Administration , Jefferson , Arkansas , USA
| |
Collapse
|
6
|
Abbat S, Jaladanki CK, Bharatam PV. Exploring PfDHFR reaction surface: A combined molecular dynamics and QM/MM analysis. J Mol Graph Model 2018; 87:76-88. [PMID: 30508692 DOI: 10.1016/j.jmgm.2018.11.007] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2018] [Revised: 11/16/2018] [Accepted: 11/19/2018] [Indexed: 11/18/2022]
Abstract
The substrate to the enzyme PfDHFR (Plasmodium falciparum Dihydrofolate Reductase) is a small molecule dihydrofolate (DHF), it gets converted to tetrahydrofolate (THF) in the active site of the enzyme. The PfDHFR reaction surface involves the protonation of DHF to DHFP as an initial step before the catalytic conversion. The binding affinities of all these species (DHF, DHFP and THF) contribute to the mechanism of DHFR catalytic action. Molecular dynamics (MD) simulations and Quantum Mechanics/Molecular Mechanics (QM/MM) analysis were performed to evaluate the binding affinity and molecular recognition interactions of the substrate DHF/DHFP and the product THF, in the active site of wild-type PfDHFR (wtPfDHFR). The binding affinities of the cofactor NADPH/NADP+ were also estimated in all the three complexes. The molecular dynamics (MD) simulations of the substrate, product and cofactor in the cavities of wtPfDHFR revealed the variation of the atomic level interactions during the course of the catalytic conversion. It was found that the DHFP binds very strongly to the PfDHFR active site and pulls the cofactor NADPH closer to itself. The QM/MM analysis revealed that the binding energy of DHFP (-59.82 kcal/mol) and NADPH (-100.24 kcal/mol) in DHFP-wtPfDHFR complex, is higher in comparison to the binding energy of DHF (-38.67 kcal/mol) and NADPH (-77.53 kcal/mol) in DHF-wtPfDHFR complex and the binding energy of THF (-30.72 kcal/mol) and NADP+ (-73.72 kcal/mol) in THF-wtPfDHFR complex. The hydride ion donor-acceptor distance (DAD) analysis was also carried out. This combined MD and QM/MM analysis revealed that the protonation of DHF increases the proximity between the substrate and the cofactor, thus facilitates the reaction profile of PfDHFR.
Collapse
Affiliation(s)
- Sheenu Abbat
- Department of Pharmacoinformatics, National Institute of Pharmaceutical Education and Research, Sector 67, S.A.S. Nagar, Punjab, 160 062, India
| | - Chaitanya K Jaladanki
- Department of Medicinal Chemistry, National Institute of Pharmaceutical Education and Research, Sector 67, S.A.S. Nagar, Punjab, 160 062, India
| | - Prasad V Bharatam
- Department of Pharmacoinformatics, National Institute of Pharmaceutical Education and Research, Sector 67, S.A.S. Nagar, Punjab, 160 062, India; Department of Medicinal Chemistry, National Institute of Pharmaceutical Education and Research, Sector 67, S.A.S. Nagar, Punjab, 160 062, India.
| |
Collapse
|
7
|
Bak A, Kozik V, Smolinski A, Jampilek J. Multidimensional (3D/4D-QSAR) probability-guided pharmacophore mapping: investigation of activity profile for a series of drug absorption promoters. RSC Adv 2016. [DOI: 10.1039/c6ra15820j] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
A hybrid approach that combines 3D and 4D-QSAR methods based on grid and neural paradigms with automated IVE-PLS procedure was examined to identify the pharmacophore pattern for cholic acid derivatives as potential drug absorption promoters.
Collapse
Affiliation(s)
- A. Bak
- Department of Organic Chemistry
- Institute of Chemistry
- University of Silesia
- Katowice
- Poland
| | - V. Kozik
- Department of Synthesis Chemistry
- Institute of Chemistry
- University of Silesia
- Katowice
- Poland
| | - A. Smolinski
- Department of Energy Saving and Air Protection
- Central Mining Institute
- Katowice
- Poland
| | - J. Jampilek
- Department of Pharmaceutical Chemistry
- Faculty of Pharmacy
- Comenius University
- Bratislava
- Slovakia
| |
Collapse
|
8
|
de Assis TM, Gajo GC, de Assis LC, Garcia LS, Silva DR, Ramalho TC, da Cunha EFF. QSAR Models Guided by Molecular Dynamics Applied to Human Glucokinase Activators. Chem Biol Drug Des 2015; 87:455-66. [DOI: 10.1111/cbdd.12683] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2015] [Revised: 10/05/2015] [Accepted: 10/22/2015] [Indexed: 11/28/2022]
Affiliation(s)
- Tamiris Maria de Assis
- Department of Chemistry; Federal University of Lavras; P.O. Box 3037 Lavras 37200-000 Brazil
| | - Giovanna Cardoso Gajo
- Department of Chemistry; Federal University of Lavras; P.O. Box 3037 Lavras 37200-000 Brazil
| | | | - Letícia Santos Garcia
- Department of Chemistry; Federal University of Lavras; P.O. Box 3037 Lavras 37200-000 Brazil
| | - Daniela Rodrigues Silva
- Department of Chemistry; Federal University of Lavras; P.O. Box 3037 Lavras 37200-000 Brazil
| | - Teodorico Castro Ramalho
- Department of Chemistry; Federal University of Lavras; P.O. Box 3037 Lavras 37200-000 Brazil
- Center for Basic and Applied Research; Faculty of Informatics and Management; University of Hradec Kralove; Hradec Kralove Czech Republic
| | | |
Collapse
|
9
|
Abbat S, Jain V, Bharatam PV. Origins of the specificity of inhibitor P218 toward wild-type and mutantPfDHFR: a molecular dynamics analysis. J Biomol Struct Dyn 2014; 33:1913-28. [DOI: 10.1080/07391102.2014.979231] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
10
|
Varnek A, Baskin I. Machine learning methods for property prediction in chemoinformatics: Quo Vadis? J Chem Inf Model 2012; 52:1413-37. [PMID: 22582859 DOI: 10.1021/ci200409x] [Citation(s) in RCA: 145] [Impact Index Per Article: 12.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
This paper is focused on modern approaches to machine learning, most of which are as yet used infrequently or not at all in chemoinformatics. Machine learning methods are characterized in terms of the "modes of statistical inference" and "modeling levels" nomenclature and by considering different facets of the modeling with respect to input/ouput matching, data types, models duality, and models inference. Particular attention is paid to new approaches and concepts that may provide efficient solutions of common problems in chemoinformatics: improvement of predictive performance of structure-property (activity) models, generation of structures possessing desirable properties, model applicability domain, modeling of properties with functional endpoints (e.g., phase diagrams and dose-response curves), and accounting for multiple molecular species (e.g., conformers or tautomers).
Collapse
Affiliation(s)
- Alexandre Varnek
- Laboratoire d'Infochimie, UMR 7177 CNRS, Université de Strasbourg, 4, rue B. Pascal, Strasbourg 67000, France.
| | | |
Collapse
|
11
|
The great descriptor melting pot: mixing descriptors for the common good of QSAR models. J Comput Aided Mol Des 2011; 26:39-43. [DOI: 10.1007/s10822-011-9511-4] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2011] [Accepted: 12/02/2011] [Indexed: 10/14/2022]
|
12
|
Adane L, Bharatam PV. Binding modes of 2,4-diaminoquinazoline and 2,4-diaminopteridine analogs to P. falciparum dihydrofolate reductase enzyme: Molecular docking studies. Indian J Pharm Sci 2011; 72:324-33. [PMID: 21188041 PMCID: PMC3003165 DOI: 10.4103/0250-474x.70478] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2009] [Revised: 02/03/2010] [Accepted: 04/20/2010] [Indexed: 11/30/2022] Open
Abstract
A molecular docking study was carried out on 28 compounds belonging to 2,4-diaminoquinazoline and 2,4-diaminopteridine analogs using Glide, FlexX and GOLD programs and the X-ray crystallographic structures of the quadruple mutant (1J3K:pdb) and wild type (1J3I:pdb) Plasmodium falciparum dihydrofolate reductase enzyme. The experimental conformation the bound ligand WR99210 was precisely reproduced by the docking procedures as demonstrated by low (<2.00 Å) root-mean-square deviations. The results indicated that most of the compounds dock into the active sites of both the wild type and quadruple mutant P. falciparum dihydrofolate reductase enzymes. Visual inspection of the binding modes also demonstrated that most of the compounds could form H-bond interactions with the key amino acid residues (Asp54, Ile14 and Leu/Ile164) and with better docking scores than the bound compound (5). Their long side chains orient in the hydrophobic portion of the active site which is occupied by trichloro aryloxy side chain of WR99210 (5). Thus, avoid potential steric clashes with Asn108 (mutated from Ser108). Such a clash is known to be responsible for the resistance of the P. falciparum to pyrimethamine and cycloguanil.
Collapse
Affiliation(s)
- L Adane
- Department of Medicinal Chemistry, National Institute of Pharmaceutical Education and Research (NIPER), S.A.S. Nagar-160 062, India
| | | |
Collapse
|
13
|
Basak SC, Mills D. Quantitative structure-activity relationships for cycloguanil analogs as PfDHFR inhibitors using mathematical molecular descriptors. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2010; 21:215-229. [PMID: 20544548 DOI: 10.1080/10629361003770951] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2023]
Abstract
Computed molecular descriptors were used to develop quantitative structure-activity relationships (QSARs) for binding affinities (K(i)) for a set of 58 cycloguanil (2,4-diamino-1,6-dihydro-1,3,5-triazine) analogues for dihydrofolate reductase (DHFR) enzyme extracted from wild and A16V+S108T mutant type (a double mutation) malaria parasite Plasmodium falciparum (Pf). High-quality models were obtained in both cases. The results of statistical analyses show that ridge regression (RR) outperformed the two other modelling methods, principal component regression (PCR) and partial least squares (PLS). For both enzymes, recognition of the inhibitors was based on four broad categories of descriptors encoding information on: (1) the electronic character of the various atoms in the molecule, (2) the size and shape of the structure, (3) the degree of branching in the molecular skeleton, and (4) two to five atom molecular fragments with aliphatic carbon at one end and aliphatic or aromatic carbon or nitrogen at the other end. The subsets of influential descriptors underlying the QSARs for the wild versus the mutant DHFR are quite non-overlapping. This indicates that the two enzymes recognize the inhibitor molecules on the basis of mutually distinct structural attributes. Such differential QSARs can be useful in the design of novel drugs active against malaria parasites which are growing in resistant to existing chemotherapeutic agents.
Collapse
Affiliation(s)
- S C Basak
- Center for Water and the Environment, Natural Resources Research Institute, University of Minnesota Duluth, Duluth, USA.
| | | |
Collapse
|
14
|
Adane L, Patel DS, Bharatam PV. Shape- and Chemical Feature-Based 3D-Pharmacophore Model Generation and Virtual Screening: Identification of Potential Leads forP. falciparumDHFR Enzyme Inhibition. Chem Biol Drug Des 2010; 75:115-26. [DOI: 10.1111/j.1747-0285.2009.00908.x] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
15
|
Adane L, Bharatam PV. 3D-QSAR analysis of cycloguanil derivatives as inhibitors of A16V+S108T mutant Plasmodium falciparum dihydrofolate reductase enzyme. J Mol Graph Model 2009; 28:357-67. [DOI: 10.1016/j.jmgm.2009.09.001] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2009] [Revised: 08/27/2009] [Accepted: 09/01/2009] [Indexed: 12/17/2022]
|
16
|
Sivaprakasam P, Tosso PN, Doerksen RJ. Structure-activity relationship and comparative docking studies for cycloguanil analogs as PfDHFR-TS inhibitors. J Chem Inf Model 2009; 49:1787-96. [PMID: 19588935 DOI: 10.1021/ci9000663] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Drug resistance acquired by Plasmodium falciparum (Pf) is a major problem in the treatment and control of malaria. One of the major examples of drug resistance is that caused by mutations in the active site of dihydrofolate reductase (DHFR) of Pf (PfDHFR-TS). A double mutation, A16V+S108T, is specific for resistance to the marketed drug cycloguanil. In this study, we used 58 cycloguanil (2,4-diamino-1,6-dihydro-1,3,5-triazine) derivatives to explore the relationship between various physicochemical properties and reported binding affinity data on wild-type and mutant-type A16V+S108T. Using the Hansch 2D-quantitative structure-activity relationship method, we obtained a parabolic relationship of hydrophobicity of substituents at the N1-phenyl ring with the wild-type binding affinity data. Hydrophobicity being a key property for wild-type binding affinity data, we found steric factors to be crucial for A16V+S108T mutant resistance. We investigated FlexX, GOLD, Glide and Molegro virtual docking programs and 13 different scoring functions on 10 of the cycloguanil derivatives to evaluate which program was best for reproducing the experimental binding mode and correlating the docking scores with the reported binding affinity data. We identified GOLD, using its GoldScore fitness function, as the most accurate docking program for predicting binding affinity data of cycloguanil derivatives to DHFR and Molegro virtual docker, with its template docking algorithm and MolDock [GRID] scoring function, as most accurate for reproducing the experimental binding mode of a reference ligand that is structurally similar to the cycloguanil derivatives studied. We also report an interaction index which best describes the structure-activity relationships exhibited by these analogs in terms of PfDHFR-TS active site interactions.
Collapse
Affiliation(s)
- Prasanna Sivaprakasam
- Department of Medicinal Chemistry, School of Pharmacy, University of Mississippi, 417 Faser Hall, University, Mississippi 38677-1848, USA
| | | | | |
Collapse
|
17
|
Hecht D, Fogel GB. A Novel In Silico Approach to Drug Discovery via Computational Intelligence. J Chem Inf Model 2009; 49:1105-21. [DOI: 10.1021/ci9000647] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- David Hecht
- Southwestern College, 900 Otay Lakes Road, Chula Vista, California 91910, and Natural Selection, Inc., 9330 Scranton Road, Suite 150, San Diego, California 92121
| | - Gary B. Fogel
- Southwestern College, 900 Otay Lakes Road, Chula Vista, California 91910, and Natural Selection, Inc., 9330 Scranton Road, Suite 150, San Diego, California 92121
| |
Collapse
|
18
|
Romeiro NC, Albuquerque MG, de Alencastro RB, Ravi M, Hopfinger AJ. Construction of 4D-QSAR models for use in the design of novel p38-MAPK inhibitors. J Comput Aided Mol Des 2008; 19:385-400. [PMID: 16231199 DOI: 10.1007/s10822-005-7927-4] [Citation(s) in RCA: 26] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2005] [Accepted: 05/22/2005] [Indexed: 12/11/2022]
Abstract
The p38-mitogen-activated protein kinase (p38-MAPK) plays a key role in lipopolysaccharide-induced tumor necrosis factor-alpha (TNF-alpha) and interleukin-1 (IL-1) release during the inflammatory process, emerging as an attractive target for new anti-inflammatory agents. Four-dimensional quantitative structure-activity relationship (4D-QSAR) analysis [Hopfinger et al., J. Am. Chem. Soc., 119 (1997) 10509] was applied to a series of 33 (a training set of 28 and a test set of 5) pyridinyl-imidazole and pyrimidinyl-imidazole inhibitors of p38-MAPK, with IC50 ranging from 0.11 to 2100 nM [Liverton et al., J. Med. Chem., 42 (1999) 2180]. Five thousand conformations of each analogue were sampled from a molecular dynamics simulation (MDS) during 50 ps at a constant temperature of 303 K. Each conformation was placed in a 2 angstroms grid cell lattice for each of three trial alignments. 4D-QSAR models were constructed by genetic algorithm (GA) optimization and partial least squares (PLS) fitting, and evaluated by leave-one-out cross-validation technique. In the best models, with three to six terms, the adjusted cross-validated squared correlation coefficients, Q2adj, ranged from 0.67 to 0.85. Model D (Q2adj = 0.84) was identified as the most robust model from alignment 1, and it is representative of the other best models. This model encompasses new molecular regions as containing pharmacophore sites, such as the amino-benzyl moiety of pyrimidine analogs and the N1-substituent in the imidazole ring. These regions of the ligands should be further explored to identify better anti-inflammatory inhibitors of p38-MAPK.
Collapse
Affiliation(s)
- Nelilma Correia Romeiro
- Departamento de Química Orgânica, Laboratório de Modelagem Molecular (LabMMol), Instituto de Química, Centro de Tecnologia, Bloco A, Ilha do Fundão, Universidade Federal do Rio de Janeiro, 21949-900, Rio de Janeiro, RJ, Brasil.
| | | | | | | | | |
Collapse
|
19
|
Fogel GB, Cheung M, Pittman E, Hecht D. In silico screening against wild-type and mutant Plasmodium falciparum dihydrofolate reductase. J Mol Graph Model 2008; 26:1145-52. [DOI: 10.1016/j.jmgm.2007.10.006] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2007] [Revised: 10/11/2007] [Accepted: 10/11/2007] [Indexed: 12/21/2022]
|
20
|
Hecht D, Cheung M, Fogel GB. QSAR using evolved neural networks for the inhibition of mutant PfDHFR by pyrimethamine derivatives. Biosystems 2008; 92:10-5. [DOI: 10.1016/j.biosystems.2007.10.005] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2007] [Revised: 10/25/2007] [Accepted: 10/29/2007] [Indexed: 10/22/2022]
|
21
|
A QSAR study and molecular design of benzothiazole derivatives as potent anticancer agents. ACTA ACUST UNITED AC 2008. [DOI: 10.1007/s11426-007-0107-8] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
22
|
Modeling the inhibition of quadruple mutant Plasmodium falciparum dihydrofolate reductase by pyrimethamine derivatives. J Comput Aided Mol Des 2007; 22:29-38. [DOI: 10.1007/s10822-007-9152-9] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2007] [Accepted: 11/15/2007] [Indexed: 11/27/2022]
|
23
|
Hilmy KMH. Synthesis of Non-Nucleosides: 7- and 1,3-Substituents of New Pyrrolo[2,3-d]pyrimidin-4-ones on Antiviral Activity. Arch Pharm (Weinheim) 2006; 339:174-81. [PMID: 16586425 DOI: 10.1002/ardp.200500162] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
A series of non-nucleosides 9-47 were synthesized. Compounds 1-4 were reacted with formic acid (85%) to afford compounds 5-8. Then, the latter compounds were reacted with alkyl halides a-f (2-bromopropane, 2-bromobutane, benzyl bromide, benzyl chloromethyl ether, chloromethyl ethyl ether, phenacyl bromide) in the presence of NaH in dry DMF to give the desired compounds 9-47, which were evaluated for activity against herpes simplex virus type-II (HSV-II).
Collapse
|
24
|
Affiliation(s)
- Ivan M Kompis
- ARPIDA Ltd, Dammstrasse 36, 4142 Münchenstein, Switzerland
| | | | | |
Collapse
|
25
|
Marrero-Ponce Y, Montero-Torres A, Zaldivar CR, Veitía MI, Peréz MM, Sánchez RNG. Non-stochastic and stochastic linear indices of the ‘molecular pseudograph’s atom adjacency matrix’: application to ‘in silico’ studies for the rational discovery of new antimalarial compounds. Bioorg Med Chem 2005; 13:1293-304. [PMID: 15670938 DOI: 10.1016/j.bmc.2004.11.008] [Citation(s) in RCA: 52] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2004] [Revised: 11/05/2004] [Accepted: 11/05/2004] [Indexed: 11/18/2022]
Abstract
Malaria is one of the most deadly diseases, affecting million of people especially in developing countries. Because of the rapidly increasing threat worldwide of malaria epidemics multidrugs resistant to therapies, there is an urgent global need to discover new classes of antimalarial compounds. In an effort to overcome this problem, we have investigated the use of structure-based classification models for the 'rational' selection/identification or design/optimization of new lead antimalarials from virtual combinatorial data sets. In this sense, TOpological MOlecular COMputer Design strategy (TOMOCOMD approach) has been introduced in order to obtain two quantitative models for the discrimination of antimalarials. A collected data set containing 597 antimalarial compounds is presented as a helpful tool not only for theoretical chemist but for other researchers in this area. The validated models (including non-stochastic and stochastic indices) classify correctly more than 90% of compounds in both training and external prediction data sets. They showed high Matthews' correlation coefficients; 0.87 and 0.82 for training and 0.86 and 0.79 for test set. The TOMOCOMD-CARDD approach implemented in this work was successfully compared with two of the most useful models for antimalarials selection reported so far. Thus we expect that these two QSAR models can be used in the identification of previously un-known antimalarials compounds.
Collapse
Affiliation(s)
- Yovani Marrero-Ponce
- Department of Pharmacy, Faculty of Chemical-Pharmacy. Central University of Las Villas, Santa Clara, 54830 Villa Clara, Cuba
| | | | | | | | | | | |
Collapse
|
26
|
Garcia-Viloca M, Truhlar DG, Gao J. Importance of substrate and cofactor polarization in the active site of dihydrofolate reductase. J Mol Biol 2003; 327:549-60. [PMID: 12628257 DOI: 10.1016/s0022-2836(03)00123-2] [Citation(s) in RCA: 73] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
By using a combined quantum-mechanical and molecular-mechanical potential in molecular dynamics simulations, we have investigated the effects of the enzyme electric field of dihydrofolate reductase on the electronic polarization of its 5-protonated dihydrofolate substrate at various stages of the catalyzed hydride transfer reaction. Energy decomposition of the total electrostatic interaction energy between the ligands and the enzyme shows that the polarization effect is 4% of the total electrostatic interaction energy, and, significantly, it accounts for 9kcal/mol of transition state stabilization relative to the reactant state. Therefore it is essential to take account of substrate polarization for quantitative interpretation of enzymatic function and for calculation of binding free energies of inhibitors to a protein. Atomic polarizations are calculated as the differences in the average atomic charges on the atoms in gas phase and in molecular simulations of the enzyme; this analysis shows that the glutamate tail and the pterin ring are the highly polarized regions of the substrate. Electron density difference plots of the reactant and product complexes at instantaneous configurations in the enzyme active center confirm the inferences made on the basis of partial atomic charges.
Collapse
Affiliation(s)
- Mireia Garcia-Viloca
- Department of Chemistry and Minnesota Supercomputing Institute, University of Minnesota, 207 Pleasant Street, SE Minneapolis 55455-0431, USA.
| | | | | |
Collapse
|
27
|
Sirawaraporn W, Sirawaraporn R, Yongkiettrakul S, Anuwatwora A, Rastelli G, Kamchonwongpaisan S, Yuthavong Y. Mutational analysis of Plasmodium falciparum dihydrofolate reductase: the role of aspartate 54 and phenylalanine 223 on catalytic activity and antifolate binding. Mol Biochem Parasitol 2002; 121:185-93. [PMID: 12034452 DOI: 10.1016/s0166-6851(02)00035-x] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
Abstract
The catalytic activity and ability to confer resistance to antifolates of Plasmodium falciparum dihydrofolate reductase (pfDHFR) through single and double mutations at Asp-54 and Phe-223 were investigated. A single Asp54Glu (D54E) mutation in the pfDHFR domain greatly decreased the catalytic activity of the enzyme and affected both the K(m) values for the substrate dihydrofolate and the K(i) values for pyrimethamine, cycloguanil and WR99210. The Phe223Ser (F223S) single mutant had unperturbed kinetics but had very poor affinity with the first two antifolates. The ability to confer high resistance to the antifolates of F223S enzyme was, however, abolished in the D54E+F223S double mutant enzyme. When D54E mutation was present together with the A16V+S108T double mutation, the effects on the K(m) values for the substrate dihydrofolate and the binding affinity of antifolates were much more pronounced. The severely impaired kinetics and poor activity observed in A16V+S108T+D54E enzyme could, however, be restored when F223S was introduced, while the binding affinities to the antifolates remained poor. The experimental findings can be explained with a model for substrate and inhibitor binding. Our data not only indicate the importance of Asp-54 of pfDHFR in catalysis and inhibitor binding, but also provide evidence that infer the potentially crucial function of the C-terminal portion of pfDHFR domain.
Collapse
Affiliation(s)
- Worachart Sirawaraporn
- Department of Biochemistry, Faculty of Science, Mahidol University, Rama 6 Rd., 10400, Bangkok 10400, Thailand.
| | | | | | | | | | | | | |
Collapse
|
28
|
Santos-Filho OA, Mishra RK, Hopfinger AJ. Free energy force field (FEFF) 3D-QSAR analysis of a set of Plasmodium falciparum dihydrofolate reductase inhibitors. J Comput Aided Mol Des 2001; 15:787-810. [PMID: 11776291 DOI: 10.1023/a:1013199108020] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Free energy force field (FEFF) 3D-QSAR analysis was used to construct ligand-receptor binding models for a set of 18 structurally diverse antifolates including pyrimethamine, cycloguanil, methotrexate, aminopterin and trimethoprim, and 13 pyrrolo[2,3-d]pyrimidines. The molecular target ('receptor') used was a 3D-homology model of a specific mutant type of Plasmodium falciparum (Pf) dihydrofolate reductase (DHFR). The dependent variable of the 3D-QSAR models is the IC50 inhibition constant for the specific mutant type of PfDHFR. The independent variables of the 3D-QSAR models (the descriptors) are scaled energy terms of a modified first-generation AMBER force field combined with a hydration shell aqueous solvation model and a collection of 2D-QSAR descriptors often used in QSAR studies. Multiple temperature molecular dynamics simulation (MDS) and the genetic function approximation (GFA) were employed using partial least square (PLS) and multidimensional linear regressions as the fitting functions to develop FEFF 3D-QSAR models for the binding process. The significant FEFF energy terms in the best 3D-QSAR models include energy contributions of the direct ligand-receptor interaction. Some changes in conformational energy terms of the ligand due to binding to the enzyme are also found to be important descriptors. The FEFF 3D-QSAR models indicate some structural features perhaps relevant to the mechanism of resistance of the PfDHFR to current antimalarials. The FEFF 3D-QSAR models are also compared to receptor-independent (RI) 4D-QSAR models developed in an earlier study and subsequently refined using recently developed generalized alignment rules.
Collapse
Affiliation(s)
- O A Santos-Filho
- Laboratory of Molecular Modeling and Design, University of Illinois at Chicago, College of Pharmacy, 60612-7231, USA
| | | | | |
Collapse
|