[Paper Review] A Nonlinear Boundary Condition for Continuum Models of Biomolecular Electrostatics

Jaydeep P. Bardhan, D. A. Tejani+3 more|arXiv (Cornell University)|May 24, 2015

Nanopore and Nanochannel Transport Studies34 references3 citations

TL;DR

This paper proposes a nonlinear boundary condition (NLBC) for continuum electrostatic models of biomolecules that captures charge-sign asymmetry—a key physical effect arising from solvent structure and finite solvent size—by modifying the standard Maxwell boundary condition. The NLBC model, validated against all-atom molecular dynamics simulations, accurately reproduces asymmetric electrostatic free energies for surface charges, especially near the solute-solvent interface, and ensures consistency with Gauss’s law in the solvent region.

ABSTRACT

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.

Motivation & Objective

• Address the long-standing limitation of continuum electrostatic models in capturing charge-sign asymmetry near biomolecular surfaces.
• Overcome the symmetric response prediction of standard Poisson and Poisson-Boltzmann theories, which fails to reflect experimental and atomistic simulation observations.
• Develop a physically grounded, nonlinear boundary condition that incorporates the finite size of solvent molecules and interfacial water structure.
• Ensure the model respects Gauss’s law in the solvent region, improving consistency with fundamental electrostatics.
• Enable accurate, computationally efficient predictions of electrostatic free energies for charged biomolecules, particularly for surface-located ions and amino acid side chains.

Proposed method

• Replace the standard Maxwell boundary condition (continuity of normal electric flux) with a nonlinear boundary condition (NLBC) at the solute-solvent interface.
• Formulate the electrostatic problem using boundary-integral equations (BIE) to solve the Poisson equation for complex molecular geometries.
• Model the solvent response as piecewise-affine in the charge magnitude, with a discontinuity at zero charge, reflecting interfacial water structure effects.
• Enforce Gauss’s law in the outer solvent region by modifying the Green’s function to account for screening in dilute electrolytes via the linearized Poisson-Boltzmann equation (LPBE).
• Implement a boundary-element method solver in MATLAB, with publicly available source code and data for reproducibility.
• Use spherical and ellipsoidal geometries as approximations for monatomic ions and amino acids to test the model on benchmark systems.

Experimental results

Research questions

• RQ1 Can a nonlinear boundary condition accurately reproduce the charge-sign asymmetry in biomolecular electrostatics observed in all-atom molecular dynamics simulations?
• RQ2 How does enforcing Gauss’s law in the solvent region affect the accuracy of continuum electrostatic models for charged biomolecules?
• RQ3 To what extent does the magnitude of electrostatic asymmetry depend on the distance of a charge from the solute-solvent interface?
• RQ4 Can the NLBC model quantitatively reproduce charging free energies for monovalent ions and amino acid side chains with high accuracy and computational efficiency?
• RQ5 What is the impact of solvent screening (via Debye-Hückel/linearized Poisson-Boltzmann) on the performance and consistency of the NLBC model?

Key findings

• The new NLBC model successfully reproduces charging free energies for monovalent ions with significantly higher accuracy than the original symmetric model and standard Poisson theory.
• The model captures substantial charge-sign asymmetry for surface charges, with the energetic difference between positive and negative charges remaining large even for protein-sized spheres.
• The magnitude of asymmetry decays rapidly with increasing distance from the solute-solvent interface, consistent with atomistic simulations and physical intuition.
• Enforcement of Gauss’s law in the solvent region improves model consistency and is particularly critical when modeling electrolyte solutions using the linearized Poisson-Boltzmann equation.
• For amino acid side chains, the model predicts that asymmetric electrostatic responses are most significant when the charge is near the protein surface, diminishing with depth.
• The NLBC model is the first asymmetric Poisson-based continuum theory that accurately captures interfacial water effects through a physically motivated boundary condition, validated against high-fidelity all-atom simulations.