A reference map of the human binary protein interactome

Curvature-enhanced graph convolutional network for biomolecular interaction prediction

Geometric deep learning has demonstrated a great potential in non-Euclidean data analysis. The incorporation of geometric insights into learning architecture is vital to its success. Here we propose a curvature-enhanced graph convolutional network (CGCN) for biomolecular interaction prediction. Our CGCN employs Ollivier-Ricci curvature (ORC) to characterize network local geometric properties and enhance the learning capability of GCNs. More specifically, ORCs are evaluated based on the local topology from node neighborhoods, and further incorporated into the weight function for the feature aggregation in message-passing procedure. Our CGCN model is extensively validated on fourteen real-world bimolecular interaction networks and analyzed in details using a series of well-designed simulated data. It has been found that our CGCN can achieve the state-of-the-art results. It outperforms all existing models, as far as we know, in thirteen out of the fourteen real-world datasets and ranks as the second in the rest one. The results from the simulated data show that our CGCN model is superior to the traditional GCN models regardless of the positive-to-negative-curvature ratios, network densities, and network sizes (when larger than 500).

RIscoper 2.0: A deep learning tool to extract RNA biomedical relation sentences from literature

RNA plays an extensive role in a multi-dimensional regulatory system, and its biomedical relationships are scattered across numerous biological studies. However, text mining works dedicated to the extraction of RNA biomedical relations remain limited. In this study, we established a comprehensive and reliable corpus of RNA biomedical relations, recruiting over 30,000 sentences manually curated from more than 15,000 biomedical literature. We also updated RIscoper 2.0, a BERT-based deep learning tool to extract RNA biomedical relation sentences from literature. Benefiting from approximately 100,000 annotated named entities, we integrated the text classification and named entity recognition tasks in this tool. Additionally, RIscoper 2.0 outperformed the original tool in both tasks and can discover new RNA biomedical relations. Additionally, we provided a user-friendly online search tool that enables rapid scanning of RNA biomedical relationships using local and online resources. Both the online tools and data resources of RIscoper 2.0 are available at http://www.rnainter.org/riscoper.

Link prediction in protein-protein interaction network: A similarity multiplied similarity algorithm with paths of length three

Protein-protein interactions (PPIs) are crucial for various biological processes, and predicting PPIs is a major challenge. To solve this issue, the most common method is link prediction. Currently, the link prediction methods based on network Paths of Length Three (L3) have been proven to be highly effective. In this paper, we propose a novel link prediction algorithm, named SMS, which is based on L3 and protein similarities. We first design a mixed similarity that combines the topological structure and attribute features of nodes. Then, we compute the predicted value by summing the product of all similarities along the L3. Furthermore, we propose the Max Similarity Multiplied Similarity (maxSMS) algorithm from the perspective of maximum impact. Our computational prediction results show that on six datasets, including S. cerevisiae, H. sapiens, and others, the maxSMS algorithm improves the precision of the top 500, area under the precision-recall curve, and normalized discounted cumulative gain by an average of 26.99%, 53.67%, and 6.7%, respectively, compared to other optimal methods.

Prioritizing drug targets by perturbing biological network response functions

Therapeutic interventions are designed to perturb the function of a biological system. However, there are many types of proteins that cannot be targeted with conventional small molecule drugs. Accordingly, many identified gene-regulatory drivers and downstream effectors are currently undruggable. Drivers and effectors are often connected by druggable signaling and regulatory intermediates. Methods to identify druggable intermediates therefore have general value in expanding the set of targets available for hypothesis-driven validation. Here we identify and prioritize potential druggable intermediates by developing a network perturbation theory, termed NetPert, for response functions of biological networks. Dynamics are defined by a network structure in which vertices represent genes and proteins, and edges represent gene-regulatory interactions and protein-protein interactions. Perturbation theory for network dynamics prioritizes targets that interfere with signaling from driver to response genes. Applications to organoid models for metastatic breast cancer demonstrate the ability of this mathematical framework to identify and prioritize druggable intermediates. While the short-time limit of the perturbation theory resembles betweenness centrality, NetPert is superior in generating target rankings that correlate with previous wet-lab assays and are more robust to incomplete or noisy network data. NetPert also performs better than a related graph diffusion approach. Wet-lab assays demonstrate that drugs for targets identified by NetPert, including targets that are not themselves differentially expressed, are active in suppressing additional metastatic phenotypes.

Democratizing protein language models with parameter-efficient fine-tuning

Proteomics has been revolutionized by large protein language models (PLMs), which learn unsupervised representations from large corpora of sequences. These models are typically fine-tuned in a supervised setting to adapt the model to specific downstream tasks. However, the computational and memory footprint of fine-tuning (FT) large PLMs presents a barrier for many research groups with limited computational resources. Natural language processing has seen a similar explosion in the size of models, where these challenges have been addressed by methods for parameter-efficient fine-tuning (PEFT). In this work, we introduce this paradigm to proteomics through leveraging the parameter-efficient method LoRA and training new models for two important tasks: predicting protein-protein interactions (PPIs) and predicting the symmetry of homooligomer quaternary structures. We show that these approaches are competitive with traditional FT while requiring reduced memory and substantially fewer parameters. We additionally show that for the PPI prediction task, training only the classification head also remains competitive with full FT, using five orders of magnitude fewer parameters, and that each of these methods outperform state-of-the-art PPI prediction methods with substantially reduced compute. We further perform a comprehensive evaluation of the hyperparameter space, demonstrate that PEFT of PLMs is robust to variations in these hyperparameters, and elucidate where best practices for PEFT in proteomics differ from those in natural language processing. All our model adaptation and evaluation code is available open-source at https://github.com/microsoft/peft_proteomics. Thus, we provide a blueprint to democratize the power of PLM adaptation to groups with limited computational resources.

Bidirectional linkage of DNA barcodes for the multiplexed mapping of higher-order protein interactions in cells

Capturing the full complexity of the diverse hierarchical interactions in the protein interactome is challenging. Here we report a DNA-barcoding method for the multiplexed mapping of pairwise and higher-order protein interactions and their dynamics within cells. The method leverages antibodies conjugated with barcoded DNA strands that can bidirectionally hybridize and covalently link to linearize closely spaced interactions within individual 3D protein complexes, encoding and decoding the protein constituents and the interactions among them. By mapping protein interactions in cancer cells and normal cells, we found that tumour cells exhibit a larger diversity and abundance of protein complexes with higher-order interactions. In biopsies of human breast-cancer tissue, the method accurately identified the cancer subtype and revealed that higher-order protein interactions are associated with cancer aggressiveness.

Structural aspects of iron‑sulfur protein biogenesis: An NMR view

Over the last decade, structural aspects involving iron‑sulfur (Fe/S) protein biogenesis have played an increasingly important role in understanding the high mechanistic complexity of mitochondrial and cytosolic machineries maturing Fe/S proteins. In this respect, solution NMR has had a significant impact because of its ability to monitor transient protein-protein interactions, which are abundant in the networks of pathways leading to Fe/S cluster biosynthesis and transfer, as well as thanks to the developments of paramagnetic NMR in both terms of new methodologies and accurate data interpretation. Here, we review the use of solution NMR in characterizing the structural aspects of human Fe/S proteins and their interactions in the framework of Fe/S protein biogenesis. We will first present a summary of the recent advances that have been achieved by paramagnetic NMR and then we will focus our attention on the role of solution NMR in the field of human Fe/S protein biogenesis.

