Conformational Preferences Research Articles

Interaction of 3,9-disubstituted acridine with single stranded poly(rA), double stranded poly(rAU) and triple stranded poly(rUAU): molecular docking – A spectroscopic tandem study

RNA plays an important role in many biological processes which are crucial for cell survival, and it has been suggested that it may be possible to inhibit individual processes involved in many diseases by targeting specific sequences of RNA. The aim of this work is to determine the affinity of novel 3,9-disubstited acridine derivative 1 with three different RNA molecules, namely single stranded poly(rA), double stranded homopolymer poly(rAU) and triple stranded poly(rUAU). The results of the absorption titration assays show that the binding constant of the novel derivative to the RNA molecules was in the range of 1.7–6.2 × 104 mol dm−3. The fluorescence and circular dichroism titration assays revealed considerable changes. The most significant results in terms of interpreting the nature of the interactions were the melting temperatures of the RNA samples in complexes with the 1. In the case of poly(rA), denaturation resulted in a self-structure formation; increased stabilization was observed for poly(rAU), while the melting points of the ligand-poly(rUAU) complex showed significant destabilization as a result of the interaction. The principles of molecular mechanics were applied to propose the non-bonded interactions within the binding complex, pentariboadenylic acid and acridine ligand as the study model. Initial molecular docking provided the input structure for advanced simulation techniques. Molecular dynamics simulation and cluster analysis reveal π – π stacking and the hydrogen bonds formation as the main forces that can stabilize the binding complex. Subsequent MM-GBSA calculations showed negative binding enthalpy accompanied the complex formation and proposed the most preferred conformation of the interaction complex.

Delineating the impact of N21D mutation on the conformational preferences and structural transitions in human islet amyloid polypeptide

The amyloid deposits of human islet amyloid polypeptide (hIAPP) in the pancreatic β-cells are linked to a chronic metabolic disease known as type 2 diabetes (T2D). Recent studies reported that N21D mutation in hIAPP inhibited the formation of fibrillar aggregates, however, the underlying molecular mechanism of fibril inhibition in the N21D mutant remains largely unclear. Thus, molecular dynamics (MD) simulations have been employed to illuminate the effect of N21D mutation on the conformational dynamics and structural transformations in hIAPP at the atomistic level. Notably, the structural stability of WT-hIAPP (wild type-hIAPP) was enhanced upon N21D mutation. Interestingly, the sampling of helical conformation in hIAPP was significantly increased from 12.60 to 19.80% with a concomitant decrease in β-sheet from 7.40 to 1.20% upon N21D mutation that corroborates with the circular dichroism (CD). The dihedral principal component analysis (dPCA) depicted a notable reduction in the conformational phase space of [N21D]hIAPP as compared to WT-hIAPP, which highlights the prevention of conformational transitions and enhanced structural stability. Remarkably, a complete loss in the β-sheet content was noted in the conformations extracted from the free energy landscape of [N21D]hIAPP as compared to 24–25% sampling of β-sheets in minimum energy conformations of WT-hIAPP consistent with the transmission electron microscopy (TEM) and thioflavin T (ThT) fluorescence assay. The present study provided key insights into the N21D mutation-induced conformational changes in hIAPP and highlighted that Asn21 plays a key role in hIAPP aggregation. The insights from this work will be beneficial in the design of drug candidates by selectively targeting Asn21 of hIAPP which could interfere with its aggregation and associated cytotoxicity.

Dissecting the mechanism of atlastin-mediated homotypic membrane fusion at the single-molecule level

Homotypic membrane fusion of the endoplasmic reticulum (ER) is mediated by dynamin-like GTPase atlastin (ATL). This fundamental process relies on GTP-dependent domain rearrangements in the N-terminal region of ATL (ATLcyto), including the GTPase domain and three-helix bundle (3HB). However, its conformational dynamics during the GTPase cycle remain elusive. Here, we combine single-molecule FRET imaging and molecular dynamics simulations to address this conundrum. Different from the prevailing model, ATLcyto can form a loose crossover dimer upon GTP binding, which is tightened by GTP hydrolysis for membrane fusion. Furthermore, the α-helical motif between the 3HB and transmembrane domain, which is embedded in the surface of the lipid bilayer and self-associates in the crossover dimer, is required for ATL function. To recycle the proteins, Pi release, which disassembles the dimer, activates frequent relative movements between the GTPase domain and 3HB, and subsequent GDP dissociation alters the conformational preference of the ATLcyto monomer for entering the next reaction cycle. Finally, we found that two disease-causing mutations affect human ATL1 activity by destabilizing GTP binding-induced loose crossover dimer formation and the membrane-embedded helix, respectively. These results provide insights into ATL-mediated homotypic membrane fusion and the pathological mechanisms of related disease.

