Different types of interactions involving cysteine sulfhydryl group in proteins

Leveraging a Phage-Encoded Noncanonical Amino Acid: A Novel Pathway to Potent and Selective Epigenetic Reader Protein Inhibitors

Epigenetic reader proteins interpret histone epigenetic marks to regulate gene expression. Given their vital roles and the link between their dysfunction and various diseases, these proteins present compelling targets for therapeutic interventions. Nevertheless, designing selective inhibitors for these proteins poses significant challenges, primarily due to their unique properties such as shallow binding sites and similarities with homologous proteins. To overcome these challenges, we propose an innovative strategy that uses phage display with a genetically encoded noncanonical amino acid (ncAA) containing an epigenetic mark. This ncAA guides binding to the reader protein's active site, allowing the identification of peptide inhibitors with enhanced affinity and selectivity. In this study, we demonstrate this novel approach's effectiveness by identifying potent inhibitors for the ENL YEATS domain that plays a critical role in leukemogenesis. Our strategy involved genetically incorporating Nε-butyryl-l-lysine (BuK), known for its binding to ENL YEATS, into a phage display library for enriching the pool of potent inhibitors. One resultant hit was further optimized by substituting BuK with other pharmacophores to exploit a unique π-π-π stacking interaction with ENL YEATS. This led to the creation of selective ENL YEATS inhibitors with a KD value of 2.0 nM and a selectivity 28 times higher for ENL YEATS than its close homologue AF9 YEATS. One such inhibitor, tENL-S1f, demonstrated robust cellular target engagement and on-target effects to inhibit leukemia cell growth and suppress the expression of ENL target genes. As a pioneering study, this work opens up extensive avenues for the development of potent and selective peptidyl inhibitors for a broad spectrum of epigenetic reader proteins.

GlyNAC (Glycine and N-Acetylcysteine) Supplementation in Old Mice Improves Brain Glutathione Deficiency, Oxidative Stress, Glucose Uptake, Mitochondrial Dysfunction, Genomic Damage, Inflammation and Neurotrophic Factors to Reverse Age-Associated Cognitive Decline: Implications for Improving Brain Health in Aging

Cognitive decline frequently occurs with increasing age, but mechanisms contributing to age-associated cognitive decline (ACD) are not well understood and solutions are lacking. Understanding and reversing mechanisms contributing to ACD are important because increased age is identified as the single most important risk factor for dementia. We reported earlier that ACD in older humans is associated with glutathione (GSH) deficiency, oxidative stress (OxS), mitochondrial dysfunction, glucose dysmetabolism and inflammation, and that supplementing GlyNAC (glycine and N-acetylcysteine) improved these defects. To test whether these defects occur in the brain in association with ACD, and could be improved/reversed with GlyNAC supplementation, we studied young (20-week) and old (90-week) C57BL/6J mice. Old mice received either regular or GlyNAC supplemented diets for 8 weeks, while young mice received the regular diet. Cognition and brain outcomes (GSH, OxS, mitochondrial energetics, autophagy/mitophagy, glucose transporters, inflammation, genomic damage and neurotrophic factors) were measured. Compared to young mice, the old-control mice had significant cognitive impairment and multiple brain defects. GlyNAC supplementation improved/corrected the brain defects and reversed ACD. This study finds that naturally-occurring ACD is associated with multiple abnormalities in the brain, and provides proof-of-concept that GlyNAC supplementation corrects these defects and improves cognitive function in aging.

Interplay between hydrogen and chalcogen bonds in cysteine

Protein structures are stabilized by several types of chemical interactions between amino acids, which can compete with each other. This is the case of chalcogen and hydrogen bonds formed by the thiol group of cysteine, which can form three hydrogen bonds with one hydrogen acceptor and two hydrogen donors and a chalcogen bond with a nucleophile along the extension of the CS bond. A survey of the Protein Data Bank shows that hydrogen bonds are about 40-50 more common than chalcogen bonds, suggesting that they are stronger and, consequently, prevail, though not always. It is also observed that frequently a thiol group that forms a chalcogen bond is also involved, as a hydrogen donor, in a hydrogen bond.

The Potential of a Protein Model Synthesized Absent of Methionine

Methionine is an amino acid long thought to be essential, but only in the case of protein synthesis initiation. In more recent years, methionine has been found to play an important role in antioxidant defense, stability, and modulation of cell and protein activity. Though these findings have expanded the previously held sentiment of methionine having a singular purpose within cells and proteins, the essential nature of methionine can still be challenged. Many of the features that give methionine its newfound functions are shared by the other sulfur-containing amino acid: cysteine. While the antioxidant, stabilizing, and cell/protein modulatory functions of cysteine have already been well established, recent findings have shown a similar hydrophobicity to methionine which suggests cysteine may be able to replace methionine in all functions outside of protein synthesis initiation with little effect on cell and protein function. Furthermore, a number of novel mechanisms for alternative initiation of protein synthesis have been identified that suggest a potential to bypass the traditional methionine-dependent initiation during times of stress. In this review, these findings are discussed with a number of examples that demonstrate a potential model for synthesizing a protein in the absence of methionine.

Unprotected peptide macrocyclization and stapling via a fluorine-thiol displacement reaction

We report the discovery of a facile peptide macrocyclization and stapling strategy based on a fluorine thiol displacement reaction (FTDR), which renders a class of peptide analogues with enhanced stability, affinity, cellular uptake, and inhibition of cancer cells. This approach enabled selective modification of the orthogonal fluoroacetamide side chains in unprotected peptides in the presence of intrinsic cysteines. The identified benzenedimethanethiol linker greatly promoted the alpha helicity of a variety of peptide substrates, as corroborated by molecular dynamics simulations. The cellular uptake of benzenedimethanethiol stapled peptides appeared to be universally enhanced compared to the classic ring-closing metathesis (RCM) stapled peptides. Pilot mechanism studies suggested that the uptake of FTDR-stapled peptides may involve multiple endocytosis pathways in a distinct pattern in comparison to peptides stapled by RCM. Consistent with the improved cell permeability, the FTDR-stapled lead Axin and p53 peptide analogues demonstrated enhanced inhibition of cancer cells over the RCM-stapled analogues and the unstapled peptides.

Strategies capable of stapling unprotected peptides in a straightforward, chemoselective, and clean manner, as well as promoting cellular uptake are of great interest. Here the authors report a peptide macrocyclization and stapling strategy which satisfies those criteria, based on a fluorine thiol displacement reaction.

GlyNAC (Glycine and N-Acetylcysteine) Supplementation Improves Impaired Mitochondrial Fuel Oxidation and Lowers Insulin Resistance in Patients with Type 2 Diabetes: Results of a Pilot Study

Patients with type 2 diabetes (T2D) are known to have mitochondrial dysfunction and increased insulin resistance (IR), but the underlying mechanisms are not well understood. We reported previously that (a) adequacy of the antioxidant glutathione (GSH) is necessary for optimal mitochondrial fatty-acid oxidation (MFO); (b) supplementing the GSH precursors glycine and N-acetylcysteine (GlyNAC) in mice corrected GSH deficiency, reversed impaired MFO, and lowered oxidative stress (OxS) and IR; and (c) supplementing GlyNAC in patients with T2D improved GSH synthesis and concentrations, and lowered OxS. However, the effect of GlyNAC on MFO, MGO (mitochondrial glucose oxidation), IR and plasma FFA (free-fatty acid) concentrations in humans with T2D remains unknown. This manuscript reports the effect of supplementing GlyNAC for 14-days on MFO, MGO, IR and FFA in 10 adults with T2D and 10 unsupplemented non-diabetic controls. Fasted T2D participants had 36% lower MFO (p < 0.001), 106% higher MGO (p < 0.01), 425% higher IR (p < 0.001) and 76% higher plasma FFA (p < 0.05). GlyNAC supplementation significantly improved fasted MFO by 30% (p < 0.001), lowered MGO by 47% (p < 0.01), decreased IR by 22% (p < 0.01) and lowered FFA by 25% (p < 0.01). These results provide proof-of-concept that GlyNAC supplementation could improve mitochondrial dysfunction and IR in patients with T2D, and warrant additional research.

Expression and Refolding of the Plant Chitinase From Drosera capensis for Applications as a Sustainable and Integrated Pest Management

Recently, the study of chitinases has become an important target of numerous research projects due to their potential for applications, such as biocontrol pest agents. Plant chitinases from carnivorous plants of the genus Drosera are most aggressive against a wide range of phytopathogens. However, low solubility or insolubility of the target protein hampered application of chitinases as biofungicides. To obtain plant chitinase from carnivorous plants of the genus Drosera in soluble form in E.coli expression strains, three different approaches including dialysis, rapid dilution, and refolding on Ni-NTA agarose to renaturation were tested. The developed « Rapid dilution » protocol with renaturation buffer supplemented by 10% glycerol and 2M arginine in combination with the redox pair of reduced/oxidized glutathione, increased the yield of active soluble protein to 9.5 mg per 1 g of wet biomass. A structure-based removal of free cysteines in the core domain based on homology modeling of the structure was carried out in order to improve the soluble of chitinase. One improved chitinase variant (C191A/C231S/C286T) was identified which shows improved expression and solubility in E. coli expression systems compared to wild type. Computational analyzes of the wild-type and the improved variant revealed overall higher fluctuations of the structure while maintaining a global protein stability. It was shown that free cysteines on the surface of the protein globule which are not involved in the formation of inner disulfide bonds contribute to the insolubility of chitinase from Drosera capensis. The functional characteristics showed that chitinase exhibits high activity against colloidal chitin (360 units/g) and high fungicidal properties of recombinant chitinases against Parastagonospora nodorum. Latter highlights the application of chitinase from D. capensis as a promising enzyme for the control of fungal pathogens in agriculture.

