Strength of NH···S Hydrogen Bonds in Methionine Residues Revealed by Gas-Phase IR/UV Spectroscopy

Molecular Dynamics Simulations Help Determine the Molecular Mechanisms of Lasioglossin-III and Its Variant Peptides' Membrane Interfacial Interactions

Lasioglossin-III (LL-III) is a potent broad-spectrum antimicrobial peptide used in diverse antimicrobial applications. In this work, coarse-grained and all-atom molecular dynamics simulation strategies were used in tandem to interpret the molecular mechanisms involved in the interfacial dynamics of LL-III and its recombinant variants during interactions with diverse cell membrane systems. Our results indicate that the membrane charges act as the driving force for initiating the membrane-peptide interactions, while the hydrophobic or van der Waals forces help to reinforce the membrane-peptide bindings. The optimized charge-hydrophobicity ratio of the LL-III peptides helps ensure their high specificity toward bacterial membranes compared to mammalian membrane systems, which also helps explain our experimental observations. Overall, we hope that our work gives new insight into the antimicrobial action of LL-III peptides and that the adopted simulation strategy will help other scientists and engineers extract maximal information from complex molecular simulations using minimal computational power.

Modulation of the strength of weak S-H⋯O hydrogen-bond: Spectroscopic study of the complexes of 2-flurothiophenol with methanol and ethanol

Spectroscopic exploration of sulfur-centered hydrogen bonding involving a thiol group (S-H) as the hydrogen bond donor is scarce in the literature. Herein, we have investigated 1:1 complexes of 2-fluorothiophenol (2-FTP) with methanol (MeOH) and ethanol (EtOH) in the gas phase to examine the physical characteristics and strength of the S-H⋯O hydrogen bond. Structures, conformations, and the strength of the S-H⋯O interaction are investigated by measuring the electronic and Infrared (IR) spectra of the two complexes employing resonant two-photon ionization, UV-UV hole-burning, and IR-UV double resonance spectroscopic techniques combined with quantum chemical calculations. Three conformers of 2-FTP⋯MeOH and two conformers of 2-FTP⋯EtOH have been detected in the experiment. A comparison of the IR spectra obtained from the experiment with those of the low-energy conformers of 2-FTP⋯MeOH and 2-FTP⋯EtOH predicted from the theory confirms that all the observed conformers of the two complexes are primarily S-H⋯O hydrogen bonded. The IR red-shifts found in the S-H stretching frequencies in 2-FTP⋯MeOH and 2-FTP⋯EtOH concerning that in 2-FTP are ∼76 and ∼88 cm-1, respectively, which are much larger than that was reported earlier in the 2-FTP⋯H2O complex (30 cm-1). The strength and physical nature of different noncovalent interactions, including the S-H⋯O hydrogen bond existing in the complexes, are further analyzed using natural bond orbital analysis, quantum theory of atoms in molecules, and localized molecular orbital-energy decomposition analysis. The current investigation reveals that the S-H⋯O hydrogen bond can be strengthened by judicial choices of the hydrogen bond acceptors of higher proton affinities.

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.

Investigating the Spectroscopy of the Gas Phase Guanine-Cytosine Pair: Keto versus Enol Configurations

We report on a vibrational study of the guanine-cytosine dimer tautomers using state-of-the-art quasiclassical trajectory and semiclassical vibrational spectroscopy. The latter includes possible quantum mechanical effects. Through an accurate comparison to the experimental spectra, we are able to shine a light on the hydrogen bond network of one of the main subunits of DNA and put the experimental assignment on a solid footing. Our calculations corroborate the experimental conclusion that the global minimum Watson-and-Crick structure is not detected in the spectra, and there is no evidence of tunnel-effect-based double proton hopping. Our accurate assignment of the spectral features may also serve as a basis for the development of precise force fields to study the guanine-cytosine dimer.

The Realm of Unconventional Noncovalent Interactions in Proteins: Their Significance in Structure and Function

Proteins and their assemblies are fundamental for living cells to function. Their complex three-dimensional architecture and its stability are attributed to the combined effect of various noncovalent interactions. It is critical to scrutinize these noncovalent interactions to understand their role in the energy landscape in folding, catalysis, and molecular recognition. This Review presents a comprehensive summary of unconventional noncovalent interactions, beyond conventional hydrogen bonds and hydrophobic interactions, which have gained prominence over the past decade. The noncovalent interactions discussed include low-barrier hydrogen bonds, C5 hydrogen bonds, C-H···π interactions, sulfur-mediated hydrogen bonds, n → π* interactions, London dispersion interactions, halogen bonds, chalcogen bonds, and tetrel bonds. This Review focuses on their chemical nature, interaction strength, and geometrical parameters obtained from X-ray crystallography, spectroscopy, bioinformatics, and computational chemistry. Also highlighted are their occurrence in proteins or their complexes and recent advances made toward understanding their role in biomolecular structure and function. Probing the chemical diversity of these interactions, we determined that the variable frequency of occurrence in proteins and the ability to synergize with one another are important not only for ab initio structure prediction but also to design proteins with new functionalities. A better understanding of these interactions will promote their utilization in designing and engineering ligands with potential therapeutic value.

Activation of the urotensin-II receptor by remdesivir induces cardiomyocyte dysfunction

Remdesivir is an antiviral drug used for COVID-19 treatment worldwide. Cardiovascular side effects have been associated with remdesivir; however, the underlying molecular mechanism remains unknown. Here, we performed a large-scale G-protein-coupled receptor screening in combination with structural modeling and found that remdesivir is a selective, partial agonist for urotensin-II receptor (UTS2R) through the Gα_i/o-dependent AKT/ERK axis. Functionally, remdesivir treatment induced prolonged field potential and APD₉₀ in human induced pluripotent stem cell (iPS)-derived cardiomyocytes and impaired contractility in both neonatal and adult cardiomyocytes, all of which mirror the clinical pathology. Importantly, remdesivir-mediated cardiac malfunctions were effectively attenuated by antagonizing UTS2R signaling. Finally, we characterized the effect of 110 single-nucleotide variants in UTS2R gene reported in genome database and found four missense variants that show gain-of-function effects in the receptor sensitivity to remdesivir. Collectively, our study illuminates a previously unknown mechanism underlying remdesivir-related cardiovascular events and that genetic variations of UTS2R gene can be a potential risk factor for cardiovascular events during remdesivir treatment, which collectively paves the way for a therapeutic opportunity to prevent such events in the future.

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.

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.

Integration of in silico methods to determine endocrine-disrupting tobacco pollutants binding potency with steroidogenic genes: comprehensive QSAR modeling and ensemble docking strategies

A myriad of tobacco-associated chemicals may have possibilities to developmental/reproductive axis and endocrine-disruption impacts. Mostly they breach the biotransformation of cholesterol in mitochondria by interfering with steroidogenic pathway genes, prompting to adverse effects in steroid biosynthesis; however, studies are scanty. The quantitative structure-activity relationship (QSAR) modeling and comparative docking strategies were used to understand structural features of dataset compounds that influence developmental/reproductive toxicity and estrogen and androgen receptor-binding abilities, and to predict binding levels of toxicants with steroidogenic acute regulatory protein (StAR) and cholesterol side-chain cleavage enzyme (CYP11A1) active sites. Developed QSAR models presented good robustness and predictive ability that were determined from the applicability domain and, clustering and classification of chemicals by performing self-organizing maps. Accordingly, the exorbitant amount of polycyclic aromatic hydrocarbons (PAHs) and a limited number of other chemicals including N-nitrosamines and nicotine was represented as potential developmental/reproductive toxicants as well as estrogen and androgen receptor binders. From the docking analysis, hydrogen bonding, nonpolar, atomic π-stacking, and π-cation interactions were found between PAHs (bay and fjord structural pockets) and functional hotspot residues of StAR and CYP11A1, which strengthened the subtle structural changes at domains. These govern barrier effects to cholesterol binding and/or locking cholesterol to complicate its ejection from the Ω1 loop of StAR, and further mitigates steroid biosynthesis through cholesterol by CYP11A1; therefore, they are presumably considered as block-cluster mechanisms. These outcomes are significant to be hopeful to estimate developmental/reproductive toxicity and endocrine-disruption activities of other environmental pollutants, and could be useful for further assessment to discover binding mechanisms of PAHs with other steroidogenesis pathway genes.

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.

