Source of Cooperativity in Halogen-Bonded Haloamine Tetramers

Properties and Stabilities of Cyclic and Open Chains of Halogen Bonds

Open and cyclic chains from two to eight units of ICl and IF are constructed and examined by density functional theory (DFT) calculations. These chains contain either I···I or I···X halogen bonds (XBs) where X refers to Cl or F. The closed rings are more stable than the open chains due to the presence of an additional XB and enhanced cooperativity. This pattern is true even for most trimers where there is sizable geometric distortion in the rings. I···F rings are generally more stable than the corresponding I···I cycles as the I···F bond is stronger than I···I even in the simple dimer. However, I···I rings are comparable in energy to I···Cl. It is possible to construct I···I rings of at least as large as eight units, which are held together exclusively by XBs. On the other hand, the maximum possible size of I···X rings is 6. Red shifts are observed in the I-X stretching frequency bands, which magnify as the chain, both cyclic and open, grows longer. The NMR chemical shielding of the I atoms increases for I···I chains but diminishes when I···Cl bonds are present.

Cooperative Self-Assembly in Linear Chains Based on Halogen Bonds

Cooperative properties of halogen bonds were investigated with computational experiments based on dispersion-corrected relativistic density functional theory. The bonding mechanism in linear chains of cyanogen halide (X-CN), halocyanoacetylene (X-CC-CN), and 4-halobenzonitrile (X-C6 H4 -CN) were examined for X = H, Cl, Br, and I. Our energy decomposition and Kohn-Sham molecular-orbital analyses revealed the bonding mechanism of the studied systems. Cyanogen halide and halocyanoacetylene chains possess an extra stabilizing effect with increasing chain size, whereas the 4-halobenzonitrile chains do not. This cooperativity can be traced back to charge separation within the σ-electronic system by charge-transfer between the lone-pair orbital of the nitrogen (σHOMO ) on one unit and the acceptor orbital of the C-X (σ*LUMO ) on the adjacent unit. As such, the HOMO-LUMO gap in the σ-system decreases, and the cooperativity increases with chain length revealing the similarity in the bonding mechanisms of hydrogen and halogen bonds.

Lack of Cooperativity in the Triangular X3 Halogen-Bonded Synthon?

We have investigated 44 crystal structures, found in the Cambridge Structural Database, containing the X₃ synthon (where X = Cl, Br, I) in order to verify whether three type II halogen-halogen contacts forming the synthon exhibit cooperativity. A hypothesis that this triangular halogen-bonded motif is stabilized by cooperative effects is postulated on the basis of structural data. However, theoretical investigations of simplified model systems in which the X₃ motif is present demonstrate that weak synergy occurs only in the case of the I₃ motif. In the present paper we computationally investigate crystal structures in which the X₃ synthon is present, including halomesitylene structures, that are usually described as being additionally stabilized by a synergic interaction. Our computations find no cooperativity for halomesitylene trimers containing the X₃ motif. Only in the case of two other structures containing the I₃ synthon a very weak or weak synergy, i.e. the cooperative effect being stronger than -0.40 kcal mol^-1, is found. The crystal structure of iodoform has the most pronounced cooperativity of all investigated systems, amounting to about 10% of the total interaction energy.

Halogen-bonded haloamine trimers - modelling the X3 synthon

Halogen-halogen bonded haloamine trimers serve as model systems for the X3 synthon present in numerous crystal structures and in two-dimensional self-assembled nanoarchitectures. Halogen bonds forming the synthon are often considered to display cooperativity. Synergy effects were previously found for halogen-halogen bonded bromoamine and iodoamine tetramers. In the present study comparison between haloamine cyclic trimers and tetramers is made. The cooperativity occurring in bromoamine and iodoamine clusters is significantly weaker in the case of the trimers. The present study demonstrates that the bromoamine and iodoamine trimers display weaker cooperativity due to a smaller number of synergy components in comparison to the corresponding tetramers of stronger cooperativity. Moreover, the halogen-halogen interactions in bromoamine and iodoamine dimers with the geometries of the corresponding trimers and tetramers are examined using energy decomposition analysis methods (supermolecular, canonical EDA and SAPT) and the Kohn-Sham molecular orbital model. The results of the analysis indicate that although the interaction energy values for the dimers of the different spatial arrangement are very close to each other, their origin is substantially different. For pairs with the geometry of the trimers orbital interactions and electrostatic attraction are both weaker than for the corresponding dimers with the geometry of the tetramers. This is especially important because both donor-acceptor interactions and electrostatic attraction were previously proven to be responsible for cooperative effects occurring in the bromoamine and iodoamine tetramers.