Microscale measurements of protein complexes from single cells

Proteins execute numerous cell functions in concert with one another in protein-protein interactions (PPI). While essential in each cell, such interactions are not identical from cell to cell. Instead, PPI heterogeneity contributes to cellular phenotypic heterogeneity in health and diseases such as cancer. Understanding cellular phenotypic heterogeneity thus requires measurements of properties of PPIs such as abundance, stoichiometry, and kinetics at the single-cell level. Here, we review recent, exciting progress in single-cell PPI measurements. Novel technology in this area is enabled by microscale and microfluidic approaches that control analyte concentration in timescales needed to outpace PPI disassembly kinetics. We describe microscale innovations, needed technical capabilities, and methods poised to be adapted for single-cell analysis in the near future.

Evaluation of antioxidant, antimicrobial, and bioactive properties and peptide sequence composition of Malatya apricot kernels

BACKGROUND

This study used four different apricot (Prunus armeniaca) kernels cultivated in Malatya during two consecutive years. The varieties were Hacihaliloglu, Hasanbey, Kabaasi, and Zerdali. The physicochemical properties of the kernels were determined, and the bioactive content of the kernels was evaluated using kernel hydrolysates prepared using trypsin.

RESULTS

With regard to the physicochemical properties of the kernels, the dry matter ratio and protein content were the highest in the Hacihaliloglu variety; the ash ratio was the highest in the Kabaasi variety, and the free oil ratio was the highest in the Hasanbey variety. The bioactive compound content changed according to kernel variety. Angiotensin-converting enzyme inhibitors activity was found to be the highest in the Hacihaliloglu and Hasanbey varieties, which had the lowest amygdalin content, and Zerdali had the highest amygdalin content. The antioxidant and antimicrobial effects of the kernels varied, with Hasanbey and Kabaasi generally having the highest content in both analyses. Moreover, a concentration of 20 mg mL-1 of the hydrolysate was determined to have a destructive effect for the microorganisms used in this study. The storage protein of the kernels, except Hacihaliloglu, was found to be Prunin 1, with the longest matching protein chain in the kernels being R.QQQGGQLMANGLEETFCSLRLK.E.

CONCLUSION

The results suggest that the peptide sequences identified in the kernels could have antihypertensive, antioxidative, and Dipeptidyl peptidase IV (DPP-IV) inhibitory effects. Consequently, apricot kernels show potential for use in the production of functional food products. Of the kernels evaluated in this study, Hacihaliloglu and Hasanbey were deemed the most suitable varieties due to their higher bioactive content and lower amygdalin content. © 2024 The Author(s). Journal of The Science of Food and Agriculture published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Emerging paradigms and recent progress in targeting ErbB in cancers

The epidermal growth factor receptor (EGFR) family is a class of transmembrane proteins, highly regarded as anticancer targets due to their pivotal role in various malignancies. Standard cancer treatments targeting the ErbB receptors include tyrosine kinase inhibitors (TKIs) and monoclonal antibodies (mAbs). Despite their substantial survival benefits, the achievement of curative outcomes is hindered by acquired resistance. Recent advancements in anti-ErbB approaches, such as inhibitory peptides, nanobodies, targeted-protein degradation strategies, and bispecific antibodies (BsAbs), aim to overcome such resistance. More recently, emerging insights into the cell surface interactome of the ErbB family open new avenues for modulating ErbB signaling by targeting specific domains of ErbB partners. Here, we review recent progress in ErbB targeting and elucidate emerging paradigms that underscore the significance of EGF domain-containing proteins (EDCPs) as new ErbB-targeting pathways.

phuEGO: A Network-Based Method to Reconstruct Active Signaling Pathways From Phosphoproteomics Datasets

Signaling networks are critical for virtually all cell functions. Our current knowledge of cell signaling has been summarized in signaling pathway databases, which, while useful, are highly biased toward well-studied processes, and do not capture context specific network wiring or pathway cross-talk. Mass spectrometry-based phosphoproteomics data can provide a more unbiased view of active cell signaling processes in a given context, however, it suffers from low signal-to-noise ratio and poor reproducibility across experiments. While progress in methods to extract active signaling signatures from such data has been made, there are still limitations with respect to balancing bias and interpretability. Here we present phuEGO, which combines up-to-three-layer network propagation with ego network decomposition to provide small networks comprising active functional signaling modules. PhuEGO boosts the signal-to-noise ratio from global phosphoproteomics datasets, enriches the resulting networks for functional phosphosites and allows the improved comparison and integration across datasets. We applied phuEGO to five phosphoproteomics data sets from cell lines collected upon infection with SARS CoV2. PhuEGO was better able to identify common active functions across datasets and to point to a subnetwork enriched for known COVID-19 targets. Overall, phuEGO provides a flexible tool to the community for the improved functional interpretation of global phosphoproteomics datasets.

Protein family annotation for the Unified Human Gastrointestinal Proteome by DPCfam clustering

Technological advances in massively parallel sequencing have led to an exponential growth in the number of known protein sequences. Much of this growth originates from metagenomic projects producing new sequences from environmental and clinical samples. The Unified Human Gastrointestinal Proteome (UHGP) catalogue is one of the most relevant metagenomic datasets with applications ranging from medicine to biology. However, the low levels of sequence annotation may impair its usability. This work aims to produce a family classification of UHGP sequences to facilitate downstream structural and functional annotation. This is achieved through the release of the DPCfam-UHGP50 dataset containing 10,778 putative protein families generated using DPCfam clustering, an unsupervised pipeline grouping sequences into single or multi-domain architectures. DPCfam-UHGP50 considerably improves family coverage at protein and residue levels compared to the manually curated repository Pfam. In the hope that DPCfam-UHGP50 will foster future discoveries in the field of metagenomics of the human gut, we release a FAIR-compliant database of our results that is easily accessible via a searchable web server and Zenodo repository.

Predicting gene disease associations with knowledge graph embeddings for diseases with curtailed information

Knowledge graph embeddings (KGE) are a powerful technique used in the biomedical domain to represent biological knowledge in a low dimensional space. However, a deep understanding of these methods is still missing, and, in particular, regarding their applications to prioritize genes associated with complex diseases with reduced genetic information. In this contribution, we built a knowledge graph (KG) by integrating heterogeneous biomedical data and generated KGE by implementing state-of-the-art methods, and two novel algorithms: Dlemb and BioKG2vec. Extensive testing of the embeddings with unsupervised clustering and supervised methods showed that KGE can be successfully implemented to predict genes associated with diseases and that our novel approaches outperform most existing algorithms in both scenarios. Our findings underscore the significance of data quality, preprocessing, and integration in achieving accurate predictions. Additionally, we applied KGE to predict genes linked to Intervertebral Disc Degeneration (IDD) and illustrated that functions pertinent to the disease are enriched within the prioritized gene set.

TransNeT-CGP: A cluster-based comorbid gene prioritization by integrating transcriptomics and network-topological features

