Theoretical and Computational Chemistry

Optimized Halogen Atomic Radii for PBSA Calculations Using Off-Center Point-Charges

Paulo J. Costa BioISI – Biosystems & Integrative Sciences Institute, Faculty of Sciences, University of Lisboa


In force field methods, the usage of off-center point-charges, also called extra-points (EPs), is a common strategy to tackle the anisotropy of the electrostatic potential of covalently-bonded halogens (X), thus allowing the description of halogen bonds (XBs) at the molecular mechanics / molecular dynamics (MM/MD) level. Diverse EP implementations exist in the literature differing on the charge sets and/or the X–EP distances. Poisson–Boltzmann and surface area (PBSA) calculations can be used to obtain solvation free energies (∆G solv ) of small molecules, often to compute binding free energies (∆G bind ) at the MM PBSA level. This method depends, among other parameters, on the empirical assignment of atomic radii (PB radii). Given the multiplicity of off-center point-charges models and the lack of specific PB radii for halogens compatible with such implementations, in this work we assessed the performance of PBSA calculations for the estimation of ∆G solv values in water (∆G hyd ), also conducting an optimization of the halogen PB radii (Cl, Br, and I) for each EP model. We not only expand the usage of EP models in the scope of the General AMBER Force Field (GAFF) but also provide the first optimized halogen PB radii in the context of the CHARMM General Force Field (CGenFF), thus contributing to improving the description of halogenated compounds in PBSA calculations.

Version notes

Draft ver1


Thumbnail image of AFortuna_PJCosta_PBSA-manuscript.pdf
download asset AFortuna_PJCosta_PBSA-manuscript.pdf 0.85 MB [opens in a new tab]

Supplementary material

Thumbnail image of AFortuna_PJCosta_PBSA-SI.pdf
download asset AFortuna_PJCosta_PBSA-SI.pdf 0.26 MB [opens in a new tab]
AFortuna PJCosta PBSA-SI