Theoretical Investigations on Heteronuclear Chalcogen−Chalcogen Interactions: On the Nature of Weak Bonds between Chalcogen Centers

Isothiourea-Catalysed Acylative Dynamic Kinetic Resolution of Tetra-substituted Morpholinone and Benzoxazinone Lactols

The development of methods to allow the selective acylative dynamic kinetic resolution (DKR) of tetra-substituted lactols is a recognised synthetic challenge. In this manuscript, a highly enantioselective isothiourea-catalysed acylative DKR of tetra-substituted morpholinone and benzoxazinone-derived lactols is reported. The scope and limitations of this methodology have been developed, with high enantioselectivity and good to excellent yields (up to 89 %, 99 : 1 er) observed across a broad range of substrate derivatives incorporating substitution at N(4) and C(2), di- and spirocyclic substitution at C(5) and C(6), as well as benzannulation (>35 examples in total). The DKR process is amenable to scale-up on a 1 g laboratory scale. The factors leading to high selectivity in this DKR process have been probed through computation, with an N-C=O⋅⋅⋅isothiouronium interaction identified as key to producing ester products in highly enantioenriched form.

Synthesis of Tetra-Substituted 3-Hydroxyphthalide Esters by Isothiourea-Catalysed Acylative Dynamic Kinetic Resolution

A general and highly enantioselective method for the preparation of tetra-substituted 3-hydroxyphthalide esters via isothiourea-catalysed acylative dynamic kinetic resolution (DKR) is reported. Using (2S,3R)-HyperBTM (5 mol %) as the catalyst, the scope and limitations of this methodology have been extensively probed, with high enantioselectivity and good to excellent yields observed (>40 examples, up to 99 %, 99 : 1 er). Substitution of the aromatic core within the 3-hydroxyphthalide skeleton, as well as aliphatic and aromatic substitution at C(3), is readily tolerated. A diverse range of anhydrides, including those from bioactive and pharmaceutically relevant acids, can also be used. The high enantioselectivity observed in this DKR process has been probed computationally, with a key substrate heteroatom donor O⋅⋅⋅acyl-isothiouronium interaction identified through DFT analysis as necessary for enantiodiscrimination.

Isothiourea-Catalysed Acylative Kinetic Resolution of Tertiary Pyrazolone Alcohols

The development of methods for the selective acylative kinetic resolution (KR) of tertiary alcohols is a recognised synthetic challenge with relatively few successful substrate classes reported to date. In this manuscript, a highly enantioselective isothiourea-catalysed acylative KR of tertiary pyrazolone alcohols is reported. The scope and limitations of this methodology have been developed, with high selectivity observed across a broad range of substrate derivatives incorporating varying substitution at N(2)-, C(4)- and C(5)-, as well as bicyclic constraints within the pyrazolone scaffold (30 examples, selectivity factors (s) typically >100) at generally low catalyst loadings (1 mol %). The application of this KR method to tertiary alcohols derived directly from a natural product (geraniol), alongside pharmaceutically relevant drug compounds (indomethacin, gemfibrozil and probenecid), with high efficiency (s >100) is also described. The KR process is readily amenable to scale up using bench grade solvents and reagents, with effective resolution on a 50 g (0.22 mol) scale demonstrated. The key structural motif leading to excellent selectivity in this KR process has been probed through computation, with an NC=O⋅⋅⋅isothiouronium interaction from substrate to acylated catalyst observed within the favoured transition state. Similarly, the effect of C(5)-aryl substitution that leads to reduced experimental selectivity is probed, with a competitive π-isothiouronium interaction identified as leading to reduced selectivity.

One-Pot Access to Functionalised Malamides via Organocatalytic Enantioselective Formation of Spirocyclic β-Lactone-Oxindoles and Double Ring-Opening

Malamides (diamide derivatives of malic acid) are prevalent in nature and of significant biological interest, yet only limited synthetic methods to access functionalised enantiopure derivatives have been established to date. Herein, an effective synthetic method to generate this molecular class is developed through in situ formation of spirocyclic β-lactone-oxindoles (employing a known enantioselective isothiourea-catalysed formal [2+2] cycloaddition of C(1)-ammonium enolates and isatin derivatives) followed by a subsequent dual ring-opening protocol (of the β-lactone and oxindole) with amine nucleophiles. The application of this protocol is demonstrated across twelve examples to give densely functionalised malamide derivatives with high enantio- and diastereo-selectivity (up to >95:5 dr and >99:1 er).

De-epimerizing DyKAT of β-lactones generated by isothiourea-catalysed enantioselective [2 + 2] cycloaddition

