Shradheya R. R. Gupta, Rashmi Rameshwari, Indrakant Kumar Singh
{"title":"Discovery of a potent ROR1 inhibitor using μs-scale MD simulations, wt-metadynamics, and absolute binding free energy calculations","authors":"Shradheya R. R. Gupta, Rashmi Rameshwari, Indrakant Kumar Singh","doi":"10.1007/s00894-025-06508-3","DOIUrl":null,"url":null,"abstract":"<div><h3>Context</h3><p>Receptor tyrosine kinase-like orphan receptor 1 (ROR1) is a cancer-associated pseudokinase with low expression in normal adult tissues but elevated levels in various malignancies, making it a promising therapeutic target. Among ~ 4 million compounds, CHEMBL3926946 emerged as the most promising candidate, demonstrating a persistent binding pose and a well-defined free energy basin. Well-tempered metadynamics (wt-MetaD) revealed a deep minimum of 26.00 ± 2.44 kcal/mol, indicating a highly stable interaction. CHEMBL3926946 exhibited a favourable Absolute Binding Free Energy Perturbation (ABFEP) of − 16.52 ± 0.37 kcal/mol, significantly outperforming the inhibitor Ponatinib (− 8.67 ± 0.94 kcal/mol), supported by persistent interactions with GLU523 and LEU479. This study highlights CHEMBL3926946 as a robust lead for ROR1-targeted cancer therapy and emphasizes the utility of combining wt-MetaD and ABFEP for reliable hit prioritization.</p><h3>Methods</h3><p>We employed a multilayered in silico pipeline integrating high-throughput virtual screening, long-timescale molecular dynamics (MD), wt-MetaD, and ABFEP. Ligands and protein were prepared using the OPLS2005 force field, and all stages up to wt-MetaD were conducted in Maestro (v12.8.117) using the same force field. A library of ~ 4 million compounds yielded 137 candidates, Further shortlisted via MD. 7 high-confidence molecules underwent 5 independent MD replicates with randomized seeds to ensure statistical robustness. The top 3 compounds were validated by 1 μs (1000 ns) simulations to assess long-term conformational stability and wt-MetaD to reveal deep minimum. ABFEP calculations were performed using the CGenFF force field in NAMD 3.0. We benchmarked ABFEP protocol against experimentally validated ligands, successfully reproducing experimental binding free energies (ΔG), confirming the protocol’s predictive reliability.</p><h3>Graphical Abstract</h3><div><figure><div><div><picture><source><img></source></picture></div></div></figure></div></div>","PeriodicalId":651,"journal":{"name":"Journal of Molecular Modeling","volume":"31 10","pages":""},"PeriodicalIF":2.5000,"publicationDate":"2025-09-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Molecular Modeling","FirstCategoryId":"92","ListUrlMain":"https://link.springer.com/article/10.1007/s00894-025-06508-3","RegionNum":4,"RegionCategory":"化学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q4","JCRName":"BIOCHEMISTRY & MOLECULAR BIOLOGY","Score":null,"Total":0}
引用次数: 0
Abstract
Context
Receptor tyrosine kinase-like orphan receptor 1 (ROR1) is a cancer-associated pseudokinase with low expression in normal adult tissues but elevated levels in various malignancies, making it a promising therapeutic target. Among ~ 4 million compounds, CHEMBL3926946 emerged as the most promising candidate, demonstrating a persistent binding pose and a well-defined free energy basin. Well-tempered metadynamics (wt-MetaD) revealed a deep minimum of 26.00 ± 2.44 kcal/mol, indicating a highly stable interaction. CHEMBL3926946 exhibited a favourable Absolute Binding Free Energy Perturbation (ABFEP) of − 16.52 ± 0.37 kcal/mol, significantly outperforming the inhibitor Ponatinib (− 8.67 ± 0.94 kcal/mol), supported by persistent interactions with GLU523 and LEU479. This study highlights CHEMBL3926946 as a robust lead for ROR1-targeted cancer therapy and emphasizes the utility of combining wt-MetaD and ABFEP for reliable hit prioritization.
Methods
We employed a multilayered in silico pipeline integrating high-throughput virtual screening, long-timescale molecular dynamics (MD), wt-MetaD, and ABFEP. Ligands and protein were prepared using the OPLS2005 force field, and all stages up to wt-MetaD were conducted in Maestro (v12.8.117) using the same force field. A library of ~ 4 million compounds yielded 137 candidates, Further shortlisted via MD. 7 high-confidence molecules underwent 5 independent MD replicates with randomized seeds to ensure statistical robustness. The top 3 compounds were validated by 1 μs (1000 ns) simulations to assess long-term conformational stability and wt-MetaD to reveal deep minimum. ABFEP calculations were performed using the CGenFF force field in NAMD 3.0. We benchmarked ABFEP protocol against experimentally validated ligands, successfully reproducing experimental binding free energies (ΔG), confirming the protocol’s predictive reliability.
期刊介绍:
The Journal of Molecular Modeling focuses on "hardcore" modeling, publishing high-quality research and reports. Founded in 1995 as a purely electronic journal, it has adapted its format to include a full-color print edition, and adjusted its aims and scope fit the fast-changing field of molecular modeling, with a particular focus on three-dimensional modeling.
Today, the journal covers all aspects of molecular modeling including life science modeling; materials modeling; new methods; and computational chemistry.
Topics include computer-aided molecular design; rational drug design, de novo ligand design, receptor modeling and docking; cheminformatics, data analysis, visualization and mining; computational medicinal chemistry; homology modeling; simulation of peptides, DNA and other biopolymers; quantitative structure-activity relationships (QSAR) and ADME-modeling; modeling of biological reaction mechanisms; and combined experimental and computational studies in which calculations play a major role.