Implementation and testing of stable, fast implicit solvation in molecular dynamics using the smooth-permittivity
Ninad V Prabhu1, Peijuan Zhu, Kim A Sharp
1Johnson Research Foundation and Department of Biochemistry and Biophysics, University of Pennsylvania, 37th and Hamilton Walk, Philadelphia, Pennsylvania 19104, USA.
Abstract:
A fast stable finite difference Poisson-Boltzmann (FDPB) model for implicit solvation in molecular dynamics simulations was developed using the smooth permittivity FDPB method implemented in the OpenEye ZAP libraries. This was interfaced with two widely used molecular dynamics packages, AMBER and CHARMM. Using the CHARMM-ZAP software combination, the implicit solvent model was tested on eight proteins differing in size, structure, and cofactors: calmodulin, horseradish peroxidase (with and without substrate analogue bound), lipid carrier protein, flavodoxin, ubiquitin, cytochrome c, and a de novo designed 3-helix bundle. The stability and accuracy of the implicit solvent simulations was assessed by examining root-mean-squared deviations from crystal structure. This measure was compared with that of a standard explicit water solvent model. In addition we compared experimental and calculated NMR order parameters to obtain a residue level assessment of the accuracy of MD-ZAP for simulating dynamic quantities. Overall, the agreement of the implicit solvent model with experiment was as good as that of explicit water simulations. The implicit solvent method was up to eight times faster than the explicit water simulations, and approximately four times slower than a vacuum simulation (i.e., with no solvent treatment).
More Related Videos
05:37Rapid in-silico Battery Electrolyte Electrochemical Reaction Generation using 3T-VASP Multi-Scale Energy Minimization
Published on: August 22, 2025
10:52Multiscale Sampling of a Heterogeneous Water/Metal Catalyst Interface using Density Functional Theory and Force-Field Molecular Dynamics
Published on: April 12, 2019
Related Concept Videos
Kohlraush’s Law and its Applications
Maxwell-Boltzmann Distribution: Problem Solving
This distribution function f(v) is defined by saying that the expected number N (v1,v2) of particles with speeds between v1 and v2 is given by
Theories of Dissolution: The Danckwerts' Model and Interfacial Barrier Model
Electrostatic Boundary Conditions in Dielectrics
Consider a case where both the mediums across a boundary are two different dielectric materials. Recall that the electric field and electric displacement are proportional and related through the material's permittivity.
The Debye–Hückel Theory of Electrolyte Solutions
Molecular and Ionic Solids
Molecular Solids
Molecular crystalline solids, such as ice, sucrose (table sugar), and iodine, are solids that are composed of neutral molecules as their constituent units. These molecules are held together by weak intermolecular forces such as London dispersion forces, dipole-dipole interactions, or hydrogen bonds, which...
