Dataset Information

Order N algorithm for computation of electrostatic interactions in biomolecular systems.

ABSTRACT: Poisson-Boltzmann electrostatics is a well established model in biophysics; however, its application to large-scale biomolecular processes such as protein-protein encounter is still limited by the efficiency and memory constraints of existing numerical techniques. In this article, we present an efficient and accurate scheme that incorporates recently developed numerical techniques to enhance our computational ability. In particular, a boundary integral equation approach is applied to discretize the linearized Poisson-Boltzmann equation; the resulting integral formulas are well conditioned and are extended to systems with arbitrary numbers of biomolecules. The solution process is accelerated by Krylov subspace methods and a new version of the fast multipole method. In addition to the electrostatic energy, fast calculations of the forces and torques are made possible by using an interpolation procedure. Numerical experiments show that the implemented algorithm is asymptotically optimal O(N) in both CPU time and required memory, and application to the acetylcholinesterase-fasciculin complex is illustrated.

SUBMITTER: Lu B

PROVIDER: S-EPMC1748223 | biostudies-literature | 2006 Dec

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Order N algorithm for computation of electrostatic interactions in biomolecular systems.

Lu Benzhuo B Cheng Xiaolin X Huang Jingfang J McCammon J Andrew JA

Proceedings of the National Academy of Sciences of the United States of America 20061205 51

Poisson-Boltzmann electrostatics is a well established model in biophysics; however, its application to large-scale biomolecular processes such as protein-protein encounter is still limited by the efficiency and memory constraints of existing numerical techniques. In this article, we present an efficient and accurate scheme that incorporates recently developed numerical techniques to enhance our computational ability. In particular, a boundary integral equation approach is applied to discretize ...[more]

PMID: 17148613

Similar Datasets

Project description:An ability to efficiently compute the electrostatic potential produced by molecular charge distributions under realistic solvation conditions is essential for a variety of applications. Here, the simple closed-form analytical approximation to the Poisson equation rigorously derived in Part I for idealized spherical geometry is tested on realistic shapes. The effects of mobile ions are included at the Debye-Huckel level. The accuracy of the resulting closed-form expressions for electrostatic potential is assessed through comparisons with numerical Poisson-Boltzmann (NPB) reference solutions on a test set of 580 representative biomolecular structures under typical conditions of aqueous solvation. For each structure, the deviation from the reference is computed for a large number of test points placed near the dielectric boundary (molecular surface). The accuracy of the approximation, averaged over all test points in each structure, is within 0.6 kcal/mol/mid R:emid R: approximately kT per unit charge for all structures in the test set. For 91.5% of the individual test points, the deviation from the NPB potential is within 0.6 kcal/mol/mid R:emid R:. The deviations from the reference decrease with increasing distance from the dielectric boundary: The approximation is asymptotically exact far away from the source charges. Deviation of the overall shape of a structure from ideal spherical does not, by itself, appear to necessitate decreased accuracy of the approximation. The largest deviations from the NPB reference are found inside very deep and narrow indentations that occur on the dielectric boundaries of some structures. The dimensions of these pockets of locally highly negative curvature are comparable to the size of a water molecule; the applicability of a continuum dielectric models in these regions is discussed. The maximum deviations from the NPB are reduced substantially when the boundary is smoothed by using a larger probe radius (3 A) to generate the molecular surface. A detailed accuracy analysis is presented for several proteins of various shapes, including lysozyme whose surface features a functionally relevant region of negative curvature. The proposed analytical model is computationally inexpensive; this strength of the approach is demonstrated by computing and analyzing the electrostatic potential generated by a full capsid of the tobacco ring spot virus at atomic resolution (500 000 atoms). An analysis of the electrostatic potential of the inner surface of the capsid reveals what might be a RNA binding pocket. These results are generated with the modest computational power of a desktop personal computer.

Dataset Information

Order N algorithm for computation of electrostatic interactions in biomolecular systems.

Publications

Order N algorithm for computation of electrostatic interactions in biomolecular systems.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets