Competition between O-H and S-H Intermolecular Interactions in Conformationally Complex Systems: The 2-Phenylethanethiol and 2-Phenylethanol Dimers

Noncovalent interactions involving sulfur centers play a relevant role in biological and chemical environments. Yet, detailed molecular descriptions are scarce and limited to very simple model systems. Here we explore the formation of the elusive S-H···S hydrogen bond and the competition between S-H···O and O-H···S interactions in pure and mixed dimers of the conformationally flexible molecules 2-phenylethanethiol (PET) and 2-phenylethanol (PEAL), using the isolated and size-controlled environment of a jet expansion. The structure of both PET-PET and PET-PEAL dimers was unraveled through a comprehensive methodology that combined rotationally resolved microwave spectroscopy, mass-resolved isomer-specific infrared laser spectroscopy, and quantum chemical calculations. This synergic experimental-computational approach offered unique insights into the potential energy surface, conformational equilibria, molecular structure, and intermolecular interactions of the dimers. The results show a preferential order for establishing hydrogen bonds following the sequence S-H···S < S-H···O ≲ O-H···S < O-H···O, despite the hydrogen bond only accounting for a fraction of the total interaction energy.

Exploring the conformational landscape, hydrogen bonding, and internal dynamics in the diallyl ether and diallyl sulfide monohydrates

The conformational spaces of the diallyl ether (DAE) and diallyl sulfide (DAS) monohydrates were explored using rotational spectroscopy from 6 to 19 GHz. Calculations at the B3LYP-D3(BJ)/aug-cc-pVTZ level suggested significant differences in their conformational behavior, with DAE-w exhibiting 22 unique conformers and DAS-w featuring three stable structures within 6 kJ mol-1. However, only transitions from the lowest energy conformer of each were experimentally observed. Spectral analysis confirmed that binding with water does not alter the conformational preference for the lowest energy structure of the monomers, but it does influence the relative stabilities of all other conformers, particularly in the case of DAE. Non-covalent interaction and quantum theory of atoms in molecules analyses showed that the observed conformer for each complex is stabilized by two intermolecular hydrogen bonds (HBs), where water primarily interacts with the central oxygen or sulfur atom of the diallyl compounds, along with secondary interactions involving the allyl groups. The nature of these interactions was further elucidated using symmetry-adapted perturbation theory, which suggests that the primary HB interaction with S in DAS is weaker and more dispersive in nature compared to the primary HB in DAE. This supports the experimental observation of a tunneling splitting exclusively in the rotational spectrum of DAS-w, as the weaker contact allows water to undergo internal motions within the complex, as shown based on calculated transition state structures for possible tunneling pathways.

Noncovalent interactions of aromatic heterocycles: rotational spectroscopy and theoretical calculations of the thiazole-CF4 and thiazole-SF6 complexes

The complexes of thiazole with CF4 and SF6 have been investigated by Fourier transform microwave spectroscopy and quantum chemical calculations. One rotational spectrum was observed for the thiazole-CF4 complex. Experiments and theoretical computations confirmed that the observed structure of thiazole-CF4 is primarily formed due to N⋯CCF4 interaction with the C atom of CF4 located in the plane of the thiazole ring. The rotational transitions of thiazole-CF4 exhibit A/E torsional splitting induced by the internal rotation of the -CF3 top. The potential barrier of the -CF3 internal rotation is 0.2411(1) kJ mol-1, consistent with the calculated value (∼0.3 kJ mol-1). For the thiazole-SF6 complex, one conformer with SF6 located above the thiazole ring is detected. The observed structure of thiazole-SF6 is mainly stabilized by van der Waals interactions. The energy decomposition analysis reveals that the electrostatics and dispersion are the dominant attractive contributions to the formation of thiazole-CF4 and thiazole-SF6 dimers, whereas the weight of the dispersion term becomes more significant in the thiazole-SF6 complex compared to that of the thiazole-CF4 complex.

Nature and Strength of Sulfur-Centered Hydrogen Bond in Methanethiol Aqueous Solutions

