Sulfur, not too far behind O, N, and C: SH...pi hydrogen bond

DFT investigations of phenyldithiafulvene dimers at different oxidation states

Oxidative dimerization of aryl-substituted dithiafulvenes (Ar-DTFs) presents an efficient C-C bond forming method for the preparation of diverse redox-active π-conjugated molecules and conductive polymers. Previous experimental data indicated a reaction pathway in which direct combination of two Ar-DTF radical cations is a key step. However, mechanistic details about how Ar-DTF dimers are formed under different oxidation states have not yet been clearly established prior to this work. The assembly of two Ar-DTF molecules generates a vast conformational and configurational landscape, which is quite complex but fundamentally important for understanding the dimerization mechanism. To cast a deep insight into this aspect, we have performed density functional theory (DFT) calculations at the M06-2X/Def2-SVP level of theory to thoroughly investigate the potential energy surfaces (PESs) of various dimers of a phenyl-substituted dithiafulvene (Ph-DTF) in the mixed-valence radical cation and dication states. Key stationary points in these PESs, including minimum-energy conformers (π-dimers and σ-dimers) as well as the transition states connected to them, were examined and compared. We have also calculated the binding energies of these dimers to evaluate the energetic driving forces for their formation. Based on our computational results, the roles that various Ph-DTF dimers play in different pathways of oxidative dimerization have been clarified.

Novel Quinoline-Based Thiosemicarbazide Derivatives: Synthesis, DFT Calculations, and Investigation of Antitubercular, Antibacterial, and Antifungal Activities

The discovery of new antimicrobial agents as a means of treating drug-resistant microbial pathogens is of utmost significance to overcome their immense risk to human well-being. The current investigation involves the development, synthesis, and assessment of the antimicrobial efficacy of novel quinoline derivatives incorporating a thiosemicarbazide functionality. To design the target compounds (QST1-QST14), we applied the molecular hybridization approach to link various thiosemicarbazides to the quinoline core with a sulfonyl group. Upon the synthesis and completion of structural characterization via spectroscopic techniques (1H NMR, 13C NMR, 15N NMR, IR, and HRMS), the title molecules were extensively evaluated for their potential antitubercular, antibacterial, and antifungal activities. N-(3-Chlorophenyl)-2-(quinolin-8-ylsulfonyl)hydrazine-1-carbothioamide (QST4), the most effective compound against Mycobacterium tuberculosis H37Rv, was also tested on isoniazid-resistant clinical isolates with katG and inhA promoter mutations. Based on molecular docking studies, QST4 was also likely to demonstrate its antimycobacterial activity through inhibition of the InhA enzyme. Furthermore, three derivatives (QST3, QST4, and QST10) with preferable antimicrobial and drug-like profiles were also shown to be nontoxic against human embryonic kidney (HEK) cells. All compounds were optimized by the density functional theory method using B3LYP with the 6-31+G(d,p) basis set. Structural analysis, natural bond orbital calculations of donor-acceptor interactions, molecular electrostatic potential analysis, and frontier molecular orbital analysis were carried out. Quantum chemical descriptors and charges on the atoms were determined to compare the strengths of the intramolecular hydrogen bonds formed and their stabilities. We determined that the sulfur atom forms a stronger intramolecular hydrogen bond than the nitrogen, oxygen, and fluorine atoms in these sulfonyl thiosemicarbazide derivatives.

Impacts of noncovalent interactions involving sulfur atoms on protein stability, structure, folding, and bioactivity

This review discusses the various types of noncovalent interactions in which sulfur atoms participate and their effects on protein stability, structure, folding and bioactivity. Current approaches and recommendations for modelling these noncovalent interactions (in terms of both geometries and interaction energies) are highlighted.

Dynamic Long-Range Interactions Influence Substrate Binding and Catalysis by Human Histidine Triad Nucleotide-Binding Proteins (HINTs), Key Regulators of Multiple Cellular Processes and Activators of Antiviral ProTides

