Boundary element solution of macromolecular electrostatics: interaction energy between two proteins.

The boundary element technique is implemented to solve for the electrostatic potential of macromolecules in an ionic solution. This technique entails solving surface integral equations that are equivalent to the Poisson and the Poisson-Boltzmann equations governing the electrostatic potential inside the macromolecules and and in the solvent. A simple but robust method is described for discretizing the macromolecular surfaces in order to approximate the integral equations by linear algebraic equations. Particular attention is paid to the interaction energy between two macromolecules, and an iterative procedure is devised to make the calculation more efficient. This iterative procedure is illustrated in the electron transfer system of cytochrome c and cytochrome c peroxidase.