Classifying protein kinase conformations with machine learning.

Protein kinases are key actors of signaling networks and important drug targets. They cycle between active and inactive conformations, distinguished by a few elements within the catalytic domain. One is the activation loop, whose conserved DFG motif can occupy DFG-in, DFG-out, and some rarer conformations. Annotation and classification of the structural kinome are important, as different conformations can be targeted by different inhibitors and activators. Valuable resources exist; however, large-scale applications will benefit from increased automation and interpretability of structural annotation. Interpretable machine learning models are described for this purpose, based on ensembles of decision trees. To train them, a set of catalytic domain sequences and structures was collected, somewhat larger and more diverse than existing resources. The structures were clustered based on the DFG conformation and manually annotated. They were then used as training input. Two main models were constructed, which distinguished active/inactive and in/out/other DFG conformations. They considered initially 1692 structural variables, spanning the whole catalytic domain, then identified ("learned") a small subset that sufficed for accurate classification. The first model correctly labeled all but 3 of 3289 structures as active or inactive, while the second assigned the correct DFG label to all but 17 of 8826 structures. The most potent classifying variables were all related to well-known structural elements in or near the activation loop and their ranking gives insights into the conformational preferences. The models were used to automatically annotate 3850 kinase structures predicted recently with the Alphafold2 tool, showing that Alphafold2 reproduced the active/inactive but not the DFG-in proportions seen in the Protein Data Bank. We expect the models will be useful for understanding and engineering kinases.

Deciphering the conformational preference and ionization dynamics of tetrahydrofuran: Insights from advanced spectroscopic techniques.

Tetrahydrofuran (THF) has garnered significant attention due to its pivotal role in biological and chemical processes. The diverse array of conformations exhibited by THF profoundly impacts its reactivity and interactions with other molecules. Understanding these conformational preferences is crucial for comprehending its molecular behavior. In this study, we utilize infrared (IR) resonant vacuum ultraviolet photoionization/mass-analyzed threshold ionization (VUV-PI/MATI) mass spectroscopies to capture distinctive vibrational spectra of individual conformers, namely, "twisted" and "bent," within THF. Our conformer-specific vibrational spectra provide valuable insights into the relative populations of these two conformers. The analysis reveals that the twisted (C2) conformer is more stable than the bent (CS) conformer by 17 ± 15cm-1. By precisely tuning the VUV photon energy to coincide with vibrational excitation via IR absorption, we selectively ionize specific conformers, yielding two-photon IR + VUV-PI/MATI spectra corresponding to the twisted and bent conformers. This investigation conclusively affirms that both the twisted and bent conformers coexist in the neutral state, while only the twisted conformer exists in the cationic state. These findings not only bridge gaps in existing knowledge but also provide profound insights into the behavior of this pivotal molecule in the realms of biology and medicine.

Exploring helix structures of γ-peptides based on 2-(aminomethyl)cyclopentanecarboxylic acid.