Human histidine triad nucleotide-binding (hHINT) proteins catalyze nucleotide phosphoramidase and acyl-phosphatase reactions that are essential for the activation of antiviral proTides, such as Sofosbuvir and Remdesivir. hHINT1 and hHINT2 are highly homologous but exhibit disparate roles as regulators of opioid tolerance (hHINT1) and mitochondrial activity (hHINT2). NMR studies of hHINT1 reveal a pair of dynamic surface residues (Q62, E100), which gate a conserved water channel leading to the active site 13 Å away. hHINT2 crystal structures identify analogous residues (R99, D137) and water channel. hHINT1 Q62 variants significantly alter the steady-state k_cat and K_m for turnover of the fluorescent substrate (TpAd), while stopped-flow kinetics indicate that K_D also changes. hHINT2, like hHINT1, exhibits a burst phase of adenylation, monitored by fluorescent tryptamine release, prior to rate-limiting hydrolysis and nucleotide release. hHINT2 exhibits a much smaller burst-phase amplitude than hHINT1, which is further diminished in hHINT2 R99Q. Kinetic simulations suggest that amplitude variations can be accounted for by a variable fluorescent yield of the E·S complex from changes in the environment of bound TpAd. Isothermal titration calorimetry measurements of inhibitor binding show that these hHINT variants also alter the thermodynamic binding profile. We propose that these altered surface residues engender long-range dynamic changes that affect the orientation of bound ligands, altering the thermodynamic and kinetic characteristics of hHINT active site function. Thus, studies of the cellular roles and proTide activation potential by hHINTs should consider the importance of long-range interactions and possible protein binding surfaces far from the active site.

Insight into Spodium–π Bonding Characteristics of the MX2···π (M = Zn, Cd and Hg; X = Cl, Br and I) Complexes—A Theoretical Study

The spodium–π bonding between MX₂ (M = Zn, Cd, and Hg; X = Cl, Br, and I) acting as a Lewis acid, and C₂H₂/C₂H₄ acting as a Lewis base was studied by ab initio calculations. Two types of structures of cross (T) and parallel (P) forms are obtained. For the T form, the X–M–X axis adopts a cross configuration with the molecular axis of C≡C or C=C, but both of them are parallel in the P form. NCI, AIM, and electron density shifts analyses further, indicating that the spodium–π bonding exists in the binary complexes. Spodium–π bonding exhibits a partially covalent nature characterized with a negative energy density and large interaction energy. With the increase of electronegativity of the substituents on the Lewis acid or its decrease in the Lewis base, the interaction energies increase and vice versa. The spodium–π interaction is dominated by electrostatic interaction in most complexes, whereas dispersion and electrostatic energies are responsible for the stability of the MX₂⋯C₂F₂ complexes. The spodium–π bonding further complements the concept of the spodium bond and provides a wider range of research on the adjustment of the strength of spodium bond.

Carbon-Centered Hydrogen Bonds in Proteins

Hydrogen bonding (H-bonding) without lone pair(s) of electrons and π-electrons is a concept developed 2-3 years ago. H-bonds involving less electronegative tetrahedral carbon are beyond the classical concept of H-bonds. Herein, we present the first report on H-bonds with tetravalent carbons in proteins. A special bonding arrangement is needed to increase the negative charge density around the sp³-hybridized carbon atom. Therefore, less electronegative elements such as As and Mg, when bonded to sp³-C, enable the C-atoms as H-bond acceptors. Careful protein structure analysis aided by several quantum chemical calculations suggests that these H-bonds are weak to moderate in strength. We developed an empirical equation to estimate the C-H···C H-bond energy in proteins from the distances between the C- and H-atoms. In proteins, the binding energies range from -5.4 to -14.0 kJ/mol. The C-H···C H-bonds assist the substrate binding in proteins. We also explored the potential role of these carbon-centered H-bonds in C-H bond activation through σ-bond metathesis. To our surprise, contribution from these H-bonds is almost of similar magnitude as that from C-H···π H-bonds for C-H bond activation.

