Critical Assessment of the Strength of Hydrogen Bonds between the Sulfur Atom of Methionine/Cysteine and Backbone Amides in Proteins

OLIVES: A Go̅-like Model for Stabilizing Protein Structure via Hydrogen Bonding Native Contacts in the Martini 3 Coarse-Grained Force Field

Coarse-grained molecular dynamics simulations enable the modeling of increasingly complex systems at millisecond timescales. The transferable coarse-grained force field Martini 3 has shown great promise in modeling a wide range of biochemical processes, yet folded proteins in Martini 3 are not stable without the application of external bias potentials, such as elastic networks or Go̅-like models. We herein develop an algorithm, called OLIVES, which identifies native contacts with hydrogen bond capabilities in coarse-grained proteins and use it to implement a novel Go̅-like model for Martini 3. We show that the protein structure instability originates in part from the lack of hydrogen bond energy in the coarse-grained force field representation. By using realistic hydrogen bond energies obtained from literature ab initio calculations, it is demonstrated that protein stability can be recovered by the reintroduction of a coarse-grained hydrogen bond network and that OLIVES removes the need for secondary structure restraints. OLIVES is validated against known protein complexes and at the same time addresses the open question of whether there is a need for protein quaternary structure bias in Martini 3 simulations. It is shown that OLIVES can reduce the number of bias terms, hereby speeding up Martini 3 simulations of proteins by up to ≈30% on a GPU architecture compared to the established Go̅MARTINI Go̅-like model.

Structure and Substrate Specificity of S-Methyl Thiourocanate Hydratase

Nicotinamide adenine dinucleotide (NAD+) is a common cofactor in enzyme-catalyzed reactions that involve hydride transfers. In contrast, urocanase and urocanase-like enzymes use NAD+ for covalent electrophilic catalysis. Deciphering avenues by which this unusual catalytic strategy has diversified by evolution may point to approaches for the design of novel enzymes. In this report, we describe the S-methyl thiourocanate hydratase (S-Me-TUC) from Variovorax sp. RA8 as a novel member of this small family of NAD+-dependent hydratases. This enzyme catalyzes the 1,4-addition of water to S-methyl thiourocanate as the second step in the catabolism of S-methyl ergothioneine. The crystal structure of this enzyme in complex with the cofactor and a product analogue identifies critical sequence motifs that explain the narrow and nonoverlapping substrate scopes of S-methyl thiourocanate-, urocanate-, thiourocanate-, and Nτ-methyl urocanate-specific hydratases. The discovery of a S-methyl ergothioneine catabolic pathway also suggests that S-methylation or alkylation may be a significant activity in the biology of ergothioneine.

The adsorption of biogenetic odorants onto activated carbon: Adsorption characteristics and impacts of algal organic matter

Powdered activated carbon (PAC) adsorption is regarded as an efficient method for removing odorants from drinking water. However, in eutrophic aquatic environments, the presence of algal organic matter (AOM) produced by cyanobacteria considerably impedes the adsorption of odorous compounds by activated carbon. This study focused on investigating the adsorption characteristics of three representative odorants: 2-methylisoborneol (2-MIB), β-cyclocitral (β-cyclo), and butyl sulfide (BS) by PAC and the effects of AOM on the PAC adsorption of odorants. The removal of the three odorants reached 83.5-97.5% at a PAC dosage of 10 mg/L after 12 h of exposure in a competition-free scenario. The adsorption kinetics demonstrated higher conformity (R2 > 0.9) with the pseudo-second-order model, whereas the adsorption capacity exhibited stronger conformity (R2 > 0.9) with the Freundlich model. The presence of AOM resulted in varying levels of competition for PAC for the adsorption of the three odorants. As the concentration of AOM increased from 0 to 5 mg C/L, the removal of 2-MIB was the most affected (from 83.5% to 10.0%), followed by β-cyclo (from 86.6% to 55.0%), and BS (from 97.5% to 92.0%). The competitive adsorption of AOM at the molecular level was studied using density functional theory (DFT). The DFT results suggested that odorants with higher and more uniformly distributed electrostatic potentials exhibited a heightened affinity for PAC adsorption and a diminished susceptibility to disruption caused by AOM. This study provides valuable insights into the mitigation of odorous compounds during drinking water purification.

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

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