Methanethiol (M) and water (W) clusters like dimers (M₁W₁, M₂, and W₂), trimers (M₁W₂, M₂W₁, M₃, and W₃), and tetramers (M₁W₃, M₂W₂, M₃W₁, M₄, and W₄) were studied to assess the strength of sulfur-centered hydrogen bonding using different levels of theories, viz, HF, MP2, MP3, MP4, B3LYP, B3LYP-D3, CCSD, CCSD(T)-F12, and CCSD(T) along with aug-cc-pVNZ (where N = D, T, and Q) basis sets. Interaction energies were found to be in the range of -3.3 to -5.3 kcal/mol for the dimers, -8.0 to -16.7 kcal/mol for the trimers, and -13.5 to -29.5 kcal/mol for the tetramers at the B3LYP-D3/CBS limit level of theory. Normal modes of vibrations computed at the B3LYP/cc-pVDZ level of theory were seen to be in good agreement with the experimental values. Local energy decomposition calculations using the DLPNO-CCSD(T) level of theory indicated the domination of electrostatic interactions' contribution to the interaction energy in all cluster systems. Furthermore, atoms in molecules and natural bond orbital calculations both carried out at the B3LYP-D3/aug-cc-pVQZ level of theory aided in visualizing the hydrogen bonds besides proving a rationale for the strength of the hydrogen bonds and thereby the stability of these cluster systems.

Hydrogen bond interactions between thioethers and amides: A joint rotational spectroscopic and theoretical study of the formamide⋯dimethyl sulfide adduct

The rotational spectrum of the binary adduct of formamide (HCONH₂) with dimethyl sulfide (DMS) has been investigated employing cavity-based Fourier transform microwave spectroscopy combined with theoretical computations. Experimentally, only one isomer of the adduct was unambiguously observed and assigned according to the theoretically predicted spectroscopic parameters, and its rotational spectrum displays the hyperfine splittings associated with the ¹⁴N nuclear quadrupole coupling effect. The observed isomer exhibits C_s symmetry, such that the ∠CSC angle of the DMS subunit is bisected by the ab-plane of the HCONH₂ moiety. The two moieties in the detected isomer are connected via one primary NH···S and two secondary CH···O hydrogen bonds. Quantum theory of atoms in molecules (QTAIM), non-covalent interaction (NCI), natural bond orbital (NBO) and symmetry-adapted perturbation theory (SAPT) approaches were utilized for characterizing the intermolecular interactions occurring in the titled adduct. Additionally, the adduct of HCONH₂ with dimethyl ether (DME) was also theoretically investigated to compare the difference in structure and energy characteristics between the NH···S and NH···O hydrogen bonds.

Hydrogen Bonding in the Dimer and Monohydrate of 2-Adamantanol: A Test Case for Dispersion-Corrected Density Functional Methods

Weakly-bound intermolecular clusters constitute reductionist physical models for non-covalent interactions. Here we report the observation of the monomer, the dimer and the monohydrate of 2-adamantanol, a secondary alcohol with a bulky ten-carbon aliphatic skeleton. The molecular species were generated in a supersonic jet expansion and characterized using broadband chirped-pulse microwave spectroscopy in the 2–8 GHz frequency region. Two different gauche-gauche O-H···O hydrogen-bonded isomers were observed for the dimer of 2-adamantanol, while a single isomer was observed for the monomer and the monohydrate. The experimental rotational parameters were compared with molecular orbital calculations using density functional theory (B3LYP-D3(BJ), B2PLYP-D3(BJ), CAM-B3LYP-D3(BJ), ωB97XD), additionally providing energetic and electron density characterization. The shallow potential energy surface makes the dimer an interesting case study to benchmark dispersion-corrected computational methods and conformational search procedures.

Torsional chirality and molecular recognition: the homo and heterochiral dimers of thenyl and furfuryl alcohol

