Towards understanding the unbound state of drug compounds: Implications for the intramolecular reorganization energy upon binding

There has been an explosion of structural information for pharmaceutical compounds bound to biological targets, but the conformations and dynamics of compounds free in solution are poorly characterized, if at all. Yet, knowledge of the unbound state is essential to understand the fundamentals of molecular recognition, including the much debated conformational intramolecular reorganization energy of a compound upon binding (Delta E-Reorg). Also, dependable observation of the unbound compounds is important for ligandbased drug discovery, e. g. with pharmacophore modelling. Here, these questions are addressed with long (>= 0.5 mu s) state-of-the-art molecular dynamics (MD) simulations of 26 compounds (including 7 approved drugs) unbound in explicit solvent. These compounds were selected to be chemically diverse, with a range of flexibility, and good quality bioactive X-ray structures. The MD-simulated free compounds are compared to their bioactive structure and conformers generated with ad hoc sampling in vacuo or with implicit generalized Born (GB) aqueous solvation models. The GB conformational models clearly depart from those obtained in explicit solvent, and suffer from conformational collapse almost as severe as in vacuo. Thus, the global energy minima in vacuo or with GB are not suitable representations of the unbound state, which can instead be extensively sampled by MD simulations. Many, but not all, MDsimulated compounds displayed some structural similarity to their bioactive structure, supporting the notion of conformational pre-organization for binding. The ligand-protein complexes were also simulated in explicit solvent, to estimate Delta E-Reorg as an enthalpic difference Delta H-Reorg between the intramolecular energies in the bound and unbound states. This fresh approach yielded Delta H-Reorg values <= 6 kcal/mol for 18 out of 26 compounds. For three particularly polar compounds 15 <= Delta H-Reorg <= 20 kcal/mol, supporting the notion that Delta H-Reorg can be substantial. Those large Delta H-Reorg values correspond to a redistribution of electrostatic interactions upon binding. Overall, the study illustrates how MD simulations offer a promising avenue to characterize the unbound state of medicinal compounds. (C) 2016 Elsevier Ltd. All rights reserved.