Principles Guiding the Square Bonding Motif Containing a Pair of Chalcogen Bonds between Chalcogenadiazoles

Electrostatic Surface Potentials and Chalcogen-Bonding Motifs of Substituted 2,1,3-Benzoselenadiazoles Probed via 77Se Solid-State NMR Spectroscopy

Chalcogen bonds (ChB) are moderately strong, directional, and specific non-covalent interactions that have garnered substantial interest over the last decades. Specifically, the presence of two σ-holes offers great potential for crystal engineering, catalysis, biochemistry, and molecular sensing. However, ChB applications are currently hampered by a lack of methods to characterize and control chalcogen bonds. Here, we report on the influence of various substituents (halogens, cyano, and methyl groups) on the observed self-complementary ChB networks of 2,1,3-benzoselenadiazoles. From molecular electrostatic potential calculations, we show that the electrostatic surface potentials (ESP) of the σ-holes on selenium are largely influenced by the electron-withdrawing character of these substituents. Structural analyses via X-ray diffraction reveal a variety of ChB geometries and binding modes that are rationalized via the computed ESP maps, although the structure of 5,6-dimethyl-2,1,3-benzoselenadiazole also demonstrates the influence of steric interactions. 77Se solid-state magic-angle spinning NMR spectroscopy, in particular the analysis of the selenium chemical shift tensors, is found to be an effective probe able to characterize both structural and electrostatic features of these self-complementary ChB systems. We find a positive correlation between the value of the ESP maxima at the σ-holes and the experimentally measured 77Se isotropic chemical shift, while the skew of the chemical shift tensor is established as a metric which is reflective of the ChB binding motif.

On-Surface Molecular Recognition Driven by Chalcogen Bonding

Chalcogen bonding interactions (ChBIs) have been widely employed to create ordered noncovalent assemblies in solids and liquids. Yet, their ability to engineer molecular self-assembly on surfaces has not been demonstrated. Here, we report the first demonstration of on-surface molecular recognition solely governed by ChBIs. Scanning tunneling microscopy and ab initio calculations reveal that a pyrenyl derivative can undergo noncovalent chiral dimerization on the Au(111) surface through double Ch···N interactions involving Te- or Se-containing chalcogenazolo pyridine motifs. In contrast, reference chalcogenazole counterparts lacking the pyridyl moiety fail to form regular self-assemblies on Au, resulting in disordered assemblies.

Influence of donor point modifications on the assembly of chalcogen-bonded organic frameworks

Incremental, single-atom substitutions of Se-based chalcogen bond (Ch-bond) donors with stronger donating Te centers were implemented in two new triptycene tris(1,2,5-chalcogenadiazole) tectons. The appreciably more favorable Ch-bonding ability of the Te-based donors promotes assembly of low-density networks and more stable Ch-bonded organic frameworks (ChOFs).

Halide Complexes of 5,6-Dicyano-2,1,3-Benzoselenadiazole with 1 : 4 Stoichiometry: Cooperativity between Chalcogen and Hydrogen Bonding

The [M4 -Hal]- (M=the title compound; Hal=Cl, Br, and I) complexes were isolated in the form of salts of [Et4 N]+ cation and characterized by XRD, NMR, UV-Vis, DFT, QTAIM, EDD, and EDA. Their stoichiometry is caused by a cooperative interplay of σ-hole-driven chalcogen (ChB) and hydrogen (HB) bondings. In the crystal, [M4 -Hal]- are connected by the π-hole-driven ChB; overall, each [Hal]- is six-coordinated. In the ChB, the electrostatic interaction dominates over orbital and dispersion interactions. In UV-Vis spectra of the M+[Hal]- solutions, ChB-typical and [Hal]- -dependent charge-transfer bands are present; they reflect orbital interactions and allow identification of the individual [Hal]- . However, the structural situation in the solutions is not entirely clear. Particularly, the UV-Vis spectra of the solutions are different from the solid-state spectra of the [Et4 N]+ [M4 -Hal]- ; very tentatively, species in the solutions are assigned [M-Hal]- . It is supposed that the formation of the [M4 -Hal]- proceeds during the crystallization of the [Et4 N]+ [M4 -Hal]- . Overall, M can be considered as a chromogenic receptor and prototype sensor of [Hal]- . The findings are also useful for crystal engineering and supramolecular chemistry.

Competition Between the Two σ-Holes in the Formation of a Chalcogen Bond

A chalcogen atom Y contains two separate σ-holes when in a R₁ YR₂ molecular bonding pattern. Quantum chemical calculations consider competition between these two σ-holes to engage in a chalcogen bond (ChB) with a NH₃ base. R groups considered include F, Br, I, and tert-butyl (tBu). Also examined is the situation where the Y lies within a chalcogenazole ring, where its neighbors are C and N. Both electron-withdrawing substituents R₁ and R₂ act cooperatively to deepen the two σ-holes, but the deeper of the two holes consistently lies opposite to the more electron-withdrawing group, and is also favored to form a stronger ChB. The formation of two simultaneous ChBs in a triad requires the Y atom to act as double electron acceptor, and so anti-cooperativity weakens each bond relative to the simple dyad. This effect is such that some of the shallower σ-holes are unable to form a ChB at all when a base occupies the other site.

