Role of Water in Ligand Binding to Maltose-Binding Protein: Insight from a New Docking Protocol Based on the 3D-RISM-KH Molecular Theory of Solvation

Maltose-binding protein is a periplasmic binding protein responsible for transport of maltooligosaccarides through the periplasmic space of Gram-negative bacteria, as a part of the ABC transport system. The molecular mechanisms of the initial ligand binding and induced large scale motion of the protein's domains still remain elusive. In this study, we use a new docking protocol that combines a recently proposed explicit water placement algorithm based on the 3D-RISM-KH molecular theory of solvation and conventional docking software (AutoDock Vina) to explain the mechanisms of maltotriose binding to the apo-open state of a maltose-binding protein. We confirm the predictions of previous NMR spectroscopic experiments on binding modes of the ligand. We provide the molecular details on the binding mode that was not previously observed in the X-ray experiments. We show that this mode, which is defined by the fine balance between the protein-ligand direct interactions and solvation effects, can trigger the protein's domain motion resulting in the holo-closed structure of the maltose-binding protein with the maltotriose ligand in excellent agreement with the experimental data. We also discuss the role of water in blocking unfavorable binding sites and water-mediated interactions contributing to the stability of observable binding modes of maltotriose.

[1]  Robin Taylor,et al.  SuperStar: a knowledge-based approach for identifying interaction sites in proteins. , 1999, Journal of molecular biology.

[2]  Yuichi Harano,et al.  A theoretical analysis on hydration thermodynamics of proteins. , 2006, The Journal of chemical physics.

[3]  Andriy Kovalenko,et al.  3D-RISM-KH approach for biomolecular modelling at nanoscale: thermodynamics of fibril formation and beyond , 2010 .

[4]  David S. Goodsell,et al.  AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility , 2009, J. Comput. Chem..

[5]  Ricardo L. Mancera,et al.  Ligand-Protein Docking with Water Molecules , 2008, J. Chem. Inf. Model..

[6]  Winfried Boos,et al.  Maltose Transport in Escherichia coli K12 , 1976 .

[7]  L. Kay,et al.  The energetic cost of domain reorientation in maltose-binding protein as studied by NMR and fluorescence spectroscopy , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[8]  Caterina Barillari,et al.  Classification of water molecules in protein binding sites. , 2007, Journal of the American Chemical Society.

[9]  Bernard Pettitt,et al.  A Site‐Site Theory for Finite Concentration Saline Solutions. , 1993 .

[10]  Benoît Roux,et al.  NUMERICAL SOLUTION OF THE HYPERNETTED CHAIN EQUATION FOR A SOLUTE OF ARBITRARY GEOMETRY IN THREE DIMENSIONS , 1995 .

[11]  M. Karplus,et al.  Functionality maps of binding sites: A multiple copy simultaneous search method , 1991, Proteins.

[12]  Gregory A Ross,et al.  Rapid and Accurate Prediction and Scoring of Water Molecules in Protein Binding Sites , 2012, PloS one.

[13]  Yuichi Harano,et al.  Theoretical study for partial molar volume of amino acids and polypeptides by the three-dimensional reference interaction site model , 2001 .

[14]  David Chandler,et al.  Density functional theory of nonuniform polyatomic systems. II: Rational closures for integral equations , 1986 .

[15]  I. R. Mcdonald,et al.  Theory of simple liquids , 1998 .

[16]  Brian K. Shoichet,et al.  Roles for Ordered and Bulk Solvent in Ligand Recognition and Docking in Two Related Cavities , 2013, PloS one.

[17]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[18]  Andriy Kovalenko,et al.  Modeling solvatochromic shifts using the orbital-free embedding potential at statistically mechanically averaged solvent density. , 2010, The journal of physical chemistry. A.

[19]  Fumio Hirata,et al.  Analysis of biomolecular solvation sites by 3D-RISM theory. , 2013, The journal of physical chemistry. B.

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

[21]  Andriy Kovalenko,et al.  Association thermodynamics and conformational stability of beta-sheet amyloid beta(17-42) oligomers: effects of E22Q (Dutch) mutation and charge neutralization. , 2010, Biophysical journal.

