Rapid approximate calculation of water binding free energies in the whole hydration domain of (bio)macromolecules

Abstract

The evaluation of water binding free energies around solute molecules is important for the thermodynamic characterization of hydration or association processes. Here, a rapid approximate method to estimate water binding free energies around (bio)macromolecules from a single molecular dynamics simulation is presented. The basic idea is that endpoint free-energy calculation methods are applied and the endpoint quantities are monitored on a three-dimensional grid around the solute. Thus, a gridded map of water binding free energies around the solute is obtained, that is, from a single short simulation, a map of favorable and unfavorable water binding sites can be constructed. Among the employed free-energy calculation methods, approaches involving endpoint information pertaining to actual thermodynamic integration calculations or endpoint information as exploited in the linear interaction energy method were examined. The accuracy of the approximate approaches was evaluated on the hydration of a cage-like molecule representing either a nonpolar, polar, or charged water binding site and on α- and β-cyclodextrin molecules. Among the tested approaches, the linear interaction energy method is considered the most viable approach. Applying the linear interaction energy method on the grid around the solute, a semi-quantitative thermodynamic characterization of hydration around the whole solute is obtained. Disadvantages are the approximate nature of the method and a limited flexibility of the solute. © 2016 Wiley Periodicals, Inc.

Insights into Noncovalent Binding Obtained from Molecular Dynamics Simulations

“…Illustration of water binding free energy in the vicinity of the a) monomeric and b) dimeric solute molecules. The calculation is based on a linear interaction energy approach . Blue color represents favorable binding of water while red color represents unfavorable interaction.…”

Human Immune Protein C1q Selectively Disaggregates Carbon Nanotubes

“…Various methods exist for the calculation of free energy differences from MD. − When only liquids and solids are involved as is the case here, we can work with the Helmholtz free energy and consider constant temperature ( T ) and volume ( V ) conditions. For the problem in hand, it is best to obtain the free energy directly from the averaged forces.…”

Moisture sorption and thermodynamic properties of Camellia oleifera seeds as influenced by oil content

“…At constant temperature, when 1 mol water is removed from the material, the additional energy is called binding energy after removing latent heat of vaporization [22,23] , which is evaluated by Equation ( 4):…”