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 条
[1]  
[Anonymous], 1980, BOUNDARY ELEMENT TEC
[2]  
Banerjee P.K., 1981, BOUNDARY ELEMENT MET
[3]   ELECTROSTATIC EFFECTS OF CHARGE PERTURBATIONS INTRODUCED BY METAL OXIDATION IN PROTEINS - A THEORETICAL-ANALYSIS [J].
BASHFORD, D ;
KARPLUS, M ;
CANTERS, GW .
JOURNAL OF MOLECULAR BIOLOGY, 1988, 203 (02) :507-510
[4]   Volumes and hydration warmth of ions [J].
Born, M .
ZEITSCHRIFT FUR PHYSIK, 1920, 1 :45-48
[5]  
BOTTA EFF, 1978, THESIS U GRONINGEN N
[6]   CHARMM - A PROGRAM FOR MACROMOLECULAR ENERGY, MINIMIZATION, AND DYNAMICS CALCULATIONS [J].
BROOKS, BR ;
BRUCCOLERI, RE ;
OLAFSON, BD ;
STATES, DJ ;
SWAMINATHAN, S ;
KARPLUS, M .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 1983, 4 (02) :187-217
[7]   MOLECULAR-SURFACE TRIANGULATION [J].
CONNOLLY, ML .
JOURNAL OF APPLIED CRYSTALLOGRAPHY, 1985, 18 (DEC) :499-505
[8]  
DIJKMAN JP, 1989, THESIS U GRONINGEN N, pCH5
[9]  
GILSON MK, 1987, J COMPUT CHEM, V9, P329
[10]  
GILSON MK, 1985, J MOL BIOL, V183, P503