[22]  W. L. Jorgensen,et al.  Energetics of displacing water molecules from protein binding sites: consequences for ligand optimization. , 2009, Journal of the American Chemical Society.

[23]  Andriy Kovalenko,et al.  Three-dimensional Rism Theory for Molecular Liquids and Solid-Liquid Interfaces , 2004 .

[24]  Florante A. Quiocho,et al.  Structural evidence for a dominant role of nonpolar interactions in the binding of a transport/chemosensory receptor to its highly polar ligands. , 2002, Biochemistry.

[25]  Fumio Hirata,et al.  Self-consistent description of a metal–water interface by the Kohn–Sham density functional theory and the three-dimensional reference interaction site model , 1999 .

[26]  W R Pitt,et al.  Modelling of solvent positions around polar groups in proteins. , 1991, Protein engineering.

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

[28]  Arthur J. Olson,et al.  AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading , 2009, J. Comput. Chem..

[29]  F A Quiocho,et al.  Extensive features of tight oligosaccharide binding revealed in high-resolution structures of the maltodextrin transport/chemosensory receptor. , 1997, Structure.

[30]  F A Quiocho,et al.  Crystal structures of the maltodextrin/maltose-binding protein complexed with reduced oligosaccharides: flexibility of tertiary structure and ligand binding. , 2001, Journal of molecular biology.

[31]  L. Kay,et al.  Ligand-induced structural changes to maltodextrin-binding protein as studied by solution NMR spectroscopy. , 2001, Journal of molecular biology.

[32]  G. Marius Clore,et al.  Open to closed transition in apo maltose-binding protein visualized by paramagnetic NMR , 2007 .

[33]  R. Doolittle,et al.  A simple method for displaying the hydropathic character of a protein. , 1982, Journal of molecular biology.

[34]  Hak-Sung Kim,et al.  A single-molecule dissection of ligand binding to a protein with intrinsic dynamics. , 2013, Nature chemical biology.

[35]  Fumio Hirata,et al.  Placevent: An algorithm for prediction of explicit solvent atom distribution—Application to HIV‐1 protease and F‐ATP synthase , 2012, J. Comput. Chem..

[36]  F A Quiocho,et al.  Rates of ligand binding to periplasmic proteins involved in bacterial transport and chemotaxis. , 1983, The Journal of biological chemistry.

[37]  F A Quiocho,et al.  Carbohydrate-binding proteins: tertiary structures and protein-sugar interactions. , 1986, Annual review of biochemistry.

[38]  F. Quiocho,et al.  Crystallographic evidence of a large ligand-induced hinge-twist motion between the two domains of the maltodextrin binding protein involved in active transport and chemotaxis. , 1992, Biochemistry.

[39]  Fumio Hirata,et al.  Ligand mapping on protein surfaces by the 3D-RISM theory: toward computational fragment-based drug design. , 2009, Journal of the American Chemical Society.

[40]  Diego F. Gauto,et al.  Solvent structure improves docking prediction in lectin-carbohydrate complexes. , 2013, Glycobiology.

[41]  Seok-Yong Lee,et al.  Detection and characterization of xenon-binding sites in proteins by 129Xe NMR spectroscopy. , 2002, Journal of molecular biology.

[42]  Makoto Taiji,et al.  Free-energy landscapes of protein domain movements upon ligand binding. , 2011, The journal of physical chemistry. B.

[43]  Christopher R. Corbeil,et al.  Towards the development of universal, fast and highly accurate docking/scoring methods: a long way to go , 2008, British journal of pharmacology.

[44]  K. Diederichs,et al.  X-ray structures of the maltose-maltodextrin-binding protein of the thermoacidophilic bacterium Alicyclobacillus acidocaldarius provide insight into acid stability of proteins. , 2003, Journal of molecular biology.

[45]  J G Pelton,et al.  Tritium NMR spectroscopy of ligand binding to maltose-binding protein. , 1991, Biochemistry.