An enantioselective isothiourea-catalysed [2 + 2] cycloaddition of C(1)-ammonium enolates with pyrazol-4,5-diones is used to construct spirocyclic β-lactones in good yields, excellent enantioselectivity (99 : 1 er) but with modest diastereocontrol (typically 70 : 30 dr). Upon ring-opening with morpholine or alternative nucleophilic amines and alcohols β-hydroxyamide and β-hydroxyester products are generated with enhanced diastereocontrol (up to >95 : 5 dr). Control experiments show that stereoconvergence is observed in the ring-opening of diastereoisomeric β-lactones, leading to a single product (>95 : 5 dr, >99 : 1 er). Mechanistic studies and DFT analysis indicate a substrate controlled Dynamic Kinetic Asymmetric Transformation (DyKAT) involving epimerisation at C(3) of the β-lactone under the reaction conditions, coupled with a hydrogen bond-assisted nucleophilic addition to the Si-face of the β-lactone and stereodetermining ring-opening. The scope and limitations of a one-pot protocol consisting of isothiourea-catalysed enantio-determining [2 + 2] cycloaddition followed by diastereo-determining ring-opening are subsequently developed. Variation within the anhydride ammonium enolate precursor, as well as N(1) and C(3) within the pyrazol-4,5-dione scaffold is demonstrated, giving a range of functionalised β-hydroxyamides with high diastereo- and enantiocontrol (>20 examples, up to >95 : 5 dr and >99 : 1 er) via this DyKAT.

Elucidating the Binding Mode of Sulfur- and Selenium-Based Cationic Chalcogen-Bond Donors

For a comparison of the interaction modes of various chalcogen-bond donors, 2-chalcogeno-imidazolium salts have been designed, synthesized, and studied by single crystal X-ray diffraction, solution NMR and DFT as well as for their ability to act as activators in an SN1-type substitution reaction. Their interaction modes in solution were elucidated based on NMR diffusion and chemical shift perturbation experiments, which were supported by DFT-calculations. Our finding is that going from lighter to the heavier chalcogens, hydrogen bonding plays a less, while chalcogen bonding an increasingly important role for the coordination of anions. Anion-π interactions also show importance, especially for the sulfur and selenium derivatives.

Unveiling the Nature and Strength of Selenium-Centered Chalcogen Bonds in Binary Complexes of SeO2 with Oxygen-/Sulfur-Containing Lewis Bases: Insights from Theoretical Calculations

Among various non-covalent interactions, selenium-centered chalcogen bonds (SeChBs) have garnered considerable attention in recent years as a result of their important contributions to crystal engineering, organocatalysis, molecular recognition, materials science, and biological systems. Herein, we systematically investigated π-hole-type Se∙∙∙O/S ChBs in the binary complexes of SeO2 with a series of O-/S-containing Lewis bases by means of high-level ab initio computations. The results demonstrate that there exists an attractive interaction between the Se atom of SeO2 and the O/S atom of Lewis bases. The interaction energies computed at the MP2/aug-cc-pVTZ level range from -4.68 kcal/mol to -10.83 kcal/mol for the Se∙∙∙O chalcogen-bonded complexes and vary between -3.53 kcal/mol and -13.77 kcal/mol for the Se∙∙∙S chalcogen-bonded complexes. The Se∙∙∙O/S ChBs exhibit a relatively short binding distance in comparison to the sum of the van der Waals radii of two chalcogen atoms. The Se∙∙∙O/S ChBs in all of the studied complexes show significant strength and a closed-shell nature, with a partially covalent character in most cases. Furthermore, the strength of these Se∙∙∙O/S ChBs generally surpasses that of the C/O-H∙∙∙O hydrogen bonds within the same complex. It should be noted that additional C/O-H∙∙∙O interactions have a large effect on the geometric structures and strength of Se∙∙∙O/S ChBs. Two subunits are connected together mainly via the orbital interaction between the lone pair of O/S atoms in the Lewis bases and the BD*(OSe) anti-bonding orbital of SeO2, except for the SeO2∙∙∙HCSOH complex. The electrostatic component emerges as the largest attractive contributor for stabilizing the examined complexes, with significant contributions from induction and dispersion components as well.

Enantioselective Synthesis in Continuous Flow: Polymer-Supported Isothiourea-Catalyzed Enantioselective Michael Addition-Cyclization with α-Azol-2-ylacetophenones

A packed reactor bed incorporating a polymer-supported isothiourea HyperBTM catalyst derivative has been used to promote the enantioselective synthesis of a range of heterocyclic products derived from α-azol-2-ylacetophenones and -acetamides combined with alkyl, aryl, and heterocyclic α,β-unsaturated homoanhydrides in continuous flow via an α,β-unsaturated acyl-ammonium intermediate. The products are generated in good to excellent yields and generally in excellent enantiopurity (up to 97:3 er). Scale-up is demonstrated on a 15 mmol scale, giving the heterocyclic product in 68% overall yield with 98:2 er after recrystallization.

Understanding divergent substrate stereoselectivity in the isothiourea-catalysed conjugate addition of cyclic α-substituted β-ketoesters to α,β-unsaturated aryl esters

The development of enantioselective synthetic methods capable of generating vicinal stereogenic centres, where one is tetrasubstituted (such as either an all-carbon quaternary centre or where one or more substituents are heteroatoms), is a recognised synthetic challenge. Herein, the enantioselective conjugate addition of a range of carbo- and heterocyclic α-substituted β-ketoesters to α,β-unsaturated aryl esters using the isothiourea HyperBTM as a Lewis base catalyst is demonstrated. Notably, divergent diastereoselectivity is observed through the use of either cyclopentanone-derived or indanone-derived substituted β-ketoesters with both generating the desired stereodefined products with high selectivity (>95 : 5 dr, up to 99 : 1 er). The scope and limitations of these processes are demonstrated, alongside application on gram scale. The origin of the divergent substrate selectivity has been probed through the use of DFT-analysis, with preferential orientation driven by dual stabilising CH⋯O interactions. The importance of solvation with strongly polar transition-states is highlighted and the SMD solvation model is demonstrated to capture solvation effects reliably.

