Spectroscopic Evidences for Strong Hydrogen Bonds with Selenomethionine in Proteins

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.

Characteristics of nonconventional hydrogen bonds and stability of dimers of chalcogenoaldehyde derivatives: a noticeable role of oxygen compared to other chalcogens

In this work, twenty-four stable dimers of RCHZ with R = H, F, Cl, Br, CH3 or NH2 and Z = O, S, Se or Te were determined. It was found that the stability of most dimers is primarily contributed by the electrostatic force, except for the dominant role of the induction term in those involving a Te atom, which has been rarely observed. Both electron-donating and -withdrawing groups in substituted formaldehyde cause an increase in the strength of nonconventional Csp2-H⋯Z hydrogen bonds, as well as the dimers, in which the electron donating effect plays a more crucial role. The strength of nonconventional hydrogen bonds decreases in the following order: Csp2-H⋯O ≫ Csp2-H⋯S > Csp2-H⋯Se > Csp2-H⋯Te. Remarkably, a highly significant role of the O atom compared to S, Se and Te in increasing the Csp2-H stretching frequency and strength of the nonconventional hydrogen bonds and dimers is found. A Csp2-H stretching frequency red-shift is observed in Csp2-H⋯S/Se/Te, while a blue-shift is obtained in Csp2-H⋯O. When Z changes from O to S to Se and to Te, the Csp2-H blue-shift tends to decrease and eventually turns to a red-shift, in agreement with the increasing order of the proton affinity at Z in the isolated monomer. The magnitude of the Csp2-H stretching frequency red-shift is larger for Csp2-H⋯Te than Csp2-H⋯S/Se, consistent with the rising trend of proton affinity at the Z site and the polarity of the Csp2-H bond in the substituted chalcogenoaldehydes. The Csp2-H blue-shifting of the Csp2-H⋯O hydrogen bonds is observed in all dimers regardless of the electron effect of the substituents. Following complexation, the electron-donating derivatives exhibit a stronger Csp2-H blue-shift compared to the electron-withdrawing ones. Notably, the stronger Csp2-H blue-shift turns out to involve a less polarized Csp2-H bond and a decrease in the occupation at the σ*(Csp2-H) antibonding orbital in the isolated monomer.

Organelle-Mediated Dissipative Self-Assembly of Peptides in Living Cells

Implementing dissipative assembly in living systems is meaningful for creation of living materials or even artificial life. However, intracellular dissipative assembly remains scarce and is significantly impeded by the challenges lying in precisely operating chemical reaction cycles under complex physiological conditions. Here, we develop organelle-mediated dissipative self-assembly of peptides in living cells fueled by GSH, via the design of a mitochondrion-targeting and redox-responsive hexapeptide. While the hexapeptide undergoes efficient redox-responsive self-assembly, the addition of GSH into the peptide solution in the presence of mitochondrion-biomimetic liposomes containing hydrogen peroxide allows for transient assembly of peptides. Internalization of the peptide by LPS-stimulated macrophages leads to the self-assembly of the peptide driven by GSH reduction and the association of the peptide assemblies with mitochondria. The association facilitates reversible oxidation of the reduced peptide by mitochondrion-residing ROS and thereby dissociates the peptide from mitochondria to re-enter the cytoplasm for GSH reduction. The metastable peptide-mitochondrion complexes prevent the thermodynamically equilibrated self-assembly, thus establishing dissipative assembly of peptides in stimulated macrophages. The entire dissipative self-assembling process allows for elimination of elevated ROS and decrease of pro-inflammatory cytokine expression. Creating dissipative self-assembling systems assisted by internal structures provides new avenues for the development of living materials or medical agents in the future.

Critical assessment of selenourea as an efficient small molecule fluorescence quenching probe to monitor protein dynamics

Organoselenium compounds have recently been the experimentalists' delight due to their broad applications in organic synthesis, medicinal chemistry, and materials science. Selenium atom replacement of the carbonyl oxygen of the urea moiety dramatically reduces the HOMO-LUMO gap and oxidation potential, which completely changes the physicochemical properties of selenocarbonyl compounds. To our surprise, the photophysics and utility of a simple molecule such as selenourea (SeU) have not been explored in detail, which persuaded us to investigate its role in excited state processes. The steady-state emission, temperature-dependent time-correlated single photon counting, and femtosecond fluorescence upconversion experimental results confirmed that SeU significantly enhances the fluorescence quenching through a photoinduced electron transfer (PET) mechanism with an ∼10 ps ultrafast intrinsic PET lifetime component which is mostly absent in thiourea (TU). A wide range of fluorophores, based on their different redox abilities and fluorescence lifetimes covering a broad spectral window (λex: 390-590 nm and λem: 490-690 nm), were chosen to validate the proof of the concept. It was extended to tetramethylrhodamine (TMR)-5-maleimide labeled lysozyme protein, where we observed significant fluorescence quenching in the presence of SeU. The present work emphasizes that the high quenching efficiency with an ultrafast PET process, reduced orbital energy gap, and higher negative free energy change of the electron transfer reaction are the representative characteristics of selenourea or selenoamides to enable them as potential surrogates of thioamides or oxoamides quenching probes to monitor protein conformational changes and dynamics.

