Application of the level-set method to the implicit solvation of nonpolar molecules

Abstract

A level-set method is developed for numerically capturing the equilibrium solute-solvent interface that is defined by the recently proposed variational implicit solvent model (Dzubiella, Swanson, and McCammon, Phys. Rev. Lett. {\bf 104}, 527 (2006) and J. Chem.\Phys. {\bf 124}, 084905 (2006)). In the level-set method, a possible solute-solvent interface is represented by the zero level-set (i.e., the zero level surface) of a level-set function and is eventually evolved into the equilibrium solute-solvent interface. The evolution law is determined by minimization of a solvation free energy {\it functional} that couples both the interfacial energy and the van der Waals type solute-solvent interaction energy. The surface evolution is thus an energy minimizing process, and the equilibrium solute-solvent interface is an output of this process. The method is implemented and applied to the solvation of nonpolar molecules such as two xenon atoms, two parallel paraffin plates, helical alkane chains, and a single fullerene $C_{60}$. The level-set solutions show good agreement for the solvation energies when compared to available molecular dynamics simulations. In particular, the method captures solvent dewetting (nanobubble formation) and quantitatively describes the interaction in the strongly hydrophobic plate system.

References in corpus (5)

• Temperature Dependence of the Hydrophobic Hydration and Interaction of Simple Solutes: An Examination of Five Popular Water Models
• Morphological Thermodynamics of Fluids: Shape Dependence of Free Energies
• Coupling hydrophobic, dispersion, and electrostatic contributions in continuum solvent models
• Coupling nonpolar and polar solvation free energies in implicit solvent models
• Interface dynamics of microscopic cavities in water

Cited by in corpus (7)

• Dewetting-controlled binding of ligands to hydrophobic pockets
• Variational implicit-solvent predictions of the dry-wet transition pathways for ligand-receptor binding and unbinding kinetics
• Parameter optimization in differential geometry based solvation models
• The Calculus of Boundary Variations and the Dielectric Boundary Force in the Poisson-Boltzmann Theory for Molecular Solvation
• Exploring A Multi-Scale Method for Molecular Simulations in Continuum Solvent Model: Explicit Simulation of Continuum Solvent As An Incompressible Fluid
• Binary Level Set Method for Variational Implicit Solvation Model
• Explicit-Solute Implicit-Solvent Molecular Simulation with Binary Level-Set, Adaptive-Mobility, and GPU