The local disruptions caused by the genes of one disease can influence the pathways associated with the other diseases resulting in comorbidity. For gene therapies, it is necessary to prioritize the key genes that regulate common biological mechanisms to tackle the issues caused by overlapping diseases. This work proposes a clustering-based computational approach for prioritising the comorbid genes within the overlapping disease modules by analyzing Protein-Protein Interaction networks. For this, a sub-network with gene interactions of the disease pair was extracted from the interactome. The edge weights are assigned by combining the pairwise gene expression correlation and betweenness centrality scores. Further, a weighted graph clustering algorithm is applied and dominant nodes of high-density clusters are ranked based on clustering coefficients and neighborhood connectivity. Case studies based on neurodegenerative diseases such as Amyotrophic Lateral Sclerosis- Spinal Muscular Atrophy (ALS-SMA) pair and cancers such as Ovarian Carcinoma-Invasive Ductal Breast Carcinoma (OC-IDBC) pair were conducted to examine the efficacy of the proposed method. To identify the mechanistic role of top-ranked genes, we used Functional and Pathway enrichment analysis, connectivity analysis with leave-one-out (LOO) method, analysis of associated disease-related protein complexes, and prioritization tools such as TOPPGENE and Heml2.0. From pathway analysis, it was observed that the top 10 genes obtained using the proposed method were associated with 10 pathways in ALS-SMA comorbidity and 15 in the case of OC-IDBC, while that in similar methods like SAPDSB and S2B were 4, 6 respectively for ALS-SMA and 9, 10 respectively for OC-IDBC. In both case studies, 70 % of the disease-specific benchmark protein complexes were linked to top-ranked genes of the proposed method while that of SAPDSB and S2B were 55 % and 60 % respectively. Additionally, it was found that the removal of the top 10 genes disconnect the network into 14 distinct components in the case of ALS-SMA and 9 in the case of OC-IDBC. The experimental results shows that the proposed method can be effectively used for identifying key genes in comorbidity and can offer insights about the intricate molecular relationship driving comorbid diseases.

Ancient eukaryotic protein interactions illuminate modern genetic traits and disorders

All eukaryotes share a common ancestor from roughly 1.5 - 1.8 billion years ago, a single-celled, swimming microbe known as LECA, the Last Eukaryotic Common Ancestor. Nearly half of the genes in modern eukaryotes were present in LECA, and many current genetic diseases and traits stem from these ancient molecular systems. To better understand these systems, we compared genes across modern organisms and identified a core set of 10,092 shared protein-coding gene families likely present in LECA, a quarter of which are uncharacterized. We then integrated >26,000 mass spectrometry proteomics analyses from 31 species to infer how these proteins interact in higher-order complexes. The resulting interactome describes the biochemical organization of LECA, revealing both known and new assemblies. We analyzed these ancient protein interactions to find new human gene-disease relationships for bone density and congenital birth defects, demonstrating the value of ancestral protein interactions for guiding functional genetics today.

Structural biology in cellulo: Minding the gap between conceptualization and realization

Recent technological advances have deepened our perception of cellular structure. However, most structural data doesn't originate from intact cells, limiting our understanding of cellular processes. Here, we discuss current and future developments that will bring us towards a structural picture of the cell. Electron cryotomography is the standard bearer, with its ability to provide in cellulo snapshots. Single-particle electron microscopy (of purified biomolecules and of complex mixtures) and covalent crosslinking combined with mass spectrometry also have significant roles to play, as do artificial intelligence algorithms in their many guises. To integrate these multiple approaches, data curation and standardisation will be critical - as is the need to expand efforts beyond our current protein-centric view to the other (macro)molecules that sustain life.

SLiMAn 2.0: meaningful navigation through peptide-protein interaction networks

Among the myriad of protein-protein interactions occurring in living organisms, a substantial amount involves small linear motifs (SLiMs) recognized by structured domains. However, predictions of SLiM-based networks are tedious, due to the abundance of such motifs and a high portion of false positive hits. For this reason, a webserver SLiMAn (Short Linear Motif Analysis) was developed to focus the search on the most relevant SLiMs. Using SLiMAn, one can navigate into a given (meta-)interactome and tune a variety of parameters associated to each type of SLiMs in attempt to identify functional ELM motifs and their recognition domains. The IntAct and BioGRID databases bring experimental information, while IUPred and AlphaFold provide boundaries of folded and disordered regions. Post-translational modifications listed in PhosphoSite+ are highlighted. Links to PubMed accelerate scrutiny into the literature, to support (or not) putative pairings. Dedicated visualization features are also incorporated, such as Cytoscape for macromolecular networks and BINANA for intermolecular contacts within structural models generated by SCWRL 3.0. The use of SLiMAn 2.0 is illustrated on a simple example. It is freely available at https://sliman2.cbs.cnrs.fr.

Protein codes promote selective subcellular compartmentalization

Cells have evolved mechanisms to distribute ~10 billion protein molecules to subcellular compartments where diverse proteins involved in shared functions must efficiently assemble. Here, we demonstrate that proteins with shared functions share amino acid sequence codes that guide them to compartment destinations. A protein language model, ProtGPS, was developed that predicts with high performance the compartment localization of human proteins excluded from the training set. ProtGPS successfully guided generation of novel protein sequences that selectively assemble in targeted subcellular compartments. ProtGPS also identified pathological mutations that change this code and lead to altered subcellular localization of proteins. Our results indicate that protein sequences contain not only a folding code, but also a previously unrecognized code governing their distribution in specific cellular compartments.

ZEPPI: Proteome-scale sequence-based evaluation of protein-protein interaction models

We introduce ZEPPI (Z-score Evaluation of Protein-Protein Interfaces), a framework to evaluate structural models of a complex based on sequence coevolution and conservation involving residues in protein-protein interfaces. The ZEPPI score is calculated by comparing metrics for an interface to those obtained from randomly chosen residues. Since contacting residues are defined by the structural model, this obviates the need to account for indirect interactions. Further, although ZEPPI relies on species-paired multiple sequence alignments, its focus on interfacial residues allows it to leverage quite shallow alignments. ZEPPI can be implemented on a proteome-wide scale and is applied here to millions of structural models of dimeric complexes in the Escherichia coli and human interactomes found in the PrePPI database. PrePPI's scoring function is based primarily on the evaluation of protein-protein interfaces, and ZEPPI adds a new feature to this analysis through the incorporation of evolutionary information. ZEPPI performance is evaluated through applications to experimentally determined complexes and to decoys from the CASP-CAPRI experiment. As we discuss, the standard CAPRI scores used to evaluate docking models are based on model quality and not on the ability to give yes/no answers as to whether two proteins interact. ZEPPI is able to detect weak signals from PPI models that the CAPRI scores define as incorrect and, similarly, to identify potential PPIs defined as low confidence by the current PrePPI scoring function. A number of examples that illustrate how the combination of PrePPI and ZEPPI can yield functional hypotheses are provided.

USP7 interacts with and destabilizes oncoprotein SET

Oncoprotein SE translocation (SET) is frequently overexpressed in different types of tumors and correlated with poor prognosis of cancer patients. Targeting SET has been considered a promising strategy for cancer intervention. However, the mechanisms by which SET is regulated under cellular conditions are largely unknown. Here, by performing a tandem affinity purification-mass spectrometry (TAP-MS), we identify that the ubiquitin-specific protease 7 (USP7) forms a stable protein complex with SET in cancer cells. Further analyses reveal that the acidic domain of SET directly binds USP7 while both catalytic domain and ubiquitin-like (UBL) domains of USP7 are required for SET binding. Knockdown of USP7 has no effect on the mRNA level of SET. However, we surprisingly find that USP7 depletion leads to a dramatic elevation of SET protein levels, suggesting that USP7 plays a key role in destabilizing oncoprotein SET, possibly through an indirect mechanism. To our knowledge, our data report the first deubiquitinase (DUB) that physically associates with oncoprotein SET and imply an unexpected regulatory effect of USP7 on SET stability.

