1
|
Knockenhauer KE, Copeland RA. The importance of binding kinetics and drug-target residence time in pharmacology. Br J Pharmacol 2024; 181:4103-4116. [PMID: 37160660 DOI: 10.1111/bph.16104] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2023] [Revised: 04/19/2023] [Accepted: 04/27/2023] [Indexed: 05/11/2023] Open
Abstract
A dominant assumption in pharmacology throughout the 20th century has been that in vivo target occupancy-and attendant pharmacodynamics-depends on the systemic concentration of drug relative to the equilibrium dissociation constant for the drug-target complex. In turn, the duration of pharmacodynamics is temporally linked to the systemic pharmacokinetics of the drug. Yet, there are many examples of drugs for which pharmacodynamic effect endures long after the systemic concentration of a drug has waned to (equilibrium) insignificant levels. To reconcile such data, the drug-target residence time model was formulated, positing that it is the lifetime (or residence time) of the binary drug-target complex, and not its equilibrium affinity per se, that determines the extent and duration of drug pharmacodynamics. Here, we review this model, its evolution over time, and its applications to natural ligand-macromolecule biology and synthetic drug-target pharmacology.
Collapse
|
2
|
Du Y. Binding Curve Viewer: Visualizing the Equilibrium and Kinetics of Protein-Ligand Binding and Competitive Binding. J Chem Inf Model 2024; 64:4180-4192. [PMID: 38720179 PMCID: PMC11134506 DOI: 10.1021/acs.jcim.4c00130] [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: 01/24/2024] [Revised: 04/21/2024] [Accepted: 04/25/2024] [Indexed: 05/28/2024]
Abstract
Understanding the thermodynamics and kinetics of the protein-ligand interaction is essential for biologists and pharmacologists. To visualize the equilibrium and kinetics of the binding reaction with 1:1 stoichiometry and no cooperativity, we obtained the exact relationship of the concentration of the protein-ligand complex and the time in the second-order binding process and numerically simulated the process of competitive binding. First, two common concerns in measuring protein-ligand interactions were focused on how to avoid the titration regime and how to establish the appropriate incubation time. Then, we gave examples of how the commonly used experimental conditions of [L]0 ≫ [P]0 and [I]0 ≫ [P]0 affected the estimation of the kinetic and thermodynamic properties. Theoretical inhibition curves were calculated, and the apparent IC50 and IC50 were estimated accordingly under predefined conditions. Using the estimated apparent IC50, we compared the apparent Ki and Ki calculated by using the Cheng-Prusoff equation, Lin-Riggs equation, and Wang's group equation. We also applied our tools to simulate high-throughput screening and compare the results of real experiments. The visualization tool for simulating the saturation experiment, kinetic experiments of binding and competitive binding, and inhibition curve, "Binding Curve Viewer," is available at www.eplatton.net/binding-curve-viewer.
Collapse
Affiliation(s)
- Yu Du
- Department
of Clinical Laboratory, The Second Affiliated
Hospital of Jiaxing University, Huancheng North Road 1518, Jiaxing, Zhejiang 314000, China
- The
Key Laboratory, The Second Affiliated Hospital
of Jiaxing University, Huancheng North Road 1518, Jiaxing, Zhejiang 314000, China
| |
Collapse
|
3
|
Xu C, Zhang X, Zhao L, Verkhivker GM, Bai F. Accurate Characterization of Binding Kinetics and Allosteric Mechanisms for the HSP90 Chaperone Inhibitors Using AI-Augmented Integrative Biophysical Studies. JACS AU 2024; 4:1632-1645. [PMID: 38665669 PMCID: PMC11040708 DOI: 10.1021/jacsau.4c00123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/09/2024] [Revised: 03/15/2024] [Accepted: 03/18/2024] [Indexed: 04/28/2024]
Abstract
The binding kinetics of drugs to their targets are gradually being recognized as a crucial indicator of the efficacy of drugs in vivo, leading to the development of various computational methods for predicting the binding kinetics in recent years. However, compared with the prediction of binding affinity, the underlying structure and dynamic determinants of binding kinetics are more complicated. Efficient and accurate methods for predicting binding kinetics are still lacking. In this study, quantitative structure-kinetics relationship (QSKR) models were developed using 132 inhibitors targeting the ATP binding domain of heat shock protein 90α (HSP90α) to predict the dissociation rate constant (koff), enabling a direct assessment of the drug-target residence time. These models demonstrated good predictive performance, where hydrophobic and hydrogen bond interactions significantly influence the koff prediction. In subsequent applications, our models were used to assist in the discovery of new inhibitors for the N-terminal domain of HSP90α (N-HSP90α), demonstrating predictive capabilities on an experimental validation set with a new scaffold. In X-ray crystallography experiments, the loop-middle conformation of apo N-HSP90α was observed for the first time (previously, the loop-middle conformation had only been observed in holo-N-HSP90α structures). Interestingly, we observed different conformations of apo N-HSP90α simultaneously in an asymmetric unit, which was also observed in a holo-N-HSP90α structure, suggesting an equilibrium of conformations between different states in solution, which could be one of the determinants affecting the binding kinetics of the ligand. Different ligands can undergo conformational selection or alter the equilibrium of conformations, inducing conformational rearrangements and resulting in different effects on binding kinetics. We then used molecular dynamics simulations to describe conformational changes of apo N-HSP90α in different conformational states. In summary, the study of the binding kinetics and molecular mechanisms of N-HSP90α provides valuable information for the development of more targeted therapeutic approaches.
Collapse
Affiliation(s)
- Chao Xu
- Shanghai
Institute for Advanced Immunochemical Studies and School of Life Science
and Technology, ShanghaiTech University, 393 Middle Huaxia Road, Shanghai 201210, China
| | - Xianglei Zhang
- Shanghai
Institute for Advanced Immunochemical Studies and School of Life Science
and Technology, ShanghaiTech University, 393 Middle Huaxia Road, Shanghai 201210, China
| | - Lianghao Zhao
- Shanghai
Institute for Advanced Immunochemical Studies and School of Life Science
and Technology, ShanghaiTech University, 393 Middle Huaxia Road, Shanghai 201210, China
| | - Gennady M. Verkhivker
- Keck
Center for Science and Engineering, Graduate Program in Computational
and Data Sciences, Schmid College of Science and Technology, Chapman University, Orange, California 92866, United States
- Department
of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Irvine, California 92618, United States
| | - Fang Bai
- Shanghai
Institute for Advanced Immunochemical Studies and School of Life Science
and Technology, ShanghaiTech University, 393 Middle Huaxia Road, Shanghai 201210, China
- School
of Information Science and Technology, ShanghaiTech
University, 393 Middle Huaxia Road, Shanghai 201210, China
- Shanghai
Clinical Research and Trial Center, Shanghai 201210, China
| |
Collapse
|
4
|
Li Z, Huang R, Xia M, Patterson TA, Hong H. Fingerprinting Interactions between Proteins and Ligands for Facilitating Machine Learning in Drug Discovery. Biomolecules 2024; 14:72. [PMID: 38254672 PMCID: PMC10813698 DOI: 10.3390/biom14010072] [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/02/2023] [Revised: 12/26/2023] [Accepted: 12/28/2023] [Indexed: 01/24/2024] Open
Abstract
Molecular recognition is fundamental in biology, underpinning intricate processes through specific protein-ligand interactions. This understanding is pivotal in drug discovery, yet traditional experimental methods face limitations in exploring the vast chemical space. Computational approaches, notably quantitative structure-activity/property relationship analysis, have gained prominence. Molecular fingerprints encode molecular structures and serve as property profiles, which are essential in drug discovery. While two-dimensional (2D) fingerprints are commonly used, three-dimensional (3D) structural interaction fingerprints offer enhanced structural features specific to target proteins. Machine learning models trained on interaction fingerprints enable precise binding prediction. Recent focus has shifted to structure-based predictive modeling, with machine-learning scoring functions excelling due to feature engineering guided by key interactions. Notably, 3D interaction fingerprints are gaining ground due to their robustness. Various structural interaction fingerprints have been developed and used in drug discovery, each with unique capabilities. This review recapitulates the developed structural interaction fingerprints and provides two case studies to illustrate the power of interaction fingerprint-driven machine learning. The first elucidates structure-activity relationships in β2 adrenoceptor ligands, demonstrating the ability to differentiate agonists and antagonists. The second employs a retrosynthesis-based pre-trained molecular representation to predict protein-ligand dissociation rates, offering insights into binding kinetics. Despite remarkable progress, challenges persist in interpreting complex machine learning models built on 3D fingerprints, emphasizing the need for strategies to make predictions interpretable. Binding site plasticity and induced fit effects pose additional complexities. Interaction fingerprints are promising but require continued research to harness their full potential.
Collapse
Affiliation(s)
- Zoe Li
- National Center for Toxicological Research, US Food and Drug Administration, Jefferson, AR 72079, USA; (Z.L.); (T.A.P.)
| | - Ruili Huang
- National Center for Advancing Translational Sciences, National Institutes of Health, Bethesda, MD 20892, USA; (R.H.); (M.X.)
| | - Menghang Xia
- National Center for Advancing Translational Sciences, National Institutes of Health, Bethesda, MD 20892, USA; (R.H.); (M.X.)
| | - Tucker A. Patterson
- National Center for Toxicological Research, US Food and Drug Administration, Jefferson, AR 72079, USA; (Z.L.); (T.A.P.)
| | - Huixiao Hong
- National Center for Toxicological Research, US Food and Drug Administration, Jefferson, AR 72079, USA; (Z.L.); (T.A.P.)
| |
Collapse
|
5
|
Zhou F, Yin S, Xiao Y, Lin Z, Fu W, Zhang YJ. Structure-Kinetic Relationship for Drug Design Revealed by a PLS Model with Retrosynthesis-Based Pre-Trained Molecular Representation and Molecular Dynamics Simulation. ACS OMEGA 2023; 8:18312-18322. [PMID: 37251166 PMCID: PMC10210189 DOI: 10.1021/acsomega.3c02294] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/05/2023] [Accepted: 04/26/2023] [Indexed: 05/31/2023]
Abstract
Drug design based on kinetic properties is growing in application. Here, we applied retrosynthesis-based pre-trained molecular representation (RPM) in machine learning (ML) to train 501 inhibitors of 55 proteins and successfully predicted the dissociation rate constant (koff) values of 38 inhibitors from an independent dataset for the N-terminal domain of heat shock protein 90α (N-HSP90). Our RPM molecular representation outperforms other pre-trained molecular representations such as GEM, MPG, and general molecular descriptors from RDKit. Furthermore, we optimized the accelerated molecular dynamics to calculate the relative retention time (RT) for the 128 inhibitors of N-HSP90 and obtained the protein-ligand interaction fingerprints (IFPs) on their dissociation pathways and their influencing weights on the koff value. We observed a high correlation among the simulated, predicted, and experimental -log(koff) values. Combining ML, molecular dynamics (MD) simulation, and IFPs derived from accelerated MD helps design a drug for specific kinetic properties and selectivity profiles to the target of interest. To further validate our koff predictive ML model, we tested our model on two new N-HSP90 inhibitors, which have experimental koff values and are not in our ML training dataset. The predicted koff values are consistent with experimental data, and the mechanism of their kinetic properties can be explained by IFPs, which shed light on the nature of their selectivity against N-HSP90 protein. We believe that the ML model described here is transferable to predict koff of other proteins and will enhance the kinetics-based drug design endeavor.
Collapse
|
6
|
Wang J, Do HN, Koirala K, Miao Y. Predicting Biomolecular Binding Kinetics: A Review. J Chem Theory Comput 2023; 19:2135-2148. [PMID: 36989090 DOI: 10.1021/acs.jctc.2c01085] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/30/2023]
Abstract
Biomolecular binding kinetics including the association (kon) and dissociation (koff) rates are critical parameters for therapeutic design of small-molecule drugs, peptides, and antibodies. Notably, the drug molecule residence time or dissociation rate has been shown to correlate with their efficacies better than binding affinities. A wide range of modeling approaches including quantitative structure-kinetic relationship models, Molecular Dynamics simulations, enhanced sampling, and Machine Learning has been developed to explore biomolecular binding and dissociation mechanisms and predict binding kinetic rates. Here, we review recent advances in computational modeling of biomolecular binding kinetics, with an outlook for future improvements.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Hung N Do
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Kushal Koirala
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| |
Collapse
|
7
|
Sohraby F, Nunes-Alves A. Advances in computational methods for ligand binding kinetics. Trends Biochem Sci 2022; 48:437-449. [PMID: 36566088 DOI: 10.1016/j.tibs.2022.11.003] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2022] [Revised: 11/16/2022] [Accepted: 11/29/2022] [Indexed: 12/24/2022]
Abstract
Binding kinetic parameters can be correlated with drug efficacy, which in recent years led to the development of various computational methods for predicting binding kinetic rates and gaining insight into protein-drug binding paths and mechanisms. In this review, we introduce and compare computational methods recently developed and applied to two systems, trypsin-benzamidine and kinase-inhibitor complexes. Methods involving enhanced sampling in molecular dynamics simulations or machine learning can be used not only to predict kinetic rates, but also to reveal factors modulating the duration of residence times, selectivity, and drug resistance to mutations. Methods which require less computational time to make predictions are highlighted, and suggestions to reduce the error of computed kinetic rates are presented.
Collapse
Affiliation(s)
- Farzin Sohraby
- Institute of Chemistry, Technische Universität Berlin, 10623 Berlin, Germany
| | - Ariane Nunes-Alves
- Institute of Chemistry, Technische Universität Berlin, 10623 Berlin, Germany.
| |
Collapse
|