Post-translational modifications (PTMs) alter the proteome in response to intra- and extracellular signals, providing fundamental information processing in development, homeostasis and disease[1][1]–[5][2]. Here, we used proteome-wide structure predictions to investigate the structural context of protein modification sites and found nearly a fifth to be buried within proteins[6][3]. These cryptic sites within folded domains showed evolutionary conservation across species, as well as significantly less variability within the human population. We also found significantly more modification sites to be affected by recurrent cancer mutations and disease variants, with some variants appearing to mimic the modified or unmodified state, supporting their functionality and disease relevance. To understand how these buried sites become modified, we studied 23 phosphorylation sites across seven proteins using all-atom molecular dynamics simulations. The simulations revealed a range of mechanisms spanning broad timescales, ranging from picosecond fluctuations to seconds for proline isomerisation-gated phosphorylation sites. Notably, some sites remained fully buried in the native structure. Some conformational states with exposed cryptic sites resembled structures of co-translational protein folding intermediates formed on the ribosome[7][4], indicating an interplay between protein biogenesis, folding, modifications, and signalling. We propose more broadly that protein conformational fluctuations may serve as structure-specific regulatory mechanisms by kinetically burying or exposing the unmodified, and in some cases even modified, sidechains, thereby expanding the functional versatility of proteins and enabling signal transduction across diverse timescales. ### Competing Interest Statement The authors have declared no competing interest. [1]: #ref-1 [2]: #ref-5 [3]: #ref-6 [4]: #ref-7
Innate immune responses to cell damage-associated molecular patterns induce a controlled degree of inflammation, ideally avoiding the promotion of intense unwanted inflammatory adverse events. When released by damaged cells, Hsp70 can stimulate different responses that range from immune activation to immune suppression. The effects of Hsp70 are mediated through innate receptors expressed primarily by myeloid cells, such as dendritic cells (DCs). The regulatory innate receptors that bind to extracellular mouse Hsp70 (mHsp70) are not fully characterized, and neither are their potential interactions with activating innate receptors. Here, we show that extracellular mHsp70 interacts with a receptor complex formed by both inhibitory Siglec-E and activating lectin-like oxidized low-density lipoprotein receptor-1 (LOX-1) on DCs. We also find that this interaction takes place in lipid microdomains within the plasma membrane, and that Siglec-E acts as a negative regulator of LOX-1-mediated innate activation upon mHsp70 or oxidized LDL binding. Thus, Hsp70 can both bind to and modulate the interaction of inhibitory and activating innate receptors on the cell surface. These findings add another dimension of regulatory mechanism to indicate how self-molecules contribute to dampening of exacerbated inflammatory responses.
Accurate prediction of protein structures is essential for understanding biological functions and guiding biomedical research. However, maintaining synchronization between structure models and rapidly expanding, continuously evolving protein sequence databases remains a major challenge. Here, we present AlphaSync ( alphasync.stjude.org ), a comprehensive resource that complements the AlphaFold Protein Structure Database. AlphaSync currently provides 2.6 million UniProt-synchronized structural models, including predictions for 40,016 updated proteins and isoforms from 925 species. AlphaSync achieves complete, up-to-date proteome coverage for 42 species, including humans, key pathogens and model organisms. It also provides residue-level annotations such as solvent accessibility, dihedral angles, intrinsic disorder status and over 4.7 billion atom-level noncovalent contacts. Its up-to-date structural models and detailed annotations will facilitate the study of protein structure-function relationships, assessment of sequence variants and machine learning tasks including protein design. With an intuitive web interface and application programming interface, AlphaSync enables protein research at scale and in detail.
Biomolecular condensates mediate dynamic compartmentalization of cellular processes. The multivalent interactions that underlie biomolecular condensation are often promoted by intrinsically disordered regions (IDRs) within proteins. Although the role of IDRs in biomolecular condensates is well appreciated, predicting whether an IDR forms condensates in cells remains challenging. Here, we developed a machine learning model to accurately predict the condensation behavior of IDRs, analyzing 215 IDRs from fusion oncoproteins in human embryonic kidney (HEK) 293T cells. We identified distinct sequence-derived physicochemical features associated with condensation. Leveraging these data, our model predicts that ~12% of the ~13,000 IDRs in the human proteome are likely to form cellular condensates, establishing a robust framework for proteome-wide analysis of IDR-mediated biomolecular condensation. Notably, proteins with condensate-forming IDRs are significantly enriched in RNA processing and splicing functions and are predominantly localized to membraneless organelles, highlighting a central role of IDR-mediated biomolecular condensation in cellular organization and RNA biology.
Transcription factors containing disordered regions can now be mapped across the genome, aiding functional studies.
This study identified tumorigenic processes most dependent on murine heat shock protein 72 (HSP72) in the mouse mammary tumor virus-PyMT mammary tumor model, which give rise to spontaneous mammary tumors that exhibit HSP72-dependent metastasis to the lung. RNA-seq expression profiling of Hspa1a/Hspa1b (Hsp72) WT and Hsp72-/- primary mammary tumors discovered significantly lower expression of genes encoding components of the extracellular matrix (ECM) in Hsp72 knockout mammary tumors compared to WT controls. In vitro studies found that genetic or chemical inhibition of HSP72 activity in cultured collagen-expressing human or murine cells also reduces mRNA and protein levels of COL1A1 and several other ECM-encoding genes. In search of a possible mechanistic basis for this relationship, we found HSP72 to support the activation of the tumor growth factor-β-suppressor of mothers against decapentaplegic-3 signaling pathway and evidence of suppressor of mothers against decapentaplegic-3 and HSP72 coprecipitation, suggesting potential complex formation. Human COL1A1 mRNA expression was found to have prognostic value for HER2+ breast tumors over other breast cancer subtypes, suggesting a possible human disease context where targeting HSP72 may have a therapeutic rationale. Analysis of human HER2+ breast tumor gene expression data using a gene set comprising ECM-related gene and protein folding-related gene as an input to the statistical learning algorithm, Galgo, found a subset of these genes that can collectively stratify patients by relapse-free survival, further suggesting a potential interplay between the ECM and protein-folding genes may contribute to tumor progression.
We evaluate deconvolution methods, which infer levels of immune infiltration from bulk expression of tumor samples, through a community-wide DREAM Challenge. We assess six published and 22 community-contributed methods using in vitro and in silico transcriptional profiles of admixed cancer and healthy immune cells. Several published methods predict most cell types well, though they either were not trained to evaluate all functional CD8+ T cell states or do so with low accuracy. Several community-contributed methods address this gap, including a deep learning-based approach, whose strong performance establishes the applicability of this paradigm to deconvolution. Despite being developed largely using immune cells from healthy tissues, deconvolution methods predict levels of tumor-derived immune cells well. Our admixed and purified transcriptional profiles will be a valuable resource for developing deconvolution methods, including in response to common challenges we observe across methods, such as sensitive identification of functional CD4+ T cell states. Deconvolution methods infer levels of immune infiltration from bulk expression of tumour samples. Here, authors assess 6 published and 22 community-contributed methods via a DREAM Challenge using in vitro and in silico transcriptional profiles of admixed cancer and healthy immune cells.
Enhanced crosslinking and immunoprecipitation (eCLIP) sequencing is a powerful method for transcriptome-wide detection of binding sites of RNA-binding proteins (RBPs). However, identified crosslink sites can profoundly deviate from experimentally established functional elements of even well-studied RBPs. Current peak-calling strategies result in low replication and high false-positive rates. Here, we present the R/Bioconductor package DEWSeq that makes full use of replicate information and size-matched input controls. We benchmarked DEWSeq on 107 RBPs for which both eCLIP data and RNA sequence motifs are available and were able to more than double the number of motif-containing binding regions relative to standard eCLIP processing (2.3-fold median). The improvement not only relates to the number of binding sites (e.g., 3.1-fold of known motifs for RBFOX2), but also their subcellular localisation (e.g., 1.9-fold of mitochondrial genes for FASTKD2) and structural targets (e.g., 2.2-fold increase of stem-loop regions for SLBP). DEWSeq therefore shows promise as an improved processing method for eCLIP protein–RNA interaction data.
Recent advances in computational approaches and their integration into structural biology enable tackling increasingly complex questions. Here, we discuss several key areas, highlighting breakthroughs and remaining challenges. Theoretical modeling has provided tools to accurately predict and design protein structures on a scale currently difficult to achieve using experimental approaches. Molecular Dynamics simulations have become faster and more precise, delivering actionable information inaccessible by current experimental methods. Virtual screening workflows allow a high-throughput approach to discover ligands that bind and modulate protein function, while Machine Learning methods enable the design of proteins with new functionalities. Integrative structural biology combines several of these approaches, pushing the frontiers of structural and functional characterization to ever larger systems, advancing towards a complete understanding of the living cell. These breakthroughs will accelerate and significantly impact diverse areas of science.
The co-chaperone p50/Cdc37 is an important partner for Hsp90, assisting in molecular chaperone activities, particularly with regard to the regulation of protein kinases. Analysis of the structure of Hsp90-Cdc37-kinase complexes demonstrates the way in which Cdc37 interacts with and controls the folding of a large proportion of intracellular protein kinases. This co-chaperone thus stands at the hub of a multitude of intracellular signaling networks. Indeed, the influence of Cdc37 reaches beyond the housekeeping pathways of protein folding into the regulation of a wide range of cellular processes. This co-chaperone has attracted attention as a potential intermediate in carcinogenesis. Cdc37 is an attractive potential target in cancer due to (1) high expression in a number of tumor types and (2) control of multiple signaling pathways. These properties indicate (3) a potential for selectivity due to its elevated expression in malignant cells and (4) robustness, as the co-chaperone may control multiple growth signaling pathways and thus be less prone to evolution of resistance than less versatile oncoproteins. Cdc37 may also be involved in other aspects of pathophysiology and has been shown to be secreted in exosomes. Protein aggregation disorders have been linked to age-related declines in molecular chaperones and co-chaperones. Cdc37 also appears to be a potential agent in longevity due to its links to protein folding and autophagy, and it will be informative to study the role of Cdc37 maintenance/decline in aging organisms.
Fusion oncoproteins (FOs) arise from chromosomal translocations in ~17% of cancers and are often oncogenic drivers. Although some FOs can promote oncogenesis by undergoing liquid-liquid phase separation (LLPS) to form aberrant biomolecular condensates, the generality of this phenomenon is unknown. We explored this question by testing 166 FOs in HeLa cells and found that 58% formed condensates. The condensate-forming FOs displayed physicochemical features distinct from those of condensate-negative FOs and segregated into distinct feature-based groups that aligned with their sub-cellular localization and biological function. Using Machine Learning, we developed a predictor of FO condensation behavior, and discovered that 67% of ~3000 additional FOs likely form condensates, with 35% of those predicted to function by altering gene expression. 47% of the predicted condensate-negative FOs were associated with cell signaling functions, suggesting a functional dichotomy between condensate-positive and -negative FOs. Our Datasets and reagents are rich resources to interrogate FO condensation in the future.
Background: The IL-36 pathway plays a key role in the pathogenesis of generalized pustular psoriasis (GPP). In a proof-of-concept clinical trial, treatment with spesolimab, an anti-IL-36 receptor antibody, resulted in rapid skin and pustular clearance in patients presenting with GPP flares. Objective: We sought to compare the molecular profiles of lesional and nonlesional skin from patients with GPP or palmoplantar pustulosis (PPP) with skin from healthy volunteers, and to investigate the molecular changes after spesolimab treatment in the skin and blood of patients with GPP flares. Methods: Pre- and post-treatment skin and blood samples were collected from patients with GPP who participated in a single-arm, phase I study (n = 7). Skin biopsies from patients with PPP (n = 8) and healthy volunteers (n = 16) were obtained for comparison at baseline. Biomarkers were assessed by RNA-sequencing, histopathology, and immunohistochemistry. Results: In GPP and PPP lesions, 1287 transcripts were commonly upregulated or downregulated. Selected transcripts from the IL-36 signaling pathway were upregulated in untreated GPP and PPP lesions. In patients with GPP, IL-36 pathway-related signatures, T(H)1/T(H)17 and innate inflammation signaling, neutrophilic mediators, and keratinocyte-driven inflammation pathways were downregulated by spesolimab as early as week 1. Spesolimab also decreased related serum biomarkers and cell populations in the skin lesions from patients with GPP, including CD3(+) T, CD11c(+), and IL-36 gamma(+) cells and lipocalin-2-expressing cells. Conclusions: In patients with GPP, spesolimab showed rapid modulation of commonly dysregulated molecular pathways in GPP and PPP, which may be associated with improved clinical outcomes.
The zinc finger transcription factor EGR4 has previously been identified as having a critical role in the proliferation of small cell lung cancer. Here, we have identified a novel, shortened splice variant of this transcription factor (EGR4-S) that is regulated by Heat Shock Factor-1 (HSF1). Our findings demonstrate that the shortened variant (EGR4-S) is upregulated with high EGFR, HER2, and H-Rasv12-expressing breast cell lines, and its expression is inhibited in response to HER pathway inhibitors. Protein and mRNA analyses of HER2+ human breast tumours indicated the novel EGR4-S splice variant to be preferentially expressed in tumour tissue and not detectable in patient-matched normal tissue. Knockdown of EGR4-S in the HER2-amplified breast cancer cell line SKBR3 reduced cell growth, suggesting that EGR4-S supports the growth of HER2+ tumour cells. In addition to chemical inhibitors of the HER2 pathway, EGR4-S expression was also found to be suppressed by chemical stressors and the overexpression of HSF1. Under these conditions, reduced EGR4-S levels were associated with the observed lower cell growth rate, but the augmentation of properties associated with higher metastatic potential. Taken together, these findings identify EGR4-S as a potential biomarker for HER2 pathway activation in human tumours that is regulated by HSF1.
Breast cancer is a heterogeneous disease classified into four main subtypes with different clinical outcomes, such as patient survival, prognosis, and relapse. Current genetic tests for the differential diagnosis of BC subtypes showed a poor reproducibility. Therefore, an early and correct diagnosis of molecular subtypes is one of the challenges in the clinic. In the present study, we identified differentially expressed genes, long non-coding RNAs and RNA binding proteins for each BC subtype from a public dataset applying bioinformatics algorithms. In addition, we investigated their interactions and we proposed interacting biomarkers as potential signature specific for each BC subtype. We found a network of only 2 RBPs (RBM20 and PCDH20) and 2 genes (HOXB3 and RASSF7) for luminal A, a network of 21 RBPs and 53 genes for luminal B, a HER2-specific network of 14 RBPs and 30 genes, and a network of 54 RBPs and 302 genes for basal BC. We validated the signature considering their expression levels on an independent dataset evaluating their ability to classify the different molecular subtypes with a machine learning approach. Overall, we achieved good performances of classification with an accuracy >0.80. In addition, we found some interesting novel prognostic biomarkers such as RASSF7 for luminal A, DCTPP1 for luminal B, DHRS11, KLC3, NAGS, and TMEM98 for HER2, and ABHD14A and ADSSL1 for basal. The findings could provide preliminary evidence to identify putative new prognostic biomarkers and therapeutic targets for individual breast cancer subtypes.
Delivery of exogenous heat shock protein 90α (Hsp90α) and/or its induced expression in neural tissues has been suggested as a potential strategy to combat neurodegenerative disease. However, within a neurodegenerative context, a pro-inflammatory response to extracellular Hsp90α (eHsp90α) could undermine strategies to use it for therapeutic intervention. The aim of this study was to investigate the biological effects of eHsp90α on microglial cells, the primary mediators of inflammatory responses in the brain. Transcriptomic profiling by RNA-seq of primary microglia and the cultured EOC2 microglial cell line treated with eHsp90α showed the chaperone to stimulate activation of innate immune responses in microglia that were characterized by an increase in NF-kB-regulated genes. Further characterization showed this response to be substantially lower in amplitude than the effects of other inflammatory stimuli such as fibrillar amyloid-β (fAβ) or lipopolysaccharide (LPS). Additionally, the toxicity of conditioned media obtained from microglia treated with fAβ was attenuated by addition of eHsp90α. Using a co-culture system of microglia and hippocampal neuronal cell line HT22 cells separated by a chamber insert, the neurotoxicity of medium conditioned by microglia treated with fAβ was reduced when eHsp90α was also added. Mechanistically, eHsp90α was shown to activate Nrf2, a response which attenuated fAβ-induced nitric oxide production. The data thus suggested that eHsp90α protects against fAβ-induced oxidative stress. We also report eHsp90α to induce expression of macrophage receptor with collagenous structure (Marco), which would permit receptor-mediated endocytosis of fAβ.
The dynamic network of chaperone interactions known as the chaperome contributes significantly to the proteotoxic cell response and the malignant phenotype. To bypass the inherent redundancy in the network, we have used a microRNA (mir) approach to target multiple members of the chaperome simultaneously. We identified a potent microRNA, miR-570 that could bind the 3 ' untranslated regions of multiple HSP mRNAs and inhibit HSP synthesis. Transfection of cells with this miR species reduced expression of multiple HSPs, inhibited the heat shock response and reduced tumor cell growth while acted additively in combination with cytotoxic drugs. As overexpression of miR-570 elicited tumor suppressive effects, we inferred that this miR could play a potential role in inhibiting tumorigenesis and cancer cell growth. In accordance with this hypothesis, we determined a significant role for miR-570 in regulating markers of mammary tumor progression, including cell motility and invasion. Our data provide a proof of the principle that the tumor chaperome can be targeted by microRNAs suggesting a potential therapeutic avenue towards cancer therapy.
Heat shock protein 90 (Hsp90), although one of the most essential intracellular chaperones, can also play key roles in the extracellular milieu. Here, we review the properties of extracellular Hsp90 in cellular homeostasis in the heat shock response (HSR), focusing on cells of the central nervous system. Hsp90 can be secreted by microglia as well as other cell types by non-canonical pathways of secretion. The chaperone may then influence the behavior of distant cells and can for instance protect neuronal cells from the oxidative burst accompanying phagocytosis by microglia of beta-amyloid fibrils. A mechanism involving activation of the transcription factor Nrf2, and induction of the antioxidant response is reported. We review the potential role of extracellular Hsp90, Nrf2 and transcellular chaperone signaling in the non-cell-intrinsic HSR.
We have investigated the role of e xtracellular H eat s hock p rotein 90 alpha (eHsp90α) in conferring protection of neuronal cells against fibrillary amyloid beta (f-Aβ 1-42 ) toxicity mediated by microglial cells. Formation of f-Aβ 1-42 plaques leads to neurotoxic inflammation, a critical pathological feature of Alzheimer’s Disease. We observed increased uptake and clearance of internalized f-Aβ 1-42 by microglial cells treated with eHsp90α, an effect associated with activation of NRF2 (NF-E2-related factor 2) - mediated autophagy. eHsp90α thus mitigated the neuronal toxicity of f-Aβ 1-42 -activated microglia. In addition, eHsp90α facilitated f-Aβ 1-42 engulfment by microglial cells in vitro . In summary, eHsp90α triggers NRF2-mediated autophagy in microglia and thus protects against the neurotoxic effects of f-Aβ 1-42 .