The conformational landscape of three FK506-related drugs with disparate inhibition constants is determined in bulk solution using a replica exchange simulation method with solute torsional tempering. Energetic fitness of important drug conformations with respect to the FKBP12 protein is evaluated by molecular mechanics. Results show that the experimental affinity toward peptidyl-prolyl cis-trans isomerase protein (FKBP12) of the analyzed ligands appears to be positively correlated to the observed population of specific chair structures of the drug piperidinic ring in bulk solution. This observation is rationalized on the basis that such structures, stabilized by stereospecific intramolecular hydrophobic interactions, allows the formation of a pair of protein-ligand hydrogen bonds upon binding.
© 2011 American Chemical Society