Could Quantum Mechanical Properties Be Reflected on Classical Molecular Dynamics? The Case of Halogenated Organic Compounds of Biological Interest

Essential to understanding life, the biomolecular phenomena have been an important subject in science, therefore a necessary path to be covered to make progress in human knowledge. To fully comprehend these processes, the non-covalent interactions are the key. In this review, we discuss how specific protein-ligand interactions can be efficiently described by low computational cost methods, such as Molecular Mechanics (MM). We have taken as example the case of the halogen bonds (XB). Albeit generally weaker than the hydrogen bonds (HB), the XBs play a key role to drug design, enhancing the affinity and selectivity toward the biological target. Along with the attraction between two electronegative atoms in XBs explained by the σ-hole model, important orbital interactions, as well as relief of Pauli repulsion take place. Nonetheless, such electronic effects can be only well-described by accurate quantum chemical methods that have strong limitations dealing with supramolecular systems due to their high computational cost. To go beyond the poor description of XBs by MM methods, reparametrizing the force-fields equations can be a way to keep the balance between accuracy and computational cost. Thus, we have shown the steps to be considered when parametrizing force-fields to achieve reliable results of complex non-covalent interactions at MM level for In Silico drug design methods.

N-Oxide–N-oxide interactions and Cl...Cl halogen bonds in pentachloropyridine N-oxide: the many-body approach to interactions in the crystal state

Pentachloropyridine N-oxide, C5Cl5NO, crystallizes in the monoclinic space group P21/c. In the crystal structure, molecules are linked by C—Cl...Cl halogen bonds into infinite ribbons extending along the crystallographic [100] direction. These molecular aggregates are further stabilized by very short intermolecular N-oxide–N-oxide interactions into herringbone motifs. Computations based on quantum chemistry methods allowed for a more detailed description of the N-oxide–N-oxide interactions and Cl...Cl halogen bonds. For this purpose, Hirshfeld surface analysis and the many-body approach to interaction energy were applied.

Insights into the self-assembly steps of cyanuric acid toward rosette motifs: a DFT study

The nature of non-covalent interactions in self-assembling systems is a topic that has aroused great attention in literature. In this field, the 1,3,5-triazinane-2,4,6-trione or cyanuric acid (CA) is one of the most widely used molecules to formulate self-assembled materials or monolayers. In the present work, a variety of molecular aggregates of CA are examined using three different DFT functionals (B3LYP, B3LYP-D3, and ω-B97XD) in the framework of the quantum theory of atoms in molecules (QTAIM) and natural bond orbital (NBO) analysis. Herein, a step by step aggregation path is proposed and the origin of cooperative effects is also examined. It is shown that a greater cooperativity is not always associated with a greater binding energy, and the greatest cooperative effect occurs with highly directional hydrogen bonds. The intramolecular charge transfers play a key role in this effect. Graphical abstract The noncovalent interactions in cyanuric acid supramolecules were analyzed. The calculations provide insights into the self-assembly steps from dimers to rosette-like motif. The complexes with collinear hydrogen bonds show positive cooperativity, while in the arrangement with double hydrogen bonds the cooperative effect is essentially zero.

NX⋯Y halogen bonds. Comparison with NH⋯Y H-bonds and CX⋯Y halogen bonds

Quantum calculations examine how the NH⋯Y H-bond compares to the equivalent NX⋯Y halogen bond, as well as to comparable CH/CX donors.