Structural Characterization of 6-Thioguanosine and Its Monohydrate in the Gas Phase

Detailed structural analysis of 6-thioguanosine (6TGs) in relation to its tautomerization and sugar conformation is performed in the gas phase using UV and IR spectroscopy combined with ab initio calculations. We have observed a thiol tautomer of 6TGs with its sugar moiety in the syn conformation that is stabilized by a strong intramolecular H-bonding between O5'H of the sugar and N3 atom of the guanine moiety. This observation is consistent with previous results for guanosine (Gs) in which the corresponding enol form is solely detected. We have also identified a monohydrate of 6TGs consisting of a thiol tautomer with the water linking guanine moiety and sugar OH group. It is demonstrated that hydration behavior of 6TGs is significantly different from that of Gs as a result of a weaker H-bonding ability of the thiol group.

Azide···Oxygen Interaction: A Crystal Engineering Tool for Conformational Locking

We have designed, synthesized and crystallized 36 compounds, each containing an azide group and an oxygen atom separated by three bonds. Crystal structure analysis revealed that each of these molecules adopts a conformation in which the azide and oxygen groups orient syn to each other with a short O ··· N b  contact. Geometry-optimized structures [using  M06-2X/6-311G(d,p) level of theory ] also showed the syn conformation in all 36 of these cases, suggesting that this not merely a crystal packing effect. Quantum topological analysis using Bader's Atoms in Molecules (AIM) theory revealed bond paths and bond critical points (BCP) in these structures suggesting its nature and energetics to be similar to weak hydrogen bonding. The NCI-RDG plot clearly revealed the attractive interaction consisting of electrostatic or dispersive components in all the 36 systems. NBO analysis suggested a weak orbital-relaxation (charge-transfer) contribution of energy for a few  (sp2) O-donor systems. Natural population analysis (NPA) and molecular electrostatic potential mapping (MESP) of these crystal structures further revealed the existence of favorable azide-oxygen interaction.   A CSD search indicated the frequent and consistent occurrence of this interaction and its role dictating the syn conformation of azide and oxygen in molecules where these groups are separated by 2-4 bonds.

A unified cost-effective method for the construction of reliable potential energy surfaces for H2S and H2O clusters

A DFT-based methodology has been used to construct the potential energy surface of H2S clusters up to pentamers. Geometrical parameters and energetics show very good agreement with the existing experimental and high-level theoretical results. Distinct stable conformers of three dimers, six trimers, eleven tetramers and twenty-three pentamers have been identified. Both S-HS H-bond and SS interactions are identified in dimers, trimers and pentamers, while no SS interactions could be found in any of the 11 tetramer conformers. The binding energies of the most stable dimer, trimer, tetramer and pentamer are -1.66, -5.21, -8.57 and -12.54 kcal mol-1, respectively. The PES has been found to be exceedingly flat and the energy gap between the most and the least stable conformers was found to be only 0.09, 2.13, 1.65 and 1.13 kcal mol-1, from the dimer to the pentamer, respectively. The proposed method has also been used for water clusters up to the pentamer. The results obtained were found to agree closely with the existing results. Only one conformer was found for the water dimer, whereas four, five and fifteen conformers were obtained for the trimer, tetramer and pentamer, respectively. Atoms in molecular calculations were found to corroborate with the geometric and energetic results for both clusters.

Contributions of methionine to recognition of trimethyllysine in aromatic cage of PHD domains: implications of polarizability, hydrophobicity, and charge on binding