Selenium in Proteins: Conformational Changes Induced by Se Substitution on Methionine, as Studied in Isolated Model Peptides by Optical Spectroscopy and Quantum Chemistry

The side-chain of methionine residues is long enough to establish NH⋯S H-bonds with neighboring carbonyl groups of the backbone, giving rise to so-called intra-residue 6^δ and inter-residue 7^δ H-bonds. The aim of the present article is to document how the substitution of sulfur with a selenium atom affects the H-bonding of the Met system. This was investigated both experimentally and theoretically by conformation-resolved optical spectroscopy, following an isolated molecule approach. The present work emphasizes the similarities of the Met and Sem residues in terms of conformational structures, energetics, NH⋯Se/S H-bond strength and NH stretch spectral shifts, but also reveals subtle behavior differences between them. It provides evidence for the sensitivity of the H-bonding network with the folding type of the Sem/Met side-chains, where a simple flip of the terminal part of the side-chain can induce an extra 50 cm⁻¹ spectral shift of the NH stretch engaged in a 7^δ NH⋯S/Se bond.

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.

Importance of water and intramolecular interaction governs substantial blue shift of Csp2–H stretching frequency in complexes between chalcogenoaldehydes and water

Geometrical structure, stability and cooperativity, and contribution of hydrogen bonds to the stability of complexes between chalcogenoaldehydes and water were thoroughly investigated using quantum chemical methods. The stability of the complexes increases significantly when one or more H₂O molecules are added to the binary system, whereas it decreases sharply going from O to S, Se, or Te substitution. The O–H⋯O H-bond is twice as stable as C_sp²–H⋯O and O–H⋯S/Se/Te H-bonds. It is found that a considerable blue-shift of C_sp²–H stretching frequency in the C_sp²–H⋯O H-bond is mainly determined by an addition of water into the complexes along with the low polarity of the C_sp²–H covalent bond in formaldehyde and acetaldehyde. The C_sp²–H stretching frequency shift as a function of net second hyperconjugative energy for the σ*(C_sp²–H) antibonding orbital is observed. Remarkably, a considerable C_sp²–H blue shift of 109 cm⁻¹ has been reported for the first time. Upon the addition of H₂O into the binary systems, halogenated complexes witness a decreasing magnitude of the C_sp²–H stretching frequency blue-shift in the C_sp²–H⋯O H-bond, whereas CH₃-substituted complexes experience the opposite trend.

The considerable blue shift of C_sp²–H stretching frequency.

Hydrogen Bonding with Polonium

Hydrogen bonding (H-bonding) with heavier chalcogens such as polonium and tellurium is almost unexplored owing to their lower electronegativities, providing us an opportunity to delve into the uncharted territory of...

Theoretical Aspects of Nonconventional Hydrogen Bonds in the Complexes of Aldehydes and Hydrogen Chalcogenides

Hydrogen bonds (H-bonds) in the complexes between aldehydes and hydrogen chalcogenides, XCHO...nH₂Z with X = H, F, Cl, Br, and CH₃, Z = O, S, Se, and Te, and n = 1,2, were investigated using high-level ab initio calculations. The C_sp2-H...O H-bonds are found to be about twice as strong as the C_sp2-H...S/Se/Te counterparts. Remarkably, the S/Se/Te-H...S/Se/Te H-bonds are 4.5 times as weak as the O-H...O ones. The addition of the second H₂Z molecule into binary systems induces stronger complexes and causes a positive cooperative effect in ternary complexes. The blue shift of C_sp2-H stretching frequency involving the C_sp2-H...Z H-bond sharply increases when replacing one H atom in HCHO by a CH₃ group. In contrast, when one H atom in HCHO is substituted with a halogen, the magnitude of blue-shifting of the C_sp2-H...Z H-bond becomes smaller. The largest blue shift up to 92 cm^-1 of C_sp2-H stretching frequency in C_sp2-H...O H-bond in CH₃CHO...2H₂O has rarely been observed and is much greater than that in the cases of the C_sp2-H...S/Se/Te ones. The C_sp2-H blue shift of C_sp2-H...Z bonds in the halogenated aldehydes is converted into a red shift when H₂O is replaced by a heavier analogue, such as H₂S, H₂Se, or H₂Te. The stability and classification of nonconventional H-bonds including C_sp2-H...Se/Te, Te-H...Te, and Se/Te-H...O have been established for the first time.

