Investigating the ibrutinib resistance mechanism of L528W mutation on Bruton's tyrosine kinase via molecular dynamics simulations

J Mol Graph Model. 2024 Jan:126:108623. doi: 10.1016/j.jmgm.2023.108623. Epub 2023 Sep 12.

Abstract

Drug resistance to Bruton's Tyrosine Kinase (BTK) inhibitors presents a challenge in treating B-cell malignancies, and the mechanism behind drug resistance remains unclear. In this study, we focused on the BTK L528W mutation and investigated the underlying mechanisms of resistance to ibrutinib (including prototype and its active metabolite from, PCI-45227) using a combination of bioinformatics analysis, and molecular dynamics (MD) simulations. Protein stability of wild type (WT) BTK and L528W mutant was predicted using DUET, PoPMuSiC, and I-Mutant2.0. We performed MD simulations of six systems, apo-WT, metabolite-WT, prototype-WT and their mutants, to analyze the significant conformational and BTK-inhibitor binding affinity changes induced by the L528W mutation. Results show that the L528W mutation reduces the conformational stability of BTK compared to the WT. Principal component analysis (PCA) based free energy landscape (FEL) analysis shows that the L528W mutant ensemble tends to form more conformation clusters and exhibit higher levels of local minima than the WT counterpart. The interaction analysis reveal that the L528W mutation disrupts the strong hydrogen bond between Cys481 and inhibitors and reduces the number of hydrogen bonds between inhibitors and BTK in the L528W mutant complex structures compared to the WT. Porcupine plot analysis in association with cross-correlation analysis show the high-intensity flexible motion exhibited by the P-loop region. MM/GBSA calculations show that the L528W mutation in metabolite-BTK and prototype-BTK complexes increases binding free energy compared to the WT, with a reduction in binding affinity confirmed by per-residue energy decomposition. Specifically, the binding free energy increases from -57.86 kcal/mol to -48.26 kcal/mol for the metabolite-BTK complex and from -62.04 kcal/mol to -50.55 kcal/mol for the prototype-BTK complex. Overall, our study finds that the L528W mutation reduces BTK stability, decreases binding affinity, and leads to drug resistance and potential disease recurrence.

Keywords: Bruton's tyrosine kinase; Ibrutinib; Molecular dynamics simulation; Mutation; Resistance mechanism.

Publication types

  • Research Support, Non-U.S. Gov't

MeSH terms

  • Agammaglobulinaemia Tyrosine Kinase / genetics
  • Drug Resistance, Neoplasm* / genetics
  • Molecular Dynamics Simulation*
  • Mutation
  • Tyrosine Kinase Inhibitors* / pharmacology

Substances

  • Agammaglobulinaemia Tyrosine Kinase
  • ibrutinib
  • Tyrosine Kinase Inhibitors