Conformational search and density functional theory calculations were performed to explore the preferences of helical structures for chiro-specific oligo-γ-peptides of 2-(aminomethyl)cyclopentanecarboxylic acid (γAmc5) with a cyclopentyl constraint on the Cα-Cβ bond in solution. The dimer and tetramer of γAmc5 (1) with homochiral (1S, 2S) configurations exhibited a strong preference for the 9-membered helix foldamer in solution, except for the tetramer in water. However, the oligomers of γAmc5 (1) longer than tetramer preferentially adopted a right-handed (P)-2.614-helix (H1-14) as the peptide sequence becomes longer and as solvent polarity increases. The high stabilities for H1-14 foldamers of γAmc5 (1) in solution were ascribed to the favored solvation free energies. The calculated mean backbone torsion angles for H1-14 helix foldamers of γAmc5 (1) were similar to those calculated for oligomers of other γ-residues with cyclopentane or cyclohexane rings. However, the substitution of cyclopentane constraints on the Cα-Cβ bond of the γAmc5 (1) residue resulted in different conformational preferences and/or handedness of helix foldamers. In particular, the pyrrolidine-substituted analogs of the H1-14 foldamers of γAmc5 (1) with adjacent amine diads substituted at a proximal distance are expected to be potential catalysts for the crossed aldol condensation in nonpolar and polar solvents.

Isocyanide-Based Multicomponent Reactions Coupled One-Pot Process: Efficient Tools to Diversity-Oriented Synthesis of Structural Peptidomimetics.

Multicomponent diversity-oriented synthesis (DOS) of conformationally anchored structural peptidomimetics like 2,5-diketopiperazines (2,5-DKP) containing heterocyclic bioisosteres of the amide bond, such as 1,2,3-triazoles and 1,5-disubstituted tetrazoles (1,5-DS-T) is described. Structural peptidomimetics are synthesized from similar available starting materials, via a strategy based on isocyanide-based multicomponent reactions (IMCRs): Ugi-4CR and Ugi-Azide (UA), followed by a one-pot process: SN2/intramolecular alkyne-azide cycloaddition (IAAC). The sequential aligning of two powerful synthetic tools (IMCR and IAAC) has parallelly contributed to generate anchored conformation and complexity in target molecules, which are considered structural peptidomimetics of 2,5-DKP. Herein, the 1,2,3-triazole ring plays a key role in the preference for the boat conformation. Furthermore, the use of UA reaction generates scaffold diversity at the N-1 α-carbon of the pyrazinone ring, replacing a linear amide bond with a heterocyclic bioisostere such as 1,5-DS-T leading to the synthesis of novel tricyclic peptidomimetics. The DFT calculations confirmed the boat conformation of the synthesized molecules.

Beyond reader proteins: RNA binding proteins and RNA modifications in conversation to regulate gene expression.

Post-transcriptional mRNA modifications play diverse roles in gene expression and RNA function. In many cases, RNA modifications function by altering how cellular machinery such as RNA binding proteins (RBPs) interact with RNA substrates. For instance, N6-methyladenosine (m6A) is recognized by the well-characterized YTH domain-containing family of "reader" proteins. For other mRNA modifications, similar global readers of modification status have not been clearly defined. Rather, most interactions between RBPs and RNA modifications have a more complicated dependence on sequence context and binding modality. The current handful of studies that demonstrate modifications impacting protein binding likely represent only a fraction of the full landscape. In this review, we dissect the known instances of RNA modifications altering RBP binding, specifically m6A, N1-methyladenosine (m1A), 5-methylcytosine (m5C), pseudouridine (Ψ), and internal N7-methylguanosine. We then review the biochemical properties of these and other identified mRNA modifications including dihydrouridine (D), N4-acetylcytosine (ac4C), and 2'-O-Methylation (Nme). We focus on how these properties would be likely to impact RNA:RBP interactions, including by changes to hydrogen bond potential, base-stacking efficiency, and RNA conformational preferences. The effects of RNA modifications on secondary structure have been well-studied, and we briefly discuss how structural effects imparted by modifications can lead to protein binding changes. Finally, we discuss strategies for uncovering as-yet-to-be identified modification-sensitive RBP:RNA Interactions. Coordinating future efforts to intersect the epitranscriptome and the RNA-protein interactome will illuminate the rules governing RNA modification recognition and the mechanisms responsible for the biological consequences of mRNA modification. This article is categorized under: RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry RNA Interactions with Proteins and Other Molecules > Protein-RNA Recognition RNA Processing > RNA Editing and Modification.

Conformational Preference of Lithium Polysulfide Clusters Li2Sx (x = 4-8) in Lithium-Sulfur Batteries.

