Halogen bond is an important non-covalent interaction which is receiving a growing attention in the study of protein-ligand complexes. Many drugs are halogenated molecules and it has been recently shown that many halogenated ligands establish halogen bonds with biomolecules. As the halogen bond nature is due to an anisotropy of the electrostatic potential around halogen atoms, it is not possible to use traditional force fields based on a set of atom-centred charges to study halogen bonds in biomolecules. We show that the introduction of pseudo-atoms on halogens permits us to correctly describe the anisotropy of the electrostatic potential and to perform molecular dynamics simulations on complexes of proteins with halogenated ligands that reproduce experimental values. The results are compared with crystallographic data and with hybrid quantum mechanics/molecular mechanics calculations.

Halogen bonding is a noncovalent interaction between a halogen atom and a nucleophilic site. Interactions involving the π electrons of aromatic rings have received, up to now, little attention, despite the large number of systems in which they are present. We report binding energies of the interaction between either NCX or PhX (X = F, Cl, Br, I) and the aromatic benzene system as determined with the coupled cluster with perturbative triple excitations method [CCSD(T)] extrapolated at the complete basis set limit. Results are compared with those obtained by Møller-Plesset perturbation theory to second order (MP2) and density functional theory (DFT) calculations by using some of the most common functionals. Results show the important role of DFT in studying this interaction.

The performance of an extensive set of density functional theory functionals has been tested against CCSD(T) and MP2 results, extrapolated to the complete basis set (CBS) limit, for the interaction of either DCl or DBr (D = H, HCC, F, and NC) with the aromatic system of benzene. It was found that double hybrid functionals explicitly including dispersion, that is, B2PLYPD and mPW2PLYPD, provide the better agreement with the CCSD(T)/CBS results on both energies and equilibrium geometry, indicating the importance of dispersive contributions in determining this interaction. Among the less expensive functionals, the better performance is provided by the ωB97X and M062X functionals, while the ωB97XD and B97D functionals are shown to work very well for bromine complexes but not so well for chlorine complexes.

The novel strigolactam derivatives described here compare favourably with the corresponding GR-28 strigolactones in terms of biological activity and physicochemical properties. However, we believe strigolactone and strigolactam derivatives require further structural optimisation to improve their soil persistence to demonstrate a potential for agronomical applications. © 2016 Society of Chemical Industry.