Non-Classical Synthons: Supramolecular Recognition by S⋅⋅⋅O Chalcogen Bonding in Molecular Complexes of Riluzole

Classical examples of supramolecular recognition units or synthons are the ones formed by hydrogen bonds. Here, we report the ubiquity of a S⋅⋅⋅O chalcogen bonded synthon observed in a series of supramolecular complexes of the amyotrophic lateral sclerosis drug riluzole. Although the potential of higher chalcogens such as Se and Te to form robust and directional chalcogen bonded motifs is known, intermolecular sulfur chalcogen bonding is considered to be weak owing to the lower polarizability of S atoms. Here, the robustness and electronic nature of a S⋅⋅⋅O chalcogen bonding non-classical synthon, and the origin of its exceptional directionality have been explored. Bond orders of the drug–coformer chalcogen bonding are found to be as high as one third of a single bond, and they are largely ionic in nature. The contribution of the S⋅⋅⋅O chalcogen bonded motifs to the lattice energies of a series of crystals from the Cambridge Structural Database has been analyzed, showing they can be indeed significant, especially in molecules devoid of strong hydrogen bond donor groups.     

The Carbonyl⋅⋅⋅Tellurazole Chalcogen Bond as a Molecular Recognition Unit: From Model Studies to Supramolecular Organic Frameworks

In the last years, chalcogen bonding, the noncovalent interaction involving chalcogen centers, has emerged as interesting alternative to the ubiquitous hydrogen bonding in many research areas. Here, we could show by means of high‐level quantum chemical calculations that the carbonyl???tellurazole chalcogen bond is at least as strong as conventional hydrogen bonds. Using the carbonyl???tellurazole binding motif, we were able to design complex supramolecular networks in solid phase starting from tellurazole‐substituted cyclic peptides. X‐ray analyses reveal that the rigid structure of the cyclic peptides is caused by hydrogen bonds, whereas the supramolecular network is held together by chalcogen bonding. The type of the supramolecular network depends on peptide used; both linear wires and a honeycomb‐like supramolecular organic framework (SOF) were observed. The unique structure of the SOF shows two channels filled with different types of solvent mixtures that are either locked or freely movable.     

Supramolecular Wiring of Benzo‐1,3‐chalcogenazoles through Programmed Chalcogen Bonding Interactions

The high‐yielding synthesis of 2‐substituted benzo‐1,3‐tellurazoles and benzo‐1,3‐selenazoles through a dehydrative cyclization reaction has been reported, giving access to a large variety of benzo‐1,3‐chalcogenazoles. Exceptionally, these aromatic heterocycles proved to be very stable and thus very handy to form controlled solid‐state organizations in which wire‐like polymeric structures are formed through secondary N???Y bonding interactions (SBIs) engaging the chalcogen (Y=Se or Te) and nitrogen atoms. In particular, it has been shown that the recognition properties of the chalcogen centre at the solid state could be programmed by selectively barring one of its σ‐holes through a combination of electronic and steric effects exerted by the substituent at the 2‐position. As predicted by the electrostatic potential surfaces calculated by quantum chemical modelling, the pyridyl groups revealed to be the stronger chalcogen bonding acceptors, and thus the best ligand candidate for programming the molecular organization at the solid state. In contrast, the thiophenyl group is an unsuitable substituent for establishing SBIs in this molecular system as it gives rise to chalcogen–chalcogen repulsion. The weaker chalcogen donor properties of the Se analogues trigger the formation of feeble N???Se contacts, which are manifested in similar solid‐state polymers featuring longer nitrogen–chalcogen distances.     

Self‐Complementary Dimers of Oxalamide‐Functionalized Resorcinarene Tetrabenzoxazines

