A knowledge-based halogen bonding scoring function for predicting protein-ligand interactions

Halogen bonding, a non-covalent interaction between the halogen σ-hole and Lewis bases, could not be properly characterized by majority of current scoring functions. In this study, a knowledge-based halogen bonding scoring function, termed XBPMF, was developed by an iterative method for predicting protein-ligand interactions. Three sets of pairwise potentials were derived from two training sets of protein-ligand complexes from the Protein Data Bank. It was found that two-dimensional pairwise potentials could characterize appropriately the distance and angle profiles of halogen bonding, which is superior to one-dimensional pairwise potentials. With comparison to six widely used scoring functions, XBPMF was evaluated to have moderate power for predicting protein-ligand interactions in terms of “docking power”, “ranking power” and “scoring power”. Especially, it has a rather satisfactory performance for the systems with typical halogen bonds. To the best of our knowledge, XBPMF is the first halogen bonding scoring function that is not dependent on any dummy atom, and is practical for high-throughput virtual screening. Therefore, this scoring function should be useful for the study and application of halogen bonding interactions like molecular docking and lead optimization.

Competition between halogen, dihalogen and hydrogen bonds in bromo- and iodomethanol dimers

O-H…X and O-H…O H-bonds as well as C-X…X dihalogen and C-X…O halogen bonds have been investigated in halomethanol dimers (bromomethanol dimer, iodomethanol dimer, difluorobromomethanol…bromomethanol complex and difluoroiodomethanol…iodomethanol complex). Structures of all complexes were optimized at the counterpoise-corrected MP2/cc-pVTZ level and single-point energies were calculated at the CCSD(T)/aug-cc-pVTZ level. Energy decomposition for the bromomethanol dimer complex was performed using the DFT-SAPT method based on the aug-cc-pVTZ basis set. OH…O and OH…X H-bonds are systematically the strongest in all complexes investigated, with the former being the strongest bond. Halogen and dihalogen bonds, being of comparable strength, are weaker than both H-bonds but are still significant. The strongest bonds were found in the difluoroiodomethanol…iodomethanol complex, where the O-H…O H-bond exceeds 7 kcal mol^-1, and the halogen and dihalogen bonds exceed 2.5 and 2.3 kcal mol^-1, respectively. Electrostatic energy is dominant for H-bonded structures, in halogen bonded structures electrostatic and dispersion energies are comparable, and, finally, for dihalogen structures the dispersion energy is clearly dominant.

Insights into the strength and nature of carbene···halogen bond interactions: a theoretical perspective

Halogen-bonding, a noncovalent interaction between a halogen atom X in one molecule and a negative site in another, plays critical roles in fields as diverse as molecular biology, drug design and material engineering. In this work, we have examined the strength and origin of halogen bonds between carbene CH₂ and XCCY molecules, where X?=?Cl, Br, I, and Y?=?H, F, COF, COOH, CF₃, NO₂, CN, NH₂, CH₃, OH. These calculations have been carried out using M06-2X, MP2 and CCSD(T) methods, through analyses of surface electrostatic potentials V _S(r) and intermolecular interaction energies. Not surprisingly, the strength of the halogen bonds in the CH₂···XCCY complexes depend on the polarizability of the halogen X and the electron-withdrawing power of the Y group. It is revealed that for a given carbene···X interaction, the electrostatic term is slightly larger (i.e., more negative) than the dispersion term. Comparing the data for the chlorine, bromine and iodine substituted CH₂···XCCY systems, it can be seen that both the polarization and dispersion components of the interaction energy increase with increasing halogen size. One can see that increasing the size and positive nature of a halogen’s σ-hole markedly enhances the electrostatic contribution of the halogen-bonding interaction.

A theoretical investigation of the characteristics of hydrogen/halogen bonding interactions in dibromo-nitroaniline