Recognition of trimethyllysine (Kme3) by reader proteins is an important regulator of gene expression. This recognition event is mediated by an aromatic cage made up of 2-4 aromatic residues in the reader proteins that bind Kme3 via cation-π interactions. A small subset of reader proteins contain a methionine (Met) residue in place of an aromatic sidechain in the binding pocket. The unique role of sulfur in molecular recognition has been demonstrated in a number of noncovalent interactions recently, including interactions of thiols, thioethers, and sulfoxides with aromatic rings. However, the interaction of a thioether with an ammonium ion has not previously been investigated and the role of Met in binding Kme3 has not yet been explored. Herein, we systematically vary the Met in two reader proteins, DIDO1 and TAF3, and the ligand, Kme3 or its neutral analog tert-butyl norleucine (tBuNle), to determine the role of Met in the recognition of the cationic Kme3. Our studies demonstrate that Met contributes to binding via dispersion forces, with about an equal contribution to binding Kme3 and tBuNle, indicating that electrostatic interactions do not play a role. During the course of these studies, we also discovered that DIDO1 exhibits equivalent binding to tBuNle and Kme3 through a change in the mechanism of binding.

Molecular Tailoring Approach for the Estimation of Intramolecular Hydrogen Bond Energy

Hydrogen bonds (HBs) play a crucial role in many physicochemical and biological processes. Theoretical methods can reliably estimate the intermolecular HB energies. However, the methods for the quantification of intramolecular HB (IHB) energy available in the literature are mostly empirical or indirect and limited only to evaluating the energy of a single HB. During the past decade, the authors have developed a direct procedure for the IHB energy estimation based on the molecular tailoring approach (MTA), a fragmentation method. This MTA-based method can yield a reliable estimate of individual IHB energy in a system containing multiple H-bonds. After explaining and illustrating the methodology of MTA, we present its use for the IHB energy estimation in molecules and clusters. We also discuss the use of this method by other researchers as a standard, state-of-the-art method for estimating IHB energy as well as those of other noncovalent interactions.

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.

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.

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.

Absorptive Desulfurization of Model Biogas Stream Using Choline Chloride-Based Deep Eutectic Solvents

The paper presents a synthesis of deep eutectic solvents (DESs) based on choline chloride (ChCl) as hydrogen bond acceptor and phenol (Ph), glycol ethylene (EG), and levulinic acid (Lev) as hydrogen bond donors in 1:2 molar ratio. DESs were successfully used as absorption solvents for removal of dimethyl disulfide (DMDS) from model biogas steam. Several parameters affecting the absorption capacity and absorption rate have been optimized including kinds of DES, temperature, the volume of absorbent, model biogas flow rate, and initial concentration of DMDS. Furthermore, reusability and regeneration of DESs by means of adsorption and nitrogen barbotage followed by the mechanism of absorptive desulfurization by means of density functional theory (DFT) as well as FT-IR analysis were investigated. Experimental results indicate that the most promising DES for biogas purification is ChCl:Ph, due to high absorption capacity, relatively long absorption rate, and easy regeneration. The research on the absorption mechanism revealed that van der Waal interaction is the main driving force for DMDS removal from model biogas.

Dissecting C-H∙∙∙π and N-H∙∙∙π Interactions in Two Proteins Using a Combined Experimental and Computational Approach

C−H∙∙∙π and N−H∙∙∙π interactions can have an important contribution for protein stability. However, direct measurements of these interactions in proteins are rarely reported. In this work, we combined the mutant cycle experiments and molecular dynamics (MD) simulations to characterize C−H∙∙∙π and N−H∙∙∙π interactions and their cooperativity in two model proteins. It is shown that the average C−H∙∙∙π interaction per residue pair is ~ −0.5 kcal/mol while the N−H∙∙∙π interaction is slightly stronger. The triple mutant box measurement indicates that N−H∙∙∙π∙∙∙C−H∙∙∙π and C−H∙∙∙π∙∙∙C−H∙∙∙π can have a positive or negative cooperativity. MD simulations suggest that the cooperativity, depending on the local environment of the interactions, mainly arises from the geometric rearrangement when the nearby interaction is perturbed.

