THE ELECTRIC-POTENTIAL OF A MACROMOLECULE IN A SOLVENT - A FUNDAMENTAL APPROACH

被引:290
作者
JUFFER, AH [1 ]
BOTTA, EFF [1 ]
VANKEULEN, BAM [1 ]
VANDERPLOEG, A [1 ]
BERENDSEN, HJC [1 ]
机构
[1] STATE UNIV GRONINGEN,DEPT MATH,9700 AB GRONINGEN,NETHERLANDS
关键词
D O I
10.1016/0021-9991(91)90043-K
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
A general numerical method is presented to compute the electric potential for a macromolecule of arbitrary shape in a solvent with nonzero ionic strength. The model is based on a continuum description of the dielectric and screening properties of the system, which consists of a bounded internal region with discrete charges and an infinite external region. The potential obeys the Poisson equation in the internal region and the linearized Poisson-Boltzmann equation in the external region, coupled through appropriate boundary conditions. It is shown how this three-dimensional problem can be presented as a pair of coupled integral equations for the potential and the normal component of the electric field at the dielectric interface. These equations can be solved by a straightforward application of boundary element techniques. The solution involves the decomposition of a matrix that depends only on the geometry of the surface and not on the positions of the charges. With this approach the number of unknowns is reduced by an order of magnitude with respect to the usual finite difference methods. Special attention is given to the numerical inaccuracies resulting from charges which are located close to the interface; an adapted formulation is given for that case. The method is tested both for a spherical geometry, for which an exact solution is available, and for a realistic problem, for which a finite difference solution and experimental verification is available. The latter concerns the shift in acid strength (pK-values) of histidines in the copper-containing protein azurin on oxidation of the copper, for various values of the ionic strength. A general method is given to triangulate a macromolecular surface. The possibility is discussed to use the method presented here for a correct treatment of long-range electrostatic interactions in simulations of solvated macromolecules, which form an essential part of correct potentials of mean force. © 1991.
引用
收藏
页码:144 / 171
页数:28
相关论文
共 34 条
[21]  
STACKGOLD I, 1968, MACMILLAN SERIES ADV, V2
[23]   THEORY OF PROTEIN TITRATION CURVES .1. GENERAL EQUATIONS FOR IMPENETTRABLE SPHERES [J].
TANFORD, C ;
KIRKWOOD, JG .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1957, 79 (20) :5333-5339
[25]  
VANDEREERDEN JP, 1983, NUCLEI ACID CONFORMA, P61
[26]  
VONGUNSTEREN WF, 1987, GRONINGEN MOL SIMULA
[27]   CALCULATIONS OF ELECTROSTATIC INTERACTIONS IN BIOLOGICAL-SYSTEMS AND IN SOLUTIONS [J].
WARSHEL, A ;
RUSSELL, ST .
QUARTERLY REVIEWS OF BIOPHYSICS, 1984, 17 (03) :283-422
[28]   CALCULATION OF THE ELECTRIC-POTENTIAL IN THE ACTIVE-SITE CLEFT DUE TO ALPHA-HELIX DIPOLES [J].
WARWICKER, J ;
WATSON, HC .
JOURNAL OF MOLECULAR BIOLOGY, 1982, 157 (04) :671-679
[29]  
Westneimer FH, 1938, J CHEM PHYS, V6, P513, DOI 10.1063/1.1750303
[30]   THE RIGOROUS COMPUTATION OF THE MOLECULAR ELECTRIC-POTENTIAL [J].
ZAUHAR, RJ ;
MORGAN, RS .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 1988, 9 (02) :171-187