In this work, computations of density functional theory (DFT) were carried out to investigate the nature of interactions in solid 2,6-dibromo-4-nitroaniline (DBNA). This system was selected to mimic the hydrogen/halogen bonding found within crystal structures as well as within biological molecules. DFT (M06-2X/6-311++G**) calculations indicated that the binding energies for different of interactions lie in the range between ?1.66 and ?9.77 kcal mol^?1. The quantum theory of atoms in molecules (QTAIM) was applied to provide more insight into the nature of these interactions. Symmetry-adapted perturbation theory (SAPT) analysis indicated that stability of the Br···Br halogen bonds is predicted to be attributable mainly to dispersion, while electrostatic forces, which have been widely believed to be responsible for these types of interactions, play a smaller role. Our results indicate that, for those nuclei participating in hydrogen/halogen bonding interactions, nuclear quadrupole resonance parameters exhibit considerable changes on going from the isolated molecule model to crystalline DBNA.

A theoretical study on the halogen bonding interactions of C6F5I with a series of group 10 metal monohalides

The halogen bonding interactions between C₆F₅I and a series of transition metal monohalides trans-[M(X)(2-C₅NF₄)-(PR₃)₂] (M = Ni, Pd, Pt; X = F, Cl, Br; R = Me, Cy) have been studied with quantum chemical calculations. Optimized geometries of the halogen bonding complexes indicate that angles C₁-I···X are basically linear (178–180°) and angles I···X-M mainly range from 90 to 150°. The strength of these metal-influenced halogen bonds alters with different metal centers, metal-bound halogen atoms and the substitutes on phosphine ligands. Electrostatic potential and natural bond orbital analysis show that both of the electrostatic and orbital interactions make a contribution to the formation of halogen bonds, while the electrostatic term plays a dominant role. AIM analysis suggests that, for trans-[M(F)(2-C₅NF₄)-(PR₃)₂] (M = Ni, Pd, Pt) monomers, the formed halogen bonding complexes are stabilized by local concentration of the charge of intermediate character, while for the metal monomers containing chlorine and bromine, a typical closed-shell interaction exist. These results prove that the structures and geometries of these halogen bonding complexes can be tuned by changing the halogen atoms and metal centers, which may provide useful information for the design and synthesis of new functional materials.

How do halogen bonds (S–O⋯I,N–O⋯I and C–O⋯I) and halogen–halogen contacts (C–I⋯I–C,C–F⋯F–C) subsist in crystal structures? A quantum chemical insight

Thirteen X-ray crystal structures containing various non-covalent interactions such as halogen bonds, halogen–halogen contacts and hydrogen bonds (I?N, I?F, I?I, F?F, I?H and F?H) were considered and investigated using the DFT-D3 method (B97D/def2-QZVP). The interaction energies were calculated at MO62X/def2-QZVP and MP2/aug-cc-pvDZ level of theories. The higher interaction and dispersion energies (2nd crystal) of ?9.58 kcal mol^?1 and ?7.10 kcal mol^?1 observed for 1,4-di-iodotetrafluorobenzene bis [bis (2-phenylethyl) sulfoxide] structure indicates the most stable geometrical arrangement in the crystal packing. The electrostatic potential values calculated for all crystal structures have a positive σ-hole, which aids understanding of the nature of σ-hole bonds. The significance of the existence of halogen bonds in crystal packing environments was authenticated by replacing iodine atoms by bromine and chlorine atoms. Nucleus independent chemical shift analysis reported on the resonance contribution to the interaction energies of halogen bonds and halogen–halogen contacts. Hirshfeld surface analysis and topological analysis (atoms in molecules) were carried out to analyze the occurrence and strength of all non-covalent interactions. These analyses revealed that halogen bond interactions were more dominant than hydrogen bonding interactions in these crystal structures.

Open image in new window

Graphical Abstract Molecluar structure of 1,4-Di-iodotetrafluorobenzene bis(thianthrene 5-oxide) moelcule and its corresponding molecular electrostatic potential map for the view of σ-hole.

     

Simulations on the possibility of formation of complexes between fluorouracil drug and cucurbit[n]urils: ab initio van der Waals DFT study

