Critical Importance of Length-Scale Dependence in Implicit Modeling of Hydrophobic Interactions

Abstract

The existence of length-scale dependence of hydrophobic solvation has important implications in the equilibrium of disordered, partially folded, and folded protein conformations. Neglecting this dependence, such as in popular solute surface-area based implicit solvent models with fixed surface tension coefficients, severely limits the ability to accurately model protein conformational equilibrium. We illustrate such fundamental limitations by examining the potentials of mean force of forming dimeric and trimeric nonpolar clusters and propose a new empirical model that effectively captures the context dependence of the local effective surface tension. Further optimization of the new model with other components of the implicit solvent force fields provides promise to significantly improve one's ability to simulate protein folding and conformational transitions. The existence of length-scale dependence of hydrophobic solvation has important implications in the equilibrium of disordered, partially folded, and folded protein conformations. Neglecting this dependence, such as in popular solute surface-area based implicit solvent models with fixed surface tension coefficients, severely limits the ability to accurately model protein conformational equilibrium. We illustrate such fundamental limitations by examining the potentials of mean force of forming dimeric and trimeric nonpolar clusters and propose a new empirical model that effectively captures the context dependence of the local effective surface tension. Further optimization of the new model with other components of the implicit solvent force fields provides promise to significantly improve one's ability to simulate protein folding and conformational transitions.