Optimizing the Cardiotoxicity Profile of Osimertinib via Scaffold Hopping: An Integrated Cheminformatics and Molecular Modeling Approach

Osimertinib, a third-generation covalent epidermal growth factor receptor (EGFR) inhibitor widely used to treat EGFR T790M-positive non-small cell lung cancer (NSCLC), is limited by dose-dependent cardiotoxicity driven by human ether-à-go-go-related gene (hERG) channel inhibition. To design safer analogues, we applied an integrated computational workflow combining ADMET filtering, covalent docking, molecular dynamics, and quantum-chemical analysis. Using ADMETopt and admetSAR 2.0, a scaffold-hopping campaign generated 390 analogs, of these 161 compounds exhibited substantially reduced predicted hERG inhibition (<0.70) and were prioritized for structure-based studies. Covalent docking against double mutant L858R/T790M EGFR TK (PDB ID: 5HCX) using Schrödinger CovDock identified 101 analogs with covalent docking scores below −5.0 kcal ⋅ mol[Formula: see text], comparable to or better than Osimertinib (−5.79 kcal ⋅ mol[Formula: see text] and retaining productive orientation toward Cys797. Top analogue compound 252 demonstrated favorable covalent docking (−9.26 kcal ⋅ mol[Formula: see text] and was subjected to a 100-ns molecular dynamics simulation. The complex exhibited lower backbone RMSD, reduced protein fluctuations, and persistent active-site interactions relative to the Osimertinib complex. PCA, DCCM, and FEL analyses confirmed stabilization of EGFR in a single dominant conformational basin, indicative of enhanced dynamic stability. DFT calculations at the B3LYP/6-311G++ ⋆⋆ level showed compound 252 possessed a favorable HOMO–LUMO gap (0.1416 eV), lower LUMO energy, and higher electrophilicity ([Formula: see text] = 0.115 eV), supporting efficient Michael acceptor reactivity. Overall, compound 252 represents a quantitatively validated cardiotoxicity-optimized EGFR inhibitor, demonstrating improved stability, reactivity, and predicted safety. This study establishes a robust computational framework for designing next-generation covalent kinase inhibitors with reduced off-target cardiotoxicity.

Authors

Institutions

Publication Details

Journal
Journal of Computational Biophysics and Chemistry
Published
2026-04-10
DOI
https://doi.org/10.1142/s2737416526500754
Primary Topic
Lung Cancer Treatments and Mutations
Type
article
Field-Weighted Citation Impact
0.00
Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

Optimizing the Cardiotoxicity Profile of Osimertinib via Scaffold Hopping: An Integrated Cheminformatics and Molecular Modeling Approach

Harun Patel, Bhatu Patil, Narendra R. Nagpure, Vaibhav Kadam et al.
Journal of Computational Biophysics and Chemistry
Lung Cancer Treatments and Mutations
article

Optimizing the Cardiotoxicity Profile of Osimertinib via Scaffold Hopping: An Integrated Cheminformatics and Molecular Modeling Approach

Harun Patel, Bhatu Patil, Narendra R. Nagpure, Vaibhav Kadam, Chandrakant Gawli, Sampath Chinnam, Iqrar Ahmad
article en

Abstract

Osimertinib, a third-generation covalent epidermal growth factor receptor (EGFR) inhibitor widely used to treat EGFR T790M-positive non-small cell lung cancer (NSCLC), is limited by dose-dependent cardiotoxicity driven by human ether-à-go-go-related gene (hERG) channel inhibition. To design safer analogues, we applied an integrated computational workflow combining ADMET filtering, covalent docking, molecular dynamics, and quantum-chemical analysis. Using ADMETopt and admetSAR 2.0, a scaffold-hopping campaign generated 390 analogs, of these 161 compounds exhibited substantially reduced predicted hERG inhibition (<0.70) and were prioritized for structure-based studies. Covalent docking against double mutant L858R/T790M EGFR TK (PDB ID: 5HCX) using Schrödinger CovDock identified 101 analogs with covalent docking scores below −5.0 kcal ⋅ mol[Formula: see text], comparable to or better than Osimertinib (−5.79 kcal ⋅ mol[Formula: see text] and retaining productive orientation toward Cys797. Top analogue compound 252 demonstrated favorable covalent docking (−9.26 kcal ⋅ mol[Formula: see text] and was subjected to a 100-ns molecular dynamics simulation. The complex exhibited lower backbone RMSD, reduced protein fluctuations, and persistent active-site interactions relative to the Osimertinib complex. PCA, DCCM, and FEL analyses confirmed stabilization of EGFR in a single dominant conformational basin, indicative of enhanced dynamic stability. DFT calculations at the B3LYP/6-311G++ ⋆⋆ level showed compound 252 possessed a favorable HOMO–LUMO gap (0.1416 eV), lower LUMO energy, and higher electrophilicity ([Formula: see text] = 0.115 eV), supporting efficient Michael acceptor reactivity. Overall, compound 252 represents a quantitatively validated cardiotoxicity-optimized EGFR inhibitor, demonstrating improved stability, reactivity, and predicted safety. This study establishes a robust computational framework for designing next-generation covalent kinase inhibitors with reduced off-target cardiotoxicity.

Journal of Computational Biophysics and Chemistry
Shri Vile Parle Kelavani Mandal (IN), Visvesvaraya Technological University (IN)
Openalex Percentile: Top 27%
Lung Cancer Treatments and Mutations
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.