A Comprehensive Review of Electrostatic Modelling in Biomolecular Simulations
DOI:
https://doi.org/10.71426/jasm.v1.i1.pp24-39Keywords:
Electrostatics, Biomolecules, Poisson-Boltzmann equation, MD SimulationAbstract
Electrostatic interactions are omnipresent in biology and play an important role in carrying out and determining the outcome of biological processes. Understanding biomolecular processes and the underlying phenomena, thus, requires an accurate representation of the system’s electrostatics. Computational methods like molecular dynamics simulations, owing to their atomistic resolution and the ability to describe biomolecular motion, interaction and structural fluctuations in atomic detail, have over time gained a lot of importance as a supporting tool for experiments in the quest to obtain mechanistic understanding of biological processes. With the important role that electrostatic interactions play in biology, accurate modelling of these interactions is important for obtaining reliable results from molecular simulations of biomolecules. In this review, we discuss the basics of biomolecular electrostatics, and state of the art and the recent developments in modelling electrostatic interaction in molecular dynamics simulations.
References
[1] Estévez-Torres A, Baigl D. DNA compaction: fundamentals and applications. Soft Matter. 2011;7(15):6746-6756. Available from: https://doi.org/10.1039/C1SM05373F
[2] Davey CA, Sargent DF, Luger K, Maeder AW, Richmond TJ. Solvent mediated interactions in the structure of the nucleosome core particle at 1.9 Å resolution. Journal of molecular biology. 2002 Jun 21;319(5):1097-1113. Available from: https://doi.org/10.1016/S0022-2836(02)00386-8
[3] Onufriev AV, Schiessel H. The nucleosome: from structure to function through physics. Current opinion in structural biology. 2019 Jun 1;56:119-130. Available from: https://doi.org/10.1016/j.sbi.2018.11.003
[4] Saurabh S, Glaser MA, Lansac Y, Maiti PK. Atomistic Simulation of Stacked Nucleosome Core Particles: Tail Bridging, the H4 Tail, and Effect of Hydrophobic Forces. The Journal of Physical Chemistry B 2016 Mar 31;120(12):3048-3060. Available from: https://pubmed.ncbi.nlm.nih.gov/26931280/
[5] Semenyuk P, Muronetz V. Protein interaction with charged macromolecules: From model polymers to unfolded proteins and post-translational modifications. International journal of molecular sciences. 2019 Mar 12;20(5):1252. Available from: https://doi.org/10.3390/ijms20051252
[6] Mukherjee A, Saurabh S, Olive E, Jang YH, Lansac Y. Protamine binding site on DNA: molecular dynamics simulations and free energy calculations with full atomistic details. The Journal of Physical Chemistry B. 2021 Mar 23;125(12):3032-3044. Available from: https://pubs.acs.org/doi/10.1021/acs.jpcb.0c09166?goto=supporting-info
[7] {Balhorn R. The protamine family of sperm nuclear proteins. Genome biology. 2007 Sep 26;8(9):227. Available from: https://doi.org/10.1186/gb-2007-8-9-227
[8] Bartlett GJ, Porter CT, Borkakoti N, Thornton JM. Analysis of catalytic residues in enzyme active sites. Journal of molecular biology. 2002 Nov 15;324(1):105-121 Available from: https://doi.org/10.1016/S0022-2836(02)01036-7
[9] Daubner SC, Fitzpatrick PF. Site-directed mutants of charged residues in the active site of tyrosine hydroxylase. Biochemistry. 1999 Apr 6;38(14):4448-4454. Available from: https://pubs.acs.org/doi/full/10.1021/bi983012u
[10] Russell AJ, Thomas PG, Fersht AR. Electrostatic effects on modification of charged groups in the active site cleft of subtilisin by protein engineering. Journal of molecular biology. 1987 Feb 20;193(4):803-813. Available from: https://doi.org/10.1016/0022-2836(87)90360-3
[11] Urry DW, Peng S, Parker T. Delineation of electrostatic-and hydrophobic-induced pka shifts in polypentapeptides: the glutamic acid residue. Journal of the American Chemical Society. 1993 Aug;115(16):7509-7510. Available from: https://pubs.acs.org/doi/pdf/10.1021/ja00069a063
[12] Davoodi J, Wakarchuk WW, Campbell RL, Carey PR, Surewicz WK. Abnormally High pKa of an Active‐Site Glutamic Acid Residue in Bacillus Circulans Xylanase: The Role of Electrostatic Interactions. European journal of biochemistry. 1995 Sep;232(3):839-843. Available from: https://doi.org/10.1111/j.1432-1033.1995.0839a.x
[13] Saurabh S, Zhang Q, Seddon JM, Lu JR, Kalonia C, Bresme F. Unraveling the microscopic mechanism of molecular ion interaction with monoclonal antibodies: impact on protein aggregation. Molecular Pharmaceutics. 2024 Feb 12;21(3):1285-1299. Available from: https://pubs.acs.org/doi/10.1021/acs.molpharmaceut.3c00963?goto=supporting-info
[14] Ren S. Current and emerging strategies for subcutaneous delivery of high-concentration and high-dose antibody therapeutics. Journal of Pharmaceutical Sciences. 2025 Aug 1;114(8):103877. Available from: https://doi.org/10.1016/j.xphs.2025.103877
[15] Saurabh S, Kalonia C, Li Z, Hollowell P, Waigh T, Li P, Webster J, Seddon JM, Lu JR, Bresme F. Understanding the stabilizing effect of histidine on mAb aggregation: a molecular dynamics study. Molecular Pharmaceutics. 2022 Aug 10;19(9):3288-3303. Available from: https://pubs.acs.org/doi/10.1021/acs.molpharmaceut.2c00453?goto=supporting-info
[16] Watson JD, Crick FH. The structure of DNA. In Cold Spring Harbor symposia on quantitative biology. 1953 Jan 1 (Vol. 18, pp. 123-131). Available from: https://symposium.cshlp.org/content/18/123.extract
[17] Alberts B, Johnson A, Lewis J, Raff M, Roberts K, Walter P. The structure and function of DNA. In Molecular Biology of the Cell. 4th edition 2002. Garland Science. Available from: https://www.ncbi.nlm.nih.gov/books/NBK26821/
[18] Schutz CN, Warshel A. What are the dielectric “constants” of proteins and how to validate electrostatic models?. Proteins: Structure, Function, and Bioinformatics 2001 Sep 1;44(4):400-417. Available from: https://doi.org/10.1002/prot.1106
[19] Saurabh S, Sivakumar PM, Perumal V, Khosravi A, Sugumaran A, Prabhawathi V. Molecular dynamics simulations in drug discovery and drug delivery. In Integrative nanomedicine for new therapies. 2020 Mar 3 (pp. 275-301). Cham: Springer International Publishing. Available from: https://doi.org/10.1007/978-3-030-36260-7_10
[20] Zhu X, Lopes PE, MacKerell Jr AD. Recent developments and applications of the CHARMM force fields. Wiley Interdisciplinary Reviews: Computational Molecular Science. 2012 Jan;2(1):167-185. Available from: https://doi.org/10.1002/wcms.74
[21] Kadaoluwa Pathirannahalage SP, Meftahi N, Elbourne A, Weiss AC, McConville CF, Padua A, Winkler DA, Costa Gomes M, Greaves TL, Le TC, Besford QA. Systematic comparison of the structural and dynamic properties of commonly used water models for molecular dynamics simulations. Journal of Chemical Information and Modeling. 2021 Aug 18;61(9):4521-4536. Available from: https://pubs.acs.org/doi/10.1021/acs.jcim.1c00794?goto=supporting-info
[22] Frenkel D, Smit B. Understanding molecular simulation: from algorithms to applications. Elsevier. 2023 Jul 13. Available from: https://doi.org/10.1016/B978-0-12-267351-1.X5000-7
[23] Luty BA, Davis ME, Tironi IG, Van Gunsteren WF. A comparison of particle-particle, particle-mesh and Ewald methods for calculating electrostatic interactions in periodic molecular systems. Molecular Simulation. 1994 Dec 1;14(1):11-20. Available from: https://doi.org/10.1080/08927029408022004
[24] Essmann U, Perera L, Berkowitz ML, Darden T, Lee H, Pedersen LG. A smooth particle mesh Ewald method. The Journal of chemical physics. 1995 Nov 15;103(19):8577-8593. Available from: https://doi.org/10.1063/1.470117
[25] Norberg J, Nilsson L. On the truncation of long-range electrostatic interactions in DNA. Biophysical journal. 2000 Sep 1;79(3):1537-1553. Available from: https://doi.org/10.1016/S0006-3495(00)76405-8
[26] Åqvist J, Hansson T. Analysis of electrostatic potential truncation schemes in simulations of polar solvents. The Journal of Physical Chemistry B. 1998 May 7;102(19):3837-3840. Available from: https://pubs.acs.org/doi/full/10.1021/jp973207w
[27] Loncharich RJ, Brooks BR. The effects of truncating long‐range forces on protein dynamics. Proteins: Structure, Function, and Bioinformatics. 1989;6(1):32-45. Available from: https://doi.org/10.1002/prot.340060104
[28] Tironi IG, Sperb R, Smith PE, van Gunsteren WF. A generalized reaction field method for molecular dynamics simulations. The Journal of chemical physics. 1995 Apr 1;102(13):5451-5459. Available from: https://doi.org/10.1063/1.469273
[29] Bergdorf M, Peter C, Hünenberger PH. Influence of cut-off truncation and artificial periodicity of electrostatic interactions in molecular simulations of solvated ions: A continuum electrostatics study. The Journal of chemical physics. 2003 Nov 1;119(17):9129-9144. Available from: https://doi.org/10.1063/1.1614202
[30] Resat H, McCammon JA. Correcting for electrostatic cutoffs in free energy simulations: Toward consistency between simulations with different cutoffs. The Journal of chemical physics. 1998 Jun 15;108(23):9617-9123. Available from: https://doi.org/10.1063/1.476437
[31] Fogolari F, Brigo A, Molinari H. The Poisson–Boltzmann equation for biomolecular electrostatics: a tool for structural biology. Journal of Molecular Recognition. 2002 Nov;15(6):377-392. Available from: https://doi.org/10.1002/jmr.577
[32] Lamm G. The poisson–boltzmann equation. Reviews in computational chemistry 2003 Oct 3;19:147-365. Available from: https://doi.org/10.1002/0471466638.ch4
[33] Matsarskaia O, Roosen‐Runge F, Schreiber F. Multivalent ions and biomolecules: Attempting a comprehensive perspective. ChemPhysChem. 2020 Aug 18;21(16):1742-1767. Available from: https://doi.org/10.1002/cphc.202000162
[34] Sparnaay MJ. On the electrostatic contribution to the interfacial tension of semiconductor/gas and semiconductor/electrolyte interfaces. Surface Science. 1964 Jul 1;1(3):213-224. Available from: https://doi.org/10.1016/0039-6028(64)90028-7
[35] Wu Y, Hu Q, Liang H, Wang A, Xu H, Wang L, He X. Electrostatic potential as solvent descriptor to enable rational electrolyte design for lithium batteries. Advanced Energy Materials. 2023 Jun;13(22):2300259. Available from: https://doi.org/10.1002/aenm.202300259
[36] Wan Y, Guo D, Hui X, Liu L, Yao Y. Studies on Hydration Swelling and Bound Water Type of Sodium‐and Polymer‐Modified Calcium Bentonite. Advances in Polymer Technology. 2020;2020(1):9361795. Available from: https://doi.org/10.1155/2020/9361795
[37] Brown MA, Goel A, Abbas Z. Effect of electrolyte concentration on the stern layer thickness at a charged interface. Angewandte Chemie. 2016 Mar 7;128(11):3854-3858. Available from: https://doi.org/10.1002/ange.201512025
[38] Mangelsdorf CS, White LR. The dynamic double layer Part 1Theory of a mobile Stern layer. Journal of the Chemical Society, Faraday Transactions. 1998;94(16):2441-2452. Available from: https://doi.org/10.1039/A803588A
[39] Hartley GS. The application of the Debye-Hückel theory to colloidal electrolytes. Transactions of the Faraday Society. 1935;31:31-50. Available from: https://pubs.rsc.org/en/content/articlelanding/1935/tf/tf9353100031/unauth
[40] Fogolari F, Brigo A, Molinari H. The Poisson–Boltzmann equation for biomolecular electrostatics: a tool for structural biology. Journal of Molecular Recognition. 2002 Nov;15(6):377-392. Available from: https://doi.org/10.1002/jmr.577
[41] Baker NA. Poisson–Boltzmann methods for biomolecular electrostatics. In Methods in enzymology ,Academic Press. 2004 Jan 1 (Vol. 383, pp. 94-118). Available from: https://doi.org/10.1016/S0076-6879(04)83005-2
[42] de Souza Gama M, Santos MS, de Almeida Lima ER, Tavares FW, Barreto Jr AG. A modified Poisson-Boltzmann equation applied to protein adsorption. Journal of Chromatography A. 2018 Jan 5;1531:74-82. Available from: https://doi.org/10.1016/j.chroma.2017.11.022
[43] Li B. Minimization of electrostatic free energy and the Poisson–Boltzmann equation for molecular solvation with implicit solvent. SIAM Journal on Mathematical Analysis. 2009;40(6):2536-2566. Available from: https://doi.org/10.1137/080712350
[44] Hu J, Zhao S, Geng W. Accurate pKa computation using matched interface and boundary (MIB) method based Poisson-Boltzmann solver. Commun. Comput. Phys. 2018 Feb 1;23(2):520-539. Available from: https://s2.smu.edu/wgeng/research/Hu_CiCP_2018.pdf
[45] Griffiths DJ. Introduction to electrodynamics. Cambridge University Press. 2023 Nov 2. Available from: https://www.cambridge.org/highereducation/books/introduction-to-electrodynamics/3AB220820DBB628E5A43D52C4B011ED4
[46] Sharp K. Incorporating solvent and ion screening into molecular dynamics using the finite‐difference Poisson–Boltzmann method. Journal of computational chemistry. 1991 May;12(4):454-468. Available from: https://doi.org/10.1002/jcc.540120407
[47] Baker NA, Sept D, Joseph S, Holst MJ, McCammon JA. Electrostatics of nanosystems: application to microtubules and the ribosome. Proceedings of the National Academy of Sciences. 2001 Aug 28;98(18):10037-10041. Available from: https://doi.org/10.1073/pnas.181342398
[48] Jurrus E, Engel D, Star K, Monson K, Brandi J, Felberg LE, Brookes DH, Wilson L, Chen J, Liles K, Chun M. Improvements to the APBS biomolecular solvation software suite. Protein science. 2018 Jan;27(1):112-128. Available from: https://doi.org/10.1002/pro.3280
[49] Dhatt G, Lefrançois E, Touzot G. Finite element method. John Wiley $&$ Sons. 2012 Dec 27. Available from: https://onlinelibrary.wiley.com/doi/book/10.1002/9781118569764
[50] Boschitsch AH, Fenley MO, Zhou HX. Fast boundary element method for the linear Poisson− Boltzmann equation. The Journal of Physical Chemistry B. 2002 Mar 14;106(10):2741-2754. Available from: https://pubs.acs.org/doi/full/10.1021/jp013607q
[51] Li Z, Saurabh S, Hollowell P, Kalonia CK, Waigh TA, Li P, Webster JR, Seddon JM, Bresme F, Lu JR. pH-Dependent Conformational Plasticity of Monoclonal Antibodies at the SiO2/Water Interface: Insights from Neutron Reflectivity and Molecular Dynamics. ACS Applied Materials $&$ Interfaces. 2024 Dec 12;16(51):70231-70241. Available from: https://pubs.acs.org/doi/10.1021/acsami.4c14407
[52] Onufriev AV, Case DA. Generalized born implicit solvent models for biomolecules. Annual review of biophysics. 2019 May 6;48(1):275-296. Available from: https://doi.org/10.1146/annurev-biophys-052118-115325
[53] Bashford D, Case DA. Generalized born models of macromolecular solvation effects. Annual review of physical chemistry. 2000 Oct;51(1):129-152. Available from: https://doi.org/10.1146/annurev.physchem.51.1.129
[54] Bauer MR, Mackey MD. Electrostatic complementarity as a fast and effective tool to optimize binding and selectivity of protein–ligand complexes. Journal of medicinal chemistry. 2019 Feb 26;62(6):3036-3050. Available from: https://pubs.acs.org/doi/full/10.1021/acs.jmedchem.8b01925
[55] Tilocca A. Structural models of bioactive glasses from molecular dynamics simulations. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences. 2009 Apr 1;465(2104):1003-27. Available from: https://doi.org/10.1098/rspa.2008.0462
[56] Gruian C, Vanea E, Steinhoff HJ, Simon S. Glass-ceramics: fundamental aspects regarding the interaction with proteins. Handbook of Bioceramics and Biocomposites. 2016:667-701. Available from: https://doi.org/10.1007/978-3-319-09230-0_20-1
[57] Rosengren Å, Oscarsson S, Mazzocchi M, Krajewski A, Ravaglioli A. Protein adsorption onto two bioactive glass-ceramics. Biomaterials. 2003 Jan 1;24(1):147-155 Available from: https://doi.org/10.1016/S0142-9612(02)00272-7
[58] Verné E, Bretcanu O, Balagna C, Bianchi CL, Cannas M, Gatti S, Vitale-Brovarone C. Early stage reactivity and in vitro behavior of silica-based bioactive glasses and glass-ceramics. Journal of Materials Science: Materials in Medicine. 2009 Jan;20(1):75-87. Available from: https://doi.org/10.1007/s10856-008-3537-8
[59] Szefczyk B, Mulholland AJ, Ranaghan KE, Sokalski WA. Differential transition-state stabilization in enzyme catalysis: quantum chemical analysis of interactions in the chorismate mutase reaction and prediction of the optimal catalytic field. Journal of the American Chemical Society. 2004 Dec 15;126(49):16148-16159. Available from: https://pubs.acs.org/doi/full/10.1021/ja049376t
[60] Søndergaard CR, Olsson MH, Rostkowski M, Jensen JH. Improved treatment of ligands and coupling effects in empirical calculation and rationalization of p K a values. Journal of chemical theory and computation. 2011 Jul 12;7(7):2284-3395. Available from: https://pubs.acs.org/doi/full/10.1021/ct200133y
[61] Olsson MH. Protein electrostatics and pKa blind predictions; contribution from empirical predictions of internal ionizable residues Proteins: Structure, Function, and Bioinformatics. 2011 Dec;79(12):3333-3345. Available from: https://doi.org/10.1002/prot.23113Digital
[62] Pierce BG, Wiehe K, Hwang H, Kim BH, Vreven T, Weng Z. ZDOCK server: interactive docking prediction of protein–protein complexes and symmetric multimers Bioinformatics. 2014 Jun 15;30(12):1771-1773. Available from: https://doi.org/10.1093/bioinformatics/btu097
[63] Cosconati S, Forli S, Perryman AL, Harris R, Goodsell DS, Olson AJ. Virtual screening with AutoDock: theory and practice. Expert opinion on drug discovery. 2010 Jun 1;5(6):597-607. Available from: https://doi.org/10.1517/17460441.2010.484460
[64] Hou T, Wang J, Li Y, Wang W. Assessing the performance of the molecular mechanics/Poisson Boltzmann surface area and molecular mechanics/generalized Born surface area methods. II. The accuracy of ranking poses generated from docking. Journal of computational chemistry. 2011 Apr 15;32(5):866-877. Available from: https://doi.org/10.1002/jcc.21666
[65] Jing Z, Liu C, Cheng SY, Qi R, Walker BD, Piquemal JP, Ren P. Polarizable force fields for biomolecular simulations: Recent advances and applications. Annual Review of biophysics. 2019 May 6;48(1):371-394. Available from: https://doi.org/10.1146/annurev-biophys-070317-033349
[66] Mauger N, Plé T, Lagardère L, Huppert S, Piquemal JP. The Q-AMOEBA (CF) Polarizable Potential. The Journal of Physical Chemistry Letters. 2025 Jun 2;16(23):5723-31. Available from: https://pubs.acs.org/doi/full/10.1021/acs.jpclett.5c00797
[67] Lemkul JA, Huang J, Roux B, MacKerell Jr AD. An empirical polarizable force field based on the classical drude oscillator model: development history and recent applications. Chemical reviews. 2016 May 11;116(9):4983-5013. Available from: https://pubs.acs.org/doi/full/10.1021/acs.chemrev.5b00505
[68] Vanommeslaeghe K, MacKerell Jr AD. CHARMM additive and polarizable force fields for biophysics and computer-aided drug design. Biochimica et Biophysica Acta (BBA)-General Subjects. 2015 May 1;1850(5):861-871. Available from: https://doi.org/10.1016/j.bbagen.2014.08.004
[69] de Oliveira VM, Liu R, Shen J. Constant pH molecular dynamics simulations: Current status and recent applications. Current opinion in structural biology. 2022 Dec 1;77:102498. Available from: https://doi.org/10.1016/j.sbi.2022.102498
[70] Henderson JA, Liu R, Harris JA, Huang Y, de Oliveira VM, Shen J. A guide to the continuous constant pH molecular dynamics methods in Amber and CHARMM [Article v1. 0]. Living journal of computational molecular science. 2022 Aug 22;4(1):1563. Available from: https://doi.org/10.33011/livecoms.4.1.1563
[71] Aho N, Buslaev P, Jansen A, Bauer P, Groenhof G, Hess B. Scalable constant pH molecular dynamics in GROMACS. Journal of Chemical Theory and Computation. 2022 Sep 21;18(10):6148-6160. Available from: https://pubs.acs.org/doi/full/10.1021/acs.jctc.2c00516
[72] Demichelis R, Noel Y, D’Arco P, Rérat M, Zicovich-Wilson CM, Dovesi R. Properties of carbon nanotubes: an ab initio study using large gaussian basis sets and various dft functionals The Journal of Physical Chemistry C. 2011 May 12;115(18):8876-85. Available from: https://pubs.acs.org/doi/10.1021/jp110704x
[73] Neese F. The ORCA program system. Wiley Interdisciplinary Reviews: Computational Molecular Science. 2012 Jan;2(1):73-8. Available from: https://doi.org/10.1002/wcms.81
[74] Hutter J, Iannuzzi M, Schiffmann F, VandeVondele J. cp2k: atomistic simulations of condensed matter systems. Wiley Interdisciplinary Reviews: Computational Molecular Science. 2014 Jan;4(1):15-25. Available from: https://doi.org/10.1002/wcms.1159
[75] Götz AW, Clark MA, Walker RC. An extensible interface for QM/MM molecular dynamics simulations with AMBER. Journal of computational chemistry. 2014 Jan 15;35(2):95-108. Available from: https://doi.org/10.1002/jcc.23444
[76] Melo MC, Bernardi RC, Rudack T, Scheurer M, Riplinger C, Phillips JC, Maia JD, Rocha GB, Ribeiro JV, Stone JE, Neese F. NAMD goes quantum: an integrative suite for hybrid simulations. Nature methods. 2018 May 1;15(5):351-354. Available from: https://doi.org/10.1038/nmeth.4638
[77] Zeng J, Zhang D, Lu D, Mo P, Li Z, Chen Y, Rynik M, Huang LA, Li Z, Shi S, Wang Y. DeePMD-kit v2: A software package for deep potential models. The Journal of Chemical Physics. 2023 Aug 7;159(5). Available from: https://doi.org/10.1063/5.0155600
[78] Xue J, Terrel N, Pickering I, Roitberg A. LAMMPS-ANI: Large Scale Molecular Dynamics Simulations with ANI Neural Network Potential. chemrxiv.org. Available from: https://doi.org/10.26434/chemrxiv-2025-8v03m
[79] Schütt KT, Sauceda HE, Kindermans PJ, Tkatchenko A, Müller KR. Schnet–a deep learning architecture for molecules and materials. The Journal of chemical physics. 2018 Jun 28;148(24). Available from: https://doi.org/10.1063/1.5019779
Downloads
Published
Issue
Section
License
Copyright (c) 2025 Adari Sasi Sahasra, Kowsar Abdul, Mukherjee Arnab, Suman Saurabh (Author)
This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.
The Journal of Applied Sciences and Modelling publishes all articles under the Creative Commons Attribution–NonCommercial 4.0 International License (CC BY-NC 4.0). This license permits others to copy, distribute, reproduce, remix, adapt, and build upon the published work for non-commercial purposes, provided appropriate credit is given to the original authors and the source. By publishing in the Journal of Applied Sciences and Modelling, all authors agree to these licensing terms as a condition of publication.
