Resistance without a detectable bindingenergy penalty: distinct dynamic defects of gatekeeper and molecular-brake mutations in FGFR1–erdafitinib recognition
- 1School of Computer Information
- 2Inner Mongolia Medical University
Abstract
The backbone Cα–Cα distance widely used to report kinase molecular-brake status gives a false negative under a clinically relevant FGFR1 resistance mutation: the brake's side-chain contact is destabilized in every replicate, yet the criterion reports it intact. Acquired FGFR1 kinase-domain mutations limit the clinical durability of erdafitinib, the first approved pan-FGFR inhibitor; loss of binding affinity is the assumed mechanism. We tested two questions in a single design: nine all-atom MD trajectories of the erdafitinib-bound kinase domain (wild type, V561M and molecular-brake mutant N546K; three replicates each; 4.491 μs). First, what a fully replicated comparison can resolve: neither MM/GBSA nor MM/PBSA detected a mutation-associated change in binding free energy; the observed mutant-versus-WT ΔΔG values lie below the detection limit of the present replicate set and within the intrinsic uncertainty of the end-point methods (GB ANOVA p = 0.27; PB p = 0.054). A large, E3810-like penalty (≈5 kcal·mol-1) is excluded. Second, the brake itself—N546K destabilizes the persistent engaged state of the residue-546–Glu562 contact: occupancy falls to 16–55% across replicates, entirely below the 95–100% of WT and V561M, at unchanged Cα–Cα geometry, dwell times fall from tens of nanoseconds to sub-nanosecond flickers, and broken periods reach 19 ns. V561M leaves the brake contact intact and instead loosens the gatekeeper-flanking hydrophobic network (Leu484/Val492/Lys514), dispersing its conformational ensemble across replicates. Side-chain proximity occupancy thus emerges as an operational order parameter for the molecular brake; the kinetic corollary, altered rates at unchanged affinity, is testable by SPR/BLI and unbinding calculations.
Related
- pH-Responsive Erdafitinib Liposomes with FGFR3 Targeting for Intravesical Bladder Cancer Therapy
- Single Optically Detectable Tumbling Spin in Silicon
- Surface-localized thermoresponsive molecularly imprinted polymer for gentamicin recognition and aqueous regeneration
- Conserved recognition and divergent peripheral dynamics of designed statin analogues at human HMG-CoA reductase: ensemble docking and replicated molecular dynamics
- Molecular recognition of Atrazine, 2,4-Dichlorophenoxyacetic acid, and Chlorpyrifos by cyclodextrin-based adsorbents: insights into selective pesticide removal
BibTeX
@misc{he2026resistance,
title = {{Resistance without a detectable bindingenergy penalty: distinct dynamic defects of gatekeeper and molecular-brake mutations in FGFR1–erdafitinib recognition}},
author = {Jia He},
howpublished = {ChemRxiv},
year = {2026},
month = oct,
doi = {10.26434/chemrxiv.15011002/v1},
url = {https://chemrxiv.org/doi/full/10.26434/chemrxiv.15011002/v1},
note = {CC BY},
}
RIS
TY - UNPB
TI - Resistance without a detectable bindingenergy penalty: distinct dynamic defects of gatekeeper and molecular-brake mutations in FGFR1–erdafitinib recognition
AU - Jia He
PY - 2026
DA - 2026/10/25
PB - ChemRxiv
DO - 10.26434/chemrxiv.15011002/v1
UR - https://chemrxiv.org/doi/full/10.26434/chemrxiv.15011002/v1
AB - The backbone Cα–Cα distance widely used to report kinase molecular-brake status gives a false negative under a clinically relevant FGFR1 resistance mutation: the brake's side-chain contact is destabilized in every replicate, yet the criterion reports it intact. Acquired FGFR1 kinase-domain mutations limit the clinical durability of erdafitinib, the first approved pan-FGFR inhibitor; loss of binding affinity is the assumed mechanism. We tested two questions in a single design: nine all-atom MD trajectories of the erdafitinib-bound kinase domain (wild type, V561M and molecular-brake mutant N546K; three replicates each; 4.491 μs). First, what a fully replicated comparison can resolve: neither MM/GBSA nor MM/PBSA detected a mutation-associated change in binding free energy; the observed mutant-versus-WT ΔΔG values lie below the detection limit of the present replicate set and within the intrinsic uncertainty of the end-point methods (GB ANOVA p = 0.27; PB p = 0.054). A large, E3810-like penalty (≈5 kcal·mol-1) is excluded. Second, the brake itself—N546K destabilizes the persistent engaged state of the residue-546–Glu562 contact: occupancy falls to 16–55% across replicates, entirely below the 95–100% of WT and V561M, at unchanged Cα–Cα geometry, dwell times fall from tens of nanoseconds to sub-nanosecond flickers, and broken periods reach 19 ns. V561M leaves the brake contact intact and instead loosens the gatekeeper-flanking hydrophobic network (Leu484/Val492/Lys514), dispersing its conformational ensemble across replicates. Side-chain proximity occupancy thus emerges as an operational order parameter for the molecular brake; the kinetic corollary, altered rates at unchanged affinity, is testable by SPR/BLI and unbinding calculations.
ER -