Self‐complementarity is a useful concept in supramolecular chemistry, molecular biology and polymeric systems. Two resorcinarene tetrabenzoxazines decorated with four oxalamide groups were synthesized and characterized. The oxalamide groups possessed self‐complementary hydrogen bonding sites between the carbonyls and amide groups. The self‐complementary nature of the oxalamide groups resulted in self‐included dimeric assemblies. The hydrogen bonding interactions within the tetrabenzoxazines gave rise to the formation of dimers, which were confirmed by single‐crystal X‐ray diffractions analysis and supported by NMR spectroscopy and mass spectrometry. The self‐included dimers were connected by numerous and strong intermolecular N?H???O and C?H???O hydrogen bonds supplemented with C?H???π interactions, forming one‐dimensional polymers, which were then further linked into three‐dimensional networks.     

Engineering Te-Containing Recognition Modules for Chalcogen Bonding: Towards Supramolecular Polymeric Materials

Aiming at the preparation of one-dimensional (1D) chalcogen-bonded supramolecular polymers at the solid state, this work describes the different syntheses which have been challenged to obtain ditopic molecular modules. At first, tellurazolopyridyl (TZP) rings have been chosen as recognition units, given their well-proven ability and persistency to self-assemble through double Te⋅⋅⋅N chalcogen bonds (ChBs). The second synthetic strategy dealt with the preparation of pyridyl-modified ebselen Te-containing analogues. By attempting several synthetic protocols, the targeted ebselen derivatives could not be obtained, whereas an unexpected Te-containing lactone as well as a spiro-type Te(IV)-containing derivatives were isolated, with the latter investigated by X-ray diffraction (XRD) analysis.     

Color Polymorphism: Understanding the Diverse Solid‐State Packing and Color in Dimethyl‐3,6‐dichloro‐2,5‐dihydroxyterephthalate

Dimethyl‐3,6‐dichloro‐2,5‐dihydroxyterephthalate (MCHT) is known to exist in three differently packed crystals having three different colors, namely yellow (Y), light yellow (LY), and white (W). Apart from the difference in their color, the molecules in the crystals also differ in their intramolecular O?H???O and O?H???Cl hydrogen bonds. Time‐dependent DFT calculations reveal the role of the various types of hydrogen bonds in controlling the color of the polymorphs. Mechanistic pathways that lead to such transformations in the crystal are elucidated by solid‐state dispersion‐corrected DFT studies. Relative stabilities of the various polymorphs rationalize the experimentally observed transformations between them. Calculations reveal that the minimum‐energy pathway for the conversion of the Y form to a W form is through stepwise disrotatory motion of the two ?OH groups through a hybrid intermediate having one intramolecular O?H???O and one O?H???Cl bond. The LY form is shown to exist on the higher‐energy pathway involving a concerted Y→W transformation.     

Micro‐ and Nanometer‐Scale Porous,Fibrous, and Sheet Architectures Constructed by Supramolecular Assemblies of Dipyrrolyldiketones

X‐ray analysis of some 1,3‐dipyrrolyl‐1,3‐propanediones synthesized from pyrroles and malonyl chloride derivatives revealed 1D supramolecular networks formed by N? H???O?C interactions in the solid state. Micro‐ and nanometer‐scale morphologies of porous, fibrous, and sheet structures were fabricated by hydrogen‐bonding interactions and determined by fine‐tuning the substituents and the solvents used. Of the unique polymorphs, ordered 2D lamellar sheet structures of the derivatives with long alkyl chains (C₁₆H₃₃, C₁₄H₂₉, and so on) were constructed by van der Waals hydrophobic effects between aliphatic chains as well as hydrogen bonding.     

Intermolecular Weak Interactions in HTeXH Dimers (X=O,S, Se,Te): Hydrogen Bonds,Chalcogen–Chalcogen Contacts and Chiral Discrimination