A molecular dynamics study of protein denaturation induced by sulfonate-based surfactants

Microsecond timescale explicit-solvent atomistic simulations were carried out to investigate how anionic surfactants modulate protein structure and dynamics. We found that lysozyme undergoes near-complete denaturation at the high concentration (> 0.1 M) of sodium pentadecyl sulfonate (SPDS), while only partial denaturation occurs at the concentration slightly below 0.1 M. In large part, protein denaturation is structurally manifested by disappearance of helical segments and loss of tertiary interactions. The computational prediction of the extent of burial of cysteine residues was experimentally validated by measuring the accessibility of the respective sulfhydryl groups. Overall, our work indicates an interesting synergy between electrostatic and hydrophobic contributions to lysozyme's denaturation process by anionic surfactants. In fact, first disulfide bridges and hydrogen bonds from protein surface to SPDS head groups loosen the protein globule followed by fuller denaturation via insertion of the surfactant's hydrophobic tails into the protein core.

Structure of Usutu virus SAAR-1776 displays fusion loop asymmetry

Usutu virus (USUV) is an emerging arbovirus in Europe that has been increasingly identified in asymptomatic humans and donated blood samples and is a cause of increased incidents of neuroinvasive human disease. Treatment or prevention options for USUV disease are currently nonexistent, the result of a lack of understanding of the fundamental elements of USUV pathogenesis. Here, we report two structures of the mature USUV virus, determined at a resolution of 2.4 Å, using single-particle cryogenic electron microscopy. Mature USUV is an icosahedral shell of 180 copies of envelope (E) and membrane (M) proteins arranged in the classic herringbone pattern. However, unlike previous reports of flavivirus structures, we observe virus subpopulations and differences in the fusion loop disulfide bond. Presence of a second, unique E glycosylation site could elucidate host interactions, contributing to the broad USUV tissue tropism. The structures provide a basis for exploring USUV interactions with glycosaminoglycans and lectins, the role of the RGD motif as a receptor, and the inability of West Nile virus therapeutic antibody E16 to neutralize the mature USUV strain SAAR-1776. Finally, we identify three lipid binding sites and predict key residues that likely participate in virus stability and flexibility during membrane fusion. Our findings provide a framework for the development of USUV therapeutics and expand the current knowledge base of flavivirus biology.

The susceptibility of disulfide bonds towards radiation damage may be explained by S⋯O interactions

Radiation-induced damage to protein crystals during X-ray diffraction data collection is a major impediment to obtaining accurate structural information on macromolecules. Some of the specific impairments that are inflicted upon highly brilliant X-ray irradiation are metal-ion reduction, disulfide-bond cleavage and a loss of the integrity of the carboxyl groups of acidic residues. With respect to disulfide-bond reduction, previous results have indicated that not all disulfide bridges are equally susceptible to damage. A careful analysis of the chemical environment of disulfide bonds in the structures of elastase, lysozyme, acetylcholinesterase and other proteins suggests that S-S bonds which engage in a close contact with a carbonyl O atom along the extension of the S-S bond vector are more susceptible to reduction than the others. Such an arrangement predisposes electron transfer to occur from the O atom to the disulfide bond, leading to its reduction. The interaction between a nucleophile and an electrophile, akin to hydrogen bonding, stabilizes protein structures, but it also provides a pathway of electron transfer to the S-S bond, leading to its reduction during exposure of the protein crystal to an intense X-ray beam. An otherwise stabilizing interaction can thus be the cause of destabilization under the condition of radiation exposure.

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.

Increasing the Affinity of an O-Antigen Polysaccharide Binding Site in Shigella flexneri Bacteriophage Sf6 Tailspike Protein

Broad and unspecific use of antibiotics accelerates spread of resistances. Sensitive and robust pathogen detection is thus important for a more targeted application. Bacteriophages contain a large repertoire of pathogen-binding proteins. These tailspike proteins (TSP) often bind surface glycans and represent a promising design platform for specific pathogen sensors. We analysed bacteriophage Sf6 TSP that recognizes the O-polysaccharide of dysentery-causing Shigella flexneri to develop variants with increased sensitivity for sensor applications. Ligand polyrhamnose backbone conformations were obtained from 2D 1 H,1 H-trNOESY NMR utilizing methine-methine and methine-methyl correlations. They agreed well with conformations obtained from molecular dynamics (MD), validating the method for further predictions. In a set of mutants, MD predicted ligand flexibilities that were in good correlation with binding strength as confirmed on immobilized S. flexneri O-polysaccharide (PS) with surface plasmon resonance. In silico approaches combined with rapid screening on PS surfaces hence provide valuable strategies for TSP-based pathogen sensor design.

Amino acid-modified PAMAM dendritic nanocarriers as effective chemotherapeutic drug vehicles in cancer treatment: a study using zebrafish as a cancer model

