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.

In this work we have investigated hydration properties of aqueous solutions up to a solute molar fraction X₂ =0.125 of two isosteric molecules – the bioprotectant trimethylamine-N-oxide (TMAO) and the denaturant tert-butyl alcohol (TBA) – using molecular dynamics simulation at 298 K. Statistical analyses of the trajectories show in particular that as the solute concentration increases the number of the water molecules in the first hydration shell decreases uniformly for TMAO, while for TBA it decreases more rapidly in a concentration range where experiments indicate that TBA starts to self-aggregate. No appreciable solute segregation occurs for TMAO even in the most concentrated solution, where on the average each water molecule is shared by two solutes. This result parallels what has been recently found for glycine betaine, an organic osmolyte closely related to TMAO.