System-level analysis of genes mutated in muscular dystrophies reveals a functional pattern associated with muscle weakness distribution

Muscular dystrophies (MDs) are inherited genetic diseases causing weakness and degeneration of muscles. The distribution of muscle weakness differs between MDs, involving distal muscles or proximal muscles. While the mutations in most of the MD-associated genes lead to either distal or proximal onset, there are also genes whose mutations can cause both types of onsets. We hypothesized that the genes associated with different MD onsets code proteins with distinct cellular functions. To investigate this, we collected the MD-associated genes and assigned them to three onset groups: genes mutated only in distal onset dystrophies, genes mutated only in proximal onset dystrophies, and genes mutated in both types of onsets. We then systematically evaluated the cellular functions of these gene sets with computational strategies based on functional enrichment analysis and biological network analysis. Our analyses demonstrate that genes mutated in either distal or proximal onset MDs code proteins linked with two distinct sets of cellular processes. Interestingly, these two sets of cellular processes are relevant for the genes that are associated with both onsets. Moreover, the genes associated with both onsets display high centrality and connectivity in the network of muscular dystrophy genes. Our findings support the hypothesis that the proteins associated with distal or proximal onsets have distinct functional characteristics, whereas the proteins associated with both onsets are multifunctional.

The lack of Tex44 causes severe subfertility with flagellar abnormalities in male mice

By analyzing a mouse Interspecific Recombinant Congenic Strain (IRCS), we previously identified a quantitative trait locus (QTL), called Mafq1 on mouse chromosome 1, that is associated with male hypofertility and ultrastructural sperm abnormalities. Within this locus, we identified a new candidate gene that could be implicated in a reproductive phenotype: Tex44 (Testis-expressed protein 44). We thus performed a CRISPR/Cas9-mediated complete deletion of this gene in mice in order to study its function. Tex44-KO males were severely hypofertile in vivo and in vitro due to a drastic reduction of sperm motility which itself resulted from important morphological sperm abnormalities. Namely, Tex44-KO sperm showed a disorganized junction between the midpiece and the principal piece of the flagellum, leading to a 180° flagellar bending in this region. In addition, the loss of some axonemal microtubule doublets and outer dense fibers in the flagellum's principal piece has been observed. Our results suggest that, in mice, TEX44 is implicated in the correct set-up of the sperm flagellum during spermiogenesis and its absence leads to flagellar abnormalities and consequently to severe male hypofertility.

Single-cell analysis reveals an antiviral network that controls Zika virus infection in human dendritic cells

Zika virus (ZIKV) is a mosquito-borne flavivirus that caused an epidemic in the Americas in 2016 and is linked to severe neonatal birth defects, including microcephaly and spontaneous abortion. To better understand the host response to ZIKV infection, we adapted the 10× Genomics Chromium single-cell RNA sequencing (scRNA-seq) assay to simultaneously capture viral RNA and host mRNA. Using this assay, we profiled the antiviral landscape in a population of human monocyte-derived dendritic cells infected with ZIKV at the single-cell level. The bystander cells, which lacked detectable viral RNA, expressed an antiviral state that was enriched for genes coinciding predominantly with a type I interferon (IFN) response. Within the infected cells, viral RNA negatively correlated with type I IFN-dependent and -independent genes (the antiviral module). We modeled the ZIKV-specific antiviral state at the protein level, leveraging experimentally derived protein interaction data. We identified a highly interconnected network between the antiviral module and other host proteins. In this work, we propose a new paradigm for evaluating the antiviral response to a specific virus, combining an unbiased list of genes that highly correlate with viral RNA on a per-cell basis with experimental protein interaction data.

IMPORTANCE

Zika virus (ZIKV) remains a public health threat given its potential for re-emergence and the detrimental fetal outcomes associated with infection during pregnancy. Understanding the dynamics between ZIKV and its host is critical to understanding ZIKV pathogenesis. Through ZIKV-inclusive single-cell RNA sequencing (scRNA-seq), we demonstrate on the single-cell level the dynamic interplay between ZIKV and the host: the transcriptional program that restricts viral infection and ZIKV-mediated inhibition of that response. Our ZIKV-inclusive scRNA-seq assay will serve as a useful tool for gaining greater insight into the host response to ZIKV and can be applied more broadly to the flavivirus field.

Mapping protein-protein interactions by mass spectrometry

Protein-protein interactions (PPIs) are essential for numerous biological activities, including signal transduction, transcription control, and metabolism. They play a pivotal role in the organization and function of the proteome, and their perturbation is associated with various diseases, such as cancer, neurodegeneration, and infectious diseases. Recent advances in mass spectrometry (MS)-based protein interactomics have significantly expanded our understanding of the PPIs in cells, with techniques that continue to improve in terms of sensitivity, and specificity providing new opportunities for the study of PPIs in diverse biological systems. These techniques differ depending on the type of interaction being studied, with each approach having its set of advantages, disadvantages, and applicability. This review highlights recent advances in enrichment methodologies for interactomes before MS analysis and compares their unique features and specifications. It emphasizes prospects for further improvement and their potential applications in advancing our knowledge of PPIs in various biological contexts.

NetREm: Network Regression Embeddings reveal cell-type transcription factor coordination for gene regulation

Transcription factor (TF) coordination plays a key role in target gene (TG) regulation via protein-protein interactions (PPIs) and DNA co-binding to regulatory elements. Single-cell technologies facilitate gene expression measurement for individual cells and cell-type identification, yet the connection between TF coordination and TG regulation of various cell types remains unclear. To address this, we have developed a novel computational approach, Network Regression Embeddings (NetREm), to reveal cell-type TF-TF coordination activities for TG regulation. NetREm leverages network-constrained regularization using prior knowledge of direct and/or indirect PPIs among TFs to analyze single-cell gene expression data. We test NetREm by simulation data and benchmark its performance in 4 real-world applications that have gold standard TF-TG networks available: mouse (mESCs) and simulated human (hESCs) embryonic stem (ESCs), human hematopoietic stem (HSCs), and mouse dendritic (mDCs) cells. Further, we use NetREm to prioritize valid novel TF-TF coordination links in human Peripheral Blood Mononuclear cell (PBMC) sub-types. We apply NetREm to analyze various cell types in both central (CNS) and peripheral (PNS) nerve system (NS) (e.g. neuronal, glial, Schwann cells (SCs)) as well as in Alzheimers disease (AD). Our findings uncover cell-type coordinating TFs and identify new TF-TG candidate links. We validate our top predictions using Cut&Run and knockout loss-of-function expression data in rat/mouse models and compare results with additional functional genomic data, including expression quantitative trait loci (eQTL) and Genome-Wide Association Studies (GWAS) to link genetic variants (single nucleotide polymorphisms (SNPs)) to TF coordination.

SUMOylated Golgin45 associates with PML-NB to transcriptionally regulate lipid metabolism genes during heat shock stress

