Solvation Structure and Thermodynamic Mapping (SSTMap): An Open-Source, Flexible Package for the Analysis of Water in Molecular Dynamics Trajectories.

We have developed SSTMap, a software package for mapping structural and thermodynamic water properties in molecular dynamics trajectories. The package introduces automated analysis and mapping of local measures of frustration and enhancement of water structure. The thermodynamic calculations are based on Inhomogeneous Fluid Solvation Theory (IST), which is implemented using both site-based and grid-based approaches. The package also extends the applicability of solvation analysis calculations to multiple molecular dynamics (MD) simulation programs by using existing cross-platform tools for parsing MD parameter and trajectory files. SSTMap is implemented in Python and contains both command-line tools and a Python module to facilitate flexibility in setting up calculations and for automated generation of large data sets involving analysis of multiple solutes. Output is generated in formats compatible with popular Python data science packages. This tool will be used by the molecular modeling community for computational analysis of water in problems of biophysical interest such as ligand binding and protein function.

[1]  B. Berne,et al.  Role of the active-site solvent in the thermodynamics of factor Xa ligand binding. , 2008, Journal of the American Chemical Society.

[2]  Gerhard Hummer,et al.  Molecular binding: Under water's influence. , 2010, Nature chemistry.

[3]  Gerhard Klebe,et al.  Methyl, Ethyl, Propyl, Butyl: Futile But Not for Water, as the Correlation of Structure and Thermodynamic Signature Shows in a Congeneric Series of Thermolysin Inhibitors , 2014, ChemMedChem.

[4]  Michael S Bodnarchuk,et al.  Water, water, everywhere… It's time to stop and think. , 2016, Drug discovery today.

[5]  G. Whitesides,et al.  Is it the shape of the cavity, or the shape of the water in the cavity? , 2014 .

[6]  J. Ulander,et al.  Ligand Binding Thermodynamics in Drug Discovery: Still a Hot Tip? , 2015, Journal of medicinal chemistry.

[7]  I. Weber,et al.  Structural and kinetic analysis of caspase-3 reveals role for s5 binding site in substrate recognition. , 2006, Journal of molecular biology.

[8]  Michael K. Gilson,et al.  Solvation thermodynamic mapping of molecular surfaces in AmberTools: GIST , 2016, J. Comput. Chem..

[9]  Michael K. Gilson,et al.  Thermodynamics of Water in an Enzyme Active Site: Grid-Based Hydration Analysis of Coagulation Factor Xa , 2014, Journal of chemical theory and computation.

[10]  Alexander D. MacKerell,et al.  Spatial analysis and quantification of the thermodynamic driving forces in protein-ligand binding: binding site variability. , 2015, Journal of the American Chemical Society.

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

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

[13]  Lars V. Schäfer,et al.  Solvent effects on ligand binding to a serine protease. , 2017, Physical chemistry chemical physics : PCCP.

[14]  G. Whitesides,et al.  Water networks contribute to enthalpy/entropy compensation in protein-ligand binding. , 2013, Journal of the American Chemical Society.

[15]  Robert Abel,et al.  Motifs for molecular recognition exploiting hydrophobic enclosure in protein–ligand binding , 2007, Proceedings of the National Academy of Sciences.

[16]  Mike C. Payne,et al.  Assessing the Accuracy of Inhomogeneous Fluid Solvation Theory in Predicting Hydration Free Energies of Simple Solutes , 2013, The journal of physical chemistry. B.

[17]  Lauren Wickstrom,et al.  Enthalpic Breakdown of Water Structure on Protein Active-Site Surfaces. , 2016, The journal of physical chemistry. B.

[18]  Travis E. Oliphant,et al.  Python for Scientific Computing , 2007, Computing in Science & Engineering.

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

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

[21]  Gerhard Klebe,et al.  Rational Design of Thermodynamic and Kinetic Binding Profiles by Optimizing Surface Water Networks Coating Protein-Bound Ligands. , 2016, Journal of medicinal chemistry.

[22]  M. Gilson,et al.  Spatial Decomposition of Translational Water–Water Correlation Entropy in Binding Pockets , 2015, Journal of chemical theory and computation.

[23]  Gerhard Klebe,et al.  Impact of Surface Water Layers on Protein-Ligand Binding: How Well Are Experimental Data Reproduced by Molecular Dynamics Simulations in a Thermolysin Test Case? , 2016, J. Chem. Inf. Model..

