A Nonlinear Boundary Condition for Continuum Models of Biomolecular Electrostatics

Understanding the behavior of biomolecules such as proteins requires understanding the critical influence of the surrounding fluid (solvent) environment--water with mobile salt ions such as sodium. Unfortunately, for many studies, fully atomistic simulations of biomolecules, surrounded by thousands...

Full description

Saved in:
Bibliographic Details
Published inarXiv.org
Main Authors Bardhan, J P, Tejani, D A, Wieckowski, N S, Ramaswamy, A, Knepley, M G
Format Paper
LanguageEnglish
Published Ithaca Cornell University Library, arXiv.org 24.05.2015
Subjects
Biomolecules
Boltzmann transport equation
Boundary conditions
Computer simulation
Continuum modeling
Debye-Huckel theory
Electrostatics
Mathematical models
Molecular dynamics
Poisson equation
Proteins
Solvents
Water chemistry
Online AccessGet full text

Cover

More Information
Summary:Understanding the behavior of biomolecules such as proteins requires understanding the critical influence of the surrounding fluid (solvent) environment--water with mobile salt ions such as sodium. Unfortunately, for many studies, fully atomistic simulations of biomolecules, surrounded by thousands of water molecules and ions are too computationally slow. Continuum solvent models based on macroscopic dielectric theory (e.g. the Poisson equation) are popular alternatives, but their simplicity fails to capture well-known phenomena of functional significance. For example, standard theories predict that electrostatic response is symmetric with respect to the sign of an atomic charge, even though response is in fact strongly asymmetric if the charge is near the biomolecule surface. In this work, we present an asymmetric continuum theory that captures the essential physical mechanism--the finite size of solvent atoms--using a nonlinear boundary condition (NLBC) at the dielectric interface between the biomolecule and solvent. Numerical calculations using boundary-integral methods demonstrate that the new NLBC model reproduces a wide range of results computed by more realistic, and expensive, all-atom molecular-dynamics (MD) simulations in explicit water. We discuss model extensions such as modeling dilute-electrolyte solvents with Debye-Huckel theory (the linearized Poisson-Boltzmann equation) and opportunities for the electromagnetics community to contribute to research in this important area of molecular nanoscience and engineering.
ISSN:2331-8422