The AGBNP2 Implicit Solvation Model

TL;DR

The AGBNP2 model enhances implicit solvation modeling by incorporating hydration site scoring and a solvent excluding volume, improving hydration free energy predictions and conformational sampling accuracy over previous models.

Contribution

It introduces a new empirical hydration site component and an analytical solvent excluding volume model to better capture hydration effects in implicit solvent simulations.

Findings

Improved agreement with experimental hydration free energies.

Better conformational ensemble accuracy in mini-protein simulations.

Enhanced modeling of first solvation shell effects.

Abstract

The AGBNP2 implicit solvent model, an evolution of the Analytical Generalized Born plus Non-Polar (AGBNP) model we have previously reported, is presented with the aim of modeling hydration effects beyond those described by conventional continuum dielectric representations. A new empirical hydration free energy component based on a procedure to locate and score hydration sites on the solute surface is introduced to model first solvation shell effects, such as hydrogen bonding, which are poorly described by continuum dielectric models. This new component is added to the Generalized Born and non-polar AGBNP models which have been improved with respect to the description of the solute volume description. We have introduced an analytical Solvent Excluding Volume (SEV) model which reduces the effect of spurious high-dielectric interstitial spaces present in conventional van der Waals representations of the solute volume. The AGBNP2 model is parametrized and tested with respect to experimental hydration free energies of small molecules and the results of explicit solvent simulations. Modeling the granularity of water is one of the main design principles employed for the the first shell solvation function and the SEV model, by requiring that water locations have a minimum available volume based on the size of a water molecule. We show that the new volumetric model produces Born radii and surface areas in good agreement with accurate numerical evaluations. The results of Molecular Dynamics simulations of a series of mini-proteins show that the new model produces conformational ensembles in substantially better agreement with reference explicit solvent ensembles than the original AGBNP model with respect to both structural and energetics measures.