Furfuryl alcohol and thenyl alcohol contain a labile torsional chiral center, producing transiently chiral enantiomers interconverting in the nanosecond time-scale. We explored chiral molecular recognition using the weakly-bound intermolecular dimers of both alcohols, freezing stereomutation. Supersonic jet broadband microwave spectroscopy revealed homo and heterochiral diastereoisomers for each alcohol dimer and the structural characteristics of the clusters. All dimers are primarily stabilized by a moderately intense O-H⋯O hydrogen bond, but differ in the secondary interactions, which introduce additional hydrogen bonds either to the ring oxygen in furfuryl alcohol or to the π ring system in thenyl alcohol. Density-functional calculations (B2PLYP-D3(BJ)/def2-TZVP) show no clear preferences for a particular stereochemistry in the dimers, with relative energies of the order 1-2 kJ mol^-1. The study suggests opportunities for the investigation of chiral recognition in molecules with torsional barriers in between transient and permanent interconversion regimes.

t BuOLi-promoted terminal alkyne functionalizations by aliphatic thiols and alcohols

Selective radical addition to terminal alkynes is always a difficult task to achieve because it gives a mixture of stereo- and regioisomers. Herein we describe the selective addition of aliphatic thiols or alcohols to N-phenylpropiolamides (terminal alkynes) using lithium tert-butoxide (^tBuOLi) in ethanol as a promoter. Mechanistically, it has been shown that the reaction proceeded through the generation of a thiyl radical intermediate, and the amide group in N-phenylpropiolamide could help in the activation of the alkyne, which led to thioacetalization via the formation of a (Z)-selective anti-Markovnikov vinyl sulfide. The (Z)-selectivity during the formation of vinyl sulfides was controlled by an intramolecular sulfur⋯oxygen interaction.

C-H⋯S hydrogen bonding interactions

The short C-H⋯S contacts found in available structural data for both small molecules and larger biomolecular systems suggest that such contacts are an often overlooked yet important stabilizing interaction. Moreover, many of these short C-H⋯S contacts meet the definition of a hydrogen bonding interaction. Using available structural data from the Cambridge Structural Database (CSD), as well as selected examples from the literature in which important C-H⋯S contacts may have been overlooked, we highlight the generality of C-H⋯S hydrogen bonding as an important stabilizing interaction. To uncover and establish the generality of these interactions, we compare C-H⋯S contacts with other traditional hydrogen bond donors and acceptors as well as investigate how coordination number and metal bonding affect the preferred geometry of interactions in the solid state. This work establishes that the C-H⋯S bond meets the definition of a hydrogen bond and serves as a guide to identify C-H⋯S hydrogen bonds in diverse systems.

Chirality-Puckering correlation and intermolecular interactions in Sphingosines: Rotational spectroscopy of jaspine B3 and its monohydrate

Chirality is determinant for sphingosine biofunctions and pharmacological activity, yet the reasons for the biological chiral selection are not well understood. Here, we characterized the intra- and intermolecular interactions at the headgroup of the cytotoxic anhydrophytosphingosine jaspine B, revealing chirality-dependent correlations between the puckering of the ring core and the formation of amino-alcohol hydrogen bond networks, both in the monomer and the monohydrate. Following the specific synthesis of a shortened 3-carbon side-chain molecule, denoted jaspine B3, six different isomers were observed in a jet expansion using broadband (chirped-pulsed) rotational spectroscopy. Additionally, a single isomer of the jaspine B3 monohydrate was observed, revealing the insertion of water in between the hydroxy and amino groups and the formation of a network of O-H···N-H···O_ring hydrogen bonds. The specific jaspine B3 stereochemistry thus creates a double-faced molecule where the exposed lone-pair electrons may easily catalyze the formation of intermolecular aggregates and determine the sphingosine biological properties.

Setting up the HyDRA blind challenge for the microhydration of organic molecules

The procedure leading to the first HyDRA blind challenge for the prediction of water donor stretching vibrations in monohydrates of organic molecules is described. A training set of 10 monohydrates...

Molecular Recognition, Transient Chirality and Sulfur Hydrogen Bonding in the Benzyl Mercaptan Dimer

The homodimers of transiently chiral molecules offer physical insight into the process of molecular recognition, the preference for homo or heterochiral aggregation and the nature of the non-covalent interactions stabilizing the adducts. We report the observation of the benzyl mercaptan dimer in the isolation conditions of a supersonic jet expansion, using broadband (chirped-pulse) microwave spectroscopy. A single homochiral isomer was observed for the dimer, stabilized by a cooperative sequence of S-H···S and S-H···π hydrogen bonds. The structural data, stabilization energies and energy decomposition describe these non-covalent interactions as weak and dispersion-controlled. A comparison is also provided with the benzyl alcohol dimer.

