The Free Energy Landscape of Small Molecule Unbinding

The spontaneous dissociation of six small ligands from the active site of FKBP (the FK506 binding protein) is investigated by explicit water molecular dynamics simulations and network analysis. The ligands have between four (dimethylsulphoxide) and eleven (5-diethylamino-2-pentanone) non-hydrogen atoms, and an affinity for FKBP ranging from 20 to 0.2 mM. The conformations of the FKBP/ligand complex saved along multiple trajectories (50 runs at 310 K for each ligand) are grouped according to a set of intermolecular distances into nodes of a network, and the direct transitions between them are the links. The network analysis reveals that the bound state consists of several subbasins, i.e., binding modes characterized by distinct intermolecular hydrogen bonds and hydrophobic contacts. The dissociation kinetics show a simple (i.e., single-exponential) time dependence because the unbinding barrier is much higher than the barriers between subbasins in the bound state. The unbinding transition state is made up of heterogeneous positions and orientations of the ligand in the FKBP active site, which correspond to multiple pathways of dissociation. For the six small ligands of FKBP, the weaker the binding affinity the closer to the bound state (along the intermolecular distance) are the transition state structures, which is a new manifestation of Hammond behavior. Experimental approaches to the study of fragment binding to proteins have limitations in temporal and spatial resolution. Our network analysis of the unbinding simulations of small inhibitors from an enzyme paints a clear picture of the free energy landscape (both thermodynamics and kinetics) of ligand unbinding.

[1]  M. Congreve,et al.  Recent developments in fragment-based drug discovery. , 2008, Journal of medicinal chemistry.

[2]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[3]  A. Caflisch,et al.  Structure-based optimization of potent and selective inhibitors of the tyrosine kinase erythropoietin producing human hepatocellular carcinoma receptor B4 (EphB4). , 2009, Journal of medicinal chemistry.

[4]  Sergei V Krivov,et al.  One-dimensional free-energy profiles of complex systems: progress variables that preserve the barriers. , 2006, The journal of physical chemistry. B.

[5]  Amedeo Caflisch,et al.  How does a simplified-sequence protein fold? , 2009, Biophysical journal.

[6]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[7]  A. Caflisch,et al.  Kinetic analysis of molecular dynamics simulations reveals changes in the denatured state and switch of folding pathways upon single‐point mutation of a β‐sheet miniprotein , 2008, Proteins.

[8]  Jean-Claude Latombe,et al.  Stochastic roadmap simulation: an efficient representation and algorithm for analyzing molecular motion , 2002, RECOMB '02.

[9]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[10]  Alexander D. MacKerell,et al.  Computational Fragment-Based Binding Site Identification by Ligand Competitive Saturation , 2009, PLoS Comput. Biol..

[11]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[12]  E. Springman,et al.  A Biacore biosensor method for detailed kinetic binding analysis of small molecule inhibitors of p38alpha mitogen-activated protein kinase. , 2004, Analytical biochemistry.

[13]  Amedeo Caflisch,et al.  Free energy surfaces from single-distance information. , 2010, The journal of physical chemistry. B.

[14]  K. Dill,et al.  Predicting absolute ligand binding free energies to a simple model site. , 2007, Journal of molecular biology.

[15]  Amedeo Caflisch,et al.  Flaviviral protease inhibitors identified by fragment-based library docking into a structure generated by molecular dynamics. , 2009, Journal of medicinal chemistry.

[16]  A. Warshel,et al.  Calculations of electrostatic interactions in biological systems and in solutions , 1984, Quarterly Reviews of Biophysics.

[17]  J. Aqvist,et al.  A new method for predicting binding affinity in computer-aided drug design. , 1994, Protein engineering.

[18]  Jeremy C. Smith,et al.  Hierarchical analysis of conformational dynamics in biomolecules: transition networks of metastable states. , 2007, The Journal of chemical physics.

[19]  David Chandler,et al.  Statistical mechanics of isomerization dynamics in liquids and the transition state approximation , 1978 .

[20]  P. Hajduk,et al.  A decade of fragment-based drug design: strategic advances and lessons learned , 2007, Nature Reviews Drug Discovery.

[21]  A Caflisch,et al.  Slow folding of cross-linked alpha-helical peptides due to steric hindrance. , 2010, The journal of physical chemistry. B.

[22]  F. Rao,et al.  The protein folding network. , 2004, Journal of molecular biology.

[23]  G Otting,et al.  Solution structure and dynamics of PEC-60, a protein of the Kazal type inhibitor family, determined by nuclear magnetic resonance spectroscopy. , 1994, Journal of molecular biology.

[24]  Amedeo Caflisch,et al.  Estimation of protein folding probability from equilibrium simulations. , 2005, The Journal of chemical physics.

[25]  S. Kearsley.  On the orthogonal transformation used for structural comparisons , 1989 .

[26]  David L Mobley,et al.  Predicting ligand binding affinity with alchemical free energy methods in a polar model binding site. , 2009, Journal of molecular biology.

[27]  Amedeo Caflisch,et al.  α-Helix folding in the presence of structural constraints , 2008, Proceedings of the National Academy of Sciences.

[28]  W. L. Jorgensen.  The Many Roles of Computation in Drug Discovery , 2004, Science.

[29]  Amedeo Caflisch,et al.  Change of the unbinding mechanism upon a mutation: A molecular dynamics study of an antibody–hapten complex , 2005, Protein science : a publication of the Protein Society.