Tyrosine, cysteine, and proton coupled electron transfer in a ribonucleotide reductase-inspired beta hairpin maquette

Tyrosine residues act as intermediates in proton coupled electron transfer reactions (PCET) in proteins. For example, in ribonucleotide reductase (RNR), a tyrosyl radical oxidizes an active site cysteine via a 35 Å pathway that contains multiple aromatic groups. When singlet tyrosine is oxidized, the radical becomes a strong acid, and proton transfer reactions, which are coupled with the redox reaction, may be used to control reaction rate. Here, we characterize a tyrosine-containing beta hairpin, Peptide O, which has a cross-strand, noncovalent interaction between its single tyrosine, Y5, and a cysteine (C14). Circular dichroism provides evidence for a thermostable beta-turn. EPR spectroscopy shows that Peptide O forms a neutral tyrosyl radical after UV photolysis at 160 K. Molecular dynamics simulations support a phenolic/SH interaction in the tyrosine singlet and radical states. Differential pulse voltammetry exhibits pH dependence consistent with the formation of a neutral tyrosyl radical and a pK_a change in two other residues. A redox-coupled decrease in cysteine pK_a from 9 (singlet) to 6.9 (radical) is assigned. At pD 11, picosecond transient absorption spectroscopy after UV photolysis monitors tyrosyl radical recombination via electron transfer (ET). The ET rate in Peptide O is indistinguishable from the ET rates observed in peptides containing a histidine and a cyclohexylalanine (Cha) at position 14. However, at pD 9, the tyrosyl radical decays via PCET, and the decay rate is slowed, when compared to the histidine 14 variant. Notably, the decay rate is accelerated, when compared to the Cha 14 variant. We conclude that redox coupling between tyrosine and cysteine can act as a PCET control mechanism in proteins.

The interplay and strength of the π⋯HF, C⋯HF, F⋯HF and F⋯HC hydrogen bonds upon the formation of multimolecular complexes based on C2H2⋯HF and C2H4⋯HF small dimers

The conception of this theoretical research was idealized aiming to unveil the intermolecular structures of complexes formed by acetylene or ethylene and hydrofluoric acid. At light of computational calculations by using the B3LYP/6-311++G(d,p) method, the geometries of the C₂H₂⋯(HF), C₂H₂⋯2(HF), C₂H₂⋯4(HF), C₂H₄⋯(HF), C₂H₄⋯2(HF) and C₂H₄⋯4(HF) hydrogen-bonded complexes were fully optimized. Moreover, the Post-Hartree-Fock calculations MP2/6-311++G(d,p), MP2/aug-cc-pVTZ, MP4(SDQ)/6-311++G(d,p) and CCSD/6-311++G(d,p) also were also used. The infrared spectra were analyzed in order to identify the new vibrational modes and frequencies of the proton donors shifted to red region. Through the modeling of charge-fluxes on the basis of the Quantum Theory of Atoms In Molecules (QTAIM) and, by contradicting the expectation of the hydrofluorination mechanisms of acetylene or ethylene, C⋯HF was recognized as a new type of hydrogen bond instead of the already well known π⋯H. The calculations of the Natural Bonding Orbital (NBO) and Charges derived from the Electrostatic Potential Grid-based (ChElPG) were also applied to interpret the shifting frequencies as well as measuring of the punctual charge-transfer after the formation of the complexes. Finally, the determination of the stabilization energy was carried out through the arguments of the Fock matrix in NBO basis and through the supermolecule approach. Also it is worthwhile to notice that some algebraic formulations were used for determining the electronic cooperative effect (CE).

Tipping the Balance between S-π and O-π Interactions

