Show simple item record

dc.contributor.authorAnandakrishnan, Ramamoorthien_US
dc.description.abstractMolecular 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 direct interaction approximation (DIA) and the hierarchical charge partitioning (HCP) approximation, 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 (2^N) microstates for a molecule, where N is the number of ionizable sites. In general this computation scales exponentially as O(N^2 2^N). The DIA is an O(N^2) approximation for calculating the average charge state of ionizable sites. For each site, the DIA treats direct interactions (interactions involving the site of interest) exactly, while using an average 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(n^2), where n is the number of atoms. For most biologically relevant timescales the simulations involve 10^12â 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(n log n) 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(n^2) 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 spherical cutoff method. And, unlike the PME, the DIA can be easily extended to implicit solvent GB models. 50 ns implicit solvent simulations for a representative set of four biomolecules suggests that the HCP could be a practical alternative for implicit solvent simulations, and preferable to the cutoff based method. The HCP is available for general use in the open source MD software, NAB within AmberTools.

dc.publisherVirginia Techen_US
dc.rightsI hereby certify that, if appropriate, I have obtained and attached hereto a written permission statement from the owner(s) of each third party copyrighted matter to be included in my thesis, dissertation, or project report, allowing distribution as specified below. I certify that the version I submitted is the same as that approved by my advisory committee. I hereby grant to Virginia Tech or its agents the non-exclusive license to archive and make accessible, under the conditions specified below, my thesis, dissertation, or project report in whole or in part in all forms of media, now or hereafter known. I retain all other ownership rights to the copyright of the thesis, dissertation or project report. I also retain the right to use in future works (such as articles or books) all or part of this thesis, dissertation, or project report.en_US
dc.subjectstatistical mechanicsen_US
dc.subjectbiomolecular electrostaticsen_US
dc.subjectmolecular dynamicsen_US
dc.titleSpeeding up electrostatic computations for molecular dynamicsen_US
dc.contributor.departmentComputer Science and Applicationsen_US
dc.description.degreePh. D.en_US D.en_US Polytechnic Institute and State Universityen_US Science and Applicationsen_US
dc.contributor.committeechairOnufriev, Alexey V.en_US
dc.contributor.committeememberBevan, David R.en_US
dc.contributor.committeememberCao, Yangen_US
dc.contributor.committeememberSandu, Adrianen_US
dc.contributor.committeememberXing, Jianhuaen_US

Files in this item


This item appears in the following Collection(s)

Show simple item record