Structures and stability of I2 and ICl complexes with pyridine: Ab initio and DFT study

Theoretical investigation of thermodynamic stability and bonding features of possible isomers of the molecular and ionic complexes of pyridine with molecular iodine and iodine monochloride IX (X = I,Cl) is presented. M06-2X DFT functional is found to provide bond distances and dissociation energies which are close to those obtained at high-level ab initio CCSD(T)/aug-cc-pvtz//CCSD/aug-cc-pvtz benchmark computations for the most stable isomers, formed via donation of a lone pair of nitrogen atom of pyridine to the iodine atom. These isomers are by 23-33 kJ mol-1 (in case of I2) and by 39-56 kJ mol-1 (in case of ICl) more stable than other molecular complexes. T-shaped π-σ* bonded isomers turn out to be energetically comparable with van der Waals bound compounds. Among the ionic isomers, structures featuring [IPy2]+ cation with I3 - or ICl2 - counterions are more stable. Oligomerization favors ionic isomers starting from the tetrameric clusters of the composition (IX)4Py4.

Substituent Effect and Its Halogen-Atom Dependence of Halogen Bonding Viewed through Electron Density Changes

Elucidating how the halogen-bonding ability and strength are controlled by the substituent effect and how this control depends on halogen atom will be essential for finely-tuned design of functionally important molecules. Here, this problem is tackled by analyzing the electron density differences/changes for variously substituted halobenzenes. It is shown that the anisotropy of the electron distribution around the halogen atom, which is an important factor for halogen-bonding ability, is not much affected by the substituent effect and rather simply depends on the halogen atom, while the partial charge on the halogen atom, which is related to the bond dipole of the C-X bond, is significantly modulated by the substituent effect and gives rise to enhancement of the electrostatic potential on the line extended from the C-X bond. The properties related to the polarization effect are also discussed.

Nature of hydrogen-bond-enhanced halogen bonding viewed through electron density changes

Elucidating the mechanism of how we can achieve fine tuning of intermolecular interaction strength will be helpful for designing functionally important molecules. In the present study, a theoretical analysis is conducted, by examining the electron density changes, for two halogen-bonding iodinated systems whose halogen-bond strengths have been considered to be enhanced by the presence of a hydrogen-bond donating group (termed hydrogen-bond-enhanced halogen bonding). It is shown that, contrary to the expectation obtained from the enhancement of electrostatic potential along the line extended from the C-I bond, the anisotropy of electron distribution on the iodine atom remains nearly the same. This means that the hydrogen bond and halogen bond contribute almost independently and additively to the enhancement of electrostatic potential, indicating the nature of this enhancement and, in a more general sense, the relationship between the strength and the extent of directionality of halogen bonding.

Halogen bond-directed self-assembly in bicomponent blends at the solid/liquid interface: Effect of the alkyl chain substitution position

The fabrication of well-organised molecular assemblies on surfaces is fundamental for the creation of functional molecular systems applicable to nanoelectronics and molecular devices. In this study, we investigated the effect...

FMODB: The World's First Database of Quantum Mechanical Calculations for Biomacromolecules Based on the Fragment Molecular Orbital Method

We developed the world's first web-based public database for the storage, management, and sharing of fragment molecular orbital (FMO) calculation data sets describing the complex interactions between biomacromolecules, named FMO Database (https://drugdesign.riken.jp/FMODB/). Each entry in the database contains relevant background information on how the data was compiled as well as the total energy of each molecular system and interfragment interaction energy (IFIE) and pair interaction energy decomposition analysis (PIEDA) values. Currently, the database contains more than 13 600 FMO calculation data sets, and a comprehensive search function implemented at the front-end. The procedure for selecting target proteins, preprocessing the experimental structures, construction of the database, and details of the database front-end were described. Then, we demonstrated a use of the FMODB by comparing IFIE value distributions of hydrogen bond, ion-pair, and XH/π interactions obtained by FMO method to those by molecular mechanics approach. From the comparison, the statistical analysis of the data provided standard reference values for the three types of interactions that will be useful for determining whether each interaction in a given system is relatively strong or weak compared to the interactions contained within the data in the FMODB. In the final part, we demonstrate the use of the database to examine the contribution of halogen atoms to the binding affinity between human cathepsin L and its inhibitors. We found that the electrostatic term derived by PIEDA greatly correlated with the binding affinities of the halogen containing cathepsin L inhibitors, indicating the importance of QM calculation for quantitative analysis of halogen interactions. Thus, the FMO calculation data in FMODB will be useful for conducting statistical analyses to drug discovery, for conducting molecular recognition studies in structural biology, and for other studies involving quantum mechanics-based interactions.

Halogen bonding interactions in the XC5H4N···YCF3 (X = CH3, H, Cl, CN, NO2; Y = Cl, Br, I) complexes

