Chalcogen Bonding in Protein-Ligand Complexes: PDB Survey and Quantum Mechanical Calculations

Reactivity of Zinc Fingers in Oxidizing Environments: Insight from Molecular Models Through Activation Strain Analysis

The reactivity of Zn2+ tetrahedral complexes with H2O2 was investigated in silico, as a first step in their disruption process. The substrates were chosen to represent the cores of three different zinc finger protein motifs, i. e., a Zn2+ ion coordinated to four cysteines (CCCC), to three cysteines and one histidine (CCCH), and to two cysteines and two histidines (CCHH). The cysteine and histidine ligands were further simplified to methyl thiolate and imidazole, respectively. H2O2 was chosen as an oxidizing agent due to its biological role as a metabolic product and species involved in signaling processes. The mechanism of oxidation of a coordinated cysteinate to sulfenate-κS and the trends for the different substrates were rationalized through activation strain analysis and energy decomposition analysis in the framework of scalar relativistic Density Functional Theory (DFT) calculations at ZORA-M06/TZ2P ae // ZORA-BLYP-D3(BJ)/TZ2P. CCCC is oxidized most easily, an outcome explained considering both electrostatic and orbital interactions. The isomerization to sulfenate-κO was attempted to assess whether this step may affect the ligand dissociation; however, it was found to introduce a kinetic barrier without improving the energetics of the dissociation. Lastly, ligand exchange with free thiolates and selenolates was investigated as a trigger for ligand dissociation, possibly leading to metal ejection; molecular docking simulations also support this hypothesis.

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.

σ-Hole Site-Based Interactions within Hypervalent Pnicogen, Halogen, and Aerogen-Bearing Molecules with Lewis Bases: A Comparative Study

σ-Hole site-based interactions in the trigonal bipyramidal geometrical structure of hypervalent pnicogen, halogen, and aerogen-bearing molecules with pyridine and NCH Lewis bases (LBs) were comparatively examined. In this respect, the ZF5···, XF3O2···, and AeF2O3···LB complexes (where Z = As, Sb; X = Br, I; Ae = Kr, Xe; and LB = pyridine and NCH) were investigated. The electrostatic potential (EP) analysis affirmations outlined the occurrence of σ-holes on the systems under consideration with disparate magnitudes that increased according to the following order: AeF2O3 < XF3O2 < ZF5. In line with EP outcomes, the proficiency of σ-hole site-based interactions increased as the atomic size of the central atom increased with a higher favorability for the pyridine-based complexes over NCH-based ones. The interaction energy showed the most favorable negative values of -35.97, -44.53, and -56.06 kcal/mol for the XeF2O3···, IF3O2···, and SbF5···pyridine complexes, respectively. The preferentiality pattern of the studied interactions could be explained as a consequence of (i) the dramatic rearrangement of ZF5 molecules from the trigonal bipyramid geometry to the square pyramidal one, (ii) the significant and tiny deformation energy in the case of the interaction of XF3O2 molecules with pyridine and NCH, respectively, and (iii) the absence of geometrical deformation within the AeF2O3···pyridine and ···NCH complexes other than the XeF2O3···pyridine one. Quantum theory of atoms in molecules and noncovalent interaction index findings reveal the partially covalent nature of most of the investigated interactions. Symmetry-adapted perturbation theory affirmations declared that the electrostatic component was the driving force beyond the occurrence of the considered interactions. The obtained findings will help in improving our understanding of the effect of geometrical deformation on intermolecular interactions.

Synthesis of a 3,7-Disubstituted Isothiazolo[4,3-b]pyridine as a Potential Inhibitor of Cyclin G-Associated Kinase

