Targeted covalent inhibitor design is gaining increasing interest and acceptance. A typical covalent kinase inhibitor design targets a reactive cysteine; however, this strategy is limited due to the low abundance of cysteine and acquired drug resistance from point mutations. Inspired by the recent development of lysine-targeted chemical probes, we asked if nucleophilic (reactive) catalytic lysines are common based on the published crystal structures of the human kinome. Using a newly developed pK a prediction tool based on continuous constant pH molecular dynamics, the catalytic lysines of 8 unique kinases from various human kinase groups were retro-and prospectively predicted to be nucleophilic, when kinase is in the rare DFG-out/αC-out type of conformation. Importantly, other reactive lysines as well as cysteines at various locations were also identified. Based on the finding, we proposed a new strategy in which selective type II reversible kinase inhibitors are modified to design highly selective, lysine-targeted covalent inhibitors. Traditional covalent drugs were discovered serendipitously; the presented tool, which can assess the reactivities of any potentially targetable residues, may accelerate the rational discovery of new covalent inhibitors. Another significant finding of the work is that lysines and cysteines in kinases may adopt neutral and charged states at physiological pH, respectively. This finding may shift the current paradigm of computational studies of kinases, which assume standard protonation states.

Recent development of titratable coions has paved the way for realizing all-atom molecular dynamics at constant pH. To further improve physical realism, here we describe a technique in which proton titration of the solute is directly coupled to the interconversion between water and hydroxide or hydronium. We test the new method in replica-exchange continuous constant pH molecular dynamics simulations of three proteins, HP36, BBL, and HEWL. The calculated pKa values based on 10-ns sampling per replica have the average absolute and root-mean-square errors of 0.7 and 0.9 pH units, respectively. Introducing titratable water in molecular dynamics offers a means to model proton exchange between solute and solvent, thus opening a door to gaining new insights into the intricate details of biological phenomena involving proton translocation.

Despite the relevance of understanding structure-function relationships, robust prediction of proton donors and nucleophiles in enzyme active sites remains challenging. Here we tested three types of state-of-the-art computational methods to calculate the p K's of the buried and hydrogen bonded catalytic dyads in five enzymes. We asked the question what determines the p K order, i.e., what makes a residue proton donor vs a nucleophile. The continuous constant pH molecular dynamics simulations captured the experimental p K orders and revealed that the negative nucleophile is stabilized by increased hydrogen bonding and solvent exposure as compared to the proton donor. Surprisingly, this simple trend is not apparent from crystal structures and the static structure-based calculations. While the generality of the findings awaits further testing via a larger set of data, they underscore the role of dynamics in bridging enzyme structures and functions.

Permeability (P m) across biological membranes is of fundamental importance and a key factor in drug absorption, distribution and development. Although the majority of drugs will be charged at some point during oral delivery, our understanding of membrane permeation by charged species is limited. The canonical model assumes that only neutral molecules partition into and passively permeate across membranes, but there is mounting evidence that these processes are also facile for certain charged species. However, it is unknown whether such ionizable permeants dynamically neutralize at the membrane surface or permeate in their charged form. To probe protonationcoupled permeation in atomic detail, we herein apply continuous constant-pH molecular dynamics along with free energy sampling to study the permeation of a weak base propranolol (PPL), and evaluate the impact of including dynamic protonation on P m. The simulations reveal that PPL dynamically neutralizes at the lipid-tail interface, which dramatically influences the permeation free energy landscape and explains why the conventional model overestimates the assigned intrinsic permeability. We demonstrate how fixed-charge-state simulations can account for this effect, and propose a revised model that better describes pH-coupled partitioning and permeation. Our results demonstrate how dynamic changes in protonation state may play a critical role in the permeation of ionizable molecules, including pharmaceuticals and drug-like molecules, thus requiring a revision of the standard picture.