Abstract
Membrane proteins represent a considerable fraction of pharmaceutical drug targets. A computational technique to identify ligand binding pockets in these proteins is therefore of great importance. We recently reported such a technique called pMD-membrane that utilizes small molecule probes to detect ligand binding sites and surface hotspots on membrane proteins based on probe-based molecular dynamics simulation. The current work extends pMD-membrane to a diverse set of small organic molecular species that can be used as cosolvents during simulation of membrane proteins. We also describe a projection technique for globally quantifying probe densities on the protein surface and introduce a technique to construct surface topography maps directly from the probe-binding propensity of surface residues. The map reveals surface patterns and geometric features that aid in filtering out high probe density hotspots lacking pocketlike characteristics. We demonstrate the applicability of the extended pMD-membrane and the new analysis tool by exploring the druggability of full-length G12D, G12V, and G13D oncogenic K-Ras mutants bound to a negatively charged lipid bilayer. Using data from 30 pMD-membrane runs conducted in the presence of a 2.8 M cosolvent made up of an equal proportion of seven small organic molecules, we show that our approach robustly identifies known allosteric ligand binding sites and other reactive regions on K-Ras. Our results also show that accessibility of some pockets is modulated by differential membrane interactions.