Disubstituted isothiazolo[4,3-b]pyridines are known inhibitors of cyclin G-associated kinase. Since 3-substituted-7-aryl-isothiazolo[4,3-b]pyridines remain elusive, a strategy was established to prepare this chemotype, starting from 2,4-dichloro-3-nitropyridine. Selective C-4 arylation using ligand-free Suzuki-Miyaura coupling and palladium-catalyzed aminocarbonylation functioned as key steps in the synthesis. The 3-N-morpholinyl-7-(3,4-dimethoxyphenyl)-isothiazolo[4,3-b]pyridine was completely devoid of GAK affinity, in contrast to its 3,5- and 3,6-disubstituted congeners. Molecular modeling was applied to rationalize its inactivity as a GAK ligand.

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 bonds formed by protein sulfur atoms in proteins. A survey of high-resolution structures deposited in the protein data bank

The presence of chalcogen bonds in native proteins was investigated on a non-redundant and high-resolution (≤ 1 Angstrom) set of protein crystal structures deposited in the Protein Data Bank. It was observed that about one half of the sulfur atoms of methionines and disulfide bridges from chalcogen bonds with nucleophiles (oxygen and sulfur atoms, and aromatic rings). This suggests that chalcogen bonds are a non-bonding interaction important for protein stability. Quite numerous chalcogen bonds involve water molecules. Interestingly, in the case of disulfide bridges, chalcogen bonds have a marked tendency to occur along the S-S bond extension rather than along the C-S bond extension. Additionally, it has been observed that closer residues have a higher probability of being connected by a chalcogen bonds, while the secondary structure of the two residues connected by a chalcogen bond do not correlate with its formation.Communicated by Ramaswamy H. Sarma.

Evidence for and evaluation of fluorine-tellurium chalcogen bonding

In the field of noncovalent interactions, chalcogen bonding (ChB) involving the tellurium atom is currently attracting much attention in supramolecular chemistry and in catalysis. However, as a prerequisite for its application, the ChB should be studied in solution to assess its formation and, if possible, to evaluate its strength. In this context, new tellurium derivatives bearing CH2F and CF3 groups were designed to exhibit Te⋯F ChB and were synthesized in good to high yields. In both types of compounds, Te⋯F interactions were characterized in solution by combining 19F, 125Te and HOESY NMR techniques. These Te⋯F ChBs were shown to contribute to the overall JTe-F coupling constants (94-170 Hz) measured in the CH2F- and CF3-based tellurium derivatives. Finally, a variable temperature NMR study allowed us to approximate the energy of the Te⋯F ChB, from 3 kJ mol-1 for the compounds with weak Te σ-holes to 11 kJ mol-1 for Te σ-holes activated by the presence of strong electron withdrawing substituents.

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.

Sulfur-mediated chalcogen versus hydrogen bonds in proteins: a see-saw effect in the conformational space

Divalent sulfur (S) forms a chalcogen bond (Ch-bond) via its σ-holes and a hydrogen bond (H-bond) via its lone pairs. The relevance of these interactions and their interplay for protein structure and function is unclear. Based on the analyses of the crystal structures of small organic/organometallic molecules and proteins and their molecular electrostatic surface potential, we show that the reciprocity of the substituent-dependent strength of the σ-holes and lone pairs correlates with the formation of either Ch-bond or H-bond. In proteins, cystines preferentially form Ch-bonds, metal-chelated cysteines form H-bonds, while methionines form either of them with comparable frequencies. This has implications for the positioning of these residues and their role in protein structure and function. Computational analyses reveal that the S-mediated interactions stabilise protein secondary structures by mechanisms such as helix capping and protecting free β-sheet edges by negative design. The study highlights the importance of S-mediated Ch-bond and H-bond for understanding protein folding and function, the development of improved strategies for protein/peptide structure prediction and design and structure-based drug discovery.

Strategies for the Design of PEDOT Analogues Unraveled: the Use of Chalcogen Bonds and σ-Holes

In this theoretical study, we set out to demonstrate the substitution effect of PEDOT analogues on planarity as an intrinsic indicator for electronic performance. We perform a quantum mechanical (DFT) study of PEDOT and analogous model systems and demonstrate the usefulness of the ωB97X-V functional to simulate chalcogen bonds and other noncovalent interactions. We confirm that the chalcogen bond stabilizes the planar conformation and further visualize its presence via the electrostatic potential surface. In comparison to the prevalent B3LYP, we gain 4-fold savings in computational time and simulate model systems of up to a dodecamer. Implications for design of conductive polymers can be drawn from the results, and an example for self-doped polymers is presented where modulation of the strength of the chalcogen bond plays a significant role.

Interplay between hydrogen and chalcogen bonds in cysteine

Protein structures are stabilized by several types of chemical interactions between amino acids, which can compete with each other. This is the case of chalcogen and hydrogen bonds formed by the thiol group of cysteine, which can form three hydrogen bonds with one hydrogen acceptor and two hydrogen donors and a chalcogen bond with a nucleophile along the extension of the CS bond. A survey of the Protein Data Bank shows that hydrogen bonds are about 40-50 more common than chalcogen bonds, suggesting that they are stronger and, consequently, prevail, though not always. It is also observed that frequently a thiol group that forms a chalcogen bond is also involved, as a hydrogen donor, in a hydrogen bond.

Chalcogen Bonds: How to Characterize Them in Solution?

Chalcogen bonds (ChBs) occur between molecules containing Lewis acidic chalcogen substituents and Lewis bases. Recently, ChB emerged as a pivotal interaction in solution-based applications such as anion recognition, anion transport and catalysis. However, before moving to applications, the involvement of ChB must be established in solution. In this Concept article, we provide a brief review of the currently available experimental investigations of ChB in solution.

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.

Recognition in the Domain of Molecular Chirality: From Noncovalent Interactions to Separation of Enantiomers

It is not a coincidence that both chirality and noncovalent interactions are ubiquitous in nature and synthetic molecular systems. Noncovalent interactivity between chiral molecules underlies enantioselective recognition as a fundamental phenomenon regulating life and human activities. Thus, noncovalent interactions represent the narrative thread of a fascinating story which goes across several disciplines of medical, chemical, physical, biological, and other natural sciences. This review has been conceived with the awareness that a modern attitude toward molecular chirality and its consequences needs to be founded on multidisciplinary approaches to disclose the molecular basis of essential enantioselective phenomena in the domain of chemical, physical, and life sciences. With the primary aim of discussing this topic in an integrated way, a comprehensive pool of rational and systematic multidisciplinary information is provided, which concerns the fundamentals of chirality, a description of noncovalent interactions, and their implications in enantioselective processes occurring in different contexts. A specific focus is devoted to enantioselection in chromatography and electromigration techniques because of their unique feature as "multistep" processes. A second motivation for writing this review is to make a clear statement about the state of the art, the tools we have at our disposal, and what is still missing to fully understand the mechanisms underlying enantioselective recognition.

Noncovalent interactions in proteins and nucleic acids: beyond hydrogen bonding and π-stacking

Understanding the noncovalent interactions (NCIs) among the residues of proteins and nucleic acids, and between drugs and proteins/nucleic acids, etc., has extraordinary relevance in biomolecular structure and function. It helps in interpreting the dynamics of complex biological systems and enzymatic activity, which is esential for new drug design and efficient drug delivery. NCIs like hydrogen bonding (H-bonding) and π-stacking have been researchers' delight for a long time. Prominent among the recently discovered NCIs are halogen, chalcogen, pnictogen, tetrel, carbo-hydrogen, and spodium bonding, and n → π* interaction. These NCIs have caught the imaginations of various research groups in recent years while explaining several chemical and biological processes. At this stage, a holistic view of these new ideas and findings lying scattered can undoubtedly trigger our minds to explore more. The present review attempts to address NCIs beyond H-bonding and π-stacking, which are mainly n → σ*, n → π* and σ → σ* type interactions. Five of the seven NCIs mentioned earlier are linked to five non-inert end groups of the modern periodic table. Halogen (group-17) bonding is one of the oldest and most explored NCIs, which finds its relevance in biomolecules due to the phase correction and inhibitory properties of halogens. Chalcogen (group 16) bonding serves as a redox-active functional group of different active sites of enzymes and acts as a nucleophile in proteases and phosphates. Pnictogen (group 15), tetrel (group 14), triel (group 13) and spodium (group 12) bonding does exist in biomolecules. The n → π* interactions are linked to backbone carbonyl groups and protein side chains. Thus, they are crucial in determining the conformational stability of the secondary structures in proteins. In addition, a more recently discovered to and fro σ → σ* type interaction, namely carbo-hydrogen bonding, is also present in protein-ligand systems. This review summarizes these grand epiphanies routinely used to elucidate the structure and dynamics of biomolecules, their enzymatic activities, and their application in drug discovery. It also briefs about the future perspectives and challenges posed to the spectroscopists and theoreticians.