Comparing C2=O and C2=S Barbiturates: Different Hydrogen-Bonding Patterns of Thiobarbiturates in Solution and the Solid State

Carbonyl-centered hydrogen bonds with various strength and geometries are often exploited in materials to embed dynamic and adaptive properties, with the use of thiocarbonyl groups as hydrogen-bonding acceptors remaining only scarcely investigated. We herein report a comparative study of C2=O and C2=S barbiturates in view of their differing hydrogen bonds, using the 5,5-disubstituted barbiturate B and the thiobarbiturate TB as model compounds. Owing to the different hydrogen-bonding strength and geometries of C2=O vs. C2=S, we postulate the formation of different hydrogen-bonding patterns in C2=S in comparison to the C2=O in conventional barbiturates. To study differences in their association in solution, we conducted concentration- and temperature-dependent NMR experiments to compare their association constants, Gibbs free energy of association ∆Gassn., and the coalescence behavior of the N-H‧‧‧S=C bonded assemblies. In Langmuir films, the introduction of C2=S suppressed 2D crystallization when comparing B and TB using Brewster angle microscopy, also revealing a significant deviation in morphology. When embedded into a hydrophobic polymer such as polyisobutylene, a largely different rheological behavior was observed for the barbiturate-bearing PB compared to the thiobarbiturate-bearing PTB polymers, indicative of a stronger hydrogen bonding in the thioanalogue PTB. We therefore prove that H-bonds, when affixed to a polymer, here the thiobarbiturate moieties in PTB, can reinforce the nonpolar PIB matrix even better, thus indicating the formation of stronger H-bonds among the thiobarbiturates in polymers in contrast to the effects observed in solution.

O-H stretching frequency red shifts do not correlate with the dissociation energies in the dimethylether and dimethylsulfide complexes of phenol derivatives

In this perspective, we present a comprehensive report on the spectroscopic and computational investigations of the hydrogen bonded (H-bonded) complexes of Me2O and Me2S with seven para-substituted H-bond donor phenols. The salient finding was that although the dissociation energies, D0, of the Me2O complexes were consistently higher than those of the analogous Me2S complexes, the red-shifts in phenolic O-H frequencies, Δν(O-H), showed the exactly opposite trend. This is in contravention of the general perception that the red shift in the X-H stretching frequency in the X-HY hydrogen bonded complexes is a reliable indicator of H-bond strength (D0), a concept popularly known as the Badger-Bauer rule. This is also in contrast to the trend reported for the H-bonded complexes of H2S/H2O with several para substituted phenols of different pKa values wherein the oxygen centered hydrogen bonded (OCHB) complexes consistently showed higher Δν(O-H) and D0 compared to those of the analogous sulfur centered hydrogen bonded (SCHB) complexes. Our effort was to understand these intriguing observations based on the spectroscopic investigations of 1 : 1 complexes in combination with a variety of high level quantum chemical calculations. Ab initio calculations at the MP2 level and the DFT calculations using various dispersion corrected density functionals (including DFT-D3) were performed on counterpoise corrected surfaces to compute the dissociation energy, D0, of the H-bonded complexes. The importance of anharmonic frequency computations is underscored as they were able to correctly reproduce the observed trend in the relative OH frequency shifts unlike the harmonic frequency computations. We have attempted to find a unified correlation that would globally fit the observed red shifts in the O-H frequency with the H-bonding strength for the four bases, namely, H2S, H2O, Me2O, and Me2S, in this set of H-bond donors. It was found that the proton affinity normalized Δν(O-H) values scale very well with the H-bond strength.

Doubling Förster Resonance Energy Transfer Efficiency in Proteins with Extrinsic Thioamide Probes: Implications for Thiomodified Nucleobases