Hydrogen bonding networks and cooperativity effects in the aqueous solvation of trimethylene oxide and sulfide rings by microwave spectroscopy and computational chemistry

The intermolecular interactions responsible for the microsolvation of the highly flexible trimethylene oxide (TMO) and trimethylene sulfide (TMS) rings with one and two water (w) molecules were investigated using rotational spectroscopy (8-22 GHz) and quantum chemical calculations. The observed patterns of transitions are consistent with the most stable geometries of the TMO-w, TMO-(w)₂, and TMS-w complexes at the B2PLYP-D3(BJ)/aug-cc-pVTZ level and were confirmed using spectra of the ¹⁸O isotopologue. Due to its effectively planar backbone, TMO offers one unique binding site for solvation, while water can bind to the puckered TMS ring in either an axial or equatorial site of the heteroatom. In all clusters, the first water molecule binds in the σ_v symmetry plane of the ring monomer and serves as a hydrogen bond donor to the heteroatom. The second water molecule is predicted to form a cooperative hydrogen bonding network between the three moieties. Secondary C-H⋯O interactions are a key stabilizing influence in trimers and also drive the preferred binding site in the TMS clusters with the axial binding site preferred in TMS-w and the equatorial form calculated to be more stable in the dihydrate. Using an energy partition scheme from the symmetry-adapted perturbation theory for the O, S, and Se containing mono- and dihydrates, the intermolecular interactions are revealed to be mainly electrostatic, but the dispersive character of the contacts is enhanced with the increasing size of the ring's heteroatom due to the key role of longer-range secondary interactions.

Sulfuroxygen interaction-controlled (Z)-selective anti-Markovnikov vinyl sulfides

The sulfur oxygen (SO) interaction was used herein to obtain (Z)-selective anti-Markovnikov vinyl sulfides from the addition of thiyl radicals to terminal alkynes. DFT calculations predicted that SO interaction originated from the delocalization of the lone-pair of the carbonyl oxygen to the adjacent σ* orbital of the S atom of C-S.

Switching Hydrogen Bonding to π-Stacking: The Thiophenol Dimer and Trimer

We used jet-cooled broadband rotational spectroscopy to explore the balance between π-stacking and hydrogen-bonding interactions in the self-aggregation of thiophenol. Two different isomers were detected for the thiophenol dimer, revealing dispersion-controlled π-stacked structures anchored by a long S-H···S sulfur hydrogen bond. The weak intermolecular forces allow for noticeable internal dynamics in the dimers, as tunneling splittings are observed for the global minimum. The large-amplitude motion is ascribed to a concerted inversion motion between the two rings, exchanging the roles of the proton donor and acceptor in the thiol groups. The determined torsional barrier of B₂ = 250.3 cm^-1 is consistent with theoretical predictions (290-502 cm^-1) and the monomer barrier of 277.1(3) cm^-1. For the thiophenol trimer, a symmetric top structure was assigned in the spectrum. The results highlight the relevance of substituent effects to modulate π-stacking geometries and the role of the sulfur-centered hydrogen bonds.

NMR quantification of H-bond donating ability for bioactive functional groups and isosteres

The H-bond donating ability for 127 compounds including drug fragments and isosteres have been quantified using a simple and rapid method with ³¹P NMR spectroscopy. Functional groups important to medicinal chemistry were evaluated including carboxylic acids, alcohols, phenols, thioic acids and nitrogen group H-bond donors. ³¹P NMR shifts for binding to a phosphine oxide probe have a higher correlation with equilibrium constants for H-bonding (log K^H_A) than acidity (pK_a), indicating that these binding experiments are representative of H-bonding ability and not proton transfer. Additionally, ³¹P NMR binding data for carboxylic acid isosteres correlates with physicochemical properties such as lipophilicity, membrane permeability and plasma protein binding. This method has been used to evaluate the H-bond donating ability of small molecule drug compounds such as NSAIDs and antimicrobials.