Golgin tethers are known to mediate vesicular transport in the secretory pathway, whereas it is relatively unknown whether they may mediate cellular stress response within the cell. Here, we describe a cellular stress response during heat shock stress via SUMOylation of a Golgin tether, Golgin45. We found that Golgin45 is a SUMOylated Golgin via SUMO1 under steady state condition. Upon heat shock stress, the Golgin enters the nucleus by interacting with Importin-β2 and gets further modified by SUMO3. Importantly, SUMOylated Golgin45 appears to interact with PML and SUMO-deficient Golgin45 mutant functions as a dominant negative for PML-NB formation during heat shock stress, suppressing transcription of lipid metabolism genes. These results indicate that Golgin45 may play a role in heat stress response by transcriptional regulation of lipid metabolism genes in SUMOylation-dependent fashion.

Functional Genomics and Insights into the Pathogenesis and Treatment of Psoriasis

Psoriasis is a lifelong, systemic, immune mediated inflammatory skin condition, affecting 1-3% of the world's population, with an impact on quality of life similar to diseases like cancer or diabetes. Genetics are the single largest risk factor in psoriasis, with Genome-Wide Association (GWAS) studies showing that many psoriasis risk genes lie along the IL-23/Th17 axis. Potential psoriasis risk genes determined through GWAS can be annotated and characterised using functional genomics, allowing the identification of novel drug targets and the repurposing of existing drugs. This review is focused on the IL-23/Th17 axis, providing an insight into key cell types, cytokines, and intracellular signaling pathways involved. This includes examination of currently available biological treatments, time to relapse post drug withdrawal, and rates of primary/secondary drug failure, showing the need for greater understanding of the underlying genetic mechanisms of psoriasis and how they can impact treatment. This could allow for patient stratification towards the treatment most likely to reduce the burden of disease for the longest period possible.

A UL26-PIAS1 complex antagonizes anti-viral gene expression during Human Cytomegalovirus infection

Viral disruption of innate immune signaling is a critical determinant of productive infection. The Human Cytomegalovirus (HCMV) UL26 protein prevents anti-viral gene expression during infection, yet the mechanisms involved are unclear. We used TurboID-driven proximity proteomics to identify putative UL26 interacting proteins during infection to address this issue. We find that UL26 forms a complex with several immuno-regulatory proteins, including several STAT family members and various PIAS proteins, a family of E3 SUMO ligases. Our results indicate that UL26 prevents STAT phosphorylation during infection and antagonizes transcriptional activation induced by either interferon α (IFNA) or tumor necrosis factor α (TNFα). Additionally, we find that the inactivation of PIAS1 sensitizes cells to inflammatory stimulation, resulting in an anti-viral transcriptional environment similar to ΔUL26 infection. Further, PIAS1 is important for HCMV cell-to-cell spread, which depends on the presence of UL26, suggesting that the UL26-PIAS1 interaction is vital for modulating intrinsic anti-viral defense.

Identification of a chromatin-bound ERRα interactome network in mouse liver

OBJECTIVES

Estrogen-related-receptor α (ERRα) plays a critical role in the transcriptional regulation of cellular bioenergetics and metabolism, and perturbations in its activity have been associated with metabolic diseases. While several coactivators and corepressors of ERRα have been identified to date, a knowledge gap remains in understanding the extent to which ERRα cooperates with coregulators in the control of gene expression. Herein, we mapped the primary chromatin-bound ERRα interactome in mouse liver.

METHODS

RIME (Rapid Immuno-precipitation Mass spectrometry of Endogenous proteins) analysis using mouse liver samples from two circadian time points was used to catalog ERRα-interacting proteins on chromatin. The genomic crosstalk between ERRα and its identified cofactors in the transcriptional control of precise gene programs was explored through cross-examination of genome-wide binding profiles from chromatin immunoprecipitation-sequencing (ChIP-seq) studies. The dynamic interplay between ERRα and its newly uncovered cofactor Host cell factor C1 (HCFC1) was further investigated by loss-of-function studies in hepatocytes.

RESULTS

Characterization of the hepatic ERRα chromatin interactome led to the identification of 48 transcriptional interactors of which 42 were previously unknown including HCFC1. Interrogation of available ChIP-seq binding profiles highlighted oxidative phosphorylation (OXPHOS) under the control of a complex regulatory network between ERRα and multiple cofactors. While ERRα and HCFC1 were found to bind to a large set of common genes, only a small fraction showed their colocalization, found predominately near the transcriptional start sites of genes particularly enriched for components of the mitochondrial respiratory chain. Knockdown studies demonstrated inverse regulatory actions of ERRα and HCFC1 on OXPHOS gene expression ultimately dictating the impact of their loss-of-function on mitochondrial respiration.

CONCLUSIONS

Our work unveils a repertoire of previously unknown transcriptional partners of ERRα comprised of chromatin modifiers and transcription factors thus advancing our knowledge of how ERRα regulates metabolic transcriptional programs.

State of the Interactomes: an evaluation of molecular networks for generating biological insights

Advancements in genomic and proteomic technologies have powered the use of gene and protein networks ("interactomes") for understanding genotype-phenotype translation. However, the proliferation of interactomes complicates the selection of networks for specific applications. Here, we present a comprehensive evaluation of 46 current human interactomes, encompassing protein-protein interactions as well as gene regulatory, signaling, colocalization, and genetic interaction networks. Our analysis shows that large composite networks such as HumanNet, STRING, and FunCoup are most effective for identifying disease genes, while smaller networks such as DIP and SIGNOR demonstrate strong interaction prediction performance. These findings provide a benchmark for interactomes across diverse network biology applications and clarify factors that influence network performance. Furthermore, our evaluation pipeline paves the way for continued assessment of emerging and updated interaction networks in the future.

Evolution of chromosome-arm aberrations in breast cancer through genetic network rewiring

The basal breast cancer subtype is enriched for triple-negative breast cancer (TNBC) and displays consistent large chromosomal deletions. Here, we characterize evolution and maintenance of chromosome 4p (chr4p) loss in basal breast cancer. Analysis of The Cancer Genome Atlas data shows recurrent deletion of chr4p in basal breast cancer. Phylogenetic analysis of a panel of 23 primary tumor/patient-derived xenograft basal breast cancers reveals early evolution of chr4p deletion. Mechanistically we show that chr4p loss is associated with enhanced proliferation. Gene function studies identify an unknown gene, C4orf19, within chr4p, which suppresses proliferation when overexpressed-a member of the PDCD10-GCKIII kinase module we name PGCKA1. Genome-wide pooled overexpression screens using a barcoded library of human open reading frames identify chromosomal regions, including chr4p, that suppress proliferation when overexpressed in a context-dependent manner, implicating network interactions. Together, these results shed light on the early emergence of complex aneuploid karyotypes involving chr4p and adaptive landscapes shaping breast cancer genomes.

PBX/Knotted 1 homeobox-2 (PKNOX2) is a novel regulator of myocardial fibrosis

