1
|
Wang J, Miao Y. Ligand Gaussian Accelerated Molecular Dynamics 3 (LiGaMD3): Improved Calculations of Binding Thermodynamics and Kinetics of Both Small Molecules and Flexible Peptides. J Chem Theory Comput 2024; 20:5829-5841. [PMID: 39002136 DOI: 10.1021/acs.jctc.4c00502] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/15/2024]
Abstract
Binding thermodynamics and kinetics play critical roles in drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics of small molecules and flexible peptides using conventional molecular dynamics (cMD), due to limited simulation time scales. Based on our previously developed ligand Gaussian accelerated molecular dynamics (LiGaMD) method, we present a new approach, termed "LiGaMD3″, in which we introduce triple boosts into three individual energy terms that play important roles in small-molecule/peptide dissociation, rebinding, and system conformational changes to improve the sampling efficiency of small-molecule/peptide interactions with target proteins. To validate the performance of LiGaMD3, MDM2 bound by a small molecule (Nutlin 3) and two highly flexible peptides (PMI and P53) were chosen as the model systems. LiGaMD3 could efficiently capture repetitive small-molecule/peptide dissociation and binding events within 2 μs simulations. The predicted binding kinetic constant rates and free energies from LiGaMD3 were in agreement with the available experimental values and previous simulation results. Therefore, LiGaMD3 provides a more general and efficient approach to capture dissociation and binding of both small-molecule ligands and flexible peptides, allowing for accurate prediction of their binding thermodynamics and kinetics.
Collapse
Affiliation(s)
- Jinan Wang
- Computational Medicine Program and Department of Pharmacology, University of North Carolina-Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Yinglong Miao
- Computational Medicine Program and Department of Pharmacology, University of North Carolina-Chapel Hill, Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
2
|
Wang J, Miao Y. Ligand Gaussian accelerated Molecular Dynamics 3 (LiGaMD3): Improved Calculations of Binding Thermodynamics and Kinetics of Both Small Molecules and Flexible Peptides. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.06.592668. [PMID: 38766067 PMCID: PMC11100592 DOI: 10.1101/2024.05.06.592668] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/22/2024]
Abstract
Binding thermodynamics and kinetics play critical roles in drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics of small molecules and flexible peptides using conventional Molecular Dynamics (cMD), due to limited simulation timescales. Based on our previously developed Ligand Gaussian accelerated Molecular Dynamics (LiGaMD) method, we present a new approach, termed "LiGaMD3", in which we introduce triple boosts into three individual energy terms that play important roles in small-molecule/peptide dissociation, rebinding and system conformational changes to improve the sampling efficiency of small-molecule/peptide interactions with target proteins. To validate the performance of LiGaMD3, MDM2 bound by a small molecule (Nutlin 3) and two highly flexible peptides (PMI and P53) were chosen as model systems. LiGaMD3 could efficiently capture repetitive small-molecule/peptide dissociation and binding events within 2 microsecond simulations. The predicted binding kinetic constant rates and free energies from LiGaMD3 agreed with available experimental values and previous simulation results. Therefore, LiGaMD3 provides a more general and efficient approach to capture dissociation and binding of both small-molecule ligand and flexible peptides, allowing for accurate prediction of their binding thermodynamics and kinetics.
Collapse
Affiliation(s)
- Jinan Wang
- Computational Medicine Program and Department of Pharmacology, University of North Carolina – Chapel Hill, Chapel Hill, North Carolina, USA 27599
| | - Yinglong Miao
- Computational Medicine Program and Department of Pharmacology, University of North Carolina – Chapel Hill, Chapel Hill, North Carolina, USA 27599
| |
Collapse
|
3
|
Ligand Gaussian Accelerated Molecular Dynamics 2 (LiGaMD2): Improved Calculations of Ligand Binding Thermodynamics and Kinetics with Closed Protein Pocket. J Chem Theory Comput 2023; 19:733-745. [PMID: 36706316 DOI: 10.1021/acs.jctc.2c01194] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Ligand binding thermodynamics and kinetics are critical parameters for drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics from molecular simulations due to limited simulation timescales. Protein dynamics, especially in the ligand binding pocket, often plays an important role in ligand binding. Based on our previously developed Ligand Gaussian accelerated molecular dynamics (LiGaMD), here we present LiGaMD2 in which a selective boost potential was applied to both the ligand and protein residues in the binding pocket to improve sampling of ligand binding and dissociation. To validate the performance of LiGaMD2, the T4 lysozyme (T4L) mutants with open and closed pockets bound by different ligands were chosen as model systems. LiGaMD2 could efficiently capture repetitive ligand dissociation and binding within microsecond simulations of all T4L systems. The obtained ligand binding kinetic rates and free energies agreed well with available experimental values and previous modeling results. Therefore, LiGaMD2 provides an improved approach to sample opening of closed protein pockets for ligand dissociation and binding, thereby allowing for efficient calculations of ligand binding thermodynamics and kinetics.
Collapse
|
4
|
Non-ionic surfactants in drug delivery vehicles: Physicochemical insights with systems of drugs, Igepal CA-630, bovine serum albumin and hen egg-white lysozyme. Colloids Surf A Physicochem Eng Asp 2023. [DOI: 10.1016/j.colsurfa.2022.130433] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
5
|
Minetti CA, Remeta DP. Forces Driving a Magic Bullet to Its Target: Revisiting the Role of Thermodynamics in Drug Design, Development, and Optimization. Life (Basel) 2022; 12:1438. [PMID: 36143474 PMCID: PMC9504344 DOI: 10.3390/life12091438] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2022] [Revised: 09/02/2022] [Accepted: 09/03/2022] [Indexed: 11/27/2022] Open
Abstract
Drug discovery strategies have advanced significantly towards prioritizing target selectivity to achieve the longstanding goal of identifying "magic bullets" amongst thousands of chemical molecules screened for therapeutic efficacy. A myriad of emerging and existing health threats, including the SARS-CoV-2 pandemic, alarming increase in bacterial resistance, and potentially fatal chronic ailments, such as cancer, cardiovascular disease, and neurodegeneration, have incentivized the discovery of novel therapeutics in treatment regimens. The design, development, and optimization of lead compounds represent an arduous and time-consuming process that necessitates the assessment of specific criteria and metrics derived via multidisciplinary approaches incorporating functional, structural, and energetic properties. The present review focuses on specific methodologies and technologies aimed at advancing drug development with particular emphasis on the role of thermodynamics in elucidating the underlying forces governing ligand-target interaction selectivity and specificity. In the pursuit of novel therapeutics, isothermal titration calorimetry (ITC) has been utilized extensively over the past two decades to bolster drug discovery efforts, yielding information-rich thermodynamic binding signatures. A wealth of studies recognizes the need for mining thermodynamic databases to critically examine and evaluate prospective drug candidates on the basis of available metrics. The ultimate power and utility of thermodynamics within drug discovery strategies reside in the characterization and comparison of intrinsic binding signatures that facilitate the elucidation of structural-energetic correlations which assist in lead compound identification and optimization to improve overall therapeutic efficacy.
Collapse
Affiliation(s)
- Conceição A. Minetti
- Department of Chemistry and Chemical Biology, Rutgers—The State University of New Jersey, Piscataway, NJ 08854, USA
| | - David P. Remeta
- Department of Chemistry and Chemical Biology, Rutgers—The State University of New Jersey, Piscataway, NJ 08854, USA
| |
Collapse
|
6
|
Schuetz DA, Richter L, Martini R, Ecker GF. A structure-kinetic relationship study using matched molecular pair analysis. RSC Med Chem 2020; 11:1285-1294. [PMID: 34085042 PMCID: PMC8126976 DOI: 10.1039/d0md00178c] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
The lifetime of a binary drug–target complex is increasingly acknowledged as an important parameter for drug efficacy and safety. With a better understanding of binding kinetics and better knowledge about kinetic parameter optimization, intentionally induced prolongation of the drug–target residence time through structural changes of the ligand could become feasible. In this study we assembled datasets from 21 publications and the K4DD (Kinetic for Drug Discovery) database to conduct large scale data analysis. This resulted in 3812 small molecules annotated to 78 different targets from five protein classes (GPCRs: 273, kinases: 3238, other enzymes: 240, HSPs: 160, ion channels: 45). Performing matched molecular pair (MMP) analysis to further investigate the structure–kinetic relationship (SKR) in this data collection allowed us to identify a fundamental contribution of a ligand's polarity to its association rate, and in selected cases, also to its dissociation rate. However, we furthermore observed that the destabilization of the transition state introduced by increased polarity is often accompanied by simultaneous destabilization of the ground state resulting in an unaffected or even worsened residence time. Supported by a set of case studies, we provide concepts on how to alter ligands in ways to trigger on-rates, off-rates, or both. A large-scale study employing matched molecular pair (MMP) analysis to uncover the contribution of a compound's polarity to its association and dissociation rates.![]()
Collapse
Affiliation(s)
- Doris A Schuetz
- Department of Pharmaceutical Chemistry, University of Vienna UZA 2, Althanstrasse 14 1090 Vienna Austria
| | - Lars Richter
- Department of Pharmaceutical Chemistry, University of Vienna UZA 2, Althanstrasse 14 1090 Vienna Austria
| | - Riccardo Martini
- Department of Pharmaceutical Chemistry, University of Vienna UZA 2, Althanstrasse 14 1090 Vienna Austria
| | - Gerhard F Ecker
- Department of Pharmaceutical Chemistry, University of Vienna UZA 2, Althanstrasse 14 1090 Vienna Austria
| |
Collapse
|
7
|
Du Y, Wang R. Revealing the Unbinding Kinetics and Mechanism of Type I and Type II Protein Kinase Inhibitors by Local-Scaled Molecular Dynamics Simulations. J Chem Theory Comput 2020; 16:6620-6632. [PMID: 32841004 DOI: 10.1021/acs.jctc.0c00342] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Protein kinase inhibitors disrupt phosphorylation of the target kinases, which are an important class of drug for treating cancer and other diseases. Conventional structure-based design methods (such as molecular docking) focus on the static binding mode of the kinase inhibitor with its target. However, dissociation kinetic properties of a drug molecule are found to correlate with its residence time in vivo and thus have drawn the attention of drug designers in recent years. In this study, we have applied the local-scaled molecular dynamics (MD) simulation enabled in GROMACS software to explore the unbinding mechanism of a total of 41 type I and type II kinase inhibitors. Our simulation considered multiple starting configurations as well as possible protonation states of kinase inhibitors. Based on our local-scaled MD results, we discovered that the integrals of the favorable binding energy during dissociation correlated well (R2 = 0.64) with the experimental dissociation rate constants of those kinase inhibitors on the entire data set. Given its accuracy and technical advantage, this method may serve as a practical option for estimating this important property in reality. Our simulation also provided a reasonable explanation of the dynamic properties of kinase and its inhibitor as well as the role of relevant water molecules in dissociation.
Collapse
Affiliation(s)
- Yu Du
- State Key Laboratory of Bioorganic and Natural Products Chemistry, Center for Excellence in Molecular Synthesis, Shanghai Institute of Organic Chemistry, Chinese Academy of Sciences, Shanghai 200032, People's Republic of China.,University of Chinese Academy of Sciences, Beijing 100049, People's Republic of China
| | - Renxiao Wang
- State Key Laboratory of Bioorganic and Natural Products Chemistry, Center for Excellence in Molecular Synthesis, Shanghai Institute of Organic Chemistry, Chinese Academy of Sciences, Shanghai 200032, People's Republic of China.,University of Chinese Academy of Sciences, Beijing 100049, People's Republic of China.,Department of Medicinal Chemistry, School of Pharmacy, Fudan University, Shanghai 201203, People's Republic of China
| |
Collapse
|
8
|
Multisecond ligand dissociation dynamics from atomistic simulations. Nat Commun 2020; 11:2918. [PMID: 32522984 PMCID: PMC7286908 DOI: 10.1038/s41467-020-16655-1] [Citation(s) in RCA: 51] [Impact Index Per Article: 12.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2020] [Accepted: 05/12/2020] [Indexed: 12/22/2022] Open
Abstract
Coarse-graining of fully atomistic molecular dynamics simulations is a long-standing goal in order to allow the description of processes occurring on biologically relevant timescales. For example, the prediction of pathways, rates and rate-limiting steps in protein-ligand unbinding is crucial for modern drug discovery. To achieve the enhanced sampling, we perform dissipation-corrected targeted molecular dynamics simulations, which yield free energy and friction profiles of molecular processes under consideration. Subsequently, we use these fields to perform temperature-boosted Langevin simulations which account for the desired kinetics occurring on multisecond timescales and beyond. Adopting the dissociation of solvated sodium chloride, trypsin-benzamidine and Hsp90-inhibitor protein-ligand complexes as test problems, we reproduce rates from molecular dynamics simulation and experiments within a factor of 2–20, and dissociation constants within a factor of 1–4. Analysis of friction profiles reveals that binding and unbinding dynamics are mediated by changes of the surrounding hydration shells in all investigated systems. Protein-ligand unbinding processes are out of reach for atomistic simulations due to time-scale involved. Here the authors demonstrate an approach relying on dissipation-corrected targeted molecular dynamics that enables to provide binding and unbinding rates with a speed-up of several orders of magnitude.
Collapse
|
9
|
Koehler M, Delguste M, Sieben C, Gillet L, Alsteens D. Initial Step of Virus Entry: Virion Binding to Cell-Surface Glycans. Annu Rev Virol 2020; 7:143-165. [PMID: 32396772 DOI: 10.1146/annurev-virology-122019-070025] [Citation(s) in RCA: 78] [Impact Index Per Article: 19.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Virus infection is an intricate process that requires the concerted action of both viral and host cell components. Entry of viruses into cells is initiated by interactions between viral proteins and cell-surface receptors. Various cell-surface glycans function as initial, usually low-affinity attachment factors, providing a first anchor of the virus to the cell surface, and further facilitate high-affinity binding to virus-specific cell-surface receptors, while other glycans function as specific entry receptors themselves. It is now possible to rapidly identify specific glycan receptors using different techniques, define atomic-level structures of virus-glycan complexes, and study these interactions at the single-virion level. This review provides a detailed overview of the role of glycans in viral infection and highlights experimental approaches to study virus-glycan binding along with specific examples. In particular, we highlight the development of the atomic force microscope to investigate interactions with glycans at the single-virion level directly on living mammalian cells, which offers new perspectives to better understand virus-glycan interactions in physiologically relevant conditions.
Collapse
Affiliation(s)
- Melanie Koehler
- Louvain Institute of Biomolecular Science and Technology, Université Catholique de Louvain, 1348 Louvain-la-Neuve, Belgium;
| | - Martin Delguste
- Louvain Institute of Biomolecular Science and Technology, Université Catholique de Louvain, 1348 Louvain-la-Neuve, Belgium;
| | - Christian Sieben
- Institute of Physics, École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| | - Laurent Gillet
- Immunology-Vaccinology Laboratory, Department of Infectious and Parasitic Diseases, Fundamental and Applied Research for Animals and Health center (FARAH), University of Liège, 4000 Liège, Belgium
| | - David Alsteens
- Louvain Institute of Biomolecular Science and Technology, Université Catholique de Louvain, 1348 Louvain-la-Neuve, Belgium; .,Walloon Excellence in Life sciences and Biotechnology (WELBIO), 1300 Wavre, Belgium
| |
Collapse
|
10
|
Glöckner S, Ngo K, Sager CP, Hüfner-Wulsdorf T, Heine A, Klebe G. Conformational Changes in Alkyl Chains Determine the Thermodynamic and Kinetic Binding Profiles of Carbonic Anhydrase Inhibitors. ACS Chem Biol 2020; 15:675-685. [PMID: 32027480 DOI: 10.1021/acschembio.9b00895] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Thermodynamics and kinetics of protein-ligand binding are both important aspects for the design of novel drug molecules. Presently, thermodynamic data are collected with isothermal titration calorimetry, while kinetic data are mostly derived from surface plasmon resonance. The new method of kinITC provides both thermodynamic and kinetic data from calorimetric titration measurements. The present study demonstrates the convenient collection of calorimetric data suitable for both thermodynamic and kinetic analysis for two series of congeneric ligands of human carbonic anhydrase II and correlates these findings with structural data obtained by macromolecular crystallography to shed light on the importance of shape complementarity for thermodynamics and kinetics governing a protein-ligand binding event. The study shows how minute chemical alterations change preferred ligand conformation and can be used to manipulate thermodynamic and kinetic signatures of binding. They give rise to the observation that analogous n-alkyl and n-alkyloxy derivatives of identical chain length swap their binding kinetic properties at unchanged binding affinity.
Collapse
Affiliation(s)
- Steffen Glöckner
- Institut für Pharmazeutische Chemie, Philipps-Universität Marburg, Marbacher Weg 6, 35032 Marburg, Germany
| | - Khang Ngo
- Institut für Pharmazeutische Chemie, Philipps-Universität Marburg, Marbacher Weg 6, 35032 Marburg, Germany
| | - Christoph P Sager
- Institut für Pharmazeutische Chemie, Philipps-Universität Marburg, Marbacher Weg 6, 35032 Marburg, Germany
| | - Tobias Hüfner-Wulsdorf
- Institut für Pharmazeutische Chemie, Philipps-Universität Marburg, Marbacher Weg 6, 35032 Marburg, Germany
| | - Andreas Heine
- Institut für Pharmazeutische Chemie, Philipps-Universität Marburg, Marbacher Weg 6, 35032 Marburg, Germany
| | - Gerhard Klebe
- Institut für Pharmazeutische Chemie, Philipps-Universität Marburg, Marbacher Weg 6, 35032 Marburg, Germany
| |
Collapse
|
11
|
Hüfner-Wulsdorf T, Klebe G. Role of Water Molecules in Protein–Ligand Dissociation and Selectivity Discrimination: Analysis of the Mechanisms and Kinetics of Biomolecular Solvation Using Molecular Dynamics. J Chem Inf Model 2020; 60:1818-1832. [DOI: 10.1021/acs.jcim.0c00156] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- Tobias Hüfner-Wulsdorf
- Institut für Pharmazeutische Chemie, Philipps Universität Marburg, Marbacher Weg 6, 35037 Marburg, Germany
| | - Gerhard Klebe
- Institut für Pharmazeutische Chemie, Philipps Universität Marburg, Marbacher Weg 6, 35037 Marburg, Germany
| |
Collapse
|
12
|
Wolf S, Amaral M, Lowinski M, Vallée F, Musil D, Güldenhaupt J, Dreyer MK, Bomke J, Frech M, Schlitter J, Gerwert K. Estimation of Protein-Ligand Unbinding Kinetics Using Non-Equilibrium Targeted Molecular Dynamics Simulations. J Chem Inf Model 2019; 59:5135-5147. [PMID: 31697501 DOI: 10.1021/acs.jcim.9b00592] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
We here report on nonequilibrium targeted molecular dynamics simulations as a tool for the estimation of protein-ligand unbinding kinetics. Correlating simulations with experimental data from SPR kinetics measurements and X-ray crystallography on two small molecule compound libraries bound to the N-terminal domain of the chaperone Hsp90, we show that the mean nonequilibrium work computed in an ensemble of trajectories of enforced ligand unbinding is a promising predictor for ligand unbinding rates. We furthermore investigate the molecular basis determining unbinding rates within the compound libraries. We propose ligand conformational changes and protein-ligand nonbonded interactions to impact on unbinding rates. Ligands may remain longer at the protein if they exhibit strong electrostatic and/or van der Waals interactions with the target. In the case of ligands with a rigid chemical scaffold that exhibit longer residence times, transient electrostatic interactions with the protein appear to facilitate unbinding. Our results imply that understanding the unbinding pathway and the protein-ligand interactions along this path is crucial for the prediction of small molecule ligands with defined unbinding kinetics.
Collapse
Affiliation(s)
- Steffen Wolf
- Department of Biophysics , Ruhr-University Bochum , 44780 Bochum , Germany.,Institute of Physics , Albert-Ludwigs-University Freiburg , 79104 Freiburg , Germany
| | - Marta Amaral
- Instituto de Biologia Experimental e Tecnológica , 2780-157 Oeiras , Portugal.,Molecular Interactions and Biophysics , Merck KGaA , 64293 Darmstadt , Germany.,Sanofi-Aventis Deutschland GmbH , Biologics Research/Protein Therapeutics , 65926 Frankfurt am Main , Germany
| | - Maryse Lowinski
- Sanofi IDD-BioStructure and Biophysics , 94400 Vitry-sur-Seine , France
| | - Francois Vallée
- Sanofi IDD-BioStructure and Biophysics , 94400 Vitry-sur-Seine , France
| | - Djordje Musil
- Molecular Interactions and Biophysics , Merck KGaA , 64293 Darmstadt , Germany
| | - Jörn Güldenhaupt
- Department of Biophysics , Ruhr-University Bochum , 44780 Bochum , Germany
| | - Matthias K Dreyer
- Sanofi-Aventis Deutschland GmbH , R&D Integrated Drug Discovery , 65926 Frankfurt am Main , Germany
| | - Jörg Bomke
- Molecular Pharmacology , Merck KGaA , 64293 Darmstadt , Germany
| | - Matthias Frech
- Molecular Interactions and Biophysics , Merck KGaA , 64293 Darmstadt , Germany
| | - Jürgen Schlitter
- Department of Biophysics , Ruhr-University Bochum , 44780 Bochum , Germany
| | - Klaus Gerwert
- Department of Biophysics , Ruhr-University Bochum , 44780 Bochum , Germany
| |
Collapse
|
13
|
Magarkar A, Schnapp G, Apel AK, Seeliger D, Tautermann CS. Enhancing Drug Residence Time by Shielding of Intra-Protein Hydrogen Bonds: A Case Study on CCR2 Antagonists. ACS Med Chem Lett 2019; 10:324-328. [PMID: 30891134 DOI: 10.1021/acsmedchemlett.8b00590] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2018] [Accepted: 02/07/2019] [Indexed: 12/15/2022] Open
Abstract
The target residence time (RT) for a given ligand is one of the important parameters that have to be optimized during drug design. It is well established that shielding the receptor-ligand hydrogen bond (H-bond) interactions from water has been one of the factors in increasing ligand RT. Building on this foundation, here we report that shielding an intra-protein H-bond, which confers rigidity to the binding pocket and which is not directly involved in drug-receptor interactions, can strongly influence RT for CCR2 antagonists. Based on our recently solved CCR2 structure with MK-0812 and molecular dynamics (MD) simulations, we show that the RT for this and structurally related ligands is directly dependent on the shielding of the Tyr120-Glu291 H-bond from the water. If solvated this H-bond is often broken, making the binding pocket flexible and leading to shorter RT.
Collapse
Affiliation(s)
- Aniket Magarkar
- Medicinal Chemistry, Boehringer Ingelheim Pharma GmbH & Co. KG, Birkendorfer Strasse 65, D-88397 Biberach a.d. Riss, Germany
| | - Gisela Schnapp
- Medicinal Chemistry, Boehringer Ingelheim Pharma GmbH & Co. KG, Birkendorfer Strasse 65, D-88397 Biberach a.d. Riss, Germany
| | - Anna-Katharina Apel
- Medicinal Chemistry, Boehringer Ingelheim Pharma GmbH & Co. KG, Birkendorfer Strasse 65, D-88397 Biberach a.d. Riss, Germany
| | - Daniel Seeliger
- Medicinal Chemistry, Boehringer Ingelheim Pharma GmbH & Co. KG, Birkendorfer Strasse 65, D-88397 Biberach a.d. Riss, Germany
| | - Christofer S. Tautermann
- Medicinal Chemistry, Boehringer Ingelheim Pharma GmbH & Co. KG, Birkendorfer Strasse 65, D-88397 Biberach a.d. Riss, Germany
| |
Collapse
|
14
|
Schuetz DA, Bernetti M, Bertazzo M, Musil D, Eggenweiler HM, Recanatini M, Masetti M, Ecker GF, Cavalli A. Predicting Residence Time and Drug Unbinding Pathway through Scaled Molecular Dynamics. J Chem Inf Model 2018; 59:535-549. [DOI: 10.1021/acs.jcim.8b00614] [Citation(s) in RCA: 44] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Doris A. Schuetz
- Department of Pharmaceutical Chemistry, University of Vienna, UZA 2, Althanstrasse 14, 1090 Vienna, Austria
| | - Mattia Bernetti
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum—Università di Bologna, via Belmeloro 6, I-40126 Bologna, Italy
| | - Martina Bertazzo
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum—Università di Bologna, via Belmeloro 6, I-40126 Bologna, Italy
- Computational Sciences, Istituto Italiano di Tecnologia, via Morego 30, 16163 Genova, Italy
| | - Djordje Musil
- Discovery Technologies, Merck KGaA, Frankfurter Straße 250, 64293 Darmstadt, Germany
| | | | - Maurizio Recanatini
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum—Università di Bologna, via Belmeloro 6, I-40126 Bologna, Italy
| | - Matteo Masetti
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum—Università di Bologna, via Belmeloro 6, I-40126 Bologna, Italy
| | - Gerhard F. Ecker
- Department of Pharmaceutical Chemistry, University of Vienna, UZA 2, Althanstrasse 14, 1090 Vienna, Austria
| | - Andrea Cavalli
- Department of Pharmacy and Biotechnology, Alma Mater Studiorum—Università di Bologna, via Belmeloro 6, I-40126 Bologna, Italy
- Computational Sciences, Istituto Italiano di Tecnologia, via Morego 30, 16163 Genova, Italy
| |
Collapse
|
15
|
Su H, Xu Y. Application of ITC-Based Characterization of Thermodynamic and Kinetic Association of Ligands With Proteins in Drug Design. Front Pharmacol 2018; 9:1133. [PMID: 30364164 PMCID: PMC6193069 DOI: 10.3389/fphar.2018.01133] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2018] [Accepted: 09/18/2018] [Indexed: 01/20/2023] Open
Abstract
A comprehensive characterization of the thermodynamic and kinetic profiling of ligands binding to a given target protein is crucial for the hit selection as well as the hit-to-lead-to-drug evolution. Isothermal titration calorimetry (ITC), widely known as an invaluable tool to measure the thermodynamic data, has recently found its way to determine the binding kinetics too. The extensive application of ITC in measurement of both thermodynamic and kinetic data manifests unique roles of ITC in drug discovery and development. This mini-review concentrates on elaborating how to gain the thermodynamic and kinetic data using ITC, highlighting the importance of these data in lead discovery and optimization, and intends to provide an overview of the technical and conceptual advances that offer unprecedented access to protein–ligand recognition by ITC measurement.
Collapse
Affiliation(s)
- Haixia Su
- CAS Key Laboratory of Receptor Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China.,School of Pharmacy, University of Chinese Academy of Sciences, Beijing, China
| | - Yechun Xu
- CAS Key Laboratory of Receptor Research, Drug Discovery and Design Center, Shanghai Institute of Materia Medica, Chinese Academy of Sciences, Shanghai, China.,School of Pharmacy, University of Chinese Academy of Sciences, Beijing, China
| |
Collapse
|
16
|
Zihlmann P, Silbermann M, Sharpe T, Jiang X, Mühlethaler T, Jakob RP, Rabbani S, Sager CP, Frei P, Pang L, Maier T, Ernst B. KinITC-One Method Supports both Thermodynamic and Kinetic SARs as Exemplified on FimH Antagonists. Chemistry 2018; 24:13049-13057. [PMID: 29939458 DOI: 10.1002/chem.201802599] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2018] [Revised: 06/21/2018] [Indexed: 11/09/2022]
Abstract
Affinity data, such as dissociation constants (KD ) or inhibitory concentrations (IC50 ), are widely used in drug discovery. However, these parameters describe an equilibrium state, which is often not established in vivo due to pharmacokinetic effects and they are therefore not necessarily sufficient for evaluating drug efficacy. More accurate indicators for pharmacological activity are the kinetics of binding processes, as they shed light on the rate of formation of protein-ligand complexes and their half-life. Nonetheless, although highly desirable for medicinal chemistry programs, studies on structure-kinetic relationships (SKR) are still rare. With the recently introduced analytical tool kinITC this situation may change, since not only thermodynamic but also kinetic information of the binding process can be deduced from isothermal titration calorimetry (ITC) experiments. Using kinITC, ITC data of 29 mannosides binding to the bacterial adhesin FimH were re-analyzed to make their binding kinetics accessible. To validate these kinetic data, surface plasmon resonance (SPR) experiments were conducted. The kinetic analysis by kinITC revealed that the nanomolar affinities of the FimH antagonists arise from both (i) an optimized interaction between protein and ligand in the bound state (reduced off-rate constant koff ) and (ii) a stabilization of the transition state or a destabilization of the unbound state (increased on-rate constant kon ). Based on congeneric ligand modifications and structural input from co-crystal structures, a strong relationship between the formed hydrogen-bond network and koff could be concluded, whereas electrostatic interactions and conformational restrictions upon binding were found to have mainly an impact on kon .
Collapse
Affiliation(s)
- Pascal Zihlmann
- Institute of Molecular Pharmacy, Pharmazentrum, University of Basel, Klingelbergstrasse 50, 4056, Basel, Switzerland
| | - Marleen Silbermann
- Institute of Molecular Pharmacy, Pharmazentrum, University of Basel, Klingelbergstrasse 50, 4056, Basel, Switzerland
| | - Timothy Sharpe
- Biophysics Facility, Biozentrum, University of Basel, Klingelbergstrasse 70, 4056, Basel, Switzerland
| | - Xiaohua Jiang
- Institute of Molecular Pharmacy, Pharmazentrum, University of Basel, Klingelbergstrasse 50, 4056, Basel, Switzerland
| | - Tobias Mühlethaler
- Institute of Molecular Pharmacy, Pharmazentrum, University of Basel, Klingelbergstrasse 50, 4056, Basel, Switzerland
| | - Roman P Jakob
- Focal Area Structural Biology, Biozentrum, University of Basel, Klingelbergstrasse 70, 4056, Basel, Switzerland
| | - Said Rabbani
- Institute of Molecular Pharmacy, Pharmazentrum, University of Basel, Klingelbergstrasse 50, 4056, Basel, Switzerland
| | - Christoph P Sager
- Institute of Molecular Pharmacy, Pharmazentrum, University of Basel, Klingelbergstrasse 50, 4056, Basel, Switzerland
| | - Priska Frei
- Institute of Molecular Pharmacy, Pharmazentrum, University of Basel, Klingelbergstrasse 50, 4056, Basel, Switzerland
| | - Lijuan Pang
- Institute of Molecular Pharmacy, Pharmazentrum, University of Basel, Klingelbergstrasse 50, 4056, Basel, Switzerland
| | - Timm Maier
- Focal Area Structural Biology, Biozentrum, University of Basel, Klingelbergstrasse 70, 4056, Basel, Switzerland
| | - Beat Ernst
- Institute of Molecular Pharmacy, Pharmazentrum, University of Basel, Klingelbergstrasse 50, 4056, Basel, Switzerland
| |
Collapse
|
17
|
De Benedetti PG, Fanelli F. Computational modeling approaches to quantitative structure-binding kinetics relationships in drug discovery. Drug Discov Today 2018; 23:1396-1406. [PMID: 29574212 DOI: 10.1016/j.drudis.2018.03.010] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2017] [Revised: 02/22/2018] [Accepted: 03/19/2018] [Indexed: 11/22/2022]
Abstract
Simple comparative correlation analyses and quantitative structure-kinetics relationship (QSKR) models highlight the interplay of kinetic rates and binding affinity as an essential feature in drug design and discovery. The choice of the molecular series, and their structural variations, used in QSKR modeling is fundamental to understanding the mechanistic implications of ligand and/or drug-target binding and/or unbinding processes. Here, we discuss the implications of linear correlations between kinetic rates and binding affinity constants and the relevance of the computational approaches to QSKR modeling.
Collapse
Affiliation(s)
- Pier G De Benedetti
- Department of Life Sciences, University of Modena and Reggio Emilia, via Campi 103, 41125 Modena, Italy.
| | - Francesca Fanelli
- Department of Life Sciences, University of Modena and Reggio Emilia, via Campi 103, 41125 Modena, Italy; Center for Neuroscience and Neurotechnology, University of Modena and Reggio Emilia, via Campi 103, 41125 Modena, Italy
| |
Collapse
|
18
|
Schuetz DA, Richter L, Amaral M, Grandits M, Grädler U, Musil D, Buchstaller HP, Eggenweiler HM, Frech M, Ecker GF. Ligand Desolvation Steers On-Rate and Impacts Drug Residence Time of Heat Shock Protein 90 (Hsp90) Inhibitors. J Med Chem 2018; 61:4397-4411. [DOI: 10.1021/acs.jmedchem.8b00080] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Doris A. Schuetz
- Department of Pharmaceutical Chemistry, University of Vienna, UZA 2, Althanstrasse 14, 1090 Vienna, Austria
| | - Lars Richter
- Department of Pharmaceutical Chemistry, University of Vienna, UZA 2, Althanstrasse 14, 1090 Vienna, Austria
| | - Marta Amaral
- Discovery Technologies, Merck KGaA, Frankfurter Straße 250, 64293 Darmstadt, Germany
- Instituto de Biologia Experimental e Tecnológica, Oeiras, Portugal
| | - Melanie Grandits
- Department of Pharmaceutical Chemistry, University of Vienna, UZA 2, Althanstrasse 14, 1090 Vienna, Austria
| | - Ulrich Grädler
- Discovery Technologies, Merck KGaA, Frankfurter Straße 250, 64293 Darmstadt, Germany
| | - Djordje Musil
- Discovery Technologies, Merck KGaA, Frankfurter Straße 250, 64293 Darmstadt, Germany
| | | | | | - Matthias Frech
- Discovery Technologies, Merck KGaA, Frankfurter Straße 250, 64293 Darmstadt, Germany
| | - Gerhard F. Ecker
- Department of Pharmaceutical Chemistry, University of Vienna, UZA 2, Althanstrasse 14, 1090 Vienna, Austria
| |
Collapse
|
19
|
Linkuvienė V, Talibov VO, Danielson UH, Matulis D. Introduction of Intrinsic Kinetics of Protein–Ligand Interactions and Their Implications for Drug Design. J Med Chem 2018; 61:2292-2302. [DOI: 10.1021/acs.jmedchem.7b01408] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Vaida Linkuvienė
- Department of Biothermodynamics and Drug Design, Institute of Biotechnology, Vilnius University, Saulėtekio 7, Vilnius, LT-10257, Lithuania
| | - Vladimir O. Talibov
- Department of Chemistry - BMC, Uppsala University, Box 576, Uppsala, SE-751 23, Sweden
| | - U. Helena Danielson
- Department of Chemistry - BMC, Uppsala University, Box 576, Uppsala, SE-751 23, Sweden
- Science for Life Laboratory, Uppsala University, Uppsala, SE-751 23, Sweden
| | - Daumantas Matulis
- Department of Biothermodynamics and Drug Design, Institute of Biotechnology, Vilnius University, Saulėtekio 7, Vilnius, LT-10257, Lithuania
| |
Collapse
|
20
|
Wang Y, Edalji RP, Panchal SC, Sun C, Djuric SW, Vasudevan A. Are We There Yet? Applying Thermodynamic and Kinetic Profiling on Embryonic Ectoderm Development (EED) Hit-to-Lead Program. J Med Chem 2017; 60:8321-8335. [DOI: 10.1021/acs.jmedchem.7b00576] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Ying Wang
- AbbVie Inc., 1 North Waukegan Road, North
Chicago, Illinois 60064, United States
| | - Rohinton P. Edalji
- AbbVie Inc., 1 North Waukegan Road, North
Chicago, Illinois 60064, United States
| | - Sanjay C. Panchal
- AbbVie Inc., 1 North Waukegan Road, North
Chicago, Illinois 60064, United States
| | - Chaohong Sun
- AbbVie Inc., 1 North Waukegan Road, North
Chicago, Illinois 60064, United States
| | - Stevan W. Djuric
- AbbVie Inc., 1 North Waukegan Road, North
Chicago, Illinois 60064, United States
| | - Anil Vasudevan
- AbbVie Inc., 1 North Waukegan Road, North
Chicago, Illinois 60064, United States
| |
Collapse
|
21
|
Appel EA, Biedermann F, Hoogland D, Del Barrio J, Driscoll MD, Hay S, Wales DJ, Scherman OA. Decoupled Associative and Dissociative Processes in Strong yet Highly Dynamic Host-Guest Complexes. J Am Chem Soc 2017; 139:12985-12993. [PMID: 28661667 DOI: 10.1021/jacs.7b04821] [Citation(s) in RCA: 45] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Kinetics and thermodynamics in supramolecular systems are intimately linked, yet both are independently important for application in sensing assays and stimuli-responsive switching/self-healing of materials. Host-guest interactions are of particular interest in many water-based materials, sensing, and drug delivery applications. Herein we investigate the binding dynamics of a variety of electron-rich aromatic moieties forming hetero-ternary complexes with the macrocycle cucurbit[8]uril (CB[8]) and an auxiliary guest, dimethyl viologen, with high selectivity and equilibrium binding constants (Keq up to 1014 M-2). Using stopped-flow spectrofluorimetry, association rate constants were observed to approach the diffusion limit and were found to be insensitive to the structure of the guest. Conversely, the dissociation rate constants of the ternary complexes varied dramatically with the guest structure and were correlated with the thermodynamic binding selectivity. Hence differing molecular features were found to contribute to the associative and dissociative processes, mimicking naturally occurring reactions and giving rise to a decoupling of these kinetic parameters. Moreover, we demonstrate the ability to exploit these phenomena and selectively perturb the associative process with external stimuli (e.g., viscosity and pressure). Significantly, these complexes exhibit increased binding equilibria with increasing pressure, with important implications for the application of the CB[8] ternary complex for the formation of hydrogels, as these gels exhibit unprecedented pressure-insensitive rheological properties. A high degree of flexibility therefore exists in the design of host-guest systems with tunable kinetic and thermodynamic parameters for tailor-made applications across a broad range of fields.
Collapse
Affiliation(s)
- Eric A Appel
- Department of Materials Science and Engineering, Stanford University , Stanford California 94305, United States.,Department of Chemistry, University of Cambridge , Lensfield Road, Cambridge CB2 1EW, U.K
| | - Frank Biedermann
- Institute of Nanotechnology, Karlsruhe Institute of Technology (KIT) , Hermann-von-Helmholtz Platz 1, 76344 Eggenstein-Leopoldshafen, Germany.,Department of Chemistry, University of Cambridge , Lensfield Road, Cambridge CB2 1EW, U.K
| | - Dominique Hoogland
- Department of Chemistry, University of Cambridge , Lensfield Road, Cambridge CB2 1EW, U.K
| | - Jesús Del Barrio
- Schlumberger Gould Research Center, High Cross , Madingley Road, Cambridge CB3 0EL, U.K
| | - Max D Driscoll
- Manchester Institute of Biotechnology and School of Chemistry, University of Manchester , 131 Princess Street, Manchester M1 7DN, U.K
| | - Sam Hay
- Manchester Institute of Biotechnology and School of Chemistry, University of Manchester , 131 Princess Street, Manchester M1 7DN, U.K
| | - David J Wales
- Department of Chemistry, University of Cambridge , Lensfield Road, Cambridge CB2 1EW, U.K
| | - Oren A Scherman
- Department of Chemistry, University of Cambridge , Lensfield Road, Cambridge CB2 1EW, U.K
| |
Collapse
|
22
|
Linkuvienė V, Krainer G, Chen WY, Matulis D. Isothermal titration calorimetry for drug design: Precision of the enthalpy and binding constant measurements and comparison of the instruments. Anal Biochem 2016; 515:61-64. [DOI: 10.1016/j.ab.2016.10.005] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2016] [Accepted: 10/03/2016] [Indexed: 10/20/2022]
|
23
|
Renaud JP, Chung CW, Danielson UH, Egner U, Hennig M, Hubbard RE, Nar H. Biophysics in drug discovery: impact, challenges and opportunities. Nat Rev Drug Discov 2016; 15:679-98. [PMID: 27516170 DOI: 10.1038/nrd.2016.123] [Citation(s) in RCA: 209] [Impact Index Per Article: 26.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Over the past 25 years, biophysical technologies such as X-ray crystallography, nuclear magnetic resonance spectroscopy, surface plasmon resonance spectroscopy and isothermal titration calorimetry have become key components of drug discovery platforms in many pharmaceutical companies and academic laboratories. There have been great improvements in the speed, sensitivity and range of possible measurements, providing high-resolution mechanistic, kinetic, thermodynamic and structural information on compound-target interactions. This Review provides a framework to understand this evolution by describing the key biophysical methods, the information they can provide and the ways in which they can be applied at different stages of the drug discovery process. We also discuss the challenges for current technologies and future opportunities to use biophysical methods to solve drug discovery problems.
Collapse
Affiliation(s)
- Jean-Paul Renaud
- NovAliX, Boulevard Sébastien Brant, 67405 Illkirch Cedex, France.,Institut de Génétique et Biologie Moléculaire et Cellulaire, CNRS UMR7104/INSERM U964/Université de Strasbourg, 1 rue Laurent Fries - BP10142, 67404 Illkirch Cedex, France.,RiboStruct, 15 rue Neuve, 67540 Ostwald, France
| | - Chun-Wa Chung
- GlaxoSmithKline R&D, Gunnels Wood Road, Stevenage, SG1 2NY, UK
| | - U Helena Danielson
- Department of Chemistry - BMC and Science for Life Laboratory, Drug Discovery &Development Platform, Uppsala University, SE-751 05 Uppsala, Sweden.,Beactica AB, Uppsala Business Park, 754 50 Uppsala, Sweden
| | - Ursula Egner
- Bayer Pharma AG, Müllerstrasse 178, 13353 Berlin, Germany
| | - Michael Hennig
- Hoffmann-La Roche Ltd, Grenzacherstrasse 124, 4070 Basel, Switzerland.,leadXpro AG, PARK INNOVAARE, CH-5234 Villigen, Switzerland
| | - Roderick E Hubbard
- University of York, Heslington, York, YO10 5DD, UK.,Vernalis (R&D), Granta Park, Cambridge, CB21 6GB, UK
| | - Herbert Nar
- Boehringer Ingelheim GmbH &Co. KG, Birkendorfer Strasse 65, 88400 Biberach, Germany
| |
Collapse
|
24
|
Impact, determination and prediction of drug-receptor residence times for GPCRs. Curr Opin Pharmacol 2016; 30:22-26. [PMID: 27428776 DOI: 10.1016/j.coph.2016.07.004] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2016] [Revised: 07/03/2016] [Accepted: 07/04/2016] [Indexed: 01/02/2023]
Abstract
The residence time of a ligand on a GPCR of interest has become an optimization parameter in many examples in drug design. Long residence times can counterbalance unfavorable pharmacokinetic parameters, contributing to compound safety, and short residence times can be a tool to avoid target related side effects. Unlike the prediction and interpretation of the structure-activity relationship (SAR) of a ligand class on a receptor, the understanding and prediction of the structure-kinetics relationship (SKR) is much more demanding. Experimental and computational approaches are described, which serve to either rationalize SKR or to predict the kinetic parameters such as on-rates and off-rates.
Collapse
|
25
|
Abstract
The dynamics of protein binding pockets are crucial for their interaction specificity. Structural flexibility allows proteins to adapt to their individual molecular binding partners and facilitates the binding process. This implies the necessity to consider protein internal motion in determining and predicting binding properties and in designing new binders. Although accounting for protein dynamics presents a challenge for computational approaches, it expands the structural and physicochemical space for compound design and thus offers the prospect of improved binding specificity and selectivity. A cavity on the surface or in the interior of a protein that possesses suitable properties for binding a ligand is usually referred to as a binding pocket. The set of amino acid residues around a binding pocket determines its physicochemical characteristics and, together with its shape and location in a protein, defines its functionality. Residues outside the binding site can also have a long-range effect on the properties of the binding pocket. Cavities with similar functionalities are often conserved across protein families. For example, enzyme active sites are usually concave surfaces that present amino acid residues in a suitable configuration for binding low molecular weight compounds. Macromolecular binding pockets, on the other hand, are located on the protein surface and are often shallower. The mobility of proteins allows the opening, closing, and adaptation of binding pockets to regulate binding processes and specific protein functionalities. For example, channels and tunnels can exist permanently or transiently to transport compounds to and from a binding site. The influence of protein flexibility on binding pockets can vary from small changes to an already existent pocket to the formation of a completely new pocket. Here, we review recent developments in computational methods to detect and define binding pockets and to study pocket dynamics. We introduce five different classes of protein pocket dynamics: (1) appearance/disappearance of a subpocket in an existing pocket; (2) appearance/disappearance of an adjacent pocket on the protein surface in the direct vicinity of an already existing pocket; (3) pocket breathing, which may be caused by side-chain fluctuations or backbone or interdomain vibrational motion; (4) opening/closing of a channel or tunnel, connecting a pocket inside the protein with solvent, including lid motion; and (5) the appearance/disappearance of an allosteric pocket at a site on a protein distinct from an already existing pocket with binding of a ligand to the allosteric binding site affecting the original pocket. We suggest that the class of pocket dynamics, as well as the type and extent of protein motion affecting the binding pocket, should be factors considered in choosing the most appropriate computational approach to study a given binding pocket. Furthermore, we examine the relationship between pocket dynamics classes and induced fit, conformational selection, and gating models of ligand binding on binding kinetics and thermodynamics. We discuss the implications of protein binding pocket dynamics for drug design and conclude with potential future directions for computational analysis of protein binding pocket dynamics.
Collapse
Affiliation(s)
- Antonia Stank
- Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
| | - Daria B. Kokh
- Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
| | - Jonathan C. Fuller
- Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
| | - Rebecca C. Wade
- Heidelberg Institute for Theoretical Studies (HITS), Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
- Center
for Molecular Biology of the University of Heidelberg (ZMBH), DKFZ-ZMBH Alliance, Im Neuenheimer Feld 282, 69120 Heidelberg, Germany
| |
Collapse
|
26
|
Lawton G, Nussbaumer P. The Evolving Role of the Medicinal Chemist. PROGRESS IN MEDICINAL CHEMISTRY 2016; 55:193-226. [PMID: 26852936 DOI: 10.1016/bs.pmch.2015.11.001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Key Words] [Subscribe] [Scholar Register] [Indexed: 10/21/2022]
Affiliation(s)
- Geoff Lawton
- Gardenfields, St Ippolyts, Hertfordshire, United Kingdom
| | | |
Collapse
|
27
|
Biggin PC, Aldeghi M, Bodkin MJ, Heifetz A. Beyond Membrane Protein Structure: Drug Discovery, Dynamics and Difficulties. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2016; 922:161-181. [PMID: 27553242 DOI: 10.1007/978-3-319-35072-1_12] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Most of the previous content of this book has focused on obtaining the structures of membrane proteins. In this chapter we explore how those structures can be further used in two key ways. The first is their use in structure based drug design (SBDD) and the second is how they can be used to extend our understanding of their functional activity via the use of molecular dynamics. Both aspects now heavily rely on computations. This area is vast, and alas, too large to consider in depth in a single book chapter. Thus where appropriate we have referred the reader to recent reviews for deeper assessment of the field. We discuss progress via the use of examples from two main drug target areas; G-protein coupled receptors (GPCRs) and ion channels. We end with a discussion of some of the main challenges in the area.
Collapse
Affiliation(s)
- Philip C Biggin
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, UK.
| | - Matteo Aldeghi
- Department of Biochemistry, University of Oxford, South Parks Road, Oxford, OX1 3QU, UK
| | - Michael J Bodkin
- Evotec Ltd, 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK
| | - Alexander Heifetz
- Evotec Ltd, 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire, OX14 4RZ, UK
| |
Collapse
|
28
|
Thermodynamics of protein–ligand interactions as a reference for computational analysis: how to assess accuracy, reliability and relevance of experimental data. J Comput Aided Mol Des 2015; 29:867-83. [DOI: 10.1007/s10822-015-9867-y] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2015] [Accepted: 09/05/2015] [Indexed: 12/11/2022]
|
29
|
Ekins S, Litterman NK, Lipinski CA, Bunin BA. Thermodynamic Proxies to Compensate for Biases in Drug Discovery Methods. Pharm Res 2015; 33:194-205. [DOI: 10.1007/s11095-015-1779-y] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2015] [Accepted: 08/13/2015] [Indexed: 11/24/2022]
|
30
|
Stockdale TP, Williams CM. Pharmaceuticals that contain polycyclic hydrocarbon scaffolds. Chem Soc Rev 2015; 44:7737-63. [DOI: 10.1039/c4cs00477a] [Citation(s) in RCA: 137] [Impact Index Per Article: 15.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
This review comprehensively explores approved pharmaceutical compounds that contain polycyclic scaffolds and the properties that these skeletons convey.
Collapse
Affiliation(s)
- Tegan P. Stockdale
- School of Chemistry and Molecular Biosciences
- University of Queensland
- St Lucia
- Australia
| | - Craig M. Williams
- School of Chemistry and Molecular Biosciences
- University of Queensland
- St Lucia
- Australia
| |
Collapse
|