[46]  Andriy Kovalenko,et al.  3D-RISM-Dock: A New Fragment-Based Drug Design Protocol. , 2012, Journal of chemical theory and computation.

[47]  Benoît Roux,et al.  An Integral Equation To Describe the Solvation of Polar Molecules in Liquid Water , 1997 .

[48]  Frederick C. Neidhardt,et al.  Escherichia coli and Salmonella :cellular and molecular biology , 2016 .

[49]  Yuichi Harano,et al.  Theoretical analysis on changes in thermodynamic quantities upon protein folding: essential role of hydration. , 2007, The Journal of chemical physics.

[50]  F A Quiocho,et al.  Refined 1.8-A structure reveals the mode of binding of beta-cyclodextrin to the maltodextrin binding protein. , 1993, Biochemistry.

[51]  Daiwen Yang,et al.  A new strategy for structure determination of large proteins in solution without deuteration , 2006, Nature Methods.

[52]  James Andrew McCammon,et al.  Extracting hydration sites around proteins from explicit water simulations , 2002, J. Comput. Chem..

[53]  P. Goodford.  A computational procedure for determining energetically favorable binding sites on biologically important macromolecules. , 1985, Journal of medicinal chemistry.

[54]  Karl Nicholas Kirschner,et al.  GLYCAM06: A generalizable biomolecular force field. Carbohydrates , 2008, J. Comput. Chem..

[55]  Masaru Yoshida,et al.  Molecular theory of solvation for supramolecules and soft matter structures: application to ligand binding, ion channels, and oligomeric polyelectrolyte gelators , 2012 .

[56]  M F Sanner,et al.  Python: a programming language for software integration and development. , 1999, Journal of molecular graphics & modelling.

[57]  F. Hirata,et al.  Three-dimensional density profiles of water in contact with a solute of arbitrary shape: a RISM approach , 1998 .

[58]  Andriy Kovalenko,et al.  Multiscale methods for nanochemistry and biophysics in solution , 2011 .

[59]  Julien Michel,et al.  Prediction of the water content in protein binding sites. , 2009, The journal of physical chemistry. B.

[60]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

[61]  Fumio Hirata,et al.  Locating missing water molecules in protein cavities by the three‐dimensional reference interaction site model theory of molecular solvation , 2006, Proteins.

[62]  J. Mccammon,et al.  Induced Fit or Conformational Selection? The Role of the Semi-closed State in the Maltose Binding Protein , 2011, Biochemistry.

[63]  Andriy Kovalenko,et al.  Calculation of local water densities in biological systems: a comparison of molecular dynamics simulations and the 3D-RISM-KH molecular theory of solvation. , 2011, The journal of physical chemistry. B.

[64]  Sandor Vajda,et al.  Algorithms for computational solvent mapping of proteins , 2003, Proteins.

[65]  P. Kollman,et al.  How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? , 2000 .

[66]  Niu Huang,et al.  Exploiting ordered waters in molecular docking. , 2008, Journal of medicinal chemistry.

[67]  Andriy Kovalenko,et al.  Hydration effects on the HET-s prion and amyloid-beta fibrillous aggregates, studied with three-dimensional molecular theory of solvation. , 2008, Biophysical journal.

[68]  W. Sherman,et al.  Thermodynamic analysis of water molecules at the surface of proteins and applications to binding site prediction and characterization , 2011, Proteins.

[69]  J M Sturtevant,et al.  A thermodynamic study of the binding of linear and cyclic oligosaccharides to the maltodextrin-binding protein of Escherichia coli. , 1998, Biophysical chemistry.

[70]  David Chandler,et al.  Density functional theory of nonuniform polyatomic systems. I. General formulation , 1986 .

[71]  Thomas Stockner,et al.  A salt-bridge motif involved in ligand binding and large-scale domain motions of the maltose-binding protein. , 2005, Biophysical journal.

[72]  Ashutosh Kumar,et al.  Investigation on the Effect of Key Water Molecules on Docking Performance in CSARdock Exercise , 2013, J. Chem. Inf. Model..

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