The binding geometry of fluorouracil/cucurbit[n]urils (CB[n]s) complexes with n?=?5–8 is investigated using the first-principles van der Waals density functional (vdW-DF) method, involving full geometry optimization. Such host-guest complexes are typically calculated using conventional DFT method, which significantly underestimates non-local dispersion forces (or vdW contributions) and therefore affects interactions between respected entities. We address here the role of vdW forces for the fluorouracil and CB[n]s molecules which can form directional hydrogen bonds with each other. It was found that the inclusion of dispersion interactions significantly affects the host-guest binding properties and the hydrogen bonding between the molecules provides the main binding mechanism, while results in the same geometries for the considered complexes. The 0.84 eV binding energy, for the thermodynamically favorable state, reveals that the interaction of fluorouracil with CB[n]s is an exothermic interaction and typical for strong hydrogen bonding. Furthermore, we have investigated the binding nature of these host-guest systems in aqueous solution with ab initio MD simulations adopting vdW-DF method. These findings afford evidence for the applicability of the vdW-DF approach and provide a realistic benchmark for the investigation of the host-guest complexes.

ETS-NOCV description of σ-hole bonding

The ETS-NOCV analysis was applied to describe the σ-hole in a systematic way in a series of halogen compounds, CF₃-X (X?=?I, Br, Cl, F), CH₃I, and C(CH3)_nH_3-n-I (n?=?1,2,3), as well as for the example germanium-based systems. GeXH₃, X?=?F, Cl, H. Further, the ETS-NOCV analysis was used to characterize bonding with ammonia for these systems. The results show that the dominating contribution to the deformation density, Δρ ₁ , exhibits the negative-value area with a minimum, corresponding to σ-hole. The “size” (spatial extension of negative value) and “depth” (minium value) of the σ-hole varies for different X in CF₃-X, and is influenced by the carbon substituents (fluorine atoms, hydrogen atoms, methyl groups). The size and depth of σ-hole decreases in the order: I, Br, Cl, F in CF₃-X. In CH₃-I and C(CH3)_nH_3-n-I, compared to CF₃-I, introduction of hydrogen atoms and their subsequent replacements by methyl groups lead to the systematic decrease in the σ-hole size and depth. The ETS-NOCV σ-hole picture is consistent with the existence the positive MEP area at the extension of σ-hole generating bond. Finally, the NOCV deformation density contours as well as by the ETS orbital-interaction energy indicate that the σ-hole-based bond with ammonia contains a degree of covalent contribution. In all analyzed systems, it was found that the electrostatic energy is approximately two times larger than the orbital-interaction term, confirming the indisputable role of the electrostatic stabilization in halogen bonding and σ-hole bonding.

Cooperativity between fluorine-centered halogen bonds: investigation of substituent effects

This article analyzes the substitution effects on cooperativity between fluorin-centered halogen bonds in NCF?·?·?·?NCF?·?·?·?NCX and CNF?·?·?·?CNF?·?·?·?CNX complexes, where X?=?H, F, Cl, CN, OH, and NH₂. These effects are investigated theoretically in terms of geometric and energetic features of the complexes, which are computed by ab initio methods. The topological analysis, based on the quantum theory of atoms in molecules (QTAIM), is used to characterize the interactions and analyze their enhancement with varying electron density at bond critical points. It is found that the complexes with electron-donating groups exhibit a strong cooperativity, while a much weaker cooperativity occurs in the NCF?·?·?·?NCF?·?·?·?NCCN and CNF?·?·?·?CNF?·?·?·?CNCN trimers. An excellent correlation is found between the cooperative energy in the ternary complexes and the calculated three-body interaction energies. The energy decomposition analysis (EDA) indicates that the electrostatic and dispersion effects play a main role in the cooperativity of fluorine-centered halogen bonding.

Chalcogen- and halogen-bonds involving SX<Subscript>2</Subscript> (X = F,Cl, and Br) with formaldehyde