Selection of regioselective DNA aptamer for detection of homocysteine in nondeproteinized human plasma

Small molecule-binding aptamers often suffer from high cross reactivity to structure analogues in biological samples, limiting their value for clinical diagnosis. Herein, we present a method to overcome this issue, by performing binding-inhibited organic reaction-based regioselective selection of aptamers against homocysteine (Hcy), which is a marker for diagnosing many disorders including stroke and Alzheimer's. This approach has led to isolation of a DNA aptamer that binds to the alkane thiol chain of Hcy with exceptional specificity against cysteine. It also binds with oxidized Hcy at weaker affinity. Using this new aptamer, we produced a reusable fluorescent optical fiber aptasensor for direct and validated detection of both free and total Hcy in nondeproteinized patient plasma in the diagnostic concentration range. The binding site-specific aptamer selection and optical-fiber-sensing strategy can expand the practical utility of aptamers in clinical diagnosis.

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.

Sulfur-mediated chalcogen versus hydrogen bonds in proteins: a see-saw effect in the conformational space

Divalent sulfur (S) forms a chalcogen bond (Ch-bond) via its σ-holes and a hydrogen bond (H-bond) via its lone pairs. The relevance of these interactions and their interplay for protein structure and function is unclear. Based on the analyses of the crystal structures of small organic/organometallic molecules and proteins and their molecular electrostatic surface potential, we show that the reciprocity of the substituent-dependent strength of the σ-holes and lone pairs correlates with the formation of either Ch-bond or H-bond. In proteins, cystines preferentially form Ch-bonds, metal-chelated cysteines form H-bonds, while methionines form either of them with comparable frequencies. This has implications for the positioning of these residues and their role in protein structure and function. Computational analyses reveal that the S-mediated interactions stabilise protein secondary structures by mechanisms such as helix capping and protecting free β-sheet edges by negative design. The study highlights the importance of S-mediated Ch-bond and H-bond for understanding protein folding and function, the development of improved strategies for protein/peptide structure prediction and design and structure-based drug discovery.

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.

Noncovalent interactions in proteins and nucleic acids: beyond hydrogen bonding and π-stacking

Understanding the noncovalent interactions (NCIs) among the residues of proteins and nucleic acids, and between drugs and proteins/nucleic acids, etc., has extraordinary relevance in biomolecular structure and function. It helps in interpreting the dynamics of complex biological systems and enzymatic activity, which is esential for new drug design and efficient drug delivery. NCIs like hydrogen bonding (H-bonding) and π-stacking have been researchers' delight for a long time. Prominent among the recently discovered NCIs are halogen, chalcogen, pnictogen, tetrel, carbo-hydrogen, and spodium bonding, and n → π* interaction. These NCIs have caught the imaginations of various research groups in recent years while explaining several chemical and biological processes. At this stage, a holistic view of these new ideas and findings lying scattered can undoubtedly trigger our minds to explore more. The present review attempts to address NCIs beyond H-bonding and π-stacking, which are mainly n → σ*, n → π* and σ → σ* type interactions. Five of the seven NCIs mentioned earlier are linked to five non-inert end groups of the modern periodic table. Halogen (group-17) bonding is one of the oldest and most explored NCIs, which finds its relevance in biomolecules due to the phase correction and inhibitory properties of halogens. Chalcogen (group 16) bonding serves as a redox-active functional group of different active sites of enzymes and acts as a nucleophile in proteases and phosphates. Pnictogen (group 15), tetrel (group 14), triel (group 13) and spodium (group 12) bonding does exist in biomolecules. The n → π* interactions are linked to backbone carbonyl groups and protein side chains. Thus, they are crucial in determining the conformational stability of the secondary structures in proteins. In addition, a more recently discovered to and fro σ → σ* type interaction, namely carbo-hydrogen bonding, is also present in protein-ligand systems. This review summarizes these grand epiphanies routinely used to elucidate the structure and dynamics of biomolecules, their enzymatic activities, and their application in drug discovery. It also briefs about the future perspectives and challenges posed to the spectroscopists and theoreticians.

Carbon-Centered Hydrogen Bonds in Proteins

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

