QM/MM Simulations of Afatinib-EGFR Addition: The Role of β-Dimethylaminomethyl Substitution

Acrylamides are the most commonly used warheads of targeted covalent inhibitors (TCIs) directed at cysteines; however, the reaction mechanisms of acrylamides in proteins remain controversial, particularly for those involving protonated or unreactive cysteines. Using the combined semiempirical quantum mechanics (QM)/molecular mechanics (MM) free energy simulations, we investigated the reaction between afatinib, the first TCI drug for cancer treatment, and Cys797 in the EGFR kinase. Afatinib contains a β-dimethylaminomethyl (β-DMAM) substitution which has been shown to enhance the intrinsic reactivity and potency against EGFR for related inhibitors. Two hypothesized reaction mechanisms were tested. Our data suggest that Cys797 becomes deprotonated in the presence of afatinib, and the reaction proceeds via a classical Michael addition mechanism, with Asp800 stabilizing the ion-pair reactant state β-DMAM⁺/C797^- and the transition state of the nucleophilic attack. Our work elucidates an important structure-activity relationship of acrylamides in proteins.