Structural study of bioisosteric derivatives of 5-(1H-indol-3-yl)-benzotriazole and their ability to form chalcogen bonds

Recently, inter-est in the isosteric replacement of a nitro-gen atom to selenium, sulfur or oxygen atoms has been highlighted in the design of potential inhibitors for cancer research. In this context, the structures of 5-(1H-indol-3-yl)-2,1,3-benzotriazole derivatives [5-(1H-indol-3-yl)-2,1,3-benzo-thia-diazole (bS, C14H9N3S) and 5-(1H-indol-3-yl)-2,1,3-benzoxa-diazole (bO, C14H9N3O)], as well as a synthesis inter-mediate of the selenated bioisostere [5-[1-(benzensulfon-yl)-1H-indol-3-yl]-2,1,3-benzoselena-diazole (p-bSe, C20H13N3O2SSe)] were determined using single-crystal X-ray diffraction (SCXRD) analyses. Despite being analogues, different crystal packing, torsion angles and supra-molecular features were observed, depending on the substitution of the central atoms of the benzotriazole. In particular, chalcogen inter-actions were described in the case of p-bSe and not in the bS and bO derivatives. An investigation by ab initio computational methods was therefore conducted to understand the effect of the substitution on the ability to form chalcogen bonds and the flexibility of the compounds.

Type I-IV Halogen⋯Halogen Interactions: A Comparative Theoretical Study in Halobenzene⋯Halobenzene Homodimers

In the current study, unexplored type IV halogen⋯halogen interaction was thoroughly elucidated, for the first time, and compared to the well-established types I−III interactions by means of the second-order Møller−Plesset (MP2) method. For this aim, the halobenzene⋯halobenzene homodimers (where halogen = Cl, Br, and I) were designed into four different types, parodying the considered interactions. From the energetic perspective, the preference of scouted homodimers was ascribed to type II interactions (i.e., highest binding energy), whereas the lowest binding energies were discerned in type III interactions. Generally, binding energies of the studied interactions were observed to decline with the decrease in the σ-hole size in the order, C6H5I⋯IC6H5 > C6H5Br⋯BrC6H5 > C6H5Cl⋯ClC6H5 homodimers and the reverse was noticed in the case of type IV interactions. Such peculiar observations were relevant to the ample contributions of negative-belt⋯negative-belt interactions within the C6H5Cl⋯ClC6H5 homodimer. Further, type IV torsional trans → cis interconversion of C6H5X⋯XC6H5 homodimers was investigated to quantify the π⋯π contributions into the total binding energies. Evidently, the energetic features illustrated the amelioration of the considered homodimers (i.e., more negative binding energy) along the prolonged scope of torsional trans → cis interconversion. In turn, these findings outlined the efficiency of the cis configuration over the trans analog. Generally, symmetry-adapted perturbation theory-based energy decomposition analysis (SAPT-EDA) demonstrated the predominance of all the scouted homodimers by the dispersion forces. The obtained results would be beneficial for the omnipresent studies relevant to the applications of halogen bonds in the fields of materials science and crystal engineering.

Unusual chalcogen⋯chalcogen interactions in like⋯like and unlike YCY⋯YCY complexes (Y = O, S, and Se)