The use of nanomaterials for drug delivery offers many advantages including the targeted delivery of drugs and their controlled release. Nonetheless, entry into the target cells remains a challenge for many nanomaterials used for drug delivery. Moreover, cellular uptake limits the therapeutic efficiency of many anticancer drugs. An important goal is to increase the specific accumulation of these nanoparticles (NPs) at the desired cancerous tissues. Notably, cancer cells show a high demand for some amino acids and we have used this knowledge to develop novel carrier systems. In this study, drug carriers were produced by the conjugation of multiple amino acids such as l-histidine (H) and l-cysteine (C) or single amino acids such as only H with the G4.5 dendrimers (G) to produce GHC aggregates and GH NP carriers, respectively. Doxorubicin was loaded into the G4.5, GH, and GHC dendrimers (G/DOX, GH/DOX and GHC/DOX, respectively) and the release mechanism was demonstrated at pH 7.4 and pH 5.0. GH/DOX and GHC/DOX showed better stability under physiological conditions than the dendrimer alone (G/DOX). GH/DOX and GHC/DOX exhibited higher inhibition of HeLa cell proliferation in in vitro and in vivo studies in zebrafish, confirming the early release of DOX by disrupting the endosomal membrane and triggering the destabilization of carriers at a lower pH of 5.0.

Discovery of Novel Druggable Sites on Zika Virus NS3 Helicase Using X-ray Crystallography-Based Fragment Screening

The flavivirus family contains several important human pathogens, such as Zika virus (ZIKV), dengue, West Nile, and Yellow Fever viruses, that collectively lead to a large, global disease burden. Currently, there are no approved medicines that can target these viruses. The sudden outbreak of ZIKV infections in 2015⁻2016 posed a serious threat to global public health. While the epidemic has receded, persistent reservoirs of ZIKV infection can cause reemergence. Here, we have used X-ray crystallography-based screening to discover two novel sites on ZIKV NS3 helicase that can bind drug-like fragments. Both sites are structurally conserved in other flaviviruses, and mechanistically significant. The binding poses of four fragments, two for each of the binding sites, were characterized at atomic precision. Site A is a surface pocket on the NS3 helicase that is vital to its interaction with NS5 polymerase and formation of the flaviviral replication complex. Site B corresponds to a flexible, yet highly conserved, allosteric site at the intersection of the three NS3 helicase domains. Saturation transfer difference nuclear magnetic resonance (NMR) experiments were additionally used to evaluate the binding strength of the fragments, revealing dissociation constants (K_D) in the lower mM range. We conclude that the NS3 helicase of flaviviruses is a viable drug target. The data obtained open opportunities towards structure-based design of first-in-class anti-ZIKV compounds, as well as pan-flaviviral therapeutics.

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

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

Structural basis for the interaction of the beta-secretase with copper

The β-secretase (BACE1) features a unique sulfur rich motif (M₄₆₂xxxC₄₆₆xxxM₄₇₀xxxC₄₇₄xxxC₄₇₈) in its transmembrane helix (BACE1-TM) which is characteristic for proteins involved in copper ion storage and transport. While this motif has been shown to promote BACE1-TM trimerization and binding of copper ions in vitro, the structural basis for the interaction of copper ions with the BACE1-TM is still not well understood. Using molecular dynamics (MD) simulations, we show that membrane embedded BACE1-TMs adopt a flexible trimeric structure that binds and conducts copper ions through variable coordination. In coarse-grained (CG) MD simulations, the spontaneous assembly of BACE1-TMs trimers results in a right-handed helix packing arrangement. In subsequent atomistic MD simulations the sulfur rich motif defines characteristic copper ion coordination sites along a constricted partially solvated axial pore. Sliding and tilting of BACE1-TMs along smooth A₄₅₉xxxA_463/464xxA₄₆₇ surfaces, facilitated by a central P472 induced kink, enables copper ions to alternate between different coordination sites, including the prominent C466 and M470. We shed light into the structural arrangement of BACE1-TM trimers and propose a mechanism for copper ion conduction that might also apply to other proteins involved in metal ion transport.

Predicting structural and energetic changes in Met–aromatic motifs on methionine oxidation to the sulfoxide and sulfone

Methionine oxidation increases its affinity for aromatics in the gas phase but lowers it for most complexes in water.

Nanodisc-to-Nanofiber Transition of Noncovalent Peptide-Phospholipid Assemblies

