Next: 3-Phase Calculation
Up: Generalized Born Implicit Solvent
Previous: Generalized Born Implicit Solvent
Contents
Index
Subsections
Water has many biologically necessary properties, one of which is as a dielectric.
As a dielectric, water screens (lessens) electrostatic interactions between charged particles.
Water can therefore be crudely modeled as a dielectric continuum.
In this manner, the electrostatic forces of a biological system can be expressed as a system of differential equations which can be solved for the electric field caused by a collection of charges.
The Poisson Boltzmann equation (PBE),
is a nonlinear equation which solves for the electrostatic field,
, based on
the position dependent dielectric,
,
the position-dependent accessibility of position
to the ions in solution,
,
the solute charge distribution,
,
and the bulk charge density,
, of ion
.
While this equation does exactly solve for the electrostic field of a charge distribution in a dielectric, it is very expensive to solve, and therefore not suitable for molecular dynamics.
The Generalized Born (GB) equation is an approximation of the PBE.
It models atoms as charged spheres whose internal dielectric is lower than that of the environment.
The screening which each atom,
, experiences is determined by the local environment;
the more atom
is surrounded by other atoms, the less it's electrostatics will be screened since it is more surrounded by low dielectric; this property is called one atom descreening another.
Different GB models calculate atomic descreening differently.
Descreening is used to calculate the Born radius,
, of each atom.
The Born radius of an atom measures the degree of descreening.
A large Born radius represents small screening (strong electric field) as if the atom were in vacuum.
A small Born radius represents large screening (weak electric field) as if the atom were in bulk water.
The next section describes how the Born radius is calculated and how it is used to calculate electrostatics.
In a GB simulation, the total electrostatic force on an atom,
, is the net Coulomb force on atom
(from nearby atoms) minus the GB force on atom
(also caused by nearby atoms):
Forces are contributed by other nearby atoms within a cutoff.
The GB force on atom
is the derivative of the total GB energy with respect to relative atom distances
,
where the partial derivatives are included since the Born radius,
, is a function of all relative atom distances.
The total GB energy of the system is
|
(11) |
where
is the Born radius dependent self energy of atom
, and the GB energy between atoms
and
is given by
|
(12) |
The dielectric term [69] is
|
(13) |
and the GB function [72] is
|
(14) |
As the Born radii of atoms
and
decrease (increasing screening), the effective distance between the atoms (
) increases.
The implicit solvent implemented in NAMD is the model of Onufriev, Bashford and Case [53,54] which calculates the Born radius as
|
(15) |
where
|
(16) |
is the piecewise descreening function [54,28,61]; the seven piecewise regimes are
|
(17) |
and the values of
are
|
(18) |
Below are defined the derivatives of the above functions which are required for force calculations.
|
(19) |
|
(20) |
|
(21) |
|
(22) |
|
(27) |
|
(28) |
|
(29) |
Other variables referenced in the above GB equations are
-
- distance between atoms i and j; calculated from atom coordinates.
-
- debye screening length; calculated from ion concentration,
;
Å for 0.1 M monovalent salt.
-
- dielectric constant of solvent.
-
- dielectric constant of protein.
-
- Born radius of atom
.
-
- intrinsic radius of atom
taken from Bondi [9].
-
- intrinsic radius offset;
Å by default [54].
-
-
-
- atom radius scaling factor [28,69].
-
- Coulomb's constant,
, 332.063711 kcal Å / e
.
-
[54]
Next: 3-Phase Calculation
Up: Generalized Born Implicit Solvent
Previous: Generalized Born Implicit Solvent
Contents
Index
http://www.ks.uiuc.edu/Research/namd/