Related Experiment Video
Updated: Jul 4, 2026
12:11
Computation of Atmospheric Concentrations of Molecular Clusters from ab initio Thermochemistry
Published on: April 8, 2020
Accurate solution of multi-region continuum biomolecule electrostatic problems using the linearized Poisson-Boltzmann
Michael D Altman1, Jaydeep P Bardhan, Jacob K White
1Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, USA.
Journal of Computational Chemistry
|June 21, 2008
Summary
This study introduces a boundary-element method (BEM) for biomolecular electrostatics, improving accuracy and efficiency for complex molecular geometries. The novel approach enhances calculations of electrostatic solvation and binding energies.
Area of Science:
- Computational Biology
- Biophysics
- Biomolecular Electrostatics
Background:
- Continuum models are crucial for understanding biological molecule electrostatics.
- Accurate representation of molecular geometries and topologies is essential.
- Existing methods face challenges with complex biomolecular structures.
Purpose of the Study:
- To develop a precise boundary-element method (BEM) solver for biomolecular electrostatics.
- To accurately model complex molecular geometries and topologies.
- To improve the calculation of electrostatic energies in biological systems.
Main Methods:
- Implemented curved boundary elements for faithful geometric representation.
- Utilized preconditioned iterative methods (GMRES) with matrix compression (FFTSVD) to solve linear systems efficiently.
- Employed robust numerical integration for singular and near-singular integrals.
- Developed a boundary-integral approach for multiple dielectric regions, salt effects, and point charges.
Main Results:
- The BEM implementation accurately solves linearized Poisson-Boltzmann equation problems.
- Demonstrated improved convergence and significant impact on computed energetics for solvation and binding free energies.
- Showcased the importance of curved-element BEM accuracy for nonrigid-binding models.
- Achieved high accuracy in comparable compute times to finite-difference methods for multiple solves.
Conclusions:
- The presented BEM solver offers enhanced accuracy and efficiency for biomolecular electrostatics.
- This method is particularly beneficial for calculations involving complex geometries and advanced modeling techniques.
- The BEM approach provides a robust tool for charge optimization and component analysis in biomolecular systems.
Related Concept Videos
01:27Electrostatic Boundary Conditions in Dielectrics
When an electric field passes from one homogeneous medium to another, crossing the boundary between the two mediums imparts a discontinuity in the electric field. This results in electrostatic boundary conditions that depend on the type of mediums the field propagates through.
Consider a case where both the mediums across a boundary are two different dielectric materials. Recall that the electric field and electric displacement are proportional and related through the material's permittivity.
Consider a case where both the mediums across a boundary are two different dielectric materials. Recall that the electric field and electric displacement are proportional and related through the material's permittivity.
01:16Electrostatic Boundary Conditions
Consider an external electric field propagating through a homogeneous medium. When the electric field crosses the surface boundary of the medium, it undergoes a discontinuity. The electric field can be resolved into normal and tangential components. The amount by which the field changes at any boundary is given by the difference between the field components above and below the surface boundary.
The surface integral of an electric field is given by Gauss's law in integral form and is related to...
The surface integral of an electric field is given by Gauss's law in integral form and is related to...
01:25Poisson's And Laplace's Equation
The electric potential of the system can be calculated by relating it to the electric charge densities that give rise to the electric potential. The differential form of Gauss's law expresses the electric field's divergence in terms of the electric charge density.
01:20Maxwell-Boltzmann Distribution: Problem Solving
Individual molecules in a gas move in random directions, but a gas containing numerous molecules has a predictable distribution of molecular speeds, which is known as the Maxwell-Boltzmann distribution, f(v).
This distribution function f(v) is defined by saying that the expected number N (v1,v2) of particles with speeds between v1 and v2 is given by
This distribution function f(v) is defined by saying that the expected number N (v1,v2) of particles with speeds between v1 and v2 is given by
01:23Equilibrium Conditions for a Particle
When an object is in equilibrium, it is either at rest or moving with a constant velocity. There are two types of equilibrium: static and dynamic. Static equilibrium occurs when an object is at rest, while dynamic equilibrium occurs when an object is moving with a constant velocity. In both cases, there must be a balance of forces acting on the object.
To understand the concept of equilibrium, let us first consider the forces acting on an object. When different forces act on an object, they can...
To understand the concept of equilibrium, let us first consider the forces acting on an object. When different forces act on an object, they can...
01:20Region of Convergence of Laplace Tarnsform
The Region of Convergence (ROC) is a fundamental concept in signal processing and system analysis, particularly associated with the Laplace transform. The ROC represents an area in the complex plane where the Laplace transform of a given signal converges, determining the transform's applicability and utility.
Consider a decaying exponential signal that begins at a specific time. When deriving its Laplace transform, the time-domain variable is replaced with a complex variable. This substitution...
Consider a decaying exponential signal that begins at a specific time. When deriving its Laplace transform, the time-domain variable is replaced with a complex variable. This substitution...