Much effort has been made to uncover the cellular heterogeneities of human hearts by single-nucleus RNA sequencing. However, the cardiac transcriptional regulation networks have not been systematically described because of the limitations in detecting transcription factors. In this study, we optimized a pipeline for isolating nuclei and conducting single-nucleus RNA sequencing targeted to detect a higher number of cell signal genes and an optimal number of transcription factors. With this unbiased protocol, we characterized the cellular composition of healthy human hearts and investigated the transcriptional regulation networks involved in determining the cellular identities and functions of the main cardiac cell subtypes. Particularly in fibroblasts, a novel regulator, PKNOX2, was identified as being associated with physiological fibroblast activation in healthy hearts. To validate the roles of these transcription factors in maintaining homeostasis, we used single-nucleus RNA-sequencing analysis of transplanted failing hearts focusing on fibroblast remodelling. The trajectory analysis suggested that PKNOX2 was abnormally decreased from fibroblast activation to pathological myofibroblast formation. Both gain- and loss-of-function in vitro experiments demonstrated the inhibitory role of PKNOX2 in pathological fibrosis remodelling. Moreover, fibroblast-specific overexpression and knockout of PKNOX2 in a heart failure mouse model induced by transverse aortic constriction surgery significantly improved and aggravated myocardial fibrosis, respectively. In summary, this study established a high-quality pipeline for single-nucleus RNA-sequencing analysis of heart muscle. With this optimized protocol, we described the transcriptional regulation networks of the main cardiac cell subtypes and identified PKNOX2 as a novel regulator in suppressing fibrosis and a potential therapeutic target for future translational studies.

BTB domain mutations perturbing KCTD15 oligomerisation cause a distinctive frontonasal dysplasia syndrome

INTRODUCTION

KCTD15 encodes an oligomeric BTB domain protein reported to inhibit neural crest formation through repression of Wnt/beta-catenin signalling, as well as transactivation by TFAP2. Heterozygous missense variants in the closely related paralogue KCTD1 cause scalp-ear-nipple syndrome.

METHODS

Exome sequencing was performed on a two-generation family affected by a distinctive phenotype comprising a lipomatous frontonasal malformation, anosmia, cutis aplasia of the scalp and/or sparse hair, and congenital heart disease. Identification of a de novo missense substitution within KCTD15 led to targeted sequencing of DNA from a similarly affected sporadic patient, revealing a different missense mutation. Structural and biophysical analyses were performed to assess the effects of both amino acid substitutions on the KCTD15 protein.

RESULTS

A heterozygous c.310G>C variant encoding p.(Asp104His) within the BTB domain of KCTD15 was identified in an affected father and daughter and segregated with the phenotype. In the sporadically affected patient, a de novo heterozygous c.263G>A variant encoding p.(Gly88Asp) was present in KCTD15. Both substitutions were found to perturb the pentameric assembly of the BTB domain. A crystal structure of the BTB domain variant p.(Gly88Asp) revealed a closed hexameric assembly, whereas biophysical analyses showed that the p.(Asp104His) substitution resulted in a monomeric BTB domain likely to be partially unfolded at physiological temperatures.

CONCLUSION

BTB domain substitutions in KCTD1 and KCTD15 cause clinically overlapping phenotypes involving craniofacial abnormalities and cutis aplasia. The structural analyses demonstrate that missense substitutions act through a dominant negative mechanism by disrupting the higher order structure of the KCTD15 protein complex.

PheSeq, a Bayesian deep learning model to enhance and interpret the gene-disease association studies

Despite the abundance of genotype-phenotype association studies, the resulting association outcomes often lack robustness and interpretations. To address these challenges, we introduce PheSeq, a Bayesian deep learning model that enhances and interprets association studies through the integration and perception of phenotype descriptions. By implementing the PheSeq model in three case studies on Alzheimer's disease, breast cancer, and lung cancer, we identify 1024 priority genes for Alzheimer's disease and 818 and 566 genes for breast cancer and lung cancer, respectively. Benefiting from data fusion, these findings represent moderate positive rates, high recall rates, and interpretation in gene-disease association studies.

Golgi-associated retrograde protein (GARP) complex-dependent endosomes to trans Golgi network retrograde trafficking is controlled by Rab4b

BACKGROUND

The trafficking of cargoes from endosomes to the trans-Golgi network requires numerous sequential and coordinated steps. Cargoes are sorted into endosomal-derived carriers that are transported, tethered, and fused to the trans-Golgi network. The tethering step requires several complexes, including the Golgi-associated retrograde protein complex, whose localization at the trans-Golgi network is determined by the activity of small GTPases of the Arl and Rab family. However, how the Golgi-associated retrograde protein complex recognizes the endosome-derived carriers that will fuse with the trans-Golgi network is still unknown.

METHODS

We studied the retrograde trafficking to the trans-Golgi network by using fluorescent cargoes in cells overexpressing Rab4b or after Rab4b knocked-down by small interfering RNA in combination with the downregulation of subunits of the Golgi-associated retrograde protein complex. We used immunofluorescence and image processing (Super Resolution Radial Fluctuation and 3D reconstruction) as well as biochemical approaches to characterize the consequences of these interventions on cargo carriers trafficking.

RESULTS

We reported that the VPS52 subunit of the Golgi-associated retrograde protein complex is an effector of Rab4b. We found that overexpression of wild type or active Rab4b increased early endosomal to trans-Golgi network retrograde trafficking of the cation-independent mannose-6-phosphate receptor in a Golgi-associated retrograde protein complex-dependent manner. Conversely, overexpression of an inactive Rab4b or Rab4b knockdown attenuated this trafficking. In the absence of Rab4b, the internalized cation-independent mannose 6 phosphate receptor did not have access to VPS52-labeled structures that look like endosomal subdomains and/or endosome-derived carriers, and whose subcellular distribution is Rab4b-independent. Consequently, the cation-independent mannose-6-phosphate receptor was blocked in early endosomes and no longer had access to the trans-Golgi network.

CONCLUSION

Our results support that Rab4b, by controlling the sorting of the cation-independent mannose-6-phosphate receptor towards VPS52 microdomains, confers a directional specificity for cargo carriers en route to the trans-Golgi network. Given the importance of the endocytic recycling in cell homeostasis, disruption of the Rab4b/Golgi-associated retrograde protein complex-dependent step could have serious consequences in pathologies.

Construction of gastric cancer prognostic signature based on the E26 transcription factor and the identification of novel oncogene ELK3

The aim of the present study was to investigate the function of 29 E26 (ETS) transcription factor families in gastric cancer (GC) and determine their association with prognosis. Our analysis of the expression of the ETS family revealed that 28 genes were dysregulated in GC, and that their expression was associated with multiple clinicopathological features (P<0.05). Based on the expression signature of the ETS family, consensus clustering was performed to generate two gastric cancer subtypes. These subtypes exhibited differences in overall survival (OS, P = 0.161), disease-free survival (DFS, P<0.05) and GC grade (P<0.01). Functional enrichment analysis of the target genes associated with the ETS family indicated that these genes primarily contribute to functions that facilitate tumor progression. A systematic statistical analysis was used to construct a prognostic model related to OS and DFS in association with the ETS family. This model demonstrated that the maximum area under the curve (AUC) values for predicting OS and DFS were 0.729 and 0.670, respectively, establishing ETS as an independent prognostic factor for GC Furthermore, a nomogram was created from the prognostic signature, and its predictive accuracy was confirmed by a calibration curve. Finally, the expression and prognostic significance of the six genes comprising the model were also examined. Among these, ELK3 was found to be significantly overexpressed in GC clinical samples. Subsequent in vitro and in vivo studies verified that ELK3 regulates GC proliferation and metastasis, highlighting its potential as a therapeutic target for gastric cancer.

Predictomes: A classifier-curated database of AlphaFold-modeled protein-protein interactions