Computational Insight into the Nature and Strength of the π-Hole Type Chalcogen∙∙∙Chalcogen Interactions in the XO2∙∙∙CH3YCH3 Complexes (X = S, Se, Te; Y = O, S, Se, Te)

In recent years, the non-covalent interactions between chalcogen centers have aroused substantial research interest because of their potential applications in organocatalysis, materials science, drug design, biological systems, crystal engineering, and molecular recognition. However, studies on π-hole-type chalcogen∙∙∙chalcogen interactions are scarcely reported in the literature. Herein, the π-hole-type intermolecular chalcogen∙∙∙chalcogen interactions in the model complexes formed between XO2 (X = S, Se, Te) and CH3YCH3 (Y = O, S, Se, Te) were systematically studied by using quantum chemical computations. The model complexes are stabilized via one primary X∙∙∙Y chalcogen bond (ChB) and the secondary C-H∙∙∙O hydrogen bonds. The binding energies of the studied complexes are in the range of -21.6~-60.4 kJ/mol. The X∙∙∙Y distances are significantly smaller than the sum of the van der Waals radii of the corresponding two atoms. The X∙∙∙Y ChBs in all the studied complexes except for the SO2∙∙∙CH3OCH3 complex are strong in strength and display a partial covalent character revealed by conducting the quantum theory of atoms in molecules (QTAIM), a non-covalent interaction plot (NCIplot), and natural bond orbital (NBO) analyses. The symmetry-adapted perturbation theory (SAPT) analysis discloses that the X∙∙∙Y ChBs are primarily dominated by the electrostatic component.

Chalcogen Bonding Catalysis: Tellurium, the Last Frontier?

Chalcogen bonding (ChB) is the non-covalent interaction occurring between chalcogen atoms as Lewis acid sites and atoms or groups of atoms able to behave as Lewis bases through their lone pair or π electrons. Analogously to its sister halogen bonding, the high directionality of this interaction was implemented for precise structural organizations in the solid state and in solution. Regarding catalysis, ChB is now accepted as a new mode of activation as demonstrated by the increased number of examples in the last five years. In the family of ChB catalysts, those based on tellurium rapidly appeared to overcome their lighter sulfur and selenium counterparts. In this review, we highlight the Lewis acid properties of tellurium-based derivatives in solution and summarize the start-of-the-art of their applications in catalysis.

Organically Templated Uranyl Sulfates and Selenates: Structural Complexity and Crystal Chemical Restrictions for Isotypic Compounds Formation

This paper reviews the state of the art in the structural chemistry of organically templated uranyl sulfates and selenates, which are considered as the most representative groups of U-bearing synthetic compounds. In total, there are 194 compounds known for both groups, the crystal structures of which include 84 various organic molecules. Structural studies and topological analysis clearly indicate complex crystal chemical limitations in terms of the isomorphic substitution implementation, since the existence of isotypic phases has to date been confirmed only for 24 compounds out of 194, which is slightly above 12%. The structural architecture of the entire compound depends on the combination of the organic and oxyanion parts, changes in which are sometimes realized even while maintaining the topology of the U-bearing complex. An increase in the size of the hydrocarbon part and number of charge functional groups of the organic cation leads to the formation of rare and more complex topologies. In addition, the crystal structures of two novel uranyl sulfates and one uranyl selenate, templated by isopropylammonium cations, are reported.

Enantioselective isothiourea-catalysed reversible Michael addition of aryl esters to 2-benzylidene malononitriles

Catalytic enantioselective transformations usually rely upon optimal enantioselectivity being observed in kinetically controlled reaction processes, with energy differences between diastereoisomeric transition state energies translating to stereoisomeric product ratios. Herein, stereoselectivity resulting from an unusual reversible Michael addition of an aryl ester to 2-benzylidene malononitrile electrophiles using an isothiourea as a Lewis base catalyst is demonstrated. Notably, the basicity of the aryloxide component and reactivity of the isothiourea Lewis base both affect the observed product selectivity, with control studies and crossover experiments indicating the feasibility of a constructive reversible Michael addition from the desired product. When this reversible addition is coupled with a crystallisation-induced diastereomer transformation (CIDT) it allows isolation of products in high yield and stereocontrol (14 examples, up to 95 : 5 dr and 99 : 1 er). Application of this process to gram scale, plus derivatisations to provide further useful products, is demonstrated.

Probing Regio- and Enantioselectivity in the Formal [2 + 2] Cycloaddition of C(1)-Alkyl Ammonium Enolates with β- and α,β-Substituted Trifluoromethylenones

The isothiourea-catalyzed regio- and enantioselective formal [2 + 2] cycloaddition of C(1)-alkyl and C(1)-unsubstituted ammonium enolates with β- and α,β-substituted trifluoromethylenones has been developed. In all cases, preferential [2 + 2]-cycloaddition over the alternative [4 + 2]-cycloaddition is observed, giving β-lactones with excellent diastereo- and enantioselectivity (34 examples, up to >95:5 dr, >99:1 er). The regioselectivity of the process was dictated by the nature of the substituents on both reaction components. Solely [2 + 2] cycloaddition products are observed when using α,β-substituted trifluoromethylenones or α-trialkylsilyl acetic acid derivatives; both [2 + 2] and [4 + 2] cycloaddition products are observed when using β-substituted trifluoromethylenones and α-alkyl-α-trialkylsilyl acetic acids as reactants, with the [2 + 2] cycloaddition as the major reaction product. The beneficial role of the α-silyl substituent within the acid component in this protocol has been demonstrated by control experiments.

Optoelectronic design and charge transport properties of Benzodifuran (BDF) isomers for organic electronic devices: DFT/TD-DFT insights

The primary goal of this work is to provide a comprehensive analysis of the charge transport and optoelectronic characteristics of all the isomers of benzodifuran (BDF) for organic electronic devices in order to suggest qualified materials/candidates for organic photovoltaic devices. Density functional theory (DFT) calculations were performed for all possible isomers of BDF and results are compared with corresponding experimental known isomers. Time Dependent-Density Functional Theory (TD-DFT) is used for the calculation of the absorption and HOMO-LUMO energy levels. To characterize the electronic charge transport state in these isomers, the ionization potentials (IP), reorganization energies (hole and electron), and electron affinities (EA) of all the isomers are investigated. Comparatively, all the BDF isomers are having low electron and hole reorganization energies and hence they can be used in the organic electronic material fabrication.

A Desilylative Approach to Alkyl Substituted C(1)-Ammonium Enolates: Application in Enantioselective [2+2] Cycloadditions

The catalytic generation of C(1)-ammonium enolates from the corresponding α-silyl-α-alkyl substituted carboxylic acids using the isothiourea HyperBTM is reported. This desilylative approach grants access to α-unsubstituted and α-alkyl substituted C(1)-ammonium enolates, which are typically difficult to access through traditional methods reliant upon deprotonation. The scope and limitations of this process is established in enantioselective [2+2]-cycloaddition processes with perfluoroalkylketones (31 examples, up to 96 % yield and >99 : 1 er), as well as selective [2+2]-cycloaddition with trifluoromethyl enones (4 examples, up to 75 % yield and >99 : 1 er). Preliminary mechanistic studies indicate this process proceeds through an initial kinetic resolution of an in situ prepared (±)-α-silyl-α-alkyl substituted anhydride, while the reaction process exhibits overall pseudo zero-order kinetics.

Catalytic Activation of Imines by Chalcogen Bond Donors in a Povarov [4+2] Cycloaddition Reaction

Recently, chalcogen bonding has been investigated in more detail in organocatalysis and the scope of activated functionalities continues to increase. Herein, the activation of imines in a Povarov [4+2] cycloaddition reaction with bidentate cationic chalcogen bond donors is presented. Tellurium-based Lewis acids show superior properties compared to selenium-based catalysts and inactive sulfur-based analogues. The catalytic activity of the chalcogen bonding donors increases with weaker binding anions. Triflate, however, is not suitable due to its participation in the catalytic pathway. A solvent screening revealed a more efficient activation in less polar solvents and a pronounced effect of solvent (and catalyst) on endo : exo diastereomeric ratio. Finally, new chiral chalcogen bonding catalysts were applied but provided only racemic mixtures of the product.

Isothiourea-Catalyzed [2 + 2] Cycloaddition of C(1)-Ammonium Enolates and N-Alkyl Isatins

Enantioselective [2 + 2] cycloaddition of C(1)-ammonium enolates generated catalytically using the isothiourea HyperBTM with N-alkyl isatins gives spirocyclic β-lactones. In situ ring opening with an amine nucleophile generates isolable highly enantioenriched products in up to 92:8 dr and in >99:1 er.

Various Sorts of Chalcogen Bonds Formed by an Aromatic System

The chalcogen Y atom in the aromatic ring of thiophene and its derivatives YC₄H₄ (Y = S, Se, Te) can engage in a number of different interactions with another such unit within the homodimer. Quantum calculations show that the two rings can be oriented perpendicular to one another in a T-shaped dimer in which the Y atom accepts electron density from the π-system of the other unit in a Y···π chalcogen bond (ChB). This geometry best takes advantage of attractions between the electrostatic potentials surrounding the two monomers. There are two other geometries in which the two Y atoms engage in a ChB with one another. However, instead of a simple interaction between a σ-hole on one Y and the lone pair of its neighbor, the interaction is better described as a pair of symmetrically equivalent Y···Y interactions, in which charge is transferred in both directions simultaneously, thereby effectively doubling the strength of the bond. These geometries differ from what might be expected based simply on the juxtaposition of the electrostatic potentials of the two monomers.

Cooperative Palladium/Isothiourea Catalyzed Enantioselective Formal (3+2) Cycloaddition of Vinylcyclopropanes and α,β-Unsaturated Esters

A protocol for the enantioselective synthesis of substituted vinylcyclopentanes has been realised using cooperative palladium and isothiourea catalysis. Treatment of vinylcyclopropanes with Pd(PPh3 )4 generates a zwitterionic π-allyl palladium intermediate that intercepts a catalytically generated α,β-unsaturated acyl ammonium species prepared from the corresponding α,β-unsaturated para-nitrophenyl ester and the isothiourea (R)-BTM. Intermolecular formal (3+2) cycloaddition between these reactive intermediates generates functionalised cyclopentanes in generally good yields and excellent diastereo- and enantiocontrol (up to >95 : 5 dr, 97 : 3 er), with the use of LiCl as an additive proving essential for optimal stereocontrol. To the best of our knowledge a dual transition metal/organocatalytic process involving α,β-unsaturated acyl ammonium intermediates has not been demonstrated previously.

Isothiourea-Catalyzed Enantioselective Michael Addition of Malonates to α,β-Unsaturated Aryl Esters

An enantioselective Michael addition of malonates to α,β-unsaturated para-nitrophenyl esters was achieved using the Lewis basic isothiourea HyperBTM, giving excellent levels of product enantioselectivity (up to >99:1 enantiomeric ratio) in good yields and with complete regioselectivity (>20:1 regioselectivity ratio) in the presence of alternative (phenyl ketone and ethyl ester) Michael acceptors. Density functional theory calculations indicate that N-acylation is rate-limiting. This constitutes a rare example of a highly enantioselective addition of simple, readily available malonates to α,β-unsaturated esters.

Simulating chalcogen bonding using molecular mechanics: a pseudoatom approach to model ebselen

The organoselenium compound ebselen has recently been investigated as a treatment for COVID-19; however, efforts to model ebselen in silico have been hampered by the lack of an efficient and accurate method to assess its binding to biological macromolecules. We present here a Generalized Amber Force Field modification which incorporates classical parameters for the selenium atom in ebselen, as well as a positively charged pseudoatom to simulate the σ-hole, a quantum mechanical phenomenon that dominates the chemistry of ebselen. Our approach is justified using an energy decomposition analysis of a number of density functional theory–optimized structures, which shows that the σ-hole interaction is primarily electrostatic in origin. Finally, our model is verified by conducting molecular dynamics simulations on a number of simple complexes, as well as the clinically relevant enzyme SOD1 (superoxide dismutase), which is known to bind to ebselen.

Graphical Abstract

Inter-anion chalcogen bonds: Are they anti-electrostatic in nature?

Inter-anion hydrogen and halogen bonds have emerged as counterintuitive linkers and inspired us to expand the range of this unconventional bonding pattern. Here, the inter-anion chalcogen bond (IAChB) was proposed and theoretically analyzed in a series of complexes formed by negatively charged bidentate chalcogen bond donors with chloride anions. The kinetic stability of IAChB was evidenced by the minima on binding energy profiles and further supported by ab initio molecular dynamic simulations. The block-localized wave function (BLW) method and its subsequent energy decomposition (BLW-ED) approach were employed to elucidate the physical origin of IAChB. While all other energy components vary monotonically as anions get together, the electrostatic interaction behaves exceptionally as it experiences a Coulombic repulsion barrier. Before reaching the barrier, the electrostatic repulsion increases with the shortening Ch⋯Cl^- distance as expected from classical electrostatics. However, after passing the barrier, the electrostatic repulsion decreases with the Ch⋯Cl^- distance shortening and subsequently turns into the most favorable trend among all energy terms at short ranges, representing a dominating force for the kinetic stability of inter-anions. For comparison, all energy components exhibit the same trends and vary monotonically in the conventional counterparts where donors are neutral. By comparing inter-anions and their conventional counterparts, we found that only the electrostatic energy term is affected by the extra negative charge. Remarkably, the distinctive (nonmonotonic) electrostatic energy profiles were reproduced using quantum mechanical-based atomic multipoles, suggesting that the crucial electrostatic interaction in IAChB can be rationalized within the classical electrostatic theory just like conventional non-covalent interactions.

Design of Azobenzene beyond Simple On-Off Behavior

Azobenzenes are without a doubt the most widely used light-induced switching units, and there is a plethora of application examples ranging from supramolecular chemistry to material science and biological chemistry. Here, we present a smart azobenzene, in which the photoswitching capability of the azobenzene moiety can be reversibly switched on and off using a second unit (redox switch). This second switching unit is based on the variation of the strength of a chalcogen bond between the azo group and a Te-Ph unit in ortho position to the azo group. This allows the selective switching of only one azobenzene unit in the presence of other azobenzene switches. The entire double-switch is a very simple, small system that can also be easily synthesized. As a result, this double-switch can be used as a smarter replacement for the established azobenzene system in the future. For example, in contrast to the latter this double-switch could be employed to store state information analogous to a flip-flop in digital electronic systems.

On the Importance of Pnictogen and Chalcogen Bonding Interactions in Supramolecular Catalysis

In this review, several examples of the application of pnictogen (Pn) (group 15) and chalcogen (Ch) bonding (group 16) interactions in organocatalytic processes are gathered, backed up with Molecular Electrostatic Potential surfaces of model systems. Despite the fact that the use of catalysts based on pnictogen and chalcogen bonding interactions is taking its first steps, it should be considered and used by the scientific community as a novel, promising tool in the field of organocatalysis.

Isothiourea-Catalyzed Enantioselective α-Alkylation of Esters via 1,6-Conjugate Addition to para-Quinone Methides

The isothiourea-catalyzed enantioselective 1,6-conjugate addition of para-nitrophenyl esters to 2,6-disubstituted para-quinone methides is reported. para-Nitrophenoxide, generated in situ from initial N-acylation of the isothiourea by the para-nitrophenyl ester, is proposed to facilitate catalyst turnover in this transformation. A range of para-nitrophenyl ester products can be isolated, or derivatized in situ by addition of benzylamine to give amides at up to 99% yield. Although low diastereocontrol is observed, the diastereoisomeric ester products are separable and formed with high enantiocontrol (up to 94:6 er).

On the Importance of σ–Hole Interactions in Crystal Structures

Elements from groups 14–18 and periods 3–6 commonly behave as Lewis acids, which are involved in directional noncovalent interactions (NCI) with electron-rich species (lone pair donors), π systems (aromatic rings, triple and double bonds) as well as nonnucleophilic anions (BF4−, PF6−, ClO4−, etc.). Moreover, elements of groups 15 to 17 are also able to act as Lewis bases (from one to three available lone pairs, respectively), thus presenting a dual character. These emerging NCIs where the main group element behaves as Lewis base, belong to the σ–hole family of interactions. Particularly (i) tetrel bonding for elements belonging to group 14, (ii) pnictogen bonding for group 15, (iii) chalcogen bonding for group 16, (iv) halogen bonding for group 17, and (v) noble gas bondings for group 18. In general, σ–hole interactions exhibit different features when moving along the same group (offering larger and more positive σ–holes) or the same row (presenting a different number of available σ–holes and directionality) of the periodic table. This is illustrated in this review by using several examples retrieved from the Cambridge Structural Database (CSD), especially focused on σ–hole interactions, complemented with molecular electrostatic potential surfaces of model systems.

Metal-Involving Chalcogen Bond. The Case of Platinum(II) Interaction with Se/Te-Based σ-Hole Donors

Platinum(II) complexes exhibiting an expressed d_z²-nucleophilicity of the positively charged metal centers, namely, [Pt(ppy)(acac)] (1; acacH is acetylacetone; ppyH is 2-Ph-pyridine) and [Pt(ppy)(tmhd)] (2; tmhdH is 2,2,6,6-tetramethylheptanedione-3,5), were cocrystallized with the chalcogen bond donors (4-NC₅F₄)₂Ch (Ch = Se, Te) to form two isostructural cocrystals 1·¹/₂(4-NC₅F₄)₂Ch, and 2·²/₃(4-NC₅F₄)₂Se and 2·(4-NC₅F₄)₂Te. The X-ray data for these cocrystals and appropriate theoretical DFT calculations (PBE0-D3BJ) allowed the recognition of the metal-involving chalcogen bond, namely, Ch···d_z²-Pt^II (its energy spans from -7 to -12 kcal/mol). In 1·¹/₂(4-NC₅F₄)₂Ch, Ch···d_z²-Pt^II bonding is accompanied by the C···d_z²-Pt^II interaction, representing a three-center bifurcate, whereas in 2·(4-NC₅F₄)₂Te the chalcogen bond Te···d_z²-Pt^II is purely two-centered and is stronger than that in 1·¹/₂(4-NC₅F₄)₂Ch because of more efficient orbital overlap. The association of 2 with (4-NC₅F₄)₂Te and the structure of the formed adduct in CDCl₃ solutions was studied by using ¹H, ¹³C, ¹⁹F, ¹⁹⁵Pt, ¹²⁵Te NMR, ¹⁹F-¹H HOESY, and diffusion NMR methods. The ¹⁹⁵Pt and ¹²⁵Te NMR titration and the isothermal titration calorimetry results revealed a 1:1 association of 2 with (4-NC₅F₄)₂Te.

Catalytic enantioselective synthesis of 1,4-dihydropyridines via the addition of C(1)-ammonium enolates to pyridinium salts

The regio- and stereoselective addition of C(1)-ammonium enolates - generated in situ from aryl esters and the isothiourea catalyst (R)-BTM - to pyridinium salts bearing an electron withdrawing substituent in the 3-position allows the synthesis of a range of enantioenriched 1,4-dihydropyridines. This represents the first organocatalytic approach to pyridine dearomatisation using pronucleophiles at the carboxylic acid oxidation level. Optimisation studies revealed a significant solvent dependency upon product enantioselectivity, with only toluene providing significant asymmetric induction. Using DABCO as a base also proved beneficial for product enantioselectivity, while investigations into the nature of the counterion showed that co-ordinating bromide or chloride substrates led to higher product er than the corresponding tetrafluoroborate or hexafluorophosphate. The scope and limitations of this process are developed, with enantioselective addition to 3-cyano- or 3-sulfonylpyridinium salts giving the corresponding 1,4-dihydropyridines (15 examples, up to 95 : 5 dr and 98 : 2 er).

Emission Intensity Enhancement for Iridium(III) Complex in Dimethyl Sulfoxide under Photoirradiation

We found emission intensity enhancement for fac-Ir(ppy)₃ (ppy = 2-(2'-phenyl)pyridine) in aerated dimethyl sulfoxide (DMSO) during photoirradiation for the first time. This phenomenon was concluded to be responsible for the consumption of ³O₂ dissolved in DMSO through dimethyl sulfone production by photosensitized reaction using fac-Ir(ppy)₃. A ³O₂ adduct of DMSO molecule was detected by UV absorption measurement and theoretical calculation. We proposed a mechanism for the emission enhancement reaction including ^1,3O₂ molecules and ^1,3O₂-DMSO adducts and validated it through a simulation of emission intensity change using an ordinary differential equation solver.

Classification of So-Called Non-Covalent Interactions Based on VSEPR Model

The variety of interactions have been analyzed in numerous studies. They are often compared with the hydrogen bond that is crucial in numerous chemical and biological processes. One can mention such interactions as the halogen bond, pnicogen bond, and others that may be classified as σ-hole bonds. However, not only σ-holes may act as Lewis acid centers. Numerous species are characterized by the occurrence of π-holes, which also may play a role of the electron acceptor. The situation is complicated since numerous interactions, such as the pnicogen bond or the chalcogen bond, for example, may be classified as a σ-hole bond or π-hole bond; it ultimately depends on the configuration at the Lewis acid centre. The disadvantage of classifications of interactions is also connected with their names, derived from the names of groups such as halogen and tetrel bonds or from single elements such as hydrogen and carbon bonds. The chaos is aggravated by the properties of elements. For example, a hydrogen atom can act as the Lewis acid or as the Lewis base site if it is positively or negatively charged, respectively. Hence names of the corresponding interactions occur in literature, namely hydrogen bonds and hydride bonds. There are other numerous disadvantages connected with classifications and names of interactions; these are discussed in this study. Several studies show that the majority of interactions are ruled by the same mechanisms related to the electron charge shifts, and that the occurrence of numerous interactions leads to specific changes in geometries of interacting species. These changes follow the rules of the valence-shell electron-pair repulsion model (VSEPR). That is why the simple classification of interactions based on VSEPR is proposed here. This classification is still open since numerous processes and interactions not discussed in this study may be included within it.

Chalcogen vs Halogen Bonding Catalysis in a Water-Bridge-Cocatalyzed Nitro-Michael Reaction

Recently, a tellurium-based chalcogen-bond-catalyzed nitro-Michael reaction was reported ( Angew. Chem. Int. Ed. 2019, 58, 16923), taking advantage of the strong Lewis acidity of the catalyst. This species was found to be more effective than an analogous iodine-based halogen bond organocatalyst. Herein, we present a detailed mechanistic and kinetic analysis of these catalytic cycles including the influence of the solvent (and the performance of different intrinsic solvation models). While the chalcogen bonding interaction is fundamental to activate the C-C bond formation, we found that the presence of a two-water molecular bridge is critical to allow the following, otherwise high-energy proton transfer step. Even though the iodine-based halogen bonding interaction is stronger than the tellurium-based chalcogen bonding one, which makes the former a stronger Lewis acid and hence in principle a more efficient catalyst, solvation effects explain the smaller energy span of the latter.

Triangular Interchalcogen Interactions: A Joint Crystallographic Data Analysis and Theoretical Study

Noncovalent chalcogen-chalcogen interactions are being actively investigated from both crystallographic and theoretical viewpoints in recent years. According to our search of the Cambridge Structural Database (CSD), a huge number of crystal structures containing triangular Ch₃ synthons were extracted. On the basis of the results of the CSD survey, a series of trimeric complexes of organic divalent chalcogen molecules were selected to model such interaction motifs. Similar to that in conventional chalcogen bonds, triangular interchalcogen interactions become gradually stronger along the sequence of Ch = S, Se, Te. Particularly, hydrogen bonds between the chalcogen centers and the H atoms in the substituents occur, which play a significant role in stabilizing the Ch₃ motifs in the trimers. Through multibody energy calculations, the complexes under study exhibit no or only weak cooperativity. Finally, the differences between the Ch₃ interaction motifs and the well-studied windmill X₃ bonding (X means halogen and this is halogen bond) patterns were elucidated.

Enantioselective Synthesis of α-Aryl-β2 -Amino-Esters by Cooperative Isothiourea and Brønsted Acid Catalysis

The synthesis of α-aryl-β2 -amino esters through enantioselective aminomethylation of an arylacetic acid ester in high yields and enantioselectivity via cooperative isothiourea and Brønsted acid catalysis is demonstrated. The scope and limitations of this process are explored (25 examples, up to 94 % yield and 96:4 er), with applications to the synthesis of (S)-Venlafaxine⋅HCl and (S)-Nakinadine B. Mechanistic studies are consistent with a C(1)-ammonium enolate pathway being followed rather than an alternative dynamic kinetic resolution process. Control studies indicate that (i) a linear effect between catalyst and product er is observed; (ii) an acyl ammonium ion can be used as a precatalyst; (iii) reversible isothiourea addition to an in situ generated iminium ion leads to an off-cycle intermediate that can be used as a productive precatalyst.

CHAL336 Benchmark Set: How Well Do Quantum-Chemical Methods Describe Chalcogen-Bonding Interactions?

We present the CHAL336 benchmark set-the most comprehensive database for the assessment of chalcogen-bonding (CB) interactions. After careful selection of suitable systems and identification of three high-level reference methods, the set comprises 336 dimers each consisting of up to 49 atoms and covers both σ- and π-hole interactions across four categories: chalcogen-chalcogen, chalcogen-π, chalcogen-halogen, and chalcogen-nitrogen interactions. In a subsequent study of DFT methods, we re-emphasize the need for using proper London dispersion corrections when treating noncovalent interactions. We also point out that the deterioration of results and systematic overestimation of interaction energies for some dispersion-corrected DFT methods does not hint at problems with the chosen dispersion correction but is a consequence of large density-driven errors. We conclude this work by performing the most detailed DFT benchmark study for CB interactions to date. We assess 109 variations of dispersion-corrected and dispersion-uncorrected DFT methods and carry out a detailed analysis of 80 of them. Double-hybrid functionals are the most reliable approaches for CB interactions, and they should be used whenever computationally feasible. The best three double hybrids are SOS0-PBE0-2-D3(BJ), revDSD-PBEP86-D3(BJ), and B2NCPLYP-D3(BJ). The best hybrids in this study are ωB97M-V, PW6B95-D3(0), and PW6B95-D3(BJ). We do not recommend using the popular B3LYP functional nor the MP2 approach, which have both been frequently used to describe CB interactions in the past. We hope to inspire a change in computational protocols surrounding CB interactions that leads away from the commonly used, popular methods to the more robust and accurate ones recommended herein. We would also like to encourage method developers to use our set for the investigation and reduction of density-driven errors in new density functional approximations.

α,β-Unsaturated acyl ammonium species as reactive intermediates in organocatalysis: an update

α,β-Unsaturated acyl ammonium species are versatile intermediates that have been applied in a variety of transformations including Michael additions, domino reactions and cycloadditions. Many of these transformations are promoted by chiral Lewis base catalysts, enabling the rapid generation of molecular complexity with high stereochemical control. This review highlights recent developments in the generation and application of α,β-unsaturated acyl ammonium intermediates reported since a previous review of this area in 2016. Particular emphasis will be placed on reports providing mechanistic insight into catalytic transformations and observed selectivities. A perspective on current challenges and potential future developments in the field of α,β-unsaturated acyl ammonium catalysis is also provided.

Quantification of Noncovalent Interactions in Azide-Pnictogen, -Chalcogen, and -Halogen Contacts

The noncovalent interactions between azides and oxygen‐containing moieties are investigated through a computational study based on experimental findings. The targeted synthesis of organic compounds with close intramolecular azide–oxygen contacts yielded six new representatives, for which X‐ray structures were determined. Two of those compounds were investigated with respect to their potential conformations in the gas phase and a possible significantly shorter azide–oxygen contact. Furthermore, a set of 44 high‐quality, gas‐phase computational model systems with intermolecular azide–pnictogen (N, P, As, Sb), –chalcogen (O, S, Se, Te), and –halogen (F, Cl, Br, I) contacts are compiled and investigated through semiempirical quantum mechanical methods, density functional approximations, and wave function theory. A local energy decomposition (LED) analysis is applied to study the nature of the noncovalent interaction. The special role of electrostatic and London dispersion interactions is discussed in detail. London dispersion is identified as a dominant factor of the azide–donor interaction with mean London dispersion energy‐interaction energy ratios of 1.3. Electrostatic contributions enhance the azide–donor coordination motif. The association energies range from −1.00 to −5.5 kcal mol⁻¹.

N-Oxide S-O chalcogen bonding in conjugated materials

Non-covalent bonding interactions, such as chalcogen bonding, can have a substantial effect on the electronic and physical properties of conjugated polymers and is largely dependent on the strength of interaction. Functional groups that are traditionally used to instill chalcogen bonding such as alkoxy or fluorine substituents can demand challenging synthetic effort, as well as have drastic effects on the electronics of a π-system. The incorporation of a N-oxide functionality into bithiazole-containing materials, a synthetically simple transformation, has been entirely overlooked until now. A systematic analysis of the effects of N-oxidation on the electronic and physical properties of bithiazole-containing materials has been undertaken. N-Oxidation has been found to affect the electronic band gap through increase of the HOMO and lowering of the LUMO. Furthermore, exceptionally strong intramolecular S-O chalcogen bonding interactions in the bithiazole core contribute to rigidification of the conjugated system. Computational analysis of this system has shown this N-oxide chalcogen bonding interaction to be significantly stronger than other chalcogen bonding interactions commonly exploited in conjugated materials.

Participation of S and Se in hydrogen and chalcogen bonds

The heavier chalcogen atoms S, Se, and Te can each participate in a range of different noncovalent interactions. They can serve as both proton donor and acceptor in H-bonds. Each atom can also act as electron acceptor in a chalcogen bond.

Geometrical and energetic characteristics of Se⋯Se interactions in crystal structures of organoselenium molecules

Combined crystallographic and quantum chemical study was performed to reveal the nature of selenium–selenium interactions in the crystal structures of organoselenium compounds.