Can two wrongs make a right? F508del-CFTR ion channel rescue by second-site mutations in its transmembrane domains

Deletion of phenylalanine 508 (F508del) in the cystic fibrosis transmembrane conductance regulator (CFTR) anion channel is the most common cause of cystic fibrosis. The F508 residue is located on nucleotide-binding domain 1 (NBD1) in contact with the cytosolic extensions of the transmembrane helices, in particular intracellular loop 4 (ICL4). To investigate how absence of F508 at this interface impacts the CFTR protein, we carried out a mutagenesis scan of ICL4 by introducing second-site mutations at 11 positions in cis with F508del. Using an image-based fluorescence assay, we measured how each mutation affected membrane proximity and ion-channel function. The scan strongly validated the effectiveness of R1070W at rescuing F508del defects. Molecular dynamics simulations highlighted two features characterizing the ICL4/NBD1 interface of F508del/R1070W-CFTR: flexibility, with frequent transient formation of interdomain hydrogen bonds, and loosely stacked aromatic sidechains (F1068, R1070W, and F1074, mimicking F1068, F508, and F1074 in WT CFTR). F508del-CFTR displayed a distorted aromatic stack, with F1068 displaced toward the space vacated by F508, while in F508del/R1070F-CFTR, which largely retained F508del defects, R1070F could not form hydrogen bonds and the interface was less flexible. Other ICL4 second-site mutations which partially rescued F508del-CFTR included F1068M and F1074M. Methionine side chains allow hydrophobic interactions without the steric rigidity of aromatic rings, possibly conferring flexibility to accommodate the absence of F508 and retain a dynamic interface. These studies highlight how both hydrophobic interactions and conformational flexibility might be important at the ICL4/NBD1 interface, suggesting possible structural underpinnings of F508del-induced dysfunction.

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.

Hydrogen Bonding with Polonium

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

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.

Structural insights into the enhanced thermostability of cysteine substitution mutants of L-histidine decarboxylase from Photobacterium phosphoreum

Enzymatic amino acid assays are important in physiological research and clinical diagnostics because abnormal amino acid concentrations in biofluids are associated with various diseases. L-histidine decarboxylase from Photobacterium phosphoreum (PpHDC) is a pyridoxal 5'-phosphate-dependent enzyme and a candidate for use in an L-histidine quantitation assay. Previous cysteine substitution experiments demonstrated that the PpHDC C57S mutant displayed improved long-term storage stability and thermostability when compared with those of the wild-type enzyme. In this study, combinational mutation experiments of single cysteine substitution mutants of PpHDC were performed, revealing that the PpHDC C57S/C101V/C282V mutant possessed the highest thermostability. The stabilizing mechanism of these mutations were elucidated by solving the structures of PpHDC C57S and C57S/C101V/C282V mutants by X-ray crystallography. In the crystal structures, two symmetry-related PpHDC molecules form a domain-swapped homodimer. The side chain of S57 is solvent exposed in the structure, indicating that the C57S mutation eliminates chemical oxidation or disulfide bond formation with a free thiol group, thereby providing greater stability. Residues 101 and 282 form hydrophobic interactions with neighboring hydrophobic residues. Mutations C101V and C282V enhanced thermostability of PpHDC by filling a cavity present in the hydrophobic core (C101V) and increasing hydrophobic interactions.

Intramolecular O-H⋯S hydrogen bonding in threefold symmetry: Line broadening dynamics from ultrafast 2DIR-spectroscopy and ab initio calculations

The dynamics of intramolecular hydrogen-bonding involving sulfur atoms as acceptors is studied using two-dimensional infrared (2DIR) spectroscopy. The molecular system is a tertiary alcohol whose donating hydroxy group is embedded in a hydrogen-bond potential with torsional C₃-symmetry about the carbon-oxygen bond. The linear and 2DIR-spectra recorded in the OH-stretching region of the alcohol can be simulated very well using Kubo's line shape theory based on the cumulant expansion for evaluating the linear and nonlinear optical response functions. The correlation function for OH-stretching frequency fluctuations reveals an ultrafast component decaying with a time constant of 700 fs, which is in line with the apparent decay of the center line slopes averaged over absorption and bleach/emission signals. In addition, a quasi-static inhomogeneity is detected, which prevents the 2DIR line shape to fully homogenize within the observation window of 4 ps. The experimental data were then analyzed in more detail using a full ab initio approach that merges time-dependent structural information from classical molecular dynamics (MD) simulations with an OH-stretching frequency map derived from density functional theory (DFT). The latter method was also used to obtain a complementary transition dipole map to account for non-Condon effects. The 2DIR-spectra obtained from the MD/DFT method are in good agreement with the experimental data at early waiting delays, thereby corroborating an assignment of the fast decay of the correlation function to the dynamics of hydrogen-bond breakage and formation.

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.

