Calculation of Relative Binding Free Energy in the Water-Filled Active Site of Oligopeptide-Binding Protein A

The periplasmic oligopeptide binding protein A (OppA) represents a well-known example of water-mediated protein-ligand interactions. Here, we perform free-energy calculations for three different ligands binding to OppA, using a thermodynamic integration approach. The tripeptide ligands share a high structural similarity (all have the sequence KXK), but their experimentally-determined binding free energies differ remarkably. Thermodynamic cycles were constructed for the ligands, and simulations conducted in the bound and (freely solvated) unbound states. In the unbound state, it was observed that the difference in conformational freedom between alanine and glycine leads to a surprisingly slow convergence, despite their chemical similarity. This could be overcome by increasing the softness parameter during alchemical transformations. Discrepancies remained in the bound state however, when comparing independent simulations of the three ligands. These difficulties could be traced to a slow relaxation of the water network within the active site. Fluctuations in the number of water molecules residing in the binding cavity occur mostly on a timescale larger than the simulation time along the alchemical path. After extensive simulations, relative binding free energies that were converged to within thermal noise could be obtained, which agree well with available experimental data.

[1]  D. Oliver,et al.  Lyme disease-causing Borrelia species encode multiple lipoproteins homologous to peptide-binding proteins of ABC-type transporters. , 1998, Infection and immunity.

[2]  C. Löw,et al.  Escherichia coli peptide binding protein OppA has a preference for positively charged peptides. , 2011, Journal of molecular biology.

[3]  C. Higgins,et al.  The crystal structures of the oligopeptide-binding protein OppA complexed with tripeptide and tetrapeptide ligands. , 1995, Structure.

[4]  Markus Christen,et al.  Architecture, implementation and parallelisation of the GROMOS software for biomolecular simulation , 2012, Comput. Phys. Commun..

[5]  Wilfred F. van Gunsteren,et al.  A generalized reaction field method for molecular dynamics simulations , 1995 .

[6]  F. Tian,et al.  Why OppA protein can bind sequence-independent peptides? A combination of QM/MM, PB/SA, and structure-based QSAR analyses , 2011, Amino Acids.

[7]  R E Hubbard,et al.  Relating structure to thermodynamics: The crystal structures and binding affinity of eight OppA‐peptide complexes , 1999, Protein science : a publication of the Protein Society.

[8]  Themis Lazaridis,et al.  Water at biomolecular binding interfaces. , 2007, Physical chemistry chemical physics : PCCP.

[9]  E J Dodson,et al.  Peptide binding in OppA, the crystal structures of the periplasmic oligopeptide binding protein in the unliganded form and in complex with lysyllysine. , 1997, Biochemistry.

[10]  J. Tame,et al.  Crystallographic and calorimetric analysis of peptide binding to OppA protein. , 1999, Journal of molecular biology.

[11]  Tim N. Heinz,et al.  Comparison of four methods to compute the dielectric permittivity of liquids from molecular dynamics simulations , 2001 .

[12]  V. Monnet Bacterial oligopeptide-binding proteins , 2003, Cellular and Molecular Life Sciences CMLS.

[13]  Ting Wang,et al.  Comparative binding energy (COMBINE) analysis of OppA-peptide complexes to relate structure to binding thermodynamics. , 2002, Journal of medicinal chemistry.

[14]  Chris Oostenbrink,et al.  Comparison of thermodynamic integration and Bennett acceptance ratio for calculating relative protein‐ligand binding free energies , 2013, J. Comput. Chem..

[15]  Themis Lazaridis,et al.  The effect of water displacement on binding thermodynamics: concanavalin A. , 2005, The journal of physical chemistry. B.

[16]  N. Vermeulen,et al.  The role of water molecules in computational drug design. , 2010, Current topics in medicinal chemistry.

[17]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[18]  A. Mark,et al.  Avoiding singularities and numerical instabilities in free energy calculations based on molecular simulations , 1994 .

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

[20]  C. Higgins,et al.  Periplasmic protein associated with the oligopeptide permeases of Salmonella typhimurium and Escherichia coli , 1983, Journal of bacteriology.

[21]  John E. Ladbury,et al.  The role of water in sequence-independent ligand binding by an oligopeptide transporter protein , 1996, Nature Structural Biology.

[22]  C. Higgins,et al.  Uptake of cell wall peptides by Salmonella typhimurium and Escherichia coli , 1987, Journal of bacteriology.

[23]  Yi Wang,et al.  Independent-Trajectory Thermodynamic Integration: a practical guide to protein-drug binding free energy calculations using distributed computing. , 2012, Methods in molecular biology.

[24]  J. Kirkwood Statistical Mechanics of Fluid Mixtures , 1935 .

[25]  Wilfred F. van Gunsteren,et al.  An improved GROMOS96 force field for aliphatic hydrocarbons in the condensed phase , 2001, J. Comput. Chem..

[26]  B. Poolman,et al.  The structural basis for peptide selection by the transport receptor OppA , 2009, The EMBO journal.

[27]  Michael R. Shirts,et al.  Optimal pairwise and non-pairwise alchemical pathways for free energy calculations of molecular transformation in solution phase. , 2012, The Journal of chemical physics.

[28]  M Levitt,et al.  Water: now you see it, now you don't. , 1993, Structure.

[29]  R. Wootton,et al.  Structural colour: Now you see it — now you don't , 2001, Nature.

[30]  Wilfred F van Gunsteren,et al.  Calculation of relative free energies for ligand-protein binding, solvation, and conformational transitions using the GROMOS software. , 2011, The journal of physical chemistry. B.

[31]  L. Larsson,et al.  A new class of synthetic antibacterials acting on lipopolysaccharide biosynthesis , 1987, Nature.

[32]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[33]  N. Panopoulos,et al.  Phaseolotoxin transport in Escherichia coli and Salmonella typhimurium via the oligopeptide permease , 1980, Journal of bacteriology.