We report a novel molecular architecture of peptide-phospholipid coassemblies. The amphiphilic peptide Ac-18A-NH₂ (18A), which was designed to mimic apolipoprotein α-helices, has been shown to form nanodisc structures with phospholipid bilayers. We show that an 18A peptide cysteine substitution at residue 11, 18A[A11C], forms fibrous assemblies with 1-palmitoyl-2-oleoyl-phosphatidylcholine at a lipid-to-peptide (L/P) molar ratio of 1, a fiber diameter of 10-20 nm, and a length of more than 1 μm. Furthermore, 18A[A11C] can form nanodiscs with these lipid bilayers at L/P ratios of 4-6. The peptide adopts α-helical structures in both the nanodisc and nanofiber assemblies, although the α-helical bundle structures were evident only in the nanofibers, and the phospholipids of the nanofibers were not lamellar. Fluorescence spectroscopic analysis revealed that the peptide and lipid molecules in the nanofibers exhibited different solvent accessibility and hydrophobicity from those of the nanodiscs. Furthermore, the cysteine substitution at residue 11 did not result in disulfide bond formation, although it was responsible for the nanofiber formation, suggesting that this free sulfhydryl group has an important functional role. Alternatively, the disulfide dimer of 18A[A11C] preferentially formed nanodiscs, even at an L/P ratio of 1. Interconversions of these discoidal and fibrous assemblies were induced by the stepwise addition of free 18A[A11C] or liposomes into the solution. Furthermore, these structural transitions could also be induced by the introduction of oxidative and reductive stresses to the assemblies. Our results demonstrate that heteromolecular lipid-peptide complexes represent a novel approach to the construction of controllable and functional nanoscale assemblies.

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

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

In situ Raman spectroscopy distinguishes between reversible and irreversible thiol modifications in l-cysteine

Reversibility of disulphide formation and breakage on l-cysteine examined through vibrational modes using in situ Raman spectroscopy.

Analysis of the interactions of sulfur-containing amino acids in membrane proteins

The interactions of Met and Cys with other amino acid side chains have received little attention, in contrast to aromatic-aromatic, aromatic-aliphatic or/and aliphatic-aliphatic interactions. Precisely, these are the only amino acids that contain a sulfur atom, which is highly polarizable and, thus, likely to participate in strong Van der Waals interactions. Analysis of the interactions present in membrane protein crystal structures, together with the characterization of their strength in small-molecule model systems at the ab-initio level, predicts that Met-Met interactions are stronger than Met-Cys ≈ Met-Phe ≈ Cys-Phe interactions, stronger than Phe-Phe ≈ Phe-Leu interactions, stronger than the Met-Leu interaction, and stronger than Leu-Leu ≈ Cys-Leu interactions. These results show that sulfur-containing amino acids form stronger interactions than aromatic or aliphatic amino acids. Thus, these amino acids may provide additional driving forces for maintaining the 3D structure of membrane proteins and may provide functional specificity.

Human Frataxin Folds Via an Intermediate State. Role of the C-Terminal Region

The aim of this study is to investigate the folding reaction of human frataxin, whose deficiency causes the neurodegenerative disease Friedreich's Ataxia (FRDA). The characterization of different conformational states would provide knowledge about how frataxin can be stabilized without altering its functionality. Wild-type human frataxin and a set of mutants, including two highly destabilized FRDA-associated variants were studied by urea-induced folding/unfolding in a rapid mixing device and followed by circular dichroism. The analysis clearly indicates the existence of an intermediate state (I) in the folding route with significant secondary structure content but relatively low compactness, compared with the native ensemble. However, at high NaCl concentrations I-state gains substantial compaction, and the unfolding barrier is strongly affected, revealing the importance of electrostatics in the folding mechanism. The role of the C-terminal region (CTR), the key determinant of frataxin stability, was also studied. Simulations consistently with experiments revealed that this stretch is essentially unstructured, in the most compact transition state ensemble (TSE2). The complete truncation of the CTR drastically destabilizes the native state without altering TSE2. Results presented here shed light on the folding mechanism of frataxin, opening the possibility of mutating it to generate hyperstable variants without altering their folding kinetics.

Molecular Mechanisms of Glutamine Synthetase Mutations that Lead to Clinically Relevant Pathologies

Glutamine synthetase (GS) catalyzes ATP-dependent ligation of ammonia and glutamate to glutamine. Two mutations of human GS (R324C and R341C) were connected to congenital glutamine deficiency with severe brain malformations resulting in neonatal death. Another GS mutation (R324S) was identified in a neurologically compromised patient. However, the molecular mechanisms underlying the impairment of GS activity by these mutations have remained elusive. Molecular dynamics simulations, free energy calculations, and rigidity analyses suggest that all three mutations influence the first step of GS catalytic cycle. The R324S and R324C mutations deteriorate GS catalytic activity due to loss of direct interactions with ATP. As to R324S, indirect, water-mediated interactions reduce this effect, which may explain the suggested higher GS residual activity. The R341C mutation weakens ATP binding by destabilizing the interacting residue R340 in the apo state of GS. Additionally, the mutation is predicted to result in a significant destabilization of helix H8, which should negatively affect glutamate binding. This prediction was tested in HEK293 cells overexpressing GS by dot-blot analysis: Structural stability of H8 was impaired through mutation of amino acids interacting with R341, as indicated by a loss of masking of an epitope in the glutamate binding pocket for a monoclonal anti-GS antibody by L-methionine-S-sulfoximine; in contrast, cells transfected with wild type GS showed the masking. Our analyses reveal complex molecular effects underlying impaired GS catalytic activity in three clinically relevant mutants. Our findings could stimulate the development of ATP binding-enhancing molecules by which the R324S mutant can be repaired extrinsically.

