The reorganization energy of compounds upon binding to proteins, from dynamic and solvated bound and unbound states. (1st December 2021)

• Record Type: Journal Article
• Authors: Foloppe, Nicolas; Chen, I-Jen
• Is Part Of: Bioorganic & medicinal chemistry. Volume 51(2021)
• Publication Date: 2021-12-01
• Subjects: Computational chemistry -- Drug discovery -- Electrostatics -- Molecular dynamics -- Molecular recognition -- Reorganization energy -- Simulation; DHFR Dihydrofolate Reductase -- FF Force Field -- GB Generalized Born solvation model -- GPU Graphical Processing Unit -- HCV hepatitis C virus -- HIV human immunodeficiency virus -- MD Molecular Dynamics -- MM Molecular Mechanics -- MOE Molecular Operating Environment -- MT/LMOD Mixed torsional/Low-mode -- MW molecular weight -- NRot Oprea Number of Rotatable bonds -- OPLS Optimized Potential for Liquid Simulations force-field -- PDB Protein Data Bank -- Rgyr Radius of Gyration -- SBDD Structure-Based Drug Design -- SD Standard Deviation -- TPSA Topological Polar Surface Area; Bioorganic chemistry -- Periodicals; Pharmaceutical chemistry -- Periodicals; Biochemistry -- Periodicals; Chemistry, Clinical -- Periodicals; Chemistry, Organic -- Periodicals; Chimie bio-organique -- Périodiques; Chimie pharmaceutique -- Périodiques; 615.19
• DOI: 10.1016/j.bmc.2021.116464
• Languages: English
• ISSNs: 0968-0896

Abstract

Graphical abstract: Abstract: The intramolecular reorganization energy (ΔEReorg ) of compounds upon binding to proteins is a component of the binding free energy, which has long received particular attention, for fundamental and practical reasons. Understanding ΔEReorg would benefit the science of molecular recognition and drug design. For instance, the tolerable strain energy of compounds upon binding has been elusive. Prior studies found some large ΔEReorg values (e.g. > 10 kcal/mol), received with skepticism since they imply excessive opposition to binding. Indeed, estimating ΔEReorg is technically difficult. Typically, ΔEReorg has been approached by taking two energy-minimized conformers representing the bound and unbound states, and subtracting their conformational energy. This is a drastic oversimplification, liable to conformational collapse of the unbound conformer. Instead, the present work applies extensive molecular dynamics (MD) and the modern OPLS3 force-field to simulate compounds bound and unbound states, in explicit solvent under physically relevant conditions. The thermalized unbound compounds populate multiple conformations, not reducible to one or a few energy-minimized conformers. The intramolecular energies in the bound and unbound states were averaged over pertinent conformational ensembles, and the reorganization enthalpy upon binding (ΔHReorg ) deduced by subtraction. This was applied to 76 systems, including 43 approved drugs, carefully selected for i) the quality of the bioactive X-ray structures and ii) the diversity of the chemotypes, their properties and protein targets. It yielded comparatively low ΔHReorg values (median = 1.4 kcal/mol, mean = 3.0 kcal/mol). A new finding is the observation of negative ΔHReorg values. Indeed, reorganization energies do not have to oppose binding, e.g. when intramolecular interactions stabilize preferentially the bound state. Conversely, even with competing water molecules, intramolecular interactions can occur predominantly in the unbound compound, and be replaced by intermolecular counterparts upon protein binding. Such disruption of intramolecular interactions upon binding gives rise to occasional larger ΔHReorg values. Such counterintuitive larger ΔHReorg values may be rationalized as a redistribution of interactions upon binding, qualitatively compatible with binding.