[30]  A. Fersht,et al.  Movement of the position of the transition state in protein folding. , 1995, Biochemistry.

[31]  Kyle A. Beauchamp,et al.  Molecular simulation of ab initio protein folding for a millisecond folder NTL9(1-39). , 2010, Journal of the American Chemical Society.

[32]  P Burkhard,et al.  X-ray structures of small ligand-FKBP complexes provide an estimate for hydrophobic interaction energies. , 2000, Journal of molecular biology.

[33]  A. Fersht,et al.  Toward the rational design of p53-stabilizing drugs: probing the surface of the oncogenic Y220C mutant. , 2010, Chemistry & biology.

[34]  Edward M. Reingold,et al.  Graph drawing by force‐directed placement , 1991, Softw. Pract. Exp..

[35]  Amedeo Caflisch,et al.  Bulky side chains and non-native salt bridges slow down the folding of a cross-linked helical peptide: a combined molecular dynamics and time-resolved infrared spectroscopy study. , 2009, The journal of physical chemistry. B.

[36]  W. L. Jorgensen,et al.  Binding affinities for sulfonamide inhibitors with human thrombin using Monte Carlo simulations with a linear response method. , 1997, Journal of medicinal chemistry.

[37]  Amedeo Caflisch,et al.  Wordom: a program for efficient analysis of molecular dynamics simulations , 2007, Bioinform..

[38]  A. Caflisch,et al.  Is quantum mechanics necessary for predicting binding free energy? , 2008, Journal of medicinal chemistry.

[39]  K. Dill,et al.  Binding of small-molecule ligands to proteins: "what you see" is not always "what you get". , 2009, Structure.

[40]  Alexander D. MacKerell,et al.  CHARMM general force field: A force field for drug‐like molecules compatible with the CHARMM all‐atom additive biological force fields , 2009, J. Comput. Chem..

[41]  L. Nieba,et al.  Competition BIAcore for measuring true affinities: large differences from values determined from binding kinetics. , 1996, Analytical biochemistry.

[42]  V. Pande,et al.  On the transition coordinate for protein folding , 1998 .

[43]  Vijay S Pande,et al.  Protein folded states are kinetic hubs , 2010, Proceedings of the National Academy of Sciences.

[44]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[45]  Amedeo Caflisch,et al.  Amyloid fibril polymorphism is under kinetic control. , 2010, Journal of the American Chemical Society.

[46]  A. Cavalli,et al.  Single-molecule pulling simulations can discern active from inactive enzyme inhibitors. , 2010, Journal of the American Chemical Society.

[47]  Amedeo Caflisch,et al.  Library screening by fragment‐based docking , 2009, Journal of molecular recognition : JMR.

[48]  Amedeo Caflisch,et al.  One-dimensional barrier-preserving free-energy projections of a beta-sheet miniprotein: new insights into the folding process. , 2008, The journal of physical chemistry. B.

[49]  G. Hummer,et al.  Coarse master equations for peptide folding dynamics. , 2008, The journal of physical chemistry. B.

[50]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

[51]  Sergei V Krivov,et al.  Diffusive reaction dynamics on invariant free energy profiles , 2008, Proceedings of the National Academy of Sciences.

[52]  Benoît Roux,et al.  Molecular basis for the Born model of ion solvation , 1990 .

[53]  Johan Åqvist,et al.  Ligand binding affinity prediction by linear interaction energy methods , 1998, J. Comput. Aided Mol. Des..

[54]  William Swope,et al.  Describing Protein Folding Kinetics by Molecular Dynamics Simulations. 1. Theory , 2004 .

[55]  George S. Hammond,et al.  A Correlation of Reaction Rates , 1955 .

[56]  A Caflisch,et al.  Identification of the protein folding transition state from molecular dynamics trajectories. , 2009, The Journal of chemical physics.

[57]  K Wüthrich,et al.  Protein hydration studied with homonuclear 3D1H NMR experiments , 1991, Journal of biomolecular NMR.

[58]  F. J. Luque,et al.  Binding site detection and druggability index from first principles. , 2009, Journal of medicinal chemistry.

[59]  K. Dill,et al.  Automatic discovery of metastable states for the construction of Markov models of macromolecular conformational dynamics. , 2007, The Journal of chemical physics.

[60]  C. Schütte,et al.  Supplementary Information for “ Constructing the Equilibrium Ensemble of Folding Pathways from Short Off-Equilibrium Simulations ” , 2009 .

[61]  Eric J. Deeds,et al.  Understanding ensemble protein folding at atomic detail , 2006, Proceedings of the National Academy of Sciences.

[62]  K Wüthrich,et al.  Interaction of urea with an unfolded protein The DNA‐binding domain of the 434‐repressor , 1995, FEBS letters.

[63]  Amedeo Caflisch,et al.  Network and graph analyses of folding free energy surfaces. , 2006, Current opinion in structural biology.

[64]  Luyuan Zhang,et al.  Mapping hydration dynamics around a protein surface , 2007, Proceedings of the National Academy of Sciences.

[65]  F. Noé,et al.  Transition networks for modeling the kinetics of conformational change in macromolecules. , 2008, Current opinion in structural biology.

[66]  A. Berezhkovskii,et al.  Reactive flux and folding pathways in network models of coarse-grained protein dynamics. , 2009, The Journal of chemical physics.