Abstract:The simulation of the interactions of proteins with charged surfaces in a condensed-phase aqueous solution containing electrolytes using empirical force field based methods is predominantly governed by nonbonded interactions between the atoms of the protein, surface, and the solvent. Electrostatic effects represent the strongest type of these interactions and the type that is most difficult to accurately represent because of their long-range influence. While many different methods have been developed to represent electrostatic interactions, the particle mesh Ewald summation (PME) method is generally considered to be the most accurate one for calculating these effects. However, the PME method was designed for systems with three-dimensional (3D) periodicity, and not for interfacial systems such as the case of protein adsorption to a charged surface. Interfacial systems such as these have only two-dimensional periodicity, which may not be appropriate for treatment with PME due to the possibility that the presence of multiple charged image surfaces parallel to the primary simulation cell's surface, may introduce nonphysical effects on the behavior of the charged molecules in the system. In an effort to address this issue, the authors have conducted a set of nanosecond-scale molecular dynamics simulations to calculate the equilibrium distribution of Na(+) and Cl(-) ions near a charged surface using PME and a range of radial cutoff methods for treating electrostatic interactions, where the cutoffs prevent interaction with the periodic images of the system. The resulting ion concentration profiles were compared to one another and to a continuum analytical solution of the theoretical ion distribution obtained from the Poisson-Boltzmann equation. Their results show that the PME method does not introduce the suspected nonphysical effects in the ion distributions due to the 3D periodic images of the system, thus indicating that it is appropriate for use for this type of molecular simulation. Although their interest is motivated by protein-surface interactions, the conclusions are applicable for the treatment of electrostatics in other aqueous systems with two-dimensional periodicity.

Electrostatic interaction of a solute with a continuum. A direct utilizaion of AB initio molecular potentials for the prevision of solvent effects

Perspective on “Electrostatic interactions of a solute with a continuum. A direct utilization of ab initio molecular potentials for the prevision of solvent effects”

Electrostatic interactions in water: nonlocal electrostatic approach

Influence of cut-off truncation and artificial periodicity of electrostatic interactions in molecular simulations of solvated ions: A continuum electrostatics study

A quantum mechanical computational method for modeling electrostatic and solvation effects of protein

Electrostatics of solvated systems in periodic boundary conditions

Electrostatic solvation free energies of charged hard spheres using molecular dynamics with density functional theory interactions

Development of molecular simulation methods to accurately represent protein-surface interactions: Method assessment for the calculation of electrostatic effects

Accurate boundary-integral formulations for the calculation of electrostatic forces with an implicit-solvent model

Molecular Theories and Simulation of Ions and Polar Molecules in Water

Some challenges of diffused interfaces in implicit-solvent models

Tackling solvent effect by coupling electronic and molecular Density Functional Theory

Implicit self-consistent electrolyte model in plane-wave density-functional theory

Dielectric continuum methods for quantum chemistry

Electrostatic interactions in charged nanoslits within an explicit solvent theory

Implicit Solvation Methods for Catalysis at Electrified Interfaces

Linear Response Properties of Solvated Systems: A Computational Study

Dielectric response with short-ranged electrostatics

Entropic Origin of Ionic Interactions in Polar Solvents

The Poisson-Boltzmann model for implicit solvation of electrolyte solutions: Quantum chemical implementation and assessment via Sechenov coefficients

Accounting for electronic polarization in nonpolarizable force fields