Solvent & geometric effects on non-covalent interactions

Transcription

1 Solvent & geometric effects on non-covalent interactions Scott L. Cockroft PhysChem Forum 10, Syngenta, Jealott s Hill, 23 rd March 11

2 QSAR & Physical Organic Chemistry Quantifiable Physicochemical Properties Non-covalent interactions (α, β, MEP) Electronic descriptors (σ m, σ p, σ*, F, R) Solubility Lipophilicity (logp/d, π) Ionisable groups (pk a ) Shape Molecular weight Administration QSAR Quantifiable Biological Properties Interactions with target Absorption Distribution Metabolism Excretion Toxicology output: e.g. identify structural features giving good potency, identify outliers, predict new structures with better biological profiles

3 Why study non-covalent interactions and solvation? DRUG Drug-solvent interactions Drug-lipid interactions Drug-solvent interactions Drug-receptor interactions (octanol?) + solvent-solvent interactions (liberated water) + Proteinsolvent interactions + solvent-solvent interactions (liberated water)

4 Non-covalent interactions in biological recognition COMPLICATED! P. A. Williams, J. Cosme et al, Nature, 2003, 424,

5 Total binding energy contributions interaction changes structural changes Binding Energy, G = H T S electrostatic rotational induction translational repulsion vibrational van der Waals interactions dispersion desolvation Energy contributions from the WHOLE system = drug, solvent, receptor

6 H-bonds: back to basics, experiment vs. theory H-bond interaction energies are highly correlated with the maxima and minima of calculated molecular electrostatic potential surfaces. Calculated elec ctrostatic potential α s β MEP MAX MEP MIN α 2 H (H-bond donors) β 2 H (H-bond acceptors) Experimental H-bond parameters measured in CCl 4 M. H. Abraham, J. A. Platts, J. Org. Chem. 2001, 66, Hunter & Cockroft unpublished results

7 H-bond scales Possible to rescale α 2Η and β 2Η H-bond constants such that: + α = H-bond donor constant β = H-bond acceptor G = αβ constant in kj mol -1 H-bond donor-acceptor interaction C. A. Hunter, Angew. Chem. Int. Ed. 2004, 43,

8 Interaction energies in kj mol-1

9 Total binding energy contributions interaction changes structural changes Binding Energy, G = H T S electrostatic rotational induction translational repulsion vibrational van der Waals interactions dispersion desolvation Energy contributions from the WHOLE system = drug, solvent, receptor

10 Solvent effects on non-covalent interactions K a α s β αβ s αβ α s β s G = α s β + αβ s αβ + α s β s + entropic term α/β obtained from experimental measurements or calculated electrostatic potentials C. A. Hunter, Angew. Chem. Int. Ed. 2004, 43,

11 Solvent effects on non-covalent interactions This simple model has been applied in a range of circumstances: (CF 3 ) 3 C OH O=P(n-Bu) 3 COM MPLEXITY Single H-bonds in pure solvents and solvent mixtures J. L. Cook, C. A. Hunter, C. M. R. Low, A. Perez-Velasco, J. Vinter, Angew. Chem. Int. Ed Single H-bonds in solvent mixtures R. Cabot, C. A. Hunter, Org. Biomol. Chem., 2010 Supramolecular complexes & Folding molecules S. L. Cockroft, C. A. Hunter. Chem. Commun & 2009

12 Solvent effects on H-bonding interactions (CF 3 ) 3 C OH O=P(n-Bu) 3 log K experiment log K predicted J. L. Cook, C. A. Hunter et al. Angew. Chem. Int. Ed. 2007, 46,

13 The puzzle of the Wilcox torsion balance The α/β scale works well for simple H-bond interactions one step up from a single H-bond, a simple supramolecular system. Is it really that simple? X G fold /kj mol -1 NO 2 CN I Br CH 3 OH NH X Me N N O O GK fold δ+ + δ X O O N N Me The experiments support a conclusion that the electrostatic potential of the aromatic ring is not a dominant aspect of the aryl-aryl interaction. The results should No obvious encourage trend increased with variation emphasis of X-substituents. on the importance of London dispersion forces in on-face aryl interactions involving neutral components >250 literature citations S. Paliwal, S. Geib, C. S. Wilcox, J. Am. Chem. Soc. 1994, 116, 4497 E. Kim, S. Paliwal, C. S. Wilcox, J. Am. Chem. Soc. 1998, 120, 11192