Chalcogen⋯chalcogen interactions were investigated within four types of like⋯like and unlike YCY⋯YCY complexes (where Y = O, S, or Se). A plethora of quantum mechanical calculations, including molecular electrostatic potential (MEP), surface electrostatic potential extrema, point-of-charge (PoC), quantum theory of atoms in molecules (QTAIM), noncovalent interaction (NCI), and symmetry-adapted perturbation theory-based energy decomposition analysis (SAPT-EDA) calculations, were executed. The energetic findings revealed a preferential tendency of the studied chalcogen-bearing molecules to engage in type I, II, III, or IV chalcogen⋯chalcogen interactions. Notably, the selenium-bearing molecules exhibited the most potent ability to favorably participate in all the explored chalcogen⋯chalcogen interactions. Among like⋯like complexes, type IV interactions showed the most favorable negative binding energies, whereas type III interactions exhibited the weakest binding energies. Unexpectedly, oxygen-containing complexes within type IV interactions showed an alien pattern of binding energies that decreased along with an increase in the chalcogen atomic size level. QTAIM analysis provided a solo BCP, via chalcogen⋯chalcogen interactions, with no clues as to any secondary ones. SAPT-EDA outlined the domination of the explored interactions by the dispersion forces and indicated the pivotal shares of the electrostatic forces, except type III σ-hole⋯σ-hole and di-σ-hole interactions. These observations demonstrate in better detail all the types of chalcogen⋯chalcogen interactions, providing persuasive reasons for their more intensive use in versatile fields related to materials science and drug design.

Photomodulation of ultrastable host–guest complexes in water and their application in light-controlled steroid release

Host–guest complexation of dithienylethene photoswitches with cucurbit[8]uril leads to photoresponsive binding pairs with picomolar affinity in water.

Enhancing chalcogen bonding by metal coordination

This manuscript shows that chalcogen bonding (ChB) interaction is enhanced by the coordination of the chalcogen atom to metal centers, as evidenced using DFT calculations (PBE0-D3/def2-TZVP level of theory). X-ray...

Chalcogen Bonds Involving Selenium in Protein Structures

Chalcogen bonds are the specific interactions involving group 16 elements as electrophilic sites. The role of chalcogen atoms as sticky sites in biomolecules is underappreciated, and the few available studies have mostly focused on S. Here, we carried out a statistical analysis over 3562 protein structures in the Protein Data Bank (PDB) containing 18 266 selenomethionines and found that Se···O chalcogen bonds are commonplace. These findings may help the future design of functional peptides and contribute to understanding the role of Se in nature.

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

In this manuscript the ability of selenium carbohydrates to undergo chalcogen bonding (ChB) interactions with protein residues has been studied at the RI-MP2/def2-TZVP level of theory. An inspection of the Protein Data Bank (PDB) revealed SeA (A = O, C and S) intermolecular contacts involving Se-pyranose ligands and ASP, TYR, SER and MET residues. Theoretical models were built to analyse the strength and directionality of the interaction together with "Atoms in Molecules" (AIM), Natural Bonding Orbital (NBO) and Non Covalent Interactions plot (NCIplot) analyses, which further assisted in the characterization of the ChBs described herein. We expect that the results from this study will be useful to expand the current knowledge regarding biological ChBs as well as to increase the visibility of the interaction among the carbohydrate chemistry community.

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

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

Supramolecular Self-Assembly Built by Weak Hydrogen, Chalcogen, and Unorthodox Nonbonded Motifs in 4-(4-Chlorophenyl)-3-[(4-fluorobenzyl)sulfanyl]-5-(thiophen-2-yl)-4H-1,2,4-triazole, a Selective COX-2 Inhibitor: Insights from X-ray and Theoretical Studies

