Enhanced sampling and ligandability assessment to expand the repertoire of potentially druggable cryptic pockets

Certain proteins known to be involved in life-threatening diseases remained challenging therapeutic targets for decades simply because a suitable binding pocket for a potent molecular inhibitor could not be identified in their ground state apo structures. In the last decade, discovery of cryptic pockets in challenging targets like KRAS and Werner helicase has proven to be a major turning point for therapeutic development. However, the alternate protein conformations required for these cryptic pockets to exist were only revealed by experiments conducted in the presence of ligands that can bind to them. Time-consuming and expensive experiments currently used to uncover these biologically rare events could be usefully complemented by a computational method capable of reliably identifying cryptic pockets. We have previously shown that aqueous and mixed-solvent Weighted Ensemble molecular dynamics (WEMD) simulations driven by normal modes representing the direction of the most-collective motion of a protein can predict known cryptic pockets in the KRAS oncoprotein. Here, we evaluate this cryptic pocket detection technique on a dataset of diverse proteins and show that it successfully identifies cryptic pockets over 92% of the time starting with just the apo structure. We also show that we can successfully rank candidate pockets from WEMD using our pocket ligandability prediction model, Target X.