Bruton's tyrosine kinase (BTK) is a key protein from the TEC family and is involved in B-cell lymphoma occurrence and development. Targeting BTK is therefore an effective strategy for B-cell lymphoma treatment. Since previous studies on BTK have been limited to structure-function analyses of static protein structures, the dynamics of conformational change of BTK upon inhibitor binding remain unclear. Here, molecular dynamics simulations were conducted to investigate the molecular mechanisms of association and dissociation of a reversible (ARQ531) and irreversible (ibrutinib) small-molecule inhibitor to/from BTK. The results indicated that the BTK kinase domain was found to be locked in an inactive state through local conformational changes in the DFG motif, and P-, A-, and gatekeeper loops. The binding of the inhibitors drove the outward rotation of the C-helix, resulting in the upfolded state of Trp395 and the formation of the salt bridge of Glu445-Arg544, which maintained the inactive conformation state. Met477 and Glu475 in the hinge region were found to be the key residues for inhibitor binding. These findings can be used to evaluate the inhibitory activity of the pharmacophore and applied to the design of effective BTK inhibitors. In addition, the drug resistance to the irreversible inhibitor Ibrutinib was mainly from the strong interaction of Cys481, which was evidenced by the mutational experiment, and further confirmed by the measurement of rupture force and rupture times from steered molecular dynamics simulation. Our results provide mechanistic insights into resistance against BTK-targeting drugs and the key interaction sites for the development of high-quality BTK inhibitors. The steered dynamics simulation also offers a means to rapidly assess the binding capacity of newly designed inhibitors.
Keywords: ARQ531; Bruton’s tyrosine kinase; ibrutinib; molecular simulation; small molecule inhibitors.