Comparison of Bifurcated Halogen with Hydrogen Bonds

Crystal Structure Survey and Theoretical Analysis of Bifurcated Halogen Bonds

The possibility that two Lewis bases can share a single halogen atom within the context of a bifurcated halogen bond (XB) is explored first by a detailed examination of the CSD. Of the more than 22,000 geometries that fit the definition of an XB (with X = Cl, Br, I), less than 2% are bifurcated. There is a heavy weighting of I in such bifurcated arrangements as opposed to Br, which prefers monofurcated bonds. The conversion from mono to bifurcated is associated with a smaller number of short contact distances, as well as a trend toward lesser linearity. The two XBs within a bifurcated system are somewhat symmetrical: the two lengths generally differ by less than 0.05 Å, and the two XB angles are within several degrees of one another. Quantum calculations of model systems reflect the patterns observed in crystals and reinforce the idea that the negative cooperativity within a bifurcated XB weakens and lengthens each individual bond.

Bifurcated Halogen Bond-Driven Supramolecular Double Helices from 1,2-Dihalotetrafluorobenzene and 2,2′-Bi(1,8-naphthyridine)

The unique enantiomeric pairs of double helices have been found in the structure of the cocrystal between 1,2-diiodotetrafluorobenzene and 2,2′-bi(1,8-naphthyridine). The formation of the supramolecular double helices is driven by the strong bifurcated iodine bonds which can force the herringbone packing arrangement of the molecules 2,2′-bi(1,8-naphthyridine) into a face-to-face π···π stacking pattern. In contrast, the cocrystal between 1,2-dibromotetrafluorobenzene (or 1,2-dichlorotetrafluorobenzene) and 2,2′-bi(1,8-naphthyridine) was not obtained under the same conditions. The interaction energies of the bifurcated halogen bonds and π···π stacking interactions were computed with the reliable dispersion-corrected density functional theory. The computational results show that the bifurcated iodine bond is much stronger than the bifurcated bromine bond and bifurcated chlorine bond, and it is the much stronger bifurcated iodine bond that makes the cocrystal of 1,2-diiodotetrafluorobenzene and 2,2′-bi(1,8-naphthyridine) much easier to be synthesized.

Anticooperativity of Multiple Halogen Bonds and Its Effect on Stoichiometry of Cocrystals of Perfluorinated Iodobenzenes

To investigate influences on the topicity of perfluorinated halobenzenes as halogen bond (XB) donors in the solid state, we have conducted a database survey and prepared 18 novel cocrystals of potentially ditopic (13ditfb, 14ditfb) and tritopic (135titfb) XB donors with 15 monotopic pyridines. 135titfb shows high tendency to be mono- or ditopic, but with strong bases it can act as a tritopic XB donor. DFT calculations have shown that binding of a single acceptor molecule on one of the iodine atoms of the XB donor reduces the ESP_max on the remaining iodine atoms and dramatically decreases their potential for forming further halogen bonds, which explains both the high occurrence of crystal structures where the donors do not achieve their maximal topicity and the observed differences in halogen bond lengths. Despite the fact that this effect increases with the basicity of the acceptor, when the increase of halogen bond energy due to the basicity of the acceptor compensates its decrease due to the reduction of the acidity of the donor, it enables strong bases to form cocrystals in which a potentially polytopic XB donor achieves its maximal topicity.

Ability of Lewis Acids with Shallow σ-Holes to Engage in Chalcogen Bonds in Different Environments

Molecules of the type XYT = Ch (T = C, Si, Ge; Ch = S, Se; X,Y = H, CH₃, Cl, Br, I) contain a σ-hole along the T = Ch bond extension. This hole can engage with the N lone pair of NCH and NCCH₃ so as to form a chalcogen bond. In the case of T = C, these bonds are rather weak, less than 3 kcal/mol, and are slightly weakened in acetone or water. They owe their stability to attractive electrostatic energy, supplemented by dispersion, and a much smaller polarization term. Immersion in solvent reverses the electrostatic interaction to repulsive, while amplifying the polarization energy. The σ-holes are smaller for T = Si and Ge, even negative in many cases. These Lewis acids can nonetheless engage in a weak chalcogen bond. This bond owes its stability to dispersion in the gas phase, but it is polarization that dominates in solution.

Assessment of Computational Tools for Predicting Supramolecular Synthons

The ability to predict the most likely supramolecular synthons in a crystalline solid is a valuable starting point for subsequently predicting the full crystal structure of a molecule with multiple competing molecular recognition sites. Energy and informatics-based prediction models based on molecular electrostatic potentials (MEPs), hydrogen-bond energies (HBE), hydrogen-bond propensity (HBP), and hydrogen-bond coordination (HBC) were applied to the crystal structures of twelve pyrazole-based molecules. HBE, the most successful method, correctly predicted 100% of the experimentally observed primary intermolecular-interactions, followed by HBP (87.5%), and HBC = MEPs (62.5%). A further HBC analysis suggested a risk of synthon crossover and synthon polymorphism in molecules with multiple binding sites. These easy-to-use models (based on just 2-D chemical structure) can offer a valuable risk assessment of potential formulation challenges.

Hydrogen Bonding in Natural and Unnatural Base Pairs-A Local Vibrational Mode Study

In this work hydrogen bonding in a diverse set of 36 unnatural and the three natural Watson Crick base pairs adenine (A)-thymine (T), adenine (A)-uracil (U) and guanine (G)-cytosine (C) was assessed utilizing local vibrational force constants derived from the local mode analysis, originally introduced by Konkoli and Cremer as a unique bond strength measure based on vibrational spectroscopy. The local mode analysis was complemented by the topological analysis of the electronic density and the natural bond orbital analysis. The most interesting findings of our study are that (i) hydrogen bonding in Watson Crick base pairs is not exceptionally strong and (ii) the N-H⋯N is the most favorable hydrogen bond in both unnatural and natural base pairs while O-H⋯N/O bonds are the less favorable in unnatural base pairs and not found at all in natural base pairs. In addition, the important role of non-classical C-H⋯N/O bonds for the stabilization of base pairs was revealed, especially the role of C-H⋯O bonds in Watson Crick base pairs. Hydrogen bonding in Watson Crick base pairs modeled in the DNA via a QM/MM approach showed that the DNA environment increases the strength of the central N-H⋯N bond and the C-H⋯O bonds, and at the same time decreases the strength of the N-H⋯O bond. However, the general trends observed in the gas phase calculations remain unchanged. The new methodology presented and tested in this work provides the bioengineering community with an efficient design tool to assess and predict the type and strength of hydrogen bonding in artificial base pairs.

Anion-Anion Interactions in Aerogen-Bonded Complexes. Influence of Solvent Environment

Ab initio calculations are applied to the question as to whether a AeX₅^- anion (Ae = Kr, Xe) can engage in a stable complex with another anion: F^-, Cl^-, or CN^-. The latter approaches the central Ae atom from above the molecular plane, along its C₅ axis. While the electrostatic repulsion between the two anions prevents their association in the gas phase, immersion of the system in a polar medium allows dimerization to proceed. The aerogen bond is a weak one, with binding energies less than 2 kcal/mol, even in highly polar aqueous solvent. The complexes are metastable in the less polar solvents THF and DMF, with dissociation opposed by a small energy barrier.