Unprecedented observation and characterization of sulfur-centred bifurcated hydrogen bonds

Owing to the small electronegativity of the sulfur atom, it is commonly supposed that at most one weak H-bond can be formed between a sulfur atom and an H-bond donor. In this paper, an unprecedented 2 : 1 binding species generated from two molecules of phenol and a molecule of thioether was observed and characterized by various nuclear magnetic resonance (NMR) techniques, Fourier transform-infrared (FT-IR) techniques and density functional theory (DFT) calculations, revealing the formation of sulfur-centred O-H⋯S⋯H-O bifurcated H-bonds. This work may provide a simple and efficient method for the quantitative analysis of weak H-bonds between small organic molecules.

Unsupervised Machine Learning Neural Gas Algorithm for Accurate Evaluations of the Hessian Matrix in Molecular Dynamics

The Hessian matrix of the potential energy of molecular systems is employed not only in geometry optimizations or high-order molecular dynamics integrators but also in many other molecular procedures, such as instantaneous normal mode analysis, force field construction, instanton calculations, and semiclassical initial value representation molecular dynamics, to name a few. Here, we present an algorithm for the calculation of the approximated Hessian in molecular dynamics. The algorithm belongs to the family of unsupervised machine learning methods, and it is based on the neural gas idea, where neurons are molecular configurations whose Hessians are adopted for groups of molecular dynamics configurations with similar geometries. The method is tested on several molecular systems of different dimensionalities both in terms of accuracy and computational time versus calculating the Hessian matrix at each time-step, that is, without any approximation, and other Hessian approximation schemes. Finally, the method is applied to the on-the-fly, full-dimensional simulation of a small synthetic peptide (the 46 atom N-acetyl-l-phenylalaninyl-l-methionine amide) at the level of DFT-B3LYP-D/6-31G* theory, from which the semiclassical vibrational power spectrum is calculated.

Double-Ligand Strategy to Construct an Inhibitor-Loaded Zn-MOF and Its Corrosion Protection Ability for Aluminum Alloy 2A12

A promising double-ligand strategy for the delivery of active corrosion inhibitors by a Zn(II)-based metal-organic framework (Zn-MOF) is developed. Zn-MOF compounds were synthesized by a facile one-pot solvothermal method and characterized. The Zn-MOF is based on the corrosion inhibitor benzotriazole (BTA) and 2,5-furandicarboxylic acid (H₂FDA) ligand, which is a promising renewable building block alternative to terephthalic or isophthalic acid. The crystal structure and morphology are characterized by single-crystal X-ray diffraction analysis, powder X-ray diffraction analysis (PXRD), Fourier transform infrared spectroscopy (FT-IR), and scanning electron microscopy (SEM). The synthesized MOF crystallites are in the trigonal space group R3c with the cell parameters in a three-dimensional (3D) anionic framework. Their ability to inhibit the corrosion process of aluminum alloy 2A12 in NaCl solution was also evaluated by immersion tests in solutions with and without a MOF. The postcorrosion analysis was performed by SEM and X-ray photoelectron spectroscopy (XPS). Additional information about the inhibition efficiency was obtained by electrochemical impedance spectroscopy (EIS). The results suggest that the as-synthesized MOF can release the inhibitors and form protective layers effectively on the surface of the aluminum alloy. The use of inhibitor-loaded MOF nanocontainers provides promising opportunities for the smart delivery of inhibitors and effective corrosion protection of 2A12 aluminum alloys.

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.

A theoretical and experimental case study of the hydrogen bonding predilection of S-methylcysteine

