Abstract

The binding affinity of a ligand to a protein is influenced by the protonation and tautomeric states of both partners. However, this relationship remains under-investigated due to the limited availability of computational tools capable of considering all charge and tautomer states in a scalable manner to study clinically relevant systems. Here, we use Multi-Conformation Continuum Electrostatics (MCCE) to calculate the protonation and tautomer distributions of nine kinase domains bound to 18 FDA-approved inhibitors while considering their Boltzmann-ensemble. Our simulations show that protein net charge and proton distribution remain largely stable even upon binding charged inhibitors. Our results find that individual inhibitor charges are dynamic, frequently increasing, or decreasing upon binding a specific protein target. Kinase-inhibitor binding significantly shifts the relative probabilities of low-energy states ({Delta}G < 2.5 kcal/mol), though it does not recruit higher-energy conformers into the bound population. Our consideration of all possible charge states and tautomers enable us to identify when tautomer have significant significant free binding energy differentials (3-6 kcal/mol). In turn, we find that minority species can become the dominant component in the bound state, emphasizing the necessity of considering ensemble-wide protonation and tautomer states to accurately predict protein-ligand binding energetics.