Proof of concept for poor inhibitor binding and efficient formation of covalent adducts of KRAS~(G12C) and ARS compounds

摘要

The use of selective covalent inhibitors with low binding affinity and high reactivity with the target enzyme is a promising way to solve a long-standing problem of the "undruggable" RAS-like proteins. Specifically, compounds of the ARS family that prevent the activation of the GDP-bound G12C mutant of Kirsten RAS (KRAS) are in the focus of recent experimental research. We report the first computational characterization of the entire reaction mechanism of the covalent binding of ARS-853 to the KRAS~(G12C). GDP complex. The application of molecular dynamics, molecular docking and quantum mechanics/molecular mechanics approaches allowed us to model the inhibitor binding to the protein and the chemical reaction of ARS-853 with Cysl2 in the enzyme binding site. We estimated a full set of kinetic constants and carried out numerical kinetic analysis of the process. Thus, we were able to compare directly the physicochemical parameters of the reaction obtained in silico and the macroscopic parameters observed in experimental studies. From our computational results, we explain the observed unusual dependence of the rate constant of covalent complex formation, k_(obs), on the ARS concentration. The latter depends both on the non-covalent binding step with the equilibrium constant, K,, and on the rate constant of covalent adduct formation, k_(inact). The calculated ratio k_(inact)/K_i = 213 M~(-1)s~(-1) reproduces the corresponding experimental value of 250 ± 30 M~(1) s~(-1) for the interaction of ARS-853 with KRAS~(G12C) Electron density analysis in the reactive region demonstrates that covalent bond formation occurs efficiently according to the Michael addition mechanism, which assumes the activation of the C=C bond of ARS-853 by a water molecule and Lysl6 in the binding site of KRAS~(G12C). We also refine the k_(inact) and K_i constants of the ARS-107 compound, which shares common features with ARS-853, and show that the decrease in the k_(inact)/K_i ratio in the case of ARS-107 is explained by changes in both K_i and k_(inact) constants.