S-containing amino acids can lead to two types of local NH···S interactions which bridge backbone NH sites to the side chain to form either intra- or inter-residue H-bonds. The present work reports on the conformational preferences of S-methyl-L-cysteine, Cys(Me), using a variety of investigating tools, ranging from quantum chemistry simulations, gas-phase UV and IR laser spectroscopy, and solution state IR and NMR spectroscopies, on model compounds comprising one or two Cys(Me) residues. We demonstrate that in gas phase and in low polarity solution, the C- and N-capped model compound for one Cys(Me) residue adopts a preferred C5-C6γ conformation which combines an intra-residue N-H···O=C backbone interaction (C5) and an inter-residue N-H···S interaction implicating the side-chain sulfur atom (C6γ). In contrast, the dominant conformation of the C- and N-capped model compound featuring two consecutive Cys(Me) residues is a regular type I β-turn. This structure is incompatible with concomitant C6γ interactions, which are no longer in evidence. Instead, C5γ interactions occur, that are fully consistent with the turn geometry and additionally stabilize the structure. Comparison with the thietane amino acid Attc, which exhibits a rigid cyclic side chain, pinpoints the significance of side chain flexibility for the specific conformational behavior of Cys(Me).

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.

Intrinsic folding of the cysteine residue: competition between folded and extended forms mediated by the -SH group

A dual microwave and optical spectroscopic study of a capped cysteine amino acid isolated in a supersonic expansion, combined with quantum chemistry modelling, enabled us to characterize the conformational preferences of Cys embedded in a protein chain. IR/UV double resonance spectroscopy provided evidence for the coexistence of two conformers, assigned to folded and extended backbones (with classical C7 and C5 backbone H-bonding respectively), each of them additionally stabilized by specific main-chain/side-chain H-bonding, where the sulfur atom essentially plays the role of H-bond acceptor. The folded structure was confirmed by microwave spectroscopy, which demonstrated the validity of the DFT-D methods currently used in the field. These structural and spectroscopic results, complemented by a theoretical Natural Bond Orbital analysis, enabled us to document the capacity of the weakly polar -CH2-SH side chain of Cys to adapt itself to the intrinsic local preferences of the peptide backbone, i.e., a γ-turn or a β-sheet extended secondary structure. The corresponding local H-bonding bridges the side chain acceptor S atom to the backbone NH donor site of the same or the next residue along the chain, through a 5- or a 6-membered ring respectively.

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.

Neutral Peptides in the Gas Phase: Conformation and Aggregation Issues

Combined IR and UV laser spectroscopic techniques in molecular beams merged with theoretical approaches have proven to be an ideal tool to elucidate intrinsic structural properties on a molecular level. It offers the possibility to analyze structural changes, in a controlled molecular environment, when successively adding aggregation partners. By this, it further makes these techniques a valuable starting point for a bottom-up approach in understanding the forces shaping larger molecular systems. This bottom-up approach was successfully applied to neutral amino acids starting around the 1990s. Ever since, experimental and theoretical methods developed further, and investigations could be extended to larger peptide systems. Against this background, the review gives an introduction to secondary structures and experimental methods as well as a summary on theoretical approaches. Vibrational frequencies being characteristic probes of molecular structure and interactions are especially addressed. Archetypal biologically relevant secondary structures investigated by molecular beam spectroscopy are described, and the influences of specific peptide residues on conformational preferences as well as the competition between secondary structures are discussed. Important influences like microsolvation or aggregation behavior are presented. Beyond the linear α-peptides, the main results of structural analysis on cyclic systems as well as on β- and γ-peptides are summarized. Overall, this contribution addresses current aspects of molecular beam spectroscopy on peptides and related species and provides molecular level insights into manifold issues of chemical and biochemical relevance.

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.

Conformation control through concurrent N-H⋯S and N-H⋯O[double bond, length as m-dash]C hydrogen bonding and hyperconjugation effects

In addition to the classical N-H⋯O[double bond, length as m-dash]C non-covalent interaction, less conventional types of hydrogen bonding, such as N-H⋯S, may play a key role in determining the molecular structure. In this work, using theoretical calculations in combination with spectroscopic analysis in both gas phase and solution phase, we demonstrate that both these H-bonding modes exist simultaneously in low-energy conformers of capped derivatives of Attc, a thietane α-amino acid. 6-Membered ring inter-residue N-H⋯S interactions (C6^γ), assisted by hyperconjugation between the thietane ring and the backbone, combine with 5-membered ring intra-residue backbone N-H⋯O[double bond, length as m-dash]C interactions (C5) to provide a C5-C6^γ feature that stabilizes a planar geometry in the monomer unit. Two contiguous C5-C6^γ features in the planar dimer implicate an unprecedented three-centre H-bond of the type C[double bond, length as m-dash]O⋯H(N)⋯SR₂, while the trimer adopts two C5-C6^γ features separated by a Ramachandran α-type backbone configuration. These low-energy conformers are fully characterized in the gas phase and support is presented for their existence in solution state.