Structures are of fundamental importance for diverse studies of lithium polysulfide clusters, which govern the performance of lithium-sulfur batteries. The ring-like geometries were regarded as the most stable structures, but their physical origin remains elusive. In this work, we systematically explored the minimal structures of Li2Sx (x = 4-8) clusters to uncover the driving force for their conformational preferences. All low-lying isomers were generated by performing global searches using the ABCluster program, and the ionic nature of the Li···S interactions was evidenced with the energy decomposition analysis based on the block-localized wave function (BLW-ED) approach and further confirmed with the quantum theory of atoms in molecule (QTAIM). By analysis of the contributions of various energy components to the relative stability with the references of the lowest-lying isomers, the controlling factor for isomer preferences was found to be the polarization interaction. Notably, although the electrostatic interaction dominates the binding energies, it contributes favorably to the relative stabilities of most isomers. The Li+···Li+ distance is identified as the key geometrical parameter that correlates with the strength of the polarization of the Sx2- fragment imposed by the Li+ cations. Further BLW-ED analyses reveal that the cooperativity of the Li+ cations primarily determines the relative strength of the polarization.

Does Inter-Residue Hydrogen Bonding in β-(1→4)-Linked Disaccharides Influence Linkage Conformation in Aqueous Solution?

Two disaccharides, methyl β-d-galactopyranosyl-(1→4)-α-d-glucopyranoside (1) and methyl β-d-galactopyranosyl-(1→4)-3-deoxy-α-d-ribo-hexopyranoside (3), were prepared with selective 13C-enrichment to allow measurement of six trans-O-glycosidic J-couplings (2JCOC, 3JCOCH, and 3JCOCC) in each compound. Density functional theory (DFT) was used to parameterize Karplus-like equations that relate these J-couplings to either ϕ or ψ. MA'AT analysis was applied to both linkages to determine mean values of ϕ and ψ in each disaccharide and their associated circular standard deviations (CSDs). Results show that deoxygenation at C3 of 1 has little effect on both the mean values and librational motions of the linkage torsion angles. This finding implies that, if inter-residue hydrogen bonding between O3H and O5' of 1 is present in aqueous solution and persistent, it plays little if any role in dictating preferred linkage conformation. Hydrogen bonding may lower the energy of the preferred linkage geometry but does not determine it to any appreciable extent. Aqueous 1-μs MD simulation supports this conclusion and also indicates greater conformational flexibility in deoxydisaccharide 3 in terms of sampling several, conformationally distinct, higher-energy conformers in solution. The populations of these latter conformers are low (3-14%) and could not be validated by MA'AT analysis. If the MD model is correct, however, C3 deoxygenation does enable conformational sampling over a wider range of ϕ/ψ values, but linkage conformation in the predominant conformer is essentially identical in both 1 and 3.

Intrinsic Lipid Curvature and Bilayer Elasticity as Regulators of Channel Function: A Comparative Single-Molecule Study.

Perturbations in bilayer material properties (thickness, lipid intrinsic curvature and elastic moduli) modulate the free energy difference between different membrane protein conformations, thereby leading to changes in the conformational preferences of bilayer-spanning proteins. To further explore the relative importance of curvature and elasticity in determining the changes in bilayer properties that underlie the modulation of channel function, we investigated how the micelle-forming amphiphiles Triton X-100, reduced Triton X-100 and the HII lipid phase promoter capsaicin modulate the function of alamethicin and gramicidin channels. Whether the amphiphile-induced changes in intrinsic curvature were negative or positive, amphiphile addition increased gramicidin channel appearance rates and lifetimes and stabilized the higher conductance states in alamethicin channels. When the intrinsic curvature was modulated by altering phospholipid head group interactions, however, maneuvers that promote a negative-going curvature stabilized the higher conductance states in alamethicin channels but destabilized gramicidin channels. Using gramicidin channels of different lengths to probe for changes in bilayer elasticity, we found that amphiphile adsorption increases bilayer elasticity, whereas altering head group interactions does not. We draw the following conclusions: first, confirming previous studies, both alamethicin and gramicidin channels are modulated by changes in lipid bilayer material properties, the changes occurring in parallel yet differing dependent on the property that is being changed; second, isolated, negative-going changes in curvature stabilize the higher current levels in alamethicin channels and destabilize gramicidin channels; third, increases in bilayer elasticity stabilize the higher current levels in alamethicin channels and stabilize gramicidin channels; and fourth, the energetic consequences of changes in elasticity tend to dominate over changes in curvature.

