Abstract:Halogenation is a commonly used strategy in drug design, not only due to the improvement of ADME properties of the compounds but also due to the ability of halogens to engage in halogen bonds (XBs) with biological targets such as proteins. The prediction of binding free energies (\dGbind{}) is paramount in drug development, and, in this context, MM-PBSA calculations are a fast and commonly used method to achieve such a task. However, the usage of parameters to describe halogen anisotropy, such as extra-points of charge (EPs) that allow for the sampling of XBs, in MM-PBSA calculations is not common. Despite the recent development of optimized radii for halogens (r$_{opt}$), compatible with EP implementations, providing accurate hydration free energies at the PBSA level, the impact of these parameters in the calculation of binding free energies (\dGbind{}) on protein-ligand complexes was yet to be tackled. In this work, we evaluated the performance of r$_{opt}$ in the estimation of \dGbind{} values using MM-PBSA calculations for three sets of CK2 inhibitors. for which the experimental binding energies are available. We tested the performance of three EP models and standard RESP charges (i.e. when no EP is added), along with the effect of the internal dielectric constant ($\epsilon$_{\textnormal{in}}$) and sampling time. Our results indicate that, while directly using the X-ray structure does not typically lead to good results, relaxation through MM minimization, especially with $\epsilon$_{\textnormal{in}}$ = 4, leads to very decent correlations with the experimental values. Better correlations were achieved when adding configurational sampling via MD simulations, though, the Pearson coefficients varied substantially with sampling time and PBSA setup used, indicating that longer sampling times do not necessarily lead to higher correlations with experimental values. For all cases, optimized halogen radii (r$_{opt}$) leads to slight improvements whereas the description of the halogen anisotropy \textit{via} and EP considerably improved the correlations for \textbf{Set A} and \textbf{Set C}. In addition, our findings suggest that $\epsilon$_{\textnormal{in}}$ has a considerable impact on the prediction of the \dGbind{} values, with the usage of an EP being enough to account for polarizability for systems with several halogens, whereas for systems with fewer halogens, a higher $\epsilon$_{\textnormal{in}}$ might be needed. Moreover, the use of an EP allowed the sampling of XBs without impairing the sampling of hydrogen bonds (HBs). To date, this study is the first to provide a comprehensive benchmark using MM-PBSA calculations, thus giving insights into which parameters contribute to a better description of halogen bonding systems.

Impact of Halogen Bonds on Protein-Peptide Binding and Protein Structural Stability Revealed by Computational Approaches

Geometric Characteristics and Energy Landscapes of Halogen–water–hydrogen Bridges at Protein–ligand Interfaces

Preferred microenvironments of halogen bonds and hydrogen bonds revealed using statistics and QM/MM calculation studies

Rediscovery of halogen bonds in protein-ligand complexes.

Halogen-water-hydrogen Bridges in Biomolecules.

Halogen Bonds Exist between Noncovalent Ligands and Natural Nucleic Acids

Halogen-Ionic Bridges: Do They Exist in the Biomolecular World?

Halogen Bonding in Nucleic Acid Complexes

Halogen Bonding on Water─A Drop in the Ocean?

Ab initio investigation of the complexes between bromobenzene and several electron donors: some insights into the magnitude and nature of halogen bonding interactions.

Data and Molecular Fingerprint-Driven Machine Learning Approaches to Halogen Bonding

Halogen Bonding in Haspin-Halogenated Tubercidin Complexes: Molecular Dynamics and Quantum Chemical Calculations

Halogen Bonding: an AIM Analysis of the Weak Interactions

Mechanism and Application of Halogen Bond Induced Fluorescence Enhancement and Iodine Molecule Cleavage in Solution

Impact of the halogen PB radii in the estimation of protein-ligand binding energies using MM-PBSA calculations

Modulation of halogen-bonded 2D self-assemblies of benzothiadiazole derivative: Concentration and solvent effects

Assessment of Halogen Off-Center Point-Charge Models Using Explicit Solvent Simulations

Ab Initio Study of the Complexes of Halogen-Containing Molecules RX (x=cl, Br, and I) and NH3: Towards Understanding the Nature of Halogen Bonding and the Electron-Accepting Propensities of Covalently Bonded Halogen Atoms.

Volumetric Assessment of Glioma Removal by Intraoperative High-field Magnetic Resonance Imaging.

Competition of Hydrogen Bonds, Halogen Bonds, and π–π Interactions in Crystal Structures. Exploring the Effect of One Atom Substitution

Halogen Bond Structure and Dynamics from Molecular Simulations