An intraresidue H-bonding motif in selenocysteine and cysteine, revealed by gas phase laser spectroscopy and quantum chemistry calculations

Models of protein chains containing a seleno-cysteine (Sec) residue have been investigated by gas phase laser spectroscopy in order to document the effect of the H-bonding properties of the SeH group in the folding of the Sec side chain, by comparison with recent data on Ser- and Cys-containing sequences. Experimental data, complemented by quantum chemistry calculations and natural bonding orbital (NBO) analyses, are interpreted in terms of the formation of a so-called 5γ intra-residue motif, which bridges the acceptor chalcogen atom of the side chain to the NH bond of the same residue. This local structure, in which the O/S/Se atom is close to the plane of the N-terminal side amide, is constrained by local backbone-side chain hyperconjugation effects involving the S and Se atoms. Theoretical investigations of the Cys/Sec side chain show that (i) this 5γ motif is an intrinsic feature of these residues, (ii) the corresponding H-bond is strongly non-linear and intrinsically weak, (iii) but enhanced by γ- and β-turn secondary structures, which promote a more favorable 5γ H-bonding approach and distance. The resulting H-bonds are slightly stronger in selenocysteine than in cysteine, but nearly inexistent in serine, whose side chain in contrast behaves as a H-bonding donor. The modest spectral shifts of the Cys/Sec NH stretches measured experimentally reflect the moderate strength of the 5γ H-bonding, in agreement with the correlation obtained with a NBO-based H-bond strength indicator. The evolution along the Ser, Cys and Sec series emphasizes the compromise between the several factors that control the H-bonding in a hyperconjugation-constrained geometry, among them the chalcogen van der Waals and covalent radii. It also illustrates the 5γ H-bond enhancements with the Sec and Cys residues favoured by the constraints imposed by the γ- and β-turn structures of the peptide chain.

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.

In silico allicin induced S-thioallylation of SARS-CoV-2 main protease

Coronavirus disease 2019 (COVID-19) is an ongoing pandemic caused due to new coronavirus infection with ³716075 deaths across the world as reported by the World Health Organization (WHO). SARS-CoV-2 main protease (M^pro) plays a vital role in the replication of coronavirus and thus an attractive target for the screening of inhibitors for the therapy of COVID-19. The preclinical drugs ebselen and PX-12 are potent inhibitors of SARS-CoV-2 M^pro and covalently modifies the active site Cys-145 residue of M^pro through selenosulfide/disulfide. In the current report, using virtual screening methods, reactive sulfur species allicin is subjecting for covalent docking at the active site of SARS-CoV-2 M^pro using PX-12 as a benchmark reference compound. The results indicate that allicin induces dual S-thioallylation of Cys-145 and Cys-85/ Cys-156 residues of SARS-CoV-2 M^pro. Using density functional theory (DFT), Gibbs free energy change (DG) is calculated for the putative reactions between N-acetylcysteine amide thiol and allicin/allyl sulfenic acid. The overall reaction is exergonic and allyl disulfide of Cys-145 residue of M^pro is involved in a sulfur mediated hydrogen bond. The results indicate that allicin causes dual S-thioallylation of SARS-CoV-2 M^pro which may be of interest for treatment and attenuation of ongoing coronavirus infection.

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.

Structural Characterization, DFT Calculation, NCI, Scan-Rate Analysis and Antifungal Activity against Botrytis cinerea of (E)-2-{[(2-Aminopyridin-2-yl)imino]-methyl}-4,6-di-tert-butylphenol (Pyridine Schiff Base)