Exploring the Dynamical Nature of Intermolecular Hydrogen Bonds in Benzamide, Quinoline and Benzoic Acid Derivatives

The hydrogen bonds properties of 2,6-difluorobenzamide, 5-hydroxyquinoline and 4-hydroxybenzoic acid were investigated by Car-Parrinello and path integral molecular dynamics (CPMD and PIMD), respectively. The computations were carried out in vacuo and in the crystalline phase. The studied complexes possess diverse networks of intermolecular hydrogen bonds (N-H…O, O-H…N and O-H…O). The time evolution of hydrogen bridges gave a deeper insight into bonds dynamics, showing that bridged protons are mostly localized on the donor side; however, the proton transfer phenomenon was registered as well. The vibrational features associated with O-H and N-H stretching were analyzed on the basis of the Fourier transform of the atomic velocity autocorrelation function. The spectroscopic effects of hydrogen bond formation were studied. The PIMD revealed quantum effects influencing the hydrogen bridges providing more accurate free energy sampling. It was found that the N…O or O…O interatomic distances decreased (reducing the length of the hydrogen bridge), while the O-H or N-H covalent bond was elongated, which led to the increase in the proton sharing. Furthermore, Quantum Theory of Atoms in Molecules (QTAIM) was used to give insight into electronic structure parameters. Finally, Symmetry-Adapted Perturbation Theory (SAPT) was employed to estimate the energy contributions to the interaction energy of the selected dimers.

Chalcogen Bond as a Factor Stabilizing Ligand Conformation in the Binding Pocket of Carbonic Anhydrase IX Receptor Mimic

It is postulated that the overexpression of Carbonic Anhydrase isozyme IX in some cancers contributes to the acidification of the extracellular matrix. It was proved that this promotes the growth and metastasis of the tumor. These observations have made Carbonic Anhydrase IX an attractive drug target. In the light of the findings and importance of the glycoprotein in the cancer treatment, we have employed quantum-chemical approaches to study non-covalent interactions in the binding pocket. As a ligand, the acetazolamide (AZM) molecule was chosen, being known as a potential inhibitor exhibiting anticancer properties. First-Principles Molecular Dynamics was performed to study the chalcogen and other non-covalent interactions in the AZM ligand and its complexes with amino acids forming the binding site. Based on Density Functional Theory (DFT) and post-Hartree-Fock methods, the metric and electronic structure parameters were described. The Non-Covalent Interaction (NCI) index and Atoms in Molecules (AIM) methods were applied for qualitative/quantitative analyses of the non-covalent interactions. Finally, the AZM-binding pocket interaction energy decomposition was carried out. Chalcogen bonding in the AZM molecule is an important factor stabilizing the preferred conformation. Free energy mapping via metadynamics and Path Integral molecular dynamics confirmed the significance of the chalcogen bond in structuring the conformational flexibility of the systems. The developed models are useful in the design of new inhibitors with desired pharmacological properties.

Noncovalent interactions between benzochalcogenadiazoles and nitrogen bases

A theoretical study has been carried out on the intermolecular interactions between tetrafluoro-benzochalcogenadiazoles (chalcogen = S, Se, Te) and a series of nitrogen bases (FCN, ClCN, NP, trans-N₂H₂, pyridine, pyrazole, imidazole) at the B97-D3/def2-TZVP level, to obtain a better insight into the nature and strength of Ch···N chalcogen bond and secondary interaction in the binary and 1:2 ternary complexes. The dispersion force plays a prominent role on the stability of the sulfur complexes, and the electrostatic effect enhanced for the heavier chalcogen complexes. Most of intermolecular bonds display the characters of closed-shell and noncovalent interaction. For the complexes involving pyridine and imidazole, chalcogen bond is stronger than hydrogen bond, while the strength of chalcogen bond is equivalent to the secondary interaction for other complexes. With the addition of nitrogen base in the 1:2 complexes, chalcogen bond is weakened, while the secondary interaction remains unchanged. In the 1:2 complexes formed by pyridine and imidazole, stronger chalcogen bond results in larger negative cooperativity than that of other complexes.

Chalcogen-bonded donor–acceptor complexes of 5,6-dicyano[1,2,5]selenadiazolo[3,4-b]pyrazine with halide ions

With halides X− (X = Cl, Br, I) 5,6-dicyano-[1,2,5]selenadiazolo[3,4-b]pyrazine 1 forms chalcogen-bonded complexes [1–X]− structurally defined by XRD. UV/Vis spectra of [1–X]− feature red-shifted charge-transfer bands in the Vis part.