Glutamine synthetase (GS) catalyzes the ATP-dependent ligation of ammonia and glutamate to glutamine, which makes the enzyme essential for human nitrogen metabolism. Three mutations in human GS, R324C, R324S, and R341C, had been identified previously that lead to a glutamine deficiency, resulting in neonatal death in the case of R324C and R341C. However, the molecular mechanisms underlying this impairment of GS activity have remained elusive. Our results from computational biophysics approaches suggest that all three mutants influence the first step of GS’ catalytic cycle, namely ATP or glutamate binding. The analyses reveal a complex set of effects including the loss of direct interactions to substrates, the involvement of water-mediated interactions that alleviate part of the mutation effect, and long-range effects between the catalytic site and structural parts distant from it. As to the latter, experimental validation is in line with our prediction of a significant destabilization of helix H8 in the R341C mutant, which should negatively affect glutamate binding. Finally, our findings could stimulate the development of ATP-binding enhancing molecules for the R324S mutant, which has been suggested to have residual activity, that way extrinsically “repairing” the mutant.

Serine and Cysteine π-Interactions in Nature: A Comparison of the Frequency, Structure, and Stability of Contacts Involving Oxygen and Sulfur

Despite many DNA–protein π-interactions in high-resolution crystal structures, only four X–H···π or X···π interactions were found between serine (Ser) or cysteine (Cys) and DNA nucleobase π-systems in over 100 DNA–protein complexes (where X = O for Ser and X = S for Cys). Nevertheless, 126 non-covalent contacts occur between Ser or Cys and the aromatic amino acids in many binding arrangements within proteins. Furthermore, Ser and Cys protein–protein π-interactions occur with similar frequencies and strengths. Most importantly, due to the great stability that can be provided to biological macromolecules (up to –20 kJ mol–1 for neutral π-systems or –40 kJ mol–1 for cationic π-systems), Ser and Cys π-interactions should be considered when analyzing protein stability and function.

Intra-residue interactions in proteins: interplay between serine or cysteine side chains and backbone conformations, revealed by laser spectroscopy of isolated model peptides

Intra-residue interactions play an important role in proteins by influencing local folding of the backbone. Taking advantage of the capability of gas phase experiments to provide relevant information on the intrinsic H-bonding pattern of isolated peptide chains, the intra-residue interactions of serine and cysteine residues, i.e., OH/SH···OC(i) C6 and NH(i···)O/S C5 interactions in Ser/Cys residues, are probed by laser spectroscopy of isolated peptides. The strength of these local side chain-main chain interactions, elegantly documented from their IR spectral features for well-defined conformations of the main chain, demonstrates that a subtle competition exists between the two types of intra-residue bond: the C6 H-bond is the major interaction with Ser, in contrast to Cys where C5 interaction takes over. The restricted number of conformers observed in the gas phase experiment with Ser compared to Cys (where both extended and folded forms are observed) also suggests a significant mediation role of these intra-residue interactions on the competition between the several main chain folding patterns.

Crystal structure of HlyU, the hemolysin gene transcription activator, from Vibrio cholerae N16961 and functional implications

HlyU in Vibrio cholerae is known to be the transcriptional activator of the hemolysin gene, HlyA and possibly a regulator of other virulence factors influencing growth, colonization and pathogenicity of this infective agent. Here we report the crystal structure of HlyU from V. cholerae N16961 (HlyU_Vc) at 1.8Å. The protein, with five α-helices and three β-strands in the topology of α1-α2-β1-α3-α4-β2-β3-α5, forms a homodimer. Helices α3-α4 and a β sheet form the winged helix-turn-helix (wHTH) DNA-binding motif common to the transcription regulators of the SmtB/ArsR family. In spite of an overall fold similar to SmtB/ArsR family, it lacks any metal binding site seen in SmtB. A comparison of the dimeric interfaces showed that the one in SmtB is much larger and have salt bridges that can be disrupted to accommodate metal ions. A model of HlyU-DNA complex suggests bending of the DNA. Cys38 in the structure was found to be modified as sulfenic acid; the oxidized form was not seen in another structure solved under reducing condition. Although devoid of any metal binding site, the presence of a Cys residue exhibiting oxidation-reduction suggests the possibility of the existence of a redox switch in transcription regulation. A structure-based phylogenetic analysis of wHTH proteins revealed the segregation of metal and non-metal binding proteins as well as those in the latter group that are under redox control.

Crystal structure of the extended-spectrum β-lactamase PER-2 and insights into the role of specific residues in the interaction with β-lactams and β-lactamase inhibitors