The noncovalent interactions between the σ-hole region outside the halogen atom and the nitrogen atom of pyridine and its para-substituted derivatives are the focus of this work. Based on the analyses of the electrostatic potentials, YCF₃ (Y = Cl, Br, I) act as halogen bond donors, XC₅H₄N (X = CH₃, H, Cl, CN, NO₂) act as halogen bond acceptors, and the binary halogen-bonded complexes XC₅H₄N···YCF₃ have been designed and investigated by B3LYP-D3/aug-cc-pVDZ calculations together with the aug-cc-pVDZ-PP basis set for iodine. When the halogen bond acceptor remains unchanged, the interactions between C₅H₅N and YCF₃ (Y = Cl, Br, I) increase with the order of Y = Cl, Br, and I. When the halogen donor ICF₃ is fixed, the halogen bonding interactions decrease along the sequence of X = CH₃, H, Cl, CN, NO₂. Therefore, the halogen bond of the CH₃C₅H₄N···ICF₃ complex is the strongest. The interactions between Lewis acid YCF₃ (Y = Cl, Br, I) and pyridine and para-substituted pyridine are closed-shell and noncovalent interactions. On the one hand, when the halogen bond acceptor XC₅H₄N is fixed, with the increase of halogen atomic number, the strength of halogen bond increases; on the other hand, when the halogen bond donor ICF₃ is fixed, as the electron-withdrawing ability of the electron-withdrawing group (X) increases, the halogen bond gradually weakens.

Dynamic host-guest behavior in halogen-bonded two-dimensional molecular networks investigated by scanning tunneling microscopy at the solid/liquid interface

The fabrication of supramolecularly engineered two-dimensional (2D) networks using simple molecular building blocks is an effective means for studying host-guest chemistry at surfaces toward the potential application of such systems in nanoelectronics and molecular devices. In this study, halogen-bonded molecular networks were constructed by the combination of linear halogen-bond donor and acceptor ligands, and their 2D structures at the highly oriented pyrolytic graphite/1-phenyloctane interface were studied by scanning tunneling microscopy. The bi-component blend of the molecular building blocks possessing tetradecyloxy chains formed a lozenge structure via halogen bonding. Upon the introduction of an appropriate guest molecule (e.g., coronene) into the system, the 2D structure transformed into a hexagonal array, and the central pore of this array was occupied by the guest molecules. Remarkably, the halogen bonding of the original structure was maintained after the introduction of the guest molecule. Thus, the halogen-bonded molecular networks are applicable for assembling guest species on the substrate without the requirement of the conventional rigid molecular building blocks with C 3 symmetry.

Hexagonal array formation by intermolecular halogen bonding using a binary blend of linear building blocks: STM study

Hexagonal arrays were fabricated via intermolecular halogen bonding between two linear molecular building blocks in a bicomponent blend. The substitution position of the pyridine N atom involved in the halogen bond plays an important role in the formation of the hexagonal structures.

Correlation of the partial charge-transfer and covalent nature of halogen bonding with the THz and IR spectral changes

Changes in the spectral intensities in the THz region are good probes for the non-electrostatic aspect of halogen bonding.

Modeling cooperative effects in halogen-bonded infinite linear chains

Non-additivity in noncovalent interactions is an important aspect of complex systems that can lead to stronger (cooperative) interactions when three or more molecular units influence each other. The halogen bond (XB) is a highly-directional noncovalent interaction that has been found to be cooperative. Here the strength and nature of cooperativity arising in X-bonded infinite linear chains of cyanogen halides and 4-halopyridines are investigated by means of density functional theory calculations. It is found that cyanogen halide chains are highly cooperative (up to 77%), whereas pyridines are only slightly cooperative (below 21%). It is demonstrated that XB and its non-additivity can be modeled as the sum of a local term, which depends on first nearest-neighbors only, and long-range effective dipole-dipole attractions. It is shown that the local term in cyanogen halides primarily accounts for repulsive short-range screened Coulomb interactions, whereas in 4-halopyridines such a term also includes attractive contributions, which are particularly sizeable in some elongated XB conformations. This outcome reveals differences in the nature of the XBs formed in these molecular systems. Nevertheless, it is shown that both systems behave as effective point dipoles regarding cooperative effects, at any point of the XB dissociation path. As such, these results are useful contributions for the understanding and modeling of non-additive effects of noncovalent interactions.

The15N NMR chemical shift in the characterization of weak halogen bonding in solution

We have studied the applicability of¹⁵N NMR spectroscopy in the characterization of the very weak halogen bonds of nonfluorinated halogen bond donors with a nitrogenous Lewis base in solution. The ability of the technique to detect the relative strength of iodine-, bromine- and chlorine-centered halogen bonds, as well as solvent and substituent effects was evaluated. Whereas computations on the DFT level indicate that¹⁵N NMR chemical shifts reflect the diamagnetic deshielding associated with the formation of a weak halogen bond, the experimentally observed chemical shift differences were on the edge of detectability due to the low molar fraction of halogen-bonded complexes in solution. The formation of the analogous yet stronger hydrogen bond of phenols have induced approximately ten times larger chemical shift changes, and could be detected and correlated to the electronic properties of substituents of the hydrogen bond donors. Overall,¹⁵N NMR is shown to be a suitable tool for the characterization of comparably strong secondary interactions in solution, but not sufficiently accurate for the detection of the formation of thermodynamically labile, weak halogen bonded complexes.