From Functional Fatty Acids to Potent and Selective Natural-Product-Inspired Mimetics via Conformational Profiling.

Fatty acids play important signaling roles in biology, albeit typically lacking potency or selectivity, due to their substantial conformational flexibility. While being recognized as having properties of potentially great value as therapeutics, it is often the case that the functionally relevant conformation of the natural fatty acid is not known, thereby complicating efforts to develop natural-product-inspired ligands that have similar functional properties along with enhanced potency and selectivity profiles. In other words, without structural information associated with a particular functional relationship and the hopelessly unbiased conformational preferences of the endogenous ligand, one is molecularly ill-informed regarding the precise ligand-receptor interactions that play a role in driving the biological activity of interest. To address this problem, a molecular strategy to query the relevance of distinct subpopulations of fatty acid conformers has been established through "conformational profiling", a process whereby a unique collection of chiral and conformationally constrained fatty acids is employed to deconvolute beneficial structural features that impart natural-product-inspired function. Using oleic acid as an example because it is known to engage a variety of receptors, including GPR40, GPR120, and TLX, a 24-membered collection of mimetics was designed and synthesized. It was then demonstrated that this collection contained members that have enhanced potency and selectivity profiles, with some being clearly biased for engagement of the GPCRs GPR40 and GPR120 while others were identified as potent and selective modulators of the nuclear receptor TLX. A chemical synthesis strategy that exploited the power of modern technology for stereoselective synthesis was critical to achieving success, establishing a common sequence of bond-forming reactions to access a disparate collection of chiral mimetics, whose conformational preferences are impacted by the nature of stereodefined moieties differentially positioned about the C18 skeleton of the parent fatty acid. Overall, this study establishes a foundation to fuel future programs aimed at developing natural-product-inspired fatty acid mimetics as valuable tools in chemical biology and potential therapeutic leads.

Unraveling Atomic Contributions to the London Dispersion Energy: Insights into Molecular Recognition and Reactivity.

We present a general framework that enables quantification with atomic resolution of the overall London dispersion energy, which can be readily integrated with currently available energy decomposition schemes. This approach can be used to determine the contribution of individual atoms and functional groups to molecular recognition, conformational preferences, molecular stability, and reactivity. Its efficacy across diverse realms of molecular chemistry and biology is demonstrated with application to molecular balances in solution, asymmetric organocatalytic transformations, and a subcomplex of the F1FO ATP synthase.

Deciphering the molecular choreography of Janus kinase 2 inhibition via Gaussian accelerated molecular dynamics simulations: a dynamic odyssey.

The Janus kinases (JAK) are crucial targets in drug development for several diseases. However, accounting for the impact of possible structural rearrangements on the binding of different kinase inhibitors is complicated by the extensive conformational variability of their catalytic kinase domain (KD). The dynamic KD contains mainly four prominent mobile structural motifs: the phosphate-binding loop (P-loop), the αC-helix within the N-lobe, the Asp-Phe-Gly (DFG) motif, and the activation loop (A-loop) within the C-lobe. These distinct structural orientations imply a complex signal transmission path for regulating the A-loop's flexibility and conformational preference for optimal JAK function. Nevertheless, the precise dynamical features of the JAK induced by different types of inhibitors still remain elusive. We performed comparative, microsecond-long, Gaussian accelerated molecular dynamics simulations in triplicate of three phosphorylated JAK2 systems: the KD alone, type-I ATP-competitive inhibitor (CI) bound KD in the catalytically active DFG-in conformation, and the type-II inhibitor (AI) bound KD in the catalytically inactive DFG-out conformation. Our results indicate significant conformational variations observed in the A-loop and αC helix motions upon inhibitor binding. Our studies also reveal that the DFG-out inactive conformation is characterized by the closed A-loop rearrangement, open catalytic cleft of N and C-lobe, the outward movement of the αC helix, and open P-loop states. Moreover, the outward positioning of the αC helix impacts the hallmark salt bridge formation between Lys882 and Glu898 in an inactive conformation. Finally, we compared their ligand binding poses and free energy by the MM/PBSA approach. The free energy calculations suggested that the AI's binding affinity is higher than CI against JAK2 due to an increased favorable contribution from the total non-polar interactions and the involvement of the αC helix. Overall, our study provides the structural and energetic insights crucial for developing more promising type I/II JAK2 inhibitors for treating JAK-related diseases.

