1
|
Pala D, Clark DE. Caught between a ROCK and a hard place: current challenges in structure-based drug design. Drug Discov Today 2024; 29:104106. [PMID: 39029868 DOI: 10.1016/j.drudis.2024.104106] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2024] [Revised: 06/27/2024] [Accepted: 07/13/2024] [Indexed: 07/21/2024]
Abstract
The discipline of structure-based drug design (SBDD) is several decades old and it is tempting to think that the proliferation of experimental structures for many drug targets might make computer-aided drug design (CADD) straightforward. However, this is far from true. In this review, we illustrate some of the challenges that CADD scientists face every day in their work, even now. We use Rho-associated protein kinase (ROCK), and public domain structures and data, as an example to illustrate some of the challenges we have experienced during our project targeting this protein. We hope that this will help to prevent unrealistic expectations of what CADD can accomplish and to educate non-CADD scientists regarding the challenges still facing their CADD colleagues.
Collapse
Affiliation(s)
- Daniele Pala
- Medicinal Chemistry and Drug Design Technologies Department, Chiesi Farmaceutici S.p.A, Research Center, Largo Belloli 11/a, 43122 Parma, Italy
| | - David E Clark
- Charles River, 6-9 Spire Green Centre, Flex Meadow, Harlow CM19 5TR, UK.
| |
Collapse
|
2
|
Zhao M, Yu W, MacKerell AD. Enhancing SILCS-MC via GPU Acceleration and Ligand Conformational Optimization with Genetic and Parallel Tempering Algorithms. J Phys Chem B 2024; 128:7362-7375. [PMID: 39031121 PMCID: PMC11294009 DOI: 10.1021/acs.jpcb.4c03045] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/22/2024]
Abstract
In the domain of computer-aided drug design, achieving precise and accurate estimates of ligand-protein binding is paramount in the context of screening extensive drug libraries and performing ligand optimization. A fundamental aspect of the SILCS (site identification by ligand competitive saturation) methodology lies in the generation of comprehensive 3D free-energy functional group affinity maps (FragMaps), encompassing the entirety of the target molecule structure. These FragMaps offer an intricate landscape of functional group affinities across the protein, bilayer, or RNA, acting as the basis for subsequent SILCS-Monte Carlo (MC) simulations wherein ligands are docked to the target molecule. To augment the efficiency and breadth of ligand sampling capabilities, we implemented an improved SILCS-MC methodology. By harnessing the parallel computing capability of GPUs, our approach facilitates concurrent calculations over multiple ligands and binding sites, markedly enhancing the computational efficiency. Moreover, the integration of a genetic algorithm (GA) with MC allows us to employ an evolutionary approach to perform ligand sampling, assuring enhanced convergence characteristics. In addition, the potential utility of parallel tempering (PT) to improve sampling was investigated. Implementation of SILCS-MC on GPU architecture is shown to accelerate the speed of SILCS-MC calculations by over 2-orders of magnitude. Use of GA and PT yield improvements over Markov-chain MC, increasing the precision of the resultant docked orientations and binding free energies, though the extent of improvements is relatively small. Accordingly, significant improvements in speed are obtained through the GPU implementation with minor improvements in the precision of the docking obtained via the tested GA and PT algorithms.
Collapse
Affiliation(s)
- Mingtian Zhao
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland, School of Pharmacy, 20 Penn St., Baltimore, Maryland 21201, USA
| | - Wenbo Yu
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland, School of Pharmacy, 20 Penn St., Baltimore, Maryland 21201, USA
| | - Alexander D. MacKerell
- Computer Aided Drug Design Center, Department of Pharmaceutical Sciences, University of Maryland, School of Pharmacy, 20 Penn St., Baltimore, Maryland 21201, USA
| |
Collapse
|
3
|
Tan S, Gong X, Liu H, Yao X. Identification of novel LRRK2 inhibitors by structure-based virtual screening and alchemical free energy calculation. Phys Chem Chem Phys 2024; 26:19775-19786. [PMID: 38984923 DOI: 10.1039/d4cp01762e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/11/2024]
Abstract
The Leucine-rich repeat kinase 2 (LRRK2) target has been identified as a promising drug target for Parkinson's disease (PD) treatment. This study focuses on optimizing the activity of LRRK2 inhibitors using alchemical relative binding free energy (RBFE) calculations. Initially, we assessed various free energy calculation methods across different LRRK2 kinase inhibitor scaffolds. The results indicate that alchemical free energy calculations are promising for prospective predictions on LRRK2 inhibitors, especially for the aminopyrimidine scaffold with an RMSE of 1.15 kcal mol-1 and Rp of 0.83. Following this, we optimized a potent LRRK2 kinase inhibitor identified from previous virtual screenings, featuring a novel scaffold. Guided by RBFE predictions using alchemical methods, this optimization led to the discovery of compound LY2023-001. This compound, with a [1,2,4]triazolo[5,6-b]indole scaffold, exhibited enhanced inhibitory activity against G2019S LRRK2 (IC50 = 12.9 nM). Molecular dynamics (MD) simulations revealed that LY2023-001 formed stable hydrogen bonds with Glu1948, and Ala1950 in the G2019S LRRK2 protein. Additionally, its phenyl substituents engage in strong electrostatic interactions with Lys1906 and van der Waals interactions with Leu1885, Phe1890, Val1893, Ile1933, Met1947, Leu1949, Leu2001, Ala2016, and Asp2017. Our findings underscore the potential of computational methods in the successful optimization of small molecules, offering important insights for the development of novel LRRK2 inhibitors.
Collapse
Affiliation(s)
- Shuoyan Tan
- College of Chemistry and Chemical Engineering, Lanzhou University, Lanzhou 730000, China
| | - Xiaoqing Gong
- Faculty of Applied Sciences, Macao Polytechnic University, Macao SAR, 999078, China.
| | - Huanxiang Liu
- Faculty of Applied Sciences, Macao Polytechnic University, Macao SAR, 999078, China.
| | - Xiaojun Yao
- Faculty of Applied Sciences, Macao Polytechnic University, Macao SAR, 999078, China.
| |
Collapse
|
4
|
Hahn DF, Gapsys V, de Groot BL, Mobley DL, Tresadern G. Current State of Open Source Force Fields in Protein-Ligand Binding Affinity Predictions. J Chem Inf Model 2024; 64:5063-5076. [PMID: 38895959 PMCID: PMC11234369 DOI: 10.1021/acs.jcim.4c00417] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2024] [Revised: 04/23/2024] [Accepted: 04/25/2024] [Indexed: 06/21/2024]
Abstract
In drug discovery, the in silico prediction of binding affinity is one of the major means to prioritize compounds for synthesis. Alchemical relative binding free energy (RBFE) calculations based on molecular dynamics (MD) simulations are nowadays a popular approach for the accurate affinity ranking of compounds. MD simulations rely on empirical force field parameters, which strongly influence the accuracy of the predicted affinities. Here, we evaluate the ability of six different small-molecule force fields to predict experimental protein-ligand binding affinities in RBFE calculations on a set of 598 ligands and 22 protein targets. The public force fields OpenFF Parsley and Sage, GAFF, and CGenFF show comparable accuracy, while OPLS3e is significantly more accurate. However, a consensus approach using Sage, GAFF, and CGenFF leads to accuracy comparable to OPLS3e. While Parsley and Sage are performing comparably based on aggregated statistics across the whole dataset, there are differences in terms of outliers. Analysis of the force field reveals that improved parameters lead to significant improvement in the accuracy of affinity predictions on subsets of the dataset involving those parameters. Lower accuracy can not only be attributed to the force field parameters but is also dependent on input preparation and sampling convergence of the calculations. Especially large perturbations and nonconverged simulations lead to less accurate predictions. The input structures, Gromacs force field files, as well as the analysis Python notebooks are available on GitHub.
Collapse
Affiliation(s)
- David F. Hahn
- Computational
Chemistry, Janssen Research & Development, Turnhoutseweg 30, Beerse 2340, Belgium
| | - Vytautas Gapsys
- Computational
Chemistry, Janssen Research & Development, Turnhoutseweg 30, Beerse 2340, Belgium
- Computational
Biomolecular Dynamics Group, Max Planck
Institute for Multidisciplinary Sciences, Am Fassberg 11, Göttingen 37077, Germany
| | - Bert L. de Groot
- Computational
Biomolecular Dynamics Group, Max Planck
Institute for Multidisciplinary Sciences, Am Fassberg 11, Göttingen 37077, Germany
| | - David L. Mobley
- Department
of Chemistry, University of California, Irvine, California 92697, United States
- Department
of Pharmaceutical Sciences, University of
California, Irvine, California 92697, United States
| | - Gary Tresadern
- Computational
Chemistry, Janssen Research & Development, Turnhoutseweg 30, Beerse 2340, Belgium
| |
Collapse
|
5
|
Singh M, Verma H, Gera N, Baddipadige R, Choudhary S, Bhandu P, Silakari O. Evaluation of Cordyceps militaris steroids as anti-inflammatory agents to combat the Covid-19 cytokine storm: a bioinformatics and structure-based drug designing approach. J Biomol Struct Dyn 2024; 42:5159-5177. [PMID: 37551029 DOI: 10.1080/07391102.2023.2245039] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2022] [Accepted: 06/08/2023] [Indexed: 08/09/2023]
Abstract
Since the SARS-CoV-2 epidemic, researchers have been working on figuring out ways to tackle multi-organ failure and hyperinflation, which are brought on by a cytokine storm. Angiotensin-converting enzyme 2 (ACE2), a SARS-CoV-2 spike glycoprotein's cellular receptor, is involved in complicated molecular processes that result in hyperinflammation. Cordyceps militaris is one of the traditional Chinese medicines that is used as an immune booster, and it has exhibited efficacy in lowering blood glucose levels, seminal emissions, and infertility. In the current study, we explored the potential of Cordyceps militaris steroids as key agents in managing the anger of cytokine storm in Covid-19 using network ethnopharmacological techniques and structure-based drug designing approaches. The steroids present in Cordyceps militaris were initially screened against the targets involved in inflammatory pathways. The results revealed that out of 16 steroids, 5 may be effective against 17 inflammatory pathways by targeting 11 pathological proteins. Among the five steroids, beta-sitosterol, Cholest-5-en-3β-ol, 3β, and 7α-Dihydroxycholest-5-ene were found to interact with thrombin (F2), an important protein reported to reduce the severity of inflammatory mediators and Cholest-4-en-3-one was found to target Glucocorticoid receptor (NR3C1). The top docked steroid displayed key interactions with both targets, which retained key interactions throughout the 100 ns simulation period. These compounds were also shown high binding free energy scores in water swap studies. Based on obtained results the current study suggests the use of Cordyceps militaris as an add-on therapy that may reduce the progression of inflammatory co-morbidities among patients infected with SARS-CoV-2.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Manmeet Singh
- Molecular Modeling Lab (MML), Department of Pharmaceutical Sciences and Drug Research, Punjabi University, Patiala, Punjab, India
| | - Himanshu Verma
- Molecular Modeling Lab (MML), Department of Pharmaceutical Sciences and Drug Research, Punjabi University, Patiala, Punjab, India
| | - Narendra Gera
- Molecular Modeling Lab (MML), Department of Pharmaceutical Sciences and Drug Research, Punjabi University, Patiala, Punjab, India
| | - Raju Baddipadige
- Molecular Modeling Lab (MML), Department of Pharmaceutical Sciences and Drug Research, Punjabi University, Patiala, Punjab, India
| | - Shalki Choudhary
- Molecular Modeling Lab (MML), Department of Pharmaceutical Sciences and Drug Research, Punjabi University, Patiala, Punjab, India
| | - Priyanka Bhandu
- Molecular Modeling Lab (MML), Department of Pharmaceutical Sciences and Drug Research, Punjabi University, Patiala, Punjab, India
| | - Om Silakari
- Molecular Modeling Lab (MML), Department of Pharmaceutical Sciences and Drug Research, Punjabi University, Patiala, Punjab, India
| |
Collapse
|
6
|
Lim CP, Leow CH, Lim HT, Kok BH, Chuah C, Oliveira JIN, Jones M, Leow CY. Insights into structural vaccinology harnessed for universal coronavirus vaccine development. Clin Exp Vaccine Res 2024; 13:202-217. [PMID: 39144127 PMCID: PMC11319108 DOI: 10.7774/cevr.2024.13.3.202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2024] [Revised: 05/10/2024] [Accepted: 05/15/2024] [Indexed: 08/16/2024] Open
Abstract
Structural vaccinology is pivotal in expediting vaccine design through high-throughput screening of immunogenic antigens. Leveraging the structural and functional characteristics of antigens and immune cell receptors, this approach employs protein structural comparison to identify conserved patterns in key pathogenic components. Molecular modeling techniques, including homology modeling and molecular docking, analyze specific three-dimensional (3D) structures and protein interactions and offer valuable insights into the 3D interactions and binding affinity between vaccine candidates and target proteins. In this review, we delve into the utilization of various immunoinformatics and molecular modeling tools to streamline the development of broad-protective vaccines against coronavirus disease 2019 variants. Structural vaccinology significantly enhances our understanding of molecular interactions between hosts and pathogens. By accelerating the pace of developing effective and targeted vaccines, particularly against the rapidly mutating severe acute respiratory syndrome coronavirus 2 and other prevalent infectious diseases, this approach stands at the forefront of advancing immunization strategies. The combination of computational techniques and structural insights not only facilitates the identification of potential vaccine candidates but also contributes to the rational design of vaccines, fostering a more efficient and targeted approach to combatting infectious diseases.
Collapse
Affiliation(s)
- Chin Peng Lim
- School of Pharmaceutical Sciences, Universiti Sains Malaysia, Gelugor, Malaysia
- Institute for Research in Molecular Medicine, Universiti Sains Malaysia, Gelugor, Malaysia
| | - Chiuan Herng Leow
- Institute for Research in Molecular Medicine, Universiti Sains Malaysia, Gelugor, Malaysia
| | - Hui Ting Lim
- Institute for Research in Molecular Medicine, Universiti Sains Malaysia, Gelugor, Malaysia
| | - Boon Hui Kok
- Institute for Research in Molecular Medicine, Universiti Sains Malaysia, Gelugor, Malaysia
| | - Candy Chuah
- Faculty of Medicine, Asian Institute of Medical Science and Technology University, Bedong, Malaysia
| | - Jonas Ivan Nobre Oliveira
- Department of Biophysics and Pharmacology, Bioscience Center, Federal University of Rio Grande do Norte, Natal, Brazil
| | - Malcolm Jones
- School of Veterinary Science, The University of Queensland, Gatton, Australia
| | - Chiuan Yee Leow
- School of Pharmaceutical Sciences, Universiti Sains Malaysia, Gelugor, Malaysia
| |
Collapse
|
7
|
Prasanna A, Karunakar P, Pillai A, Mukundan S, Y V M, Balaji R, Niranjan V, Skariyachan S, Narayanappa R. Screening of bioactive compounds from selected mushroom species against putative drug targets in Mycobacterium tuberculosis: a multi-target approach. J Biomol Struct Dyn 2024:1-16. [PMID: 38895953 DOI: 10.1080/07391102.2024.2335292] [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: 07/31/2023] [Accepted: 03/20/2024] [Indexed: 06/21/2024]
Abstract
Mycobacterium tuberculosis (Mtb) is a notorious pathogen that causes one of the highest mortalities globally. Due to a pressing demand to identify novel therapeutic alternatives, the present study aims to focus on screening the putative drug targets and prioritizing their role in antibacterial drug development. The most vital proteins involved in the Biotin biosynthesis pathway and the Lipoarabinomannan (LAM) pathway such as biotin synthase (bioB) and alpha-(1->6)-mannopyranosyltransferase A (mptA) respectively, along with other essential virulence proteins of Mtb were selected as drug targets. Among these, the ones without native structures were modelled and validated using standard bioinformatics tools. Further, the interactions were performed with naturally available lead molecules present in selected mushroom species such as Agaricus bisporus, Pleurotus djamor, Hypsizygus ulmarius. Through Gas Chromatography-Mass Spectrometry (GC-MS), 15 bioactive compounds from the methanolic extract of mushrooms were identified. Further, 4 were selected based on drug-likeness and pharmacokinetic screening for molecular docking analysis against our prioritized targets wherein Benz[e]azulene from Pleurotus djamor illustrated a good binding affinity with a LF rank score of -9.036 kcal mol -1 against nuoM (NADH quinone oxidoreductase subunit M) and could be used as a prospective candidate in order to combat Tuberculosis (TB). Furthermore, the stability of the complex are validated using MD Simulations and subsequently, the binding free energy was calculated using MM-GBSA analysis. Thus, the current in silico analysis suggests a promising role of compounds extracted from mushrooms in tackling the TB burden.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Akshatha Prasanna
- Department of Biotechnology, Dayananda Sagar College of Engineering (Affiliated to Visvesvaraya Technological University, Belagavi), Bengaluru, Karnataka, India
| | - Prashantha Karunakar
- Department of Biotechnology, Dayananda Sagar College of Engineering (Affiliated to Visvesvaraya Technological University, Belagavi), Bengaluru, Karnataka, India
| | - Anushka Pillai
- Department of Biotechnology, Dayananda Sagar College of Engineering (Affiliated to Visvesvaraya Technological University, Belagavi), Bengaluru, Karnataka, India
| | - Shreyashree Mukundan
- Department of Biotechnology, Dayananda Sagar College of Engineering (Affiliated to Visvesvaraya Technological University, Belagavi), Bengaluru, Karnataka, India
| | - Mansi Y V
- Department of Biotechnology, Dayananda Sagar College of Engineering (Affiliated to Visvesvaraya Technological University, Belagavi), Bengaluru, Karnataka, India
| | - Renu Balaji
- Department of Biotechnology, Dayananda Sagar College of Engineering (Affiliated to Visvesvaraya Technological University, Belagavi), Bengaluru, Karnataka, India
| | - Vidya Niranjan
- Department of Biotechnology, RV College of Engineering, Bengaluru, Karnataka, India
| | - Sinosh Skariyachan
- Department of Microbiology, St. Pius X College Rajapuram, Kasaragod, Kerala, India
| | - Rajeswari Narayanappa
- Department of Biotechnology, Dayananda Sagar College of Engineering (Affiliated to Visvesvaraya Technological University, Belagavi), Bengaluru, Karnataka, India
| |
Collapse
|
8
|
Zampieri D, Romano M, Fortuna S, Amata E, Dichiara M, Cosentino G, Marrazzo A, Mamolo MG. Design, Synthesis, and Cytotoxic Assessment of New Haloperidol Analogues as Potential Anticancer Compounds Targeting Sigma Receptors. Molecules 2024; 29:2697. [PMID: 38893570 PMCID: PMC11173765 DOI: 10.3390/molecules29112697] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2024] [Revised: 05/31/2024] [Accepted: 06/03/2024] [Indexed: 06/21/2024] Open
Abstract
Sigma receptors (SRs), including SR1 and SR2 subtypes, have attracted increasing interest in recent years due to their involvement in a wide range of activities, including the modulation of opioid analgesia, neuroprotection, and potential anticancer activity. In this context, haloperidol (HAL), a commonly used antipsychotic drug, also possesses SR activity and cytotoxic effects. Herein, we describe the identification of novel SR ligands, obtained by a chemical hybridization approach. There wereendowed with pan-affinity for both SR subtypes and evaluated their potential anticancer activity against SH-SY5Y and HUH-7 cancer cell lines. Through a chemical hybridization approach, we identified novel compounds (4d, 4e, 4g, and 4j) with dual affinity for SR1 and SR2 receptors. These compounds were subjected to cytotoxicity testing using a resazurin assay. The results revealed potent cytotoxic effects against both cancer cell lines, with IC50 values comparable to HAL. Interestingly, the cytotoxic potency of the novel compounds resembled that of the SR1 antagonist HAL rather than the SR2 agonist siramesine (SRM), indicating the potential role of SR1 antagonism in their mechanism of action. The further exploration of their structure-activity relationships and their evaluation in additional cancer cell lines will elucidate their therapeutic potential and may pave the way for the development of novel anticancer agents that target SRs.
Collapse
Affiliation(s)
- Daniele Zampieri
- Department of Chemical and Pharmaceutical Sciences, University of Trieste, Via Giorgieri 1, 34127 Trieste, Italy;
| | - Maurizio Romano
- Department of Life Sciences, University of Trieste, Via Valerio 28, 34127 Trieste, Italy;
| | - Sara Fortuna
- Cresset-New Cambridge House, Bassingbourn Road, Litlington, Cambridge SG8 0SS, UK;
| | - Emanuele Amata
- Department of Drug and Health Sciences, University of Catania, Viale Doria 6, 95125 Catania, Italy; (E.A.); (M.D.); (G.C.); (A.M.)
| | - Maria Dichiara
- Department of Drug and Health Sciences, University of Catania, Viale Doria 6, 95125 Catania, Italy; (E.A.); (M.D.); (G.C.); (A.M.)
| | - Giuseppe Cosentino
- Department of Drug and Health Sciences, University of Catania, Viale Doria 6, 95125 Catania, Italy; (E.A.); (M.D.); (G.C.); (A.M.)
| | - Agostino Marrazzo
- Department of Drug and Health Sciences, University of Catania, Viale Doria 6, 95125 Catania, Italy; (E.A.); (M.D.); (G.C.); (A.M.)
| | - Maria Grazia Mamolo
- Department of Chemical and Pharmaceutical Sciences, University of Trieste, Via Giorgieri 1, 34127 Trieste, Italy;
| |
Collapse
|
9
|
Woods CJ, Hedges LO, Mulholland AJ, Malaisree M, Tosco P, Loeffler HH, Suruzhon M, Burman M, Bariami S, Bosisio S, Calabro G, Clark F, Mey ASJS, Michel J. Sire: An interoperability engine for prototyping algorithms and exchanging information between molecular simulation programs. J Chem Phys 2024; 160:202503. [PMID: 38814008 DOI: 10.1063/5.0200458] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2024] [Accepted: 04/29/2024] [Indexed: 05/31/2024] Open
Abstract
Sire is a Python/C++ library that is used both to prototype new algorithms and as an interoperability engine for exchanging information between molecular simulation programs. It provides a collection of file parsers and information converters that together make it easier to combine and leverage the functionality of many other programs and libraries. This empowers researchers to use sire to write a single script that can, for example, load a molecule from a PDBx/mmCIF file via Gemmi, perform SMARTS searches via RDKit, parameterize molecules using BioSimSpace, run GPU-accelerated molecular dynamics via OpenMM, and then display the resulting dynamics trajectory in a NGLView Jupyter notebook 3D molecular viewer. This functionality is built on by BioSimSpace, which uses sire's molecular information engine to interconvert with programs such as GROMACS, NAMD, Amber, and AmberTools for automated molecular parameterization and the running of molecular dynamics, metadynamics, and alchemical free energy workflows. Sire comes complete with a powerful molecular information search engine, plus trajectory loading and editing, analysis, and energy evaluation engines. This, when combined with an in-built computer algebra system, gives substantial flexibility to researchers to load, search for, edit, and combine molecular information from multiple sources and use that to drive novel algorithms by combining functionality from other programs. Sire is open source (GPL3) and is available via conda and at a free Jupyter notebook server at https://try.openbiosim.org. Sire is supported by the not-for-profit OpenBioSim community interest company.
Collapse
Affiliation(s)
- Christopher J Woods
- Advanced Computing Research Centre, University of Bristol, Bristol, United Kingdom
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, United Kingdom
- OpenBioSim Community Interest Company, Edinburgh, United Kingdom
| | - Lester O Hedges
- Advanced Computing Research Centre, University of Bristol, Bristol, United Kingdom
- OpenBioSim Community Interest Company, Edinburgh, United Kingdom
| | - Adrian J Mulholland
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, United Kingdom
| | - Maturos Malaisree
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, United Kingdom
| | | | | | | | - Matthew Burman
- OpenBioSim Community Interest Company, Edinburgh, United Kingdom
| | - Sofia Bariami
- EaStCHEM School of Chemistry, University of Edinburgh, Edinburgh, United Kingdom
| | - Stefano Bosisio
- EaStCHEM School of Chemistry, University of Edinburgh, Edinburgh, United Kingdom
| | - Gaetano Calabro
- EaStCHEM School of Chemistry, University of Edinburgh, Edinburgh, United Kingdom
| | - Finlay Clark
- EaStCHEM School of Chemistry, University of Edinburgh, Edinburgh, United Kingdom
| | - Antonia S J S Mey
- EaStCHEM School of Chemistry, University of Edinburgh, Edinburgh, United Kingdom
| | - Julien Michel
- OpenBioSim Community Interest Company, Edinburgh, United Kingdom
- EaStCHEM School of Chemistry, University of Edinburgh, Edinburgh, United Kingdom
| |
Collapse
|
10
|
Zeng J, Qian Y. Adaptive lambda schemes for efficient relative binding free energy calculation. J Comput Chem 2024; 45:855-862. [PMID: 38153254 DOI: 10.1002/jcc.27287] [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: 07/19/2023] [Revised: 11/13/2023] [Accepted: 12/03/2023] [Indexed: 12/29/2023]
Abstract
The relative free energy perturbation (RFEP) calculation is one of the most theoretically sound computational chemistry approaches for the binding affinity prediction. However, its application is often hindered by the complexity of the calculation choices and the requirement of expertise in the field. Improper lambda scheme of RFEP may result in deviations from an accurate description of the perturbation process and is prone to erroneous affinity predictions. To address such challenges, an automated adaptive lambda method is proposed where the adaptive lambda schemes are obtained through a split-and-merge algorithm based on the pilot runs. The newly established workflow along with a series of improvements to the perturbation settings increases the consistency of the RFEP calculation results. Comparing the pilot and adaptive lambda schemes, the latter demonstrated improvements in convergence and reproducibility and lowered the mean unsigned error and the root-mean-square error. Overall, the adaptive lambda method is a reliable and robust choice to predict small molecule relative binding free energy and can be capitalized to benefit routine RFEP calculations for drug discovery projects.
Collapse
Affiliation(s)
- Jin Zeng
- AIxplorerBio Biotech Co., Ltd., Jiaxing, Zhejiang Province, China
| | - Yue Qian
- Viva Biotech (Shanghai) Ltd., Shanghai, China
| |
Collapse
|
11
|
Falsaperla R, Sortino V, Marino SD, Collotta AD, Gammeri C, Sipala FM, Volti GL, Ruggieri M, Ronsisvalle S. Molecular Dynamic Simulations to Determine Individualized Therapy: Tetrabenazine for the GNAO1 Encephalopathy E246K Variant. Mol Diagn Ther 2024; 28:329-337. [PMID: 38581611 DOI: 10.1007/s40291-024-00706-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/17/2024] [Indexed: 04/08/2024]
Abstract
INTRODUCTION GNAO1 encephalopathy is characterized by severe hypotonia, psychomotor retardation, epilepsy, and movement disorders. Genetic variations in GNAO1 have been linked to neurological symptoms including movement disorders like dystonia. The correlation between the E246K mutation in the Gα subunit and aberrant signal transduction of G proteins has been established but no data are reported regarding the efficacy of medical treatment with tetrabenazine. METHODS Molecular modeling studies were performed to elucidate the molecular mechanisms underlying this mutation. We developed drug efficacy models using molecular dynamic simulations that replicated the behavior of wild-type and mutated proteins in the presence or absence of ligands. RESULTS AND DISCUSSION We demonstrated that the absence of the mutation leads to normal signal transduction upon receptor activation by the endogenous ligand, but not in the presence of tetrabenazine. In contrast, the presence of the mutation resulted in abnormal signal transduction in the presence of the endogenous ligand, which was corrected by the drug tetrabenazine. Tetrabenazine was identified as a promising therapeutic option for pediatric patients suffering from encephalopathy due to an E246K mutation in the GNAO1 gene validated through molecular dynamics. This is a potential first example of the use of this technique in a rare neurological pediatric disease.
Collapse
Affiliation(s)
- Raffaele Falsaperla
- Neonatal Intensive Care Unit and Neonatal Accompaniment Unit, Azienda Ospedaliero-Universitaria Policlinico "Rodolico-San Marco", San Marco Hospital, University of Catania, Catania, Italy.
- Unit of Clinical Paediatrics, Azienda Ospedaliero-Universitaria Policlinico, "Rodolico-San Marco", San Marco Hospital, Catania, Italy.
- Pediatric Clinic, University of Ferrara, Ferrara, Italy.
| | - Vincenzo Sortino
- Unit of Clinical Paediatrics, Azienda Ospedaliero-Universitaria Policlinico, "Rodolico-San Marco", San Marco Hospital, Catania, Italy
| | - Simona Domenica Marino
- Unit of Clinical Paediatrics, Azienda Ospedaliero-Universitaria Policlinico, "Rodolico-San Marco", San Marco Hospital, Catania, Italy
| | - Ausilia Desiree Collotta
- Unit of Clinical Paediatrics, Azienda Ospedaliero-Universitaria Policlinico, "Rodolico-San Marco", San Marco Hospital, Catania, Italy
- Postgraduate Training Program in Pediatrics, Department of Clinical and Experimental Medicine, University of Catania, Catania, Italy
| | - Carmela Gammeri
- Postgraduate Training Program in Pediatrics, Department of Clinical and Experimental Medicine, University of Catania, Catania, Italy
| | - Federica Maria Sipala
- Department of Drug and Health Sciences, University of Catania, Viale A. Doria 6, 95125, Catania, Italy
| | - Giovanni Li Volti
- Department of Biomedical and Biotechnological Sciences, University of Catania, Via S. Sofia 87, 95125, Catania, Italy
| | - Martino Ruggieri
- Unit of Clinical Pediatrics, Unit of Rare Diseases, AOU "Policlinico", PO "G. Rodolico", University of Catania, Catania, Italy
| | - Simone Ronsisvalle
- Department of Drug and Health Sciences, University of Catania, Viale A. Doria 6, 95125, Catania, Italy
| |
Collapse
|
12
|
Tan LH, Kwoh CK, Mu Y. RmsdXNA: RMSD prediction of nucleic acid-ligand docking poses using machine-learning method. Brief Bioinform 2024; 25:bbae166. [PMID: 38695120 PMCID: PMC11063749 DOI: 10.1093/bib/bbae166] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2023] [Revised: 03/15/2024] [Accepted: 03/19/2024] [Indexed: 05/04/2024] Open
Abstract
Small molecule drugs can be used to target nucleic acids (NA) to regulate biological processes. Computational modeling methods, such as molecular docking or scoring functions, are commonly employed to facilitate drug design. However, the accuracy of the scoring function in predicting the closest-to-native docking pose is often suboptimal. To overcome this problem, a machine learning model, RmsdXNA, was developed to predict the root-mean-square-deviation (RMSD) of ligand docking poses in NA complexes. The versatility of RmsdXNA has been demonstrated by its successful application to various complexes involving different types of NA receptors and ligands, including metal complexes and short peptides. The predicted RMSD by RmsdXNA was strongly correlated with the actual RMSD of the docked poses. RmsdXNA also outperformed the rDock scoring function in ranking and identifying closest-to-native docking poses across different structural groups and on the testing dataset. Using experimental validated results conducted on polyadenylated nuclear element for nuclear expression triplex, RmsdXNA demonstrated better screening power for the RNA-small molecule complex compared to rDock. Molecular dynamics simulations were subsequently employed to validate the binding of top-scoring ligand candidates selected by RmsdXNA and rDock on MALAT1. The results showed that RmsdXNA has a higher success rate in identifying promising ligands that can bind well to the receptor. The development of an accurate docking score for a NA-ligand complex can aid in drug discovery and development advancements. The code to use RmsdXNA is available at the GitHub repository https://github.com/laiheng001/RmsdXNA.
Collapse
Affiliation(s)
- Lai Heng Tan
- Interdisciplinary Graduate School, Nanyang Technological University, 61 Nanyang Drive, 637335 Singapore, Singapore
| | - Chee Keong Kwoh
- School of Computer Science and Engineering, Nanyang Technological University, 50 Nanyang Avenue, 639798 Singapore, Singapore
| | - Yuguang Mu
- School of Biological Sciences, Nanyang Technological University, 60 Nanyang Drive, 637551 Singapore, Singapore
| |
Collapse
|
13
|
Karrenbrock M, Rizzi V, Procacci P, Gervasio FL. Addressing Suboptimal Poses in Nonequilibrium Alchemical Calculations. J Phys Chem B 2024; 128:1595-1605. [PMID: 38323915 DOI: 10.1021/acs.jpcb.3c06516] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2024]
Abstract
Alchemical transformations can be used to quantitatively estimate absolute binding free energies at a reasonable computational cost. However, most of the approaches currently in use require knowledge of the correct (crystallographic) pose. In this paper, we present a combined Hamiltonian replica exchange nonequilibrium alchemical method that allows us to reliably calculate absolute binding free energies, even when starting from suboptimal initial binding poses. Performing a preliminary Hamiltonian replica exchange enhances the sampling of slow degrees of freedom of the ligand and the target, allowing the system to populate the correct binding pose when starting from an approximate docking pose. We apply the method on 6 ligands of the first bromodomain of the BRD4 bromodomain-containing protein. For each ligand, we start nonequilibrium alchemical transformations from both the crystallographic pose and the top-scoring docked pose that are often significantly different. We show that the method produces statistically equivalent binding free energies, making it a useful tool for computational drug discovery pipelines.
Collapse
Affiliation(s)
- Maurice Karrenbrock
- School of Pharmaceutical Sciences, University of Geneva, Rue Michel-Servet 1, CH-1206 Geneva, Switzerland
| | - Valerio Rizzi
- School of Pharmaceutical Sciences, University of Geneva, Rue Michel-Servet 1, CH-1206 Geneva, Switzerland
| | - Piero Procacci
- Chemistry Department, University of Florence, Via della Lastruccia 3-13, 50019 Sesto Fiorentino, Italy
| | - Francesco Luigi Gervasio
- School of Pharmaceutical Sciences, University of Geneva, Rue Michel-Servet 1, CH-1206 Geneva, Switzerland
- Institute of Pharmaceutical Sciences of Western Switzerland, University of Geneva, CH-1206 Geneva, Switzerland
- Chemistry Department, University College London (UCL), WC1E 6BT London, U.K
- Swiss Bioinformatics Institute, University of Geneva, CH-1206 Geneva, Switzerland
| |
Collapse
|
14
|
Vu NP, Ali L, Chua TL, Barr DA, Hendrickson HP, Trivedi DJ. Computational Insights into Prostaglandin E 2 Ligand Binding and Activation of G-Protein-Coupled Receptors. ACS APPLIED BIO MATERIALS 2024; 7:579-587. [PMID: 37058420 DOI: 10.1021/acsabm.2c01049] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/15/2023]
Abstract
G-protein coupled receptors (GPCRs) are eukaryotic integral membrane proteins that regulate signal transduction cascade pathways implicated in a variety of human diseases and are consequently of interest as drug targets. For this reason, it is of interest to investigate the way in which specific ligands bind and trigger conformational changes in the receptor during activation and how this in turn modulates intracellular signaling. In the present study, we investigate the way in which the ligand Prostaglandin E2 interacts with three GPCRs in the E-prostanoid family: EP1, EP2, and EP3. We examine information transfer pathways based on long-time scale molecular dynamics simulations using transfer entropy and betweenness centrality to measure the physical transfer of information among residues in the system. We monitor specific residues involved in binding to the ligand and investigate how the information transfer behavior of these residues changes upon ligand binding. Our results provide key insights that enable a deeper understanding of EP activation and signal transduction functioning pathways at the molecular level, as well as enabling us to make some predictions about the activation pathway for the EP1 receptor, for which little structural information is currently available. Our results should advance ongoing efforts in the development of potential therapeutics targeting these receptors.
Collapse
Affiliation(s)
- Nam P Vu
- Department of Chemistry, Lafayette College, Easton, Pennsylvania 18042, United States
| | - Luke Ali
- Department of Physics, Clarkson University, Potsdam, New York 13699, United States
| | - Theresa L Chua
- Department of Chemistry, Lafayette College, Easton, Pennsylvania 18042, United States
| | - Daniel A Barr
- Department of Chemistry, University of Mary, Bismarck, North Dakota 58504, United States
| | - Heidi P Hendrickson
- Department of Chemistry, Lafayette College, Easton, Pennsylvania 18042, United States
| | - Dhara J Trivedi
- Department of Physics, Clarkson University, Potsdam, New York 13699, United States
| |
Collapse
|
15
|
Khuttan S, Gallicchio E. What to Make of Zero: Resolving the Statistical Noise from Conformational Reorganization in Alchemical Binding Free Energy Estimates with Metadynamics Sampling. J Chem Theory Comput 2024; 20:1489-1501. [PMID: 38252868 PMCID: PMC10867849 DOI: 10.1021/acs.jctc.3c01250] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2023] [Revised: 01/03/2024] [Accepted: 01/03/2024] [Indexed: 01/24/2024]
Abstract
We introduce the self-relative binding free energy (self-RBFE) approach to evaluate the intrinsic statistical variance of dual-topology alchemical binding free energy estimators. The self-RBFE is the relative binding free energy between a ligand and a copy of the same ligand, and its true value is zero. Nevertheless, because the two copies of the ligand move independently, the self-RBFE value produced by a finite-length simulation fluctuates and can be used to measure the variance of the model. The results of this validation provide evidence that a significant fraction of the errors observed in benchmark studies reflect the statistical fluctuations of unconverged estimates rather than the models' accuracy. Furthermore, we find that ligand reorganization is a significant contributing factor to the statistical variance of binding free energy estimates and that metadynamics-accelerated conformational sampling of the torsional degrees of freedom of the ligand can drastically reduce the time to convergence.
Collapse
Affiliation(s)
- Sheenam Khuttan
- Department
of Chemistry and Biochemistry, Brooklyn
College of the City University of New York, New York, New York 11210, United States
- Ph.D.
Program in Biochemistry, The Graduate Center
of the City University of New York, New York, New York 10016, United States
| | - Emilio Gallicchio
- Department
of Chemistry and Biochemistry, Brooklyn
College of the City University of New York, New York, New York 11210, United States
- Ph.D.
Program in Biochemistry, The Graduate Center
of the City University of New York, New York, New York 10016, United States
- Ph.D.
Program in Chemistry, The Graduate Center
of the City University of New York, New York, New York 10016, United States
| |
Collapse
|
16
|
Georgiou K, Konstantinidi A, Hutterer J, Freudenberger K, Kolarov F, Lambrinidis G, Stylianakis I, Stampelou M, Gauglitz G, Kolocouris A. Accurate calculation of affinity changes to the close state of influenza A M2 transmembrane domain in response to subtle structural changes of adamantyl amines using free energy perturbation methods in different lipid bilayers. BIOCHIMICA ET BIOPHYSICA ACTA. BIOMEMBRANES 2024; 1866:184258. [PMID: 37995846 DOI: 10.1016/j.bbamem.2023.184258] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/06/2023] [Revised: 10/18/2023] [Accepted: 11/15/2023] [Indexed: 11/25/2023]
Abstract
Experimental binding free energies of 27 adamantyl amines against the influenza M2(22-46) WT tetramer, in its closed form at pH 8, were measured by ITC in DPC micelles. The measured Kd's range is ~44 while the antiviral potencies (IC50) range is ~750 with a good correlation between binding free energies computed with Kd and IC50 values (r = 0.76). We explored with MD simulations (ff19sb, CHARMM36m) the binding profile of complexes with strong, moderate and weak binders embedded in DMPC, DPPC, POPC or a viral mimetic membrane and using different experimental starting structures of M2. To predict accurately differences in binding free energy in response to subtle changes in the structure of the ligands, we performed 18 alchemical perturbative single topology FEP/MD NPT simulations (OPLS2005) using the BAR estimator (Desmond software) and 20 dual topology calculations TI/MD NVT simulations (ff19sb) using the MBAR estimator (Amber software) for adamantyl amines in complex with M2(22-46) WT in DMPC, DPPC, POPC. We observed that both methods with all lipids show a very good correlation between the experimental and calculated relative binding free energies (r = 0.77-0.87, mue = 0.36-0.92 kcal mol-1) with the highest performance achieved with TI/MBAR and lowest performance with FEP/BAR in DMPC bilayers. When antiviral potencies are used instead of the Kd values for computing the experimental binding free energies we obtained also good performance with both FEP/BAR (r = 0.83, mue = 0.75 kcal mol-1) and TI/MBAR (r = 0.69, mue = 0.77 kcal mol-1).
Collapse
Affiliation(s)
- Kyriakos Georgiou
- Laboratory of Medicinal Chemistry, Department of Pharmaceutical Chemistry, Faculty of Pharmacy, National and Kapodistrian University of Athens (NKUA), Panepistimiopolis-Zografou, 15771 Athens, Greece
| | - Athina Konstantinidi
- Laboratory of Medicinal Chemistry, Department of Pharmaceutical Chemistry, Faculty of Pharmacy, National and Kapodistrian University of Athens (NKUA), Panepistimiopolis-Zografou, 15771 Athens, Greece
| | - Johanna Hutterer
- Institut für Physikalische und Theoretische Chemie, Eberhard-Karls-Universität, D-72076 Tübingen, Germany
| | - Kathrin Freudenberger
- Institut für Physikalische und Theoretische Chemie, Eberhard-Karls-Universität, D-72076 Tübingen, Germany
| | - Felix Kolarov
- Institut für Physikalische und Theoretische Chemie, Eberhard-Karls-Universität, D-72076 Tübingen, Germany; Roche, Penzberg, Bavaria, Germany
| | - George Lambrinidis
- Laboratory of Medicinal Chemistry, Department of Pharmaceutical Chemistry, Faculty of Pharmacy, National and Kapodistrian University of Athens (NKUA), Panepistimiopolis-Zografou, 15771 Athens, Greece
| | - Ioannis Stylianakis
- Laboratory of Medicinal Chemistry, Department of Pharmaceutical Chemistry, Faculty of Pharmacy, National and Kapodistrian University of Athens (NKUA), Panepistimiopolis-Zografou, 15771 Athens, Greece
| | - Margarita Stampelou
- Laboratory of Medicinal Chemistry, Department of Pharmaceutical Chemistry, Faculty of Pharmacy, National and Kapodistrian University of Athens (NKUA), Panepistimiopolis-Zografou, 15771 Athens, Greece
| | - Günter Gauglitz
- Institut für Physikalische und Theoretische Chemie, Eberhard-Karls-Universität, D-72076 Tübingen, Germany
| | - Antonios Kolocouris
- Laboratory of Medicinal Chemistry, Department of Pharmaceutical Chemistry, Faculty of Pharmacy, National and Kapodistrian University of Athens (NKUA), Panepistimiopolis-Zografou, 15771 Athens, Greece.
| |
Collapse
|
17
|
Napolitano C, Benfatti F, Hamdan FB, Bristow JA, Dapiaggi F, Firth LC, Guest M, Saunders HA, Hall RG, Monaco MR, Quetglas V, Rendine S, Eterovic M. Synthesis and insecticidal activity of N-(5-phenylpyrazin-2-yl)-benzamide derivatives: Elucidation of mode of action on chitin biosynthesis through symptomology and genetic studies. PESTICIDE BIOCHEMISTRY AND PHYSIOLOGY 2024; 199:105771. [PMID: 38458679 DOI: 10.1016/j.pestbp.2024.105771] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/26/2023] [Revised: 12/17/2023] [Accepted: 01/08/2024] [Indexed: 03/10/2024]
Abstract
Among the six-membered heterocycles, the pyrazine ring is poorly explored in crop protection and does not feature in any product listed in the current IRAC MoA classification. In an effort to identify new leads for internal research, we synthesized a series of N-(5-phenylpyrazin-2-yl)-benzamide derivatives and evaluated them for their insecticidal activity. N-(5-phenylpyrazin-2-yl)-benzamide derivatives 3 were prepared using an automated two-step synthesis protocol. These compounds were tested for their initial biological activity against a wide range of sucking and chewing insect pests and found to be active against lepidopterans only. More detailed experiments, including symptomology studies on the diamondback moth, Plutella xylostella (L.) and the Egyptian cotton leafworm, Spodoptera littoralis (Boisduval) showed that analog 3q causes severe abnormalities in the lepidopteran cuticle leading to larval mortality. Compound 3q shows strong potency against both P. xylostella and S. littoralis, whereas analog 3i shows better potency against S. littoralis causing also impaired cuticular structure and death of the larvae. Additionally, P. xylostella genetic studies showed that compound 3q resistance is linked to Chitin Synthase 1. Our studies show that N-(5-phenylpyrazin-2-yl)-benzamide derivatives 3, and in particular analogs 3i and 3q, act as insect growth modulator insecticides. Conformational similarities with lufenuron are discussed.
Collapse
Affiliation(s)
- Carmela Napolitano
- Syngenta Crop Protection AG, Schaffhauserstrasse 101, Stein CH-4332, Switzerland
| | - Fides Benfatti
- Syngenta Crop Protection AG, Schaffhauserstrasse 101, Stein CH-4332, Switzerland
| | - Farhan Bou Hamdan
- Syngenta Crop Protection AG, Schaffhauserstrasse 101, Stein CH-4332, Switzerland
| | - Julia A Bristow
- Syngenta Crop Protection, Jealott's Hill International Research Centre, Bracknell, Berkshire RG42 6EY, UK
| | - Federico Dapiaggi
- Syngenta Crop Protection AG, Schaffhauserstrasse 101, Stein CH-4332, Switzerland
| | - Lucy C Firth
- Syngenta Crop Protection, Jealott's Hill International Research Centre, Bracknell, Berkshire RG42 6EY, UK
| | - Marcus Guest
- Syngenta Crop Protection, Jealott's Hill International Research Centre, Bracknell, Berkshire RG42 6EY, UK.
| | - Helena A Saunders
- Syngenta Crop Protection, Jealott's Hill International Research Centre, Bracknell, Berkshire RG42 6EY, UK
| | - Roger G Hall
- Syngenta Crop Protection AG, Schaffhauserstrasse 101, Stein CH-4332, Switzerland
| | - Mattia R Monaco
- Syngenta Crop Protection AG, Schaffhauserstrasse 101, Stein CH-4332, Switzerland
| | - Vincent Quetglas
- Syngenta Crop Protection AG, Schaffhauserstrasse 101, Stein CH-4332, Switzerland
| | - Stefano Rendine
- Syngenta Crop Protection AG, Schaffhauserstrasse 101, Stein CH-4332, Switzerland
| | - Marisa Eterovic
- Syngenta Crop Protection AG, Schaffhauserstrasse 101, Stein CH-4332, Switzerland.
| |
Collapse
|
18
|
Narendra G, Raju B, Verma H, Kumar M, Jain SK, Tung GK, Thakur S, Kaur R, Kaur S, Sapra B, Silakari O. Scaffold hopping based designing of selective ALDH1A1 inhibitors to overcome cyclophosphamide resistance: synthesis and biological evaluation. RSC Med Chem 2024; 15:309-321. [PMID: 38283216 PMCID: PMC10809718 DOI: 10.1039/d3md00543g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2023] [Accepted: 11/27/2023] [Indexed: 01/30/2024] Open
Abstract
Aldehyde dehydrogenase 1A1 (ALDH1A1) is an isoenzyme that catalyzes the conversion of aldehydes to acids. However, the overexpression of ALDH1A1 in a variety of malignancies is the major cause of resistance to an anti-cancer drug, cyclophosphamide (CP). CP is a prodrug that is initially converted into 4-hydroxycyclophosphamide and its tautomer aldophosphamide, in the liver. These compounds permeate into the cell and are converted as active metabolites, i.e., phosphoramide mustard (PM), through spontaneous beta-elimination. On the other hand, the conversion of CP to PM is diverted at the level of aldophosphamide by converting it into inactive carboxyphosphamide using ALDH1A1, which ultimately leads to high drug inactivation and CP resistance. Hence, in combination with our earlier work on the target of resistance, i.e., ALDH1A1, we hereby report selective ALDH1A1 inhibitors. Herein, we selected a lead molecule from our previous virtual screening and implemented scaffold hopping analysis to identify a novel scaffold that can act as an ALDH1A1 inhibitor. This results in the identification of various novel scaffolds. Among these, on the basis of synthetic feasibility, the benzimidazole scaffold was selected for the design of novel ALDH1A1 inhibitors, followed by machine learning-assisted structure-based virtual screening. Finally, the five best compounds were selected and synthesized. All synthesized compounds were evaluated using in vitro enzymatic assay against ALDH1A1, ALDH2, and ALDH3A1. The results disclosed that three molecules A1, A2, and A3 showed significant selective ALDH1A1 inhibitory potential with an IC50 value of 0.32 μM, 0.55 μM, and 1.63 μM, respectively, and none of the compounds exhibits potency towards the other two ALDH isoforms i.e. ALDH2 and ALDH3A1. Besides, the potent compounds (A1, A2, and A3) have been tested for in vitro cell line assay in combination with mafosfamide (analogue of CP) on two cell lines i.e. A549 and MIA-PaCa-2. All three compounds show significant potency to reverse mafosfamide resistance by inhibiting ALDH1A1 against these cell lines.
Collapse
Affiliation(s)
- Gera Narendra
- Department of Pharmaceutical Sciences and Drug Research, Punjabi University Patiala Punjab 147002 India +91 17522 83075 +91 95015 42696
| | - Baddipadige Raju
- Department of Pharmaceutical Sciences and Drug Research, Punjabi University Patiala Punjab 147002 India +91 17522 83075 +91 95015 42696
| | - Himanshu Verma
- Department of Pharmaceutical Sciences and Drug Research, Punjabi University Patiala Punjab 147002 India +91 17522 83075 +91 95015 42696
| | - Manoj Kumar
- Department of Pharmaceutical Sciences and Drug Research, Punjabi University Patiala Punjab 147002 India +91 17522 83075 +91 95015 42696
| | - Subheet Kumar Jain
- Department of Pharmaceutical Sciences, Guru Nanak Dev University Amritsar India
| | - Gurleen Kaur Tung
- Centre for Basic and Translational Research in Health Sciences, Guru Nanak Dev University Amritsar India
| | - Shubham Thakur
- Department of Pharmaceutical Sciences, Guru Nanak Dev University Amritsar India
| | - Rasdeep Kaur
- Department of Botany and Environmental Sciences, Guru Nanak Dev University Amritsar India
| | - Satwinderjeet Kaur
- Department of Botany and Environmental Sciences, Guru Nanak Dev University Amritsar India
| | - Bharti Sapra
- Department of Pharmaceutical Sciences and Drug Research, Punjabi University Patiala Punjab 147002 India +91 17522 83075 +91 95015 42696
| | - Om Silakari
- Department of Pharmaceutical Sciences and Drug Research, Punjabi University Patiala Punjab 147002 India +91 17522 83075 +91 95015 42696
| |
Collapse
|
19
|
Chen L, Wu Y, Wu C, Silveira A, Sherman W, Xu H, Gallicchio E. Performance and Analysis of the Alchemical Transfer Method for Binding-Free-Energy Predictions of Diverse Ligands. J Chem Inf Model 2024; 64:250-264. [PMID: 38147877 DOI: 10.1021/acs.jcim.3c01705] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2023]
Abstract
The Alchemical Transfer Method (ATM) is herein validated against the relative binding-free energies (RBFEs) of a diverse set of protein-ligand complexes. We employed a streamlined setup workflow, a bespoke force field, and AToM-OpenMM software to compute the RBFEs of the benchmark set prepared by Schindler and collaborators at Merck KGaA. This benchmark set includes examples of standard small R-group ligand modifications as well as more challenging scenarios, such as large R-group changes, scaffold hopping, formal charge changes, and charge-shifting transformations. The novel coordinate perturbation scheme and a dual-topology approach of ATM address some of the challenges of single-topology alchemical RBFE methods. Specifically, ATM eliminates the need for splitting electrostatic and Lennard-Jones interactions, atom mapping, defining ligand regions, and postcorrections for charge-changing perturbations. Thus, ATM is simpler and more broadly applicable than conventional alchemical methods, especially for scaffold-hopping and charge-changing transformations. Here, we performed well over 500 RBFE calculations for eight protein targets and found that ATM achieves accuracy comparable to that of existing state-of-the-art methods, albeit with larger statistical fluctuations. We discuss insights into the specific strengths and weaknesses of the ATM method that will inform future deployments. This study confirms that ATM can be applied as a production tool for RBFE predictions across a wide range of perturbation types within a unified, open-source framework.
Collapse
Affiliation(s)
- Lieyang Chen
- Roivant Sciences, 151 W 42nd Street, 15th Floor, New York, New York 10036, United States
| | - Yujie Wu
- Roivant Sciences, 151 W 42nd Street, 15th Floor, New York, New York 10036, United States
- Atommap Corporation, New York, New York 10017, United States
| | - Chuanjie Wu
- Roivant Sciences, 151 W 42nd Street, 15th Floor, New York, New York 10036, United States
| | - Ana Silveira
- Psivant Therapeutics, 451 D Street, Boston, Massachusetts 02210, United States
| | - Woody Sherman
- Psivant Therapeutics, 451 D Street, Boston, Massachusetts 02210, United States
| | - Huafeng Xu
- Roivant Sciences, 151 W 42nd Street, 15th Floor, New York, New York 10036, United States
- Atommap Corporation, New York, New York 10017, United States
| | - Emilio Gallicchio
- Department of Chemistry and Biochemistry, Brooklyn College of the City University of New York, New York, New York 11210, United States
- Ph.D. Program in Chemistry, The Graduate Center of the City University of New York, New York, New York 10016, United States
- Ph.D. Program in Biochemistry, The Graduate Center of the City University of New York, New York, New York 10016, United States
| |
Collapse
|
20
|
Valdivia A, Luque FJ, Llabrés S. Binding of Cholesterol to the N-Terminal Domain of the NPC1L1 Transporter: Analysis of the Epimerization-Related Binding Selectivity and Loop Mutations. J Chem Inf Model 2024; 64:189-204. [PMID: 38152929 PMCID: PMC10777396 DOI: 10.1021/acs.jcim.3c01319] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2023] [Revised: 12/14/2023] [Accepted: 12/15/2023] [Indexed: 12/29/2023]
Abstract
Cholesterol is a fat-like substance with a pivotal physiological relevance in humans, and its homeostasis is tightly regulated by various cellular processes, including the import in the small intestine and the reabsorption in the biliary ducts by the Niemann-Pick C1 Like 1 (NPC1L1) importer. NPC1L1 can mediate the absorption of a variety of sterols but strikingly exhibits a large sensitivity to cholesterol epimerization. This study examines the molecular basis of the epimerization-related selective binding of cholesterol by combining extended unbiased molecular dynamics simulations of the apo and holo species of the N-terminal domain of wild-type NPC1L1, in conjunction with relative binding free energy, umbrella sampling, and well-tempered metadynamics calculations. The analysis of the results discloses the existence of two distinct binding modes for cholesterol and epi-cholesterol. The former binds deeper in the cavity, forming key hydrogen-bond interactions with Q95, S56, and a water molecule. In contrast, epi-cholesterol is shifted ca. 3 Å to the mouth of the cavity and the transition to the Q95 site is prevented by an energetic barrier of 4.1 kcal·mol-1. Thus, the configuration of the hydroxyl group of cholesterol, together with the presence of a structural water molecule, is a key feature for effective absorption. Finally, whereas these findings may seemingly be challenged by single-point mutations that impair cholesterol transport but have a mild impact on the binding of cholesterol to the Q95 binding site, our results reveal that they have a drastic influence on the conformational landscape of the α8/β7 loop in the apo species compared to the wild-type protein. Overall, the results give support to the functional role played by the α8/β7 loop in regulating the access of ligands to NPC1L1, and hence to interpreting the impact of these mutations on diseases related to disruption of sterol absorption, paving the way to understanding certain physiological dysfunctions.
Collapse
Affiliation(s)
- Aitor Valdivia
- Departament
de Nutrició, Ciències de l′Alimentació
i Gastronomia, Facultat de Farmàcia
i Ciències de l′Alimentació—Campus Torribera,
Universitat de Barcelona, Prat de la Riba 171, 08921 Santa Coloma de Gramenet, Spain
- Institut
de Biomedicina (IBUB), Universitat de Barcelona, 08028 Barcelona, Spain
| | - F. Javier Luque
- Departament
de Nutrició, Ciències de l′Alimentació
i Gastronomia, Facultat de Farmàcia
i Ciències de l′Alimentació—Campus Torribera,
Universitat de Barcelona, Prat de la Riba 171, 08921 Santa Coloma de Gramenet, Spain
- Institut
de Biomedicina (IBUB), Universitat de Barcelona, 08028 Barcelona, Spain
- Institut
de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, 08921 Barcelona, Spain
| | - Salomé Llabrés
- Departament
de Nutrició, Ciències de l′Alimentació
i Gastronomia, Facultat de Farmàcia
i Ciències de l′Alimentació—Campus Torribera,
Universitat de Barcelona, Prat de la Riba 171, 08921 Santa Coloma de Gramenet, Spain
- Institut
de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona, 08921 Barcelona, Spain
| |
Collapse
|
21
|
Melancon K, Pliushcheuskaya P, Meiler J, Künze G. Targeting ion channels with ultra-large library screening for hit discovery. Front Mol Neurosci 2024; 16:1336004. [PMID: 38249296 PMCID: PMC10796734 DOI: 10.3389/fnmol.2023.1336004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2023] [Accepted: 12/05/2023] [Indexed: 01/23/2024] Open
Abstract
Ion channels play a crucial role in a variety of physiological and pathological processes, making them attractive targets for drug development in diseases such as diabetes, epilepsy, hypertension, cancer, and chronic pain. Despite the importance of ion channels in drug discovery, the vastness of chemical space and the complexity of ion channels pose significant challenges for identifying drug candidates. The use of in silico methods in drug discovery has dramatically reduced the time and cost of drug development and has the potential to revolutionize the field of medicine. Recent advances in computer hardware and software have enabled the screening of ultra-large compound libraries. Integration of different methods at various scales and dimensions is becoming an inevitable trend in drug development. In this review, we provide an overview of current state-of-the-art computational chemistry methodologies for ultra-large compound library screening and their application to ion channel drug discovery research. We discuss the advantages and limitations of various in silico techniques, including virtual screening, molecular mechanics/dynamics simulations, and machine learning-based approaches. We also highlight several successful applications of computational chemistry methodologies in ion channel drug discovery and provide insights into future directions and challenges in this field.
Collapse
Affiliation(s)
- Kortney Melancon
- Department of Chemistry, Vanderbilt University, Nashville, TN, United States
- Center for Structural Biology, Vanderbilt University, Nashville, TN, United States
| | | | - Jens Meiler
- Department of Chemistry, Vanderbilt University, Nashville, TN, United States
- Center for Structural Biology, Vanderbilt University, Nashville, TN, United States
- Medical Faculty, Institute for Drug Discovery, Leipzig University, Leipzig, Germany
- Center for Scalable Data Analytics and Artificial Intelligence, Leipzig University, Leipzig, Germany
| | - Georg Künze
- Medical Faculty, Institute for Drug Discovery, Leipzig University, Leipzig, Germany
- Center for Scalable Data Analytics and Artificial Intelligence, Leipzig University, Leipzig, Germany
- Interdisciplinary Center for Bioinformatics, Leipzig University, Leipzig, Germany
| |
Collapse
|
22
|
Herz AM, Kellici T, Morao I, Michel J. Alchemical Free Energy Workflows for the Computation of Protein-Ligand Binding Affinities. Methods Mol Biol 2024; 2716:241-264. [PMID: 37702943 DOI: 10.1007/978-1-0716-3449-3_11] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/14/2023]
Abstract
Alchemical free energy methods can be used for the efficient computation of relative binding free energies during preclinical drug discovery stages. In recent years, this has been facilitated further by the implementation of workflows that enable non-experts to quickly and consistently set up the required simulations. Given the correct input structures, workflows handle the difficult aspects of setting up perturbations, including consistently defining the perturbable molecule, its atom mapping and topology generation, perturbation network generation, running of the simulations via different sampling methods, and analysis of the results. Different academic and commercial workflows are discussed, including FEW, FESetup, FEPrepare, CHARMM-GUI, Transformato, PMX, QLigFEP, TIES, ProFESSA, PyAutoFEP, BioSimSpace, FEP+, Flare, and Orion. These workflows differ in various aspects, such as mapping algorithms or enhanced sampling methods. Some workflows can accommodate more than one molecular dynamics (MD) engine and use external libraries for tasks. Differences between workflows can present advantages for different use cases, however a lack of interoperability of the workflows' components hinders systematic comparisons.
Collapse
Affiliation(s)
- Anna M Herz
- EaStChem School of Chemistry, Joseph Black Building, University of Edinburgh, Edinburgh, UK
| | - Tahsin Kellici
- Evotec (UK) Ltd., In Silico Research and Development, Abingdon, Oxfordshire, UK
- Merck & Co., Inc., Modelling and Informatics, West Point, PA, USA
| | - Inaki Morao
- Evotec (UK) Ltd., In Silico Research and Development, Abingdon, Oxfordshire, UK
| | - Julien Michel
- EaStChem School of Chemistry, Joseph Black Building, University of Edinburgh, Edinburgh, UK.
| |
Collapse
|
23
|
Robo MT, Hayes RL, Ding X, Pulawski B, Vilseck JZ. Fast free energy estimates from λ-dynamics with bias-updated Gibbs sampling. Nat Commun 2023; 14:8515. [PMID: 38129400 PMCID: PMC10740020 DOI: 10.1038/s41467-023-44208-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/12/2022] [Accepted: 12/04/2023] [Indexed: 12/23/2023] Open
Abstract
Relative binding free energy calculations have become an integral computational tool for lead optimization in structure-based drug design. Classical alchemical methods, including free energy perturbation or thermodynamic integration, compute relative free energy differences by transforming one molecule into another. However, these methods have high operational costs due to the need to perform many pairwise perturbations independently. To reduce costs and accelerate molecular design workflows, we present a method called λ-dynamics with bias-updated Gibbs sampling. This method uses dynamic biases to continuously sample between multiple ligand analogues collectively within a single simulation. We show that many relative binding free energies can be determined quickly with this approach without compromising accuracy. For five benchmark systems, agreement to experiment is high, with root mean square errors near or below 1.0 kcal mol-1. Free energy results are consistent with other computational approaches and within statistical noise of both methods (0.4 kcal mol-1 or less). Notably, large efficiency gains over thermodynamic integration of 18-66-fold for small perturbations and 100-200-fold for whole aromatic ring substitutions are observed. The rapid determination of relative binding free energies will enable larger chemical spaces to be more readily explored and structure-based drug design to be accelerated.
Collapse
Affiliation(s)
- Michael T Robo
- Department of Biochemistry and Molecular Biology, Indiana University School of Medicine, Indianapolis, IN, 46202, USA
- Center for Computational Biology and Bioinformatics, Indiana University School of Medicine, Indianapolis, IN, 46202, USA
- Indiana Biosciences Research Institute, 1210 Waterway Blvd Ste. 2000, Indianapolis, IN, 46202, USA
| | - Ryan L Hayes
- Chemical and Biomolecular Engineering, University of California, Irvine, California, 92617, USA
- Pharmaceutical Sciences, University of California, Irvine, CA, 92617, USA
| | - Xinqiang Ding
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
- Department of Chemistry, Tufts University, Medford, MA, 02144, USA
| | - Brian Pulawski
- Department of Biochemistry and Molecular Biology, Indiana University School of Medicine, Indianapolis, IN, 46202, USA
| | - Jonah Z Vilseck
- Department of Biochemistry and Molecular Biology, Indiana University School of Medicine, Indianapolis, IN, 46202, USA.
- Center for Computational Biology and Bioinformatics, Indiana University School of Medicine, Indianapolis, IN, 46202, USA.
| |
Collapse
|
24
|
Fu H, Chipot C, Shao X, Cai W. Standard Binding Free-Energy Calculations: How Far Are We from Automation? J Phys Chem B 2023; 127:10459-10468. [PMID: 37824848 DOI: 10.1021/acs.jpcb.3c04370] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2023]
Abstract
Recent success stories suggest that in silico protein-ligand binding free-energy calculations are approaching chemical accuracy. However, their widespread application remains limited by the extensive human intervention required, posing challenges for the neophyte. As such, it is critical to develop automated workflows for estimating protein-ligand binding affinities with minimum personal involvement. Key human efforts include setting up and tuning enhanced-sampling or alchemical-transformation algorithms as a preamble to computational binding free-energy estimations. Additionally, preparing input files, bookkeeping, and postprocessing represent nontrivial tasks. In this Perspective, we discuss recent progress in automating standard binding free-energy calculations, featuring the development of adaptive or parameter-free algorithms, standardization of binding free-energy calculation workflows, and the implementation of user-friendly software. We also assess the current state of automated standard binding free-energy calculations and evaluate the limitations of existing methods. Last, we outline the requirements for future algorithms and workflows to facilitate automated free-energy calculations for diverse protein-ligand complexes.
Collapse
Affiliation(s)
- Haohao Fu
- State Key Laboratory of Medicinal Chemical Biology, Tianjin Key Laboratory of Biosensing and Molecular Recognition, Research Center for Analytical Sciences, Frontiers Science Center for New Organic Matter, College of Chemistry, Nankai University, Tianjin 300071, China
- Haihe Laboratory of Sustainable Chemical Transformations, Tianjin 300192, China
| | - Christophe Chipot
- Laboratoire International Associé CNRS and University of Illinois at Urbana-Champaign, UMR no. 7019, Université de Lorraine, BP 70239, F-54506 Vandoeuvre-lès-Nancy, France
- Department of Physics, University of Illinois at Urbana-Champaign, 1110 West Green Street, Urbana, Illinois 61801, United States
- Department of Chemistry, The University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637, United States
- Department of Chemistry, The University of Hawai'i at Ma̅noa, 2545 McCarthy Mall, Honolulu, Hawaii 96822, United States
| | - Xueguang Shao
- State Key Laboratory of Medicinal Chemical Biology, Tianjin Key Laboratory of Biosensing and Molecular Recognition, Research Center for Analytical Sciences, Frontiers Science Center for New Organic Matter, College of Chemistry, Nankai University, Tianjin 300071, China
- Haihe Laboratory of Sustainable Chemical Transformations, Tianjin 300192, China
| | - Wensheng Cai
- State Key Laboratory of Medicinal Chemical Biology, Tianjin Key Laboratory of Biosensing and Molecular Recognition, Research Center for Analytical Sciences, Frontiers Science Center for New Organic Matter, College of Chemistry, Nankai University, Tianjin 300071, China
- Haihe Laboratory of Sustainable Chemical Transformations, Tianjin 300192, China
| |
Collapse
|
25
|
Champion C, Gall R, Ries B, Rieder SR, Barros EP, Riniker S. Accelerating Alchemical Free Energy Prediction Using a Multistate Method: Application to Multiple Kinases. J Chem Inf Model 2023; 63:7133-7147. [PMID: 37948537 PMCID: PMC10685456 DOI: 10.1021/acs.jcim.3c01469] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2023] [Revised: 10/23/2023] [Accepted: 10/23/2023] [Indexed: 11/12/2023]
Abstract
Alchemical free-energy methods based on molecular dynamics (MD) simulations have become important tools to identify modifications of small organic molecules that improve their protein binding affinity during lead optimization. The routine application of pairwise free-energy methods to rank potential binders from best to worst is impacted by the combinatorial increase in calculations to perform when the number of molecules to assess grows. To address this fundamental limitation, our group has developed replica-exchange enveloping distribution sampling (RE-EDS), a pathway-independent multistate method, enabling the calculation of alchemical free-energy differences between multiple ligands (N > 2) from a single MD simulation. In this work, we apply the method to a set of four kinases with diverse binding pockets and their corresponding inhibitors (42 in total), chosen to showcase the general applicability of RE-EDS in prospective drug design campaigns. We show that for the targets studied, RE-EDS is able to model up to 13 ligands simultaneously with high sampling efficiency, leading to a substantial decrease in computational cost when compared to pairwise methods.
Collapse
Affiliation(s)
- Candide Champion
- Department of Chemistry and
Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - René Gall
- Department of Chemistry and
Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | | | - Salomé R. Rieder
- Department of Chemistry and
Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Emilia P. Barros
- Department of Chemistry and
Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Sereina Riniker
- Department of Chemistry and
Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
26
|
Schuurs ZP, Martyn AP, Soltau CP, Beard S, Shah ET, Adams MN, Croft LV, O’Byrne KJ, Richard DJ, Gandhi NS. An Exploration of Small Molecules That Bind Human Single-Stranded DNA Binding Protein 1. BIOLOGY 2023; 12:1405. [PMID: 37998004 PMCID: PMC10669474 DOI: 10.3390/biology12111405] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/27/2023] [Revised: 10/24/2023] [Accepted: 10/30/2023] [Indexed: 11/25/2023]
Abstract
Human single-stranded DNA binding protein 1 (hSSB1) is critical to preserving genome stability, interacting with single-stranded DNA (ssDNA) through an oligonucleotide/oligosaccharide binding-fold. The depletion of hSSB1 in cell-line models leads to aberrant DNA repair and increased sensitivity to irradiation. hSSB1 is over-expressed in several types of cancers, suggesting that hSSB1 could be a novel therapeutic target in malignant disease. hSSB1 binding studies have focused on DNA; however, despite the availability of 3D structures, small molecules targeting hSSB1 have not been explored. Quinoline derivatives targeting hSSB1 were designed through a virtual fragment-based screening process, synthesizing them using AlphaLISA and EMSA to determine their affinity for hSSB1. In parallel, we further screened a structurally diverse compound library against hSSB1 using the same biochemical assays. Three compounds with nanomolar affinity for hSSB1 were identified, exhibiting cytotoxicity in an osteosarcoma cell line. To our knowledge, this is the first study to identify small molecules that modulate hSSB1 activity. Molecular dynamics simulations indicated that three of the compounds that were tested bound to the ssDNA-binding site of hSSB1, providing a framework for the further elucidation of inhibition mechanisms. These data suggest that small molecules can disrupt the interaction between hSSB1 and ssDNA, and may also affect the ability of cells to repair DNA damage. This test study of small molecules holds the potential to provide insights into fundamental biochemical questions regarding the OB-fold.
Collapse
Affiliation(s)
- Zachariah P. Schuurs
- Centre for Genomics and Personalised Health, School of Chemistry and Physics, Queensland University of Technology (QUT), Brisbane, QLD 4000, Australia; (Z.P.S.); (A.P.M.)
- Cancer and Ageing Research Program, Centre for Genomics and Personalised Health, Queensland University of Technology (QUT), Translational Research Institute (TRI), Woolloongabba, QLD 4102, Australia; (S.B.); (M.N.A.); (L.V.C.); (K.J.O.); (D.J.R.)
| | - Alexander P. Martyn
- Centre for Genomics and Personalised Health, School of Chemistry and Physics, Queensland University of Technology (QUT), Brisbane, QLD 4000, Australia; (Z.P.S.); (A.P.M.)
- Cancer and Ageing Research Program, Centre for Genomics and Personalised Health, Queensland University of Technology (QUT), Translational Research Institute (TRI), Woolloongabba, QLD 4102, Australia; (S.B.); (M.N.A.); (L.V.C.); (K.J.O.); (D.J.R.)
| | - Carl P. Soltau
- School of Chemistry and Physics, Centre for Materials Science, Queensland University of Technology, Brisbane, QLD 4000, Australia;
| | - Sam Beard
- Cancer and Ageing Research Program, Centre for Genomics and Personalised Health, Queensland University of Technology (QUT), Translational Research Institute (TRI), Woolloongabba, QLD 4102, Australia; (S.B.); (M.N.A.); (L.V.C.); (K.J.O.); (D.J.R.)
| | - Esha T. Shah
- Cancer and Ageing Research Program, Centre for Genomics and Personalised Health, Queensland University of Technology (QUT), Translational Research Institute (TRI), Woolloongabba, QLD 4102, Australia; (S.B.); (M.N.A.); (L.V.C.); (K.J.O.); (D.J.R.)
- School of Biomedical Sciences, Faculty of Health, Queensland University of Technology, Translational Research Institute, Woolloongabba, QLD 4102, Australia
| | - Mark N. Adams
- Cancer and Ageing Research Program, Centre for Genomics and Personalised Health, Queensland University of Technology (QUT), Translational Research Institute (TRI), Woolloongabba, QLD 4102, Australia; (S.B.); (M.N.A.); (L.V.C.); (K.J.O.); (D.J.R.)
- School of Biomedical Sciences, Faculty of Health, Queensland University of Technology, Translational Research Institute, Woolloongabba, QLD 4102, Australia
| | - Laura V. Croft
- Cancer and Ageing Research Program, Centre for Genomics and Personalised Health, Queensland University of Technology (QUT), Translational Research Institute (TRI), Woolloongabba, QLD 4102, Australia; (S.B.); (M.N.A.); (L.V.C.); (K.J.O.); (D.J.R.)
- School of Biomedical Sciences, Faculty of Health, Queensland University of Technology, Translational Research Institute, Woolloongabba, QLD 4102, Australia
| | - Kenneth J. O’Byrne
- Cancer and Ageing Research Program, Centre for Genomics and Personalised Health, Queensland University of Technology (QUT), Translational Research Institute (TRI), Woolloongabba, QLD 4102, Australia; (S.B.); (M.N.A.); (L.V.C.); (K.J.O.); (D.J.R.)
- Cancer Services, Princess Alexandra Hospital—Metro South Health, Woolloongabba, QLD 4102, Australia
| | - Derek J. Richard
- Cancer and Ageing Research Program, Centre for Genomics and Personalised Health, Queensland University of Technology (QUT), Translational Research Institute (TRI), Woolloongabba, QLD 4102, Australia; (S.B.); (M.N.A.); (L.V.C.); (K.J.O.); (D.J.R.)
- School of Biomedical Sciences, Faculty of Health, Queensland University of Technology, Translational Research Institute, Woolloongabba, QLD 4102, Australia
| | - Neha S. Gandhi
- Centre for Genomics and Personalised Health, School of Chemistry and Physics, Queensland University of Technology (QUT), Brisbane, QLD 4000, Australia; (Z.P.S.); (A.P.M.)
- Cancer and Ageing Research Program, Centre for Genomics and Personalised Health, Queensland University of Technology (QUT), Translational Research Institute (TRI), Woolloongabba, QLD 4102, Australia; (S.B.); (M.N.A.); (L.V.C.); (K.J.O.); (D.J.R.)
- Department of Computer Science and Engineering, Manipal Institute of Technology, Manipal Academy of Higher Education, Manipal 576104, India
| |
Collapse
|
27
|
Braconi L, Dei S, Contino M, Riganti C, Bartolucci G, Manetti D, Romanelli MN, Perrone MG, Colabufo NA, Guglielmo S, Teodori E. Tetrazole and oxadiazole derivatives as bioisosteres of tariquidar and elacridar: New potent P-gp modulators acting as MDR reversers. Eur J Med Chem 2023; 259:115716. [PMID: 37573829 DOI: 10.1016/j.ejmech.2023.115716] [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: 05/11/2023] [Revised: 06/22/2023] [Accepted: 08/07/2023] [Indexed: 08/15/2023]
Abstract
New 2,5- and 1,5-disubstituted tetrazoles, and 2,5-disubstituted-1,3,4-oxadiazoles were synthesized as tariquidar and elacridar derivatives and studied as multidrug resistance (MDR) reversers. Their behaviour on the three ABC transporters P-gp, MRP1 and BCRP was investigated. All compounds inhibited the P-gp transport activity in MDCK-MDR1 cells overexpressing P-gp, showing EC50 values even in the low nanomolar range (compounds 15, 22). Oxadiazole derivatives were able to increase the antiproliferative effect of doxorubicin in MDCK-MDR1 and in HT29/DX cells confirming their nature of P-gp modulators, with derivative 15 being the most potent in these assays. Compound 15 also displayed a dual inhibitory effect showing good activities towards both P-gp and BCRP. A computational study suggested a common interaction pattern on P-gp for most of the potent compounds. The bioisosteric substitution of the amide group of lead compounds allowed identifying a new set of potent oxadiazole derivatives that modulate MDR through inhibition of the P-gp efflux activity. If compared to previous amide derivatives, the introduction of the heterocycle rings greatly enhances the activity on P-gp, introduces in two compounds a moderate inhibitory activity on MRP1 and maintains in some cases the effect on BCRP, leading to the unveiling of dual inhibitor 15.
Collapse
Affiliation(s)
- Laura Braconi
- Department of Neuroscience, Psychology, Drug Research and Child Health - Section of Pharmaceutical and Nutraceutical Sciences, University of Florence, Via Ugo Schiff 6, 50019, Sesto Fiorentino, Italy
| | - Silvia Dei
- Department of Neuroscience, Psychology, Drug Research and Child Health - Section of Pharmaceutical and Nutraceutical Sciences, University of Florence, Via Ugo Schiff 6, 50019, Sesto Fiorentino, Italy.
| | - Marialessandra Contino
- Department of Pharmacy - Drug Sciences, University of Bari "A. Moro", via Orabona 4, 70125, Bari, Italy
| | - Chiara Riganti
- Department of Oncology, University of Turin, Via Santena 5/bis, 10126, Torino, Italy
| | - Gianluca Bartolucci
- Department of Neuroscience, Psychology, Drug Research and Child Health - Section of Pharmaceutical and Nutraceutical Sciences, University of Florence, Via Ugo Schiff 6, 50019, Sesto Fiorentino, Italy
| | - Dina Manetti
- Department of Neuroscience, Psychology, Drug Research and Child Health - Section of Pharmaceutical and Nutraceutical Sciences, University of Florence, Via Ugo Schiff 6, 50019, Sesto Fiorentino, Italy
| | - Maria Novella Romanelli
- Department of Neuroscience, Psychology, Drug Research and Child Health - Section of Pharmaceutical and Nutraceutical Sciences, University of Florence, Via Ugo Schiff 6, 50019, Sesto Fiorentino, Italy
| | - Maria Grazia Perrone
- Department of Pharmacy - Drug Sciences, University of Bari "A. Moro", via Orabona 4, 70125, Bari, Italy
| | - Nicola Antonio Colabufo
- Department of Pharmacy - Drug Sciences, University of Bari "A. Moro", via Orabona 4, 70125, Bari, Italy
| | - Stefano Guglielmo
- Department of Drug Science and Technology, University of Turin, Via P. Giuria 9, 10125, Torino, Italy
| | - Elisabetta Teodori
- Department of Neuroscience, Psychology, Drug Research and Child Health - Section of Pharmaceutical and Nutraceutical Sciences, University of Florence, Via Ugo Schiff 6, 50019, Sesto Fiorentino, Italy
| |
Collapse
|
28
|
Schindler CEM, Kuhn D, Hartung IV. The experiment is the limit. Nat Rev Chem 2023; 7:752-753. [PMID: 37880428 DOI: 10.1038/s41570-023-00552-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2023]
Affiliation(s)
| | - Daniel Kuhn
- Discovery & Development Technologies, Merck KGaA, Darmstadt, Germany
| | - Ingo V Hartung
- Discovery & Development Technologies, Merck KGaA, Darmstadt, Germany
| |
Collapse
|
29
|
Ross GA, Lu C, Scarabelli G, Albanese SK, Houang E, Abel R, Harder ED, Wang L. The maximal and current accuracy of rigorous protein-ligand binding free energy calculations. Commun Chem 2023; 6:222. [PMID: 37838760 PMCID: PMC10576784 DOI: 10.1038/s42004-023-01019-9] [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: 10/18/2022] [Accepted: 10/02/2023] [Indexed: 10/16/2023] Open
Abstract
Computational techniques can speed up the identification of hits and accelerate the development of candidate molecules for drug discovery. Among techniques for predicting relative binding affinities, the most consistently accurate is free energy perturbation (FEP), a class of rigorous physics-based methods. However, uncertainty remains about how accurate FEP is and can ever be. Here, we present what we believe to be the largest publicly available dataset of proteins and congeneric series of small molecules, and assess the accuracy of the leading FEP workflow. To ascertain the limit of achievable accuracy, we also survey the reproducibility of experimental relative affinity measurements. We find a wide variability in experimental accuracy and a correspondence between binding and functional assays. When careful preparation of protein and ligand structures is undertaken, FEP can achieve accuracy comparable to experimental reproducibility. Throughout, we highlight reliable protocols that can help maximize the accuracy of FEP in prospective studies.
Collapse
Affiliation(s)
- Gregory A Ross
- Schrödinger Inc, New York, NY, USA.
- Isomorphic Labs, London, UK.
| | - Chao Lu
- Schrödinger Inc, New York, NY, USA
| | | | | | | | | | | | | |
Collapse
|
30
|
Rathje OH, Perryman L, Payne RJ, Hamprecht DW. PROTACs Targeting MLKL Protect Cells from Necroptosis. J Med Chem 2023; 66:11216-11236. [PMID: 37535857 DOI: 10.1021/acs.jmedchem.3c00665] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/05/2023]
Abstract
Mixed Lineage Kinase domain-Like pseudokinase (MLKL) is implicated in a broad range of diseases due to its role as the ultimate effector of necroptosis and has therefore emerged as an attractive drug target. Here, we describe the development of PROteolysis TArgeting Chimeras (PROTACs) as a novel approach to knock down MLKL through chemical means. A series of candidate degraders were synthesized from a high-affinity pyrazole carboxamide-based MLKL ligand leading to the identification of a PROTAC molecule that effectively degraded MLKL and completely abrogated cell death in a TSZ model of necroptosis. By leveraging the innate ability of these PROTACs to degrade MLKL in a dose-dependent manner, the quantitative relationship between MLKL levels and necroptosis was interrogated. This work demonstrates the feasibility of targeting MLKL using a PROTAC approach and provides a powerful tool to further our understanding of the role of MLKL within the necroptotic pathway.
Collapse
Affiliation(s)
- Oliver H Rathje
- Pharmaxis Ltd., 20 Rodborough Road, Frenchs Forest, NSW 2086, Australia
- School of Chemistry, The University of Sydney, Sydney, New South Wales 2006, Australia
- Australian Research Council Centre of Excellence for Innovations in Peptide and Protein Science, The University of Sydney, Sydney, New South Wales 2006, Australia
| | - Lara Perryman
- Pharmaxis Ltd., 20 Rodborough Road, Frenchs Forest, NSW 2086, Australia
| | - Richard J Payne
- School of Chemistry, The University of Sydney, Sydney, New South Wales 2006, Australia
- Australian Research Council Centre of Excellence for Innovations in Peptide and Protein Science, The University of Sydney, Sydney, New South Wales 2006, Australia
| | | |
Collapse
|
31
|
Zhang D, Qiao L, Lei X, Dong X, Tong Y, Wang J, Wang Z, Zhou R. Mutagenesis and structural studies reveal the basis for the specific binding of SARS-CoV-2 SL3 RNA element with human TIA1 protein. Nat Commun 2023; 14:3715. [PMID: 37349329 PMCID: PMC10287707 DOI: 10.1038/s41467-023-39410-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2022] [Accepted: 06/12/2023] [Indexed: 06/24/2023] Open
Abstract
Viral RNA-host protein interactions are indispensable during RNA virus transcription and replication, but their detailed structural and dynamical features remain largely elusive. Here, we characterize the binding interface for the SARS-CoV-2 stem-loop 3 (SL3) cis-acting element to human TIA1 protein with a combined theoretical and experimental approaches. The highly structured SARS-CoV-2 SL3 has a high binding affinity to TIA1 protein, in which the aromatic stacking, hydrogen bonds, and hydrophobic interactions collectively direct this specific binding. Further mutagenesis studies validate our proposed 3D binding model and reveal two SL3 variants have enhanced binding affinities to TIA1. And disruptions of the identified RNA-protein interactions with designed antisense oligonucleotides dramatically reduce SARS-CoV-2 infection in cells. Finally, TIA1 protein could interact with conserved SL3 RNA elements within other betacoronavirus lineages. These findings open an avenue to explore the viral RNA-host protein interactions and provide a pioneering structural basis for RNA-targeting antiviral drug design.
Collapse
Affiliation(s)
- Dong Zhang
- Institute of Quantitative Biology, College of Life Sciences, Zhejiang University, Hangzhou, Zhejiang, 310058, China
| | - Lulu Qiao
- State Key Laboratory of Plant Physiology and Biochemistry, College of Life Sciences, Zhejiang University, Hangzhou, Zhejiang, 310058, China
| | - Xiaobo Lei
- NHC Key Laboratory of Systems Biology of Pathogens and Christophe Mérieux Laboratory, Institute of Pathogen Biology, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing, 100730, China
| | - Xiaojing Dong
- NHC Key Laboratory of Systems Biology of Pathogens and Christophe Mérieux Laboratory, Institute of Pathogen Biology, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing, 100730, China
| | - Yunguang Tong
- College of Life Sciences, China Jiliang University, Hangzhou, Zhejiang, 310018, China
- Department of Pharmacy, China Jiliang University, Hangzhou, Zhejiang, 310018, China
| | - Jianwei Wang
- NHC Key Laboratory of Systems Biology of Pathogens and Christophe Mérieux Laboratory, Institute of Pathogen Biology, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing, 100730, China.
| | - Zhiye Wang
- State Key Laboratory of Plant Physiology and Biochemistry, College of Life Sciences, Zhejiang University, Hangzhou, Zhejiang, 310058, China.
- The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou, Zhejiang, 310058, China.
| | - Ruhong Zhou
- Institute of Quantitative Biology, College of Life Sciences, Zhejiang University, Hangzhou, Zhejiang, 310058, China.
- The First Affiliated Hospital, College of Medicine, Zhejiang University, Hangzhou, Zhejiang, 310058, China.
| |
Collapse
|
32
|
Bastos RS, de Lima LR, Neto MFA, Yousaf N, Cruz JN, Campos JM, Kimani NM, Ramos RS, Santos CBR. Design and Identification of Inhibitors for the Spike-ACE2 Target of SARS-CoV-2. Int J Mol Sci 2023; 24:ijms24108814. [PMID: 37240165 DOI: 10.3390/ijms24108814] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2022] [Revised: 11/23/2022] [Accepted: 11/28/2022] [Indexed: 05/28/2023] Open
Abstract
When an epidemic started in the Chinese city of Wuhan in December 2019, coronavirus was identified as the cause. Infection by the virus occurs through the interaction of viral S protein with the hosts' angiotensin-converting enzyme 2. By leveraging resources such as the DrugBank database and bioinformatics techniques, ligands with potential activity against the SARS-CoV-2 spike protein were designed and identified in this investigation. The FTMap server and the Molegro software were used to determine the active site of the Spike-ACE2 protein's crystal structure. Virtual screening was performed using a pharmacophore model obtained from antiparasitic drugs, obtaining 2000 molecules from molport®. The ADME/Tox profiles were used to identify the most promising compounds with desirable drug characteristics. The binding affinity investigation was then conducted with selected candidates. A molecular docking study showed five structures with better binding affinity than hydroxychloroquine. Ligand_003 showed a binding affinity of -8.645 kcal·mol-1, which was considered an optimal value for the study. The values presented by ligand_033, ligand_013, ligand_044, and ligand_080 meet the profile of novel drugs. To choose compounds with favorable potential for synthesis, synthetic accessibility studies and similarity analyses were carried out. Molecular dynamics and theoretical IC50 values (ranging from 0.459 to 2.371 µM) demonstrate that these candidates are promising for further tests. Chemical descriptors showed that the candidates had strong molecule stability. Theoretical analyses here show that these molecules have potential as SARS-CoV-2 antivirals and therefore warrant further investigation.
Collapse
Affiliation(s)
- Ruan S Bastos
- Graduate Program in Medicinal Chemistry and Molecular Modeling, Federal University of Pará, Belem 66075-110, PA, Brazil
- Laboratory of Modeling and Computational Chemistry, Department of Biological and Health Sciences, Federal University of Amapá, Macapa 68903-419, AP, Brazil
| | - Lúcio R de Lima
- Graduate Program in Medicinal Chemistry and Molecular Modeling, Federal University of Pará, Belem 66075-110, PA, Brazil
- Laboratory of Modeling and Computational Chemistry, Department of Biological and Health Sciences, Federal University of Amapá, Macapa 68903-419, AP, Brazil
| | - Moysés F A Neto
- Laboratory of Molecular Modeling, State University of Feira de Santana, Feira de Santana 44036-900, BA, Brazil
| | - Numan Yousaf
- Department of Biosciences, COMSATS University Islamabad, Park Road, Islamabad 45550, Pakistan
| | - Jorddy N Cruz
- Laboratory of Modeling and Computational Chemistry, Department of Biological and Health Sciences, Federal University of Amapá, Macapa 68903-419, AP, Brazil
| | - Joaquín M Campos
- Department of Pharmaceutical and Organic Chemistry, Faculty of Pharmacy, Campus of Cartuja, University of Granada, 18071 Granada, Spain
- Biosanitary Institute of Granada (ibs.GRANADA), University of Granada, 18071 Granada, Spain
| | - Njogu M Kimani
- Department of Physical Sciences, University of Embu, Embu 6-60100, Kenya
| | - Ryan S Ramos
- Laboratory of Modeling and Computational Chemistry, Department of Biological and Health Sciences, Federal University of Amapá, Macapa 68903-419, AP, Brazil
| | - Cleydson B R Santos
- Graduate Program in Medicinal Chemistry and Molecular Modeling, Federal University of Pará, Belem 66075-110, PA, Brazil
- Laboratory of Modeling and Computational Chemistry, Department of Biological and Health Sciences, Federal University of Amapá, Macapa 68903-419, AP, Brazil
| |
Collapse
|
33
|
Chen W, Cui D, Jerome SV, Michino M, Lenselink EB, Huggins DJ, Beautrait A, Vendome J, Abel R, Friesner RA, Wang L. Enhancing Hit Discovery in Virtual Screening through Absolute Protein-Ligand Binding Free-Energy Calculations. J Chem Inf Model 2023; 63:3171-3185. [PMID: 37167486 DOI: 10.1021/acs.jcim.3c00013] [Citation(s) in RCA: 27] [Impact Index Per Article: 27.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
In the hit identification stage of drug discovery, a diverse chemical space needs to be explored to identify initial hits. Contrary to empirical scoring functions, absolute protein-ligand binding free-energy perturbation (ABFEP) provides a theoretically more rigorous and accurate description of protein-ligand binding thermodynamics and could, in principle, greatly improve the hit rates in virtual screening. In this work, we describe an implementation of an accurate and reliable ABFEP method in FEP+. We validated the ABFEP method on eight congeneric compound series binding to eight protein receptors including both neutral and charged ligands. For ligands with net charges, the alchemical ion approach is adopted to avoid artifacts in electrostatic potential energy calculations. The calculated binding free energies correlate with experimental results with a weighted average of R2 = 0.55 for the entire dataset. We also observe an overall root-mean-square error (RMSE) of 1.1 kcal/mol after shifting the zero-point of the simulation data to match the average experimental values. Through ABFEP calculations using apo versus holo protein structures, we demonstrated that the protein conformational and protonation state changes between the apo and holo proteins are the main physical factors contributing to the protein reorganization free energy manifested by the overestimation of raw ABFEP calculated binding free energies using the holo structures of the proteins. Furthermore, we performed ABFEP calculations in three virtual screening applications for hit enrichment. ABFEP greatly improves the hit rates as compared to docking scores or other methods like metadynamics. The good performance of ABFEP in rank ordering compounds demonstrated in this work confirms it as a useful tool to improve the hit rates in virtual screening, thus facilitating hit discovery.
Collapse
Affiliation(s)
- Wei Chen
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - Di Cui
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - Steven V Jerome
- Schrödinger, Inc., 10201 Wateridge Circle, Suite 220, San Diego, California 92121, United States
| | - Mayako Michino
- Tri-Institutional Therapeutics Discovery Institute, 413 E. 69th Street, New York, New York 10065, United States
| | | | - David J Huggins
- Tri-Institutional Therapeutics Discovery Institute, 413 E. 69th Street, New York, New York 10065, United States
- Department of Physiology and Biophysics, Weill Cornell Medical College of Cornell University, New York, New York 10065, United States
| | - Alexandre Beautrait
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - Jeremie Vendome
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - Robert Abel
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Lingle Wang
- Schrödinger, Inc., 1540 Broadway, 24th Floor, New York, New York 10036, United States
| |
Collapse
|
34
|
Sipala F, Cavallaro G, Forte G, Satriano C, Giuffrida A, Fraix A, Spadaro A, Petralia S, Bonaccorso C, Fortuna CG, Ronsisvalle S. Different In Silico Approaches Using Heterocyclic Derivatives against the Binding between Different Lineages of SARS-CoV-2 and ACE2. Molecules 2023; 28:molecules28093908. [PMID: 37175318 PMCID: PMC10180195 DOI: 10.3390/molecules28093908] [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: 03/22/2023] [Revised: 04/24/2023] [Accepted: 05/04/2023] [Indexed: 05/15/2023] Open
Abstract
Over the last few years, the study of the SARS-CoV-2 spike protein and its mutations has become essential in understanding how it interacts with human host receptors. Since the crystallized structure of the spike protein bound to the angiotensin-converting enzyme 2 (ACE2) receptor was released (PDB code 6M0J), in silico studies have been performed to understand the interactions between these two proteins. Specifically, in this study, heterocyclic compounds with different chemical characteristics were examined to highlight the possibility of interaction with the spike protein and the disruption of the interaction between ACE2 and the spike protein. Our results showed that these compounds interacted with the spike protein and interposed in the interaction zone with ACE2. Although further studies are needed, this work points to these heterocyclic push-pull compounds as possible agents capable of interacting with the spike protein, with the potential for the inhibition of spike protein-ACE2 binding.
Collapse
Affiliation(s)
- Federica Sipala
- Department of Drug and Health Sciences, University of Catania, Viale A. Doria 6, 95125 Catania, Italy
| | - Gianfranco Cavallaro
- Department of Chemical Science, University of Catania, Viale A. Doria 6, 95125 Catania, Italy
| | - Giuseppe Forte
- Department of Drug and Health Sciences, University of Catania, Viale A. Doria 6, 95125 Catania, Italy
| | - Cristina Satriano
- Department of Chemical Science, University of Catania, Viale A. Doria 6, 95125 Catania, Italy
| | - Alessandro Giuffrida
- Department of Chemical Science, University of Catania, Viale A. Doria 6, 95125 Catania, Italy
| | - Aurore Fraix
- Department of Drug and Health Sciences, University of Catania, Viale A. Doria 6, 95125 Catania, Italy
| | - Angelo Spadaro
- Department of Drug and Health Sciences, University of Catania, Viale A. Doria 6, 95125 Catania, Italy
| | - Salvatore Petralia
- Department of Drug and Health Sciences, University of Catania, Viale A. Doria 6, 95125 Catania, Italy
| | - Carmela Bonaccorso
- Department of Chemical Science, University of Catania, Viale A. Doria 6, 95125 Catania, Italy
| | - Cosimo Gianluca Fortuna
- Department of Chemical Science, University of Catania, Viale A. Doria 6, 95125 Catania, Italy
| | - Simone Ronsisvalle
- Department of Drug and Health Sciences, University of Catania, Viale A. Doria 6, 95125 Catania, Italy
| |
Collapse
|
35
|
Sabanés Zariquiey F, Pérez A, Majewski M, Gallicchio E, De Fabritiis G. Validation of the Alchemical Transfer Method for the Estimation of Relative Binding Affinities of Molecular Series. J Chem Inf Model 2023; 63:2438-2444. [PMID: 37042797 PMCID: PMC10577236 DOI: 10.1021/acs.jcim.3c00178] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/13/2023]
Abstract
The accurate prediction of protein-ligand binding affinities is crucial for drug discovery. Alchemical free energy calculations have become a popular tool for this purpose. However, the accuracy and reliability of these methods can vary depending on the methodology. In this study, we evaluate the performance of a relative binding free energy protocol based on the alchemical transfer method (ATM), a novel approach based on a coordinate transformation that swaps the positions of two ligands. The results show that ATM matches the performance of more complex free energy perturbation (FEP) methods in terms of Pearson correlation but with marginally higher mean absolute errors. This study shows that the ATM method is competitive compared to more traditional methods in speed and accuracy and offers the advantage of being applicable with any potential energy function.
Collapse
Affiliation(s)
- Francesc Sabanés Zariquiey
- Computational Science Laboratory, Universitat Pompeu Fabra, Barcelona Biomedical Research Park (PRBB), C Dr. Aiguader 88, 08003 Barcelona, Spain
| | - Adrià Pérez
- Acellera Labs, C Dr Trueta 183, 08005 Barcelona, Spain
| | | | - Emilio Gallicchio
- Department of Chemistry, Brooklyn College of the City University of New York, New York, New York 11210, United States
- PhD Program in Chemistry Graduate Center of the City University of New York, New York, New York 10016, United States
- PhD Program in Biochemistry, Graduate Center of the City University of New York, New York, New York 10016, United States
| | - Gianni De Fabritiis
- Computational Science Laboratory, Universitat Pompeu Fabra, Barcelona Biomedical Research Park (PRBB), C Dr. Aiguader 88, 08003 Barcelona, Spain
- Acellera, Devonshire House 582 Honeypot Lane, Stanmore, Middlesex HA7 1JS, United Kingdom
- Institució Catalana de Recerca i Estudis Avançats (ICREA), Passeig Lluis Companys 23, 08010 Barcelona, Spain
| |
Collapse
|
36
|
Parui S, Robertson JC, Somani S, Tresadern G, Liu C, Dill KA. MELD-Bracket Ranks Binding Affinities of Diverse Sets of Ligands. J Chem Inf Model 2023; 63:2857-2865. [PMID: 37093848 DOI: 10.1021/acs.jcim.3c00243] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2023]
Abstract
Affinity ranking of structurally diverse small-molecule ligands is a challenging problem with important applications in structure-based drug discovery. Absolute binding free energy methods can model diverse ligands, but the high computational cost of the current methods limits application to data sets with few ligands. We recently developed MELD-Bracket, a Molecular Dynamics method for efficient affinity ranking of ligands [ JCTC 2022, 18 (1), 374-379]. It utilizes a Bayesian framework to guide sampling to relevant regions of phase space, and it couples this with a bracket-like competition on a pool of ligands. Here we find that 6-competitor MELD-Bracket can rank dozens of diverse ligands that have low structural similarity and different net charges. We benchmark it on four protein systems─PTB1B, Tyk2, BACE, and JAK3─having varied modes of interactions. We also validated 8-competitor and 12-competitor protocols. The MELD-Bracket protocols presented here may have the appropriate balance of accuracy and computational efficiency to be suitable for ranking diverse ligands from typical drug discovery campaigns.
Collapse
Affiliation(s)
- Sridip Parui
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, United States
| | - James C Robertson
- Janssen Research and Development, Spring House, Pennsylvania 19477, United States
| | - Sandeep Somani
- Janssen Research and Development, Spring House, Pennsylvania 19477, United States
| | - Gary Tresadern
- Janssen Research and Development, Turnhoutseweg 30, Beerse B-2340, Belgium
| | - Cong Liu
- Center for the Development of Therapeutics, Broad Institute of MIT and Harvard, Cambridge, Massachusetts 02142, United States
| | - Ken A Dill
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, United States
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794, United States
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794, United States
| |
Collapse
|
37
|
Fernández-Sánchez SY, Cerón-Carrasco JP, Risco C, Fernández de Castro I. Antiviral Activity of Acetylsalicylic Acid against Bunyamwera Virus in Cell Culture. Viruses 2023; 15:v15040948. [PMID: 37112928 PMCID: PMC10141918 DOI: 10.3390/v15040948] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2023] [Revised: 04/05/2023] [Accepted: 04/06/2023] [Indexed: 04/29/2023] Open
Abstract
The Bunyavirales order is a large group of RNA viruses that includes important pathogens for humans, animals and plants. With high-throughput screening of clinically tested compounds we have looked for potential inhibitors of the endonuclease domain of a bunyavirus RNA polymerase. From a list of fifteen top candidates, five compounds were selected and their antiviral properties studied with Bunyamwera virus (BUNV), a prototypic bunyavirus widely used for studies about the biology of this group of viruses and to test antivirals. Four compounds (silibinin A, myricetin, L-phenylalanine and p-aminohippuric acid) showed no antiviral activity in BUNV-infected Vero cells. On the contrary, acetylsalicylic acid (ASA) efficiently inhibited BUNV infection with a half maximal inhibitory concentration (IC50) of 2.02 mM. In cell culture supernatants, ASA reduced viral titer up to three logarithmic units. A significant dose-dependent reduction of the expression levels of Gc and N viral proteins was also measured. Immunofluorescence and confocal microscopy showed that ASA protects the Golgi complex from the characteristic BUNV-induced fragmentation in Vero cells. Electron microscopy showed that ASA inhibits the assembly of Golgi-associated BUNV spherules that are the replication organelles of bunyaviruses. As a consequence, the assembly of new viral particles is also significantly reduced. Considering its availability and low cost, the potential usability of ASA to treat bunyavirus infections deserves further investigation.
Collapse
Affiliation(s)
| | - José P Cerón-Carrasco
- Centro Universitario de la Defensa, Universidad Politécnica de Cartagena, C/Coronel López Peña s/n, Base Aérea de San Javier, Santiago de la Ribera, 30720 Murcia, Spain
| | - Cristina Risco
- Cell Structure Laboratory, Centro Nacional de Biotecnología, CSIC, Campus de Cantoblanco, 28049 Madrid, Spain
| | - Isabel Fernández de Castro
- Cell Structure Laboratory, Centro Nacional de Biotecnología, CSIC, Campus de Cantoblanco, 28049 Madrid, Spain
| |
Collapse
|
38
|
Gul I, Hassan A, Haq E, Ahmad SM, Shah RA, Ganai NA, Chikan NA, Abdul-Careem MF, Shabir N. An Investigation of the Antiviral Potential of Phytocompounds against Avian Infectious Bronchitis Virus through Template-Based Molecular Docking and Molecular Dynamics Simulation Analysis. Viruses 2023; 15:v15040847. [PMID: 37112828 PMCID: PMC10144825 DOI: 10.3390/v15040847] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2023] [Revised: 03/17/2023] [Accepted: 03/24/2023] [Indexed: 03/29/2023] Open
Abstract
Vaccination is widely used to control Infectious Bronchitis in poultry; however, the limited cross-protection and safety issues associated with these vaccines can lead to vaccination failures. Keeping these limitations in mind, the current study explored the antiviral potential of phytocompounds against the Infectious Bronchitis virus using in silico approaches. A total of 1300 phytocompounds derived from fourteen botanicals were screened for their potential ability to inhibit the main protease, papain-like protease or RNA-dependent RNA–polymerase of the virus. The study identified Methyl Rosmarinate, Cianidanol, Royleanone, and 6,7-Dehydroroyleanone as dual-target inhibitors against any two of the key proteins. At the same time, 7-alpha-Acetoxyroyleanone from Rosmarinus officinalis was found to be a multi-target protein inhibitor against all three proteins. The potential multi-target inhibitor was subjected to molecular dynamics simulations to assess the stability of the protein–ligand complexes along with the corresponding reference ligands. The findings specified stable interactions of 7-alpha-Acetoxyroyleanone with the protein targets. The results based on the in silico study indicate that the phytocompounds can potentially inhibit the essential proteins of the Infectious Bronchitis virus; however, in vitro and in vivo studies are required for validation. Nevertheless, this study is a significant step in exploring the use of botanicals in feed to control Infectious Bronchitis infections in poultry.
Collapse
Affiliation(s)
- Irfan Gul
- Division of Animal Biotechnology, Faculty of Veterinary Sciences and Animal Husbandry, Shuhama, Sher-e-Kashmir University of Agricultural Sciences and Technology of Kashmir, Srinagar 190006, India; (I.G.)
- Department of Biotechnology, University of Kashmir, Srinagar 190006, India
| | - Amreena Hassan
- Division of Animal Biotechnology, Faculty of Veterinary Sciences and Animal Husbandry, Shuhama, Sher-e-Kashmir University of Agricultural Sciences and Technology of Kashmir, Srinagar 190006, India; (I.G.)
- Department of Biotechnology, University of Kashmir, Srinagar 190006, India
| | - Ehtishamul Haq
- Department of Biotechnology, University of Kashmir, Srinagar 190006, India
| | - Syed Mudasir Ahmad
- Division of Animal Biotechnology, Faculty of Veterinary Sciences and Animal Husbandry, Shuhama, Sher-e-Kashmir University of Agricultural Sciences and Technology of Kashmir, Srinagar 190006, India; (I.G.)
| | - Riaz Ahmad Shah
- Division of Animal Biotechnology, Faculty of Veterinary Sciences and Animal Husbandry, Shuhama, Sher-e-Kashmir University of Agricultural Sciences and Technology of Kashmir, Srinagar 190006, India; (I.G.)
| | - Nazir Ahmad Ganai
- Division of Animal Biotechnology, Faculty of Veterinary Sciences and Animal Husbandry, Shuhama, Sher-e-Kashmir University of Agricultural Sciences and Technology of Kashmir, Srinagar 190006, India; (I.G.)
| | - Naveed Anjum Chikan
- Division of Computational Biology, Daskdan Innovations, Pvt. Ltd., Kashmir 190006, India
| | - Mohamed Faizal Abdul-Careem
- Faculty of Veterinary Medicine, University of Calgary, 3330 Hospital Drive NW, Calgary, AB T2N 4N1, Canada
- Correspondence: (M.F.A.-C.); (N.S.)
| | - Nadeem Shabir
- Division of Animal Biotechnology, Faculty of Veterinary Sciences and Animal Husbandry, Shuhama, Sher-e-Kashmir University of Agricultural Sciences and Technology of Kashmir, Srinagar 190006, India; (I.G.)
- Correspondence: (M.F.A.-C.); (N.S.)
| |
Collapse
|
39
|
Nakamura S, Akaki T, Nishiwaki K, Nakatani M, Kawase Y, Takahashi Y, Nakanishi I. System truncation accelerates binding affinity calculations with the fragment molecular orbital method: A benchmark study. J Comput Chem 2023; 44:824-831. [PMID: 36444861 DOI: 10.1002/jcc.27044] [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: 08/07/2022] [Revised: 11/01/2022] [Accepted: 11/03/2022] [Indexed: 11/30/2022]
Abstract
The fragment molecular orbital (FMO) method is a fast quantum-mechanical method that divides systems into pieces of fragments and performs ab initio calculations. The system truncation enables further speed improvement. In this article, we systematically study the effects of system truncations on binding affinity calculations obtained with FMO in combination with either the polarizable continuum model (FMO/PCM) or in combination with the Møller-Plesset method (FMO-MP2). We have used five protein complexes with ligands of several charged states. The calculated binding energies of the size variants of the truncated system, including only a restricted number of atoms around the ligand, are compared to the energy obtained from a full system. The result shows that the systems could be truncated to a radius of 8 Å from neutral ligands within an error of 0.7 kcal/mol, and 12 Å from charged ligands within an error of 1.1 kcal/mol for calculating the binding energy in solution.
Collapse
Affiliation(s)
- Shinya Nakamura
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Tatsuo Akaki
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan.,Chemical Research Laboratories, Central Pharmaceutical Research Institute, Japan Tobacco Inc., Osaka, Japan
| | - Keiji Nishiwaki
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Midori Nakatani
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Yuji Kawase
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Yuki Takahashi
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| | - Isao Nakanishi
- Computational Drug Design and Discovery, Department of Pharmaceutical Sciences, Kindai University, Osaka, Japan
| |
Collapse
|
40
|
Barros EP, Ries B, Champion C, Rieder SR, Riniker S. Accounting for Solvation Correlation Effects on the Thermodynamics of Water Networks in Protein Cavities. J Chem Inf Model 2023; 63:1794-1805. [PMID: 36917685 PMCID: PMC10052353 DOI: 10.1021/acs.jcim.2c01610] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/15/2023]
Abstract
Macromolecular recognition and ligand binding are at the core of biological function and drug discovery efforts. Water molecules play a significant role in mediating the protein-ligand interaction, acting as more than just the surrounding medium by affecting the thermodynamics and thus the outcome of the binding process. As individual water contributions are impossible to measure experimentally, a range of computational methods have emerged to identify hydration sites in protein pockets and characterize their energetic contributions for drug discovery applications. Even though several methods model solvation effects explicitly, they focus on determining the stability of specific water sites independently and neglect solvation correlation effects upon replacement of clusters of water molecules, which typically happens in hit-to-lead optimization. In this work, we rigorously determine the conjoint effects of replacing all combinations of water molecules in protein binding pockets through the use of the RE-EDS multistate free-energy method, which combines Hamiltonian replica exchange (RE) and enveloping distribution sampling (EDS). Applications on the small bovine pancreatic trypsin inhibitor and four proteins of the bromodomain family illustrate the extent of solvation correlation effects on water thermodynamics, with the favorability of replacement of the water sites by pharmacophore probes highly dependent on the composition of the water network and the pocket environment. Given the ubiquity of water networks in biologically relevant protein targets, we believe our approach can be helpful for computer-aided drug discovery by providing a pocket-specific and a priori systematic consideration of solvation effects on ligand binding and selectivity.
Collapse
Affiliation(s)
- Emilia P Barros
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Benjamin Ries
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Candide Champion
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Salomé R Rieder
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Sereina Riniker
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
41
|
Borkakoti N, Thornton JM. AlphaFold2 protein structure prediction: Implications for drug discovery. Curr Opin Struct Biol 2023; 78:102526. [PMID: 36621153 PMCID: PMC7614146 DOI: 10.1016/j.sbi.2022.102526] [Citation(s) in RCA: 14] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2022] [Revised: 12/01/2022] [Accepted: 12/03/2022] [Indexed: 01/09/2023]
Abstract
The drug discovery process involves designing compounds to selectively interact with their targets. The majority of therapeutic targets for low molecular weight (small molecule) drugs are proteins. The outstanding accuracy with which recent artificial intelligence methods compile the three-dimensional structure of proteins has made protein targets more accessible to the drug design process. Here, we present our perspective of the significance of accurate protein structure prediction on various stages of the small molecule drug discovery life cycle focusing on current capabilities and assessing how further evolution of such predictive procedures can have a more decisive impact in the discovery of new medicines.
Collapse
Affiliation(s)
- Neera Borkakoti
- European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SD, UK
| | - Janet M Thornton
- European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SD, UK.
| |
Collapse
|
42
|
Xu H. The slow but steady rise of binding free energy calculations in drug discovery. J Comput Aided Mol Des 2023; 37:67-74. [PMID: 36469232 DOI: 10.1007/s10822-022-00494-x] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2022] [Accepted: 11/28/2022] [Indexed: 12/12/2022]
Abstract
Binding free energy calculations are increasingly used in drug discovery research to predict protein-ligand binding affinities and to prioritize candidate drug molecules accordingly. It has taken decades of collective effort to transform this academic concept into a technology adopted by the pharmaceutical and biotech industry. Having personally witnessed and taken part in this transformation, here I recount the (incomplete) list of problems that had to be solved to make this computational tool practical and suggest areas of future development.
Collapse
Affiliation(s)
- Huafeng Xu
- Roivant Discovery, 151 West 42nd Street, New York, NY, 10036, USA.
| |
Collapse
|
43
|
Gusev F, Gutkin E, Kurnikova MG, Isayev O. Active Learning Guided Drug Design Lead Optimization Based on Relative Binding Free Energy Modeling. J Chem Inf Model 2023; 63:583-594. [PMID: 36599125 DOI: 10.1021/acs.jcim.2c01052] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
In silico identification of potent protein inhibitors commonly requires prediction of a ligand binding free energy (BFE). Thermodynamics integration (TI) based on molecular dynamics (MD) simulations is a BFE calculation method capable of acquiring accurate BFE, but it is computationally expensive and time-consuming. In this work, we have developed an efficient automated workflow for identifying compounds with the lowest BFE among thousands of congeneric ligands, which requires only hundreds of TI calculations. Automated machine learning (AutoML) orchestrated by active learning (AL) in an AL-AutoML workflow allows unbiased and efficient search for a small set of best-performing molecules. We have applied this workflow to select inhibitors of the SARS-CoV-2 papain-like protease and were able to find 133 compounds with improved binding affinity, including 16 compounds with better than 100-fold binding affinity improvement. We obtained a hit rate that outperforms that expected of traditional expert medicinal chemist-guided campaigns. Thus, we demonstrate that the combination of AL and AutoML with free energy simulations provides at least 20× speedup relative to the naïve brute force approaches.
Collapse
Affiliation(s)
- Filipp Gusev
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania15213, United States
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, Pennsylvania15213, United States
| | - Evgeny Gutkin
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania15213, United States
| | - Maria G Kurnikova
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania15213, United States
| | - Olexandr Isayev
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania15213, United States
- Computational Biology Department, Carnegie Mellon University, Pittsburgh, Pennsylvania15213, United States
| |
Collapse
|
44
|
Breznik M, Ge Y, Bluck JP, Briem H, Hahn DF, Christ CD, Mortier J, Mobley DL, Meier K. Prioritizing Small Sets of Molecules for Synthesis through in-silico Tools: A Comparison of Common Ranking Methods. ChemMedChem 2023; 18:e202200425. [PMID: 36240514 PMCID: PMC9868080 DOI: 10.1002/cmdc.202200425] [Citation(s) in RCA: 20] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2022] [Revised: 10/10/2022] [Indexed: 01/26/2023]
Abstract
Prioritizing molecules for synthesis is a key role of computational methods within medicinal chemistry. Multiple tools exist for ranking molecules, from the cheap and popular molecular docking methods to more computationally expensive molecular-dynamics (MD)-based methods. It is often questioned whether the accuracy of the more rigorous methods justifies the higher computational cost and associated calculation time. Here, we compared the performance on ranking the binding of small molecules for seven scoring functions from five docking programs, one end-point method (MM/GBSA), and two MD-based free energy methods (PMX, FEP+). We investigated 16 pharmaceutically relevant targets with a total of 423 known binders. The performance of docking methods for ligand ranking was strongly system dependent. We observed that MD-based methods predominantly outperformed docking algorithms and MM/GBSA calculations. Based on our results, we recommend the application of MD-based free energy methods for prioritization of molecules for synthesis in lead optimization, whenever feasible.
Collapse
Affiliation(s)
- Marko Breznik
- Computational Molecular Design, Pharmaceuticals, R&D, Bayer AG, 13342 Berlin, Germany
| | - Yunhui Ge
- Department of Pharmaceutical Sciences, University of California, Irvine, CA 92697, USA
| | - Joseph P. Bluck
- Computational Molecular Design, Pharmaceuticals, R&D, Bayer AG, 13342 Berlin, Germany
| | - Hans Briem
- Computational Molecular Design, Pharmaceuticals, R&D, Bayer AG, 13342 Berlin, Germany
| | - David F. Hahn
- Computational Chemistry, Janssen Research & Development, Turnhoutseweg 30, Beerse B-2340, Belgium
| | - Clara D. Christ
- Molecular Design, Pharmaceuticals, R&D, Bayer AG, 13342 Berlin, Germany
| | - Jérémie Mortier
- Computational Molecular Design, Pharmaceuticals, R&D, Bayer AG, 13342 Berlin, Germany
| | - David L. Mobley
- Department of Pharmaceutical Sciences, University of California, Irvine, CA 92697, USA,Department of Chemistry, University of California, Irvine, CA 92697, USA
| | - Katharina Meier
- Computational Life Science Technology Functions, Crop Science, R&D, Bayer AG, 40789 Monheim, Germany
| |
Collapse
|
45
|
Schuurs ZP, McDonald JP, Croft LV, Richard DJ, Woodgate R, Gandhi NS. Integration of molecular modelling and in vitro studies to inhibit LexA proteolysis. Front Cell Infect Microbiol 2023; 13:1051602. [PMID: 36936756 PMCID: PMC10020695 DOI: 10.3389/fcimb.2023.1051602] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2022] [Accepted: 02/14/2023] [Indexed: 03/06/2023] Open
Abstract
Introduction As antibiotic resistance has become more prevalent, the social and economic impacts are increasingly pressing. Indeed, bacteria have developed the SOS response which facilitates the evolution of resistance under genotoxic stress. The transcriptional repressor, LexA, plays a key role in this response. Mutation of LexA to a non-cleavable form that prevents the induction of the SOS response sensitizes bacteria to antibiotics. Achieving the same inhibition of proteolysis with small molecules also increases antibiotic susceptibility and reduces drug resistance acquisition. The availability of multiple LexA crystal structures, and the unique Ser-119 and Lys-156 catalytic dyad in the protein enables the rational design of inhibitors. Methods We pursued a binary approach to inhibit proteolysis; we first investigated β-turn mimetics, and in the second approach we tested covalent warheads targeting the Ser-119 residue. We found that the cleavage site region (CSR) of the LexA protein is a classical Type II β-turn, and that published 1,2,3-triazole compounds mimic the β-turn. Generic covalent molecule libraries and a β-turn mimetic library were docked to the LexA C-terminal domain using molecular modelling methods in FlexX and CovDock respectively. The 133 highest-scoring molecules were screened for their ability to inhibit LexA cleavage under alkaline conditions. The top molecules were then tested using a RecA-mediated cleavage assay. Results The β-turn library screen did not produce any hit compounds that inhibited RecA-mediated cleavage. The covalent screen discovered an electrophilic serine warhead that can inhibit LexA proteolysis, reacting with Ser-119 via a nitrile moiety. Discussion This research presents a starting point for hit-to-lead optimisation, which could lead to inhibition of the SOS response and prevent the acquisition of antibiotic resistance.
Collapse
Affiliation(s)
- Zachariah P. Schuurs
- Cancer and Ageing Research Program, Centre for Genomics and Personalised Health, Queensland University of Technology (QUT), Translational Research Institute (TRI), Brisbane, QLD, Australia
- School of Chemistry and Physics, Queensland University of Technology (QUT), Brisbane, QLD, Australia
| | - John P. McDonald
- Laboratory of Genomic Integrity, National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD, United States
| | - Laura V. Croft
- Cancer and Ageing Research Program, Centre for Genomics and Personalised Health, Queensland University of Technology (QUT), Translational Research Institute (TRI), Brisbane, QLD, Australia
| | - Derek J. Richard
- Cancer and Ageing Research Program, Centre for Genomics and Personalised Health, Queensland University of Technology (QUT), Translational Research Institute (TRI), Brisbane, QLD, Australia
| | - Roger Woodgate
- Laboratory of Genomic Integrity, National Institute of Child Health and Human Development, National Institutes of Health, Bethesda, MD, United States
- *Correspondence: Neha S. Gandhi, ; Roger Woodgate,
| | - Neha S. Gandhi
- Cancer and Ageing Research Program, Centre for Genomics and Personalised Health, Queensland University of Technology (QUT), Translational Research Institute (TRI), Brisbane, QLD, Australia
- School of Chemistry and Physics, Queensland University of Technology (QUT), Brisbane, QLD, Australia
- *Correspondence: Neha S. Gandhi, ; Roger Woodgate,
| |
Collapse
|
46
|
Roja R, Kalakotla S, Ravula AR, Boyina HK, Navanita SK, Vallika PBS, Gangarapu K, Devarakonda KP, Bakshi V. Insilico Screening of Pentacyclic Triterpenoids against Vascular Dementia Target's. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2023; 1423:237-243. [PMID: 37525050 DOI: 10.1007/978-3-031-31978-5_23] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/02/2023]
Abstract
Vascular dementia (VaD) accounts to 30% of cases and is predicted as second most common form of dementia after Alzheimer's disease by WHO. Earlier studies reported that plant-derived pentacyclic triterpenoids possess a wide range of pharmacological activities but these compounds are not extensively studied for their neuroprotective potential against VaD. This in silico approach was designed to screen 20 pentacyclic triterpenoid plant compounds against known targets of VaD using Flare software. S-Adenyl homocysteine hydrolase, Acetylcholinesterase, and Butyrylcholinesterase were selected as important VaD targets, and various parameters like intermolecular interaction energies, binding energy, and dock scores were analyzed and compared between selected ligands. Our results showed that Ursolic acid has lowest binding energy when docked with most of the target proteins, and among all 20 pentacyclic triterpenoids studied, only three ligands Betulinic acid, Ambolic acid, and Madecassic acid, showed better binding energy scores, and they can be shortlisted as lead compounds to further study their therapeutic potential against VaD using in vitro and in vivo animal models.
Collapse
Affiliation(s)
- Rathna Roja
- School of Pharmacy, Department of Pharmacology, Anurag University, Hyderabad, Telangana, India
| | - Shanker Kalakotla
- Department of Pharmacognosy & Phyto-Pharmacy, JSS College of Pharmacy, JSS Academy of Higher Education & Research, Nilgiris, Tamil Nadu, India
| | - Arun Reddy Ravula
- Department of Biomedical Engineering, Center for Injury Biomechanics, Materials and Medicine, New Jersey Institute of Technology, Newark, NJ, USA
| | - Hemanth Kumar Boyina
- School of Pharmacy, Department of Pharmacology, Anurag University, Hyderabad, Telangana, India
| | - S K Navanita
- Department of Pharmacognosy & Phyto-Pharmacy, JSS College of Pharmacy, JSS Academy of Higher Education & Research, Nilgiris, Tamil Nadu, India
| | | | - Kiran Gangarapu
- School of Pharmacy, Department of Pharmaceutical Analysis, Anurag University, Hyderabad, Telangana, India
| | | | - Vasudha Bakshi
- School of Pharmacy, Department of Pharmacology, Anurag University, Hyderabad, Telangana, India
| |
Collapse
|
47
|
Andole S, Sd H, Sudhula S, Vislavath L, Boyina HK, Gangarapu K, Bakshi V, Devarakonda KP. 3D QSAR based Virtual Screening of Flavonoids as Acetylcholinesterase Inhibitors. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2023; 1424:233-240. [PMID: 37486499 DOI: 10.1007/978-3-031-31982-2_26] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/25/2023]
Abstract
In an attempt to develop therapeutic agents to treat Alzheimer's disease, a series of flavonoid analogues were collected, which already had established acetylcholinesterase (AChE) enzyme inhibition activity. For each molecule we also collected biological activity data (Ki). Then, 3D-QSAR (quantitative structure-activity relationship model) was developed which showed acceptable predictive and descriptive capability as represented by standard statistical parameters r2 and q2. This SAR data can explain the key descriptors which can be related to AChE inhibitory activity. Using the QSAR model, pharmacophores were developed based on which, virtual screening was done and a dataset was obtained which loaded as a prediction set to fit the developed QSAR model. Top 10 compounds fitting the QSAR model were subjected to molecular docking. CHEMBL1718051 was found to be the lead compound. This study is offering an example of a computationally-driven tool for prioritisation and discovery of probable AChE inhibitors. Further, in vivo and in vitro testing will show its therapeutic potential.
Collapse
Affiliation(s)
- Sowmya Andole
- School of Pharmacy, Department of Pharmacology, Anurag University, Hyderabad, Telangana, India
| | - Husna Sd
- School of Pharmacy, Department of Pharmacology, Anurag University, Hyderabad, Telangana, India
| | - Srija Sudhula
- School of Pharmacy, Department of Pharmacology, Anurag University, Hyderabad, Telangana, India
| | - Lavanya Vislavath
- School of Pharmacy, Department of Pharmacology, Anurag University, Hyderabad, Telangana, India
| | - Hemanth Kumar Boyina
- School of Pharmacy, Department of Pharmacology, Anurag University, Hyderabad, Telangana, India
| | - Kiran Gangarapu
- School of Pharmacy, Department of Pharmaceutical Analysis, Anurag University, Hyderabad, Telangana, India
| | - Vasudha Bakshi
- School of Pharmacy, Department of Pharmacology, Anurag University, Hyderabad, Telangana, India
| | | |
Collapse
|
48
|
Chang Y, Hawkins BA, Du JJ, Groundwater PW, Hibbs DE, Lai F. A Guide to In Silico Drug Design. Pharmaceutics 2022; 15:pharmaceutics15010049. [PMID: 36678678 PMCID: PMC9867171 DOI: 10.3390/pharmaceutics15010049] [Citation(s) in RCA: 27] [Impact Index Per Article: 13.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2022] [Revised: 12/16/2022] [Accepted: 12/17/2022] [Indexed: 12/28/2022] Open
Abstract
The drug discovery process is a rocky path that is full of challenges, with the result that very few candidates progress from hit compound to a commercially available product, often due to factors, such as poor binding affinity, off-target effects, or physicochemical properties, such as solubility or stability. This process is further complicated by high research and development costs and time requirements. It is thus important to optimise every step of the process in order to maximise the chances of success. As a result of the recent advancements in computer power and technology, computer-aided drug design (CADD) has become an integral part of modern drug discovery to guide and accelerate the process. In this review, we present an overview of the important CADD methods and applications, such as in silico structure prediction, refinement, modelling and target validation, that are commonly used in this area.
Collapse
Affiliation(s)
- Yiqun Chang
- Sydney Pharmacy School, Faculty of Medicine and Health, The University of Sydney, Camperdown, NSW 2006, Australia
| | - Bryson A. Hawkins
- Sydney Pharmacy School, Faculty of Medicine and Health, The University of Sydney, Camperdown, NSW 2006, Australia
| | - Jonathan J. Du
- Department of Biochemistry, Emory University School of Medicine, Atlanta, GA 30322, USA
| | - Paul W. Groundwater
- Sydney Pharmacy School, Faculty of Medicine and Health, The University of Sydney, Camperdown, NSW 2006, Australia
| | - David E. Hibbs
- Sydney Pharmacy School, Faculty of Medicine and Health, The University of Sydney, Camperdown, NSW 2006, Australia
| | - Felcia Lai
- Sydney Pharmacy School, Faculty of Medicine and Health, The University of Sydney, Camperdown, NSW 2006, Australia
- Correspondence:
| |
Collapse
|
49
|
Gizzio J, Thakur A, Haldane A, Levy RM. Evolutionary divergence in the conformational landscapes of tyrosine vs serine/threonine kinases. eLife 2022; 11:83368. [PMID: 36562610 PMCID: PMC9822262 DOI: 10.7554/elife.83368] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/09/2022] [Accepted: 12/22/2022] [Indexed: 12/24/2022] Open
Abstract
Inactive conformations of protein kinase catalytic domains where the DFG motif has a "DFG-out" orientation and the activation loop is folded present a druggable binding pocket that is targeted by FDA-approved 'type-II inhibitors' in the treatment of cancers. Tyrosine kinases (TKs) typically show strong binding affinity with a wide spectrum of type-II inhibitors while serine/threonine kinases (STKs) usually bind more weakly which we suggest here is due to differences in the folded to extended conformational equilibrium of the activation loop between TKs vs. STKs. To investigate this, we use sequence covariation analysis with a Potts Hamiltonian statistical energy model to guide absolute binding free-energy molecular dynamics simulations of 74 protein-ligand complexes. Using the calculated binding free energies together with experimental values, we estimated free-energy costs for the large-scale (~17-20 Å) conformational change of the activation loop by an indirect approach, circumventing the very challenging problem of simulating the conformational change directly. We also used the Potts statistical potential to thread large sequence ensembles over active and inactive kinase states. The structure-based and sequence-based analyses are consistent; together they suggest TKs evolved to have free-energy penalties for the classical 'folded activation loop' DFG-out conformation relative to the active conformation, that is, on average, 4-6 kcal/mol smaller than the corresponding values for STKs. Potts statistical energy analysis suggests a molecular basis for this observation, wherein the activation loops of TKs are more weakly 'anchored' against the catalytic loop motif in the active conformation and form more stable substrate-mimicking interactions in the inactive conformation. These results provide insights into the molecular basis for the divergent functional properties of TKs and STKs, and have pharmacological implications for the target selectivity of type-II inhibitors.
Collapse
Affiliation(s)
- Joan Gizzio
- Center for Biophysics and Computational Biology, Temple University, Philadelphia, United States.,Department of Chemistry, Temple University, Philadelphia, United States
| | - Abhishek Thakur
- Center for Biophysics and Computational Biology, Temple University, Philadelphia, United States.,Department of Chemistry, Temple University, Philadelphia, United States
| | - Allan Haldane
- Center for Biophysics and Computational Biology, Temple University, Philadelphia, United States.,Department of Physics, Temple University, Philadelphia, United States
| | - Ronald M Levy
- Center for Biophysics and Computational Biology, Temple University, Philadelphia, United States.,Department of Chemistry, Temple University, Philadelphia, United States
| |
Collapse
|
50
|
Dolezal R. Accuracy and precision of binding free energy prediction for a tacrine related lead inhibitor of acetylcholinesterase with an arsenal of supercomputerized molecular modelling methods: a comparative study. J Biomol Struct Dyn 2022; 40:11291-11319. [PMID: 34323654 DOI: 10.1080/07391102.2021.1957716] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Abstract
Nowadays, advanced computational chemistry methods offer various strategies for revealing prospective hit structures in drug development essentially through accurate binding free energy predictions. After the era of molecular docking and quantitative structure-activity relationships, much interest has been lately oriented to perturbed molecular dynamic approaches like replica exchange with solute tempering and free energy perturbation (REST/FEP) and the potential of the mean force with adaptive biasing and accelerated weight histograms (PMF/AWH). Both of these receptor-based techniques can exploit exascale CPU&GPU supercomputers to achieve high throughput performance. In this fundamental study, we have compared the predictive power of a panel of supercomputerized molecular modelling methods to distinguish the major binding modes and the corresponding binding free energies of a promising tacrine related potential antialzheimerics in human acetylcholinesterase. The binding free energies were estimated using flexible molecular docking, molecular mechanics/generalized Born surface area/Poisson-Boltzmann surface area (MM/GBSA/PBSA), transmutation REST/FEP with 12 x 5 ns/λ windows, annihilation FEP with 20 x 5 ns/λ steps, PMF with weight histogram analysis method (WHAM) and 40 x 5 ns samples, and PMF/AWH with 10 x 100 ns replicas. Confrontation of the classical approaches such as canonical molecular dynamics and molecular docking with alchemical calculations and steered molecular dynamics enabled us to show how large errors in ΔG predictions can be expected if these in silico methods are employed in the elucidation of a common case of enzyme inhibition.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Rafael Dolezal
- Department of Chemistry, Faculty of Science, University of Hradec Kralove, Hradec Kralove, Czech Republic.,Biomedical Research Center, University Hospital Hradec Kralove, Hradec Kralove, Czech Republic
| |
Collapse
|