Designing a potential protein-ligand pair is pivotal, not only to track the protein structure dynamics, but also to assist in an atomistic understanding of drug delivery. Herein, the potential of a small model thioamide probe being used to study albumin proteins is reported. By monitoring the Förster resonance energy transfer (FRET) dynamics with the help of fluorescence spectroscopic techniques, a twofold enhancement in the FRET efficiency of 2-thiopyridone (2TPY), relative to that of its amide analogue, is observed. Molecular dynamics simulations depict the relative position of the free energy minimum to be quite stable in the case of 2TPY through noncovalent interactions with sulfur, which help to enhance the FRET efficiency. Finally, its application is shown by pairing thiouracils with protein. It is found that the site-selective sulfur atom substitution approach and noncovalent interactions with sulfur can substantially enhance the FRET efficiency, which could be a potential avenue to explore in the design of FRET probes to study the structure and dynamics of biomolecules.

Observation of an Unusually Large IR Red-Shift in an Unconventional S-H···S Hydrogen-Bond

The S-H···S non-covalent interaction is generally known as an extremely unconventional weak hydrogen-bond in the literature. The present gas-phase spectroscopic investigation shows that the S-H···S hydrogen-bond can be as strong as any conventional hydrogen-bond in terms of the IR red-shift in the stretching frequency of the hydrogen-bond donor group. Herein, the strength of the S-H···S hydrogen-bond has been determined by measuring the red-shift (∼150 cm^-1) of the S-H stretching frequency in a model complex of 2-chlorothiophenol and dimethyl sulfide using isolated gas-phase IR spectroscopy coupled with quantum chemistry calculations. The observation of an unusually large IR red-shift in the S-H···S hydrogen-bond is explained in terms of the presence of a significant amount of charge-transfer interactions in addition to the usual electrostatic interactions. The existence of ∼750 S-H···S interactions between the cysteine and methionine residues in 642 protein structures determined from an extensive Protein Data Bank analysis also indicates that this interaction is important for the structures of proteins.

Comparison of Bifurcated Halogen with Hydrogen Bonds

Bifurcated halogen bonds are constructed with FBr and FI as Lewis acids, paired with NH₃ and NCH bases. The first type considered places two bases together with a single acid, while the reverse case of two acids sharing a single base constitutes the second type. These bifurcated systems are compared with the analogous H-bonds wherein FH serves as the acid. In most cases, a bifurcated system is energetically inferior to a single linear bond. There is a larger energetic cost to forcing the single σ-hole of an acid to interact with a pair of bases, than the other way around where two acids engage with the lone pair of a single base. In comparison to FBr and FI, the H-bonding FH acid is better able to participate in a bifurcated sharing with two bases. This behavior is traced to the properties of the monomers, in particular the specific shape of the molecular electrostatic potential, the anisotropy of the orbitals of the acid and base that interact directly with one another, and the angular extent of the total electron density of the two molecules.

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.

Adsorption by Granular Activated Carbon and Nano Zerovalent Iron from Wastewater: A Study on Removal of Selenomethionine and Selenocysteine

Selenomethionine (SeMet) and selenocysteine (SeCys) are the most common forms of organic selenium, which is often found in the effluent of industrial wastewater. These organic selenium compounds are toxic, bioavailable and most likely to bioaccumulate in aquatic organisms. This study investigated the use of two adsorbent candidates (granular activated carbon (GAC) and nano zerovalent iron (nZVI)) as treatment technologies for SeMet and SeCys removal. Batch experiments were performed and inductively coupled plasma optical emission spectrometer (ICP-OES) was used for sample analysis. Experimental data showed GAC demonstrated a higher affinity towards the removal of SeMet and SeCys compared to nZVI. The removal efficiency of SeCys and SeMet by GAC was 96.1% and 86.7%, respectively. NZVI adsorption capacity for SeCys was 39.4% and SeMet < 1.1%. Irrespective of the adsorbent, SeMet is more refractory to be adsorbed compared to SeCys. Kinetics data of GAC and nZVI agreed well with the pseudo-second-order model (R2 > 0.990). The experimental data of SeCys was characterized by Langmuir model, indicating monolayer adsorption. The adsorption capacity of nZVI for SeCys increased significantly by about 35%, with a decrease in pH from 9.0 to 4.0, indicating that SeCy removal by nZVI is pH dependent. While electrostatic attraction is considered the driving mechanism for nZVI adsorption, GAC uptake capacity is controlled by weak van der Waal forces. The adsorption of binary adsorbates (SeMet and SeCys) exhibited an inhibitory effect due to the competitive interaction between contaminant molecules.

Bis(silanetellurone) with C-H···Te Interaction

