Generalized Born (GB) models play an important role in biomolecular dynamics simulations. The issue of the accuracy of solvation forces computed according to these models has not been investigated in much detail to the best of our knowledge. Optimal atomic radii have been derived for most used force fields in order to reproduce molecular dynamics results under the Poisson–Boltzmann (PB) continuum model. In principle, if optimized radii are used, the PB model reproduces well MD simulations. For this reason in this contribution, we consider the most used
GB model implementations in publicly available molecular dynamics simulation softwares and assess the accuracy of GB solvation forces compared to the reference PB ones. The results show that current implementations that compute forces by pairwise summations correlate very well with PB forces if judicious choice of model parameters is performed. The latter can be obtained by linear fitting of solvation self-energies computed by the most accurate GBR6 GB model.