A theoretical study of the HTeXH (X=O, S, Se and Te) monomers and homodimers was carried out by means of second‐order Møller‐Plesset perturbation theory (MP2) computational methods. In the case of monomers, the isomerization energy from HTeXH to H₂Te=X and H₂X=Te (X=O, S, Se, and Te) and the rotational transition‐state barriers were obtained. Due to the chiral nature of these compounds, homo and heterochiral dimers were found. The electron density of the complexes was characterized with the atoms‐in‐molecules (AIM) methodology, finding a large variety of interactions. The charge transfer within the dimers was analyzed by means of natural bond orbitals (NBO). The density functional theory‐symmetry adapted perturbation theory (DFT‐SAPT) method was used to compute the components of the interaction energies. Hydrogen bonds and chalcogen–chalcogen interactions were characterized and their influence analyzed concerning the stability and chiral discrimination of the dimers.     

A Combined Experimental and Theoretical Electron Density Study of Intra‐ and Intermolecular Interactions in Thiourea S,S‐Dioxide

The thiourea S,S‐dioxide molecule is recognized as a zwitterion with a high dipole moment and an unusually long C? S bond. The molecule has a most interesting set of intermolecular interactions in the crystalline state—a relatively strong O???H? N hydrogen bond and very weak intermolecular C???S and N???O interactions. The molecule has C_s symmetry, and each oxygen atom is hydrogen‐bonded to two hydrogen atoms with O???H? N distances of 2.837 and 2.826 Å and angles of 176.61 and 158.38°. The electron density distribution is obtained both from Xray diffraction data at 110 K and from a periodic density functional theory (DFT) calculation. Bond characterization is made in terms of the analysis of topological properties. The covalent characters of the C? N, N? H, C? S, and S? O bonds are apparent, and the agreement on the topological properties between experiment and theory is adequate. The features of the Laplacian distributions, bond paths, and atomic domains are comparable. In a systematic approach, DFT calculations are performed based on a monomer, a dimer, a heptamer, and a crystal to see the effect on the electron density distribution due to the intermolecular interactions. The dipole moment of the molecule is enhanced in the solid state. The typical values of ρ_b and H_b of the hydrogen bonds and weak intermolecular C???S and N???O interactions are given. All the interactions are verified by the location of the bond critical point and its associated topological properties. The isovalue surface of Laplacian charge density and the detailed atomic graph around each atomic site reveal the shape of the valence‐shell charge concentration and provide a reasonable interpretation of the bonding of each atom.     

Synthetic and Structural Studies of 1,8‐Chalcogen Naphthalene Derivatives

Four novel 1,8‐disubstituted naphthalene derivatives 4 – 7 that contain chalcogen atoms occupying the peri positions have been prepared and fully characterised by using X‐ray crystallography, multinuclear NMR spectroscopy, IR spectroscopy and MS. Molecular distortion due to noncovalent substituent interactions was studied as a function of the bulk of the interacting chalcogen atoms and the size and nature of the alkyl group attached to them. X‐ray data for 4 – 7 was compared to the series of known 1,8‐bis(phenylchalcogeno)naphthalenes 1 – 3 , which were themselves prepared from novel synthetic routes. A general increase in the E???E′ distance was observed for molecules containing bulkier atoms at the peri positions. The decreased S???S distance from phenyl‐ 1 and ethyl‐ 4 analogues is ascribed to a weaker chalcogen lone pair–lone pair repulsion acting in the ethyl analogue due to the presence of two equatorial S(naphthyl) ring conformations. Two novel peri‐substituted naphthalene sulfoxides of 1 , Nap(O?SPh)(SPh) 8 and Nap(O?SPh)₂ 9 , which contain different valence states of sulfur, were prepared and fully characterised by using X‐ray crystallography and multinuclear NMR spectroscopy, IR spectroscopy and MS. Molecular structures were analysed by using naphthalene ring torsions, peri‐atom displacement, splay angle magnitude, S???S interactions, aromatic ring orientations and quasi‐linear O?S???S arrangements. The axial S(naphthyl) rings in 8 and 9 are unfavourable for S???S contacts due to stronger chalcogen lone pair–lone pair repulsion. Although quasi‐linear O?S???S alignments suggest attractive interaction is conceivable, analysis of the B3LYP wavefunctions affords no evidence for direct bonding interactions between the S atoms.     