14 Aromatic rings as H-bond acceptors kj mol -1 β 0 kj mol face kj mol -1 (H-bond acceptor constant)

15 Benzene vs. Chloroform -5-4 in C 6 D 6 kj mol G fold / in CDCl 3 0 δ β face (H-bond acceptor constant) 2 3 δ - F. Hof, D. M. Scofield, W. B. Schweizer, F. Diederich, Angew. Chem. Int. Ed. 2004, 43, S. L. Cockroft, C. A. Hunter. Chem. Commun. 2006, 36,

16 G sol-face G sol-edge G edge-face G sol-sol α sol x β face constant α edge x β face constant S. L. Cockroft, C. A. Hunter. Chem. Commun. 2006, 36,

17 Aromatic solvation 1 x α sol 2 x α sol S. L. Cockroft, C. A. Hunter. Chem. Commun. 2006, 36,

18 -5-4 BENZENE CF CHLOROFORM G fold / kj mol H β face G fold / kj mol β face α sol x β face α edge x β face Chloroform 2.2 sl. > 1.0 x 2 Benzene 1.0 x 2 << 1.6 x 2 S. L. Cockroft, C. A. Hunter. Chem. Commun. 2006, 36,

19 -5-4 BENZENE CF CHLOROFORM G fold / kj mol H 0 G fold / kj mol CF 3 H β face β face α sol x β face α edge x β face Chloroform 2.2 << sl. > 2.0 Benzene 2.0 << = 2.0 Cockroft & Hunter, Chem. Comm. 2009

20 Solvent effects on non-covalent interactions This simple model has been applied in a range of circumstances: (CF 3 ) 3 C OH O=P(n-Bu) 3 COM MPLEXITY Single H-bonds in pure solvents and solvent mixtures J. L. Cook, C. A. Hunter, C. M. R. Low, A. Perez-Velasco, J. Vinter, Angew. Chem. Int. Ed Single H-bonds in solvent mixtures R. Cabot, C. A. Hunter, Org. Biomol. Chem., 2010 Supramolecular complexes & Folding molecules S. L. Cockroft, C. A. Hunter. Chem. Commun & 2009 BIOLOGY???

21 Range of aromatic interactions quantified Y = NMe 2, H, NO 2 X = NMe 2, t-bu, H, NO 2, F, I, CF 3 27 edge-face interactions (+2 to 7 kj mol -1 ) X analysis Y S Y Cockroft & Hunter, Chem. Comm. 2009

22 from electrostatic potentials to α edge & β face +100 kj mol -1 0 kj mol kj mol -1 N H Cockroft & Hunter, Chem. Comm. 2009

23 Electrostatic control of aromatic interactions Same model as used to analyse Diederich data! interaction energy / kj mol -1 Experimental i G X Y H N O H N O N H N H O O Y X G / kj mol -1 Predicted interaction energy Cockroft & Hunter, Chem. Comm. 2009

24 Non-covalent interactions: non-equilibrium geometry NON-EQUILIBRIUM GEOMETRY & BRIDGING SOLVENT MOLECULES! P. A. Williams, J. Cosme et al, Nature, 2003, 424,

25 Separation-interaction profiles G

26 Separation-interaction profiles Implications: understand interplay of solvent and geometry & conformation of small molecules determine importance of dispersion interactions-what happens when solvent excluded? quantify significance of long-range solvent-mediated (bridging) interactions predict solubility better force-fields for molecular modelling/docking, crystal structure prediction etc Many computational docking simulations use solvent fields because adding explicit solvent molecules is computationally expensive, is this a valid approach?

27 Solvent effects on molecular folding Folded Unfolded + K fold transition state G folded G G fold unfolded = RT lnk fold Folding coordinate I. K. Mati & S. L. Cockroft. Chem. Soc. Rev. 2010

28 Foldamers for quantifying non-covalent interactions Y Z K fold Y Z Wilcox, Diederich Motherwell Oki, Gung Gellman I. K. Mati & S. L. Cockroft. Chem. Soc. Rev. 2010

29 Non-covalent interactions summary α β K a X O H F N X K fold dispersion interactions non-equilibrium geometries solvent effects

30 March 2010 Sheep Heid, Edinburgh (Oldest pub in Scotland)

31