[24]  Michael K Gilson,et al.  Grid inhomogeneous solvation theory: hydration structure and thermodynamics of the miniature receptor cucurbit[7]uril. , 2012, The Journal of chemical physics.

[25]  Michael K Gilson,et al.  Testing inhomogeneous solvation theory in structure-based ligand discovery , 2017, Proceedings of the National Academy of Sciences.

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

[27]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[28]  Gerhard Klebe,et al.  Dissecting the hydrophobic effect on the molecular level: the role of water, enthalpy, and entropy in ligand binding to thermolysin. , 2013, Angewandte Chemie.

[29]  Gaël Varoquaux,et al.  The NumPy Array: A Structure for Efficient Numerical Computation , 2011, Computing in Science & Engineering.

[30]  C. Simmerling,et al.  ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. , 2015, Journal of chemical theory and computation.

[31]  Themis Lazaridis,et al.  Orientational correlations and entropy in liquid water , 1996 .

[32]  Federico D. Sacerdoti,et al.  Scalable Algorithms for Molecular Dynamics Simulations on Commodity Clusters , 2006, ACM/IEEE SC 2006 Conference (SC'06).

[33]  Felice C. Lightstone,et al.  Accounting for water molecules in drug design , 2011, Expert opinion on drug discovery.

[34]  Marcelo A. Marti,et al.  WATCLUST: a tool for improving the design of drugs based on protein-water interactions , 2015, Bioinform..

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

[36]  Philip Ball,et al.  Water as an Active Constituent in Cell Biology , 2008 .

[37]  George M. Whitesides,et al.  Mechanism of the hydrophobic effect in the biomolecular recognition of arylsulfonamides by carbonic anhydrase , 2011, Proceedings of the National Academy of Sciences.

[38]  D. Huggins.  Estimating Translational and Orientational Entropies Using the k-Nearest Neighbors Algorithm. , 2014, Journal of chemical theory and computation.

[39]  Robert Abel,et al.  Thermodynamic properties of liquid water: an application of a nonparametric approach to computing the entropy of a neat fluid. , 2009, Journal of chemical theory and computation.

[40]  Ricardo L Mancera.  Molecular modeling of hydration in drug design. , 2007, Current opinion in drug discovery & development.

[41]  Thomas J Lane,et al.  MDTraj: a modern, open library for the analysis of molecular dynamics trajectories , 2014, bioRxiv.

[42]  Diwakar Shukla,et al.  OpenMM 4: A Reusable, Extensible, Hardware Independent Library for High Performance Molecular Simulation. , 2013, Journal of chemical theory and computation.

[43]  J. Ladbury.  Just add water! The effect of water on the specificity of protein-ligand binding sites and its potential application to drug design. , 1996, Chemistry & biology.

[44]  W. Sherman,et al.  Calculating Water Thermodynamics in the Binding Site of Proteins - Applications of WaterMap to Drug Discovery. , 2017, Current topics in medicinal chemistry.

[45]  J. Onuchic,et al.  Water mediation in protein folding and molecular recognition. , 2006, Annual review of biophysics and biomolecular structure.

[46]  Felice C Lightstone,et al.  Approaches to efficiently estimate solvation and explicit water energetics in ligand binding: the use of WaterMap , 2013, Expert opinion on drug discovery.

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

[48]  K. Sharp,et al.  Water in the half shell: structure of water, focusing on angular structure and solvation. , 2010, Accounts of chemical research.

[49]  R. Baron,et al.  Water in Cavity−Ligand Recognition , 2010, Journal of the American Chemical Society.

[50]  G. Whitesides,et al.  Water-Restructuring Mutations Can Reverse the Thermodynamic Signature of Ligand Binding to Human Carbonic Anhydrase. , 2017, Angewandte Chemie.

[51]  Themis Lazaridis,et al.  Inhomogeneous Fluid Approach to Solvation Thermodynamics. 2. Applications to Simple Fluids , 1998 .

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

[53]  Andrew T. Fenley,et al.  Bridging Calorimetry and Simulation through Precise Calculations of Cucurbituril–Guest Binding Enthalpies , 2014, Journal of chemical theory and computation.

[54]  T. Lazaridis.  Inhomogeneous Fluid Approach to Solvation Thermodynamics. 1. Theory , 1998 .

[55]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..