Extracting ligands from receptors by reversed targeted molecular dynamics

Short targeted MD trajectories are used to expel ligands from binding sites. The expulsion is governed by a linear increase of the target RMSD value, growing from zero to an arbitrary chosen final RMSD that forces the ligand to a selected distance outside of the receptor. The RMSD lag (i.e., the difference between the imposed and the actual RMSD) can be used to follow barriers encountered by the ligand during its way out of the receptor. The force constant used for the targeted MD can transform the RMSD lag into a strain energy. Integration of the (time-dependent) strain energy over time yields a value with the dimensions of “action” (i.e, energy multiplied by time) and can serve as a measure for the overall effort required to extract the ligand from its binding site. Possibilities to compare (numerically and graphically) the randomly detected exit pathways are discussed. As an example, the method is tested on the exit of bisphenol A from the human estrogen-related receptor $$\gamma$$γ and of GW0072 from the peroxysome proliferator activated receptor.

[1]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[2]  Araz Jakalian,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: I. Method , 2000 .

[3]  Helmut Grubmüller,et al.  Ligand-release pathways in the pheromone-binding protein of Bombyx mori. , 2006, Structure.

[4]  Woody Sherman,et al.  Contributions of water transfer energy to protein‐ligand association and dissociation barriers: Watermap analysis of a series of p38α MAP kinase inhibitors , 2013, Proteins.

[5]  C. Abrams,et al.  Ligand escape pathways and (un)binding free energy calculations for the hexameric insulin-phenol complex. , 2008, Biophysical journal.

[6]  Christopher I. Bayly,et al.  Fast, efficient generation of high‐quality atomic charges. AM1‐BCC model: II. Parameterization and validation , 2002, J. Comput. Chem..

[7]  Lennart Nilsson,et al.  Unbinding of retinoic acid from the retinoic acid receptor by random expulsion molecular dynamics. , 2006, Biophysical journal.

[8]  Arne Traulsen,et al.  Evolutionary Dynamics of Strategic Behavior in a Collective-Risk Dilemma , 2012, PLoS Comput. Biol..

[9]  Markus A. Lill,et al.  Including ligand‐induced protein flexibility into protein tunnel prediction , 2014, J. Comput. Chem..

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

[11]  P. Kollman,et al.  Automatic atom type and bond type perception in molecular mechanical calculations. , 2006, Journal of molecular graphics & modelling.

[12]  R. Wade,et al.  Comparison of the dynamics of substrate access channels in three cytochrome P450s reveals different opening mechanisms and a novel functional role for a buried arginine , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[13]  M. Peräkylä.  Ligand unbinding pathways from the vitamin D receptor studied by molecular dynamics simulations , 2009, European Biophysics Journal.

[14]  W. Delano.  The PyMOL Molecular Graphics System , 2002 .

[15]  Antonín Pavelka,et al.  CAVER 3.0: A Tool for the Analysis of Transport Pathways in Dynamic Protein Structures , 2012, PLoS Comput. Biol..

[16]  M. Kimura,et al.  Structural evidence for endocrine disruptor bisphenol A binding to human nuclear receptor ERR gamma. , 2007, Journal of biochemistry.

[17]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[18]  N. Garnier,et al.  Ligand-escape pathways from the ligand-binding domain of PPARγ receptor as probed by molecular dynamics simulations , 2008, European Biophysics Journal.

[19]  R. Wade,et al.  How do substrates enter and products exit the buried active site of cytochrome P450cam? 1. Random expulsion molecular dynamics investigation of ligand access channels and mechanisms. , 2000, Journal of molecular biology.

[20]  Gabriele Costantino,et al.  Unbinding Pathways of VEGFR2 Inhibitors Revealed by Steered Molecular Dynamics , 2014, J. Chem. Inf. Model..

[21]  R C Wade,et al.  How do substrates enter and products exit the buried active site of cytochrome P450cam? 2. Steered molecular dynamics and adiabatic mapping of substrate pathways. , 2000, Journal of molecular biology.

[22]  Matteo Masetti,et al.  Insights into Ligand–Protein Binding from Local Mechanical Response , 2011, Journal of chemical theory and computation.

[23]  Daniel R Roe,et al.  PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. , 2013, Journal of chemical theory and computation.

[24]  R. Copeland,et al.  Drug–target residence time and its implications for lead optimization , 2007, Nature Reviews Drug Discovery.

[25]  Leandro Martínez,et al.  Molecular dynamics simulations of ligand dissociation from thyroid hormone receptors: evidence of the likeliest escape pathway and its implications for the design of novel ligands. , 2006, Journal of medicinal chemistry.

[26]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..