Combined Free-Energy Calculation and Machine Learning Methods for Understanding Ligand Unbinding Kinetics
- PMID: 35195418
- PMCID: PMC9097281
- DOI: 10.1021/acs.jctc.1c00924
Combined Free-Energy Calculation and Machine Learning Methods for Understanding Ligand Unbinding Kinetics
Abstract
The determination of drug residence times, which define the time an inhibitor is in complex with its target, is a fundamental part of the drug discovery process. Synthesis and experimental measurements of kinetic rate constants are, however, expensive and time consuming. In this work, we aimed to obtain drug residence times computationally. Furthermore, we propose a novel algorithm to identify molecular design objectives based on ligand unbinding kinetics. We designed an enhanced sampling technique to accurately predict the free-energy profiles of the ligand unbinding process, focusing on the free-energy barrier for unbinding. Our method first identifies unbinding paths determining a corresponding set of internal coordinates (ICs) that form contacts between the protein and the ligand; it then iteratively updates these interactions during a series of biased molecular dynamics (MD) simulations to reveal the ICs that are important for the whole of the unbinding process. Subsequently, we performed finite-temperature string simulations to obtain the free-energy barrier for unbinding using the set of ICs as a complex reaction coordinate. Importantly, we also aimed to enable the further design of drugs focusing on improved residence times. To this end, we developed a supervised machine learning (ML) approach with inputs from unbiased "downhill" trajectories initiated near the transition state (TS) ensemble of the string unbinding path. We demonstrate that our ML method can identify key ligand-protein interactions driving the system through the TS. Some of the most important drugs for cancer treatment are kinase inhibitors. One of these kinase targets is cyclin-dependent kinase 2 (CDK2), an appealing target for anticancer drug development. Here, we tested our method using two different CDK2 inhibitors for the potential further development of these compounds. We compared the free-energy barriers obtained from our calculations with those observed in available experimental data. We highlighted important interactions at the distal ends of the ligands that can be targeted for improved residence times. Our method provides a new tool to determine unbinding rates and to identify key structural features of the inhibitors that can be used as starting points for novel design strategies in drug discovery.
Conflict of interest statement
The authors declare no competing financial interest.
Figures
![Figure 1](https://www.ncbi.nlm.nih.gov/pmc/articles/instance/9097281/bin/ct1c00924_0002.gif)
![Figure 2](https://www.ncbi.nlm.nih.gov/pmc/articles/instance/9097281/bin/ct1c00924_0003.gif)
![Figure 3](https://www.ncbi.nlm.nih.gov/pmc/articles/instance/9097281/bin/ct1c00924_0004.gif)
![Figure 4](https://www.ncbi.nlm.nih.gov/pmc/articles/instance/9097281/bin/ct1c00924_0005.gif)
![Figure 5](https://www.ncbi.nlm.nih.gov/pmc/articles/instance/9097281/bin/ct1c00924_0006.gif)
![Figure 6](https://www.ncbi.nlm.nih.gov/pmc/articles/instance/9097281/bin/ct1c00924_0007.gif)
![Figure 7](https://www.ncbi.nlm.nih.gov/pmc/articles/instance/9097281/bin/ct1c00924_0008.gif)
![Figure 8](https://www.ncbi.nlm.nih.gov/pmc/articles/instance/9097281/bin/ct1c00924_0009.gif)
Similar articles
-
Predicting Protein-Ligand Binding and Unbinding Kinetics with Biased MD Simulations and Coarse-Graining of Dynamics: Current State and Challenges.J Chem Inf Model. 2023 May 22;63(10):2902-2910. doi: 10.1021/acs.jcim.3c00151. Epub 2023 May 3. J Chem Inf Model. 2023. PMID: 37133392 Review.
-
The prediction of protein-ligand unbinding for modern drug discovery.Expert Opin Drug Discov. 2022 Feb;17(2):191-205. doi: 10.1080/17460441.2022.2002298. Epub 2021 Nov 12. Expert Opin Drug Discov. 2022. PMID: 34731059 Review.
-
Contact Map Fingerprints of Protein-Ligand Unbinding Trajectories Reveal Mechanisms Determining Residence Times Computed from Scaled Molecular Dynamics.J Chem Theory Comput. 2021 Oct 12;17(10):6522-6535. doi: 10.1021/acs.jctc.1c00453. Epub 2021 Sep 8. J Chem Theory Comput. 2021. PMID: 34494849
-
Estimation of Protein-Ligand Unbinding Kinetics Using Non-Equilibrium Targeted Molecular Dynamics Simulations.J Chem Inf Model. 2019 Dec 23;59(12):5135-5147. doi: 10.1021/acs.jcim.9b00592. Epub 2019 Nov 22. J Chem Inf Model. 2019. PMID: 31697501
-
Variational implicit-solvent predictions of the dry-wet transition pathways for ligand-receptor binding and unbinding kinetics.Proc Natl Acad Sci U S A. 2019 Jul 23;116(30):14989-14994. doi: 10.1073/pnas.1902719116. Epub 2019 Jul 3. Proc Natl Acad Sci U S A. 2019. PMID: 31270236 Free PMC article.
Cited by
-
Surface-Accelerated String Method for Locating Minimum Free Energy Paths.J Chem Theory Comput. 2024 Mar 12;20(5):2058-2073. doi: 10.1021/acs.jctc.3c01401. Epub 2024 Feb 17. J Chem Theory Comput. 2024. PMID: 38367218
-
How Robust Is the Ligand Binding Transition State?J Am Chem Soc. 2023 Nov 22;145(46):25318-25331. doi: 10.1021/jacs.3c08940. Epub 2023 Nov 9. J Am Chem Soc. 2023. PMID: 37943667 Free PMC article.
-
Perspectives on Ligand/Protein Binding Kinetics Simulations: Force Fields, Machine Learning, Sampling, and User-Friendliness.J Chem Theory Comput. 2023 Sep 26;19(18):6047-6061. doi: 10.1021/acs.jctc.3c00641. Epub 2023 Sep 1. J Chem Theory Comput. 2023. PMID: 37656199 Free PMC article. Review.
-
Investigating the Unbinding of Muscarinic Antagonists from the Muscarinic 3 Receptor.J Chem Theory Comput. 2023 Aug 8;19(15):5260-5272. doi: 10.1021/acs.jctc.3c00023. Epub 2023 Jul 17. J Chem Theory Comput. 2023. PMID: 37458730 Free PMC article.
-
Simulation of ligand dissociation kinetics from the protein kinase PYK2.J Comput Chem. 2022 Oct 30;43(28):1911-1922. doi: 10.1002/jcc.26991. Epub 2022 Sep 8. J Comput Chem. 2022. PMID: 36073605 Free PMC article.
References
-
- Copeland R. A.Evaluation of Enzyme Inhibitors in Drug Discovery: A Guide for Medicinal Chemists and Pharmacologists, 2nd ed.; John Wiley and Sons, 2013. - PubMed
MeSH terms
Substances
LinkOut - more resources
Full Text Sources