PER-2 belongs to a small (7 members to date) group of extended-spectrum β-lactamases. It has 88% amino acid identity with PER-1 and both display high catalytic efficiencies toward most β-lactams. In this study, we determined the X-ray structure of PER-2 at 2.20 Å and evaluated the possible role of several residues in the structure and activity toward β-lactams and mechanism-based inhibitors. PER-2 is defined by the presence of a singular trans bond between residues 166 to 167, which generates an inverted Ω loop, an expanded fold of this domain that results in a wide active site cavity that allows for efficient hydrolysis of antibiotics like the oxyimino-cephalosporins, and a series of exclusive interactions between residues not frequently involved in the stabilization of the active site in other class A β-lactamases. PER β-lactamases might be included within a cluster of evolutionarily related enzymes harboring the conserved residues Asp136 and Asn179. Other signature residues that define these enzymes seem to be Gln69, Arg220, Thr237, and probably Arg/Lys240A ("A" indicates an insertion according to Ambler's scheme for residue numbering in PER β-lactamases), with structurally important roles in the stabilization of the active site and proper orientation of catalytic water molecules, among others. We propose, supported by simulated models of PER-2 in combination with different β-lactams, the presence of a hydrogen-bond network connecting Ser70-Gln69-water-Thr237-Arg220 that might be important for the proper activity and inhibition of the enzyme. Therefore, we expect that mutations occurring in these positions will have impacts on the overall hydrolytic behavior.

Stabilizing effect of electrostatic vs. aromatic interactions in diproline nucleated peptide β-hairpins

The contribution of Tyr-His vs. Cys-His interacting pairs to the scaffold stability of (D)Pro-(L)Pro nucleated peptide β-hairpins has been examined. We present direct evidence for the superiority of the Cys-His pairs, mediated by sulphur-imidazole interactions, as added stabilizing agents of the β-hairpin scaffold.

Weak hydrogen bonds formed by thiol groups in N-acetyl-(L)-cysteine and their response to the crystal structure distortion on increasing pressure

The effect of hydrostatic pressure on single crystals of N-acetyl-l-cysteine was followed at multiple pressure points from 10(-4) to 6.2 GPa with a pressure step of 0.2-0.3 GPa by Raman spectroscopy and X-ray diffraction. Since in the crystals of N-acetyl-l-cysteine the thiol group is involved in intermolecular hydrogen bonds not as a donor only (bonds S-H···O) but also as an acceptor (bonds N-H···S), increasing the pressure does not result in phase transitions. This makes a contrast with the polymorphs of l- and dl-cysteine, in which multiple phase transitions are observed already at relatively low hydrostatic pressures and are related to the changes in the conformation of the thiol side chains only weakly bound to the neighboring molecules in the structure and thus easily switching over the weak S-H···O and S-H···S hydrogen bonds. No phase transitions occur in N-acetyl-l-cysteine with increasing pressure, and changes in cell parameters and volume vs pressure do not reveal any peculiar features. Nevertheless, a more detailed analysis of the changes in intermolecular distances, in particular, of the geometric parameters of the hydrogen bonds based on X-ray single crystal diffraction analysis, complemented by an equally detailed study of the positions of all the significant bands in Raman spectra, allowed us to study the fine details of subtle changes in the hydrogen bond network. Thus, as pressure increases, a continuous shift of the hydrogen atom of the thiol group from one acceptor (a carboxyl group) to another acceptor (a carbonyl group) is observed. Precise single-crystal X-ray diffraction and polarized Raman spectroscopy structural data reveal the formation of a bifurcated S-H···O hydrogen bond with increasing pressure starting with ∼1.5 GPa. The analysis of the vibrational bands in Raman spectra has shown that different donor and acceptor groups start "feeling" the formation of the bifurcated S-H···O hydrogen bond in different pressure ranges. The results are discussed in relation to some of the previously published data on the effect of high pressure on the polymorphs of l-cysteine, dl-cysteine, and glutathione, that show similarity with the effects reported here for N-acetyl-l-cysteine. The results obtained in this work allow one to suggest new models for the pressure-induced structural rearrangements in the whole family of cysteine-containing crystals.

Evolution of the genetic code by incorporation of amino acids that improved or changed protein function

Fifty years have passed since the genetic code was deciphered, but how the genetic code came into being has not been satisfactorily addressed. It is now widely accepted that the earliest genetic code did not encode all 20 amino acids found in the universal genetic code as some amino acids have complex biosynthetic pathways and likely were not available from the environment. Therefore, the genetic code evolved as pathways for synthesis of new amino acids became available. One hypothesis proposes that early in the evolution of the genetic code four amino acids-valine, alanine, aspartic acid, and glycine-were coded by GNC codons (N = any base) with the remaining codons being nonsense codons. The other sixteen amino acids were subsequently added to the genetic code by changing nonsense codons into sense codons for these amino acids. Improvement in protein function is presumed to be the driving force behind the evolution of the code, but how improved function was achieved by adding amino acids has not been examined. Based on an analysis of amino acid function in proteins, an evolutionary mechanism for expansion of the genetic code is described in which individual coded amino acids were replaced by new amino acids that used nonsense codons differing by one base change from the sense codons previously used. The improved or altered protein function afforded by the changes in amino acid function provided the selective advantage underlying the expansion of the genetic code. Analysis of amino acid properties and functions explains why amino acids are found in their respective positions in the genetic code.

Molecular basis of antiangiogenic thrombospondin-1 type 1 repeat domain interactions with CD36

OBJECTIVE

Antiangiogenic activity of thrombospondin-1 and related proteins is mediated by interactions between thrombospondin type 1 repeat (TSR) domains and the CD36, LIMP-2, Emp sequence homology (CLESH) domain of the endothelial cell receptor CD36. We sought to characterize key molecular determinants of the interaction between thrombospondin-1 TSR domains and the CD36 CLESH domain.

APPROACH AND RESULTS

Recombinant thrombospondin-1 TSR2 and TSR(2,3) constructs inhibited microvascular endothelial cell migration, microvascular endothelial cell tube formation, and vessel sprouting in aortic ring assays. Interaction with CD36 CLESH decoy peptides negated these effects. Mutational analyses identified a cluster of residues that confer positive charge to the TSR2 surface and mediate interaction with CD36 CLESH. Antiangiogenic activity was significantly reduced by charge-neutralizing mutations of the Arg-Trp ladder in TSR2, but not TSR3. Additionally, I438 and K464 of TSR2 were shown to be required for CD36 CLESH binding to TSR2. A complementary acidic cluster within CD36 CLESH is also required for antiangiogenic activity.

CONCLUSIONS

Thrombospondin-1 interacts with CD36 CLESH through electrostatic interactions mediated by a positively charged TSR2 surface and multiple negatively charged CD36 CLESH residues. Two key residues serve as specificity determinants that identify other TSR domains that interact with CD36 CLESH.

The role of Cysteine 6.47 in class A GPCRs

Background

The CWxP motif of transmembrane helix 6 (x: any residue) is highly conserved in class A GPCRs. Within this motif, W6.48 is a big star in the theory of the global “toggle switch” because of its key role in the activation mechanism of GPCRs upon ligand binding. With all footlights focused on W6.48, the reason why the preceding residue, C6.47, is largely conserved is still unknown. The present study is aimed to fill up this lack of knowledge by characterizing the role of C6.47 of the CWxP motif.

Results

A complete analysis of available crystal structures has been made alongside with molecular dynamics simulations of model peptides to explore a possible structural role for C6.47.

Conclusions

We conclude that C6.47 does not modulate the conformation of the TM6 proline kink and propose that C6.47 participates in the rearrangement of the TM6 and TM7 interface accompanying activation.

Hyperdimensional analysis of amino acid pair distributions in proteins

Our manuscript presents a novel approach to protein structure analyses. We have organized an 8-dimensional data cube with protein 3D-structural information from 8706 high-resolution non-redundant protein-chains with the aim of identifying packing rules at the amino acid pair level. The cube contains information about amino acid type, solvent accessibility, spatial and sequence distance, secondary structure and sequence length. We are able to pose structural queries to the data cube using program ProPack. The response is a 1, 2 or 3D graph. Whereas the response is of a statistical nature, the user can obtain an instant list of all PDB-structures where such pair is found. The user may select a particular structure, which is displayed highlighting the pair in question. The user may pose millions of different queries and for each one he will receive the answer in a few seconds. In order to demonstrate the capabilities of the data cube as well as the programs, we have selected well known structural features, disulphide bridges and salt bridges, where we illustrate how the queries are posed, and how answers are given. Motifs involving cysteines such as disulphide bridges, zinc-fingers and iron-sulfur clusters are clearly identified and differentiated. ProPack also reveals that whereas pairs of Lys residues virtually never appear in close spatial proximity, pairs of Arg are abundant and appear at close spatial distance, contrasting the belief that electrostatic repulsion would prevent this juxtaposition and that Arg-Lys is perceived as a conservative mutation. The presented programs can find and visualize novel packing preferences in proteins structures allowing the user to unravel correlations between pairs of amino acids. The new tools allow the user to view statistical information and visualize instantly the structures that underpin the statistical information, which is far from trivial with most other SW tools for protein structure analysis.

Aromatic rings in chemical and biological recognition: energetics and structures

This review describes a multidimensional treatment of molecular recognition phenomena involving aromatic rings in chemical and biological systems. It summarizes new results reported since the appearance of an earlier review in 2003 in host-guest chemistry, biological affinity assays and biostructural analysis, data base mining in the Cambridge Structural Database (CSD) and the Protein Data Bank (PDB), and advanced computational studies. Topics addressed are arene-arene, perfluoroarene-arene, S⋅⋅⋅aromatic, cation-π, and anion-π interactions, as well as hydrogen bonding to π systems. The generated knowledge benefits, in particular, structure-based hit-to-lead development and lead optimization both in the pharmaceutical and in the crop protection industry. It equally facilitates the development of new advanced materials and supramolecular systems, and should inspire further utilization of interactions with aromatic rings to control the stereochemical outcome of synthetic transformations.