1
|
Bhatnagar A, Nath V, Kumar N, Kumar V. Discovery of novel PARP-1 inhibitors using tandem in silico studies: integrated docking, e-pharmacophore, deep learning based de novo and molecular dynamics simulation approach. J Biomol Struct Dyn 2024; 42:3396-3409. [PMID: 37216358 DOI: 10.1080/07391102.2023.2214223] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2023] [Accepted: 05/05/2023] [Indexed: 05/24/2023]
Abstract
Cancer accounts for the majority of deaths worldwide, and the increasing incidence of breast cancer is a matter of grave concern. Poly (ADP-ribose) polymerase-1 (PARP-1) has emerged as an attractive target for the treatment of breast cancer as it has an important role in DNA repair. The focus of the study was to identify novel PARP-1 inhibitors using a blend of tandem structure-based screening (Docking and e-pharmacophore-based screening) and artificial intelligence (deep learning)-based de novo approaches. The scrutiny of compounds having good binding characteristics for PARP-1 was carried out using a tandem mode of screening along with parameters such as binding energy and ADME analysis. The efforts afforded compound Vab1 (PubChem ID 129142036), which was chosen as a seed for obtaining novel compounds through a trained artificial intelligence (AI)-based model. Resultant compounds were assessed for PARP-1 inhibition; binding affinity prediction and interaction pattern analysis were carried out using the extra precision (XP) mode of docking. Two best hits, Vab1-b and Vab1-g, exhibiting good dock scores and suitable interactions, were subjected to 100 nanoseconds (ns) of molecular dynamics simulation in the active site of PARP-1 and compared with the reference Protein-Ligand Complex. The stable nature of PARP-1 upon binding to these compounds was revealed through MD simulation.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Aayushi Bhatnagar
- Department of Pharmacy, School of Chemical Sciences and Pharmacy, Central University of Rajasthan, Ajmer, India
| | - Virendra Nath
- Department of Pharmacy, School of Chemical Sciences and Pharmacy, Central University of Rajasthan, Ajmer, India
| | - Neeraj Kumar
- Bhupal Nobles' College of Pharmacy, Bhupal Nobles' University, Udaipur, India
| | - Vipin Kumar
- Department of Pharmacy, School of Chemical Sciences and Pharmacy, Central University of Rajasthan, Ajmer, India
| |
Collapse
|
2
|
Morales JF, Chuguransky S, Alberca LN, Alice JI, Goicoechea S, Ruiz ME, Bellera CL, Talevi A. Positive Predictive Value Surfaces as a Complementary Tool to Assess the Performance of Virtual Screening Methods. Mini Rev Med Chem 2021; 20:1447-1460. [PMID: 32072906 DOI: 10.2174/1871525718666200219130229] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2019] [Revised: 10/28/2019] [Accepted: 10/29/2019] [Indexed: 11/22/2022]
Abstract
BACKGROUND Since their introduction in the virtual screening field, Receiver Operating Characteristic (ROC) curve-derived metrics have been widely used for benchmarking of computational methods and algorithms intended for virtual screening applications. Whereas in classification problems, the ratio between sensitivity and specificity for a given score value is very informative, a practical concern in virtual screening campaigns is to predict the actual probability that a predicted hit will prove truly active when submitted to experimental testing (in other words, the Positive Predictive Value - PPV). Estimation of such probability is however, obstructed due to its dependency on the yield of actives of the screened library, which cannot be known a priori. OBJECTIVE To explore the use of PPV surfaces derived from simulated ranking experiments (retrospective virtual screening) as a complementary tool to ROC curves, for both benchmarking and optimization of score cutoff values. METHODS The utility of the proposed approach is assessed in retrospective virtual screening experiments with four datasets used to infer QSAR classifiers: inhibitors of Trypanosoma cruzi trypanothione synthetase; inhibitors of Trypanosoma brucei N-myristoyltransferase; inhibitors of GABA transaminase and anticonvulsant activity in the 6 Hz seizure model. RESULTS Besides illustrating the utility of PPV surfaces to compare the performance of machine learning models for virtual screening applications and to select an adequate score threshold, our results also suggest that ensemble learning provides models with better predictivity and more robust behavior. CONCLUSION PPV surfaces are valuable tools to assess virtual screening tools and choose score thresholds to be applied in prospective in silico screens. Ensemble learning approaches seem to consistently lead to improved predictivity and robustness.
Collapse
Affiliation(s)
- Juan F Morales
- Laboratory of Bioactive Research and Development (LIDeB), Department of Biological Sciences, Faculty of Exact Sciences, University of La Plata (UNLP) - 47 & 115, La Plata (1900), Buenos Aires, Argentina
| | - Sara Chuguransky
- Laboratory of Bioactive Research and Development (LIDeB), Department of Biological Sciences, Faculty of Exact Sciences, University of La Plata (UNLP) - 47 & 115, La Plata (1900), Buenos Aires, Argentina
| | - Lucas N Alberca
- Laboratory of Bioactive Research and Development (LIDeB), Department of Biological Sciences, Faculty of Exact Sciences, University of La Plata (UNLP) - 47 & 115, La Plata (1900), Buenos Aires, Argentina
| | - Juan I Alice
- Laboratory of Bioactive Research and Development (LIDeB), Department of Biological Sciences, Faculty of Exact Sciences, University of La Plata (UNLP) - 47 & 115, La Plata (1900), Buenos Aires, Argentina
| | - Sofía Goicoechea
- Laboratory of Bioactive Research and Development (LIDeB), Department of Biological Sciences, Faculty of Exact Sciences, University of La Plata (UNLP) - 47 & 115, La Plata (1900), Buenos Aires, Argentina
| | - María E Ruiz
- Laboratory of Bioactive Research and Development (LIDeB), Department of Biological Sciences, Faculty of Exact Sciences, University of La Plata (UNLP) - 47 & 115, La Plata (1900), Buenos Aires, Argentina
| | - Carolina L Bellera
- Laboratory of Bioactive Research and Development (LIDeB), Department of Biological Sciences, Faculty of Exact Sciences, University of La Plata (UNLP) - 47 & 115, La Plata (1900), Buenos Aires, Argentina
| | - Alan Talevi
- Laboratory of Bioactive Research and Development (LIDeB), Department of Biological Sciences, Faculty of Exact Sciences, University of La Plata (UNLP) - 47 & 115, La Plata (1900), Buenos Aires, Argentina
| |
Collapse
|
3
|
Selvaraj C, Panwar U, Dinesh DC, Boura E, Singh P, Dubey VK, Singh SK. Microsecond MD Simulation and Multiple-Conformation Virtual Screening to Identify Potential Anti-COVID-19 Inhibitors Against SARS-CoV-2 Main Protease. Front Chem 2021; 8:595273. [PMID: 33585398 PMCID: PMC7873971 DOI: 10.3389/fchem.2020.595273] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2020] [Accepted: 11/19/2020] [Indexed: 12/14/2022] Open
Abstract
The recent pandemic outbreak of COVID-19, caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), raised global health and economic concerns. Phylogenetically, SARS-CoV-2 is closely related to SARS-CoV, and both encode the enzyme main protease (Mpro/3CLpro), which can be a potential target inhibiting viral replication. Through this work, we have compiled the structural aspects of Mpro conformational changes, with molecular modeling and 1-μs MD simulations. Long-scale MD simulation resolves the mechanism role of crucial amino acids involved in protein stability, followed by ensemble docking which provides potential compounds from the Traditional Chinese Medicine (TCM) database. These lead compounds directly interact with active site residues (His41, Gly143, and Cys145) of Mpro, which plays a crucial role in the enzymatic activity. Through the binding mode analysis in the S1, S1′, S2, and S4 binding subsites, screened compounds may be functional for the distortion of the oxyanion hole in the reaction mechanism, and it may lead to the inhibition of Mpro in SARS-CoV-2. The hit compounds are naturally occurring compounds; they provide a sustainable and readily available option for medical treatment in humans infected by SARS-CoV-2. Henceforth, extensive analysis through molecular modeling approaches explained that the proposed molecules might be promising SARS-CoV-2 inhibitors for the inhibition of COVID-19, subjected to experimental validation.
Collapse
Affiliation(s)
- Chandrabose Selvaraj
- Computer Aided Drug Design and Molecular Modeling Lab, Department of Bioinformatics, Alagappa University, Karaikudi, India
| | - Umesh Panwar
- Computer Aided Drug Design and Molecular Modeling Lab, Department of Bioinformatics, Alagappa University, Karaikudi, India
| | - Dhurvas Chandrasekaran Dinesh
- Section of Molecular Biology and Biochemistry, Institute of Organic Chemistry and Biochemistry AS CR, v.v.i., Prague, Czechia
| | - Evzen Boura
- Section of Molecular Biology and Biochemistry, Institute of Organic Chemistry and Biochemistry AS CR, v.v.i., Prague, Czechia
| | - Poonam Singh
- Corrosion and Materials Protection Division, Council of Scientific and Industrial Research (CSIR)-Central Electrochemical Research Institute, Karaikudi, India
| | - Vikash Kumar Dubey
- School of Biochemical Engineering, Indian Institute of Technology (BHU), Varanasi, India
| | - Sanjeev Kumar Singh
- Computer Aided Drug Design and Molecular Modeling Lab, Department of Bioinformatics, Alagappa University, Karaikudi, India
| |
Collapse
|
4
|
Narayanan D, Gani OABSM, Gruber FXE, Engh RA. Data driven polypharmacological drug design for lung cancer: analyses for targeting ALK, MET, and EGFR. J Cheminform 2017; 9:43. [PMID: 29086093 PMCID: PMC5496928 DOI: 10.1186/s13321-017-0229-8] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2016] [Accepted: 06/18/2017] [Indexed: 12/14/2022] Open
Abstract
Drug design of protein kinase inhibitors is now greatly enabled by thousands of publicly available X-ray structures, extensive ligand binding data, and optimized scaffolds coming off patent. The extensive data begin to enable design against a spectrum of targets (polypharmacology); however, the data also reveal heterogeneities of structure, subtleties of chemical interactions, and apparent inconsistencies between diverse data types. As a result, incorporation of all relevant data requires expert choices to combine computational and informatics methods, along with human insight. Here we consider polypharmacological targeting of protein kinases ALK, MET, and EGFR (and its drug resistant mutant T790M) in non small cell lung cancer as an example. Both EGFR and ALK represent sources of primary oncogenic lesions, while drug resistance arises from MET amplification and EGFR mutation. A drug which inhibits these targets will expand relevant patient populations and forestall drug resistance. Crizotinib co-targets ALK and MET. Analysis of the crystal structures reveals few shared interaction types, highlighting proton-arene and key CH–O hydrogen bonding interactions. These are not typically encoded into molecular mechanics force fields. Cheminformatics analyses of binding data show EGFR to be dissimilar to ALK and MET, but its structure shows how it may be co-targeted with the addition of a covalent trap. This suggests a strategy for the design of a focussed chemical library based on a pan-kinome scaffold. Tests of model compounds show these to be compatible with the goal of ALK, MET, and EGFR polypharmacology.
Collapse
Affiliation(s)
- Dilip Narayanan
- The Norwegian Structural Biology Center, Department of Chemistry, Faculty of Science, UiT The Arctic University of Norway, Tromsø, Norway
| | - Osman A B S M Gani
- The Norwegian Structural Biology Center, Department of Chemistry, Faculty of Science, UiT The Arctic University of Norway, Tromsø, Norway
| | - Franz X E Gruber
- The Norwegian Structural Biology Center, Department of Chemistry, Faculty of Science, UiT The Arctic University of Norway, Tromsø, Norway
| | - Richard A Engh
- The Norwegian Structural Biology Center, Department of Chemistry, Faculty of Science, UiT The Arctic University of Norway, Tromsø, Norway.
| |
Collapse
|
5
|
Assessing protein kinase target similarity: Comparing sequence, structure, and cheminformatics approaches. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2015; 1854:1605-16. [PMID: 26001898 DOI: 10.1016/j.bbapap.2015.05.004] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/15/2015] [Revised: 05/08/2015] [Accepted: 05/11/2015] [Indexed: 11/22/2022]
Abstract
In just over two decades, structure based protein kinase inhibitor discovery has grown from trial and error approaches, using individual target structures, to structure and data driven approaches that may aim to optimize inhibition properties across several targets. This is increasingly enabled by the growing availability of potent compounds and kinome-wide binding data. Assessing the prospects for adapting known compounds to new therapeutic uses is thus a key priority for current drug discovery efforts. Tools that can successfully link the diverse information regarding target sequence, structure, and ligand binding properties now accompany a transformation of protein kinase inhibitor research, away from single, block-buster drug models, and toward "personalized medicine" with niche applications and highly specialized research groups. Major hurdles for the transformation to data driven drug discovery include mismatches in data types, and disparities of methods and molecules used; at the core remains the problem that ligand binding energies cannot be predicted precisely from individual structures. However, there is a growing body of experimental data for increasingly successful focussing of efforts: focussed chemical libraries, drug repurposing, polypharmacological design, to name a few. Protein kinase target similarity is easily quantified by sequence, and its relevance to ligand design includes broad classification by key binding sites, evaluation of resistance mutations, and the use of surrogate proteins. Although structural evaluation offers more information, the flexibility of protein kinases, and differences between the crystal and physiological environments may make the use of crystal structures misleading when structures are considered individually. Cheminformatics may enable the "calibration" of sequence and crystal structure information, with statistical methods able to identify key correlates to activity but also here, "the devil is in the details." Examples from specific repurposing and polypharmacology applications illustrate these points. This article is part of a Special Issue entitled: Inhibitors of Protein Kinases.
Collapse
|
6
|
Perspective on computational and structural aspects of kinase discovery from IPK2014. BIOCHIMICA ET BIOPHYSICA ACTA-PROTEINS AND PROTEOMICS 2015; 1854:1595-604. [PMID: 25861861 DOI: 10.1016/j.bbapap.2015.03.014] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/20/2015] [Revised: 03/29/2015] [Accepted: 03/30/2015] [Indexed: 01/16/2023]
Abstract
Recent advances in understanding the activity and selectivity of kinase inhibitors and their relationships to protein structure are presented. Conformational selection in kinases is studied from empirical, data-driven and simulation approaches. Ligand binding and its affinity are, in many cases, determined by the predetermined active and inactive conformation of kinases. Binding affinity and selectivity predictions highlight the current state of the art and advances in computational chemistry as it applies to kinase inhibitor discovery. Kinome wide inhibitor profiling and cell panel profiling lead to a better understanding of selectivity and allow for target validation and patient tailoring hypotheses. This article is part of a Special Issue entitled: Inhibitors of Protein Kinases.
Collapse
|