Investigations into the Nature of Halogen Bonding Including Symmetry Adapted Perturbation Theory Analyses

In recent years it has been recognized that, because of their unique properties, halogen bonds have tremendous potential in the development of new pharmaceutical compounds and materials. In this study we investigate the phenomenon of halogen bonding by carrying out ab initio calculations on the halomethane-formaldehyde complexes as well as the fluorine substituted FnH3-nCX···OCH2 dimers, where the halogen bonding halogens (X) are chlorine, bromine, and iodine. Coupled cluster (CCSD(T)/aug-cc-pVTZ) calculations indicate that the binding energies for these type of interactions lie in the range between -1.05 kcal/mol (H3CCl···OCH2) and -3.72 kcal/mol (F3CI···OCH2). One of the most important findings in this study is that, according to symmetry adapted perturbation theory (SAPT) analyses, halogen bonds are largely dependent on both electrostatic and dispersion type interactions. As the halogen atom involved in halogen bonding becomes larger the interaction strength for this type of interaction also gets larger and, interestingly, more electrostatic (and less dispersive) in character. Halogen bonding interactions also become stronger and more electrostatic upon substitution of (the very electronegative) fluorines onto the halomethane molecule.

Energy Matrix of Structurally Important Side-Chain/Side-Chain Interactions in Proteins

The interactions between amino acid side chains in proteins are generally considered to be the most important stabilizing factor controlling the precise arrangement of the polypeptide chain into a well-defined spatial structure. We used the RI-DFT-D method to calculate the full 20 × 20 matrix of interaction energies between all pairs of amino acid side chains. For each pair, we used a representative 3D conformation extracted from an analysis of known protein structures from Protein Data Bank (PDB). The representative comes from the largest cluster of relative orientations of the two side chains. We find that all of the calculated interaction energies between selected pairs of amino acids are attractive in the gas phase with the exception of side chain pairs having the same total charge. We compared these data with those calculated by the parm03 and OPLS-AA/L force fields to investigate the reliability of simple methods in modeling biomolecules and their behavior. The force fields yield good overall interaction energies for our set but have problems in evaluation of some particular interactions which could be of principal importance for protein stability. We then looked in detail at the 20 side chain interactions involving tryptophan. The histograms of interaction energies showed that the distributions of the interaction energies are neither normal nor Boltzmann-like and that our representative geometries correspond mostly to the minimum energy geometry which is rather poorly populated in the whole pairwise energy distribution. We concluded that cluster representatives obtained by the clusterization algorithm based on geometry criteria cannot be considered as a typical interaction for the whole side chain/side chain interaction distribution. They seem to epitomize the strongest interactions in a protein and are often functionally or structurally important.

Representative Amino Acid Side Chain Interactions in Proteins. A Comparison of Highly Accurate Correlated ab Initio Quantum Chemical and Empirical Potential Procedures

Interactions between amino acid side chains play a crucial role both within a folded protein and between the interacting protein molecules. Here we have selected a representative set of 24 of the 400 (20 × 20) possible interacting side chain pairs based on data from Atlas of Protein Side-Chain Interactions. For each pair, we obtained its most favorable interaction geometry from the structural data and computed the interaction energy in the gas phase using several different, commonly used, ab initio and force field methods, namely Møller-Plesset perturbation theory (MP2), density functional theory combined with symmetry-adapted perturbation theory (DFT-SAPT), density functional theory empirically augmented with an empirical dispersion term (DFT-D), and empirical potentials using the OPLS-AA/L and Amber03 force fields. All the methods were compared against a reference method taken to be the CCSD(T) level of theory extrapolated to the complete basis set limit. We found a high degree of agreement between the different methods, even though the range of binding energies obtained was extremely large. The most computationally intensive methods yielded the best results. Among the less computationally time-consuming methods, the DFT-D method as well as parm03 force field provided consistently good results when compared to the reference values. We also tested how representative the chosen geometries of the side chains were and investigated the effect on the binding energies of the dielectric constant of the surrounding medium.

Extensions and applications of the A24 data set of accurate interaction energies

The A24 data set (Řezáč and Hobza, J. Chem. Theory Comput. 2013, 9, 2151-2155) is a set of noncovalent complexes large enough to showcase various types of interactions yet small enough to make highly accurate calculations possible. It is intended for the testing of accurate computational methods which are then used as a benchmark in larger model systems. In this work, we improve the best estimate of the interaction energies in the set by updating their CCSD(T)/CBS and CCSDT(Q) components with calculations in larger basis sets. The data set is then used to test a large number of composite CCSD(T) approaches. Special attention is paid to the use of the explicitly correlated MP2-F12 method in these composite calculations. It is shown that an accuracy of 1-2% can be achieved with setups applicable to larger molecules. The effect of frozen natural orbital approximation on the accuracy of composite CCSD(T)/CBS calculations is also quantified. In four trivially saturated complexes where CCSDT(Q)/CBS data are now available, the convergence of the many-body correlation effects is assessed by fixed-node diffusion Monte Carlo (FN-DMC) calculations. A good agreement is achieved between FN-DMC and high-level coupled-cluster which represents an important cross-check of both approaches.

