Speeding up electrostatic computations for molecular dynamics.
Molecular dynamics (MD) simulations are routinely used to study the structure and function of biological molecules. However the accuracy and duration of these simulations are constrained by their computational costs, thus limiting the ability to accurately simulate systems of realistic sizes over biologically relevant time periods. The two most computationally demanding steps in these simulations are (1) determining the charge state of ionizable sites in biomolecules, which is a key input to the simulation, and (2) calculating long range electrostatic interactions during the simulation. Presented here are two novel methods, the <i>direct interaction approximation (DIA)</i> and the <i>hierarchical charge partitioning (HCP) approximation</i>, for speeding up each of these two computations.
The average charge state of ionizable sites in biomolecules can be calculated as the statistical average over all possible (2N
) microstates for a molecule, where N is the number of ionizable sites. In general this computation scales exponentially as O(N² 2N
). The DIA is an O(²) approximation for calculating the average charge state of ionizable sites. For each site, the DIA treats direct interactions (interactions involving the site of interest) <i>exactly</i>, while using an <i>average</i> value for indirect interactions (interactions not involving the site of interest). The DIA was tested on two problems. The computation of thermal average properties for the 2-D Ising model of ferromagnetism, and the average charge state of ionizable residues in biomolecules. Compared to the commonly used non-deterministic Monte Carlo method, for the same computational cost, the deterministic DIA was found to be at least as accurate, as measured by RMS error relative to the exact computation. Thus, the DIA may be a practical alternative to the Monte Carlo method for some problems.
In atomistic MD simulations, the computation of long range electrostatic interactions, scale as O(<i>n</i>²), where <i>n</i> is the number of atoms. For most biologically relevant timescales the simulations involve 1012–16
simulation steps. Thus, the computational cost of long range interactions become the limiting factor in the size and duration of MD simulations. The HCP is an O(<i>n</i> log <i>n</i>) approximation for computing long range electrostatic interactions. The approximation is based on multiple levels of natural partitioning of biomolecular structures into a hierarchical set of components. For components that are far from the point of interest, the charge distribution for each component is approximated by a much smaller number of charges. For nearby components, the HCP uses the full set of atomic charges. For large structures the HCP can be several orders of magnitude faster than the exact pairwise O(<i>n</i>²) all-atom computation. For a representative set of structures, the accuracy of the HCP is comparable to the industry standard explicit solvent particle mesh Ewald (PME), and is in general more accurate than the…
Advisors/Committee Members: Onufriev, Alexey V. (committeechair), Bevan, David R. (committee member), Cao, Yang (committee member), Sandu, Adrian (committee member), Xing, Jianhua (committee member).
to Zotero / EndNote / Reference
APA (6th Edition):
Anandakrishnan, R. (2011). Speeding up electrostatic computations for molecular dynamics. (Doctoral Dissertation). Virginia Tech. Retrieved from http://hdl.handle.net/10919/40262
Chicago Manual of Style (16th Edition):
Anandakrishnan, Ramamoorthi. “Speeding up electrostatic computations for molecular dynamics.” 2011. Doctoral Dissertation, Virginia Tech. Accessed March 04, 2021.
MLA Handbook (7th Edition):
Anandakrishnan, Ramamoorthi. “Speeding up electrostatic computations for molecular dynamics.” 2011. Web. 04 Mar 2021.
Anandakrishnan R. Speeding up electrostatic computations for molecular dynamics. [Internet] [Doctoral dissertation]. Virginia Tech; 2011. [cited 2021 Mar 04].
Available from: http://hdl.handle.net/10919/40262.
Council of Science Editors:
Anandakrishnan R. Speeding up electrostatic computations for molecular dynamics. [Doctoral Dissertation]. Virginia Tech; 2011. Available from: http://hdl.handle.net/10919/40262