Botrytis cinerea is a ubiquitous necrotrophic filamentous fungal phytopathogen that lacks host specificity and can affect more than 1000 different plant species. In this work, we explored L1 [(E)-2-{[(2-aminopyridin-2-yl)imino]-methyl}-4,6-di-tert-butylphenol], a pyridine Schiff base harboring an intramolecular bond (IHB), regarding their antifungal activity against Botrytis cinerea. Moreover, we present a full characterization of the L1 by NMR and powder diffraction, as well as UV–vis, in the presence of previously untested different organic solvents. Complementary time-dependent density functional theory (TD-DFT) calculations were performed, and the noncovalent interaction (NCI) index was determined. Moreover, we obtained a scan-rate study on cyclic voltammetry of L1. Finally, we tested the antifungal activity of L1 against two strains of Botrytis cinerea (B05.10, a standard laboratory strain; and A1, a wild type strains isolated from Chilean blueberries). We found that L1 acts as an efficient antifungal agent against Botrytis cinerea at 26 °C, even better than the commercial antifungal agent fenhexamid. Although the antifungal activity was also observed at 4 °C, the effect was less pronounced. These results show the high versatility of this kind of pyridine Schiff bases in biological applications.

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

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

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

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

Sulfur hydrogen bonding and internal dynamics in the monohydrates of thenyl mercaptan and thenyl alcohol

Water forms weak H-bonds with thenyl compounds, simultaneously retaining internal mobility in the dimer.

Internal Motions and Sulfur Hydrogen Bonding in Methyl 3-Mercaptopropionate

The effect of sulfur hydrogen bonding on the conformational equilibrium of methyl 3-mercaptopropionate was investigated using microwave spectroscopy in a supersonic jet expansion. The two most stable conformers (I and II) were assigned in the rotational spectra, and complex splitting patterns owing to the methyl internal rotation and SH tunneling motion were resolved and analyzed in detail. For both conformers, the experimental torsional barriers for the methyl top are similar and about 5.1 kJ mol^-1, revealing that their geometrical differences do not affect the methyl internal rotation. The experimentally derived rotational and centrifugal distortion constants, along with the methyl internal rotation barriers, are discussed and compared with results from density functional theory and ab initio calculations. Quantum theory of atoms in molecules, noncovalent interactions, and natural bond orbital analyses show that the global minimum geometry (I), which has the thiol hydrogen oriented toward the carbonyl of the ester, is stabilized by an SH···O=C hydrogen bond. The presence of a hydrogen bond is confirmed by the derivation of an accurate experimental geometry that reveals a hydrogen bond distance and S-H-O angle of 2.515(4) Å and 117.4(1)°, respectively. These results are key benchmarks to expand the current knowledge of sulfur hydrogen bonds and the relationship between internal motions and conformational preferences in esters.

A mutation in MAP2 is associated with prenatal hair follicle density

Hairlessness is usually a rare trait in pigs; however, in this study, we found hairless (HR) pigs at a relatively high frequency in 1 pig herd. We observed that, the lower hair shaft density of HR pigs could be mainly attributed to the lower hair follicle density, and during the embryonic period, d 39-45 were a critical stage for the formation of the hair follicle. In this regard, d 41 during gestation was a particularly important point. Hair follicle morphogenesis occurring at an early stage of embryo development is similar to humans and mice. Further analyses of association studies based on single-nucleotide polymorphism chip as well as sequence data, mRNA sequencing, immunohistochemistry, and comparative genomics demonstrated that microtubule-associated protein 2 (MAP2) is a key gene responsible for hair follicle density and 1 missense mutation of A-to-G at rs328005415 in MAP2, causing a valine-to-methionine substitution leads to the HR phenotype. Considering the high homology between pigs and humans, our research has some significance for the study of the mechanisms of skin development, hair morphogenesis, and hair loss in humans by showing that the pig may be a more appropriate model in which to study these processes.-Jiang, Y., Jiang, Y., Zhang, H., Mei, M., Song, H., Ma, X., Jiang, L., Yu, Z., Zhang, Q., Ding, X. A mutation in MAP2 is associated with prenatal hair follicle density.

Understanding Curly Hair Mechanics: Fiber Strength

