Dataset Information

Bind3P: Optimization of a Water Model Based on Host-Guest Binding Data.

ABSTRACT: We report a water model, Bind3P (Version 0.1), which was obtained by using sensitivity analysis to readjust the Lennard-Jones parameters of the TIP3P model against experimental binding free energies for six host-guest systems, along with pure liquid properties. Tests of Bind3P against >100 experimental binding free energies and enthalpies for host-guest systems distinct from the training set show a consistent drop in the mean signed error, relative to matched calculations with TIP3P. Importantly, Bind3P also yields some improvement in the hydration free energies of small organic molecules and preserves the accuracy of bulk water properties, such as density and the heat of vaporization. The same approach can be applied to more sophisticated water models that can better represent pure water properties. These results lend further support to the concept of integrating host-guest binding data into force field parametrization.

SUBMITTER: Yin J

PROVIDER: S-EPMC6039243 | biostudies-literature | 2018 Jul

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Bind3P: Optimization of a Water Model Based on Host-Guest Binding Data.

Yin Jian J Henriksen Niel M NM Muddana Hari S HS Gilson Michael K MK

Journal of chemical theory and computation 20180622 7

We report a water model, Bind3P (Version 0.1), which was obtained by using sensitivity analysis to readjust the Lennard-Jones parameters of the TIP3P model against experimental binding free energies for six host-guest systems, along with pure liquid properties. Tests of Bind3P against >100 experimental binding free energies and enthalpies for host-guest systems distinct from the training set show a consistent drop in the mean signed error, relative to matched calculations with TIP3P. Importantly ...[more]

PMID: 29874074

Similar Datasets

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.

Project description:Accurately predicting the binding affinities of small organic molecules to biological macromolecules can greatly accelerate drug discovery by reducing the number of compounds that must be synthesized to realize desired potency and selectivity goals. Unfortunately, the process of assessing the accuracy of current computational approaches to affinity prediction against binding data to biological macromolecules is frustrated by several challenges, such as slow conformational dynamics, multiple titratable groups, and the lack of high-quality blinded datasets. Over the last several SAMPL blind challenge exercises, host-guest systems have emerged as a practical and effective way to circumvent these challenges in assessing the predictive performance of current-generation quantitative modeling tools, while still providing systems capable of possessing tight binding affinities. Here, we present an overview of the SAMPL6 host-guest binding affinity prediction challenge, which featured three supramolecular hosts: octa-acid (OA), the closely related tetra-endo-methyl-octa-acid (TEMOA), and cucurbit[8]uril (CB8), along with 21 small organic guest molecules. A total of 119 entries were received from ten participating groups employing a variety of methods that spanned from electronic structure and movable type calculations in implicit solvent to alchemical and potential of mean force strategies using empirical force fields with explicit solvent models. While empirical models tended to obtain better performance than first-principle methods, it was not possible to identify a single approach that consistently provided superior results across all host-guest systems and statistical metrics. Moreover, the accuracy of the methodologies generally displayed a substantial dependence on the system considered, emphasizing the need for host diversity in blind evaluations. Several entries exploited previous experimental measurements of similar host-guest systems in an effort to improve their physical-based predictions via some manner of rudimentary machine learning; while this strategy succeeded in reducing systematic errors, it did not correspond to an improvement in statistical correlation. Comparison to previous rounds of the host-guest binding free energy challenge highlights an overall improvement in the correlation obtained by the affinity predictions for OA and TEMOA systems, but a surprising lack of improvement regarding root mean square error over the past several challenge rounds. The data suggests that further refinement of force field parameters, as well as improved treatment of chemical effects (e.g., buffer salt conditions, protonation states), may be required to further enhance predictive accuracy.

Dataset Information

Bind3P: Optimization of a Water Model Based on Host-Guest Binding Data.

Publications

Bind3P: Optimization of a Water Model Based on Host-Guest Binding Data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure