Selenium chalcogen bonds are involved in protein-carbohydrate recognition: a combined PDB and theoretical study

Organoselenium Compounds in Medicinal Chemistry

The chemical and biological interest in this element and the molecules bearing selenium has been exponentially growing over the years. Selenium, formerly designated as a toxin, becomes a vital trace element for life that appears as selenocysteine and its dimeric form, selenocystine, in the active sites of selenoproteins, which catalyze a wide variety of reactions, including the detoxification of reactive oxygen species and modulation of redox activities. From the point of view of drug developments, organoselenium drugs are isosteres of sulfur-containing and oxygen-containing drugs with the advantage that the presence of the selenium atom confers antioxidant properties and high lipophilicity, which would increase cell membrane permeation leading to better oral bioavailability. This statement is the paramount relevance considering the big number of clinically employed compounds bearing sulfur or oxygen atoms in their structures including nucleosides and carbohydrates. Thus, in this article we have focused on the relevant features of the application of selenium in medicinal chemistry. With the increasing interest in selenium chemistry, we have attempted to highlight the most significant published data on this subject, mainly concentrating the analysis on the last years. In consequence, the recent advances of relevant pharmacological organoselenium compounds are discussed.

Competing electrophilic substitution and oxidative polymerization of arylamines with selenium dioxide

This article describes the detailed analysis of the reaction between arylamines, such as aniline, o-anisidine, and methyl anthranilate, with selenium dioxide in acetonitrile. A systematic analysis of the reaction products with the help of 77Se NMR and single-crystal X-ray crystallography revealed that the reaction progress follows three major reaction pathways, electrophilic selenation, oxidative polymerization, and solvent oxidation. For aniline and o-anisidine, predominant oxidative polymerization occurred, leading to the formation of the respective polyaniline polymers as major products. For methyl anthranilate, the oxidative polymerization was suppressed due to the delocalization of amine lone pair electrons over the adjacent carboxylate function, which prompted the selenation pathway, leading to the formation of two of the isomeric diorganyl selenides of methyl anthranilate. The diaryl selenides were structurally characterized using single-crystal X-ray diffraction. Density functional theory calculations suggest that the highest occupied molecular orbital of methyl anthranilate was deeply buried, which suppressed the oxidative polymerization pathway. Due to solvent oxidation, oxamide formation was also noticed to a considerable extent. This study provides that utmost care must be exercised while using SeO2 as an electrophile source in aromatic electrophilic substitution reactions.

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.

Exploiting π and Chalcogen Interactions for the β-Selective Glycosylation of Indoles through Glycal Conformational Distortion

Harnessing unconventional noncovalent interactions (NCIs) is emerging as a formidable synthetic approach in difficult-to-access glycosidic chemical space. C-Glycosylation, in particular, has gained a flurry of recent attention. However, most reported methods are restricted to the relatively facile access to α-C-glycosides. Herein, we disclose a β-stereoselective glycosylation of indoles by employing a phosphonoselenide catalyst. The robustness of this protocol is exemplified by its amenability for reaction at both the indolyl C- and N- reactivity sites. In contrast to previous reports, in which the chalcogens were solely involved in Lewis acidic activation, our mechanistic investigation unraveled that the often neglected flanking aromatic substituents of phosphonoselenides can substantially contribute to catalysis by engaging in π-interactions. Computations and NMR spectroscopy indicated that the chalcogenic and aromatic components of the catalyst can be collectively exploited to foster conformational distortion of the glycal away from the usual half-chair to the boat conformation, which liberates the convex β-face for nucleophilic attack.

Chalcogen bonding catalysis

Catalysts play a major role in chemical synthesis, and catalysis is considered to be a green and economic process. Catalysis is dominated by covalent interactions between the catalyst and substrate. The design of non-covalent catalysts came into limelight only recently. Hydrogen bonding (HB) catalysts are well established among non-covalent catalysts, including asymmetric HB catalysts. Though halogen bonding (XB) catalysis and its asymmetric version are gaining admiration, non-covalent chalcogen bonding catalysis (ChB) is in the budding stage. This tutorial review will focus on the recently evolved chalcogen bonding catalysis and emphasis will be given to the chalcogen bonding of chiral molecules. Since successful enantioselective chalcogen bonding catalysis is yet to be reported, this review will focus on the basics of non-covalent bonding catalysis, chalcogen bonding catalysis, chiral chalcogenide synthesis, rigidification of transition states by ChB, stabilization of cations by chiral chalcogens, details of unsuccessful asymmetric chalcogen bonding catalysis, enantioseparation of racemic molecules using ChB, and the existence of ChB in chiral biomolecules.

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.

Pd and Pt metal atoms as electron donors in σ-hole bonded complexes

Quantum calculations provide a systematic assessment of the ability of Group 10 transition metals M = Pd and Pt to act as an electron donor within the context of pnicogen, chalcogen, and halogen bonds. These M atoms are coordinated in a square planar geometry, attached to two N atoms of a modified phenanthrene unit, as well as two ligand atoms Cl, Br, or I. As the Lewis acid, a series of AFn molecules were chosen, which could form a pnicogen bond (A = P, As, Sb), chalcogen bond (A = S, Se, Te) or halogen bond (A = Cl, Br, I) with M. These noncovalent bonds are fairly strong, varying between 6 and 20 kcal mol-1, with the occupied dz2 orbital of M acting as the origin of charge transferred to the acid. Pt forms somewhat stronger bonds than Pd, and the bond strength rises with the size of the A atom of the acid. Within the context of smaller A atoms, the bond strength rises in the order pnicogen < chalcogen < halogen, but this distinction vanishes for the fifth-row A atoms. The nature of the ligand atoms on M has little bearing on the bond strength. Based on the Harmonic Oscillator Model of Aromaticity (HOMA) index, the ZB, YB and XB bonds were shown to have only a subtle effect on the ring electronic structures.

On the Importance of Halogen and Chalcogen Bonds in the Solid State of Nucleic Acids: A Combined Crystallographic and Theoretical Perspective

In this work, intra- and intermolecular halogen and chalcogen bonds (HlgBs and ChBs, respectively) present in the solid state of nucleic acids (NAs) have been studied at the RI-MP2/def2-TZVP level of theory. To achieve this, a Protein Data Bank (PDB) survey was carried out, revealing a series of structures in which Br/I or S/Se/Te atoms belonging to nucleobases or pentose rings were involved in noncovalent interactions (NCIs) with electron-rich species. The energetics and directionality of these NCIs were rationalized through a computational study, which included the use of Molecular Electrostatic Potential (MEP) surfaces, the Quantum Theory of Atoms in Molecules (QTAIM), and Non Covalent Interaction plot (NCIplot) and Natural Bonding Orbital (NBO) techniques.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Adjusting the balance between hydrogen and chalcogen bonds

A complex is assembled which pairs a carboxyl group of X₁COOH with a 1,2,5-chalcogenadiazole ring containing substituents on its C atoms. The OH of the carboxyl group donates a proton to a N atom of the ring to form a OH⋯N H-bond (HB), while its carbonyl O engages in a Y⋯O chalcogen bond (ChB) with the ring in which Y = S, Se, Te. The ChB is strengthened by enlarging the size of the Y atom from S to Se to Te. Placement of an electron-withdrawing group (EWG) X₁ on the acid strengthens the HB while weakening the ChB; the reverse occurs when EWGs are placed on the ring. By selection of the proper substituents on the two units, it is possible to achieve a near perfect balance between the strengths of these two bonds. These bond strengths are also reflected in the NMR spectroscopic properties of the chemical shielding of the various atoms and the coupling between the nuclei directly involved in each bond.

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.

Non-covalent interactions atlas benchmark data sets 4: σ-hole interactions

The SH250×10 dataset presented here extends the Non-Covalent Interactions Atlas database (https://www.nciatlas.org) to complexes bound by σ-hole interactions - halogen, chalcogen and pnictogen bonds. It comprises 250 complexes where Cl, Br, I, S, Se, P and As interact with diverse electron donors. An accurate CCSD(T)/CBS benchmark is provided for ten points along a dissociation curve of each complex. The SH250×10 set is used in testing a wide variety of DFT functionals and semiempirical quantum-mechanical methods. In DFT calculations, the new data set exposes large errors of some functionals related to exaggerated charge transfer. The size and diversity of the data set have also been exploited in the reparametrization of a halogen-bond correction for the PM6 semiempirical method.

Host-guest complexes vs. supramolecular polymers in chalcogen bonding receptors: an experimental and theoretical study

We describe here a comparative study between two tripodal anion receptors based on selenophene as the binding motif. The receptors use benzene or perfluorobenzene as a spacer. The presence of the electron-withdrawing ring activates the selenium atom for anion recognition inducing the formation of self-assembled supramolecular structures in the presence of chloride or bromide anions, which are bonded by the cooperative action of hydrogen and chalcogen bonding interactions. DOSY NMR and DLS experiments provided evidence for the formation of the supramolecular structures only in the presence of a perfluorobenzene based anion receptor while the analogous benzene one shows the classical anion/receptor complex without the participation of the selenium atom. The energetic and geometric features of the complexes of both receptors with the Cl and Br anions have been studied in solution. These results combined with the molecular electrostatic potential (MEP) surface plots allow us to rationalize the quite different behaviors of both receptors observed experimentally.

Triel bonds within anionanion complexes

The ability of two anions to interact with one another is tested in the context of pairs of TrX₄^- homodimers, where Tr represents any of the triel atoms B, Al, Ga, In, or Tl, and X refers to a halogen substituent F, Cl, or Br. None of these pairs engage in a stable complex in the gas phase, but the situation reverses in water where the two monomers are held together by Tr⋯X triel bonds, complemented by stabilizing interactions between X atoms. Some of these bonds are quite strong, notably those involving TrF₄^-, with interaction energies surpassing 30 kcal mol^-1. Others are very much weaker, with scarcely exothermic binding energies. The highly repulsive electrostatic interactions are counteracted by large polarization energies.

Se⋅⋅⋅O/S and S⋅⋅⋅O Chalcogen Bonds in Small Molecules and Proteins: A Combined CSD and PDB Study

The importance of selenium-centered noncovalent chalcogen bonds represented as Se⋅⋅⋅A (A=O/S) has been explored for short directional contacts in small molecules and proteins. In addition, S⋅⋅⋅O centered contacts have been analyzed. Computational analyses involving the quantitative assessment of the associated energetics, the molecular electrostatic potentials (MEP), and electron density derived topological parameters, namely, quantum theory of atom in molecules (QTAIM) analyses, and NBO (natural bond orbital) based calculations, have been performed to unequivocally establish the strength, stability, and attractive role of chalcogen bonds in the solid-state. This investigation has been performed in molecules from both the Cambridge Structural Database (CSD) and Protein Data Bank (PDB). Thus futuristic materials may be designed keeping in mind the significance of these interactions, including their relevance in biology.