Herein, we report the synthesis of a series of bis(silanechalcogenones) [Ch = Te (2), S (3), or Se (4)] using an N-heterocyclic silylene-based SiCSi pincer ligand (1). 2 is the first example of a bis(silanetellurone) derivative. The bonding patterns of 2-4 were extensively studied by natural bond orbital, quantum theory of atoms in molecules, and noncovalent interaction index analyses, and these exhibit weak C-H···Ch interaction. The analogous reaction of 1 with trimethyl N-oxide produced a novel bis(cyclosiloxane) derivative (5). All of the complexes are duly characterized by single-crystal X-ray diffraction studies, multinuclear nuclear magnetic resonance (¹H, ¹³C, and ²⁹Si) spectroscopy, and high-resolution mass spectrometry.

Implication of Threonine-Based Ionic Liquids on the Structural Stability, Binding and Activity of Cytochrome c

Ionic liquids (ILs) are useful in pharmaceutical industries and biotechnology as alternative solvents or sources for protein extraction and purification, preservation of biomolecules and for regulating the catalytic activity of enzymes. However, the binding mechanism, the non-covalent forces responsible for protein-IL interactions and dynamics of proteins in IL need to be investigated in depth for the effective use of ILs as alternatives. Herein, we disclose the molecular level understanding of the structural intactness and reactivity of a model protein cytochrome c (Cyt c) in biocompatible threonine-based ILs with the help of experimental techniques such as isothermal titration calorimetry (ITC), fluorescence spectroscopy, transmission electron microscopy (TEM) as well as molecular docking. Hydrophobic and electrostatic forces are responsible for the structural and conformational integrity of Cyt c in IL. The ITC experiments revealed the Cyt c-IL binding free energies are in the range of 10-14 kJ/mol and the molecular docking studies demonstrated that ILs interact at the surfaces of Cyt c. The results look promising as the ILs used here are non-toxic and biocompatible, and thus may find potential applications in structural biology and biotechnology.

Non-conventional Hydrogen Bonding and Aromaticity: A Systematic Study on Model Nucleobases and Their Solvated Clusters

The conceptual development of aromaticity is essential to rationalize and understand the structure and behavior of aromatic heterocycles. This work addresses for the first time, the interconnection between aromaticity and sulfur/selenium centered hydrogen bonds (S/SeCHBs) involved in representative heterocycle models of canonical nucleobases (2-Pyridone; 2PY) and its sulfur (2-Thiopyridone; 2TPY) and selenium (2-Selenopyridone; 2SePY) analogs. The nucleus-independent chemical shift (NICS) and gauge induced magnetic current density (GIMIC) values suggested significant reduction of aromaticity upon replacement of exocyclic carbonyl oxygen with sulfur and selenium. However, we observed two-fold (57 %) and three-fold (80 %) enhancement in the aromaticity for 2TPY dimer, and 2SePY dimer, respectively which are connected through S/SeCHBs. Aromaticity enhancement was also noticed in 1 : 1 H-bonded complexes (heterodimers), micro hydrated clusters and for bulk hydration. It is expected that exocyclic S and Se incorporation into heterocycles without compromising aromatic loss would definitely reinforce to design new supramolecular building blocks via S/SeCH-bonded complexes.

The Prodigious Hydrogen Bonds with Sulfur and Selenium in Molecular Assemblies, Structural Biology, and Functional Materials

