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

详细信息    查看全文

• 作者：Nicolas Foloppe ; ^{n.foloppe@vernalis.com" class="auth_mail" title="E-mail the corresponding author} ; I-Jen Chen ; ^{i.chen@vernalis.com" class="auth_mail" title="E-mail the corresponding author}
• 关键词：ACE ; angiotensin converting enzyme ; DHFR ; dihydrofolate reductase ; Diel ; distance-dependent dielectric ; FF ; force field ; GB ; generalized Born solvation model ; GlobMin_Diel ; global energy minimum obtained with a distant-dependent dielectric constant ; GlobMin_GB ; global energy minimum obtained with a GB solvation model ; GlobMin_Vac ; global energy minimum obtained in vacuo ; GPU ; Graphical Processing Unit ; GUI ; graphical user interface ; LowModeMD ; search method combining low-mode moves and molecular
• 刊名：Bioorganic & Medicinal Chemistry
• 出版年：2016
• 出版时间：15 May 2016
• 年：2016
• 卷：24
• 期：10
• 页码：2159-2189
• 全文大小：8965 K

文摘

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 (ΔE_Reorg). Also, dependable observation of the unbound compounds is important for ligand-based drug discovery, e.g. with pharmacophore modelling. Here, these questions are addressed with long (⩾0.5 μ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, MD-simulated 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 ΔE_Reorg as an enthalpic difference ΔH_Reorg between the intramolecular energies in the bound and unbound states. This fresh approach yielded ΔH_Reorg values ⩽ 6 kcal/mol for 18 out of 26 compounds. For three particularly polar compounds 15 ⩽ ΔH_Reorg ⩽ 20 kcal/mol, supporting the notion that ΔH_Reorg can be substantial. Those large Δ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.