1
|
Huang X, Chorianopoulou A, Kalkounou P, Georgiou M, Pousias A, Davies A, Pearce A, Harris M, Lambrinidis G, Marakos P, Pouli N, Kolocouris A, Lougiakis N, Ladds G. Hit-to-Lead Optimization of Heterocyclic Carbonyloxycarboximidamides as Selective Antagonists at Human Adenosine A3 Receptor. J Med Chem 2024; 67:13117-13146. [PMID: 39073853 PMCID: PMC11320584 DOI: 10.1021/acs.jmedchem.4c01092] [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: 05/08/2024] [Revised: 06/18/2024] [Accepted: 07/01/2024] [Indexed: 07/30/2024]
Abstract
Antagonism of the human adenosine A3 receptor (hA3R) has potential therapeutic application. Alchemical relative binding free energy calculations of K18 and K32 suggested that the combination of a 3-(2,6-dichlorophenyl)-isoxazolyl group with 2-pyridinyl at the ends of a carbonyloxycarboximidamide group should improve hA3R affinity. Of the 25 new analogues synthesized, 37 and 74 showed improved hA3R affinity compared to K18 (and K32). This was further improved through the addition of a bromine group to the 2-pyridinyl at the 5-position, generating compound 39. Alchemical relative binding free energy calculations, mutagenesis studies and MD simulations supported the compounds' binding pattern while suggesting that the bromine of 39 inserts deep into the hA3R orthosteric pocket, so highlighting the importance of rigidification of the carbonyloxycarboximidamide moiety. MD simulations highlighted the importance of rigidification of the carbonyloxycarboximidamide, while suggesting that the bromine of 39 inserts deep into the hA3R orthosteric pocket, which was supported through mutagenesis studies 39 also selectively antagonized endogenously expressed hA3R in nonsmall cell lung carcinoma cells, while pharmacokinetic studies indicated low toxicity enabling in vivo evaluation. We therefore suggest that 39 has potential for further development as a high-affinity hA3R antagonist.
Collapse
Affiliation(s)
- Xianglin Huang
- Department
of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, U.K.
| | - Anna Chorianopoulou
- Laboratory
of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department
of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou 15771, Athens, Greece
| | - Panagoula Kalkounou
- Laboratory
of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department
of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou 15771, Athens, Greece
| | - Maria Georgiou
- Laboratory
of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department
of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou 15771, Athens, Greece
| | - Athanasios Pousias
- Laboratory
of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department
of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou 15771, Athens, Greece
| | - Amy Davies
- Department
of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, U.K.
| | - Abigail Pearce
- Department
of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, U.K.
| | - Matthew Harris
- Department
of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, U.K.
| | - George Lambrinidis
- Laboratory
of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department
of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou 15771, Athens, Greece
| | - Panagiotis Marakos
- Laboratory
of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department
of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou 15771, Athens, Greece
| | - Nicole Pouli
- Laboratory
of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department
of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou 15771, Athens, Greece
| | - Antonios Kolocouris
- Laboratory
of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department
of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou 15771, Athens, Greece
| | - Nikolaos Lougiakis
- Laboratory
of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department
of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou 15771, Athens, Greece
| | - Graham Ladds
- Department
of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, U.K.
| |
Collapse
|
2
|
Jiang T, Zhang Y, Yu S, Hu B. Discovering potential WRN inhibitors from natural product database through computational methods. J Mol Graph Model 2024; 129:108758. [PMID: 38507856 DOI: 10.1016/j.jmgm.2024.108758] [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/15/2023] [Revised: 02/29/2024] [Accepted: 03/10/2024] [Indexed: 03/22/2024]
Abstract
Microsatellite instability (MSI) is a relatively common feature associated with multiple cancers, and Werner syndrome (WRN) ATP-dependent helicase has been recognized as a novel target for treating MSI cancers, such as colorectal cancer. A small-molecule inhibitor targeting WRN would be a promising strategy for treating colorectal cancer with high MSI expression. In this study, we employed a computer-assisted drug discovery strategy to screen over 30,000 natural product molecules. By using a combination of docking, ligand efficiency, Molecular Mechanics/Generalized Born Surface Area (MM/GBSA), and thermodynamic integration (TI) calculations, we identified MOL008980, MOL010740, MOL011832, T4743, TN1166, and TNP-002173 as potential WRN inhibitors. Subsequent molecular dynamics simulation revealed that these screened natural products possessed better binding dynamic characteristics than ATP substrate and were capable of inhibiting the dynamic process of WRN, making them potential strong ATP competitive inhibitors. In conclusion, our computational approach revealed potential WRN inhibitors from a natural product database, providing a theoretical basis for future research.
Collapse
Affiliation(s)
| | | | - Shuihong Yu
- College of Pharmacy, Anqing Medical College, Anqing, China
| | - Bingde Hu
- Navy Anqing Hospital, Anqing, China.
| |
Collapse
|
3
|
Stampelou M, Ladds G, Kolocouris A. Computational Workflow for Refining AlphaFold Models in Drug Design Using Kinetic and Thermodynamic Binding Calculations: A Case Study for the Unresolved Inactive Human Adenosine A 3 Receptor. J Phys Chem B 2024; 128:914-936. [PMID: 38236582 DOI: 10.1021/acs.jpcb.3c05986] [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: 01/19/2024]
Abstract
A structure-based drug design pipeline that considers both thermodynamic and kinetic binding data of ligands against a receptor will enable the computational design of improved drug molecules. For unresolved GPCR-ligand complexes, a workflow that can apply both thermodynamic and kinetic binding data in combination with alpha-fold (AF)-derived or other homology models and experimentally resolved binding modes of relevant ligands in GPCR-homologs needs to be tested. Here, as test case, we studied a congeneric set of ligands that bind to a structurally unresolved G protein-coupled receptor (GPCR), the inactive human adenosine A3 receptor (hA3R). We tested three available homology models from which two have been generated from experimental structures of hA1R or hA2AR and one model was a multistate alphafold 2 (AF2)-derived model. We applied alchemical calculations with thermodynamic integration coupled with molecular dynamics (TI/MD) simulations to calculate the experimental relative binding free energies and residence time (τ)-random accelerated MD (τ-RAMD) simulations to calculate the relative residence times (RTs) for antagonists. While the TI/MD calculations produced, for the three homology models, good Pearson correlation coefficients, correspondingly, r = 0.74, 0.62, and 0.67 and mean unsigned error (mue) values of 0.94, 1.31, and 0.81 kcal mol-1, the τ-RAMD method showed r = 0.92 and 0.52 for the first two models but failed to produce accurate results for the multistate AF2-derived model. With subsequent optimization of the AF2-derived model by reorientation of the side chain of R1735.34 located in the extracellular loop 2 (EL2) that blocked ligand's unbinding, the computational model showed r = 0.84 for kinetic data and improved performance for thermodynamic data (r = 0.81, mue = 0.56 kcal mol-1). Overall, after refining the multistate AF2 model with physics-based tools, we were able to show a strong correlation between predicted and experimental ligand relative residence times and affinities, achieving a level of accuracy comparable to an experimental structure. The computational workflow used can be applied to other receptors, helping to rank candidate drugs in a congeneric series and enabling the prioritization of leads with stronger binding affinities and longer residence times.
Collapse
Affiliation(s)
- Margarita Stampelou
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771 Athens, Greece
| | - Graham Ladds
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, U.K
| | - Antonios Kolocouris
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771 Athens, Greece
| |
Collapse
|
4
|
Xia S, Chen E, Zhang Y. Integrated Molecular Modeling and Machine Learning for Drug Design. J Chem Theory Comput 2023; 19:7478-7495. [PMID: 37883810 PMCID: PMC10653122 DOI: 10.1021/acs.jctc.3c00814] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2023] [Revised: 10/10/2023] [Accepted: 10/11/2023] [Indexed: 10/28/2023]
Abstract
Modern therapeutic development often involves several stages that are interconnected, and multiple iterations are usually required to bring a new drug to the market. Computational approaches have increasingly become an indispensable part of helping reduce the time and cost of the research and development of new drugs. In this Perspective, we summarize our recent efforts on integrating molecular modeling and machine learning to develop computational tools for modulator design, including a pocket-guided rational design approach based on AlphaSpace to target protein-protein interactions, delta machine learning scoring functions for protein-ligand docking as well as virtual screening, and state-of-the-art deep learning models to predict calculated and experimental molecular properties based on molecular mechanics optimized geometries. Meanwhile, we discuss remaining challenges and promising directions for further development and use a retrospective example of FDA approved kinase inhibitor Erlotinib to demonstrate the use of these newly developed computational tools.
Collapse
Affiliation(s)
- Song Xia
- Department
of Chemistry, New York University, New York, New York 10003, United States
| | - Eric Chen
- Department
of Chemistry, New York University, New York, New York 10003, United States
| | - Yingkai Zhang
- Department
of Chemistry, New York University, New York, New York 10003, United States
- Simons
Center for Computational Physical Chemistry at New York University, New York, New York 10003, United States
- NYU-ECNU
Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
5
|
Millán-Pacheco C, Rios-Soto L, Corral-Rodríguez N, Sierra-Campos E, Valdez-Solana M, Téllez-Valencia A, Avitia-Domínguez C. Discovery of Potential Noncovalent Inhibitors of Dehydroquinate Dehydratase from Methicillin-Resistant Staphylococcus aureus through Computational-Driven Drug Design. Pharmaceuticals (Basel) 2023; 16:1148. [PMID: 37631063 PMCID: PMC10458038 DOI: 10.3390/ph16081148] [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: 06/30/2023] [Revised: 08/10/2023] [Accepted: 08/10/2023] [Indexed: 08/27/2023] Open
Abstract
Bacteria resistance to antibiotics is a concerning global health problem; in this context, methicillin-resistant Staphylococcus aureus (MRSA) is considered as a high priority by the World Health Organization. Furthermore, patients with a positive result for COVID-19 received early antibiotic treatment, a fact that potentially encourages the increase in antibiotic resistance. Therefore, there is an urgency to develop new drugs with molecular mechanisms different from those of the actual treatments. In this context, enzymes from the shikimate pathway, a route absent in humans, such as dehydroquinate dehydratase (DHQD), are considered good targets. In this work, a computer-aided drug design strategy, which involved exhaustive virtual screening and molecular dynamics simulations with MM-PBSA analysis, as well as an in silico ADMETox characterization, was performed to find potential noncovalent inhibitors of DHQD from MRSA (SaDHQD). After filtering the 997 million compounds from the ZINC database, 6700 compounds were submitted to an exhaustive virtual screening protocol. From these data, four molecules were selected and characterized (ZINC000005753647 (1), ZINC000001720488 (2), ZINC000082049768 (3), and ZINC000644149506 (4)). The results indicate that the four potential inhibitors interacted with residues important for substrate binding and catalysis, with an estimated binding free energy like that of the enzyme's substrate. Their ADMETox-predicted properties suggest that all of them support the structural characteristics to be considered good candidates. Therefore, the four compounds reported here are excellent option to be considered for future in vitro studies to design new SaDHQD noncovalent inhibitors and contribute to the search for new drugs against MRSA.
Collapse
Affiliation(s)
- César Millán-Pacheco
- Facultad de Farmacia, Universidad Autónoma del Estado de Morelos, Cuernavaca, Morelos 62209, Mexico;
| | - Lluvia Rios-Soto
- Facultad de Medicina y Nutrición, Universidad Juárez del Estado de Durango, Av. Universidad y Fanny Anitua S/N, Durango 34000, Mexico; (L.R.-S.); (N.C.-R.)
| | - Noé Corral-Rodríguez
- Facultad de Medicina y Nutrición, Universidad Juárez del Estado de Durango, Av. Universidad y Fanny Anitua S/N, Durango 34000, Mexico; (L.R.-S.); (N.C.-R.)
| | - Erick Sierra-Campos
- Facultad de Ciencias Químicas, Universidad Juárez del Estado de Durango Campus Gómez Palacio, Avenida Artículo 123 S/N, Fracc. Filadelfia, Gómez Palacio 35010, Mexico; (E.S.-C.); (M.V.-S.)
| | - Mónica Valdez-Solana
- Facultad de Ciencias Químicas, Universidad Juárez del Estado de Durango Campus Gómez Palacio, Avenida Artículo 123 S/N, Fracc. Filadelfia, Gómez Palacio 35010, Mexico; (E.S.-C.); (M.V.-S.)
| | - Alfredo Téllez-Valencia
- Facultad de Medicina y Nutrición, Universidad Juárez del Estado de Durango, Av. Universidad y Fanny Anitua S/N, Durango 34000, Mexico; (L.R.-S.); (N.C.-R.)
| | - Claudia Avitia-Domínguez
- Facultad de Medicina y Nutrición, Universidad Juárez del Estado de Durango, Av. Universidad y Fanny Anitua S/N, Durango 34000, Mexico; (L.R.-S.); (N.C.-R.)
| |
Collapse
|
6
|
Silvestri G, Arrigoni F, Persico F, Bertini L, Zampella G, De Gioia L, Vertemara J. Assessing the Performance of Non-Equilibrium Thermodynamic Integration in Flavodoxin Redox Potential Estimation. Molecules 2023; 28:6016. [PMID: 37630271 PMCID: PMC10459689 DOI: 10.3390/molecules28166016] [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: 06/30/2023] [Revised: 08/07/2023] [Accepted: 08/09/2023] [Indexed: 08/27/2023] Open
Abstract
Flavodoxins are enzymes that contain the redox-active flavin mononucleotide (FMN) cofactor and play a crucial role in numerous biological processes, including energy conversion and electron transfer. Since the redox characteristics of flavodoxins are significantly impacted by the molecular environment of the FMN cofactor, the evaluation of the interplay between the redox properties of the flavin cofactor and its molecular surroundings in flavoproteins is a critical area of investigation for both fundamental research and technological advancements, as the electrochemical tuning of flavoproteins is necessary for optimal interaction with redox acceptor or donor molecules. In order to facilitate the rational design of biomolecular devices, it is imperative to have access to computational tools that can accurately predict the redox potential of both natural and artificial flavoproteins. In this study, we have investigated the feasibility of using non-equilibrium thermodynamic integration protocols to reliably predict the redox potential of flavodoxins. Using as a test set the wild-type flavodoxin from Clostridium Beijerinckii and eight experimentally characterized single-point mutants, we have computed their redox potential. Our results show that 75% (6 out of 8) of the calculated reaction free energies are within 1 kcal/mol of the experimental values, and none exceed an error of 2 kcal/mol, confirming that non-equilibrium thermodynamic integration is a trustworthy tool for the quantitative estimation of the redox potential of this biologically and technologically significant class of enzymes.
Collapse
Affiliation(s)
| | | | | | | | | | - Luca De Gioia
- Department of Biotechnology and Biosciences BtBs, University of Milano-Bicocca, Piazza dell’Ateneo Nuovo 1, 20126 Milan, Italy
| | - Jacopo Vertemara
- Department of Biotechnology and Biosciences BtBs, University of Milano-Bicocca, Piazza dell’Ateneo Nuovo 1, 20126 Milan, Italy
| |
Collapse
|
7
|
Stampolaki M, Stylianakis I, Zgurskaya HI, Kolocouris A. Study of SQ109 analogs binding to mycobacterium MmpL3 transporter using MD simulations and alchemical relative binding free energy calculations. J Comput Aided Mol Des 2023; 37:245-264. [PMID: 37129848 DOI: 10.1007/s10822-023-00504-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2022] [Accepted: 04/03/2023] [Indexed: 05/03/2023]
Abstract
N-geranyl-N΄-(2-adamantyl)ethane-1,2-diamine (SQ109) is a tuberculosis drug that has high potency against Mycobacterium tuberculosis (Mtb) and may function by blocking cell wall biosynthesis. After the crystal structure of MmpL3 from Mycobacterium smegmatis in complex with SQ109 became available, it was suggested that SQ109 inhibits Mmpl3 mycolic acid transporter. Here, we showed using molecular dynamics (MD) simulations that the binding profile of nine SQ109 analogs with inhibitory potency against Mtb and alkyl or aryl adducts at C-2 or C-1 adamantyl carbon to MmpL3 was consistent with the X-ray structure of MmpL3 - SQ109 complex. We showed that rotation of SQ109 around carbon-carbon bond in the monoprotonated ethylenediamine unit favors two gauche conformations as minima in water and lipophilic solvent using DFT calculations as well as inside the transporter's binding area using MD simulations. The binding assays in micelles suggested that the binding affinity of the SQ109 analogs was increased for the larger, more hydrophobic adducts, which was consistent with our results from MD simulations of the SQ109 analogues suggesting that sizeable C-2 adamantyl adducts of SQ109 can fill a lipophilic region between Y257, Y646, F260 and F649 in MmpL3. This was confirmed quantitatively by our calculations of the relative binding free energies using the thermodynamic integration coupled with MD simulations method with a mean assigned error of 0.74 kcal mol-1 compared to the experimental values.
Collapse
Affiliation(s)
- Marianna Stampolaki
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771, Athens, Greece
- Department of NMR-Based Structural Biology, Max Planck Institute for Multidisciplinary Sciences, Am Faßberg 11, 37077, Göttingen, Germany
| | - Ioannis Stylianakis
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771, Athens, Greece
| | - Helen I Zgurskaya
- Department of Chemistry and Biochemistry, University of Oklahoma, Stephenson Life Sciences Research Center, 101 Stephenson Parkway, Norman, OK, 73019-5251, USA
| | - Antonios Kolocouris
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771, Athens, Greece.
| |
Collapse
|
8
|
Wang S, Sun X, Cui W, Yuan S. MM/PB(GB)SA benchmarks on soluble proteins and membrane proteins. Front Pharmacol 2022; 13:1018351. [PMID: 36532746 PMCID: PMC9751045 DOI: 10.3389/fphar.2022.1018351] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2022] [Accepted: 11/17/2022] [Indexed: 08/27/2023] Open
Abstract
Predicting protein-ligand binding free energy rapidly and accurately remains a challenging question in modern drug discovery. Molecular mechanics/Poisson-Boltzmann (Generalized Born) surface area (MM/PB(GB)SA) has emerged as an essential tool for accelerating cost-efficient binding free energy calculation. This study presents benchmarks with three membrane-bound protein systems and six soluble protein systems. Different parameters were sampled for different benchmarks to explore the highest accuracy. These include ligand charges, protein force fields, extra points, GB models, nonpolar optimization methods, internal dielectric constants and membrane dielectric constants. Comparisons of accuracy were made between MM/PB(GB)SA, docking and free energy perturbation (FEP). The results reveal a competitive performance between MM/PB(GB)SA and FEP. In summary, MM/PB(GB)SA is a powerful approach to predict ligand binding free energy rapidly and accurately. Parameters of MM/PB(GB)SA calculations, such as the GB models and membrane dielectric constants, need to be optimized for different systems. This method can be served as a powerful tool for drug design.
Collapse
Affiliation(s)
- Shiyu Wang
- Research Center for Computer-Aided Drug Discovery, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China
- College of Chemical Science, University of Chinese Academy of Sciences, Beijing, China
- AlphaMol-SIAT Joint Laboratory, Shenzhen, China
| | - Xiaolin Sun
- Research Center for Computer-Aided Drug Discovery, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China
- AlphaMol-SIAT Joint Laboratory, Shenzhen, China
| | - Wenqiang Cui
- Research Center for Computer-Aided Drug Discovery, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China
- AlphaMol-SIAT Joint Laboratory, Shenzhen, China
| | - Shuguang Yuan
- Research Center for Computer-Aided Drug Discovery, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, China
- AlphaMol-SIAT Joint Laboratory, Shenzhen, China
- Faculty of Chemistry, University of Warsaw, Warsaw, Poland
| |
Collapse
|
9
|
Stampelou M, Suchankova A, Tzortzini E, Dhingra L, Barkan K, Lougiakis N, Marakos P, Pouli N, Ladds G, Kolocouris A. Dual A1/A3 Adenosine Receptor Antagonists: Binding Kinetics and Structure-Activity Relationship Studies Using Mutagenesis and Alchemical Binding Free Energy Calculations. J Med Chem 2022; 65:13305-13327. [PMID: 36173355 DOI: 10.1021/acs.jmedchem.2c01123] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Drugs targeting adenosine receptors (AR) can provide treatment for diseases. We report the identification of 7-(phenylamino)-pyrazolo[3,4-c]pyridines L2-L10, A15, and A17 as low-micromolar to low-nanomolar A1R/A3R dual antagonists, with 3-phenyl-5-cyano-7-(trimethoxyphenylamino)-pyrazolo[3,4-c]pyridine (A17) displaying the highest affinity at both receptors with a long residence time of binding, as determined using a NanoBRET-based assay. Two binding orientations of A17 produce stable complexes inside the orthosteric binding area of A1R in molecular dynamics (MD) simulations, and we selected the most plausible orientation based on the agreement with alanine mutagenesis supported by affinity experiments. Interestingly, for drug design purposes, the mutation of L2506.51 to alanine increased the binding affinity of A17 at A1R. We explored the structure-activity relationships against A1R using alchemical binding free energy calculations with the thermodynamic integration coupled with the MD simulation (TI/MD) method, applied on the whole G-protein-coupled receptor-membrane system, which showed a good agreement (r = 0.73) between calculated and experimental relative binding free energies.
Collapse
Affiliation(s)
- Margarita Stampelou
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771 Athens, Greece
| | - Anna Suchankova
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, U.K
| | - Efpraxia Tzortzini
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771 Athens, Greece
| | - Lakshiv Dhingra
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, U.K
| | - Kerry Barkan
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, U.K
| | - Nikolaos Lougiakis
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771 Athens, Greece
| | - Panagiotis Marakos
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771 Athens, Greece
| | - Nicole Pouli
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771 Athens, Greece
| | - Graham Ladds
- Department of Pharmacology, University of Cambridge, Tennis Court Road, Cambridge CB2 1PD, U.K
| | - Antonios Kolocouris
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, School of Health Sciences, National and Kapodistrian University of Athens, Panepistimiopolis-Zografou, 15771 Athens, Greece
| |
Collapse
|
10
|
Yang C, Chen EA, Zhang Y. Protein-Ligand Docking in the Machine-Learning Era. Molecules 2022; 27:4568. [PMID: 35889440 PMCID: PMC9323102 DOI: 10.3390/molecules27144568] [Citation(s) in RCA: 32] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2022] [Accepted: 07/14/2022] [Indexed: 11/16/2022] Open
Abstract
Molecular docking plays a significant role in early-stage drug discovery, from structure-based virtual screening (VS) to hit-to-lead optimization, and its capability and predictive power is critically dependent on the protein-ligand scoring function. In this review, we give a broad overview of recent scoring function development, as well as the docking-based applications in drug discovery. We outline the strategies and resources available for structure-based VS and discuss the assessment and development of classical and machine learning protein-ligand scoring functions. In particular, we highlight the recent progress of machine learning scoring function ranging from descriptor-based models to deep learning approaches. We also discuss the general workflow and docking protocols of structure-based VS, such as structure preparation, binding site detection, docking strategies, and post-docking filter/re-scoring, as well as a case study on the large-scale docking-based VS test on the LIT-PCBA data set.
Collapse
Affiliation(s)
- Chao Yang
- Department of Chemistry, New York University, New York, NY 10003, USA; (C.Y.); (E.A.C.)
| | - Eric Anthony Chen
- Department of Chemistry, New York University, New York, NY 10003, USA; (C.Y.); (E.A.C.)
| | - Yingkai Zhang
- Department of Chemistry, New York University, New York, NY 10003, USA; (C.Y.); (E.A.C.)
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| |
Collapse
|
11
|
Suruzhon M, Bodnarchuk MS, Ciancetta A, Viner R, Wall ID, Essex JW. Sensitivity of Binding Free Energy Calculations to Initial Protein Crystal Structure. J Chem Theory Comput 2021; 17:1806-1821. [PMID: 33534995 DOI: 10.1021/acs.jctc.0c00972] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Binding free energy calculations using alchemical free energy (AFE) methods are widely considered to be the most rigorous tool in the computational drug discovery arsenal. Despite this, the calculations suffer from accuracy, precision, and reproducibility issues. In this publication, we perform a high-throughput study of more than a thousand AFE calculations, utilizing over 220 μs of total sampling time, on three different protein systems to investigate the impact of the initial crystal structure on the resulting binding free energy values. We also consider the influence of equilibration time and discover that the initial crystal structure can have a significant effect on free energy values obtained at short timescales that can manifest itself as a free energy difference of more than 1 kcal/mol. At longer timescales, these differences are largely overtaken by important rare events, such as torsional ligand motions, typically resulting in a much higher uncertainty in the obtained values. This work emphasizes the importance of rare event sampling and long-timescale dynamics in free energy calculations even for routinely performed alchemical perturbations. We conclude that an optimal protocol should not only concentrate computational resources on achieving convergence in the alchemical coupling parameter (λ) space but also on longer simulations and multiple repeats.
Collapse
Affiliation(s)
- Miroslav Suruzhon
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, U.K
| | | | | | - Russell Viner
- Syngenta, Jealott's Hill International Research Centre, Bracknell RG42 6EY, U.K
| | - Ian D Wall
- GSK Medicines Research Centre, Gunnels Wood Road, Stevenage SG1 2NY, U.K
| | - Jonathan W Essex
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, U.K
| |
Collapse
|
12
|
Wang X, Li W. Comparative Study of Interactions between Human cGAS and Inhibitors: Insights from Molecular Dynamics and MM/PBSA Studies. Int J Mol Sci 2021; 22:ijms22031164. [PMID: 33503880 PMCID: PMC7865699 DOI: 10.3390/ijms22031164] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2020] [Revised: 01/20/2021] [Accepted: 01/21/2021] [Indexed: 11/20/2022] Open
Abstract
Recent studies have identified cyclic GMP-AMP synthase (cGAS) as an important target for treating autoimmune diseases, and several inhibitors of human cGAS (hcGAS) and their structures in complexation with hcGAS have been reported. However, the mechanisms via which these inhibitors interact with hcGAS are not completely understood. Here, we aimed to assess the performance of molecular mechanics/Poisson–Boltzmann solvent-accessible surface area (MM/PBSA) in evaluating the binding affinity of various hcGAS inhibitors and to elucidate their detailed interactions with hcGAS from an energetic viewpoint. Using molecular dynamics (MD) simulation and MM/PBSA approaches, the estimated free energies were in good agreement with the experimental ones, with a Pearson’s correlation coefficient and Spearman’s rank coefficient of 0.67 and 0.46, respectively. In per-residue energy decomposition analysis, four residues, K362, R376, Y436, and K439 in hcGAS were found to contribute significantly to the binding with inhibitors via hydrogen bonding, salt bridges, and various π interactions, such as π· · ·π stacking, cation· · ·π, hydroxyl· · ·π, and alkyl· · ·π interactions. In addition, we discussed other key interactions between specific residues and ligands, in particular, between H363 and JUJ, F379 and 9BY, and H437 and 8ZM. The sandwiched structures of the inhibitor bound to the guanidinium group of R376 and the phenyl ring of Y436 were also consistent with the experimental data. The results indicated that MM/PBSA in combination with other virtual screening methods, could be a reliable approach to discover new hcGAS inhibitors and thus is valuable for potential treatments of cGAS-dependent inflammatory diseases.
Collapse
Affiliation(s)
- Xiaowen Wang
- Institute for Advanced Study, Shenzhen University, Shenzhen 518060, China;
- College of Physics and Optoelectronic Engineering, Shenzhen University, Shenzhen 518060, China
| | - Wenjin Li
- Institute for Advanced Study, Shenzhen University, Shenzhen 518060, China;
- Correspondence: ; Tel.: +86-755-2694-2336
| |
Collapse
|
13
|
Decherchi S, Cavalli A. Thermodynamics and Kinetics of Drug-Target Binding by Molecular Simulation. Chem Rev 2020; 120:12788-12833. [PMID: 33006893 PMCID: PMC8011912 DOI: 10.1021/acs.chemrev.0c00534] [Citation(s) in RCA: 125] [Impact Index Per Article: 31.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2020] [Indexed: 12/19/2022]
Abstract
Computational studies play an increasingly important role in chemistry and biophysics, mainly thanks to improvements in hardware and algorithms. In drug discovery and development, computational studies can reduce the costs and risks of bringing a new medicine to market. Computational simulations are mainly used to optimize promising new compounds by estimating their binding affinity to proteins. This is challenging due to the complexity of the simulated system. To assess the present and future value of simulation for drug discovery, we review key applications of advanced methods for sampling complex free-energy landscapes at near nonergodicity conditions and for estimating the rate coefficients of very slow processes of pharmacological interest. We outline the statistical mechanics and computational background behind this research, including methods such as steered molecular dynamics and metadynamics. We review recent applications to pharmacology and drug discovery and discuss possible guidelines for the practitioner. Recent trends in machine learning are also briefly discussed. Thanks to the rapid development of methods for characterizing and quantifying rare events, simulation's role in drug discovery is likely to expand, making it a valuable complement to experimental and clinical approaches.
Collapse
Affiliation(s)
- Sergio Decherchi
- Computational
and Chemical Biology, Fondazione Istituto
Italiano di Tecnologia, 16163 Genoa, Italy
| | - Andrea Cavalli
- Computational
and Chemical Biology, Fondazione Istituto
Italiano di Tecnologia, 16163 Genoa, Italy
- Department
of Pharmacy and Biotechnology, University
of Bologna, 40126 Bologna, Italy
| |
Collapse
|
14
|
Liu K, Kokubo H. Prediction of ligand binding mode among multiple cross-docking poses by molecular dynamics simulations. J Comput Aided Mol Des 2020; 34:1195-1205. [DOI: 10.1007/s10822-020-00340-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2019] [Accepted: 08/21/2020] [Indexed: 01/18/2023]
|
15
|
Nakajima R, Midorikawa N. Topic extraction to provide an overview of research activities: The case of the high-temperature superconductor and simulation and modelling. J Inf Sci 2020. [DOI: 10.1177/0165551520920794] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
For those who are not experts in a particular scientific field, it is difficult to understand scientific research trends. Although studies on the extraction of research trends have been conducted, most focus on extracting global trends from large-scale data, and the methods are often complicated. The purpose of this study is to develop a method of obtaining overviews of a scientific field for non-experts by capturing research trends simply and then to verify the method. To extract research topics which should express research trends, text analysis was performed using abstracts over 12 years of articles on high-temperature superconductors. We characterised three topics for the extracted word groups that frequently occurred. For these topics, we studied their appropriateness using a method that has been little used: examining research articles, review literature and co-citations among research articles used to extract the words, comparisons with controlled index terms assigned to the articles and confirming that there were no contradictions. Based on the established method, we have also applied this method to another research field: ‘simulation and modelling’. Although the method used in this article is simple, important topics were extracted, and the relations with the original articles are clear, which can lead to further investigation of the extracted topics.
Collapse
Affiliation(s)
- Ritsuko Nakajima
- Graduate School of Library, Information and Media Studies, University of Tsukuba, Japan
| | - Nobuyuki Midorikawa
- Faculty of Library, Information and Media Science, University of Tsukuba, Japan
| |
Collapse
|
16
|
Eberlein L, Beierlein FR, van Eikema Hommes NJR, Radadiya A, Heil J, Benner SA, Clark T, Kast SM, Richards NGJ. Tautomeric Equilibria of Nucleobases in the Hachimoji Expanded Genetic Alphabet. J Chem Theory Comput 2020; 16:2766-2777. [PMID: 32125859 DOI: 10.1021/acs.jctc.9b01079] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
Evolution has yielded biopolymers that are constructed from exactly four building blocks and are able to support Darwinian evolution. Synthetic biology aims to extend this alphabet, and we recently showed that 8-letter (hachimoji) DNA can support rule-based information encoding. One source of replicative error in non-natural DNA-like systems, however, is the occurrence of alternative tautomeric forms, which pair differently. Unfortunately, little is known about how structural modifications impact free-energy differences between tautomers of the non-natural nucleobases used in the hachimoji expanded genetic alphabet. Determining experimental tautomer ratios is technically difficult, and so, strategies for improving hachimoji DNA replication efficiency will benefit from accurate computational predictions of equilibrium tautomeric ratios. We now report that high-level quantum-chemical calculations in aqueous solution by the embedded cluster reference interaction site model, benchmarked against free-energy molecular simulations for solvation thermodynamics, provide useful quantitative information on the tautomer ratios of both Watson-Crick and hachimoji nucleobases. In agreement with previous computational studies, all four Watson-Crick nucleobases adopt essentially only one tautomer in water. This is not the case, however, for non-natural nucleobases and their analogues. For example, although the enols of isoguanine and a series of related purines are not populated in water, these heterocycles possess N1-H and N3-H keto tautomers that are similar in energy, thereby adversely impacting accurate nucleobase pairing. These robust computational strategies offer a firm basis for improving experimental measurements of tautomeric ratios, which are currently limited to studying molecules that exist only as two tautomers in solution.
Collapse
Affiliation(s)
- Lukas Eberlein
- Physikalische Chemie III, Technische Universität Dortmund, Dortmund 44227, Germany
| | - Frank R Beierlein
- Computer-Chemistry-Centre and Interdisciplinary Centre for Molecular Materials, Department of Chemistry & Pharmacy, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen 91054, Germany
| | - Nico J R van Eikema Hommes
- Computer-Chemistry-Centre and Interdisciplinary Centre for Molecular Materials, Department of Chemistry & Pharmacy, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen 91054, Germany
| | - Ashish Radadiya
- School of Chemistry, Cardiff University, Cardiff CF10 3AT, U.K
| | - Jochen Heil
- Physikalische Chemie III, Technische Universität Dortmund, Dortmund 44227, Germany
| | - Steven A Benner
- Foundation for Applied Molecular Evolution, Alachua, Florida 32615, United States
| | - Timothy Clark
- Computer-Chemistry-Centre and Interdisciplinary Centre for Molecular Materials, Department of Chemistry & Pharmacy, Friedrich-Alexander-Universität Erlangen-Nürnberg, Erlangen 91054, Germany
| | - Stefan M Kast
- Physikalische Chemie III, Technische Universität Dortmund, Dortmund 44227, Germany
| | - Nigel G J Richards
- School of Chemistry, Cardiff University, Cardiff CF10 3AT, U.K.,Foundation for Applied Molecular Evolution, Alachua, Florida 32615, United States
| |
Collapse
|
17
|
Cong Y, Duan L, Huang K, Bao J, Zhang JZH. Alanine scanning combined with interaction entropy studying the differences of binding mechanism on HIV-1 and HIV-2 proteases with inhibitor. J Biomol Struct Dyn 2020; 39:1588-1599. [DOI: 10.1080/07391102.2020.1734488] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Yalong Cong
- School of Physics and Electronics, Shandong Normal University, Jinan, China
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, China
| | - Lili Duan
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Kaifang Huang
- School of Physics and Electronics, Shandong Normal University, Jinan, China
| | - Jinxiao Bao
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, China
| | - John Z. H. Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, China
- Department of Chemistry, New York University, NY, NY, USA
| |
Collapse
|
18
|
Leonhard AC, Whitmer JK. Accurate Determination of Cavitand Binding Free Energies via Unrestrained Advanced Sampling. J Chem Theory Comput 2019; 15:5761-5768. [DOI: 10.1021/acs.jctc.9b00348] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Anne C. Leonhard
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana 46556, United States
| | - Jonathan K. Whitmer
- Department of Chemical and Biomolecular Engineering, University of Notre Dame, Notre Dame, Indiana 46556, United States
| |
Collapse
|
19
|
Shinkafi TS, Kaushik A, Mahmood A, Tiwari AK, Alam MM, Akhter M, Gupta D, Ali S. Computational prediction and experimental validation of the activator function of C2-β-D-glucopyranosyl-1,3,6,7-tetrahydroxyxanthone on pancreatic and hepatic hexokinase. J Biomol Struct Dyn 2019; 38:2976-2987. [DOI: 10.1080/07391102.2019.1650829] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Tijjani Salihu Shinkafi
- Department of Biochemistry, School of Chemical and Life Sciences, Jamia Hamdard (Deemed University), Hamdard Nagar, New Delhi, India
- Department of Biochemistry, Faculty of Science, Usmanu Danfodiyo University, Sokoto, Nigeria
| | - Abhinav Kaushik
- Bioinformatics Infrastructure Facility, Jamia Hamdard, New Delhi, India
- International Center for Genetic Engineering and Biotechnology, New Delhi, India
| | - Amena Mahmood
- Department of Biochemistry, School of Chemical and Life Sciences, Jamia Hamdard (Deemed University), Hamdard Nagar, New Delhi, India
| | | | - Mohammad Mumtaz Alam
- School of Pharmaceutical Education and Research, Jamia Hamdard, New Delhi, India
| | - Mymoona Akhter
- Bioinformatics Infrastructure Facility, Jamia Hamdard, New Delhi, India
- School of Pharmaceutical Education and Research, Jamia Hamdard, New Delhi, India
| | - Dinesh Gupta
- International Center for Genetic Engineering and Biotechnology, New Delhi, India
| | - Shakir Ali
- Department of Biochemistry, School of Chemical and Life Sciences, Jamia Hamdard (Deemed University), Hamdard Nagar, New Delhi, India
- Bioinformatics Infrastructure Facility, Jamia Hamdard, New Delhi, India
| |
Collapse
|
20
|
Berishvili VP, Perkin VO, Voronkov AE, Radchenko EV, Syed R, Venkata Ramana Reddy C, Pillay V, Kumar P, Choonara YE, Kamal A, Palyulin VA. Time-Domain Analysis of Molecular Dynamics Trajectories Using Deep Neural Networks: Application to Activity Ranking of Tankyrase Inhibitors. J Chem Inf Model 2019; 59:3519-3532. [PMID: 31276400 DOI: 10.1021/acs.jcim.9b00135] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Molecular dynamics simulations provide valuable insights into the behavior of molecular systems. Extending the recent trend of using machine learning techniques to predict physicochemical properties from molecular dynamics data, we propose to consider the trajectories as multidimensional time series represented by 2D tensors containing the ligand-protein interaction descriptor values for each time step. Similar in structure to the time series encountered in modern approaches for signal, speech, and natural language processing, these time series can be directly analyzed using long short-term memory (LSTM) recurrent neural networks or convolutional neural networks (CNNs). The predictive regression models for the ligand-protein affinity were built for a subset of the PDBbind v.2017 database and applied to inhibitors of tankyrase, an enzyme of the poly(ADP-ribose)-polymerase (PARP) family that can be used in the treatment of colorectal cancer. As an additional test set, a subset of the Community Structure-Activity Resource (CSAR) data set was used. For comparison, the random forest and simple neural network models based on the crystal pose or the trajectory-averaged descriptors were used, as well as the commonly employed docking and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) scores. Convolutional neural networks based on the 2D tensors of ligand-protein interaction descriptors for short (2 ns) trajectories provide the best accuracy and predictive power, reaching the Spearman rank correlation coefficient of 0.73 and Pearson correlation coefficient of 0.70 for the tankyrase test set. Taking into account the recent increase in computational power of modern GPUs and relatively low computational complexity of the proposed approach, it can be used as an advanced virtual screening filter for compound prioritization.
Collapse
Affiliation(s)
- Vladimir P Berishvili
- Department of Chemistry , Lomonosov Moscow State University , Moscow 119991 , Russia
| | - Valentin O Perkin
- Department of Chemistry , Lomonosov Moscow State University , Moscow 119991 , Russia
| | - Andrew E Voronkov
- Department of Chemistry , Lomonosov Moscow State University , Moscow 119991 , Russia.,Digital BioPharm Ltd. , Hovseterveien 42 A, H0301 , Oslo 0768 , Norway
| | - Eugene V Radchenko
- Department of Chemistry , Lomonosov Moscow State University , Moscow 119991 , Russia
| | - Riyaz Syed
- Department of Chemistry , Jawaharlal Nehru Technological University , Kukatpally, Hyderabad 500 085 , India
| | | | - Viness Pillay
- Wits Advanced Drug Delivery Platform Research Unit, Faculty of Health Sciences, School of Therapeutic Sciences, Department of Pharmacy and Pharmacology , University of the Witwatersrand, Johannesburg , 7 York Road , Parktown 2193 , South Africa
| | - Pradeep Kumar
- Wits Advanced Drug Delivery Platform Research Unit, Faculty of Health Sciences, School of Therapeutic Sciences, Department of Pharmacy and Pharmacology , University of the Witwatersrand, Johannesburg , 7 York Road , Parktown 2193 , South Africa
| | - Yahya E Choonara
- Wits Advanced Drug Delivery Platform Research Unit, Faculty of Health Sciences, School of Therapeutic Sciences, Department of Pharmacy and Pharmacology , University of the Witwatersrand, Johannesburg , 7 York Road , Parktown 2193 , South Africa
| | - Ahmed Kamal
- School of Pharmaceutical Education and Research , Jamia Hamdard , New Delhi 110 062 , India
| | - Vladimir A Palyulin
- Department of Chemistry , Lomonosov Moscow State University , Moscow 119991 , Russia
| |
Collapse
|
21
|
Wang E, Sun H, Wang J, Wang Z, Liu H, Zhang JZH, Hou T. End-Point Binding Free Energy Calculation with MM/PBSA and MM/GBSA: Strategies and Applications in Drug Design. Chem Rev 2019; 119:9478-9508. [DOI: 10.1021/acs.chemrev.9b00055] [Citation(s) in RCA: 578] [Impact Index Per Article: 115.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Ercheng Wang
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Huiyong Sun
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Junmei Wang
- Department of Pharmaceutical Sciences, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Zhe Wang
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - Hui Liu
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| | - John Z. H. Zhang
- Shanghai Engineering Research Center of Molecular Therapeutics & New Drug Development, Shanghai Key Laboratory of Green Chemistry & Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- NYU−ECNU Center for Computational Chemistry, NYU Shanghai, Shanghai 200122, China
- Department of Chemistry, New York University, New York, New York 10003, United States
- Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| | - Tingjun Hou
- Hangzhou Institute of Innovative Medicine, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou, Zhejiang 310058, China
| |
Collapse
|
22
|
Du Q, Qian Y, Yao X, Xue W. Elucidating the tight-binding mechanism of two oral anticoagulants to factor Xa by using induced-fit docking and molecular dynamics simulation. J Biomol Struct Dyn 2019; 38:625-633. [PMID: 30806177 DOI: 10.1080/07391102.2019.1583605] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Qingqing Du
- Department of Pharmacy, the Second Affiliated Hospital of Chongqing Medical University, Chongqing, China
| | - Yan Qian
- Department of Pharmacy, the Second Affiliated Hospital of Chongqing Medical University, Chongqing, China
| | - Xiaojun Yao
- State Key Laboratory of Applied Organic Chemistry and Department of Chemistry, Lanzhou University, Lanzhou, China
| | - Weiwei Xue
- School of Pharmaceutical Sciences and Chongqing Key Laboratory of Natural Drug Research, Chongqing University, Chongqing, China
| |
Collapse
|
23
|
Caldararu O, Olsson MA, Misini Ignjatović M, Wang M, Ryde U. Binding free energies in the SAMPL6 octa-acid host-guest challenge calculated with MM and QM methods. J Comput Aided Mol Des 2018; 32:1027-1046. [PMID: 30203229 DOI: 10.1007/s10822-018-0158-2] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2018] [Accepted: 08/31/2018] [Indexed: 01/15/2023]
Abstract
We have estimated free energies for the binding of eight carboxylate ligands to two variants of the octa-acid deep-cavity host in the SAMPL6 blind-test challenge (with or without endo methyl groups on the four upper-rim benzoate groups, OAM and OAH, respectively). We employed free-energy perturbation (FEP) for relative binding energies at the molecular mechanics (MM) and the combined quantum mechanical (QM) and MM (QM/MM) levels, the latter obtained with the reference-potential approach with QM/MM sampling for the MM → QM/MM FEP. The semiempirical QM method PM6-DH+ was employed for the ligand in the latter calculations. Moreover, binding free energies were also estimated from QM/MM optimised structures, combined with COSMO-RS estimates of the solvation energy and thermostatistical corrections from MM frequencies. They were performed at the PM6-DH+ level of theory with the full host and guest molecule in the QM system (and also four water molecules in the geometry optimisations) for 10-20 snapshots from molecular dynamics simulations of the complex. Finally, the structure with the lowest free energy was recalculated using the dispersion-corrected density-functional theory method TPSS-D3, for both the structure and the energy. The two FEP approaches gave similar results (PM6-DH+/MM slightly better for OAM), which were among the five submissions with the best performance in the challenge and gave the best results without any fit to data from the SAMPL5 challenge, with mean absolute deviations (MAD) of 2.4-5.2 kJ/mol and a correlation coefficient (R2) of 0.77-0.93. This is the first time QM/MM approaches give binding free energies that are competitive to those obtained with MM for the octa-acid host. The QM/MM-optimised structures gave somewhat worse performance (MAD = 3-8 kJ/mol and R2 = 0.1-0.9), but the results were improved compared to previous studies of this system with similar methods.
Collapse
Affiliation(s)
- Octav Caldararu
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Martin A Olsson
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Majda Misini Ignjatović
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Meiting Wang
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden.
| |
Collapse
|
24
|
Salawu EO. The Impairment of TorsinA's Binding to and Interactions With Its Activator: An Atomistic Molecular Dynamics Study of Primary Dystonia. Front Mol Biosci 2018; 5:64. [PMID: 30042949 PMCID: PMC6048259 DOI: 10.3389/fmolb.2018.00064] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2018] [Accepted: 06/19/2018] [Indexed: 01/23/2023] Open
Abstract
Primary dystonia's prolonged muscle contractions and the associated abnormal postures and twisting movements remain incurable. Genetic mutation/deletion of GAG from TorsonA's gene resulting in ΔE303 (which weakens the binding between TorsinA and its activator, such as LULL1) primarily cause this neurodegenerative disorder. We studied TorsinA-LULL1 (or TorsinAΔE303-LULL1) bindings and interactions. For the first time, we show the atomic details of TorsinA-LULL1 dynamic interactions and TorsinAΔE303-LULL1 dynamic interactions and their binding affinities. Our results show extensive effects of ΔE303 on TorsinAΔE303-LULL1 interactions, and suggest that the differences between TorsinA-LULL1 interactions and TorsinAΔE303-LULL1 interactions are non-subtle. ΔE303 significantly weakens TorsinAΔE303-LULL1's binding affinity. We present pieces of evidence proving that the effects of ΔE303 (on the differences between TorsinA-LULL1 interactions and TorsinAΔE303-LULL1 interactions) are more pronounced than previously suggested, and that the nanobody used for achieving the X-ray crystallization in the previous study attenuated the differences between TorsinA-LULL1 and TorsinAΔE303-LULL1 interactions. Our accounts of the dynamic interactions between “TorsinA and LULL1” and between “TorsinAΔE303 and LULL1” and the detailed effects of ΔE303 on TorsinA-/TorsinAΔE303-LULL1 build on previous findings and offer new insights for a better understanding of the molecular basis of Primary Dystonia. Our results have long-term potentials of guiding the development of medications for the disease.
Collapse
Affiliation(s)
- Emmanuel O Salawu
- TIGP Bioinformatics Program, Academia Sinica, Taipei, Taiwan.,Institute of Bioinformatics and Structural Biology, National Tsing Hua University, Hsinchu, Taiwan.,School of Computer Science, University of Hertfordshire, Hertfordshire, United Kingdom.,Bioinformatics Center, Sheridan, WY, United States
| |
Collapse
|
25
|
Bhati AP, Wan S, Hu Y, Sherborne B, Coveney PV. Uncertainty Quantification in Alchemical Free Energy Methods. J Chem Theory Comput 2018; 14:2867-2880. [PMID: 29678106 PMCID: PMC6095638 DOI: 10.1021/acs.jctc.7b01143] [Citation(s) in RCA: 43] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Alchemical
free energy methods have gained much importance recently
from several reports of improved ligand–protein binding affinity
predictions based on their implementation using molecular dynamics
simulations. A large number of variants of such methods implementing
different accelerated sampling techniques and free energy estimators
are available, each claimed to be better than the others in its own
way. However, the key features of reproducibility and quantification
of associated uncertainties in such methods have barely been discussed.
Here, we apply a systematic protocol for uncertainty quantification
to a number of popular alchemical free energy methods, covering both
absolute and relative free energy predictions. We show that a reliable
measure of error estimation is provided by ensemble simulation—an
ensemble of independent MD simulations—which applies irrespective
of the free energy method. The need to use ensemble methods is fundamental
and holds regardless of the duration of time of the molecular dynamics
simulations performed.
Collapse
Affiliation(s)
- Agastya P Bhati
- Centre for Computational Science, Department of Chemistry , University College London , 20 Gordon Street , London WC1H 0AJ , United Kingdom
| | - Shunzhou Wan
- Centre for Computational Science, Department of Chemistry , University College London , 20 Gordon Street , London WC1H 0AJ , United Kingdom
| | - Yuan Hu
- Modeling and Informatics , Merck & Co., Inc. , 2000 Galloping Hill Road , Kenilworth , New Jersey 07033 , United States
| | - Brad Sherborne
- Modeling and Informatics , Merck & Co., Inc. , 2000 Galloping Hill Road , Kenilworth , New Jersey 07033 , United States
| | - Peter V Coveney
- Centre for Computational Science, Department of Chemistry , University College London , 20 Gordon Street , London WC1H 0AJ , United Kingdom
| |
Collapse
|
26
|
Ren J, Mistry TL, Su PC, Mehboob S, Demissie R, Fung LWM, Ghosh AK, Johnson ME. Determination of absolute configuration and binding efficacy of benzimidazole-based FabI inhibitors through the support of electronic circular dichroism and MM-GBSA techniques. Bioorg Med Chem Lett 2018; 28:2074-2079. [PMID: 29730028 DOI: 10.1016/j.bmcl.2018.04.052] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2018] [Revised: 04/19/2018] [Accepted: 04/20/2018] [Indexed: 01/12/2023]
Abstract
We have previously reported benzimidazole-based compounds to be potent inhibitors of FabI for Francisella tularensis (FtFabI), making them promising antimicrobial hits. Optically active enantiomers exhibit markedly differing affinities toward FtFabI. The IC50 of benzimidazole (-)-1 is ∼100× lower than the (+)-enantiomer, with similar results for the 2 enantiomers. Determining the absolute configuration for these optical compounds and elucidating their binding modes is important for further design. Electronic circular dichroism (ECD) quantum calculations have become important in determining absolute configurations of optical compounds. We determined the absolute configuration of (-)/(+)-1 and (-)/(+)-2 by comparing experimental spectra and theoretical density functional theory (DFT) simulations of ECD spectra at the B3LYP/6-311+G(2d, p) level using Gaussian09. Comparison of experimental and calculated ECD spectra indicates that the S configuration corresponds to the (-)-rotation for both compounds 1 and 2, while the R configuration corresponds to the (+)-rotation. Further, molecular dynamics simulations and MM-GBSA binding energy calculations for these two pairs of enantiomers with FtFabI show much tighter binding MM-GBSA free energies for S-1 and S-2 than for their enantiomers, R-1 and R-2, consistent with the S configuration being the more active one, and with the ECD determination of the S configuration corresponding to (-) and the R configuration corresponding to (+). Thus, our computational studies allow us to assign (-) to (S)- and (+) to (R)- for compounds 1 and 2, and to further evaluate structural changes to improve efficacy.
Collapse
Affiliation(s)
- Jinhong Ren
- Center for Biomolecular Sciences, University of Illinois at Chicago, 900 S. Ashland Ave, Chicago, IL 60607, USA
| | - Tina L Mistry
- Center for Biomolecular Sciences, University of Illinois at Chicago, 900 S. Ashland Ave, Chicago, IL 60607, USA
| | - Pin-Chih Su
- Center for Biomolecular Sciences, University of Illinois at Chicago, 900 S. Ashland Ave, Chicago, IL 60607, USA
| | - Shahila Mehboob
- Novalex Therapeutics, Inc., 2242 W Harrison, Chicago, IL 60612, USA
| | - Robel Demissie
- Department of Chemistry, University of Illinois at Chicago, 845 W. Taylor St, Chicago, IL 60607, USA
| | - Leslie Wo-Mei Fung
- Department of Chemistry, University of Illinois at Chicago, 845 W. Taylor St, Chicago, IL 60607, USA
| | - Arun K Ghosh
- Department of Chemistry and Department of Medicinal Chemistry, Purdue University, 560 Oval Drive, West Lafayette, IN 47907, USA
| | - Michael E Johnson
- Center for Biomolecular Sciences, University of Illinois at Chicago, 900 S. Ashland Ave, Chicago, IL 60607, USA; Novalex Therapeutics, Inc., 2242 W Harrison, Chicago, IL 60612, USA.
| |
Collapse
|
27
|
Morrison EA, Bowerman S, Sylvers KL, Wereszczynski J, Musselman CA. The conformation of the histone H3 tail inhibits association of the BPTF PHD finger with the nucleosome. eLife 2018; 7:31481. [PMID: 29648537 PMCID: PMC5953545 DOI: 10.7554/elife.31481] [Citation(s) in RCA: 92] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2017] [Accepted: 04/11/2018] [Indexed: 01/08/2023] Open
Abstract
Histone tails harbor a plethora of post-translational modifications that direct the function of chromatin regulators, which recognize them through effector domains. Effector domain/histone interactions have been broadly studied, but largely using peptide fragments of histone tails. Here, we extend these studies into the nucleosome context and find that the conformation adopted by the histone H3 tails is inhibitory to BPTF PHD finger binding. Using NMR spectroscopy and MD simulations, we show that the H3 tails interact robustly but dynamically with nucleosomal DNA, substantially reducing PHD finger association. Altering the electrostatics of the H3 tail via modification or mutation increases accessibility to the PHD finger, indicating that PTM crosstalk can regulate effector domain binding by altering nucleosome conformation. Together, our results demonstrate that the nucleosome context has a dramatic impact on signaling events at the histone tails, and highlights the importance of studying histone binding in the context of the nucleosome.
Collapse
Affiliation(s)
- Emma A Morrison
- Department of Biochemistry, Carver College of Medicine, University of Iowa, Iowa City, United States
| | - Samuel Bowerman
- Department of Physics, Illinois Institute of Technology, Chicago, Illinois.,Center for Molecular Study of Condensed Soft Matter, Illinois Institute of Technology, Chicago, Illinois
| | - Kelli L Sylvers
- Department of Biochemistry, Carver College of Medicine, University of Iowa, Iowa City, United States
| | - Jeff Wereszczynski
- Department of Physics, Illinois Institute of Technology, Chicago, Illinois.,Center for Molecular Study of Condensed Soft Matter, Illinois Institute of Technology, Chicago, Illinois
| | - Catherine A Musselman
- Department of Biochemistry, Carver College of Medicine, University of Iowa, Iowa City, United States
| |
Collapse
|
28
|
Manzoni F, Ryde U. Assessing the stability of free-energy perturbation calculations by performing variations in the method. J Comput Aided Mol Des 2018. [PMID: 29536221 PMCID: PMC5889414 DOI: 10.1007/s10822-018-0110-5] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Abstract
We have calculated relative binding affinities for eight tetrafluorophenyl-triazole-thiogalactoside inhibitors of galectin-3 with the alchemical free-energy perturbation approach. We obtain a mean absolute deviation from experimental estimates of only 2-3 kJ/mol and a correlation coefficient (R2) of 0.5-0.8 for seven relative affinities spanning a range of up to 11 kJ/mol. We also studied the effect of using different methods to calculate the charges of the inhibitor and different sizes of the perturbed group (the atoms that are described by soft-core potentials and are allowed to have differing coordinates). However, the various approaches gave rather similar results and it is not possible to point out one approach as consistently and significantly better than the others. Instead, we suggest that such small and reasonable variations in the computational method can be used to check how stable the calculated results are and to obtain a more accurate estimate of the uncertainty than if performing only one calculation with a single computational setup.
Collapse
Affiliation(s)
- Francesco Manzoni
- Theoretical Chemistry, Department of Chemistry, Chemical Centre, Lund University, P. O. Box 124, 221 00, Lund, Sweden
| | - Ulf Ryde
- Theoretical Chemistry, Department of Chemistry, Chemical Centre, Lund University, P. O. Box 124, 221 00, Lund, Sweden.
| |
Collapse
|
29
|
Olsson MA, García-Sosa AT, Ryde U. Binding affinities of the farnesoid X receptor in the D3R Grand Challenge 2 estimated by free-energy perturbation and docking. J Comput Aided Mol Des 2018; 32:211-224. [PMID: 28879536 PMCID: PMC5767205 DOI: 10.1007/s10822-017-0056-z] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2017] [Accepted: 08/29/2017] [Indexed: 01/07/2023]
Abstract
We have studied the binding of 102 ligands to the farnesoid X receptor within the D3R Grand Challenge 2016 blind-prediction competition. First, we employed docking with five different docking software and scoring functions. The selected docked poses gave an average root-mean-squared deviation of 4.2 Å. Consensus scoring gave decent results with a Kendall's τ of 0.26 ± 0.06 and a Spearman's ρ of 0.41 ± 0.08. For a subset of 33 ligands, we calculated relative binding free energies with free-energy perturbation. Five transformations between the ligands involved a change of the net charge and we implemented and benchmarked a semi-analytic correction (Rocklin et al., J Chem Phys 139:184103, 2013) for artifacts caused by the periodic boundary conditions and Ewald summation. The results gave a mean absolute deviation of 7.5 kJ/mol compared to the experimental estimates and a correlation coefficient of R 2 = 0.1. These results were among the four best in this competition out of 22 submissions. The charge corrections were significant (7-8 kJ/mol) and always improved the results. By employing 23 intermediate states in the free-energy perturbation, there was a proper overlap between all states and the precision was 0.1-0.7 kJ/mol. However, thermodynamic cycles indicate that the sampling was insufficient in some of the perturbations.
Collapse
Affiliation(s)
- Martin A Olsson
- Department of Theoretical Chemistry, Chemical Centre, Lund University, P. O. Box 124, 221 00, Lund, Sweden
| | | | - Ulf Ryde
- Department of Theoretical Chemistry, Chemical Centre, Lund University, P. O. Box 124, 221 00, Lund, Sweden.
| |
Collapse
|
30
|
Manzoni F, Uranga J, Genheden S, Ryde U. Can System Truncation Speed up Ligand-Binding Calculations with Periodic Free-Energy Simulations? J Chem Inf Model 2017; 57:2865-2873. [PMID: 29076739 DOI: 10.1021/acs.jcim.7b00324] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We have investigated whether alchemical free-energy perturbation calculations of relative binding energies can be sped up by simulating a truncated protein. Previous studies with spherical nonperiodic systems showed that the number of simulated atoms could be reduced by a factor of 26 without affecting the calculated binding free energies by more than 0.5 kJ/mol on average ( Genheden, S.; Ryde, U. J. Chem. Theory Comput. 2012 , 8 , 1449 ), leading to a 63-fold decrease in the time consumption. However, such simulations are rather slow, owing to the need of a large cutoff radius for the nonbonded interactions. Periodic simulations with the electrostatics treated by Ewald summation are much faster. Therefore, we have investigated if a similar speed-up can be obtained also for periodic simulations. Unfortunately, our results show that it is harder to truncate periodic systems and that the truncation errors are larger for these systems. In particular, residues need to be removed from the calculations, which means that atoms have to be restrained to avoid that they move in an unrealistic manner. The results strongly depend on the strength on this restraint. For the binding of seven ligands to dihydrofolate reductase and ten inhibitors of blood-clotting factor Xa, the best results are obtained with a small restraining force constant. However, the truncation errors were still significant (e.g., 1.5-2.9 kJ/mol at a truncation radius of 10 Å). Moreover, the gain in computer time was only modest. On the other hand, if the snapshots are truncated after the MD simulations, the truncation errors are small (below 0.9 kJ/mol even for a truncation radius of 10 Å). This indicates that postprocessing with a more accurate energy function (e.g., with quantum chemistry) on truncated snapshots may be a viable approach.
Collapse
Affiliation(s)
- Francesco Manzoni
- Department of Theoretical Chemistry, Lund University, Chemical Centre , P.O. Box 124, SE-221 00 Lund, Sweden
| | - Jon Uranga
- Department of Theoretical Chemistry, Lund University, Chemical Centre , P.O. Box 124, SE-221 00 Lund, Sweden
| | - Samuel Genheden
- Department of Theoretical Chemistry, Lund University, Chemical Centre , P.O. Box 124, SE-221 00 Lund, Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre , P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
31
|
Misini Ignjatović M, Mikulskis P, Söderhjelm P, Ryde U. Can MM/GBSA calculations be sped up by system truncation? J Comput Chem 2017; 39:361-372. [DOI: 10.1002/jcc.25120] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2017] [Revised: 09/29/2017] [Accepted: 11/06/2017] [Indexed: 11/09/2022]
Affiliation(s)
- Majda Misini Ignjatović
- Department of Theoretical Chemistry; Lund University, Chemical Centre, P. O. Box 124; Lund SE-221 00 Sweden
| | - Paulius Mikulskis
- Department of Theoretical Chemistry; Lund University, Chemical Centre, P. O. Box 124; Lund SE-221 00 Sweden
| | - Pär Söderhjelm
- Department of Biophysical Chemistry; Lund University, Chemical Centre, P. O. Box 124; Lund SE-221 00 Sweden
| | - Ulf Ryde
- Department of Theoretical Chemistry; Lund University, Chemical Centre, P. O. Box 124; Lund SE-221 00 Sweden
| |
Collapse
|
32
|
Aldeghi M, Bodkin MJ, Knapp S, Biggin PC. Statistical Analysis on the Performance of Molecular Mechanics Poisson-Boltzmann Surface Area versus Absolute Binding Free Energy Calculations: Bromodomains as a Case Study. J Chem Inf Model 2017; 57:2203-2221. [PMID: 28786670 PMCID: PMC5615372 DOI: 10.1021/acs.jcim.7b00347] [Citation(s) in RCA: 90] [Impact Index Per Article: 12.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Binding free energy calculations that make use of alchemical pathways are becoming increasingly feasible thanks to advances in hardware and algorithms. Although relative binding free energy (RBFE) calculations are starting to find widespread use, absolute binding free energy (ABFE) calculations are still being explored mainly in academic settings due to the high computational requirements and still uncertain predictive value. However, in some drug design scenarios, RBFE calculations are not applicable and ABFE calculations could provide an alternative. Computationally cheaper end-point calculations in implicit solvent, such as molecular mechanics Poisson-Boltzmann surface area (MMPBSA) calculations, could too be used if one is primarily interested in a relative ranking of affinities. Here, we compare MMPBSA calculations to previously performed absolute alchemical free energy calculations in their ability to correlate with experimental binding free energies for three sets of bromodomain-inhibitor pairs. Different MMPBSA approaches have been considered, including a standard single-trajectory protocol, a protocol that includes a binding entropy estimate, and protocols that take into account the ligand hydration shell. Despite the improvements observed with the latter two MMPBSA approaches, ABFE calculations were found to be overall superior in obtaining correlation with experimental affinities for the test cases considered. A difference in weighted average Pearson ([Formula: see text]) and Spearman ([Formula: see text]) correlations of 0.25 and 0.31 was observed when using a standard single-trajectory MMPBSA setup ([Formula: see text] = 0.64 and [Formula: see text] = 0.66 for ABFE; [Formula: see text] = 0.39 and [Formula: see text] = 0.35 for MMPBSA). The best performing MMPBSA protocols returned weighted average Pearson and Spearman correlations that were about 0.1 inferior to ABFE calculations: [Formula: see text] = 0.55 and [Formula: see text] = 0.56 when including an entropy estimate, and [Formula: see text] = 0.53 and [Formula: see text] = 0.55 when including explicit water molecules. Overall, the study suggests that ABFE calculations are indeed the more accurate approach, yet there is also value in MMPBSA calculations considering the lower compute requirements, and if agreement to experimental affinities in absolute terms is not of interest. Moreover, for the specific protein-ligand systems considered in this study, we find that including an explicit ligand hydration shell or a binding entropy estimate in the MMPBSA calculations resulted in significant performance improvements at a negligible computational cost.
Collapse
Affiliation(s)
- Matteo Aldeghi
- Structural Bioinformatics and Computational Biochemistry, Department of Biochemistry, University of Oxford , South Parks Road, Oxford OX1 3QU, United Kingdom
| | - Michael J Bodkin
- Evotec (U.K.) Ltd. , 114 Innovation Drive, Milton Park, Abingdon OX14 4RZ, United Kingdom
| | - Stefan Knapp
- Structural Genomics Consortium, Nuffield Department of Clinical Medicine, University of Oxford , Old Road Campus Research Building, Roosevelt Drive, Oxford OX3 7DQ, United Kingdom.,Institute for Pharmaceutical Chemistry and Buchmann Institute for Life Sciences, Johann Wolfgang Goethe-University , D-60438 Frankfurt am Main, Germany
| | - Philip C Biggin
- Structural Bioinformatics and Computational Biochemistry, Department of Biochemistry, University of Oxford , South Parks Road, Oxford OX1 3QU, United Kingdom
| |
Collapse
|
33
|
Hu G, Ma A, Wang J. Ligand Selectivity Mechanism and Conformational Changes in Guanine Riboswitch by Molecular Dynamics Simulations and Free Energy Calculations. J Chem Inf Model 2017; 57:918-928. [PMID: 28345904 DOI: 10.1021/acs.jcim.7b00139] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Riboswitches regulate gene expression through direct and specific interactions with small metabolite molecules. Binding of a ligand to its RNA target is high selectivity and affinity and induces conformational changes of the RNA's secondary and tertiary structure. The structural difference of two purine riboswitches aptamers is caused by only one single mutation, where cytosine 74 in the guanine riboswitch is corresponding to a uracil 74 in adenine riboswitch. Here we employed molecular dynamics (MD) simulation, molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) and thermodynamic integration computational methodologies to evaluate the energetic and conformational changes of ligands binding to purine riboswitches. The snapshots used in MM-PBSA calculation were extracted from ten 50 ns MD simulation trajectories for each complex. These free energy results are in consistent with the experimental data and rationalize the selectivity of the riboswitches for different ligands. In particular, it is found that the loss in binding free energy upon mutation is mainly electrostatic in guanine (GUA) and riboswitch complex. Furthermore, new hydrogen bonds are found in mutated complexes. To reveal the conformational properties of guanine riboswitch, we performed a total of 6 μs MD simulations in both the presence and the absence of the ligand GUA. The MD simulations suggest that the conformation of guanine riboswitch depends on the distance of two groups in the binding pocket of ligand. The conformation is in a close conformation when U51-A52 is close to C74-U75.
Collapse
Affiliation(s)
- Guodong Hu
- Shandong Provincial Key Laboratory of Biophysics, College of Physics and Electronic Information, Dezhou University , Dezhou 253023, China
| | - Aijing Ma
- Shandong Provincial Key Laboratory of Biophysics, College of Physics and Electronic Information, Dezhou University , Dezhou 253023, China
| | - Jihua Wang
- Shandong Provincial Key Laboratory of Biophysics, College of Physics and Electronic Information, Dezhou University , Dezhou 253023, China
| |
Collapse
|
34
|
Ben-Shalom IY, Pfeiffer-Marek S, Baringhaus KH, Gohlke H. Efficient Approximation of Ligand Rotational and Translational Entropy Changes upon Binding for Use in MM-PBSA Calculations. J Chem Inf Model 2017; 57:170-189. [DOI: 10.1021/acs.jcim.6b00373] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Ido Y. Ben-Shalom
- Institute
for Pharmaceutical and Medicinal Chemistry, Department of Mathematics
and Natural Sciences, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| | - Stefania Pfeiffer-Marek
- LGCR/Pharmaceutical
Sciences Operations, Sanofi-Aventis Deutschland GmbH, Industriepark Höchst, 65926 Frankfurt am Main, Germany
| | - Karl-Heinz Baringhaus
- R&D Resources/Site Direction, Sanofi-Aventis Deutschland GmbH, Industriepark Höchst, 65926 Frankfurt am Main, Germany
| | - Holger Gohlke
- Institute
for Pharmaceutical and Medicinal Chemistry, Department of Mathematics
and Natural Sciences, Heinrich Heine University Düsseldorf, 40225 Düsseldorf, Germany
| |
Collapse
|
35
|
Barra PA, Ribeiro AJM, Ramos MJ, Jiménez VA, Alderete JB, Fernandes PA. Binding free energy calculations on E-selectin complexes with sLex
oligosaccharide analogs. Chem Biol Drug Des 2016; 89:114-123. [DOI: 10.1111/cbdd.12837] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2016] [Revised: 08/01/2016] [Accepted: 08/06/2016] [Indexed: 02/03/2023]
Affiliation(s)
- Pabla A. Barra
- Departamento de Química Orgánica; Facultad de Ciencias Químicas; Universidad de Concepción; Concepción Chile
| | | | - Maria J. Ramos
- Faculdade de Ciencias; Universidad do Porto; Porto Portugal
| | - Verónica A. Jiménez
- Departamento de Ciencias Químicas; Facultad de Ciencias Exactas; Universidad Andres Bello Sede Concepción; Talcahuano Chile
| | - Joel B. Alderete
- Departamento de Química Orgánica; Facultad de Ciencias Químicas; Universidad de Concepción; Concepción Chile
| | | |
Collapse
|
36
|
Misini Ignjatović M, Caldararu O, Dong G, Muñoz-Gutierrez C, Adasme-Carreño F, Ryde U. Binding-affinity predictions of HSP90 in the D3R Grand Challenge 2015 with docking, MM/GBSA, QM/MM, and free-energy simulations. J Comput Aided Mol Des 2016; 30:707-730. [PMID: 27565797 PMCID: PMC5078160 DOI: 10.1007/s10822-016-9942-z] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2016] [Accepted: 08/17/2016] [Indexed: 11/25/2022]
Abstract
We have estimated the binding affinity of three sets of ligands of the heat-shock protein 90 in the D3R grand challenge blind test competition. We have employed four different methods, based on five different crystal structures: first, we docked the ligands to the proteins with induced-fit docking with the Glide software and calculated binding affinities with three energy functions. Second, the docked structures were minimised in a continuum solvent and binding affinities were calculated with the MM/GBSA method (molecular mechanics combined with generalised Born and solvent-accessible surface area solvation). Third, the docked structures were re-optimised by combined quantum mechanics and molecular mechanics (QM/MM) calculations. Then, interaction energies were calculated with quantum mechanical calculations employing 970-1160 atoms in a continuum solvent, combined with energy corrections for dispersion, zero-point energy and entropy, ligand distortion, ligand solvation, and an increase of the basis set to quadruple-zeta quality. Fourth, relative binding affinities were estimated by free-energy simulations, using the multi-state Bennett acceptance-ratio approach. Unfortunately, the results were varying and rather poor, with only one calculation giving a correlation to the experimental affinities larger than 0.7, and with no consistent difference in the quality of the predictions from the various methods. For one set of ligands, the results could be strongly improved (after experimental data were revealed) if it was recognised that one of the ligands displaced one or two water molecules. For the other two sets, the problem is probably that the ligands bind in different modes than in the crystal structures employed or that the conformation of the ligand-binding site or the whole protein changes.
Collapse
Affiliation(s)
- Majda Misini Ignjatović
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Octav Caldararu
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Geng Dong
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden
| | - Camila Muñoz-Gutierrez
- Centro de Bioinformática y Simulación Molecular, Facultad de Ingeniería, Universidad de Talca, 2 Norte 685, Talca, Chile
| | - Francisco Adasme-Carreño
- Centro de Bioinformática y Simulación Molecular, Facultad de Ingeniería, Universidad de Talca, 2 Norte 685, Talca, Chile
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, 221 00, Lund, Sweden.
| |
Collapse
|
37
|
Bauß A, Langenmaier M, Strittmatter E, Plattner DA, Koslowski T. Thermodynamic Integration Networks and Their Application to Charge Transfer Reactions within the AauDyPI Fungal Peroxidase. J Phys Chem B 2016; 120:4937-44. [DOI: 10.1021/acs.jpcb.6b03327] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
- Anna Bauß
- Institut
für Physikalische Chemie, Universität Freiburg, Albertstraße
23a, D-79104 Freiburg
im Breisgau, Germany
| | - Michael Langenmaier
- Institut
für Physikalische Chemie, Universität Freiburg, Albertstraße
23a, D-79104 Freiburg
im Breisgau, Germany
| | - Eric Strittmatter
- Institut
für Organische Chemie, Universität Freiburg, Albertstraße
21, D-79104 Freiburg
im Breisgau, Germany
| | - Dietmar A. Plattner
- Institut
für Organische Chemie, Universität Freiburg, Albertstraße
21, D-79104 Freiburg
im Breisgau, Germany
| | - Thorsten Koslowski
- Institut
für Physikalische Chemie, Universität Freiburg, Albertstraße
23a, D-79104 Freiburg
im Breisgau, Germany
| |
Collapse
|
38
|
Su PC, Johnson ME. Evaluating thermodynamic integration performance of the new amber molecular dynamics package and assess potential halogen bonds of enoyl-ACP reductase (FabI) benzimidazole inhibitors. J Comput Chem 2016; 37:836-47. [PMID: 26666582 PMCID: PMC4769659 DOI: 10.1002/jcc.24274] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2015] [Revised: 10/23/2015] [Accepted: 11/16/2015] [Indexed: 12/17/2022]
Abstract
Thermodynamic integration (TI) can provide accurate binding free energy insights in a lead optimization program, but its high computational expense has limited its usage. In the effort of developing an efficient and accurate TI protocol for FabI inhibitors lead optimization program, we carefully compared TI with different Amber molecular dynamics (MD) engines (sander and pmemd), MD simulation lengths, the number of intermediate states and transformation steps, and the Lennard-Jones and Coulomb Softcore potentials parameters in the one-step TI, using eleven benzimidazole inhibitors in complex with Francisella tularensis enoyl acyl reductase (FtFabI). To our knowledge, this is the first study to extensively test the new AMBER MD engine, pmemd, on TI and compare the parameters of the Softcore potentials in the one-step TI in a protein-ligand binding system. The best performing model, the one-step pmemd TI, using 6 intermediate states and 1 ns MD simulations, provides better agreement with experimental results (RMSD = 0.52 kcal/mol) than the best performing implicit solvent method, QM/MM-GBSA from our previous study (RMSD = 3.00 kcal/mol), while maintaining similar efficiency. Briefly, we show the optimized TI protocol to be highly accurate and affordable for the FtFabI system. This approach can be implemented in a larger scale benzimidazole scaffold lead optimization against FtFabI. Lastly, the TI results here also provide structure-activity relationship insights, and suggest the parahalogen in benzimidazole compounds might form a weak halogen bond with FabI, which is a well-known halogen bond favoring enzyme.
Collapse
Affiliation(s)
- Pin-Chih Su
- Center for Pharmaceutical Biotechnology, College of Pharmacy, University of Illinois at Chicago, Chicago, IL, U.S.A., 60607
| | - Michael E. Johnson
- Center for Pharmaceutical Biotechnology, College of Pharmacy, University of Illinois at Chicago, Chicago, IL, U.S.A., 60607
| |
Collapse
|
39
|
Martins SA, Sousa SF, Ramos MJ, Fernandes PA. Prediction of Solvation Free Energies with Thermodynamic Integration Using the General Amber Force Field. J Chem Theory Comput 2015; 10:3570-7. [PMID: 26588320 DOI: 10.1021/ct500346y] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Computer-aided drug design (CADD) techniques can be very effective in reducing costs and speeding up drug discovery. The determination of binding and solvation free energies is pivotal for this process and is, therefore, the subject of many studies. In this work, the solvation free energy change (ΔΔGsolv) for a total of 92 transformations in small molecules was predicted using Thermodynamic Integration (TI). It was our aim to compare experimental and calculated solvation free energies for typical and prime additions considered in drug optimizations, analyzing trends, and optimizing a TI protocol. The results showed a good agreement between experimental and predicted values, with an overestimation of the predicted values for CH3, halogens, and NH2, as well as an underestimation for CONH2, but all fall within ±1 kcal/mol. NO2 addition showed a larger and systematic underestimation of the predicted ΔΔGsolv, indicating the need for special attention in these cases. For small molecules, if no experimental data is available, using TI as a theoretical strategy thus appears to be a suitable choice in CADD. It provides a good compromise between time and accuracy.
Collapse
Affiliation(s)
- Silvia A Martins
- REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto , Rua do Campo Alegre, s/n, 4169-007 Porto, Portugal
| | - Sergio F Sousa
- REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto , Rua do Campo Alegre, s/n, 4169-007 Porto, Portugal
| | - Maria João Ramos
- REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto , Rua do Campo Alegre, s/n, 4169-007 Porto, Portugal
| | - Pedro A Fernandes
- REQUIMTE, Departamento de Química e Bioquímica, Faculdade de Ciências, Universidade do Porto , Rua do Campo Alegre, s/n, 4169-007 Porto, Portugal
| |
Collapse
|
40
|
Prathipati P, Mizuguchi K. Integration of Ligand and Structure Based Approaches for CSAR-2014. J Chem Inf Model 2015; 56:974-87. [PMID: 26492437 DOI: 10.1021/acs.jcim.5b00477] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
The prediction of binding poses and affinities is an area of active interest in computer-aided drug design (CADD). Given the documented limitations with either ligand or structure based approaches, we employed an integrated approach and developed a rapid protocol for binding mode and affinity predictions. This workflow was applied to the three protein targets of Community Structure-Activity Resource-2014 (CSAR-2014) exercise: Factor Xa (FXa), Spleen Tyrosine Kinase (SYK), and tRNA (guanine-N(1))-methyltransferase (TrmD). Our docking and scoring workflow incorporates compound clustering and ligand and protein structure based pharmacophore modeling, followed by local docking, minimization, and scoring. While the former part of the protocol ensures high-quality ligand alignments and mapping, the subsequent minimization and scoring provides the predicted binding modes and affinities. We made blind predictions of docking pose for 1, 5, and 14 ligands docked into 1, 2, and 12 crystal structures of FXa, SYK, and TrmD, respectively. The resulting 174 poses were compared with cocrystallized structures (1, 5, and 14 complexes) made available at the end of CSAR. Our predicted poses were related to the experimentally determined structures with a mean root-mean-square deviation value of 3.4 Å. Further, we were able to classify high and low affinity ligands with the area under the curve values of 0.47, 0.60, and 0.69 for FXa, SYK, and TrmD, respectively, indicating the validity of our approach in at least two of the three systems. Detailed critical analysis of the results and CSAR methodology ranking procedures suggested that a straightforward application of our workflow has limitations, as some of the performance measures do not reflect the actual utility of pose and affinity predictions in the biological context of individual systems.
Collapse
Affiliation(s)
- Philip Prathipati
- National Institutes of Biomedical Innovation, Health and Nutrition , 7-6-8 Saito-Asagi, Ibaraki City, Osaka 567-0085, Japan
| | - Kenji Mizuguchi
- National Institutes of Biomedical Innovation, Health and Nutrition , 7-6-8 Saito-Asagi, Ibaraki City, Osaka 567-0085, Japan
| |
Collapse
|
41
|
Ucisik MN, Hammes-Schiffer S. Relative Binding Free Energies of Adenine and Guanine to Damaged and Undamaged DNA in Human DNA Polymerase η: Clues for Fidelity and Overall Efficiency. J Am Chem Soc 2015; 137:13240-3. [PMID: 26434497 PMCID: PMC4617179 DOI: 10.1021/jacs.5b08451] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Human DNA polymerase η (Pol η) plays an essential protective role against skin cancer caused by cyclobutane thymine-thymine dimers (TTDs), a frequent form of DNA damage arising from exposure to the sun. This enzyme rescues stalled replication forks at the TTDs by inserting bases opposite these DNA defects. Herein we calculate binding free energies for a free deoxyribose nucleotide triphosphate, dATP or dGTP, to Pol η complexed with undamaged or damaged DNA. The calculations indicate that the binding of dATP to the enzyme-DNA complex is thermodynamically favored for TTD-containing DNA over undamaged DNA, most likely because of more extensive hydrogen-bonding interactions between the TTD and the enzyme that hold the TTD more rigidly in place. The calculations also illustrate that dATP binding is thermodynamically favored over dGTP binding at both thymine positions of the TTD, most likely due to more persistent and stable hydrogen-bonding interactions between the TTD and dATP than between the TTD and dGTP. This free energy difference is slightly greater for binding at the 5' thymine position than at the 3' thymine position, presumably because of stabilization arising from the A:T base pair formed at the 3' position of the TTD in the previous step of Pol η function. All of these trends in binding free energies are consistent with experimental measurements of binding strength, fidelity, processivity, and overall efficiency. The insights gained from this analysis have implications for drug design efforts aimed at modifying the binding properties of this enzyme for improving cancer chemotherapy treatments.
Collapse
Affiliation(s)
- Melek N Ucisik
- Department of Chemistry, University of Illinois at Urbana-Champaign , 600 South Mathews Avenue, Urbana, Illinois 61801-3364, United States
| | - Sharon Hammes-Schiffer
- Department of Chemistry, University of Illinois at Urbana-Champaign , 600 South Mathews Avenue, Urbana, Illinois 61801-3364, United States
| |
Collapse
|
42
|
Su PC, Tsai CC, Mehboob S, Hevener KE, Johnson ME. Comparison of radii sets, entropy, QM methods, and sampling on MM-PBSA, MM-GBSA, and QM/MM-GBSA ligand binding energies of F. tularensis enoyl-ACP reductase (FabI). J Comput Chem 2015; 36:1859-73. [PMID: 26216222 PMCID: PMC4688044 DOI: 10.1002/jcc.24011] [Citation(s) in RCA: 71] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2015] [Revised: 04/20/2015] [Accepted: 06/22/2015] [Indexed: 01/10/2023]
Abstract
To validate a method for predicting the binding affinities of FabI inhibitors, three implicit solvent methods, MM-PBSA, MM-GBSA, and QM/MM-GBSA were carefully compared using 16 benzimidazole inhibitors in complex with Francisella tularensis FabI. The data suggests that the prediction results are sensitive to radii sets, GB methods, QM Hamiltonians, sampling protocols, and simulation length, if only one simulation trajectory is used for each ligand. In this case, QM/MM-GBSA using 6 ns MD simulation trajectories together with GB(neck2) , PM3, and the mbondi2 radii set, generate the closest agreement with experimental values (r(2) = 0.88). However, if the three implicit solvent methods are averaged from six 1 ns MD simulations for each ligand (called "multiple independent sampling"), the prediction results are relatively insensitive to all the tested parameters. Moreover, MM/GBSA together with GB(HCT) and mbondi, using 600 frames extracted evenly from six 0.25 ns MD simulations, can also provide accurate prediction to experimental values (r(2) = 0.84). Therefore, the multiple independent sampling method can be more efficient than a single, long simulation method. Since future scaffold expansions may significantly change the benzimidazole's physiochemical properties (charges, etc.) and possibly binding modes, which may affect the sensitivities of various parameters, the relatively insensitive "multiple independent sampling method" may avoid the need of an entirely new validation study. Moreover, due to large fluctuating entropy values, (QM/)MM-P(G)BSA were limited to inhibitors' relative affinity prediction, but not the absolute affinity. The developed protocol will support an ongoing benzimidazole lead optimization program.
Collapse
Affiliation(s)
- Pin-Chih Su
- Center for Pharmaceutical Biotechnology, College of Pharmacy, University of Illinois at Chicago, Chicago, IL, U.S.A., 60607
| | - Cheng-Chieh Tsai
- Center for Pharmaceutical Biotechnology, College of Pharmacy, University of Illinois at Chicago, Chicago, IL, U.S.A., 60607
| | - Shahila Mehboob
- Center for Pharmaceutical Biotechnology, College of Pharmacy, University of Illinois at Chicago, Chicago, IL, U.S.A., 60607
| | - Kirk E. Hevener
- Center for Pharmaceutical Biotechnology, College of Pharmacy, University of Illinois at Chicago, Chicago, IL, U.S.A., 60607
| | - Michael E. Johnson
- Center for Pharmaceutical Biotechnology, College of Pharmacy, University of Illinois at Chicago, Chicago, IL, U.S.A., 60607
| |
Collapse
|
43
|
Genheden S, Ryde U, Söderhjelm P. Binding affinities by alchemical perturbation using QM/MM with a large QM system and polarizable MM model. J Comput Chem 2015; 36:2114-24. [PMID: 26280564 DOI: 10.1002/jcc.24048] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2015] [Revised: 06/13/2015] [Accepted: 07/17/2015] [Indexed: 12/19/2022]
Abstract
The most general way to improve the accuracy of binding-affinity calculations for protein-ligand systems is to use quantum-mechanical (QM) methods together with rigorous alchemical-perturbation (AP) methods. We explore this approach by calculating the relative binding free energy of two synthetic disaccharides binding to galectin-3 at a reasonably high QM level (dispersion-corrected density functional theory with a triple-zeta basis set) and with a sufficiently large QM system to include all short-range interactions with the ligand (744-748 atoms). The rest of the protein is treated as a collection of atomic multipoles (up to quadrupoles) and polarizabilities. Several methods for evaluating the binding free energy from the 3600 QM calculations are investigated in terms of stability and accuracy. In particular, methods using QM calculations only at the endpoints of the transformation are compared with the recently proposed non-Boltzmann Bennett acceptance ratio (NBB) method that uses QM calculations at several stages of the transformation. Unfortunately, none of the rigorous approaches give sufficient statistical precision. However, a novel approximate method, involving the direct use of QM energies in the Bennett acceptance ratio method, gives similar results as NBB but with better precision, ∼3 kJ/mol. The statistical error can be further reduced by performing a greater number of QM calculations.
Collapse
Affiliation(s)
- Samuel Genheden
- School of Chemistry, University of Southampton, Highfield, Southampton, SO17 1BJ, United Kingdom
| | - Ulf Ryde
- Department of Theoretical Chemistry, Lund University, Chemical Centre, P. O. Box 124, Lund, SE-221 00, Sweden
| | - Pär Söderhjelm
- Department of Biophysical Chemistry, Lund University, Chemical Centre, P. O. Box 124, Lund, SE-221 00, Sweden
| |
Collapse
|
44
|
Sun H, Li Y, Tian S, Xu L, Hou T. Assessing the performance of MM/PBSA and MM/GBSA methods. 4. Accuracies of MM/PBSA and MM/GBSA methodologies evaluated by various simulation protocols using PDBbind data set. Phys Chem Chem Phys 2015; 16:16719-29. [PMID: 24999761 DOI: 10.1039/c4cp01388c] [Citation(s) in RCA: 553] [Impact Index Per Article: 61.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
By using different evaluation strategies, we systemically evaluated the performance of Molecular Mechanics/Generalized Born Surface Area (MM/GBSA) and Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) methodologies based on more than 1800 protein-ligand crystal structures in the PDBbind database. The results can be summarized as follows: (1) for the one-protein-family/one-binding-ligand case which represents the unbiased protein-ligand complex sampling, both MM/GBSA and MM/PBSA methodologies achieve approximately equal accuracies at the interior dielectric constant of 4 (with rp = 0.408 ± 0.006 of MM/GBSA and rp = 0.388 ± 0.006 of MM/PBSA based on the minimized structures); while for the total dataset (1864 crystal structures), the overall best Pearson correlation coefficient (rp = 0.579 ± 0.002) based on MM/GBSA is better than that of MM/PBSA (rp = 0.491 ± 0.003), indicating that biased sampling may significantly affect the accuracy of the predicted result (some protein families contain too many instances and can bias the overall predicted accuracy). Therefore, family based classification is needed to evaluate the two methodologies; (2) the prediction accuracies of MM/GBSA and MM/PBSA for different protein families are quite different with rp ranging from 0 to 0.9, whereas the correlation and ranking scores (an averaged rp/rs over a list of protein folds and also representing the unbiased sampling) given by MM/PBSA (rp-score = 0.506 ± 0.050 and rs-score = 0.481 ± 0.052) are comparable to those given by MM/GBSA (rp-score = 0.516 ± 0.047 and rs-score = 0.463 ± 0.047) at the fold family level; (3) for the overall prediction accuracies, molecular dynamics (MD) simulation may not be quite necessary for MM/GBSA (rp-minimized = 0.579 ± 0.002 and rp-1ns = 0.564 ± 0.002), but is needed for MM/PBSA (rp-minimized = 0.412 ± 0.003 and rp-1ns = 0.491 ± 0.003). However, for the individual systems, whether to use MD simulation is depended. (4) both MM/GBSA and MM/PBSA may be unable to give successful predictions for the ligands with high formal charges, with the Pearson correlation coefficient ranging from 0.621 ± 0.003 (neutral ligands) to 0.125 ± 0.142 (ligands with a formal charge of 5). Therefore, it can be summarized that, although MM/GBSA and MM/PBSA perform similarly in the unbiased dataset, for the currently available crystal structures in the PDBbind database, compared with MM/GBSA, which may be used in multi-target comparisons, MM/PBSA is more sensitive to the investigated systems, and may be more suitable for individual-target-level binding free energy ranking. This study may provide useful guidance for the post-processing of docking based studies.
Collapse
Affiliation(s)
- Huiyong Sun
- Institute of Functional Nano & Soft Materials (FUNSOM) and Collaborative Innovation Center of Suzhou Nano Science and Technology, Soochow University, Suzhou, Jiangsu 215123, China.
| | | | | | | | | |
Collapse
|
45
|
Abstract
INTRODUCTION The molecular mechanics energies combined with the Poisson-Boltzmann or generalized Born and surface area continuum solvation (MM/PBSA and MM/GBSA) methods are popular approaches to estimate the free energy of the binding of small ligands to biological macromolecules. They are typically based on molecular dynamics simulations of the receptor-ligand complex and are therefore intermediate in both accuracy and computational effort between empirical scoring and strict alchemical perturbation methods. They have been applied to a large number of systems with varying success. AREAS COVERED The authors review the use of MM/PBSA and MM/GBSA methods to calculate ligand-binding affinities, with an emphasis on calibration, testing and validation, as well as attempts to improve the methods, rather than on specific applications. EXPERT OPINION MM/PBSA and MM/GBSA are attractive approaches owing to their modular nature and that they do not require calculations on a training set. They have been used successfully to reproduce and rationalize experimental findings and to improve the results of virtual screening and docking. However, they contain several crude and questionable approximations, for example, the lack of conformational entropy and information about the number and free energy of water molecules in the binding site. Moreover, there are many variants of the method and their performance varies strongly with the tested system. Likewise, most attempts to ameliorate the methods with more accurate approaches, for example, quantum-mechanical calculations, polarizable force fields or improved solvation have deteriorated the results.
Collapse
Affiliation(s)
- Samuel Genheden
- University of Southampton, School of Chemistry, Highfield, SO17 1BJ, Southampton, UK
| | - Ulf Ryde
- Lund University, Chemical Centre, Department of Theoretical Chemistry, P. O. Box 124, SE-221 00 Lund, Sweden+46 46 2224502; +46 46 2228648;
| |
Collapse
|
46
|
Chen J, Liang Z, Wang W, Yi C, Zhang S, Zhang Q. Revealing origin of decrease in potency of darunavir and amprenavir against HIV-2 relative to HIV-1 protease by molecular dynamics simulations. Sci Rep 2014; 4:6872. [PMID: 25362963 PMCID: PMC4217091 DOI: 10.1038/srep06872] [Citation(s) in RCA: 63] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2014] [Accepted: 10/13/2014] [Indexed: 12/30/2022] Open
Abstract
Clinical inhibitors Darunavir (DRV) and Amprenavir (APV) are less effective on HIV-2 protease (PR2) than on HIV-1 protease (PR1). To identify molecular basis associated with the lower inhibition, molecular dynamics (MD) simulations and molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) calculations were performed to investigate the effectiveness of the PR1 inhibitors DRV and APV against PR1/PR2. The rank of predicted binding free energies agrees with the experimental determined one. Moreover, our results show that two inhibitors bind less strongly to PR2 than to PR1, again in agreement with the experimental findings. The decrease in binding free energies for PR2 relative to PR1 is found to arise from the reduction of the van der Waals interactions induced by the structural adjustment of the triple mutant V32I, I47V and V82I. This result is further supported by the difference between the van der Waals interactions of inhibitors with each residue in PR2 and in PR1. The results from the principle component analysis suggest that inhibitor binding tends to make the flaps of PR2 close and the one of PR1 open. We expect that this study can theoretically provide significant guidance and dynamics information for the design of potent dual inhibitors targeting PR1/PR2.
Collapse
Affiliation(s)
- Jianzhong Chen
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Zhiqiang Liang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Wei Wang
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Changhong Yi
- School of Science, Shandong Jiaotong University, Jinan 250357, China
| | - Shaolong Zhang
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| | - Qinggang Zhang
- College of Physics and Electronics, Shandong Normal University, Jinan, 250014, China
| |
Collapse
|
47
|
Mikulskis P, Genheden S, Ryde U. A large-scale test of free-energy simulation estimates of protein-ligand binding affinities. J Chem Inf Model 2014; 54:2794-806. [PMID: 25264937 DOI: 10.1021/ci5004027] [Citation(s) in RCA: 52] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
We have performed a large-scale test of alchemical perturbation calculations with the Bennett acceptance-ratio (BAR) approach to estimate relative affinities for the binding of 107 ligands to 10 different proteins. Employing 20-Å truncated spherical systems and only one intermediate state in the perturbations, we obtain an error of less than 4 kJ/mol for 54% of the studied relative affinities and a precision of 0.5 kJ/mol on average. However, only four of the proteins gave acceptable errors, correlations, and rankings. The results could be improved by using nine intermediate states in the simulations or including the entire protein in the simulations using periodic boundary conditions. However, 27 of the calculated affinities still gave errors of more than 4 kJ/mol, and for three of the proteins the results were not satisfactory. This shows that the performance of BAR calculations depends on the target protein and that several transformations gave poor results owing to limitations in the molecular-mechanics force field or the restricted sampling possible within a reasonable simulation time. Still, the BAR results are better than docking calculations for most of the proteins.
Collapse
Affiliation(s)
- Paulius Mikulskis
- Department of Theoretical Chemistry, Lund University, Chemical Centre , P.O. Box 124, SE-221 00 Lund, Sweden
| | | | | |
Collapse
|
48
|
Gadhe CG, Balupuri A, Kothandan G, Cho SJ. Theoretical Characterization of Binding Mode of Organosilicon Inhibitor with p38: Docking, MD Simulation and MM/GBSA Free Energy Approach. B KOREAN CHEM SOC 2014. [DOI: 10.5012/bkcs.2014.35.8.2494] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
49
|
Homeyer N, Stoll F, Hillisch A, Gohlke H. Binding Free Energy Calculations for Lead Optimization: Assessment of Their Accuracy in an Industrial Drug Design Context. J Chem Theory Comput 2014; 10:3331-44. [PMID: 26588302 DOI: 10.1021/ct5000296] [Citation(s) in RCA: 119] [Impact Index Per Article: 11.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
Correctly ranking compounds according to their computed relative binding affinities will be of great value for decision making in the lead optimization phase of industrial drug discovery. However, the performance of existing computationally demanding binding free energy calculation methods in this context is largely unknown. We analyzed the performance of the molecular mechanics continuum solvent, the linear interaction energy (LIE), and the thermodynamic integration (TI) approach for three sets of compounds from industrial lead optimization projects. The data sets pose challenges typical for this early stage of drug discovery. None of the methods was sufficiently predictive when applied out of the box without considering these challenges. Detailed investigations of failures revealed critical points that are essential for good binding free energy predictions. When data set-specific features were considered accordingly, predictions valuable for lead optimization could be obtained for all approaches but LIE. Our findings lead to clear recommendations for when to use which of the above approaches. Our findings also stress the important role of expert knowledge in this process, not least for estimating the accuracy of prediction results by TI, using indicators such as the size and chemical structure of exchanged groups and the statistical error in the predictions. Such knowledge will be invaluable when it comes to the question which of the TI results can be trusted for decision making.
Collapse
Affiliation(s)
- Nadine Homeyer
- Institute for Pharmaceutical and Medicinal Chemistry, Department of Mathematics and Natural Sciences, Heinrich-Heine-University Düsseldorf , Universitätsstr.1, 40225 Düsseldorf, Germany
| | - Friederike Stoll
- Global Drug Discovery, Medicinal Chemistry, Bayer Pharma AG , Aprather Weg 18A, 42113 Wuppertal, Germany
| | - Alexander Hillisch
- Global Drug Discovery, Medicinal Chemistry, Bayer Pharma AG , Aprather Weg 18A, 42113 Wuppertal, Germany
| | - Holger Gohlke
- Institute for Pharmaceutical and Medicinal Chemistry, Department of Mathematics and Natural Sciences, Heinrich-Heine-University Düsseldorf , Universitätsstr.1, 40225 Düsseldorf, Germany
| |
Collapse
|
50
|
Lee HC, Hsu WC, Liu AL, Hsu CJ, Sun YC. Using thermodynamic integration MD simulation to compute relative protein–ligand binding free energy of a GSK3β kinase inhibitor and its analogs. J Mol Graph Model 2014; 51:37-49. [DOI: 10.1016/j.jmgm.2014.04.010] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2013] [Revised: 04/20/2014] [Accepted: 04/22/2014] [Indexed: 01/15/2023]
|