Abstract
A general approach has been developed to allow accurate simulations of a small region part of a large macromolecular system while incorporating the influence of the remaining distant atoms with an effective boundary potential. The method is called the Generalized Solvent Boundary Potential (GSBP). By representing the surrounding solvent as a continuum dielectric, both the solvent-shielded static field from the distant atoms of the macromolecule and the reaction field from the dielectric solvent acting on the atoms in the region of interest are included. The static field is calculated once, using the finite-difference Poisson–Boltzmann (PB) equation, and the result is stored on a discrete grid for efficient simulations. The solvent reaction field is developed using a basis-set expansion whose coefficients correspond to generalized electrostatic multipoles. A matrix representing the reaction field Green’s function between those generalized multipoles is calculated only once using the PB equation and stored for efficient simulations. In the present work, the formalism is applied to both spherical and orthorhombic simulation regions for which orthonormal basis-sets exist based on spherical harmonics or cartesian Legendre polynomials. The GSBP method is also tested and illustrated with simple model systems and two detailed atomic systems: the active site region of aspartyl-tRNA synthetase (spherical region) and the interior of the KcsA potassium channel (orthorhombic region). Comparison with numerical finite-difference PB calculations shows that GSBP can accurately describe all long-range electrostatic interactions and remain computationally inexpensive.
Keywords
Spherical harmonicsBoundary value problemHamiltonian (control theory)Solvent modelsPoisson's equationBasis functionPhysicsComputational physicsChemistryComputational chemistryStatistical physicsQuantum mechanicsSolvationMathematicsMolecule
Affiliated Institutions
Related Publications
A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules
Abstract: We present the derivation of a new molecular mechanical force field for simulating the structures, conformational energies, and interaction energies of proteins, nucle...
CHARMM fluctuating charge force field for proteins: II Protein/solvent properties from molecular dynamics simulations using a nonadditive electrostatic model
Abstract A fluctuating charge (FQ) force field is applied to molecular dynamics simulations for six small proteins in explicit polarizable solvent represented by the TIP4P‐FQ po...
<i>A</i> <i>b</i> <i>i</i> <i>n</i> <i>i</i> <i>t</i> <i>i</i> <i>o</i> effective core potentials including relativistic effects. I. Formalism and applications to the Xe and Au atoms
An effective core potential system has been developed for heavy atoms in which relativistic effects are included in the effective potentials (EP). The EP’s are based on numerica...
Dispersion and Absorption in Dielectrics I. Alternating Current Characteristics
The dispersion and absorption of a considerable number of liquid and dielectrics are represented by the empirical formula ε*−ε∞=(ε0−ε∞)/[1+(iωτ0)1−α]. In this equation, ε* is th...
A point‐charge force field for molecular mechanics simulations of proteins based on condensed‐phase quantum mechanical calculations
Abstract Molecular mechanics models have been applied extensively to study the dynamics of proteins and nucleic acids. Here we report the development of a third‐generation point...
Publication Info
- Year
- 2001
- Type
- article
- Volume
- 114
- Issue
- 7
- Pages
- 2924-2937
- Citations
- 226
- Access
- Closed
External Links
Social Impact
Altmetric
PlumX Metrics
Social media, news, blog, policy document mentions
Citation Metrics
226
OpenAlex
Cite This
Wonpil Im,
Simon Bernèche,
Benoı̂t Roux
(2001).
Generalized solvent boundary potential for computer simulations.
The Journal of Chemical Physics
, 114
(7)
, 2924-2937.
https://doi.org/10.1063/1.1336570
Identifiers
- DOI
- 10.1063/1.1336570