Deciphering the Complexity of Ligand-Protein Recognition Pathways Using Supervised Molecular Dynamics (SuMD) Simulations

Molecular recognition is a crucial issue when aiming to interpret the mechanism of known active substances as well as to develop novel active candidates. Unfortunately, simulating the binding process is still a challenging task because it requires classical MD experiments in a long microsecond time scale that are affordable only with a high-level computational capacity. In order to overcome this limiting factor, we have recently implemented an alternative MD approach, named supervised molecular dynamics (SuMD), and successfully applied it to G protein-coupled receptors (GPCRs). SuMD enables the investigation of ligand-receptor binding events independently from the starting position, chemical structure of the ligand, and also from its receptor binding affinity. In this article, we present an extension of the SuMD application domain including different types of proteins in comparison with GPCRs. In particular, we have deeply analyzed the ligand-protein recognition pathways of six different case studies that we grouped into two different classes: globular and membrane proteins. Moreover, we introduce the SuMD-Analyzer tool that we have specifically implemented to help the user in the analysis of the SuMD trajectories. Finally, we emphasize the limit of the SuMD applicability domain as well as its strengths in analyzing the complexity of ligand-protein recognition pathways.

[1]  Francesca Fanelli,et al.  Wordom: A User-Friendly Program for the Analysis of Molecular Structures, Trajectories, and Free Energy Surfaces , 2010, J. Comput. Chem..

[2]  Gisbert Schneider,et al.  Protein-ligand interactions from molecular recognition to drug design , 2003 .

[3]  W. V. van Gunsteren,et al.  A fast SHAKE algorithm to solve distance constraint equations for small molecules in molecular dynamics simulations , 2001 .

[4]  Harini Krishnamurthy,et al.  X-ray structures of LeuT in substrate-free outward-open and apo inward-open states , 2012, Nature.

[5]  I. Sjöholm,et al.  Binding of drugs to human serum albumin:XI. The specificity of three binding sites as studied with albumin immobilized in microparticles. , 1979, Molecular pharmacology.

[6]  Conrad C. Huang,et al.  UCSF Chimera—A visualization system for exploratory research and analysis , 2004, J. Comput. Chem..

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

[8]  G. de Fabritiis,et al.  Complete reconstruction of an enzyme-inhibitor binding process by molecular dynamics simulations , 2011, Proceedings of the National Academy of Sciences.

[9]  S. Moro,et al.  Understanding allosteric interactions in G protein-coupled receptors using Supervised Molecular Dynamics: A prototype study analysing the human A3 adenosine receptor positive allosteric modulator LUF6000. , 2015, Bioorganic & medicinal chemistry.

[10]  B. Knoops,et al.  Crystal structure of human peroxiredoxin 5, a novel type of mammalian peroxiredoxin at 1.5 A resolution. , 2001, Journal of molecular biology.

[11]  Albert C. Pan,et al.  Molecular determinants of drug-receptor binding kinetics. , 2013, Drug discovery today.

[12]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[13]  Alexander D. MacKerell,et al.  Development and current status of the CHARMM force field for nucleic acids , 2000, Biopolymers.

[14]  Klaus Schulten,et al.  Rapid parameterization of small molecules using the force field toolkit , 2013, J. Comput. Chem..

[15]  J. Guichou,et al.  Comparing Binding Modes of Analogous Fragments Using NMR in Fragment-Based Drug Design: Application to PRDX5 , 2014, PloS one.

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

[17]  Andrei L. Lomize,et al.  OPM: Orientations of Proteins in Membranes database , 2006, Bioinform..

[18]  P. Labute.  proteins STRUCTURE O FUNCTION O BIOINFORMATICS Protonate3D: Assignment of ionization , 2013 .

[19]  Sara Lejon,et al.  Structural basis for the binding of naproxen to human serum albumin in the presence of fatty acids and the GA module. , 2008, Acta crystallographica. Section F, Structural biology and crystallization communications.

[20]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[21]  E. Laborde.  Glutathione transferases as mediators of signaling pathways involved in cell proliferation and cell death , 2010, Cell Death and Differentiation.

