Dissertation/Thesis Abstract

Fast direct methods for molecular electrostatics
by Ho, Kenneth L., Ph.D., New York University, 2012, 180; 3524158
Abstract (Summary)

Electrostatic interactions are vital for many aspects of biomolecular structure and function, including stability, activity, and specificity. Thus, a central problem in biophysical modeling is the electrostatic analysis of biomolecular systems. Here, we consider the numerical solution of the linearized Poisson-Boltzman equation (LPBE), a widely used model for the electrostatic potential of a large, solvated biomolecule. By combining boundary integral techniques with new multilevel matrix compression algorithms, we develop a fast direct solver for the LPBE that is accurate, robust, and can be more efficient than current methods by several orders of magnitude.

The fast direct solver is general and applies to a wide range of integral operators based on non-oscillatory Green's functions, including those for the Laplace, low-frequency Helmholtz, Stokes, and LPBEs. The core algorithm uses the interpolative decomposition to compress the matrix discretizations of such operators, producing highly efficient representations that facilitate fast inversion. For boundary integral equations in 2D, the solver has complexity O(N), where N is the number of discretization elements; in 3D, it incurs an O(N 3/2) cost for precomputation followed by O(N log N) solves. As is typical of direct methods, each solve can be performed extremely rapidly, though the cost of precomputation can be high. Thus, the solver is particularly suited to problems where the precomputation time can be amortized, e.g., systems with ill-conditioned matrices or involving multiple right-hand sides.

We demonstrate our solver on a number of examples and discuss various useful extensions. Furthermore, we apply our methods to the calculation of protein pKa values, which requires the computation of all pairwise titrating site energies. This corresponds to solving the LPBE on the same molecular geometry with many different boundary conditions on the protein surface, each manifesting as a different right-hand side, and hence presents a prime candidate for acceleration using our direct solver. Preliminary results are favorable and show the viability of our techniques for molecular electrostatics.

Such fast direct methods could well have broad impact on many areas of computational science and engineering. We describe further applications in biology, chemistry, and physics, and outline some directions for future work.

Indexing (document details)
Advisor: Greengard, Leslie
Commitee: Bonneau, Richard A., Peskin, Charles S., Tranchina, Daniel, Zhang, Yingkai
School: New York University
Department: Program in Computational Biology
School Location: United States -- New York
Source: DAI-B 74/01(E), Dissertation Abstracts International
Source Type: DISSERTATION
Subjects: Applied Mathematics, Biophysics, Computer science
Keywords: Fast direct solvers, Integral equations, Molecular electrostatics, Multilevel matrix compression, Poisson-boltzmann equation, Potential theory, Protein pka
Publication Number: 3524158
ISBN: 978-1-267-58399-4
Copyright © 2019 ProQuest LLC. All rights reserved. Terms and Conditions Privacy Policy Cookie Policy
ProQuest