Biological Effects of Single-Nucleotide Polymorphisms in the Drosophila melanogaster Malic Enzyme Locus

A pair of amino acid polymorphisms within the Drosophila melanogaster Malic enzyme (Men) locus presents an interesting case of genetic variation that appears to be under selection. The two alleles at each site are biochemically distinct, but their biological effects are unknown. One polymorphic site is near the active site and the other is buried within the protein. Strikingly, in twelve different populations, the first polymorphism is always found at approximately a 50:50 allelic frequency, whereas the second polymorphism is always found at approximately 90:10. The consistency of the frequencies between populations suggests that the polymorphisms are under selection and it is possible that balancing selection is at play. We used 16 lines of flies to create the nine genotypes needed to quantify both effects of the polymorphic sites and possible genetic background effects, which we found to be widespread. The alleles at each site differ, but in different biochemical characteristics. The first site significantly influences MEN K_m and V_max, whereas the second site affects the K_m and the V_max/K_m ratio (relative activity). Interestingly, the rarest allele is the most biochemically distinct. We also assayed three more distal phenotypes, triglyceride concentration, carbohydrate concentration, and longevity. In all cases, the phenotypes of the heterozygous genotypes are intermediate between those of the respective homozygotes suggesting that if balancing selection is maintaining the observed allele frequencies it is not through non-linear combinations of the biochemical phenotypes.

Semiclassical vibrational spectroscopy with Hessian databases

We report on a new approach to ease the computational overhead of ab initio "on-the-fly" semiclassical dynamics simulations for vibrational spectroscopy. The well known bottleneck of such computations lies in the necessity to estimate the Hessian matrix for propagating the semiclassical pre-exponential factor at each step along the dynamics. The procedure proposed here is based on the creation of a dynamical database of Hessians and associated molecular geometries able to speed up calculations while preserving the accuracy of results at a satisfactory level. This new approach can be interfaced to both analytical potential energy surfaces and on-the-fly dynamics, allowing one to study even large systems previously not achievable. We present results obtained for semiclassical vibrational power spectra of methane, glycine, and N-acetyl-L-phenylalaninyl-L-methionine-amide, a molecule of biological interest made of 46 atoms.

X-ray radiolytic labeling reveals the molecular basis of orange carotenoid protein photoprotection and its interactions with fluorescence recovery protein

In cyanobacterial photoprotection, the orange carotenoid protein (OCP) is photoactivated under excess light conditions and binds to the light-harvesting antenna, triggering the dissipation of captured light energy. In low light, the OCP relaxes to the native state, a process that is accelerated in the presence of fluorescence recovery protein (FRP). Despite the importance of the OCP in photoprotection, the precise mechanism of photoactivation by this protein is not well-understood. Using time-resolved X-ray-mediated in situ hydroxyl radical labeling, we probed real-time solvent accessibility (SA) changes at key OCP residues during photoactivation and relaxation. We observed a biphasic photoactivation process in which carotenoid migration preceded domain dissociation. We also observed a multiphasic relaxation process, with collapsed domain association preceding the final conformational rearrangement of the carotenoid. Using steady-state hydroxyl radical labeling, we identified sites of interaction between the FRP and OCP. In combination, the findings in this study provide molecular-level insights into the factors driving structural changes during OCP-mediated photoprotection in cyanobacteria, and furnish a basis for understanding the physiological relevance of the FRP-mediated relaxation process.

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.

Cytidine Diphosphoramidate Kinase: An Enzyme Required for the Biosynthesis of the O-Methyl Phosphoramidate Modification in the Capsular Polysaccharides of Campylobacter jejuni