The capacity of SX₂ (X = F, Cl, and Br) to engage in different kinds of noncovalent bonds was investigated by ab initio calculations. SCl₂ (SBr₂) has two σ-holes upon extension of Cl (Br)?S bonds, and two σ-holes upon extension of S?Cl (Br) bonds. SF₂ contains only two σ-holes upon extension of the F?S bond. Consequently, SCl₂ and SBr₂ form chalcogen and halogen bonds with the electron donor H₂CO while SF₂ forms only a chalcogen bond, i.e., no F···O halogen bond was found in the SF₂:H₂CO complex. The S···O chalcogen bond between SF₂ and H₂CO is the strongest, while the strongest halogen bond is Br···O between SBr₂ and H₂CO. The nature of these two types of noncovalent interaction was probed by a variety of methods, including molecular electrostatic potentials, QTAIM, energy decomposition, and electron density shift maps. Termolecular complexes X₂S···H₂CO···SX′₂ (X = F, Cl, Br, and X′ = Cl, Br) were constructed to study the interplay between chalcogen bonds and halogen bonds. All these complexes contained S···O and Cl (Br)···O bonds, with longer intermolecular distances, smaller values of electron density, and more positive three-body interaction energies, indicating negative cooperativity between the chalcogen bond and the halogen bond. In addition, for all complexes studied, interactions involving chalcogen bonds were more favorable than those involving halogen bonds.

Open image in new window

Graphical Abstract Molecular electrostatic potential and contour map of the Laplacian of the electron density in Cl₂S···H₂CO···SCl₂ complex

     

Theoretical description of halogen bonding – an insight based on the natural orbitals for chemical valence combined with the extended-transition-state method (ETS-NOCV)

In the present study we have characterized the halogen bonding in selected molecules H₃N–ICF₃ (1-NH ₃ ), (PH₃)₂C–ICF₃ (1-CPH ₃ ), C₃H₇Br–(IN₂H₂C₃)₂C₆H₄ (2-Br), H₂–(IN₂H₂C₃)₂C₆H₄ (2-H ₂ ) and Cl–(IC₆F₅)₂C₇H₁₀N₂O₅ (3-Cl), containing from one halogen bond (1-NH ₃ , 1-CPH ₃ ) up to four connections in 3-Cl (the two Cl–HN and two Cl–I), based on recently proposed ETS-NOCV analysis. It was found based on the NOCV-deformation density components that the halogen bonding C–X ^… B (X-halogen atom, B-Lewis base), contains a large degree of covalent contribution (the charge transfer to X ^… B inter-atomic region) supported further by the electron donation from base atom B to the empty σ*(C–X) orbital. Such charge transfers can be of similar importance compared to the electrostatic stabilization. Further, the covalent part of halogen bonding is due to the presence of σ-hole at outer part of halogen atom (X). ETS-NOCV approach allowed to visualize formation of the σ-hole at iodine atom of CF₃I molecule. It has also been demonstrated that strongly electrophilic halogen bond donor, [C₆H₄(C₃H₂N₂I)₂][OTf]₂, can activate chemically inert isopropyl bromide (2-Br) moiety via formation of Br–I bonding and bind the hydrogen molecule (2-H ₂ ). Finally, ETS-NOCV analysis performed for 3-Cl leads to the conclusion that, in terms of the orbital-interaction component, the strength of halogen (Cl–I) bond is roughly three times more important than the hydrogen bonding (Cl–HN).

Theoretical study of the triangular bonding complex formed by carbon tetrabromide, a halide, and a solvent molecule in the gas phase