The properties of substituted 3D-aromatic neutral carboranes: the potential for σ-hole bonding

The calculated properties of substituted carboranes such as dipole moment, polarisability, the magnitude of the σ-hole and the desolvation free energy are compared with these properties in comparable aromatic and cyclic aliphatic organic compounds. Dispersion and charge transfer energies are similar. However, the predicted strength of the halogen bonds with the same electron donor (based on the magnitude of the σ-hole) is larger for neutral C-vertex halogen-substituted carboranes than for their organic counterparts. Furthermore, the desolvation penalties of substituted carboranes are smaller than those of the corresponding organic compounds, which should further strengthen the halogen bonds of the former in the solvent. It is predicted that substituted carboranes have the potential to form stronger halogen bonds than comparable aromatic hydrocarbons, which will be even more pronounced in the medium. This theoretical study thus lays ground for the rational engineering of halogen bonding in inorganic crystals as well as in biomolecular complexes.

The Effect of Halogen-to-Hydrogen Bond Substitution on Human Aldose Reductase Inhibition

The effect of halogen-to-hydrogen bond substitution on the binding energetics and biological activity of a human aldose reductase inhibitor has been studied using X-ray crystallography, IC50 measurements, advanced binding free energy calculations, and simulations. The replacement of Br or I atoms by an amine (NH2) group has not induced changes in the original geometry of the complex, which made it possible to study the isolated features of selected noncovalent interactions in a biomolecular complex.

Correction: The strength and directionality of a halogen bond are co-determined by the magnitude and size of the σ-hole

Correction for 'The strength and directionality of a halogen bond are co-determined by the magnitude and size of the σ-hole' by Michal Kolář et al., Phys. Chem. Chem. Phys., 2014, 16, 9987–9996.

Structure and energetics of the anisole–Arn (n = 1, 2, 3) complexes: high-resolution resonant two-photon and threshold ionization experiments, and quantum chemical calculations

We present a concerted experimental and theoretical study of the anisole⋯Ar_n complexes with n = 1–3.

A comparison of ab initio quantum-mechanical and experimental D0 binding energies of eleven H-bonded and eleven dispersion-bound complexes

Dissociation energies (D₀) of 11 H-bonded and 11 dispersion-bound complexes were calculated as the sum of interaction energies and the change of zero-point vibrational energies (ΔZPVE).

Influence of hydrophobic residues on the binding of CB[7] toward diammonium ions of common ammonium⋯ammonium distance

We report that diammonium ion guests with common N⁺⋯N⁺ spacings but different core hydrophobic moieties exhibit remarkably different affinities toward CB[7].

Malonate-based inhibitors of mammalian serine racemase: kinetic characterization and structure-based computational study

Overactivation of NMDA receptors has been implicated in various neuropathological conditions, including brain ischaemia, neurodegenerative disorders and epilepsy. Production of d-serine, an NMDA receptor co-agonist, from l-serine is catalyzed in vivo by the pyridoxal-5'-phosphate (PLP)-dependent enzyme serine racemase. Specific inhibition of this enzyme has been proposed as a promising strategy for treatment of neurological conditions caused by NMDA receptor dysfunction. Here we present the synthesis and activity analysis of a series of malonate-based inhibitors of mouse serine racemase (mSR). The compounds possessed IC50 values ranging from 40 ± 11 mM for 2,2-bis(hydroxymethyl)malonate down to 57 ± 1 μM for 2,2-dichloromalonate, the most effective competitive mSR inhibitor known to date. The structure-activity relationship of the whole series in the human orthologue (hSR) was interpreted using Glide docking, WaterMap analysis of hydration and quantum mechanical calculations based on the X-ray structure of the hSR/malonate complex. Docking into the hSR active site with three thermodynamically favourable water molecules was able to discern qualitatively between good and weak inhibitors. Further improvement in ranking was obtained using advanced PM6-D3H4X/COSMO semiempirical quantum mechanics-based scoring which distinguished between the compounds with IC50 better/worse than 2 mM. We have thus not only found a new potent hSR inhibitor but also worked out a computer-assisted protocol to rationalize the binding affinity which will thus aid in search for more effective SR inhibitors. Novel, potent hSR inhibitors may represent interesting research tools as well as drug candidates for treatment of diseases associated with NMDA receptor overactivation.

Carborane-based carbonic anhydrase inhibitors: insight into CAII/CAIX specificity from a high-resolution crystal structure, modeling, and quantum chemical calculations