Campylobacter jejuni, a leading cause of gastroenteritis, produces a capsular polysaccharide that is derivatized with a unique O-methyl phosphoramidate (MeOPN) modification. This modification contributes to serum resistance and invasion of epithelial cells. Previously, the first three biosynthetic steps for the formation of MeOPN were elucidated. The first step is catalyzed by a novel glutamine kinase (Cj1418), which catalyzes the adenosine triphosphate (ATP)-dependent phosphorylation of the amide nitrogen of l-glutamine. l-Glutamine phosphate is used by cytidine triphosphate (CTP):phosphoglutamine cytidylyltransferase (Cj1416) to displace pyrophosphate from CTP to generate cytidine diphosphate (CDP)-l-glutamine, which is then hydrolyzed by γ-glutamyl-CDP-amidate hydrolase (Cj1417) to form cytidine diphosphoramidate (CDP-NH₂). Here, we show that Cj1415 catalyzes the ATP-dependent phosphorylation of CDP-NH₂ to form 3'-phospho-cytidine-5'-diphosphoramidate. Cj1415 will also catalyze the phosphorylation of adenosine diphosphoramidate (ADP-NH₂) and uridine diphosphoramidate (UDP-NH₂) but at significantly reduced rates. It is proposed that Cj1415 be named cytidine diphosphoramidate kinase.

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.

On the turn-inducing properties of asparagine: the structuring role of the amide side chain, from isolated model peptides to crystallized proteins

The anchoring properties of an asparagine (Asn) residue to its local backbone environment in turn model peptides is characterized using gas phase laser spectroscopy and compared to crystallized protein structures.

Stereochemistry-dependent hydrogen bonds stabilise stacked conformations in jet-cooled cyclic dipeptides: (LD) vs. (LL) cyclo tyrosine-tyrosine

Tyrosine-containing cyclic dipeptides based on a diketopiperazine (DKP) ring are studied under jet-cooled conditions using resonance-enhanced multi-photon ionisation (REMPI), conformer-selective IR-UV double resonance vibrational spectroscopy and quantum chemical calculations. The conformational landscape of the dipeptide containing natural L tyrosine (Tyr), namely c-LTyr-LTyr strongly differs from that of its diastereomer c-LTyr-DTyr. A similar family of conformers exists in both systems, with one aromatic ring folded on the dipeptide DKP ring and the other one extended. Weak NHπ and CHπ interactions are observed, which are slightly different in c-LTyr-LTyr and c-LTyr-DTyr. These structures are identical to those of LL and LD cyclo diphenylalanine, which only differ from c-Tyr-Tyr by the absence of hydroxyl on the benzene rings. While this is the only conformation observed for c-LTyr-DTyr, c-LTyr-LTyr exhibits an additional form stabilised by the interaction of the two hydroxyls, in which the two aromatic rings are in a stacked geometry. Stereochemical effects are still visible in the radical cation, for which one structure is observed for c-LTyr-DTyr, while the spectrum of the c-LTyr-LTyr radical cation is explained in terms of two co-existing structures.

CO2 Complexes with Five-Membered Heterocycles: Structure, Topology, and Spectroscopic Characterization

In a first step toward the rational design of macrocyclic structures optimized for CO₂ capture, we systematically explored the potential of 30 five-membered aromatic heterocycles to establish coordinating complexes with this pollutant. The interactions between the two moieties were studied in several orientations, and the obtained complexes were analyzed in terms of electron density and vibrational fingerprint. The former is an aid to provide an in-depth knowledge of the interaction, whereas the latter should help to select structural motifs that have not only good complexation properties but also diagnostic spectroscopic signals.

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.

Chiroptical properties of dithia[3.3]cyclophanes composed of anthracene and pyridine/pyridinium moieties: A combined experimental and theoretical study