Experimental and Computational Probes of the Nature of Halogen Bonding: Complexes of Bromine‐Containing Molecules with Bromide Anions

The nature of halogen bonding is examined via experimental and computational characterizations of a series of associates between electrophilic bromocarbons R? Br (R? Br=CBr₃F, CBr₃NO₂, CBr₃COCBr₃, CBr₃CONH₂, CBr₃CN, etc.) and bromide anions. The [R? Br, Br^?] complexes show intense absorption bands in the 200–350 nm range which follow the same Mulliken correlation as those observed for the charge‐transfer associates of bromide anions with common organic π‐acceptors. For a wide range of the associates, intermolecular R? Br???Br^? separations decrease and intramolecular C? Br bond lengths increase proportionally to the Br^?→R? Br charge transfer; and the energies of R? Br???Br^? bonds are correlated with the linear combination of orbital (charge‐transfer) and electrostatic interactions. On the whole, spectral, structural and thermodynamic characteristics of the [R? Br, Br^?] complexes indicate that besides electrostatics, the orbital (charge‐transfer) interactions play a vital role in the R? Br???Br^? halogen bonding. This indicates that in addition to controlling the geometries of supramolecular assemblies, halogen bonding leads to electronic coupling between interacting species, and thus affects reactivity of halogenated molecules, as well as conducting and magnetic properties of their solid‐state materials.     

Orthogonal Halogen‐Bonding‐Driven 3D Supramolecular Assembly of Right‐Handed Synthetic Helical Peptides

Peptide‐mediated self‐assembly is a prevalent method for creating highly ordered supramolecular architectures. Herein, we report the first example of orthogonal C?X???X?C/C?X???π halogen bonding and hydrogen bonding driven crystalline architectures based on synthetic helical peptides bearing hybrids of l ‐sulfono‐γ‐AApeptides and natural amino acids. The combination of halogen bonding, intra‐/intermolecular hydrogen bonding, and intermolecular hydrophobic interactions enabled novel 3D supramolecular assembly. The orthogonal halogen bonding in the supramolecular architecture exerts a novel mechanism for the self‐assembly of synthetic peptide foldamers and gives new insights into molecular recognition, supramolecular design, and rational design of biomimetic structures.     

A Crystallographic Charge Density Study of the Partial Covalent Nature of Strong N⋅⋅⋅Br Halogen Bonds