A comprehensive experimental survey consisting of 36 molecular balances was conducted to compare 18 pairs of S-π versus O-π interactions over a wide range of structural, geometric, and solvent parameters. A strong linear correlation was observed between the folding energies of the sulfur and oxygen balances across the entire library of balance pairs. The more stable interaction systematically switched from the O-π to S-π interaction. Computational studies of bimolecular PhSCH₃-arene and PhOCH₃-arene complexes were able to replicate the experimental trends in the molecular balances. The change in preference for the O-π to S-π interaction was due to the interplay of stabilizing (dispersion and solvophobic) and destabilizing (exchange-repulsion) terms arising from the differences in size and polarizability of the oxygen and sulfur atoms.

Thiols as Hydrogen Bond Acceptors and Donors: Spectroscopy of 2-Phenylethanethiol Complexes

Evidence and understanding of sulfur-centered hydrogen bonding, especially where the donor is a thiol, lags far behind that for conventional OH interactions. To help address this deficiency, conformer specific IR spectra of 2-phenylethanethiol (PET) and associated 1:1 solvent complexes have been measured in SH, OH, and CH stretch regions using resonant-two-photon-ionization (R2PI) and IR-UV ion dip spectroscopic techniques. The aromatic and aliphatic CH stretch regions show signature differences between anti and gauche conformers. Supported by ab initio calculations, a PET-water cluster with an OH···S arrangement and a PET-diethyl ether cluster expressing an SH···O interaction were identified. The SH stretch band of the SH···O complex is red-shifted and undergoes significant intensity enhancement compared to the bare molecule, which is characteristic of hydrogen bonding. These findings offer insight into the nature of the thiol functional group as a potential hydrogen bond donor and acceptor.

Modeling Protein S-Aromatic Motifs Reveals Their Structural and Redox Flexibility

S-aromatic motifs are important noncovalent forces for protein stability and function but remain poorly understood. Hence, we performed quantum calculations at the MP2(full)/6-311++G(d,p) level on complexes between Cys (H₂S, MeSH) and Met (Me₂S) models with models of Phe (benzene, toluene), Trp (indole, 3-methylindole), Tyr (phenol, 4-methylphenol), and His (imidazole, 4-methylimidazole). The most stable gas-phase conformers exhibit binding energies of -2 to -6 kcal/mol, and the S atom lies perpendicular to the ring plane. This reveals preferential interaction with the ring π-system, except in the imidazoles where S binds edge-on to an N atom. Complexation tunes the gas-phase vertical ionization potentials of the ligands over as much as 1 eV, and strong σ- or π-type H-bonding supports charge transfer to the H-bond donor, rendering it more oxidizable. When the S atom acts as an H-bond acceptor (N/O-H_ar···S), calibration of the CHARMM36 force field (by optimizing pair-specific Lennard-Jones parameters) is required. Implementing the optimized parameters in molecular dynamics simulations in bulk water, we find stable S-aromatic complexes with binding free energies of -0.6 to -1.1 kcal/mol at ligand separations up to 8 Å. The aqueous S-aromatics exhibit flexible binding conformations, but edge-on conformers are less stable in water. Reflecting this, only 0.3 to 10% of the S-indole, S-phenol, and S-imidazole structures are stabilized by N/O-H_ar···S or S-H···O_ar/N_ar σ-type H-bonding. The wide range of energies and geometries found for S-aromatic interactions and their tunable redox properties expose the versatility and variability of the S-aromatic motif in proteins and allow us to predict a number of their reported properties.

Unexpected solvent effects on the UV/Vis absorption spectra of o-cresol in toluene and benzene: in contrast with non-aromatic solvents