Neutral iodotriazoles as scaffolds for stable halogen-bonded assemblies in solution

The halogen bond (XB) donor properties of neutral 1,4-diaryl-5-iodo-1,2,3-triazoles are explored using a combination of computational and experimental results and are shown to be competitive in halogen bonding efficiency with the classic pentafluoroiodobenzene XB donor. The SNAr reactivity of these donors permits the facile assembly of an iodotriazole functionalised with a 3-oxypyridine XB acceptor, thus generating a molecular scaffold capable of undergoing dimerisation through the formation of two halogen bonds. The formation of this halogen-bonded dimer is demonstrated by 1H and DOSY NMR experiments and a plausible structure generated using DFT calculations.

Magnitude and Directionality of Halogen Bond of Benzene with C6F5X, C6H5X, and CF3X (X = I, Br, Cl, and F)

Geometries of benzene complexes with C6F5X, C6H5X, and CF3X (X is I, Br, Cl, and F) were optimized, and their interaction energies were evaluated. The CCSD(T) interaction energies at the basis set limit (Eint) of C6F5X (X is I, Br, Cl, and F) with benzene were -3.24, -2.88, -2.31, and -0.92 kcal mol(-1). Eint of C6H5X (X is I, Br, and Cl) with benzene were -2.31, -1.97, and -1.48 kcal mol(-1). The fluorination of halobenzenes slightly enhances the attraction. Eint of CF3X (X is I, Br, Cl, and F) with benzene (-3.11, -2.74, -2.22, and -0.71 kcal mol(-1)) were very close to Eint of corresponding C6F5X with benzene. In contrast to the halogen bond of iodine and bromine with pyridine (n-type halogen bond acceptor) where the main cause of the attraction is the electrostatic interactions, that of halogen bond with benzene (p-type acceptor) is dispersion interaction. In the halogen bonds with p-type acceptors (halogen-π interactions), the electrostatic interactions and induction interactions are small. The overall orbital-orbital interactions are repulsive. The directionality of halogen bonds with p-type acceptors is very weak, owing to the weak electrostatic interactions, in contrast to the strong directionality of the halogen bonds with n-type acceptors and hydrogen bonds.

Roles of the scalar and vector components of the solvation effects on the vibrational properties of hydrogen- or halogen-bond accepting stretching modes

Solvation-induced vibrational frequency shifts and infrared (IR) intensity changes of the hydrogen- or halogen-bond accepting stretching modes, especially their dependence on the angular position of the hydrogen- or halogen-bond donating molecule, are examined theoretically. Calculations are carried out for some modes of hydrogen- or halogen-bonding molecular complexes, including the S[double bond, length as m-dash]O stretch of dimethyl sulfoxide-(13)C2H2O, the C[triple bond, length as m-dash]N stretch of acetonitrileH2O, and the amide I' mode of the N-methylacetamide-d1BrNC 1 : 1 complex. It is shown that, in all the example cases dealt with in this study, the frequency shift depends rather strongly on the hydrogen- or halogen-bond angle (e.g., S[double bond, length as m-dash]OH angle), with a larger low-frequency shift as the hydrogen or halogen bond becomes more bent, indicating the generality of the results obtained for the amide I' mode of the N-methylacetamide-d1(2)H2O 1 : 1 complex in a previous study. Contrary to our vague expectation, the frequency shift is not well correlated to the hydrogen- or halogen-bond distance or strength, but nevertheless, it is well reproduced by an electrostatic interaction model if it is carefully constructed by considering the scalar and vector components separately in a reasonable way. On the basis of this electrostatic interaction model, the reason why our vague expectation is not realized is clarified, and a unified understanding is achieved on the hydration-induced high-frequency shift of the C[triple bond, length as m-dash]N stretch and the low-frequency shifts of the S[double bond, length as m-dash]O stretch and amide I'. With regard to the IR intensity, it is shown that, in some of the example cases, it also has rather strong angular position dependence. The mechanism of the IR intensity changes is estimated by analyzing the dipole derivative vector, especially its angular relation with the hydrogen or halogen bond.

Characteristic redshift and intensity enhancement as far-IR fingerprints of the halogen bond involving aromatic donors

Characteristic redshift and intensity enhancement of the C–I stretching band have been proven to be distinct signatures of the halogen bond involving iodopentafluorobenzene.

Halogen bonding in drug-like molecules: a computational and systematic study of the substituent effect

Systematic study of the effect of fourteen chemical groups at the ortho, para and meta positions of NMA⋯halobenzene complexes showed a significant influence on halogen bonding, and also non-additive effects. A comprehensive description is reported.