In the past several years, halogen bonds have been shown to be relevant in crystal engineering and biomedical applications. One of the reasons for the utility of these types of noncovalent interactions in the development of, for example, pharmaceutical ligands is that their strengths and geometric properties are very tunable. That is, substitution of atoms or chemical groups in the vicinity of a halogen can have a very strong effect on the strength of the halogen bond. In this study we investigate halogen-bonding interactions involving aromatically-bound halogens (Cl, Br, and I) and a carbonyl oxygen. The properties of these halogen bonds are modulated by substitution of aromatic hydrogens with fluorines, which are very electronegative. It is found that these types of substitutions have dramatic effects on the strengths of the halogen bonds, leading to interactions that can be up to 100% stronger. Very good correlations are obtained between the interaction energies and the magnitudes of the positive electrostatic potentials (σ-holes) on the halogens. Interestingly, it is seen that the substitution of fluorines in systems containing smaller halogens results in electrostatic potentials resembling those of systems with larger halogens, with correspondingly stronger interaction energies. It is also shown that aromatic fluorine substitutions affect the optimal geometries of the halogen-bonded complexes, often as the result of secondary interactions.

Noncovalent interactions of the polyhedral carborane 1-carba-closo-dodecaborane (CB(11)H(12))(-) with building blocks of biomolecules, modelled by glycine (GLY), serine (SER), phenylalanine (PHE), glutamic acid (GLU), lysine (LYS) and arginine (ARG), were investigated in vacuo by molecular dynamics simulations with the UFF empirical potential. Selected structures were further studied by accurate ab initio quantum chemical procedures. Interactions with a peptide bond (GLY-SER dipeptide) and a nucleic acid building block (guanine) were also considered. The RESP and NPA charges of carboranes and small model systems are compared and their use is discussed. The dominant interaction between carboranes and biomolecules is the formation of unconventional proton-hydride hydrogen bonds (dihydrogen bonds) characterized by a short distance between hydrogen atoms (as close as 1.8 A) and an average strength in the range of 4.2-5.8 kcal mol(-1). The total stabilization energy of complexes investigated is rather large, and the largest value (approximately 15 kcal mol(-1)) was found for the carborane complexes with ARG and the GLY-SER dipeptide. These interactions are ubiquitous under geometrical constraints influencing the strength of the interaction. The carborane forms dihydrogen bonds with biomolecules preferably with the hydrogen atoms of its lower hemisphere (i.e. the part of the cage opposite to the carbon atom). These two geometrical factors can be used to explain the specificity of inhibition of HIV protease by carboranes.

In this study, we introduce a fast and reliable rescoring scheme for docked complexes based on a semiempirical quantum mechanical PM6-DH2 method. The method utilizes a PM6-based Hamiltonian with corrections for dispersion energy and hydrogen bonds. The total score is constructed as the sum of the PM6-DH2 interaction enthalpy, the empirical force field (AMBER) interaction entropy, and the sum of the deformation (PM6-DH2, SMD) and the desolvation (SMD) energies of the ligand. The main advantage of the procedure is the fact that we do not add any empirical parameter for either an individual component of the total score or an individual protein-ligand complex. This rescoring method is applied to a very challenging system, namely, the HIV-1 protease with a set of ligands. As opposed to the conventional DOCK procedure, the PM6-DH2 rescoring based on all of the terms distinguishes between binders and nonbinders and provides a reliable correlation of the theoretical and experimental binding free energies. Such a dramatic improvement, resulting from the PM6-DH2 rescoring of all the complexes, provides a valuable yet inexpensive tool for rational drug discovery and de novo ligand design.