1
|
Rosado-Quiñones AM, Colón-Lorenzo EE, Pala ZR, Bosch J, Kudyba K, Kudyba H, Leed SE, Roncal N, Baerga-Ortiz A, Roche-Lima A, Gerena Y, Fidock DA, Roth A, Vega-Rodríguez J, Serrano AE. Novel hydrazone compounds with broad-spectrum antiplasmodial activity and synergistic interactions with antimalarial drugs. Antimicrob Agents Chemother 2024:e0164323. [PMID: 38639491 DOI: 10.1128/aac.01643-23] [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: 12/14/2023] [Accepted: 03/20/2024] [Indexed: 04/20/2024] Open
Abstract
The development of novel antiplasmodial compounds with broad-spectrum activity against different stages of Plasmodium parasites is crucial to prevent malaria disease and parasite transmission. This study evaluated the antiplasmodial activity of seven novel hydrazone compounds (referred to as CB compounds: CB-27, CB-41, CB-50, CB-53, CB-58, CB-59, and CB-61) against multiple stages of Plasmodium parasites. All CB compounds inhibited blood stage proliferation of drug-resistant or sensitive strains of Plasmodium falciparum in the low micromolar to nanomolar range. Interestingly, CB-41 exhibited prophylactic activity against hypnozoites and liver schizonts in Plasmodium cynomolgi, a primate model for Plasmodium vivax. Four CB compounds (CB-27, CB-41, CB-53, and CB-61) inhibited P. falciparum oocyst formation in mosquitoes, and five CB compounds (CB-27, CB-41, CB-53, CB-58, and CB-61) hindered the in vitro development of Plasmodium berghei ookinetes. The CB compounds did not inhibit the activation of P. berghei female and male gametocytes in vitro. Isobologram assays demonstrated synergistic interactions between CB-61 and the FDA-approved antimalarial drugs, clindamycin and halofantrine. Testing of six CB compounds showed no inhibition of Plasmodium glutathione S-transferase as a putative target and no cytotoxicity in HepG2 liver cells. CB compounds are promising candidates for further development as antimalarial drugs against multidrug-resistant parasites, which could also prevent malaria transmission.
Collapse
Affiliation(s)
- Angélica M Rosado-Quiñones
- Department of Microbiology and Medical Zoology, University of Puerto Rico School of Medicine, San Juan, Puerto Rico
| | - Emilee E Colón-Lorenzo
- Department of Microbiology and Medical Zoology, University of Puerto Rico School of Medicine, San Juan, Puerto Rico
| | - Zarna Rajeshkumar Pala
- Laboratory of Malaria and Vector Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rockville, Maryland, USA
| | - Jürgen Bosch
- Center for Global Health and Diseases, Case Western Reserve University, Cleveland, Ohio, USA
- InterRayBio, LLC, Cleveland, Ohio, USA
| | - Karl Kudyba
- Department of Drug Discovery, Experimental Therapeutics Branch, Walter Reed Army Institute of Research, Silver Spring, Maryland, USA
| | - Heather Kudyba
- Laboratory of Malaria and Vector Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rockville, Maryland, USA
| | - Susan E Leed
- Department of Drug Discovery, Experimental Therapeutics Branch, Walter Reed Army Institute of Research, Silver Spring, Maryland, USA
| | - Norma Roncal
- Department of Drug Discovery, Experimental Therapeutics Branch, Walter Reed Army Institute of Research, Silver Spring, Maryland, USA
| | - Abel Baerga-Ortiz
- Department of Biochemistry, University of Puerto Rico School of Medicine, San Juan, Puerto Rico
| | - Abiel Roche-Lima
- RCMI Program, Medical Science Campus, University of Puerto Rico, San Juan, Puerto Rico
| | - Yamil Gerena
- Department of Pharmacology and Toxicology, University of Puerto Rico School of Medicine, San Juan, Puerto Rico
| | - David A Fidock
- Department of Microbiology and Immunology, Columbia University, New York, New York, USA
- Division of Infectious Diseases, Department of Medicine, Center for Malaria Therapeutics and Antimicrobial Resistance, Columbia University Medical Center, New York, New York, USA
| | - Alison Roth
- Department of Drug Discovery, Experimental Therapeutics Branch, Walter Reed Army Institute of Research, Silver Spring, Maryland, USA
| | - Joel Vega-Rodríguez
- Laboratory of Malaria and Vector Research, National Institute of Allergy and Infectious Diseases, National Institutes of Health, Rockville, Maryland, USA
| | - Adelfa E Serrano
- Department of Microbiology and Medical Zoology, University of Puerto Rico School of Medicine, San Juan, Puerto Rico
| |
Collapse
|
2
|
Creanza TM, Delre P, Ancona N, Lentini G, Saviano M, Mangiatordi GF. Structure-Based Prediction of hERG-Related Cardiotoxicity: A Benchmark Study. J Chem Inf Model 2021; 61:4758-4770. [PMID: 34506150 PMCID: PMC9282647 DOI: 10.1021/acs.jcim.1c00744] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
![]()
Drug-induced blockade of the human
ether-à-go-go-related
gene (hERG) channel is today considered the main
cause of cardiotoxicity in postmarketing surveillance. Hence, several
ligand-based approaches were developed in the last years and are currently
employed in the early stages of a drug discovery process for in silico cardiac safety assessment of drug candidates.
Herein, we present the first structure-based classifiers able to discern hERG binders from nonbinders. LASSO regularized support
vector machines were applied to integrate docking scores and protein–ligand
interaction fingerprints. A total of 396 models were trained and validated
based on: (i) high-quality experimental bioactivity information returned
by 8337 curated compounds extracted from ChEMBL (version 25) and (ii)
structural predictor data. Molecular docking simulations were performed
using GLIDE and GOLD software programs and four different hERG structural models, namely, the recently published structures
obtained by cryoelectron microscopy (PDB codes: 5VA1 and 7CN1) and
two published homology models selected for comparison. Interestingly,
some classifiers return performances comparable to ligand-based models
in terms of area under the ROC curve (AUCMAX = 0.86 ±
0.01) and negative predictive values (NPVMAX = 0.81 ±
0.01), thus putting forward the herein proposed computational workflow
as a valuable tool for predicting hERG-related cardiotoxicity
without the limitations of ligand-based models, typically affected
by low interpretability and a limited applicability domain. From a
methodological point of view, our study represents the first example
of a successful integration of docking scores and protein–ligand
interaction fingerprints (IFs) through a support vector machine (SVM)
LASSO regularized strategy. Finally, the study highlights the importance
of using hERG structural models accounting for ligand-induced
fit effects and allowed us to select the best-performing protein conformation
(made available in the Supporting Information, SI) to be employed
for a reliable structure-based prediction of hERG-related cardiotoxicity.
Collapse
Affiliation(s)
- Teresa Maria Creanza
- CNR-Institute of Intelligent Industrial Technologies and Systems for Advanced Manufacturing, Via Amendola 122/o, 70126 Bari, Italy
| | - Pietro Delre
- Chemistry Department, University of Bari "Aldo Moro", via E. Orabona, 4, I-70125 Bari, Italy.,CNR-Institute of Crystallography, Via Amendola 122/o, 70126 Bari, Italy
| | - Nicola Ancona
- CNR-Institute of Intelligent Industrial Technologies and Systems for Advanced Manufacturing, Via Amendola 122/o, 70126 Bari, Italy
| | - Giovanni Lentini
- Department of Pharmacy-Pharmaceutical Sciences, University of Bari "Aldo Moro", via E. Orabona, 4, I-70125 Bari, Italy
| | - Michele Saviano
- CNR-Institute of Crystallography, Via Amendola 122/o, 70126 Bari, Italy
| | | |
Collapse
|
3
|
Vázquez J, López M, Gibert E, Herrero E, Luque FJ. Merging Ligand-Based and Structure-Based Methods in Drug Discovery: An Overview of Combined Virtual Screening Approaches. Molecules 2020; 25:E4723. [PMID: 33076254 PMCID: PMC7587536 DOI: 10.3390/molecules25204723] [Citation(s) in RCA: 77] [Impact Index Per Article: 19.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2020] [Revised: 10/06/2020] [Accepted: 10/11/2020] [Indexed: 12/20/2022] Open
Abstract
Virtual screening (VS) is an outstanding cornerstone in the drug discovery pipeline. A variety of computational approaches, which are generally classified as ligand-based (LB) and structure-based (SB) techniques, exploit key structural and physicochemical properties of ligands and targets to enable the screening of virtual libraries in the search of active compounds. Though LB and SB methods have found widespread application in the discovery of novel drug-like candidates, their complementary natures have stimulated continued efforts toward the development of hybrid strategies that combine LB and SB techniques, integrating them in a holistic computational framework that exploits the available information of both ligand and target to enhance the success of drug discovery projects. In this review, we analyze the main strategies and concepts that have emerged in the last years for defining hybrid LB + SB computational schemes in VS studies. Particularly, attention is focused on the combination of molecular similarity and docking, illustrating them with selected applications taken from the literature.
Collapse
Affiliation(s)
- Javier Vázquez
- Pharmacelera, Plaça Pau Vila, 1, Sector C 2a, Edificio Palau de Mar, 08039 Barcelona, Spain;
- Department of Nutrition, Food Science and Gastronomy, Faculty of Pharmacy and Food Sciences, Institute of Biomedicine (IBUB), and Institute of Theoretical and Computational Chemistry (IQTC-UB), University of Barcelona, Av. Prat de la Riba 171, E-08921 Santa Coloma de Gramanet, Spain
| | - Manel López
- AB Science, Parc Scientifique de Luminy, Zone Luminy Enterprise, Case 922, 163 Av. de Luminy, 13288 Marseille, France;
| | - Enric Gibert
- Pharmacelera, Plaça Pau Vila, 1, Sector C 2a, Edificio Palau de Mar, 08039 Barcelona, Spain;
| | - Enric Herrero
- Pharmacelera, Plaça Pau Vila, 1, Sector C 2a, Edificio Palau de Mar, 08039 Barcelona, Spain;
| | - F. Javier Luque
- Department of Nutrition, Food Science and Gastronomy, Faculty of Pharmacy and Food Sciences, Institute of Biomedicine (IBUB), and Institute of Theoretical and Computational Chemistry (IQTC-UB), University of Barcelona, Av. Prat de la Riba 171, E-08921 Santa Coloma de Gramanet, Spain
| |
Collapse
|
4
|
Colón-Lorenzo EE, Colón-López DD, Vega-Rodríguez J, Dupin A, Fidock DA, Baerga-Ortiz A, Ortiz JG, Bosch J, Serrano AE. Structure-Based Screening of Plasmodium berghei Glutathione S-Transferase Identifies CB-27 as a Novel Antiplasmodial Compound. Front Pharmacol 2020; 11:246. [PMID: 32256353 PMCID: PMC7090221 DOI: 10.3389/fphar.2020.00246] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2019] [Accepted: 02/24/2020] [Indexed: 11/13/2022] Open
Abstract
Plasmodium falciparum parasites are increasingly drug-resistant, requiring the search for novel antimalarials with distinct modes of action. Enzymes in the glutathione pathway, including glutathione S-transferase (GST), show promise as novel antimalarial targets. This study aims to better understand the biological function of Plasmodium GST, assess its potential as a drug target, and identify novel antiplasmodial compounds using the rodent model P. berghei. By using reverse genetics, we provided evidence that GST is essential for survival of P. berghei intra-erythrocytic stages and is a valid target for drug development. A structural model of the P. berghei glutathione S-transferase (PbGST) protein was generated and used in a structure-based screening of 900,000 compounds from the ChemBridge Hit2Lead library. Forty compounds were identified as potential inhibitors and analyzed in parasite in vitro drug susceptibility assays. One compound, CB-27, exhibited antiplasmodial activity with an EC50 of 0.5 μM toward P. berghei and 0.9 μM toward P. falciparum multidrug-resistant Dd2 clone B2 parasites. Moreover, CB-27 showed a concentration-dependent inhibition of the PbGST enzyme without inhibiting the human ortholog. A shape similarity screening using CB-27 as query resulted in the identification of 24 novel chemical scaffolds, with six of them showing antiplasmodial activity ranging from EC50 of 0.6-4.9 μM. Pharmacokinetic and toxicity predictions suggest that the lead compounds have drug-likeness properties. The antiplasmodial potency, the absence of hemolytic activity, and the predicted drug-likeness properties position these compounds for lead optimization and further development as antimalarials.
Collapse
Affiliation(s)
- Emilee E. Colón-Lorenzo
- Department of Microbiology and Medical Zoology, University of Puerto Rico School of Medicine, San Juan, PR, United States
| | - Daisy D. Colón-López
- Department of Biochemistry and Molecular Biology, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, United States
| | - Joel Vega-Rodríguez
- Department of Microbiology and Medical Zoology, University of Puerto Rico School of Medicine, San Juan, PR, United States
| | - Alice Dupin
- Department of Microbiology and Immunology, Columbia University Medical Center, New York, NY, United States
| | - David A. Fidock
- Department of Microbiology and Immunology, Columbia University Medical Center, New York, NY, United States
- Division of Infectious Diseases, Department of Medicine, Columbia University Medical Center, New York, NY, United States
| | - Abel Baerga-Ortiz
- Department of Biochemistry, University of Puerto Rico School of Medicine, San Juan, PR, United States
| | - José G. Ortiz
- Department of Pharmacology and Toxicology, University of Puerto Rico School of Medicine, San Juan, PR, United States
| | - Jürgen Bosch
- Department of Biochemistry and Molecular Biology, Johns Hopkins Bloomberg School of Public Health, Baltimore, MD, United States
- Division of Pediatric Pulmonology and Allergy/Immunology, Case Western Reserve University, Cleveland, OH, United States
- InterRayBio, LLC, Baltimore, MD, United States
| | - Adelfa E. Serrano
- Department of Microbiology and Medical Zoology, University of Puerto Rico School of Medicine, San Juan, PR, United States
| |
Collapse
|
5
|
Morrone JA, Weber JK, Huynh T, Luo H, Cornell WD. Combining Docking Pose Rank and Structure with Deep Learning Improves Protein-Ligand Binding Mode Prediction over a Baseline Docking Approach. J Chem Inf Model 2020; 60:4170-4179. [PMID: 32077698 DOI: 10.1021/acs.jcim.9b00927] [Citation(s) in RCA: 52] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
We present a simple, modular graph-based convolutional neural network that takes structural information from protein-ligand complexes as input to generate models for activity and binding mode prediction. Complex structures are generated by a standard docking procedure and fed into a dual-graph architecture that includes separate subnetworks for the ligand bonded topology and the ligand-protein contact map. Recent work has indicated that data set bias drives many past promising results derived from combining deep learning and docking. Our dual-graph network allows contributions from ligand identity that give rise to such biases to be distinguished from effects of protein-ligand interactions on classification. We show that our neural network is capable of learning from protein structural information when, as in the case of binding mode prediction, an unbiased data set is constructed. We next develop a deep learning model for binding mode prediction that uses docking ranking as input in combination with docking structures. This strategy mirrors past consensus models and outperforms a baseline docking program (AutoDock Vina) in a variety of tests, including on cross-docking data sets that mimic real-world docking use cases. Furthermore, the magnitudes of network predictions serve as reliable measures of model confidence.
Collapse
Affiliation(s)
- Joseph A Morrone
- Healthcare & Life Sciences Research, IBM TJ Watson Research Center, 1101 Kitchawan Road, Yorktown Heights, New York 10598, United States
| | - Jeffrey K Weber
- Healthcare & Life Sciences Research, IBM TJ Watson Research Center, 1101 Kitchawan Road, Yorktown Heights, New York 10598, United States
| | - Tien Huynh
- Healthcare & Life Sciences Research, IBM TJ Watson Research Center, 1101 Kitchawan Road, Yorktown Heights, New York 10598, United States
| | - Heng Luo
- Healthcare & Life Sciences Research, IBM TJ Watson Research Center, 1101 Kitchawan Road, Yorktown Heights, New York 10598, United States
| | - Wendy D Cornell
- Healthcare & Life Sciences Research, IBM TJ Watson Research Center, 1101 Kitchawan Road, Yorktown Heights, New York 10598, United States
| |
Collapse
|
6
|
Molecular Docking: Shifting Paradigms in Drug Discovery. Int J Mol Sci 2019; 20:ijms20184331. [PMID: 31487867 PMCID: PMC6769923 DOI: 10.3390/ijms20184331] [Citation(s) in RCA: 742] [Impact Index Per Article: 148.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2019] [Revised: 09/02/2019] [Accepted: 09/02/2019] [Indexed: 12/11/2022] Open
Abstract
Molecular docking is an established in silico structure-based method widely used in drug discovery. Docking enables the identification of novel compounds of therapeutic interest, predicting ligand-target interactions at a molecular level, or delineating structure-activity relationships (SAR), without knowing a priori the chemical structure of other target modulators. Although it was originally developed to help understanding the mechanisms of molecular recognition between small and large molecules, uses and applications of docking in drug discovery have heavily changed over the last years. In this review, we describe how molecular docking was firstly applied to assist in drug discovery tasks. Then, we illustrate newer and emergent uses and applications of docking, including prediction of adverse effects, polypharmacology, drug repurposing, and target fishing and profiling, discussing also future applications and further potential of this technique when combined with emergent techniques, such as artificial intelligence.
Collapse
|
7
|
Kumar A, Zhang KYJ. Improving ligand 3D shape similarity-based pose prediction with a continuum solvent model. J Comput Aided Mol Des 2019; 33:1045-1055. [PMID: 31463704 DOI: 10.1007/s10822-019-00220-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2019] [Accepted: 08/17/2019] [Indexed: 10/26/2022]
Abstract
In order to improve the pose prediction performance of docking methods, we have previously developed the pose prediction using shape similarity (PoPSS) method. It identifies a ligand conformation of the highest shape similarity with target protein crystal ligands. The identified ligand conformation is then placed into the target protein binding pocket and refined using side-chain repacking and Monte Carlo energy minimization. Subsequently, we have reported a modification to PoPSS, named as PoPSS-Lite, using a simple grid-based energy minimization for side-chain repacking and Tversky correlation coefficient as the similarity metric. This modification has improved the pose prediction performance and PoPSS-Lite was one of the top performers in D3R GC3. Here we report a further modification to PoPSS that utilizes a continuum solvent model to account for water mediated protein ligand interactions. In this approach, named as PoPSS-PB, the ligand conformation of the highest shape similarity with crystal ligands is refined along with the target protein binding site by incorporating the Poisson-Boltzmann electrostatics. The performance of PoPSS-PB along with PoPSS and PoPSS-Lite was prospectively evaluated in D3R GC4. PoPSS-PB not only demonstrated excellent performance with mean and median RMSDs of 1.20 and 1.13 Å but also achieved improved performance over PoPSS and PoPSS-Lite. Furthermore, the comparison with other D3R GC4 pose prediction submissions revealed admirable performance. Our results showed that the binding poses of ligands with unknown binding modes can be successfully predicted by utilizing ligand 3D shape similarity with known crystallographic ligands and that taking the solvation into consideration improves pose prediction.
Collapse
Affiliation(s)
- Ashutosh Kumar
- Laboratory for Structural Bioinformatics, Center for Biosystems Dynamics Research, RIKEN, 1-7-22 Suehiro, Tsurumi, Yokohama, Kanagawa, 230-0045, Japan
| | - Kam Y J Zhang
- Laboratory for Structural Bioinformatics, Center for Biosystems Dynamics Research, RIKEN, 1-7-22 Suehiro, Tsurumi, Yokohama, Kanagawa, 230-0045, Japan.
| |
Collapse
|
8
|
Kurkinen ST, Lätti S, Pentikäinen OT, Postila PA. Getting Docking into Shape Using Negative Image-Based Rescoring. J Chem Inf Model 2019; 59:3584-3599. [PMID: 31290660 PMCID: PMC6750746 DOI: 10.1021/acs.jcim.9b00383] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
The failure of default scoring functions to ensure virtual screening enrichment is a persistent problem for the molecular docking algorithms used in structure-based drug discovery. To remedy this problem, elaborate rescoring and postprocessing schemes have been developed with a varying degree of success, specificity, and cost. The negative image-based rescoring (R-NiB) has been shown to improve the flexible docking performance markedly with a variety of drug targets. The yield improvement is achieved by comparing the alternative docking poses against the negative image of the target protein's ligand-binding cavity. In other words, the shape and electrostatics of the binding pocket is directly used in the similarity comparison to rank the explicit docking poses. Here, the PANTHER/ShaEP-based R-NiB methodology is tested with six popular docking softwares, including GLIDE, PLANTS, GOLD, DOCK, AUTODOCK, and AUTODOCK VINA, using five validated benchmark sets. Overall, the results indicate that R-NiB outperforms the default docking scoring consistently and inexpensively, demonstrating that the methodology is ready for wide-scale virtual screening usage.
Collapse
Affiliation(s)
- Sami T Kurkinen
- Institute of Biomedicine, Kiinamyllynkatu 10, Integrative Physiology and Pharmacy , University of Turku , FI-20520 Turku , Finland
| | - Sakari Lätti
- Institute of Biomedicine, Kiinamyllynkatu 10, Integrative Physiology and Pharmacy , University of Turku , FI-20520 Turku , Finland
| | - Olli T Pentikäinen
- Institute of Biomedicine, Kiinamyllynkatu 10, Integrative Physiology and Pharmacy , University of Turku , FI-20520 Turku , Finland.,Aurlide Ltd. , FI-21420 Lieto , Finland
| | - Pekka A Postila
- Department of Biological and Environmental Science , University of Jyvaskyla , P.O. Box 35, FI-40014 Jyvaskyla , Finland
| |
Collapse
|
9
|
Kumar A, Zhang KYJ. Shape similarity guided pose prediction: lessons from D3R Grand Challenge 3. J Comput Aided Mol Des 2018; 33:47-59. [PMID: 30084081 DOI: 10.1007/s10822-018-0142-x] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2018] [Accepted: 08/01/2018] [Indexed: 12/15/2022]
Abstract
To extend the utility of ligand 3D shape similarity into pose prediction and virtual screening, we have previously developed CDVS and PoPSS methods. Both of them utilize ligand 3D shape similarity with the crystallographic ligands to improve pose prediction. While CDVS utilizes shape similarity to select suitable receptor structures for molecular docking, PoPSS places a ligand conformation of the highest shape similarity with crystal ligands into the target protein binding pocket which is then refined by side-chain repacking and Monte Carlo energy minimization. Analyses of PoPSS revealed some drawbacks in ligand conformation generation and the scoring scheme used. Moreover, as PoPSS does not sample the ligand conformation after placing it in the binding pocket, it relies solely on conformation generation methods to produce native like conformations. To address these limitations of PoPSS method, we report here a modified approach named as PoPSS-Lite, where side-chain repacking was replaced by a simple grid-based energy minimization. This modification also allowed the sampling of terminal functional groups while keeping the core scaffold fixed. Furthermore, shape similarity calculations were improved by increasing the number of ligand conformations and using a different similarity metric. The performance of PoPSS-Lite was prospectively evaluated in D3R GC3. Comparison of PoPSS-Lite demonstrated superior performance over PoPSS and CDVS with lower mean and median RMSDs. Furthermore, comparison with other D3R GC3 pose prediction submissions revealed top performance for PoPSS-Lite. Our D3R GC3 result extends our perspective that ligand 3D shape similarity with known crystallographic information can be successfully used to predict the binding pose of ligands with unknown binding modes. Our D3R GC3 results further highlight the necessity for improvement in conformer generation methods in order to improve shape similarity guided pose prediction.
Collapse
Affiliation(s)
- Ashutosh Kumar
- Laboratory for Structural Bioinformatics, Center for Biosystems Dynamics Research, RIKEN, 1-7-22 Suehiro, Tsurumi, Yokohama, Kanagawa, 230-0045, Japan
| | - Kam Y J Zhang
- Laboratory for Structural Bioinformatics, Center for Biosystems Dynamics Research, RIKEN, 1-7-22 Suehiro, Tsurumi, Yokohama, Kanagawa, 230-0045, Japan.
| |
Collapse
|
10
|
Kumar A, Zhang KYJ. Advances in the Development of Shape Similarity Methods and Their Application in Drug Discovery. Front Chem 2018; 6:315. [PMID: 30090808 PMCID: PMC6068280 DOI: 10.3389/fchem.2018.00315] [Citation(s) in RCA: 81] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2018] [Accepted: 07/09/2018] [Indexed: 12/21/2022] Open
Abstract
Molecular similarity is a key concept in drug discovery. It is based on the assumption that structurally similar molecules frequently have similar properties. Assessment of similarity between small molecules has been highly effective in the discovery and development of various drugs. Especially, two-dimensional (2D) similarity approaches have been quite popular due to their simplicity, accuracy and efficiency. Recently, the focus has been shifted toward the development of methods involving the representation and comparison of three-dimensional (3D) conformation of small molecules. Among the 3D similarity methods, evaluation of shape similarity is now gaining attention for its application not only in virtual screening but also in molecular target prediction, drug repurposing and scaffold hopping. A wide range of methods have been developed to describe molecular shape and to determine the shape similarity between small molecules. The most widely used methods include atom distance-based methods, surface-based approaches such as spherical harmonics and 3D Zernike descriptors, atom-centered Gaussian overlay based representations. Several of these methods demonstrated excellent virtual screening performance not only retrospectively but also prospectively. In addition to methods assessing the similarity between small molecules, shape similarity approaches have been developed to compare shapes of protein structures and binding pockets. Additionally, shape comparisons between atomic models and 3D density maps allowed the fitting of atomic models into cryo-electron microscopy maps. This review aims to summarize the methodological advances in shape similarity assessment highlighting advantages, disadvantages and their application in drug discovery.
Collapse
Affiliation(s)
| | - Kam Y. J. Zhang
- Laboratory for Structural Bioinformatics, Center for Biosystems Dynamics Research, RIKEN, Yokohama, Japan
| |
Collapse
|
11
|
Kadukova M, Grudinin S. Docking of small molecules to farnesoid X receptors using AutoDock Vina with the Convex-PL potential: lessons learned from D3R Grand Challenge 2. J Comput Aided Mol Des 2017; 32:151-162. [DOI: 10.1007/s10822-017-0062-1] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2017] [Accepted: 09/08/2017] [Indexed: 10/18/2022]
|
12
|
da Silva Figueiredo Celestino Gomes P, Da Silva F, Bret G, Rognan D. Ranking docking poses by graph matching of protein-ligand interactions: lessons learned from the D3R Grand Challenge 2. J Comput Aided Mol Des 2017; 32:75-87. [PMID: 28766097 DOI: 10.1007/s10822-017-0046-1] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2017] [Accepted: 07/27/2017] [Indexed: 12/11/2022]
Abstract
A novel docking challenge has been set by the Drug Design Data Resource (D3R) in order to predict the pose and affinity ranking of a set of Farnesoid X receptor (FXR) agonists, prior to the public release of their bound X-ray structures and potencies. In a first phase, 36 agonists were docked to 26 Protein Data Bank (PDB) structures of the FXR receptor, and next rescored using the in-house developed GRIM method. GRIM aligns protein-ligand interaction patterns of docked poses to those of available PDB templates for the target protein, and rescore poses by a graph matching method. In agreement with results obtained during the previous 2015 docking challenge, we clearly show that GRIM rescoring improves the overall quality of top-ranked poses by prioritizing interaction patterns already visited in the PDB. Importantly, this challenge enables us to refine the applicability domain of the method by better defining the conditions of its success. We notably show that rescoring apolar ligands in hydrophobic pockets leads to frequent GRIM failures. In the second phase, 102 FXR agonists were ranked by decreasing affinity according to the Gibbs free energy of the corresponding GRIM-selected poses, computed by the HYDE scoring function. Interestingly, this fast and simple rescoring scheme provided the third most accurate ranking method among 57 contributions. Although the obtained ranking is still unsuitable for hit to lead optimization, the GRIM-HYDE scoring scheme is accurate and fast enough to post-process virtual screening data.
Collapse
Affiliation(s)
- Priscila da Silva Figueiredo Celestino Gomes
- Laboratoire d'Innovation Thérapeutique, UMR 7200 CNRS-Université de Strasbourg, 67400, Illkirch, France.,Instituto de Biofísica Carlos Chagas Filho, Universidade Federal do Rio de Janeiro, Rio de Janeiro, Brazil
| | - Franck Da Silva
- Laboratoire d'Innovation Thérapeutique, UMR 7200 CNRS-Université de Strasbourg, 67400, Illkirch, France
| | - Guillaume Bret
- Laboratoire d'Innovation Thérapeutique, UMR 7200 CNRS-Université de Strasbourg, 67400, Illkirch, France
| | - Didier Rognan
- Laboratoire d'Innovation Thérapeutique, UMR 7200 CNRS-Université de Strasbourg, 67400, Illkirch, France.
| |
Collapse
|
13
|
Drwal MN, Jacquemard C, Perez C, Desaphy J, Kellenberger E. Do Fragments and Crystallization Additives Bind Similarly to Drug-like Ligands? J Chem Inf Model 2017; 57:1197-1209. [PMID: 28414463 DOI: 10.1021/acs.jcim.6b00769] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The success of fragment-based drug design (FBDD) hinges upon the optimization of low-molecular-weight compounds (MW < 300 Da) with weak binding affinities to lead compounds with high affinity and selectivity. Usually, structural information from fragment-protein complexes is used to develop ideas about the binding mode of similar but drug-like molecules. In this regard, crystallization additives such as cryoprotectants or buffer components, which are highly abundant in crystal structures, are frequently ignored. Thus, the aim of this study was to investigate the information present in protein complexes with fragments as well as those with additives and how they relate to the binding modes of their drug-like counterparts. We present a thorough analysis of the binding modes of crystallographic additives, fragments, and drug-like ligands bound to four diverse targets of wide interest in drug discovery and highly represented in the Protein Data Bank: cyclin-dependent kinase 2, β-secretase 1, carbonic anhydrase 2, and trypsin. We identified a total of 630 unique molecules bound to the catalytic binding sites, among them 31 additives, 222 fragments, and 377 drug-like ligands. In general, we observed that, independent of the target, protein-fragment interaction patterns are highly similar to those of drug-like ligands and mostly cover the residues crucial for binding. Crystallographic additives are also able to show conserved binding modes and recover the residues important for binding in some of the cases. Moreover, we show evidence that the information from fragments and drug-like ligands can be applied to rescore docking poses in order to improve the prediction of binding modes.
Collapse
Affiliation(s)
- Malgorzata N Drwal
- Laboratoire d'Innovation Thérapeutique UMR 7200, CNRS-Université de Strasbourg , 74 Route du Rhin, 674000 Illkirch, France
| | - Célien Jacquemard
- Laboratoire d'Innovation Thérapeutique UMR 7200, CNRS-Université de Strasbourg , 74 Route du Rhin, 674000 Illkirch, France
| | - Carlos Perez
- Eli Lilly Research Laboratories , Avenida de la Industria 30, 28108 Alcobendas, Madrid, Spain
| | - Jérémy Desaphy
- Lilly Research Laboratories, Eli Lilly and Company , Lilly Corporate Center, Indianapolis, Indiana 46285, United States
| | - Esther Kellenberger
- Laboratoire d'Innovation Thérapeutique UMR 7200, CNRS-Université de Strasbourg , 74 Route du Rhin, 674000 Illkirch, France
| |
Collapse
|
14
|
Kumar M, Kaur T, Sharma A. Role of computational efficiency indices and pose clustering in effective decision making: An example of annulated furanones in Pf-DHFR space. Comput Biol Chem 2017; 67:48-61. [PMID: 28049061 DOI: 10.1016/j.compbiolchem.2016.12.007] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2016] [Revised: 11/04/2016] [Accepted: 12/22/2016] [Indexed: 10/20/2022]
|
15
|
Jiang X, Kumar A, Liu T, Zhang KYJ, Yang Q. A Novel Scaffold for Developing Specific or Broad-Spectrum Chitinase Inhibitors. J Chem Inf Model 2016; 56:2413-2420. [PMID: 28024404 DOI: 10.1021/acs.jcim.6b00615] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Chitinases play important roles in pathogen invasion, arthropod molting, plant defense, and human inflammation. Inhibition of the activity of a typical chitinase by small molecules is of significance in drug development and biological research. On the basis of a recent reported crystal structure of OfChtI, the insect chitinase derived from the pest Ostrinia furnacalis, we computationally identified 17 compounds from a library of over 4 million chemicals by two rounds virtual screening. Among these, three compounds from one chemical class inhibited the activity of OfChtI with single-digit-micromolar IC50 values, and one compound from another chemical class exhibited a broad inhibitory activity not only toward OfChtI but also toward bacterial, fungal, and human chitinases. A new scaffold was discovered, and a structure-inhibitory activity relationship was proposed. This work may provide a novel starting point for the development of specific or broad-spectrum chitinase inhibitors.
Collapse
Affiliation(s)
- Xi Jiang
- State Key Laboratory of Fine Chemicals and School of Life Science and Biotechnology, Dalian University of Technology , No. 2 Linggong Road, Dalian 116024, China
| | - Ashutosh Kumar
- Structural Bioinformatics Team, Center for Life Science Technologies, RIKEN , 1-7-22 Suehiro, Tsurumi, Yokohama, Kanagawa 230-0045, Japan
| | - Tian Liu
- State Key Laboratory of Fine Chemicals and School of Life Science and Biotechnology, Dalian University of Technology , No. 2 Linggong Road, Dalian 116024, China
| | - Kam Y J Zhang
- Structural Bioinformatics Team, Center for Life Science Technologies, RIKEN , 1-7-22 Suehiro, Tsurumi, Yokohama, Kanagawa 230-0045, Japan
| | - Qing Yang
- State Key Laboratory of Fine Chemicals and School of Life Science and Biotechnology, Dalian University of Technology , No. 2 Linggong Road, Dalian 116024, China.,Institute of Plant Protection, Chinese Academy of Agricultural Sciences , Beijing 100193, China
| |
Collapse
|
16
|
Kumar A, Zhang KYJ. Prospective evaluation of shape similarity based pose prediction method in D3R Grand Challenge 2015. J Comput Aided Mol Des 2016; 30:685-693. [DOI: 10.1007/s10822-016-9931-2] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2016] [Accepted: 07/25/2016] [Indexed: 01/23/2023]
|
17
|
Slynko I, Da Silva F, Bret G, Rognan D. Docking pose selection by interaction pattern graph similarity: application to the D3R grand challenge 2015. J Comput Aided Mol Des 2016; 30:669-683. [DOI: 10.1007/s10822-016-9930-3] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2016] [Accepted: 07/25/2016] [Indexed: 10/21/2022]
|
18
|
A pose prediction approach based on ligand 3D shape similarity. J Comput Aided Mol Des 2016; 30:457-69. [DOI: 10.1007/s10822-016-9923-2] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2016] [Accepted: 07/01/2016] [Indexed: 11/27/2022]
|
19
|
Carlson HA, Smith RD, Damm-Ganamet KL, Stuckey JA, Ahmed A, Convery MA, Somers DO, Kranz M, Elkins PA, Cui G, Peishoff CE, Lambert MH, Dunbar JB. CSAR 2014: A Benchmark Exercise Using Unpublished Data from Pharma. J Chem Inf Model 2016; 56:1063-77. [PMID: 27149958 DOI: 10.1021/acs.jcim.5b00523] [Citation(s) in RCA: 74] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The 2014 CSAR Benchmark Exercise was the last community-wide exercise that was conducted by the group at the University of Michigan, Ann Arbor. For this event, GlaxoSmithKline (GSK) donated unpublished crystal structures and affinity data from in-house projects. Three targets were used: tRNA (m1G37) methyltransferase (TrmD), Spleen Tyrosine Kinase (SYK), and Factor Xa (FXa). A particularly strong feature of the GSK data is its large size, which lends greater statistical significance to comparisons between different methods. In Phase 1 of the CSAR 2014 Exercise, participants were given several protein-ligand complexes and asked to identify the one near-native pose from among 200 decoys provided by CSAR. Though decoys were requested by the community, we found that they complicated our analysis. We could not discern whether poor predictions were failures of the chosen method or an incompatibility between the participant's method and the setup protocol we used. This problem is inherent to decoys, and we strongly advise against their use. In Phase 2, participants had to dock and rank/score a set of small molecules given only the SMILES strings of the ligands and a protein structure with a different ligand bound. Overall, docking was a success for most participants, much better in Phase 2 than in Phase 1. However, scoring was a greater challenge. No particular approach to docking and scoring had an edge, and successful methods included empirical, knowledge-based, machine-learning, shape-fitting, and even those with solvation and entropy terms. Several groups were successful in ranking TrmD and/or SYK, but ranking FXa ligands was intractable for all participants. Methods that were able to dock well across all submitted systems include MDock,1 Glide-XP,2 PLANTS,3 Wilma,4 Gold,5 SMINA,6 Glide-XP2/PELE,7 FlexX,8 and MedusaDock.9 In fact, the submission based on Glide-XP2/PELE7 cross-docked all ligands to many crystal structures, and it was particularly impressive to see success across an ensemble of protein structures for multiple targets. For scoring/ranking, submissions that showed statistically significant achievement include MDock1 using ITScore1,10 with a flexible-ligand term,11 SMINA6 using Autodock-Vina,12,13 FlexX8 using HYDE,14 and Glide-XP2 using XP DockScore2 with and without ROCS15 shape similarity.16 Of course, these results are for only three protein targets, and many more systems need to be investigated to truly identify which approaches are more successful than others. Furthermore, our exercise is not a competition.
Collapse
Affiliation(s)
- Heather A Carlson
- Department of Medicinal Chemistry, College of Pharmacy, University of Michigan , 428 Church St., Ann Arbor, Michigan 48109-1065, United States
| | - Richard D Smith
- Department of Medicinal Chemistry, College of Pharmacy, University of Michigan , 428 Church St., Ann Arbor, Michigan 48109-1065, United States
| | - Kelly L Damm-Ganamet
- Department of Medicinal Chemistry, College of Pharmacy, University of Michigan , 428 Church St., Ann Arbor, Michigan 48109-1065, United States
| | - Jeanne A Stuckey
- Center for Structural Biology, University of Michigan , 3358E Life Sciences Institute, 210 Washtenaw Ave., Ann Arbor, Michigan 48109-2216, United States
| | - Aqeel Ahmed
- Department of Medicinal Chemistry, College of Pharmacy, University of Michigan , 428 Church St., Ann Arbor, Michigan 48109-1065, United States
| | - Maire A Convery
- Computational and Structural Sciences, Medicines Research Centre, GlaxoSmithKline Research & Development , Gunnels Wood Road, Stevenage, Hertfordshire SG1 2NY, United Kingdom
| | - Donald O Somers
- Computational and Structural Sciences, Medicines Research Centre, GlaxoSmithKline Research & Development , Gunnels Wood Road, Stevenage, Hertfordshire SG1 2NY, United Kingdom
| | - Michael Kranz
- Computational and Structural Sciences, Medicines Research Centre, GlaxoSmithKline Research & Development , Gunnels Wood Road, Stevenage, Hertfordshire SG1 2NY, United Kingdom
| | - Patricia A Elkins
- Computational and Structural Sciences, GlaxoSmithKline Research & Development , 1250 South Collegeville Road, Collegeville, Pennsylvania 19426, United States
| | - Guanglei Cui
- Computational and Structural Sciences, GlaxoSmithKline Research & Development , 1250 South Collegeville Road, Collegeville, Pennsylvania 19426, United States
| | - Catherine E Peishoff
- Computational and Structural Sciences, GlaxoSmithKline Research & Development , 1250 South Collegeville Road, Collegeville, Pennsylvania 19426, United States
| | - Millard H Lambert
- Computational and Structural Sciences, GlaxoSmithKline Research & Development , 1250 South Collegeville Road, Collegeville, Pennsylvania 19426, United States
| | - James B Dunbar
- Department of Medicinal Chemistry, College of Pharmacy, University of Michigan , 428 Church St., Ann Arbor, Michigan 48109-1065, United States
| |
Collapse
|
20
|
Anighoro A, Bajorath J. Three-Dimensional Similarity in Molecular Docking: Prioritizing Ligand Poses on the Basis of Experimental Binding Modes. J Chem Inf Model 2016; 56:580-7. [DOI: 10.1021/acs.jcim.5b00745] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Affiliation(s)
- Andrew Anighoro
- Department of Life Science
Informatics, B-IT, LIMES Program Unit Chemical Biology and Medicinal
Chemistry, Rheinische Friedrich-Wilhelms-Universität, Dahlmannstrasse 2, D-53113 Bonn, Germany
| | - Jürgen Bajorath
- Department of Life Science
Informatics, B-IT, LIMES Program Unit Chemical Biology and Medicinal
Chemistry, Rheinische Friedrich-Wilhelms-Universität, Dahlmannstrasse 2, D-53113 Bonn, Germany
| |
Collapse
|
21
|
Zhu X, Shin WH, Kim H, Kihara D. Combined Approach of Patch-Surfer and PL-PatchSurfer for Protein-Ligand Binding Prediction in CSAR 2013 and 2014. J Chem Inf Model 2015; 56:1088-99. [PMID: 26691286 DOI: 10.1021/acs.jcim.5b00625] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The Community Structure-Activity Resource (CSAR) benchmark exercise provides a unique opportunity for researchers to objectively evaluate the performance of protein-ligand docking methods. Patch-Surfer and PL-PatchSurfer, molecular surface-based methods for predicting binding ligands of proteins developed in our group, were tested on both CSAR 2013 and 2014 benchmark exercises in combination with an empirical scoring function-based method, AutoDock, while we only participated in CSAR 2013 using Patch-Surfer. The prediction results for Phase 1 task in CSAR 2013 showed that Patch-Surfer was able to rank all the four designed binding proteins within top ranks, outperforming AutoDock Vina. In Phase 2 of 2013, PL-PatchSurfer correctly selected the correct ligand pose for two target proteins. PL-PatchSurfer performed reasonably well in ranking ligands according to their binding affinity and in selecting near-native ligand poses in 2013 Phase 3 and 2014 Phase 1, respectively, although AutoDock Vina showed better performance. Lastly, in the 2014 Phase 2 exercise, the PL-PatchSurfer scores computed for ligands to target protein pairs correlated well with their pIC50 values, which was better or comparable to results by other participants. Overall, our methods showed fairly good performance in CSAR 2013 and 2014. Unique characteristics of the methods are discussed in comparison with AutoDock.
Collapse
Affiliation(s)
- Xiaolei Zhu
- School of Life Science, Anhui University , Hefei, Anhui 230601, China.,Department of Biology Science, Purdue University , West Lafayette, Indiana 47907, United States
| | - Woong-Hee Shin
- Department of Biology Science, Purdue University , West Lafayette, Indiana 47907, United States
| | - Hyungrae Kim
- Department of Biology Science, Purdue University , West Lafayette, Indiana 47907, United States
| | - Daisuke Kihara
- Department of Biology Science, Purdue University , West Lafayette, Indiana 47907, United States.,Department of Computer Science, Purdue University , West Lafayette, Indiana 47907, United States
| |
Collapse
|