Fast Methods for Simulation of Biomolecule of Electrostatics

Biomolecular structure and interactions in aqueous environment are determined by a complicated interplay between physical and chemical forces including solvation, electrostatics, van der Waals forces, the hydrophobic effect and covalent bonding. Among them, electrostatics has been of particular inte...

Full description

Bibliographic Details
Main Authors: Kuo, Shihhsien, Altman, Michael D., Bardhan, Jaydeep P., Tidor, Bruce, White, Jacob K.
Format: Article
Language:en_US
Published: 2003
Subjects:
Online Access:http://hdl.handle.net/1721.1/3717
_version_ 1810982335660687360
author Kuo, Shihhsien
Altman, Michael D.
Bardhan, Jaydeep P.
Tidor, Bruce
White, Jacob K.
author_facet Kuo, Shihhsien
Altman, Michael D.
Bardhan, Jaydeep P.
Tidor, Bruce
White, Jacob K.
author_sort Kuo, Shihhsien
collection MIT
description Biomolecular structure and interactions in aqueous environment are determined by a complicated interplay between physical and chemical forces including solvation, electrostatics, van der Waals forces, the hydrophobic effect and covalent bonding. Among them, electrostatics has been of particular interest due to its long-range nature and the tradeoff between desolvation and interaction effects [1]. In addition, electrostatic interactions play a significant role within a biomolecule as well as between biomolecules, making the balance between the two vital to the understanding of macromolecular systems. As a result, much effort has been devoted to accurate modeling and simulation of biomolecule electrostatics. One important application of this work is to compute the structure of electrostatic interactions for a biomolecule in an electrolyte solution, as well as the potential that the molecule generates in space. There are two valuable uses for these simulations. First, it provides a full picture of the electrostatic energetics of a biomolecular system, improving our understanding of how electrostatics contributes to stability, specificity, function, and molecular interaction [2]. Second, these simulations serve as a tool for molecular design, since electrostatic complementarity is an important feature of interacting molecules. Through examination of the electrostatics and potential field generated by a protein molecule, for example, it may be possible to suggest improvements to other proteins or drug molecules that interact with it, or perhaps even design new interacting molecules de novo [3]. There are two approaches in simulating a protein macromolecule in an aqueous solution with nonzero ionic strength. Discrete/atomistic approaches based on Monte-Carlo or molecular dynamics simulations treat the macromolecule and solvent explicitly at the atomic level. Therefore, an enormous number of solvent molecules are required to provide reasonable accuracy, especially when electric fields far away from macroscopic surface are of interest, leading to computational infeasibility. In this work, we adopt instead an approach based on a continuum description of the macromolecule and solvent. Although the continuum model of biomolecule electrostatics is widely used, the numerical techniques used to evaluate the model do not exploit fast solver approaches developed for analyzing integrated circuit interconnect. I will describe the formulation used for analyzing biomolecule electrostatics, and then derive an integral formulation of the problem that can be rapidly solved with precorrected-FFT method [4].
first_indexed 2024-09-23T10:29:00Z
format Article
id mit-1721.1/3717
institution Massachusetts Institute of Technology
language en_US
last_indexed 2024-09-23T10:29:00Z
publishDate 2003
record_format dspace
spelling mit-1721.1/37172019-04-12T08:09:11Z Fast Methods for Simulation of Biomolecule of Electrostatics Kuo, Shihhsien Altman, Michael D. Bardhan, Jaydeep P. Tidor, Bruce White, Jacob K. biomolecule electrostatics continuum model precorrected-FFT method Biomolecular structure and interactions in aqueous environment are determined by a complicated interplay between physical and chemical forces including solvation, electrostatics, van der Waals forces, the hydrophobic effect and covalent bonding. Among them, electrostatics has been of particular interest due to its long-range nature and the tradeoff between desolvation and interaction effects [1]. In addition, electrostatic interactions play a significant role within a biomolecule as well as between biomolecules, making the balance between the two vital to the understanding of macromolecular systems. As a result, much effort has been devoted to accurate modeling and simulation of biomolecule electrostatics. One important application of this work is to compute the structure of electrostatic interactions for a biomolecule in an electrolyte solution, as well as the potential that the molecule generates in space. There are two valuable uses for these simulations. First, it provides a full picture of the electrostatic energetics of a biomolecular system, improving our understanding of how electrostatics contributes to stability, specificity, function, and molecular interaction [2]. Second, these simulations serve as a tool for molecular design, since electrostatic complementarity is an important feature of interacting molecules. Through examination of the electrostatics and potential field generated by a protein molecule, for example, it may be possible to suggest improvements to other proteins or drug molecules that interact with it, or perhaps even design new interacting molecules de novo [3]. There are two approaches in simulating a protein macromolecule in an aqueous solution with nonzero ionic strength. Discrete/atomistic approaches based on Monte-Carlo or molecular dynamics simulations treat the macromolecule and solvent explicitly at the atomic level. Therefore, an enormous number of solvent molecules are required to provide reasonable accuracy, especially when electric fields far away from macroscopic surface are of interest, leading to computational infeasibility. In this work, we adopt instead an approach based on a continuum description of the macromolecule and solvent. Although the continuum model of biomolecule electrostatics is widely used, the numerical techniques used to evaluate the model do not exploit fast solver approaches developed for analyzing integrated circuit interconnect. I will describe the formulation used for analyzing biomolecule electrostatics, and then derive an integral formulation of the problem that can be rapidly solved with precorrected-FFT method [4]. Singapore-MIT Alliance (SMA) 2003-11-20T22:30:58Z 2003-11-20T22:30:58Z 2003-01 Article http://hdl.handle.net/1721.1/3717 en_US High Performance Computation for Engineered Systems (HPCES); 15813 bytes application/pdf application/pdf
spellingShingle biomolecule electrostatics
continuum model
precorrected-FFT method
Kuo, Shihhsien
Altman, Michael D.
Bardhan, Jaydeep P.
Tidor, Bruce
White, Jacob K.
Fast Methods for Simulation of Biomolecule of Electrostatics
title Fast Methods for Simulation of Biomolecule of Electrostatics
title_full Fast Methods for Simulation of Biomolecule of Electrostatics
title_fullStr Fast Methods for Simulation of Biomolecule of Electrostatics
title_full_unstemmed Fast Methods for Simulation of Biomolecule of Electrostatics
title_short Fast Methods for Simulation of Biomolecule of Electrostatics
title_sort fast methods for simulation of biomolecule of electrostatics
topic biomolecule electrostatics
continuum model
precorrected-FFT method
url http://hdl.handle.net/1721.1/3717
work_keys_str_mv AT kuoshihhsien fastmethodsforsimulationofbiomoleculeofelectrostatics
AT altmanmichaeld fastmethodsforsimulationofbiomoleculeofelectrostatics
AT bardhanjaydeepp fastmethodsforsimulationofbiomoleculeofelectrostatics
AT tidorbruce fastmethodsforsimulationofbiomoleculeofelectrostatics
AT whitejacobk fastmethodsforsimulationofbiomoleculeofelectrostatics