The relationship between the geometric and mechanical profiles of hair fibers has been studied, with special focus on curly samples. Incidental observations pointed to a significantly different viscoelastic character with varying curliness. Further investigations confirmed initial observations, showing an initial distinct toe region behavior for curly fibers on the stress-strain plot, which is absent for straight fibers. This behavior suggested a difference in the viscoelastic nature of the curly fiber that is linked to mechanical energy stored in the fiber. Results also suggest that the strength of hair depends on two main components, and further pointed out that de facto methods of tensile testing may erode curly fiber strength during preparation. The main outcome of this study is that the tensile strength (σ_T) of hair fibers is composed of two (rather than one main) components, namely the toe region (σ_t) and the elastic region (σ_ε), so that: σ_T=σ_t+σ_ε. For noncurly fibers, the greatest part of fiber strength is derived from σ_ε, while σ_t ≈ 0. For curly fibers, σ_t (i.e., springiness) adds significantly to the overall strength, even though σ_ε remains the major contributor. Although these results require validation in larger studies, they are significant in the current understanding of curly hair. Also, they may represent a fundamental shift from the current understanding of tensile testing of human hair in general.

Shining light on cysteine modification: connecting protein conformational dynamics to catalysis and regulation

Cysteine is a rare but functionally important amino acid that is often subject to covalent modification. Cysteine oxidation plays an important role in many human disease processes, and basal levels of cysteine oxidation are required for proper cellular function. Because reactive cysteine residues are typically ionized to the thiolate anion (Cys-S-), their formation of a covalent bond alters the electrostatic and steric environment of the active site. X-ray-induced photo-oxidation to sulfenic acids (Cys-SOH) can recapitulate some aspects of the changes that occur under physiological conditions. Here we propose how site-specific cysteine photo-oxidation can be used to interrogate ensuing changes in protein structure and dynamics at atomic resolution. Although this powerful approach can connect cysteine covalent modification to global protein conformational changes and function, careful biochemical validation must accompany all such studies to exclude misleading artifacts. New types of X-ray crystallography experiments and powerful computational methods are creating new opportunities to connect conformational dynamics to catalysis for the large class of systems that use covalently modified cysteine residues for catalysis or regulation.

A local fingerprint for hydrophobicity and hydrophilicity: From methane to peptides

An important characteristic that determines the behavior of a solute in water is whether it is hydrophobic or hydrophilic. The traditional classification is based on chemical experience and heuristics. However, this does not reveal how the local environment modulates this important property. We present a local fingerprint for hydrophobicity and hydrophilicity inspired by the two body contribution to the entropy. This fingerprint is an inexpensive, quantitative, and physically meaningful way of studying hydrophilicity and hydrophobicity that only requires as input the water-solute radial distribution functions. We apply our fingerprint to octanol, benzene, and 20 proteinogenic amino acids. Our measure of hydrophilicity is coherent with chemical experience, and moreover, it also shows how the character of an atom can change as its environment is changed. Finally, we use the fingerprint as a collective variable in a funnel metadynamics simulation of a host-guest system. The fingerprint serves as a desolvation collective variable that enhances transitions between the bound and unbound states.

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.

Biochemical characterization of three midgut chitin deacetylases of the Lepidopteran insect Bombyx mori

Peritrophic membrane (PM) is a chitin and protein-containing extracellular matrix that lines the midgut in most insect species, functioning as a barrier to exogenous toxins and pathogens. Midgut chitin deacetylases (CDAs) are chitin-modifying enzymes known to alter the mechanical property and permeability of PM. However, biochemical properties and specific roles of these enzymes remain elusive. In this study, the midgut-expressed CDAs (BmCDA6, BmCDA7 and BmCDA8) from Bombyx mori were cloned, recombinantly expressed and purified and their enzymatic activities toward PM chitin were determined. Of the three enzymes, BmCDA7 exhibited the highest activity (0.284 μmol/min/μmol), while BmCDA8 showed lower activity of 0.061 μmol/min/μmol. BmCDA6 was inactive towards PM chitin. Gene expression patterns indicated that although all three CDA genes were specifically expressed in the anterior midgut, they differed in their temporal expression patterns. BmCDA6 was expressed almost exclusively at the mid-molt stage, the stage when the PM was thick and with multiple chitin layers. Unlike BmCDA6, high expression levels of BmCDA7 and BmCDA8 were observed only at the feeding stage, the stage when the PM is thin and with fewer chitin layers. The different gene expression patterns and biochemical characteristics provide new information about the functional specialization among BmCDA6, BmCDA7 and BmCDA8 proteins.