A selective triazole-based COX-2 inhibitor, 4-(4-chlorophenyl)-3-[(4-fluorobenzyl)sulfanyl]-5-(thiophen-2-yl)-4H-1,2,4-triazole, C₁₉H₁₃ClFN₃S₂, has been synthesized, and its crystal structure was determined at 150 K. Single-crystal X-ray diffraction analysis revealed that the thiophene ring was disordered over two orientations. The crystal structure is stabilized by weak hydrogen and chalcogen bonds and unorthodox F···π and S···C(π) contacts. These noncovalent interactions cooperatively generate the supramolecular self-assembly in the crystalline state. The Hirshfeld surface and its associated two-dimensional (2D)-fingerprint plots were obtained to analyze the role of different noncovalent interactions in the crystal packing. Further, the enrichment ratio was obtained from different atom···atom pairs to calculate the propensity of these pairs to form noncovalent interactions. The strength of different dimeric motifs formed in the crystal structure and lattice energies was calculated by the PIXEL method. Furthermore, the topological analysis of the charge density of intermolecular interactions was described. A CSD survey of C–H···F hydrogen bond, C–S···Cl chalcogen bond, and unorthodox nonbonded contacts (F···π and S···C(π)) is presented. The title compound possesses selective inhibitory activity against human COX-2 enzyme rather than COX-1. The quantum mechanics (QM) polarized ligand docking analysis was used to predict the binding pose and study the title compound’s selectivity against COX-1/2 enzymes.

The Characteristics of Disulfide-Centered Hydrogen Bonds

The disulfide-centered hydrogen bonds in the three different model systems of diethyl disulfide⋅⋅⋅H₂ O/H₂ CO/HCONH₂ clusters were characterized by high-resolution Fourier transform microwave spectroscopy and quantum chemical computations. The global minimum energy structures for each cluster are experimentally observed and are characterized by one of the three different S-S⋅⋅⋅H-C/N/O disulfide-centered hydrogen bonds and two O⋅⋅⋅H-C hydrogen bonds. Non-covalent interaction and natural bond orbital analyses further confirm the experimental observations. The symmetry-adapted perturbation theory (SAPT) analysis reveals that electrostatic is dominant in diethyl disulfide⋅⋅⋅H₂ O/HCONH₂ clusters being consistent with normal hydrogen bonds, whilst dispersion takes over in diethyl disulfide⋅⋅⋅H₂ CO cluster. Our study gives accurate structural parameters for the disulfide bond involved non-covalent clusters providing important benchmarking data for the theoretical evaluation of more complex systems.

Comparison of ±σ-hole and ±R˙-hole interactions formed by tetrel-containing complexes: a computational study

For the first time, unconventional ±R˙-hole interactions were unveiled in tetrel-containing complexes. The nature and characteristics of ±R˙-hole interactions were explored relative to their ±σ-hole counterparts for ˙TF3⋯ and W-T-F3⋯B/R˙/A complexes (where T = C, Si, and Ge, W = H and F, B = Lewis bases, R˙ = free radicals, and A = Lewis acids). In an effort to thoroughly investigate such interactions, a plethora of quantum mechanical calculations, including molecular electrostatic potential (MEP), maximum positive electrostatic potential (V s,max), point-of-charge (PoC), interaction energy, symmetry adapted perturbation theory (SAPT), and reduced density gradient-noncovalent interaction (RDG-NCI) calculations, were applied. The most notable findings to emerge from this study are that (i) from the electrostatic perspective, the molecular stabilization energies of ˙TF3 and W-T-F3 monomers became more negative as the Lewis basicity increased, (ii) the most stable complexes were observed for the ones containing Lewis bases, forming -σ-hole and -R˙-hole interactions, and the interaction energies systematically increased in the order H-T-F3⋯B < ˙TF3⋯B < F-T-F3⋯B, (iii) contrariwise, the +σ-hole and +R˙-hole interactions with Lewis acids are more energetically favorable in the order F-T-F3⋯A < ˙TF3⋯A < H-T-F3⋯A, and (iv) generally, the dispersion force plays a key role in stabilizing the tetrel-containing complexes, jointly with the electrostatic and induction forces for the interactions with Lewis bases and acids, respectively. Concretely, the findings presented in this paper add to our understanding of the characteristics and nature of such intriguing interactions.