Circular dichroisms (CDs) of neutral and protonated [3.3]anthracenopyridinophane (1 and 1-H⁺ ) were investigated experimentally and theoretically. Introducing an anthracene moiety with extended conjugation affected the cyclophane structure with the bent angles being appreciably  reduced from those of parent [3.3]pyridinophane. The Cotton effects (CEs) observed at the ¹ B_b band for both 1 and 1-H⁺ were fairly strong and apparently bisignate, which, however, turned out not to be a simple exciton couplet but to be composed of multiple transitions. In contrast, the CEs were much weaker in the ¹ L_a band region. The spectral changes upon protonation were less significant compared with the parent pyridinophane, being dominated by the local transitions of anthracene. Nevertheless, the CD spectra of 1 and 1-H⁺ were well reproduced by theoretical calculations to allow us an unambiguous absolute configuration determination of the first high-performance liquid chromatography (HPLC) elute (from Chiralcel IB column) as S_p . The transannular interactions between the anthracene and pyridine/pyridinium units were examined by UV-vis and fluorescence spectroscopy to reveal a charge-transfer (CT) band in the low-energy region, particularly for 1-H⁺ . Despite the comparable CT interactions, the CE at the CT band was much stronger for the anthracenopyridinophane than for the parent pyridinophane, affording an anisotropy (g) factor as large as 4 × 10^-3 .

Improved Description of Sulfur Charge Anisotropy in OPLS Force Fields: Model Development and Parameterization

The atomic point-charge model used in most molecular mechanics force fields does not represent well the electronic anisotropy that is featured in many directional noncovalent interactions. Sulfur participates in several types of such interactions with its lone pairs and σ-holes. The current study develops a new model, via the addition of off-atom charged sites, for a variety of sulfur compounds in the OPLS-AA and OPLS/CM5 force fields to address the lack of charge anisotropy. Parameter optimization is carried out to reproduce liquid-state properties, torsional and noncovalent energetics from reliable quantum mechanical calculations, and electrostatic potentials. Significant improvements are obtained for computed free energies of hydration, reducing the mean unsigned errors from ca. 1.4 to 0.4-0.7 kcal/mol. Enhanced accuracy in directionality and energetics is also obtained for molecular complexes with sulfur-containing hydrogen and halogen bonds. Moreover, the new model reproduces the unusual conformational preferences of sulfur-containing compounds with 1,4-intramolecular chalcogen bonds. Transferability of the new force field parameters to cysteine and methionine is verified via molecular dynamic simulations of blocked dipeptides. The study demonstrates the effectiveness of using off-atom charge sites to address electronic anisotropy, and provides a parametrization methodology that can be applied to other systems.

The dual modulatory effects of efavirenz on GABAA receptors are mediated via two distinct sites

Efavirenz is a widely prescribed medicine used to treat type 1 human immunodeficiency virus (HIV-1), the most prevalent pathogenic strain of the virus responsible for the acquired immune deficiency syndrome (AIDS) pandemic. Under prescribed dosing conditions, either alone or in combination therapy, efavirenz-induced CNS disturbances are frequently reported. Efavirenz was recently reported to interact in a similar concentration range with a number of receptors, transporters and ion channels including recombinant rat α1β2γ2 GABA_A receptors whose actions were potentiated (Gatch et al., 2013; Dalwadi et al., 2016). Now we report on the molecular mechanism of efavirenz on GABA_A receptors as a function of concentration and subunit composition via whole-cell recordings of GABA-activated currents from HEK293 cells expressing varying subunit configurations of GABA_A receptors. Efavirenz elicited dual effects on the GABA response; it allosterically potentiated currents at low concentrations, whereas it inhibited currents at higher concentrations. The allosteric potentiating action on GABA_A receptors was pronounced in the α1β2γ2, α2β2γ2 and α4β2γ2 configurations, greatly diminished in the α6β2γ2 configuration, and completely absent in the α3β2γ2 or α5β2γ2 configuration. In stark contrast, the inhibitory modulation of efavirenz at higher concentrations was evident in all subunit configurations examined. Moreover, efavirenz-induced modulatory effects were dependent on GABA concentration ([GABA]), with a pronounced impact on currents activated by low [GABA] but little effect at saturating [GABA]. Mutation of a highly-conserved threonine to phenylalanine in transmembrane domain 2 of the α1 subunit abolished the inhibitory effect of efavirenz in α1β2 receptors. Finally, mutations of any of the three conserved extracellular residues in α1/2/4 subunits to the conserved residues at the corresponding positions in α3/5 subunits (i.e., R84P, M89L or I120L) completely eliminated the potentiating effect of efavirenz in α1β2γ2 configuration. These findings demonstrate that efavirenz's positive allosteric modulation of the GABA_A receptor is mediated via a novel allosteric site associated with the extracellular domain of the receptor.