Transfer of chirality in N-triphenylacetylamino acids and chiral derivatives of N-triphenylacetyl Gly–Gly dipeptide and control of their assembly with steric constraints

The N-triphenylacetyl group is utilized as a reporter of chirality and as a supramolecular protecting group for α-amino acid and peptide derivatives.

Hydrogen bond-linked pathways of peptide units and polar groups of amino acid residues suitable for electron transfer in cytochrome c proteins

Electron transfer occurs through heme-Fe across the cytochrome c protein. The current models of long range electron transfer pathways in proteins include covalent σ-bonds, van der Waals forces, and through space jump. Hydrogen-bond-linked pathways of delocalized electron units in peptide units and polar side chains of amino acid residues in proteins and internal water molecules are better suited for intramolecular atom-to-atom electron transfer in proteins. Crystal structures of cytochrome c proteins from horse (1HRC), tuna (3CYT), rice (1CCR), and yeast (3CX5) were analyzed using pymol software for 'Hydrogen Bonds' marking the polar atoms within the distance of 2.6-3.3 Å and tracing the atom-to-atom pathways linked by hydrogen bonds. Pathways of hydrogen-bond-linked peptide units, polar side chains of the amino acid residues, and buried water molecules connect heme-Fe through axially coordinated Met80-S and His18-N have been traced in cytochrome c proteins obtained from horse, tuna, rice and yeast with an identical hydrogen-bonded sequence around the heme-Fe: Asn-N-water-O-Tyr-O-Met-S-heme-Fe-His (HN-C=N)-Pro-Asn-Pro-Gly (peptide unit, HN-C=O)-water-O. More than half of the amino acid residues in these pathways are among the conserved list and delocalized electron units, internal water molecules and hydrogen bonds are conspicuous by their presence.

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.

Exome Sequencing Identifies a Novel MAP3K14 Mutation in Recessive Atypical Combined Immunodeficiency

Primary immunodeficiency disorders (PIDs) render patients vulnerable to infection with a wide range of microorganisms and thus provide good in vivo models for the assessment of immune responses during infectious challenges. Priming of the immune system, especially in infancy, depends on different environmental exposures and medical practices. This may determine the timing and phenotype of clinical appearance of immune deficits as exemplified with early exposure to Bacillus Calmette-Guérin (BCG) vaccination and dissemination in combined immunodeficiencies. Varied phenotype expression poses a challenge to identification of the putative immune deficit. Without the availability of genomic diagnosis and data analysis resources and with limited capacity for functional definition of immune pathways, it is difficult to establish a definitive diagnosis and to decide on appropriate treatment. This study describes the use of exome sequencing to identify a homozygous recessive variant in MAP3K14, NIK^Val345Met, in a patient with combined immunodeficiency, disseminated BCG-osis, and paradoxically elevated lymphocytes. Laboratory testing confirmed hypogammaglobulinemia with normal CD19, but failed to confirm a definitive diagnosis for targeted treatment decisions. NIK^Val345Met is predicted to be deleterious and pathogenic by two in silico prediction tools and is situated in a gene crucial for effective functioning of the non-canonical nuclear factor-kappa B signaling pathway. Functional analysis of NIK^Val345Met- versus NIK^WT-transfected human embryonic kidney-293T cells showed that this mutation significantly affects the kinase activity of NIK leading to decreased levels of phosphorylated IkappaB kinase-alpha (IKKα), the target of NIK. BCG-stimulated RAW264.7 cells transfected with NIK^Val345Met also presented with reduced levels of phosphorylated IKKα, significantly increased p100 levels and significantly decreased p52 levels compared to cells transfected with NIK^WT. Ideally, these experiments would have been conducted in patient-derived immune cells, but we were unable to source these cells from the patient. The functional analysis described in this paper supports previous illustrations of the importance of NIK in human immune responses and demonstrates the involvement of function-altering mutations in MAP3K14 in PIDs. The genomic approach used for this patient demonstrates its value in the diagnosis of an unusual PID and as a tool for detecting rarer mutations to help guide treatment approaches.