Experimental and Theoretical Studies of Dimers Stabilized by Two Chalcogen Bonds in the Presence of a N···N Pnicogen Bond

The structure of the 5,6-dichloro-2,1,3-benzoselenadiazole homodimer, obtained by adding the ligand, 4,5-dichloro-o-phenylenediamine, to the methanolic solution of SeCl₄, was determined by X-ray crystallography, augmented by Fourier transform infrared, Raman, and NMR spectroscopy. The binding motif involves a pair of Se···N chalcogen bonds, with a supplementary N···N pnicogen bond. Quantum calculations provide assessments of the strengths of the individual interactions as well as their contributing factors. All together, these three bonds compose a total interaction energy between 5.4 and 16.8 kcal/mol, with the larger chalcogen atom associated with the strongest interactions. Replacement of the Se atoms by S and Te analogues allows analysis of the dependence of these forces on the nature of the chalcogen atom. Calculations also measure the importance to the binding of the presence of a second N atom on each diazole unit as well as the substituted phenyl ring to which it is fused.

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.

σ-Hole and Lone-Pair Hole Interactions in Chalcogen-Containing Complexes: A Comparative Study

The potentiality of sp3-hybridized chalcogen-containing molecules to participate in lone-pair (lp) hole interactions was reported for the first time. lp hole interactions were characterized and compared to σ-hole ones for OF2 and SF2 molecules as a case study. Various quantum mechanical calculations, including molecular electrostatic potential (MEP), maximum positive electrostatic potential (V s,max), point of charge (PoC), symmetry-adapted perturbation theory (SAPT), quantum theory of atoms in molecule (QTAIM), and reduced density gradient-noncovalent interaction (RDG-NCI) calculations, were carried out. The more significant findings to emerge from this study are the following: (i) the V s,max calculation was proved to be an unreliable method to determine the precise σ-hole and lp hole locations. (ii) The maximum positive electrostatic potential of the σ hole and lp hole was found to be at the F-Chal···PoC angle (θ) of 180° and at the centroid of XYlp plane, respectively. (iii) Lewis basicity has a significant effect on the strength of σ-hole and lp hole interactions. (iv) The studied molecules more favorably interact with Lewis bases via the σ hole compared to the lp hole, and (v) stabilization of the σ-hole and lp hole interactions stems from the electrostatic and dispersion forces, respectively.

In Silico Identification of Potential Natural Product Inhibitors of Human Proteases Key to SARS-CoV-2 Infection

Presently, there are no approved drugs or vaccines to treat COVID-19, which has spread to over 200 countries and at the time of writing was responsible for over 650,000 deaths worldwide. Recent studies have shown that two human proteases, TMPRSS2 and cathepsin L, play a key role in host cell entry of SARS-CoV-2. Importantly, inhibitors of these proteases were shown to block SARS-CoV-2 infection. Here, we perform virtual screening of 14,011 phytochemicals produced by Indian medicinal plants to identify natural product inhibitors of TMPRSS2 and cathepsin L. AutoDock Vina was used to perform molecular docking of phytochemicals against TMPRSS2 and cathepsin L. Potential phytochemical inhibitors were filtered by comparing their docked binding energies with those of known inhibitors of TMPRSS2 and cathepsin L. Further, the ligand binding site residues and non-covalent interactions between protein and ligand were used as an additional filter to identify phytochemical inhibitors that either bind to or form interactions with residues important for the specificity of the target proteases. This led to the identification of 96 inhibitors of TMPRSS2 and 9 inhibitors of cathepsin L among phytochemicals of Indian medicinal plants. Further, we have performed molecular dynamics (MD) simulations to analyze the stability of the protein-ligand complexes for the three top inhibitors of TMPRSS2 namely, qingdainone, edgeworoside C and adlumidine, and of cathepsin L namely, ararobinol, (+)-oxoturkiyenine and 3α,17α-cinchophylline. Interestingly, several herbal sources of identified phytochemical inhibitors have antiviral or anti-inflammatory use in traditional medicine. Further in vitro and in vivo testing is needed before clinical trials of the promising phytochemical inhibitors identified here.

