Affinity and Selectivity Assessment of Covalent Inhibitors by Free Energy Calculations

Covalent inhibitors have been gaining increased attention in drug discovery due to their beneficial properties such as long residence time, high biochemical efficiency, and specificity. Optimization of covalent inhibitors is a complex task that involves parallel monitoring of the noncovalent recognition elements and the covalent reactivity of the molecules to avoid potential idiosyncratic side effects. This challenge calls for special design protocols, including a variety of computational chemistry methods. Covalent inhibition proceeds through multiple steps, and calculating free energy changes of the subsequent binding events along the overall binding process would help us to better control the design of drug candidates. Inspired by the recent success of free energy calculations on reversible binders, we developed a complex protocol to compute free energies related to the noncovalent and covalent binding steps with thermodynamic integration and hybrid quantum mechanical/molecular mechanical (QM/MM) potential of mean force (PMF) calculations, respectively. In optimization settings, we examined two therapeutically relevant proteins complexed with congeneric sets of irreversible cysteine targeting covalent inhibitors. In the selectivity paradigm, we studied the irreversible binding of covalent inhibitors to phylogenetically close targets by a mutational approach. The results of the calculations are in good agreement with the experimental free energy values derived from the inhibition and kinetic constants (K_i and k_inact) of the enzyme-inhibitor binding. The proposed method might be a powerful tool to predict the potency, selectivity, and binding mechanism of irreversible covalent inhibitors.

One- and Two-Proton Transfer Mechanisms Coexist in One Active Site

Acibenzolar-S-methyl (ASM) is one of the most successfully commercialized plant activators of the systemic acquired resistance (SAR). However, its activation (hydrolysis) mechanism catalyzed by the salicylic acid binding protein 2 (SABP2) remains elusive. The fundamental catalytic mechanism of the SABP2-catalyzed hydrolysis of the ASM had been investigated by extensive computational and experimental studies, including QM/MM simulations, charge transfer analysis, small-molecule synthesis, and biochemical assays. Here we report that the promiscuous SABP2 shows different catalytic mechanisms toward different substrates. To catalyze the ASM hydrolysis, the SABP2 uses a two-proton transfer mechanism, and the key intermediate is stabilized by the charge transfer effect; to catalyze the ethyl 1,2,3-benzothiadiazole-7-carboxylate (BTM, an ASM analogue) hydrolysis, the SABP2 applies the one-proton transfer mechanism, and the classic tetrahedral intermediate is stabilized by the electrostatic effect. The HPLC analyses of the SABP2 esterase activities toward the ASM and the BTM show comparable results with our computaional results, suggesting that the obtained computational mechanism insights are reasonable. The obtained mechanism is not only an important supplement to the theory of enzymes' catalytic promiscuity, but it also contributes a possible strategy for the design of next generation plant SAR activators.