Hydrogen bonds (H-bonds) play important roles in imparting functionality to the basic molecules of life by stabilizing their structures and directing their interactions. Numerous studies have been devoted to understanding H-bonds involving highly electronegative atoms like nitrogen, oxygen, and halogens and consequences of those H-bonds in chemical reactions, catalysis, and structure and function of biomolecules; but the involvement of less electronegative atoms like sulfur and selenium in H-bond formation establishes the concept of noncanonical H-bonds. Initially belittled for the "weak" nature of their interactions, these perceptions have gradually evolved over time through dedicated efforts by several research groups. This has been facilitated by advancements in experimental methods for their detection through gas-phase laser spectroscopy and solution NMR spectroscopy, as well as through theoretical predictions from high level quantum chemical calculations.In this Account, we present insights into the versatility of the sulfur and selenium centered H-bonds (S/SeCHBs) by highlighting their multifarious applications in various fields from chemical reactions to optoelectronic properties to structural biology. Our group has highlighted the significance and strength of such H-bonds in natural and modified biomolecules. Here, we have reviewed several molecular assemblies, biomolecules, and functional materials, where the role of these H-bonds is pivotal in influencing biological functions. It is worth mentioning here that the precise experimental data obtained from gas-phase laser spectroscopy have contributed considerably to changing the existing perceptions toward S/SeCHBs. Thus, molecular beam experiments, though difficult to perform on smaller model thio- or seleno-substituted Molecules, etc. (amides, nucleobases, drug molecules), are inevitable to gather elementary knowledge and convincing concepts on S/SeCHBs that can be extended from a small four-atom sulfanyl dimer to a large 14 kDa iron-sulfur protein, ferredoxin. These H-bonds can also tailor a fascinating array of molecular frameworks and design supramolecular assemblies by inter- and intralinking of individual "molecular Lego-like" units.The discussion is indeed intriguing when it turns to the usage of S/SeCHBs in facile synthetic strategies like tuning regioselectivity in reactions, as well as invoking phenomena like dual phosphorescence and chemiluminescence. This is in addition to our investigations of the dispersive nature of the hydrogen bond between metal hydrides and sulfur or selenium as acceptor, which we anticipate would lead to progress in the areas of proton and hydride transfer, as well as force-field design. This Account demonstrates how ease of fabrication, enhanced efficiency, and alteration of physicochemical properties of several functional materials is facilitated owing to the presence of S/SeCHBs. Our efforts have been instrumental in the evaluation of various S/SeCHBs in flue gas capture, as well as design of organic energy harvesting materials, where dipole moment and polarizability have important roles to play. We hope this Account invokes newer perspectives with regard to how H-bonds with sulfur and selenium can be adequately adopted for crystal engineering, for more photo- and biophysical studies with different spectroscopic methods, and for developing next-generation field-effect transistors, batteries, superconductors, and organic thin-film transistors, among many other multifunctional materials for the future.

C-HS interaction exhibits all the characteristics of conventional hydrogen bonds

This is a tale of a pair of a hydrogen bond donor and acceptor, namely the CH donor and sulphur acceptor, neither of which is a conventional hydrogen bond participant. Sulfur (S), being less electronegative (2.58) compared to its first row analogue oxygen (3.44), has not been considered as a potential HB acceptor for a long time. The C-HY (Y = HB acceptor) interaction has its own history of exhibiting omnidirectional shifts in the CH stretching frequency upon complex formation. Therefore, a systematic investigation of the C-HS interaction was the primary goal of the work presented here. Together with gas-phase vibrational spectroscopy and ab initio quantum chemical calculations, the nature and strength of the C-HS hydrogen bond (HB) have been investigated in the complexes of 1,2,4,5-tetracyanobenzene (TCNB) with various sulfur containing solvents. Despite the unconventional nature of both HB donor and HB acceptor (C-H and S, respectively), it was found that the C-HS hydrogen bond exhibits all the characteristics of the conventional hydrogen bond. The binding strength of the C-HS H-bond in these complexes was found to be comparable to that of the conventional hydrogen bonds. The unusual stabilities of these HBs have been mainly attributed to the attractive dispersion interaction.

The Se…S/N interactions as a possible mechanism of δ-aminolevulinic acid dehydratase enzyme inhibition by organoselenium compounds: A computational study

•

DPDS and PSA interacts with cysteine residues from AlaD active site.

•

The Se…S interactions could be involved in the δ-AlaD inhibition.

•

δ-AlaD from Cucumis sativus does not present cysteine residues in the active site.

•

Se…N interactions could be involved in the organoselenium action.

Organoselenium compounds present many pharmacological properties and are promising drugs. However, toxicological effects associated with inhibition of thiol-containing enzymes, such as the δ-aminolevulinic acid dehydratase (δ-AlaD), have been described. The molecular mechanism(s) by which they inhibit thiol-containing enzymes at the atomic level, is still not well known. The use of computational methods to understand the physical–chemical properties and biological activity of chemicals is essential to the rational design of new drugs. In this work, we propose an in silico study to understand the δ-AlaD inhibition mechanism by diphenyl diselenide (DPDS) and its putative metabolite, phenylseleninic acid (PSA), using δ-AlaD enzymes from Homo sapiens (Hsδ-AlaD), Drosophila melanogaster (Dmδ-AlaD) and Cucumis sativus (Csδ-AlaD). Protein modeling homology, molecular docking, and DFT calculations are combined in this study. According to the molecular docking, DPDS and PSA might bind in the Hsδ-AlaD and Dmδ-AlaD active sites interacting with the cysteine residues by Se^…S interactions. On the other hand, the DPDS does not access the active site of the Csδ-AlaD (a non-thiol protein), while the PSA interacts with the amino acids residues from the active site, such as the Lys291. These interactions might lead to the formation of a covalent bond, and consequently, to the enzyme inhibition. In fact, DFT calculations (mPW1PW91/def2TZVP) demonstrated that the selenylamide bond formation is energetically favored. The in silico data showed here are in accordance with previous experimental studies, and help us to understand the reactivity and biological activity of organoselenium compounds.