The covalent nature of strong N?Br???N halogen bonds in a cocrystal ( 2 ) of N‐bromosuccinimide ( NBS ) with 3,5‐dimethylpyridine ( lut ) was determined from X‐ray charge density studies and compared to a weak N?Br???O halogen bond in pure crystalline NBS ( 1 ) and a covalent bond in bis(3‐methylpyridine)bromonium cation (in its perchlorate salt ( 3 ). In 2 , the donor N?Br bond is elongated by 0.0954 Å, while the Br???acceptor distance of 2.3194(4) is 1.08 Å shorter than the sum of the van der Waals radii. A maximum electron density of 0.38 e Å^?3 along the Br???N halogen bond indicates a considerable covalent contribution to the total interaction. This value is intermediate to 0.067 e Å^?3 for the Br???O contact in 1 , and approximately 0.7 e Å^?3 in both N?Br bonds of the bromonium cation in 3 . A calculation of the natural bond order charges of the contact atoms, and the σ*(N1?Br) population of NBS as a function of distance between NBS and lut , have shown that charge transfer becomes significant at a Br???N distance below about 3 Å.     

The Importance of 1,5‐Oxygen⋅⋅⋅Chalcogen Interactions in Enantioselective Isochalcogenourea Catalysis

The importance of 1,5‐O???chalcogen (Ch) interactions in isochalcogenourea catalysis (Ch=O, S, Se) is investigated. Conformational analyses of N‐acyl isochalcogenouronium species and comparison with kinetic data demonstrate the significance of 1,5‐O???Ch interactions in enantioselective catalysis. Importantly, the selenium analogue demonstrates enhanced rate and selectivity profiles across a range of reaction processes including nitronate conjugate addition and formal [4+2] cycloadditions. A gram‐scale synthesis of the most active selenium analogue was developed using a previously unreported seleno‐Hugerschoff reaction, allowing the challenging kinetic resolutions of tertiary alcohols to be performed at 500 ppm catalyst loading. Density functional theory (DFT) and natural bond orbital (NBO) calculations support the role of orbital delocalization (occurring by intramolecular chalcogen bonding) in determining the conformation, equilibrium population, and reactivity of N‐acylated intermediates.     

Double Chalcogen Bonds: Crystal Engineering Stratagems via Diffraction and Multinuclear Solid-State Magnetic Resonance Spectroscopy

Group 16 chalcogens potentially provide Lewis-acidic σ-holes, which are able to form attractive supramolecular interactions with electron rich partners through chalcogen bonds. Here, a multifaceted experimental and computational study of a large series of novel chalcogen-bonded cocrystals, prepared using the principles of crystal engineering, is presented. Single-crystal X-ray diffraction studies reveal that dicyanoselenadiazole and dicyanotelluradiazole derivatives work as promising supramolecular synthons with the ability to form double chalcogen bonds with a wide range of electron donors including halides and oxygen- and nitrogen-containing heterocycles. Extensive ⁷⁷Se and ¹²⁵Te solid-state nuclear magnetic resonance spectroscopic investigations of cocrystals establish correlations between the NMR parameters of selenium and tellurium and the local chalcogen bonding geometry. The relationships between the electronic environment of the chalcogen bond and the ⁷⁷Se and ¹²⁵Te chemical shift tensors were elucidated through a natural localized molecular orbital density functional theory analysis. This systematic study of chalcogen-bond-based crystal engineering lays the foundations for the preparation of the various multicomponent systems and establishes solid-state NMR protocols to detect these interactions in powdered materials.     

Sterically Crowded peri‐Substituted Naphthalene Phosphines and their PV Derivatives

Three sterically crowded peri‐substituted naphthalene phosphines, Nap[PPh₂][ER] (Nap=naphthalene‐1,8‐diyl; ER=SEt, SPh, SePh) 1–3 , which contain phosphorus and chalcogen functional groups at the peri positions have been prepared. Each phosphine reacts to form a complete series of P^V chalcogenides Nap[P(E′)(Ph₂)(ER)] (E′=O, S, Se). The novel compounds were fully characterised by using X‐ray crystallography and multinuclear NMR spectroscopy, IR spectroscopy and MS. X‐ray data for 1 , 2 , n O , n S , n Se (n=1–3) are compared. Eleven molecular structures have been analysed by naphthalene ring torsions, peri‐atom displacement, splay angle magnitude, X???E interactions, aromatic ring orientations and quasi‐linear arrangements. An increase in the congestion of the peri region following the introduction of heavy chalcogen atoms is accompanied by a general increase in naphthalene distortion. P???E distances increase for molecules that contain bulkier atoms at the peri positions and also when larger chalcogen atoms are bound to phosphorus. The chalcogenides adopt similar conformations that contain a quasi‐linear E???P? C fragment, except for 3 O , which displays a twist‐axial‐twist conformation resulting in the formation of a linear O???Se? C alignment. Ab initio MO calculations performed on 2 O , 3 O , 3 S and 3 Se reveal Wiberg bond index values of 0.02 to 0.04, which indicates only minor non‐bonded interactions; however, calculations on radical cations of 3 O , 3 S and 3 Se reveal increased values (0.14–0.19).     

Cation–Cation Pairing by NCH⋅⋅⋅O Hydrogen Bonds

The pairing of ions of opposite charge is a fundamental principle in chemistry, and is widely applied in synthesis and catalysis. In contrast, cation–cation association remains an elusive concept, lacking in supporting experimental evidence. While studying the structure and properties of 4‐oxopiperidinium salts [OC₅H₈NH₂]X for a series of anions X^? of decreasing basicity, we observed a gradual self‐association of the cations, concluding in the formation of an isolated dicationic pair. In 4‐oxopiperidinium bis(trifluoromethylsulfonyl)amide, the cations are linked by N? H???O?C hydrogen bonds to form chains, flanked by hydrogen bonds to the anions. In the tetra(perfluoro‐tert‐butoxy)aluminate salt, the anions are fully separated from the cations, and the cations associate pairwise by N? C? H???O?C hydrogen bonds. The compounds represent the first genuine examples of self‐association of simple organic cations based merely on hydrogen bonding as evidenced by X‐ray structure analysis, and provide a paradigm for an extension of this class of compounds.     

The Double‐Faced Nature of Hydrogen Bonding in Hydroxy‐Functionalized Ionic Liquids Shown by Neutron Diffraction and Molecular Dynamics Simulations

We characterize the double‐faced nature of hydrogen bonding in hydroxy‐functionalized ionic liquids by means of neutron diffraction with isotopic substitution (NDIS), molecular dynamics (MD) simulations, and quantum chemical calculations. NDIS data are fit using the empirical potential structure refinement technique (EPSR) to elucidate the nearest neighbor H???O and O???O pair distribution functions for hydrogen bonds between ions of opposite charge and the same charge. Despite the presence of repulsive Coulomb forces, the cation–cation interaction is stronger than the cation–anion interaction. We compare the hydrogen‐bond geometries of both “doubly charged hydrogen bonds” with those reported for molecular liquids, such as water and alcohols. In combination, the NDIS measurements and MD simulations reveal the subtle balance between the two types of hydrogen bonds: The small transition enthalpy suggests that the elusive like‐charge attraction is almost competitive with conventional ion‐pair formation.     

α‐Halogenoacetanilides as Hydrogen‐Bonding Organocatalysts that Activate Carbonyl Bonds: Fluorine versus Chlorine and Bromine

α‐Halogenoacetanilides (X=F, Cl, Br) were examined as H‐bonding organocatalysts designed for the double activation of C?O bonds through NH and CH donor groups. Depending on the halide substituents, the double H‐bond involved a nonconventional C?H???O interaction with either a H?CX_n (n=1–2, X=Cl, Br) or a H?C_Ar bond (X=F), as shown in the solid‐state crystal structures and by molecular modeling. In addition, the catalytic properties of α‐halogenoacetanilides were evaluated in the ring‐opening polymerization of lactide, in the presence of a tertiary amine as cocatalyst. The α‐dichloro‐ and α‐dibromoacetanilides containing electron‐deficient aromatic groups afforded the most attractive double H‐bonding properties towards C?O bonds, with a N?H???O???H?CX₂ interaction.     

Theoretical Density Functional Theory insights into the nature of chalcogen bonding between CX2 (X = S,Se, Te) and diazine from monomer to supramolecular complexes

Chalcogen bonding is a noncovalent interaction, highly similar to halogen and hydrogen bonding, occurring between a chalcogen atom and a nucleophilic region. Two density functional theory (DFT) approaches B3LY-D3 and B97-D3 were performed on a series of complexes formed between CX₂ (X = S, Se, Te) and diazine (pyridazine, pyrimidine and pyrazine). Chalcogen atoms prefer interacting with the lone pair of a nitrogen atom rather than with the π-cloud of an aromatic ring. CTe₂ and CSe₂ form a stronger chalcogen bond than CS₂. The electrostatic potential of CX₂ (X = S, Se and Te) reveals the presence of two equivalent σ-holes, one on each chalcogen atom. These CX₂ molecules interact with diazine giving rise to supramolecular interactions. Wiberg bond index and second-order perturbation theory analysis in NBO were performed to better understand the nature of the chalcogen bond interaction.