Cresol is a prototype molecule in understanding intermolecular interactions in material and biological systems, because it offers different binding sites with various solvents and protonation states under different pH values. It is found that the UV/Vis absorption spectra of o-cresol in aromatic solvents (benzene, toluene) are characterized by a sharp peak, unlike the broad double-peaks in 11 non-aromatic solvents. Both molecular dynamics simulations and electronic structure calculations revealed the formation of intermolecular π-complexation between o-cresol and aromatic solvents. The thermal movements of solvent and solute molecules render the conformations of o-cresol changing between trans and cis isomers. The π-interaction makes the cis configuration a dominant isomer, hence leading to the single keen-edged UV/Vis absorption peak at approximately 283 nm. The free conformation changes between trans and cis in aqueous solution rationalize the broader absorption peaks in the range of 260-280 nm. The pH dependence of the UV/Vis absorption spectra in aqueous solutions is also rationalized by different protonation states of o-cresol. The explicit solvent model with long-ranged interactions is vital to describe the effects of π-complexation and electrostatic interaction on the UV/Vis absorption spectra of o-cresol in toluene and alkaline aqueous (pH > 10.3) solutions, respectively.

Dithiocarboxylic Acids: An Old Theme Revisited and Augmented by New Preparative, Spectroscopic and Structural Facts

The unstable dithiocarboxylic acids dithioacetic acid, 2-methyl-dithiopropionic acid, 2,2-dimethyl-dithiopropionic acid and dithiobenzoic acid were synthesized and characterized by NMR spectroscopy and GC/MS. The stable dithiocarboxylic acids 2,4,6-trimethyl benzoic acid, 2,4,6-tri-iso-propylbenzoic acid and 2,6-dimesityl benzoic acid were synthesized, isolated and characterized by spectroscopic methods and in parts by mass spectrometry and X-ray crystallography. The new data were used to re-evaluate literature data on the synthesis, spectroscopy and structural data of dithiocarboxylic acids as a fundamental class of organic compounds in general.

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.

Exploring the rare S—H...S hydrogen bond using charge density analysis in isomers of mercaptobenzoic acid

Experimental and theoretical charge density analyses on isomers of mercaptobenzoic acid have been carried out to quantify the hydrogen bonding of the hitherto less explored thiols, to assess the strength of the interactions using the topological features of the electron density. The electron density study offers interesting insights into the nature of the S—H...S interaction. The interaction energy is comparable with that of a weak hydrogen bond. The strength and directionality of the S—H...S hydrogen bond is demonstrated to be mainly due to the conformation locking potential of the intramolecular S...O chalcogen bond in 2-mercaptobenzoic acid and is stronger than in 3-mercaptobenzoic acid, which lacks the intramolecular S...O bond. Thepara-substituted mercaptobenzoic acid depicts a type I S...S interaction.

Insights into Thiol-Aromatic Interactions: A Stereoelectronic Basis for S-H/π Interactions

Thiols can engage favorably with aromatic rings in S-H/π interactions, within abiological systems and within proteins. However, the underlying bases for S-H/π interactions are not well understood. The crystal structure of Boc-l-4-thiolphenylalanine tert-butyl ester revealed crystal organization centered on the interaction of the thiol S-H with the aromatic ring of an adjacent molecule, with a through-space Hthiol···Caromatic distance of 2.71 Å, below the 2.90 Å sum of the van der Waals radii of H and C. The nature of this interaction was further examined by DFT calculations, IR spectroscopy, solid-state NMR spectroscopy, and analysis of the Cambridge Structural Database. The S-H/π interaction was found to be driven significantly by favorable molecular orbital interactions, between an aromatic π donor orbital and the S-H σ* acceptor orbital (a π → σ* interaction). For comparison, a structural analysis of O-H/π interactions and of cation/π interactions of alkali metal cations with aromatic rings was conducted. Na+ and K+ exhibit a significant preference for the centroid of the aromatic ring and distances near the sum of the van der Waals and ionic radii, as expected for predominantly electrostatic interactions. Li+ deviates substantially from Na+ and K+. The S-H/π interaction differs from classical cation/π interactions by the preferential alignment of the S-H σ* toward the ring carbons and an aromatic π orbital rather than toward the aromatic centroid. These results describe a potentially broadly applicable approach to understanding the interactions of weakly polar bonds with π systems.