77Se NMR Probes the Protein Environment of Selenomethionine

Sulfur is critical for the correct structure and proper function of proteins. Yet, lacking a sensitive enough isotope, nuclear magnetic resonance (NMR) experiments are unable to deliver for sulfur in proteins the usual wealth of chemical, dynamic, and structural information. This limitation can be circumvented by substituting sulfur with selenium, which has similar physicochemical properties and minimal impact on protein structures but possesses an NMR compatible isotope (⁷⁷Se). Here we exploit the sensitivity of ⁷⁷Se NMR to the nucleus' chemical milieu and use selenomethionine as a probe for its proteinaceous environment. However, such selenium NMR spectra of proteins currently resist a reliable interpretation because systematic connections between variations of system variables and changes in ⁷⁷Se NMR parameters are still lacking. To start narrowing this knowledge gap, we report here on a biological ⁷⁷Se magnetic resonance data bank based on a systematically designed library of GB1 variants in which a single selenomethionine was introduced at different locations within the protein. We recorded the resulting isotropic ⁷⁷Se chemical shifts and relaxation times for six GB1 variants by solution-state ⁷⁷Se NMR. For four of the GB1 variants we were also able to determine the chemical shift anisotropy tensor of SeM by solid-state ⁷⁷Se NMR. To enable interpretation of the NMR data, the structures of five of the GB1 variants were solved by X-ray crystallography to a resolution of 1.2 Å, allowing us to unambiguously determine the conformation of the selenomethionine. Finally, we combine our solution- and solid-state NMR data with the structural information to arrive at general insights regarding the execution and interpretation of ⁷⁷Se NMR experiments that exploit selenomethionine to probe proteins.

Non-covalent interactions with inverted carbon: a carbo-hydrogen bond or a new type of hydrogen bond?

Crystal structure analysis and quantum chemical calculations enabled us to discover a new non-covalent interaction, coined as carbo-hydrogen bond (C_H-bond).

Selenoamides modulate dipole-dipole interactions in hydrogen bonded supramolecular polymers of 1,3,5-substituted benzenes

We report the synthesis and self-assembly behavior of a chiral C₃-symmetrical benzene-tricarboselenoamide. The introduction of the selenoamide moiety enhances the dipolar character of the supramolecular interaction and confers a remarkable thermal stability to the supramolecular polymers obtained.

Water-Mediated Selenium Hydrogen-Bonding in Proteins: PDB Analysis and Gas-Phase Spectroscopy of Model Complexes

High-resolution X-ray crystallography and two-dimensional NMR studies demonstrate that water-mediated conventional hydrogen-bonding interactions (N-H···N, O-H···N, etc.) bridging two or more amino acid residues contribute to the stability of proteins and protein-ligand complexes. In this work, we have investigated single water-mediated selenium hydrogen-bonding interactions (unconventional hydrogen-bonding) between amino acid residues in proteins through extensive protein data bank (PDB) analysis coupled with gas-phase spectroscopy and quantum chemical calculation of a model complex consisting of indole, dimethyl selenide, and water. Here, indole and dimethyl selenide represent the amino acid residues tryptophan and selenomethionine, respectively. The current investigation demonstrates that the most stable structure of the model complex observed in the IR spectroscopy mimics single water-mediated selenium hydrogen-bonded structural motifs present in the crystal structures of proteins. The present work establishes that water-mediated Se hydrogen-bonding interactions are ubiquitous in proteins and the number of these interactions observed in the PDB is more than that of direct Se hydrogen-bonds present there.

Superstrong Noncovalent Interface between Melamine and Graphene Oxide