MP2(full)/aug-cc-pVDZ(-PP) computations predict that new triangular bonding complexes (where X^? is a halide and H–C refers to a protic solvent molecule) consist of one halogen bond and two hydrogen bonds in the gas phase. Carbon tetrabromide acts as the donor in the halogen bond, while it acts as an acceptor in the hydrogen bond. The halide (which commonly acts as an acceptor) can interact with both carbon tetrabromide and solvent molecule (CH₃CN, CH₂Cl₂, CHCl₃) to form a halogen bond and a hydrogen bond, respectively. The strength of the halogen bond obeys the order CBr₄???Cl^? > CBr₄???Br^? > CBr₄???I^?. For the hydrogen bonds formed between various halides and the same solvent molecule, the strength of the hydrogen bond obeys the order C-H???Cl^? > C-H???Br^? > C-H???I^?. For the hydrogen bonds formed between the same halide and various solvent molecules, the interaction strength is proportional to the acidity of the hydrogen in the solvent molecule. The diminutive effect is present between the hydrogen bonds and the halogen bond in chlorine and bromine triangular bonding complexes. Complexes containing iodide ion show weak cooperative effects.

Interplay between halogen bonds and hydrogen bonds in OH/SH···HOX···HY (X = Cl, Br; Y = F, Cl, Br) complexes

The character of the cooperativity between the HOX···OH/SH halogen bond (XB) and the Y―H···(H)OX hydrogen bond (HB) in OH/SH···HOX···HY (X = Cl, Br; Y = F, Cl, Br) complexes has been investigated by means of second-order Møller?Plesset perturbation theory (MP2) calculations and “quantum theory of atoms in molecules” (QTAIM) studies. The geometries of the complexes have been determined from the most negative electrostatic potentials (V _S,min) and the most positive electrostatic potentials (V _S,max) on the electron density contours of the individual species. The greater the V _S,max values of HY, the larger the interaction energies of halogen-bonded HOX···OH/SH in the termolecular complexes, indicating that the ability of cooperative effect of hydrogen bond on halogen bond are determined by V _S,max of HY. The interaction energies, binding distances, infrared vibrational frequencies, and electron densities ρ at the BCPs of the hydrogen bonds and halogen bonds prove that there is positive cooperativity between these bonds. The potentiation of hydrogen bonds on halogen bonds is greater than that of halogen bonds on hydrogen bonds. QTAIM studies have shown that the halogen bonds and hydrogen bonds are closed-shell noncovalent interactions, and both have greater electrostatic character in the termolecular species compared with the bimolecular species.

Theoretical study on cooperative effects between X⋯N and X⋯Carbene halogen bonds (X = F,Cl,Br and I)

Quantum chemical calculations are performed to study the interplay between halogen?nitrogen and halogen?carbene interactions in NCX?NCX?CH₂ complexes, where X?=?F, Cl, Br and I. Molecular geometries and interaction energies of dyads and triads are investigated at the MP2/aug-cc-pVTZ level of theory. It is found that the X?N and X?C_carbene interaction energies in the triads are larger than those in the dyads, indicating that both the halogen bonding interactions are enhanced. The estimated values of cooperative energy E _coop are all negative with much larger E _coop in absolute value for the systems including iodine. The nature of halogen bond interactions of the complexes is analyzed using parameters derived from the quantum theory atoms in molecules methodology and energy decomposition analysis.

Revealing substitution effects on the strength and nature of halogen-hydride interactions: a theoretical study

A quantum chemistry study was carried out to investigate the strength and nature of halogen bond interactions in HXeH···XCCY complexes, where X = Cl, Br and Y = H, F, Cl, Br, CN, NC, C₂H, CH₃, OH, SH, NH₂. Examination of the electrostatic potentials V(r) of the XCCY molecules reveals that the addition of substituents has a significant effect upon the most positive electrostatic potential on the surface of the interacting halogen atom. We found that the magnitude of atomic charges and multipole moments depends upon the halogen atom X and is rather sensitive to the electron-withdrawing/donating power of the remainder of the molecule. An excellent correlation was found between the most positive electrostatic potentials on the halogen atom and the interaction energies. For either HXeH···ClCCY or HXeH···BrCCY complexes, an approximate linear correlation between the interaction energies and halogens multipole moments are established, indicating that the electrostatic and polarization interactions are responsible for the stability of the complexes. According to energy decomposition analysis, it is revealed that the electrostatic interactions are the major source of the attraction in the HXeH···XCCY complexes. Furthermore, the changes in the electrostatic term are mainly responsible for the dependence of interaction energy on the halogen atom.