SQM/COSMO Scoring Function: Reliable Quantum-Mechanical Tool for Sampling and Ranking in Structure-Based Drug Design

Quantum mechanical (QM) methods have been gaining importance in structure-based drug design where a reliable description of protein-ligand interactions is of utmost significance. However, strategies i. e. QM/MM, fragmentation or semiempirical (SQM) methods had to be pursued to overcome the unfavorable scaling of QM methods. Various SQM-based approaches have significantly contributed to the accuracy of docking and improvement of lead compounds. Parametrizations of SQM and implicit solvent methods in our laboratory have been instrumental to obtain a reliable SQM-based scoring function. The experience gained in its application for activity ranking of ligands binding to tens of protein targets resulted in setting up a faster SQM/COSMO scoring approach, which outperforms standard scoring methods in native pose identification for two dozen protein targets with ten thousand poses. Recently, SQM/COSMO was effectively applied in a proof-of-concept study of enrichment in virtual screening. Due to its superior performance, feasibility and chemical generality, we propose the SQM/COSMO approach as an efficient tool in structure-based drug design.

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.

Direct investigation of chalcogen bonds by multinuclear solid-state magnetic resonance and vibrational spectroscopy

We report a multifaceted experimental and computational study of three self-complementary chalcogen-bond donors as well as a series of seven chalcogen bonded cocrystals.

Mechanochemical disulfide reduction reveals imprints of noncovalent sulfur⋯oxygen chalcogen bonds in protein-inspired mimics in aqueous solution

Chalcogen bonds in proteins are found to impact on the mechanochemical reduction of disulfide bridges in aqueous environments.

Chalcogen bonding of two ligands to hypervalent YF4 (Y = S, Se, Te, Po)

The ability of two NH₃ ligands to engage in simultaneous chalcogen bonds to a hypervalent YF₄ molecule, with Y = S, Se, Te, Po, is assessed via quantum calculations. The complex can take on one of two different geometries. The cis structure places the two ligands adjacent to one another in a pseudo-octahedral geometry, held there by a pair of σ-hole chalcogen bonds. The bases can also lie nearly opposite one another, in a distorted octahedron containing one π-hole and one strained σ-hole bond. The cis geometry is favored for Y = S, while Te, and Po tend toward the trans structure; they are nearly equally stable for Se. In either case, the binding energy rises rapidly with the size of the Y atom, exceeding 30 kcal mol^-1 for PoF₄.

Secondary Forces in Protein Folding

A complete inventory of the forces governing protein folding is critical for productive protein modeling, including structure prediction and de novo design, as well as understanding protein misfolding diseases of clinical significance. The dominant contributors to protein folding include the hydrophobic effect and conventional hydrogen bonding, along with Coulombic and van der Waals interactions. Over the past few decades, important additional contributors have been identified, including C-H···O hydrogen bonding, n→π* interactions, C5 hydrogen bonding, chalcogen bonding, and interactions involving aromatic rings (cation-π, X-H···π, π-π, anion-π, and sulfur-arene). These secondary contributions fall into two general classes: (1) weak but abundant interactions of the protein main chain and (2) strong but less frequent interactions involving protein side chains. Though interactions with high individual energies play important roles in specifying nonlocal molecular contacts and ligand binding, we estimate that weak but abundant interactions are likely to make greater overall contributions to protein folding, particularly at the level of secondary structure. Further research is likely to illuminate additional roles of these noncanonical interactions and could also reveal contributions yet unknown.

On the capability of metal–halogen groups to participate in halogen bonds

Halogens in a M–X bond are inhibited from forming a halogen bond but can do so in certain circumstances, with or without a σ-hole.

Nature and strength of chalcogen–π bonds

We have analyzed the chalcogen–π bonding mechanism in a systematic series of model systems through Kohn–Sham molecular orbital theory and a quantitative energy decomposition scheme.