Protein-protein interactions (PPIs) are ubiquitous in biology, yet a comprehensive structural characterization of the PPIs underlying biochemical processes is lacking. Although AlphaFold-Multimer (AF-M) has the potential to fill this knowledge gap, standard AF-M confidence metrics do not reliably separate relevant PPIs from an abundance of false positive predictions. To address this limitation, we used machine learning on well curated datasets to train a Structure Prediction and Omics informed Classifier called SPOC that shows excellent performance in separating true and false PPIs, including in proteome-wide screens. We applied SPOC to an all-by-all matrix of nearly 300 human genome maintenance proteins, generating ~40,000 predictions that can be viewed at predictomes.org, where users can also score their own predictions with SPOC. High confidence PPIs discovered using our approach suggest novel hypotheses in genome maintenance. Our results provide a framework for interpreting large scale AF-M screens and help lay the foundation for a proteome-wide structural interactome.

Elaboration of the Homer1 Recognition Landscape Reveals Incomplete Divergence of Paralogous EVH1 Domains

Short sequences that mediate interactions with modular binding domains are ubiquitous throughout eukaryotic proteomes. Networks of Short Linear Motifs (SLiMs) and their corresponding binding domains orchestrate many cellular processes, and the low mutational barrier to evolving novel interactions provides a way for biological systems to rapidly sample selectable phenotypes. Mapping SLiM binding specificity and the rules that govern SLiM evolution is fundamental to uncovering the pathways regulated by these networks and developing the tools to manipulate them. We used high-throughput screening of the human proteome to identify sequences that bind to the Enabled/VASP homology 1 (EVH1) domain of the postsynaptic density scaffolding protein Homer1. In doing so, we expanded current understanding of the determinants of Homer EVH1 binding preferences and defined a new motif that can facilitate the discovery of additional Homer-mediated interactions. Interestingly, the Homer1 EVH1 domain preferentially binds to sequences containing an N-terminally overlapping motif that is bound by the paralogous family of Ena/VASP actin polymerases, and many of these sequences can bind to EVH1 domains from both protein families. We provide evidence from orthologous EVH1 domains in pre-metazoan organisms that the overlap in human Ena/VASP and Homer binding preferences corresponds to an incomplete divergence from a common Ena/VASP ancestor. Given this overlap in binding profiles, promiscuous sequences that can be recognized by both families either achieve specificity through extrinsic regulatory strategies or may provide functional benefits via multi-specificity. This may explain why these paralogs incompletely diverged despite the accessibility of further diverged isoforms.

Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

Computational Approaches to Predict Protein-Protein Interactions in Crowded Cellular Environments

Investigating protein-protein interactions is crucial for understanding cellular biological processes because proteins often function within molecular complexes rather than in isolation. While experimental and computational methods have provided valuable insights into these interactions, they often overlook a critical factor: the crowded cellular environment. This environment significantly impacts protein behavior, including structural stability, diffusion, and ultimately the nature of binding. In this review, we discuss theoretical and computational approaches that allow the modeling of biological systems to guide and complement experiments and can thus significantly advance the investigation, and possibly the predictions, of protein-protein interactions in the crowded environment of cell cytoplasm. We explore topics such as statistical mechanics for lattice simulations, hydrodynamic interactions, diffusion processes in high-viscosity environments, and several methods based on molecular dynamics simulations. By synergistically leveraging methods from biophysics and computational biology, we review the state of the art of computational methods to study the impact of molecular crowding on protein-protein interactions and discuss its potential revolutionizing effects on the characterization of the human interactome.

Systematic HOIP interactome profiling reveals critical roles of linear ubiquitination in tissue homeostasis

Linear ubiquitination catalyzed by HOIL-1-interacting protein (HOIP), the key component of the linear ubiquitination assembly complex, plays fundamental roles in tissue homeostasis by executing domain-specific regulatory functions. However, a proteome-wide analysis of the domain-specific interactome of HOIP across tissues is lacking. Here, we present a comprehensive mass spectrometry-based interactome profiling of four HOIP domains in nine mouse tissues. The interaction dataset provides a high-quality HOIP interactome resource with an average of approximately 90 interactors for each bait per tissue. HOIP tissue interactome presents a systematic understanding of linear ubiquitination functions in each tissue and also shows associations of tissue functions to genetic diseases. HOIP domain interactome characterizes a set of previously undefined linear ubiquitinated substrates and elucidates the cross-talk among HOIP domains in physiological and pathological processes. Moreover, we show that linear ubiquitination of Integrin-linked protein kinase (ILK) decreases focal adhesion formation and promotes the detachment of Shigella flexneri-infected cells. Meanwhile, Hoip deficiency decreases the linear ubiquitination of Smad ubiquitination regulatory factor 1 (SMURF1) and enhances its E3 activity, finally causing a reduced bone mass phenotype in mice. Overall, our work expands the knowledge of HOIP-interacting proteins and provides a platform for further discovery of linear ubiquitination functions in tissue homeostasis.

Profiling Protein-Protein Interactions in the Human Brain by Refined Cofractionation Mass Spectrometry

Proteins usually execute their biological functions through interactions with other proteins and by forming macromolecular complexes, but global profiling of protein complexes directly from human tissue samples has been limited. In this study, we utilized cofractionation mass spectrometry (CF-MS) to map protein complexes within the postmortem human brain with experimental replicates. First, we used concatenated anion and cation Ion Exchange Chromatography (IEX) to separate native protein complexes in 192 fractions and then proceeded with Data-Independent Acquisition (DIA) mass spectrometry to analyze the proteins in each fraction, quantifying a total of 4,804 proteins with 3,260 overlapping in both replicates. We improved the DIA's quantitative accuracy by implementing a constant amount of bovine serum albumin (BSA) in each fraction as an internal standard. Next, advanced computational pipelines, which integrate both a database-based complex analysis and an unbiased protein-protein interaction (PPI) search, were applied to identify protein complexes and construct protein-protein interaction networks in the human brain. Our study led to the identification of 486 protein complexes and 10054 binary protein-protein interactions, which represents the first global profiling of human brain PPIs using CF-MS. Overall, this study offers a resource and tool for a wide range of human brain research, including the identification of disease-specific protein complexes in the future.

Key Proteomics Tools for Fundamental and Applied Microalgal Research

Microscopic, photosynthetic prokaryotes and eukaryotes, collectively referred to as microalgae, are widely studied to improve our understanding of key metabolic pathways (e.g., photosynthesis) and for the development of biotechnological applications. Omics technologies, which are now common tools in biological research, have been shown to be critical in microalgal research. In the past decade, significant technological advancements have allowed omics technologies to become more affordable and efficient, with huge datasets being generated. In particular, where studies focused on a single or few proteins decades ago, it is now possible to study the whole proteome of a microalgae. The development of mass spectrometry-based methods has provided this leap forward with the high-throughput identification and quantification of proteins. This review specifically provides an overview of the use of proteomics in fundamental (e.g., photosynthesis) and applied (e.g., lipid production for biofuel) microalgal research, and presents future research directions in this field.

A proteomics study of rheumatoid arthritis patients on etanercept identifies putative biomarkers associated with clinical outcome measures

OBJECTIVES

Biologic DMARDs (bDMARDs) are widely used in patients with RA, but response to bDMARDs is heterogeneous. The objective of this work was to identify pretreatment proteomic biomarkers associated with RA clinical outcome measures in patients starting bDMARDs.

METHODS