There have been growing academic interests in the study of strong organic molecule-graphene [or graphene oxide (GO)] systems, owing to their essential noncovalent nature and the consequent chemomechanical behavior within the interface. A more recent experimental measurement [ Chem 2018, 4, 896-910] reported that the melamine-GO interface exhibits a remarkable noncovalent binding strength up to ∼1 nN, even comparable with typical covalent bonds. But the poor understanding on the complex noncovalent nature in particular makes it challenging to unveil the mystery of this high-performance interface. Herein, we carry out first-principles calculations to investigate the atomistic origin of ultrastrong noncovalent interaction between the melamine molecule and the GO sheet, as well as the chemomechanical synergy in interfacial behavior. The anomalous O-H···N hydrogen bonding, formed between the triazine moiety of melamine and the -OH in GO, is found cooperatively enhanced by the pin-like NH₂-π interaction, which is responsible for the strong interface. Following static pulling simulations validates the 1 nN level rupture strength and the contribution of each noncovalent interaction within the interface. Moreover, our results show that the -OH hydrogen bonding will mainly augment the interfacial adhesion strength, whereas the -NH₂ group cooperating with the -OH hydrogen bonding and conjugating with the GO surface will greatly improve the interfacial shear performance. Our work deepens the understanding on the chemomechanical behaviors within the noncovalent interface, which is expected to provide new potential strategies in designing high-performance graphene-based artificial nacreous materials.

Nature and Strength of M-H···S and M-H···Se (M = Mn, Fe, & Co) Hydrogen Bond

The significance of dispersion contribution in the formation of strong hydrogen bonds (H-bonds) can no more be ignored. It was illustrated that less electronegative and electropositive H-bond acceptors such as S, Se, and Te are also capable of forming strong N-H···Y H-bonds, mostly due to the high polarizabilities of H-bond acceptor atoms. Herein, for the first time, we report the evidence of formation of nonconventional M-H···Y H-bonds between metal hydrides (M-H, M = Mn, Fe, Co) and chalcogen H-bond acceptors (Y = O, S, or Se). The nature and the strength of unusual M-H···Y H-bonds were revealed by several quantum chemical calculations and H-bond descriptors. The structural parameters, electron density topology, donor-acceptor natural bond orbital (NBO) interaction energies, and spectroscopic observables such as M-H stretching frequencies and ¹H chemical shifts are well-correlated to manifest the existence and strength of M-H···Y H-bonding. The M-H···Y H-bonds are dispersive in nature, and the computed H-bond energies are found to be in the range from ∼5 to 30 kJ/mol, which can be compared to those of the conventional H-bonds such as O-H···O, N-H···O, and N-H···O═C H-bonds, etc.

The nature of selenium hydrogen bonding: gas phase spectroscopy and quantum chemistry calculations

Subsequent to the recent re-definition of hydrogen bonding by the IUPAC committee, there has been a growing search for finding the presence of this ever interesting non-covalent interaction between a hydrogen atom in an X-H group and any other atom in the periodic table. In recent gas phase experiments, it has been observed that hydrogen bonding interactions involving S and Se are of similar strength to those with an O atom. However, there is no clear explanation for the unusual strength of this interaction in the case of hydrogen bond acceptors which are not conventional electronegative atoms. In this work, we have explored the nature of Se hydrogen bonding by studying indoledimethyl selenide (indmse) and phenoldimethyl selenide (phdmse) complexes using gas phase IR spectroscopy and quantum chemistry calculations. We have found through various energy decomposition analysis (EDA) methods and natural bond orbital (NBO) calculations that, along with electrostatics and polarization, charge transfer interactions are important to understand Se/S hydrogen bonding and there is a delicate balance between the various interactions that plays the crucial role rather than a single component of the interaction energy. An in-depth understanding of this type of non-covalent interaction has immense significance in biology as amino acids containing S and Se are widely present in proteins and hence hydrogen bonding interactions involving S and Se atoms contribute to the folding of proteins.

Thioamide, a Hydrogen Bond Acceptor in Proteins and Nucleic Acids

Thioamides are used as potential surrogates of amides to study the structure and dynamics of proteins and nucleic acids. However, incorporation of thioamides in biomolecules leads to changes in their structures and conformations mostly attributed to the strength of the amide-N-H···S═C hydrogen bond. In most cases, it is considered weak owing to the small electronegativity of sulfur, and in some cases, it is as strong as conventional H-bonds. Herein, adopting PDB structure analysis, NMR spectroscopy, and quantum chemistry calculations, we have shown that thioamides in a geometrical and structural constraint-free environment are capable of forming strong H-bonds like their amide counterparts. These studies also enabled us to determine the amide-N-H···S═C H-bond enthalpy (ΔH) very precisely. The estimated ΔH for the amide-N-H···S═C H-bond is ∼-30 kJ/mol, which suggests that the amide-N-H···S═C H-bond is a strong H-bond and merits its inclusion in computational force fields for biomolecular structure simulations to explore the role of amide-N-H···S═C H-bonds in nucleobase pairing and protein folding.