An n log n Generalized Born Approximation.

Area of Science:

• Computational Chemistry
• Biophysics
• Molecular Dynamics

Background:

• Generalized Born (GB) implicit solvation models accelerate molecular dynamics (MD) simulations by approximating solvent effects.
• The computational cost of GB models scales quadratically with system size (∼n²), limiting their application to large biomolecules.
• Existing approximations, like spherical cutoffs, can oversimplify interactions and reduce accuracy.

Purpose of the Study:

• To develop and validate an efficient and accurate approximation for the Generalized Born (GB) implicit solvent model.
• To reduce the computational scaling of GB models from quadratic (∼n²) to near-linear (∼n log n).
• To improve the accuracy of implicit solvent models in capturing biomolecular dynamics.

Main Methods:

• Implementation of the Hierarchical Charge Partitioning (HCP) method for GB (HCP-GB) in MD software.
• HCP approximates charge distributions of molecular components at different hierarchical levels.
• HCP-GB defines component-based effective Born radii for approximating distant interactions.

Main Results:

• The HCP-GB method achieves a computational scaling of ∼n log n, significantly outperforming the standard GB model.
• HCP-GB demonstrated 1.1–390 times speedup over reference GB computations across various biomolecular sizes.
• HCP-GB showed higher accuracy than cutoff-GB, with better reproduction of electrostatic forces and molecular dynamics characteristics.

Conclusions:

• The HCP-GB method offers a computationally efficient and accurate alternative to standard GB models for large-scale MD simulations.
• HCP-GB provides a more physically realistic treatment of distant interactions compared to cutoff-GB, avoiding potential inaccuracies.
• This method is recommended for MD simulations employing pairwise implicit solvent GB models, enhancing the study of complex biomolecular systems.