Sequential window acquisition of all theoretical fragment ion spectra mass spectrometry (SWATH-MS) was used to generate spectral maps of sera from patients with RA before and after 3 months of treatment with the bDMARD etanercept. Protein levels were regressed against RA clinical outcome measures, i.e. 28-joint DAS (DAS28) and its subcomponents and DAS28 <2.6 (i.e. remission). The proteins with the strongest evidence for association were analysed in an independent, replication dataset. Finally, subnetwork analysis was carried out using the Disease Module Detection algorithm and biological plausibility of identified proteins was assessed by enrichment analysis.

RESULTS

A total of 180 patients with RA were included in the discovery dataset and 58 in the validation dataset from a UK-based prospective multicentre study. Ten individual proteins were found to be significantly associated with RA clinical outcome measures. The association of T-complex protein 1 subunit η with DAS28 remission was replicated in an independent cohort. Subnetwork analysis of the 10 proteins from the regression analysis identified the ontological theme, with the strongest associations being with acute phase and acute inflammatory responses.

CONCLUSION

This longitudinal study of 180 patients with RA commencing etanercept has identified several putative protein biomarkers of treatment response to this drug, one of which was replicated in an independent cohort.

AI-guided pipeline for protein-protein interaction drug discovery identifies a SARS-CoV-2 inhibitor

Protein-protein interactions (PPIs) offer great opportunities to expand the druggable proteome and therapeutically tackle various diseases, but remain challenging targets for drug discovery. Here, we provide a comprehensive pipeline that combines experimental and computational tools to identify and validate PPI targets and perform early-stage drug discovery. We have developed a machine learning approach that prioritizes interactions by analyzing quantitative data from binary PPI assays or AlphaFold-Multimer predictions. Using the quantitative assay LuTHy together with our machine learning algorithm, we identified high-confidence interactions among SARS-CoV-2 proteins for which we predicted three-dimensional structures using AlphaFold-Multimer. We employed VirtualFlow to target the contact interface of the NSP10-NSP16 SARS-CoV-2 methyltransferase complex by ultra-large virtual drug screening. Thereby, we identified a compound that binds to NSP10 and inhibits its interaction with NSP16, while also disrupting the methyltransferase activity of the complex, and SARS-CoV-2 replication. Overall, this pipeline will help to prioritize PPI targets to accelerate the discovery of early-stage drug candidates targeting protein complexes and pathways.

HHIP protein interactions in lung cells provide insight into COPD pathogenesis

Chronic obstructive pulmonary disease (COPD) is the third leading cause of death worldwide. The primary causes of COPD are environmental, including cigarette smoking; however, genetic susceptibility also contributes to COPD risk. Genome-Wide Association Studies (GWASes) have revealed more than 80 genetic loci associated with COPD, leading to the identification of multiple COPD GWAS genes. However, the biological relationships between the identified COPD susceptibility genes are largely unknown. Genes associated with a complex disease are often in close network proximity, i.e. their protein products often interact directly with each other and/or similar proteins. In this study, we use affinity purification mass spectrometry (AP-MS) to identify protein interactions with HHIP , a well-established COPD GWAS gene which is part of the sonic hedgehog pathway, in two disease-relevant lung cell lines (IMR90 and 16HBE). To better understand the network neighborhood of HHIP , its proximity to the protein products of other COPD GWAS genes, and its functional role in COPD pathogenesis, we create HUBRIS, a protein-protein interaction network compiled from 8 publicly available databases. We identified both common and cell type-specific protein-protein interactors of HHIP. We find that our newly identified interactions shorten the network distance between HHIP and the protein products of several COPD GWAS genes, including DSP, MFAP2, TET2 , and FBLN5 . These new shorter paths include proteins that are encoded by genes involved in extracellular matrix and tissue organization. We found and validated interactions to proteins that provide new insights into COPD pathobiology, including CAVIN1 (IMR90) and TP53 (16HBE). The newly discovered HHIP interactions with CAVIN1 and TP53 implicate HHIP in response to oxidative stress.

TM4SF19-mediated control of lysosomal activity in macrophages contributes to obesity-induced inflammation and metabolic dysfunction

Adipose tissue (AT) adapts to overnutrition in a complex process, wherein specialized immune cells remove and replace dysfunctional and stressed adipocytes with new fat cells. Among immune cells recruited to AT, lipid-associated macrophages (LAMs) have emerged as key players in obesity and in diseases involving lipid stress and inflammation. Here, we show that LAMs selectively express transmembrane 4 L six family member 19 (TM4SF19), a lysosomal protein that represses acidification through its interaction with Vacuolar-ATPase. Inactivation of TM4SF19 elevates lysosomal acidification and accelerates the clearance of dying/dead adipocytes in vitro and in vivo. TM4SF19 deletion reduces the LAM accumulation and increases the proportion of restorative macrophages in AT of male mice fed a high-fat diet. Importantly, male mice lacking TM4SF19 adapt to high-fat feeding through adipocyte hyperplasia, rather than hypertrophy. This adaptation significantly improves local and systemic insulin sensitivity, and energy expenditure, offering a potential avenue to combat obesity-related metabolic dysfunction.

3Dmapper: a command line tool for BioBank-scale mapping of variants to protein structures

MOTIVATION

The interpretation of genomic data is crucial to understand the molecular mechanisms of biological processes. Protein structures play a vital role in facilitating this interpretation by providing functional context to genetic coding variants. However, mapping genes to proteins is a tedious and error-prone task due to inconsistencies in data formats. Over the past two decades, numerous tools and databases have been developed to automatically map annotated positions and variants to protein structures. However, most of these tools are web-based and not well-suited for large-scale genomic data analysis.

RESULTS

To address this issue, we introduce 3Dmapper, a stand-alone command-line tool developed in Python and R. It systematically maps annotated protein positions and variants to protein structures, providing a solution that is both efficient and reliable.

AVAILABILITY AND IMPLEMENTATION

https://github.com/vicruiser/3Dmapper.

Population-enriched innate immune variants may identify candidate gene targets at the intersection of cancer and cardio-metabolic disease

Both cancer and cardio-metabolic disease disparities exist among specific populations in the US. For example, African Americans experience the highest rates of breast and prostate cancer mortality and the highest incidence of obesity. Native and Hispanic Americans experience the highest rates of liver cancer mortality. At the same time, Pacific Islanders have the highest death rate attributed to type 2 diabetes (T2D), and Asian Americans experience the highest incidence of non-alcoholic fatty liver disease (NAFLD) and cancers induced by infectious agents. Notably, the pathologic progression of both cancer and cardio-metabolic diseases involves innate immunity and mechanisms of inflammation. Innate immunity in individuals is established through genetic inheritance and external stimuli to respond to environmental threats and stresses such as pathogen exposure. Further, individual genomes contain characteristic genetic markers associated with one or more geographic ancestries (ethnic groups), including protective innate immune genetic programming optimized for survival in their corresponding ancestral environment(s). This perspective explores evidence related to our working hypothesis that genetic variations in innate immune genes, particularly those that are commonly found but unevenly distributed between populations, are associated with disparities between populations in both cancer and cardio-metabolic diseases. Identifying conventional and unconventional innate immune genes that fit this profile may provide critical insights into the underlying mechanisms that connect these two families of complex diseases and offer novel targets for precision-based treatment of cancer and/or cardio-metabolic disease.