The competition of C-X⋯O=P halogen bond and π-hole⋯O=P bond between halopentafluorobenzenes C6F5X (X=F, Cl, Br, I) and triethylphosphine oxide

Calculation predicted the interacting forms of halopentafluorobenzene C₆F₅X (X=F, Cl, Br, I) with triethylphosphine oxide which is biologically interested and easily detected by ³¹P NMR. The interaction energy and geometric parameters of resultant halogen or π-hole bonding complexes were estimated and compared. Moreover, the bonding constants were determined by ³¹P NMR. Both theory and experiments indicated the C₆F₆ and C₆F₅Cl interact with triethylphosphine oxide by π-hole bonding pattern, while C₆F₅I by halogen/σ-hole bonding form. For C₆F₅Br, two interactions are comparative and should coexist competitively. The calculated interaction energies of σ-hole bonding complexes, ?5.07 kcal mol^?1 for C₆F₅Br?O=P and ?8.25 kcal mol^?1 for C₆F₅I?O=P, and π-hole bonding complexes, ?7.29 kcal mol^?1 for C₆F₆?O=P and ?7.24 kcal mol^?1 for C₆F₅Cl?O=P, are consistent with the changing tendency of bonding constants measured by ³¹P NMR, 4.37, 19.7, 2.42 and 2.23 M^?1, respectively.

Effects of bidentate coordination on the molecular properties rapta-C based complex using theoretical approach

In this work several quantum properties including the NEDA and QTAIM are computed on three models of rapta-C complexes using DFT with hybrid functional and basis set with ECP and without ECP. Several interesting correlations within the observed properties and also with the reported experimental behaviors of these complexes including their biological activities are presented. The study shows that the stability of the two complexes with bidentate ligands is associated with their high hydrogen bonding stability and existence of stronger non-covalent metal-ligand bonds. The energy decomposition analysis indicated that inter-atomic interactions in the three forms of rapta-C complexes and their stability are governed by the charge transfer term with significant contributions from polarization and electrostatic terms. The higher stability of complex 1 and 2 over 3 comes from the lower exchange repulsion and higher polarization contributions to their stability which agrees perfectly with the experimental observation. Our results provide insight into the nature of intramolecular forces that influence the structural stability of the three complexes.     

Insight on the interaction of polychlorobiphenyl with nucleic acid–base

The interaction between one polychlorobiphenyl (3,3′,4,4′,-tetrachlorobiphenyl, coded PCB77) and the four DNA nucleic acid–base is studied by means of quantum mechanics calculations in stacked conformations. It is shown that even if the intermolecular dispersion energy is the largest component of the total interaction energy, some other contributions play a non negligible role. In particular the electrostatic dipole-dipole interaction and the charge transfer from the nucleobase to the PCB are responsible for the relative orientation of the monomers in the complexes. In addition, the charge transfer tends to flatten the PCB, which could therefore intercalate more easily between DNA base pairs. From these seminal results, we predict that PCB could intercalate completely between two base pairs, preferably between Guanine:Cytosine pairs.

First-principles vdW-DF investigation on the interaction between the oxazepam molecule and C60 fullerene

The interaction between oxazepam and C₆₀ fullerene was explored using first-principles vdW-DF calculations. It was found that oxazepam binds weakly to the fullerene cage via its carbonyl group. The binding of oxazepam to C₆₀ is affected drastically by nonlocal dispersion interactions, while vdW forces affect the corresponding geometries only a little. Furthermore, aqueous solution affects the geometries of the oxazepam approaching to fullerene slightly, while oxazepam binds slightly farther away from the nanocage. The results presented provide evidence for the applicability of the vdW-DF method and serve as a practical benchmark for the investigation of host–guest interactions in biological systems.