[22]  J. Stewart.  Optimization of parameters for semiempirical methods V: Modification of NDDO approximations and application to 70 elements , 2007, Journal of molecular modeling.

[23]  K. Kitaura,et al.  Structural insight into human CK2alpha in complex with the potent inhibitor ellagic acid. , 2009, Bioorganic & medicinal chemistry letters.

[24]  Michael J. Frisch,et al.  MP2 energy evaluation by direct methods , 1988 .

[25]  S. Moro,et al.  Role of the extracellular loops of G protein-coupled receptors in ligand recognition: a molecular modeling study of the human P2Y1 receptor. , 1999, Biochemistry.

[26]  M. Filizola,et al.  Beyond standard molecular dynamics: investigating the molecular mechanisms of G protein-coupled receptors with enhanced molecular dynamics methods. , 2014, Advances in experimental medicine and biology.

[27]  E. Zomot,et al.  Sodium-coupled neurotransmitter transporters. , 2008, Chemical reviews.

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

[29]  M. Reith,et al.  Antidepressant specificity of serotonin transporter suggested by three LeuT–SSRI structures , 2009, Nature Structural &Molecular Biology.

[30]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[31]  Alexander D. MacKerell,et al.  Automation of the CHARMM General Force Field (CGenFF) I: Bond Perception and Atom Typing , 2012, J. Chem. Inf. Model..

[32]  M J Harvey,et al.  ACEMD: Accelerating Biomolecular Dynamics in the Microsecond Time Scale. , 2009, Journal of chemical theory and computation.

[33]  B. Knoops,et al.  Discovery of Fragment Molecules That Bind the Human Peroxiredoxin 5 Active Site , 2010, PloS one.

[34]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

[35]  Jonathan S. Mason,et al.  Structures of G protein-coupled receptors reveal new opportunities for drug discovery. , 2015, Drug discovery today.

[36]  S. Moro,et al.  How druggable is protein kinase CK2? , 2010, Medicinal research reviews.

[37]  S. Moro,et al.  Exploring the recognition pathway at the human A2A adenosine receptor of the endogenous agonist adenosine using supervised molecular dynamics simulations , 2015 .

[38]  B. Brooks,et al.  Langevin dynamics of peptides: The frictional dependence of isomerization rates of N‐acetylalanyl‐N′‐methylamide , 1992, Biopolymers.

[39]  Hugo Gutiérrez-de-Terán,et al.  Structural and Energetic Effects of A2A Adenosine Receptor Mutations on Agonist and Antagonist Binding , 2014, PloS one.

[40]  Stefano Moro,et al.  Supervised Molecular Dynamics (SuMD) as a Helpful Tool To Depict GPCR-Ligand Recognition Pathway in a Nanosecond Time Scale , 2014, J. Chem. Inf. Model..

[41]  A. Leslie,et al.  Agonist-bound adenosine A2A receptor structures reveal common features of GPCR activation , 2011, Nature.

[42]  S. Moro,et al.  Advances in Computational Techniques to Study GPCR-Ligand Recognition. , 2015, Trends in pharmacological sciences.

[43]  M. Parker,et al.  Structure and function of glutathione S-transferases. , 1994, Biochimica et biophysica acta.

[44]  K. Jacobson,et al.  Adenosine receptors as therapeutic targets , 2006, Nature Reviews Drug Discovery.

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

[46]  M W Parker,et al.  The ligandin (non-substrate) binding site of human Pi class glutathione transferase is located in the electrophile binding site (H-site). , 1999, Journal of molecular biology.

[47]  Alexander D. MacKerell,et al.  Automation of the CHARMM General Force Field (CGenFF) II: Assignment of Bonded Parameters and Partial Atomic Charges , 2012, J. Chem. Inf. Model..

[48]  Qiang Zhao,et al.  Modeling ligand recognition at the P2Y12 receptor in light of X-ray structural information , 2015, Journal of Computer-Aided Molecular Design.

[49]  R. Morita,et al.  Peroxiredoxin family proteins are key initiators of post-ischemic inflammation in the brain , 2012, Nature Medicine.