ZNL0325, a Pyrazolopyrimidine-Based Covalent Probe, Demonstrates an Alternative Binding Mode for Kinases.

The pyrazolopyrimidine (PP) heterocycle is a versatile and widely deployed core scaffold for the development of kinase inhibitors. Typically, a 4-amino-substituted pyrazolopyrimidine binds in the ATP-binding pocket in a conformation analogous to the 6-aminopurine of ATP. Here, we report the discovery of ZNL0325 which exhibits a flipped binding mode where the C3 position is oriented toward the ribose binding pocket. ZNL0325 and its analogues feature an acrylamide side chain at the C3 position which is capable of forming a covalent bond with multiple kinases that possess a cysteine at the αD-1 position including BTK, EGFR, BLK, and JAK3. These findings suggest that the ability to form a covalent bond can override the preferred noncovalent binding conformation of the heterocyclic core and provides an opportunity to create structurally distinct covalent kinase inhibitors.

Sequence-dependent conformational preferences of disordered single-stranded RNA

Sequence-dependent conformational preferences of disordered single-stranded RNA

Conformational preferences of disordered proteins in biomolecular condensates

Conformational preferences of disordered proteins in biomolecular condensates

Effect of Solvents on the Conformational Preferences of Oxy[bis(methylene)]bis[(1,3-dioxan-5-yl)methanol

Effect of Solvents on the Conformational Preferences of Oxy[bis(methylene)]bis[(1,3-dioxan-5-yl)methanol

Conformational Preferences at the Glycosidic Linkage of Saccharides in Solution as Deduced from NMR Experiments and MD Simulations: Comparison to Crystal Structures.

Glycans are central to information content and regulation in biological systems. These carbohydrate molecules are active either as oligo- or polysaccharides, often in the form of glycoconjugates. The monosaccharide entities are joined by glycosidic linkages and stereochemical arrangements are of utmost importance in determining conformation and flexibility of saccharides. The conformational preferences and population distributions at the glycosidic torsion angles φ and ψ have been investigated for O-methyl glycosides of three disaccharides where the substitution takes place at a secondary alcohol, viz., in α-l-Fucp-(1→3)-β-d-Glcp-OMe, α-l-Fucp-(1→3)-α-d-Galp-OMe and α-d-Glcp-(1→4)-α-d-Galp-OMe, corresponding to disaccharide structural elements present in bacterial polysaccharides. Stereochemical differences at or adjacent to the glycosidic linkage were explored by solution state NMR spectroscopy using one-dimensional 1 H,1 H-NOESY NMR experiments to obtain transglycosidic proton-proton distances and one- and two-dimensional heteronuclear NMR experiments to obtain 3 JCH transglycosidic coupling constants related to torsion angles φ and ψ. Computed effective proton-proton distances from molecular dynamics (MD) simulations showed excellent agreement to experimentally derived distances for the α-(1→3)-linked disaccharides and revealed that for the bimodal distribution at the ψ torsion angle for the α-(1→4)-linked disaccharide experiment and simulation were at variance with each other, calling for further force field developments. The MD simulations disclosed a highly intricate inter-residue hydrogen bonding pattern for the α-(1→4)-linked disaccharide, including a nonconventional hydrogen bond between H5' in the glucosyl residue and O3 in the galactosyl residue, supported by a large downfield 1 H NMR chemical shift displacement compared to α-d-Glcp-OMe. Comparison of population distributions of the glycosidic torsion angles φ and ψ in the disaccharide entities to those of corresponding crystal structures highlighted the potential importance of solvation on the preferred conformation.

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 19GHz. 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.