Abstract:The free energetics of water density fluctuations near a surface, and the rare low-density fluctuations in particular, serve as reliable indicators of surface hydrophobicity; the easier it is to displace the interfacial waters, the more hydrophobic the underlying surface. However, characterizing the free energetics of such rare fluctuations requires computationally expensive, non-Boltzmann sampling methods like umbrella sampling. This inherent computational expense associated with umbrella sampling makes it challenging to investigate the role of polarizability or electronic structure effects in influencing interfacial fluctuations. Importantly, it also limits the size of the volume, which can be used to probe interfacial fluctuations. The latter can be particularly important in characterizing the hydrophobicity of large surfaces with molecular-level heterogeneities, such as those presented by proteins. To overcome these challenges, here we present a method for the sparse sampling of water density fluctuations, which is roughly two orders of magnitude more efficient than umbrella sampling. We employ thermodynamic integration to estimate the free energy differences between biased ensembles, thereby circumventing the umbrella sampling requirement of overlap between adjacent biased distributions. Further, a judicious choice of the biasing potential allows such free energy differences to be estimated using short simulations, so that the free energetics of water density fluctuations are obtained using only a few, short simulations. Leveraging the efficiency of the method, we characterize water density fluctuations in the entire hydration shell of the protein, ubiquitin; a large volume containing an average of more than six hundred waters.

Successive Umbrella Sampling

Umbrella sampling: a powerful method to sample tails of distributions

Eigenvector method for umbrella sampling enables error analysis

Buckets Instead of Umbrellas for Enhanced Sampling and Free Energy Calculations

An Efficient Umbrella Potential for the Accurate Calculation of Free Energies by Molecular Simulation.

Variational Umbrella Seeding for Calculating Nucleation Barriers

Stratified sampling and bootstrapping for approximate Bayesian computation

Self-Healing Umbrella Sampling: Convergence and efficiency

Multi-point nonequilibrium umbrella sampling and associated fluctuation relations

A Unified Approach to Enhanced Sampling

Free energy surface reconstruction from umbrella samples using Gaussian process regression II: Multiple collective variables

Superposition Enhanced Nested Sampling

Extension to the weighted histogram analysis method: combining umbrella sampling with free energy calculations

Adaptive Multilevel Subset Simulation with Selective Refinement

Sparse Sampling of Water Density Fluctuations in Interfacial Environments

Patchwork sampling of stochastic differential equations

Sequential Sampling for Optimal Weighted Least Squares Approximations in Hierarchical Spaces

Ensemble Rejection Sampling

Sampling free energy surfaces as slices by combining umbrella sampling and metadynamics

Snowballing Nested Sampling