Carborane-based compounds are promising lead structures for development of inhibitors of carbonic anhydrases (CAs). Here, we report structural and computational analysis applicable to structure-based design of carborane compounds with selectivity toward the cancer-specific CAIX isoenzyme. We determined the crystal structure of CAII in complex with 1-methylenesulfamide-1,2-dicarba-closo-dodecaborane at 1.0 Å resolution and used this structure to model the 1-methylenesulfamide-1,2-dicarba-closo-dodecaborane interactions with CAIX. A virtual glycine scan revealed the contributions of individual residues to the energy of binding of 1-methylenesulfamide-1,2-dicarba-closo-dodecaborane to CAII and CAIX, respectively.

Halogen bonds in crystal TTF derivatives: an ab initio quantum mechanical study

The stabilisation energies of five ionic and neutral organic crystal structures containing various halogen bonds (I···I, Br···Br, I···Br, I···S and Br···S) were calculated using the DFT-D3 method (B97D/def2-QZVP). Besides them, the ionic I3(-)···I2 and neutral I2···I2, complexes (in the crystal geometries) were also studied. The nature of the bonds was deduced from the electrostatic potential evaluated for all subsystems. In almost all the cases, the σ-hole was positive; it was negative only for the ionic I3(-) system (although more positive than the respective belt value). The strongest halogen bonds were those that involved iodine as a halogen-bond donor and acceptor. Among ionic X···I3(-) and neutral X···I2 and X···Y dimers, the neutral X···I2 complexes were, surprisingly enough, the most stable; the highest stabilisation energy of 13.8 kcal mol(-1) was found for the I2···1,3-dithiole-2-thione-4-carboxylic acid complex. The stabilisation energies of the ionic I3(-)···I2 and neutral I2···1,3-dithiole-2-thione-4-carboxylic acid (20.2 and 20.42 kcal mol(-1), respectively) complexes are very high, which is explained by the favourable geometrical arrangement, allowing the formation of a strong halogen bond. An I···I halogen bond also exists in the neutral I2···I2 complex, having only moderate stabilisation energy (3.9 kcal mol(-1)). This stabilisation energy was, however, shown to be close to that in the optimal gas-phase L-shaped I2···I2 complex. In all the cases, the dispersion energy is important and comparable to electrostatic energy. Only in strong halogen bonds (e.g. I3(-)···I2), the electrostatic energy becomes dominant.

Ab Initio Quantum Mechanical Description of Noncovalent Interactions at Its Limits: Approaching the Experimental Dissociation Energy of the HF Dimer

Hydrogen fluoride dimer is a perfect model system for studying hydrogen bonding. Its size makes it possible to apply the most advanced theoretical methods available, yet it is a full-featured complex of molecules with nontrivial electronic structure and dynamic properties. Moreover, the dissociation energy of the HF dimer has been measured experimentally with an unparalleled accuracy of ±1 cm(-1)(Bohac et al. J. Chem. Phys. 1992, 9, 6681). In this work, we attempt to reproduce it by purely ab initio means, using advanced quantum-mechanical computational methods free of any empiricism. The purpose of this study is to demonstrate the capabilities of today's computational chemistry and to point out its limitations by identifying the contributions that introduce the largest uncertainty into the result. The dissociation energy is calculated using a composite scheme including large basis set CCSD(T) calculations, contributions of higher excitations up to CCSDTQ, relativistic and diagonal Born-Oppenheimer corrections and anharmonic vibrational calculations. The error of the calculated dissociation energy is 0.07 kcal/mol (25 cm(-1), 2.5%) when compared to the experiment. The major part of this error can be attributed to the inaccuracy of the calculations of the zero-point vibrational energy.

The relative roles of electrostatics and dispersion in the stabilization of halogen bonds

In this work we highlight recent work aimed at the characterization of halogen bonds. Here we discuss the origins of the σ-hole, the modulation of halogen bond strength by changing of neighboring chemical groups (i.e. halogen bond tuning), the performance of various computational methods in treating halogen bonds, and the strength and character of the halogen bond, the dihalogen bond, and two hydrogen bonds in bromomethanol dimers (which serve as model complexes) are compared. Symmetry adapted perturbation theory analysis of halogen bonding complexes indicates that halogen bonds strongly depend on both dispersion and electrostatics. The electrostatic interaction that occurs between the halogen σ-hole and the electronegative halogen bond donor is responsible for the high degree of directionality exhibited by halogen bonds. Because these noncovalent interactions have a strong dispersion component, it is important that the computational method used to treat a halogen bonding system be chosen very carefully, with correlated methods (such as CCSD(T)) being optimal. It is also noted here that most forcefield-based molecular mechanics methods do not describe the halogen σ-hole, and thus are not suitable for treating systems with halogen bonds. Recent attempts to improve the molecular mechanics description of halogen bonds are also discussed.