Dataset Information

Binding enthalpy calculations for a neutral host-guest pair yield widely divergent salt effects across water models.

ABSTRACT: Dissolved salts are a part of the physiological milieu and can significantly influence the kinetics and thermodynamics of various biomolecular processes, such as binding and catalysis; thus, it is important for molecular simulations to reliably describe their effects. The present study uses a simple, nonionized host-guest model system to study the sensitivity of computed binding enthalpies to the choice of water and salt models. Molecular dynamics simulations of a cucurbit[7]uril host with a neutral guest molecule show striking differences in the salt dependency of the binding enthalpy across four water models, TIP3P, SPC/E, TIP4P-Ew, and OPC, with additional sensitivity to the choice of parameters for sodium and chloride. In particular, although all of the models predict that binding will be less exothermic with increasing NaCl concentration, the strength of this effect varies by 7 kcal/mol across models. The differences appear to result primarily from differences in the number of sodium ions displaced from the host upon binding the guest rather than from differences in the enthalpy associated with this displacement, and it is the electrostatic energy that contributes most to the changes in enthalpy with increasing salt concentration. That a high sensitivity of salt affecting the choice of water model, as observed for the present host-guest system despite it being nonionized, raises issues regarding the selection and adjustment of water models for use with biological macromolecules, especially as these typically possess multiple ionized groups that can interact relatively strongly with ions in solution.

SUBMITTER: Gao K

PROVIDER: S-EPMC4668924 | biostudies-literature | 2015 Oct

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Binding enthalpy calculations for a neutral host-guest pair yield widely divergent salt effects across water models.

Gao Kaifu K Yin Jian J Henriksen Niel M NM Fenley Andrew T AT Gilson Michael K MK

Journal of chemical theory and computation 20150918 10

Dissolved salts are a part of the physiological milieu and can significantly influence the kinetics and thermodynamics of various biomolecular processes, such as binding and catalysis; thus, it is important for molecular simulations to reliably describe their effects. The present study uses a simple, nonionized host-guest model system to study the sensitivity of computed binding enthalpies to the choice of water and salt models. Molecular dynamics simulations of a cucurbit[7]uril host with a neu ...[more]

PMID: 26574247

Similar Datasets

Project description:Calculations using the density-functional theory (DFT) in combination with the single defect method were carried out to determine the heat of mixing behaviour of mineral solid solution phases. The accuracy of this method was tested on the halite-sylvite (NaCl-KCl) binary, pyrope-grossular garnets (Mg3Al2Si3O12-Ca3Al2Si3O12), MgO-CaO (halite structure) binary, and on Al/Si ordered alkali feldspars (NaAlSi3O8-KAlSi3O8); as members for coupled substitutions, the diopside-jadeite pyroxenes (CaMgSi2O6-NaAlSi2O6) and diopside-CaTs pyroxenes (CaMgSi2O6-CaAlAlSiO6) were chosen for testing and, as an application, the heat of mixing of the tremolite-glaucophane amphiboles (Ca2Mg5Si8O22(OH)2-Na2Mg3Al2Si8O22(OH)2) was computed. Six of these binaries were selected because of their experimentally well-known thermodynamic mixing behaviours. The comparison of the calculated heat of mixing data with calorimetric data showed good agreement for halite-sylvite, pyrope-grossular, and diopside-jadeite binaries and small differences for the Al/Si ordered alkali feldspar solid solution. In the case of the diopside-CaTs binary, the situation is more complex because CaTs is an endmember with disordered cation distributions. Good agreement with the experimental data could be, however, achieved assuming a reasonable disordered state. The calculated data for the Al/Si ordered alkali feldspars were applied to phase equilibrium calculations, i.e. calculating the Al/Si ordered alkali feldspar solvus. This solvus was then compared to the experimentally determined solvus finding good agreement. The solvus of the MgO-CaO binary was also constructed from DFT-based data and compared to the experimentally determined solvus, and the two were also in good agreement. Another application was the determination of the solvus in tremolite-glaucophane amphiboles (Ca2Mg5Si8O22(OH)2-Na2Mg3Al2Si8O22(OH)2). It was compared to solvi based on coexisting amphiboles found in eclogites and phase equilibrium experiments.

Project description:Computational prediction of noncovalent binding free energies with methods based on molecular mechanical force fields has become increasingly routine in drug discovery projects, where they promise to speed the discovery of small molecule ligands to bind targeted proteins with high affinity. Because the reliability of free energy methods still has significant room for improvement, new force fields, or modifications of existing ones, are regularly introduced with the aim of improving the accuracy of molecular simulations. However, comparatively little work has been done to systematically assess how well force fields perform, particularly in relation to the calculation of binding affinities. Hardware advances have made these calculations feasible, but comprehensive force field assessments for protein-ligand sized systems still remain costly. Here, we turn to cyclodextrin host-guest systems, which feature many hallmarks of protein-ligand binding interactions but are generally much more tractable due to their small size. We present absolute binding free energy and enthalpy calculations, using the attach-pull-release (APR) approach, on a set of 43 cyclodextrin-guest pairs for which experimental ITC data are available. The test set comprises both α- and β-cyclodextrin hosts binding a series of small organic guests, each with one of three functional groups: ammonium, alcohol, or carboxylate. Four water models are considered (TIP3P, TIP4Pew, SPC/E, and OPC), along with two partial charge assignment procedures (RESP and AM1-BCC) and two cyclodextrin host force fields. The results suggest a complex set of considerations when choosing a force field for biomolecular simulations. For example, some force field combinations clearly outperform others at the binding enthalpy calculations but not for the binding free energy. Additionally, a force field combination which we expected to be the worst performer gave the most accurate binding free energies - but the least accurate binding enthalpies. The results have implications for the development of improved force fields, and we propose this test set, and potential future elaborations of it, as a powerful validation suite to evaluate new force fields and help guide future force field development.

Dataset Information

Binding enthalpy calculations for a neutral host-guest pair yield widely divergent salt effects across water models.

Publications

Binding enthalpy calculations for a neutral host-guest pair yield widely divergent salt effects across water models.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets