Interrogation of RNA-protein interaction dynamics in bacterial growth.

Characterising RNA-protein interaction dynamics is fundamental to understand how bacteria respond to their environment. In this study, we have analysed the dynamics of 91% of the Escherichia coli expressed proteome and the RNA-interaction properties of 271 RNA-binding proteins (RBPs) at different growth phases. We find that 68% of RBPs differentially bind RNA across growth phases and characterise 17 previously unannotated proteins as bacterial RBPs including YfiF, a ncRNA-binding protein. While these new RBPs are mostly present in Proteobacteria, two of them are orthologs of human mitochondrial proteins associated with rare metabolic disorders. Moreover, we reveal novel RBP functions for proteins such as the chaperone HtpG, a new stationary phase tRNA-binding protein. For the first time, the dynamics of the bacterial RBPome have been interrogated, showcasing how this approach can reveal the function of uncharacterised proteins and identify critical RNA-protein interactions for cell growth which could inform new antimicrobial therapies.

N1-methylpseudouridylation of mRNA causes +1 ribosomal frameshifting.

In vitro-transcribed (IVT) mRNAs are modalities that can combat human disease, exemplified by their use as vaccines for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). IVT mRNAs are transfected into target cells, where they are translated into recombinant protein, and the biological activity or immunogenicity of the encoded protein exerts an intended therapeutic effect1,2. Modified ribonucleotides are commonly incorporated into therapeutic IVT mRNAs to decrease their innate immunogenicity3-5, but their effects on mRNA translation fidelity have not been fully explored. Here we demonstrate that incorporation of N1-methylpseudouridine into mRNA results in +1 ribosomal frameshifting in vitro and that cellular immunity in mice and humans to +1 frameshifted products from BNT162b2 vaccine mRNA translation occurs after vaccination. The +1 ribosome frameshifting observed is probably a consequence of N1-methylpseudouridine-induced ribosome stalling during IVT mRNA translation, with frameshifting occurring at ribosome slippery sequences. However, we demonstrate that synonymous targeting of such slippery sequences provides an effective strategy to reduce the production of frameshifted products. Overall, these data increase our understanding of how modified ribonucleotides affect the fidelity of mRNA translation, and although there are no adverse outcomes reported from mistranslation of mRNA-based SARS-CoV-2 vaccines in humans, these data highlight potential off-target effects for future mRNA-based therapeutics and demonstrate the requirement for sequence optimization.

System-wide analysis of RNA and protein subcellular localization dynamics.

Although the subcellular dynamics of RNA and proteins are key determinants of cell homeostasis, their characterization is still challenging. Here we present an integrative framework to simultaneously interrogate the dynamics of the transcriptome and proteome at subcellular resolution by combining two methods: localization of RNA (LoRNA) and a streamlined density-based localization of proteins by isotope tagging (dLOPIT) to map RNA and protein to organelles (nucleus, endoplasmic reticulum and mitochondria) and membraneless compartments (cytosol, nucleolus and cytosolic granules). Interrogating all RNA subcellular locations at once enables system-wide quantification of the proportional distribution of RNA. We obtain a cell-wide overview of localization dynamics for 31,839 transcripts and 5,314 proteins during the unfolded protein response, revealing that endoplasmic reticulum-localized transcripts are more efficiently recruited to cytosolic granules than cytosolic RNAs, and that the translation initiation factor eIF3d is key to sustaining cytoskeletal function. Overall, we provide the most comprehensive overview so far of RNA and protein subcellular localization dynamics.

A Bioconductor workflow for processing, evaluating, and interpreting expression proteomics data.

Background: Expression proteomics involves the global evaluation of protein abundances within a system. In turn, differential expression analysis can be used to investigate changes in protein abundance upon perturbation to such a system. Methods: Here, we provide a workflow for the processing, analysis and interpretation of quantitative mass spectrometry-based expression proteomics data. This workflow utilizes open-source R software packages from the Bioconductor project and guides users end-to-end and step-by-step through every stage of the analyses. As a use-case we generated expression proteomics data from HEK293 cells with and without a treatment. Of note, the experiment included cellular proteins labelled using tandem mass tag (TMT) technology and secreted proteins quantified using label-free quantitation (LFQ). Results: The workflow explains the software infrastructure before focusing on data import, pre-processing and quality control. This is done individually for TMT and LFQ datasets. The application of statistical differential expression analysis is demonstrated, followed by interpretation via gene ontology enrichment analysis. Conclusions: A comprehensive workflow for the processing, analysis and interpretation of expression proteomics is presented. The workflow is a valuable resource for the proteomics community and specifically beginners who are at least familiar with R who wish to understand and make data-driven decisions with regards to their analyses.

PLK1 inhibition dampens NLRP3 inflammasome-elicited response in inflammatory disease models.

Unabated activation of the NLR family pyrin domain-containing 3 (NLRP3) inflammasome is linked with the pathogenesis of various inflammatory disorders. Polo-like kinase 1 (PLK1) has been widely studied for its role in mitosis. Here, using both pharmacological and genetic approaches, we demonstrate that PLK1 promoted NLRP3 inflammasome activation at cell interphase. Using an unbiased proximity-dependent biotin identification (Bio-ID) screen for the PLK1 interactome in macrophages, we show an enhanced proximal association of NLRP3 with PLK1 upon NLRP3 inflammasome activation. We further confirmed the interaction between PLK1 and NLRP3 and identified the interacting domains. Mechanistically, we show that PLK1 orchestrated the microtubule-organizing center (MTOC) structure and NLRP3 subcellular positioning upon inflammasome activation. Treatment with a selective PLK1 kinase inhibitor suppressed IL-1ß production in in vivo inflammatory models, including LPS-induced endotoxemia and monosodium urate-induced peritonitis in mice. Our results uncover a role of PLK1 in regulating NLRP3 inflammasome activation during interphase and identify pharmacological inhibition of PLK1 as a potential therapeutic strategy for inflammatory diseases with excessive NLRP3 inflammasome activation.

Mechanisms and physiological function of daily haemoglobin oxidation rhythms in red blood cells.

Cellular circadian rhythms confer temporal organisation upon physiology that is fundamental to human health. Rhythms are present in red blood cells (RBCs), the most abundant cell type in the body, but their physiological function is poorly understood. Here, we present a novel biochemical assay for haemoglobin (Hb) oxidation status which relies on a redox-sensitive covalent haem-Hb linkage that forms during SDS-mediated cell lysis. Formation of this linkage is lowest when ferrous Hb is oxidised, in the form of ferric metHb. Daily haemoglobin oxidation rhythms are observed in mouse and human RBCs cultured in vitro, or taken from humans in vivo, and are unaffected by mutations that affect circadian rhythms in nucleated cells. These rhythms correlate with daily rhythms in core body temperature, with temperature lowest when metHb levels are highest. Raising metHb levels with dietary sodium nitrite can further decrease daytime core body temperature in mice via nitric oxide (NO) signalling. These results extend our molecular understanding of RBC circadian rhythms and suggest they contribute to the regulation of body temperature.

Enrichment of Membrane Proteins for Downstream Analysis Using Styrene Maleic Acid Lipid Particles (SMALPs) Extraction.

Integral membrane proteins are an important class of cellular proteins. These take part in key cellular processes such as signaling transducing receptors to transporters, many operating within the plasma membrane. More than half of the FDA-approved protein-targeting drugs operate via interaction with proteins that contain at least one membrane-spanning region, yet the characterization and study of their native interactions with therapeutic agents remains a significant challenge. This challenge is due in part to such proteins often being present in small quantities within a cell. Effective solubilization of membrane proteins is also problematic, with the detergents typically employed in solubilizing membranes leading to a loss of functional activity and key interacting partners. In recent years, alternative methods to extract membrane proteins within their native lipid environment have been investigated, with the aim of producing functional nanodiscs, maintaining protein-protein and protein-lipid interactions. A promising approach involves extracting membrane proteins in the form of styrene maleic acid lipid particles (SMALPs) that allow the retention of their native conformation. This extraction method offers many advantages for further protein analysis and allows the study of the protein interactions with other molecules, such as drugs. Here, we describe a protocol for efficient SMALP extraction of functionally active membrane protein complexes within nanodiscs. We showcase the method on the isolation of a low copy number plasma membrane receptor complex, the nicotinic acetylcholine receptor (nAChR), from adult Drosophila melanogaster heads. We demonstrate that these nanodiscs can be used to study native receptor-ligand interactions. This protocol can be applied across many biological scenarios to extract the native conformations of low copy number integral membrane proteins.

Mapping diversity in African trypanosomes using high resolution spatial proteomics.

African trypanosomes are dixenous eukaryotic parasites that impose a significant human and veterinary disease burden on sub-Saharan Africa. Diversity between species and life-cycle stages is concomitant with distinct host and tissue tropisms within this group. Here, the spatial proteomes of two African trypanosome species, Trypanosoma brucei and Trypanosoma congolense, are mapped across two life-stages. The four resulting datasets provide evidence of expression of approximately 5500 proteins per cell-type. Over 2500 proteins per cell-type are classified to specific subcellular compartments, providing four comprehensive spatial proteomes. Comparative analysis reveals key routes of parasitic adaptation to different biological niches and provides insight into the molecular basis for diversity within and between these pathogen species.

Multi-Attribute Monitoring Method for Process Development of Engineered Antibody for Site-Specific Conjugation.

Antibody drug conjugates, a class of biotherapeutic proteins, have been extensively developed in recent years, resulting in new approvals and improved standard of care for cancer patients. Among the numerous strategies of conjugating cytotoxic payloads to monoclonal antibodies, insertion of a cysteine residue achieves a tightly controlled, site-specific drug to antibody ratio. Tailored analytical tools are required to direct the development of processes capable of manufacturing novel antibody scaffolds with the desired product quality. Here, we describe the development of a 12 min, mass-spectrometry-based method capable of monitoring four distinct quality attributes simultaneously: variations in the thiol state of the inserted cysteines, N-linked glycosylation, reduction of interchain disulfide bonds, and polypeptide fragmentation. This method provides new insight into the properties of the antibody intermediate and associated manufacturing processes. Oxidized thiol states are formed within the bioreactor, of which a variant containing an additional disulfide bond was produced and remained relatively constant throughout the fed-batch process; reduced thiol variants were introduced upon harvest. Nearly 20 percent of N-linked glycans contained sialic acid, substantially higher than anticipated for wildtype IgG1. Lastly, previously unreported polypeptide fragmentation sites were identified in the C239i constant domain, and the relationship between fragmentation and glycoform were explored. This work illustrates the utility of applying a high-throughput liquid chromatography-mass spectrometry multi-attribute monitoring method to support the development of engineered antibody scaffolds.

The proteomic landscape of genome-wide genetic perturbations.

Functional genomic strategies have become fundamental for annotating gene function and regulatory networks. Here, we combined functional genomics with proteomics by quantifying protein abundances in a genome-scale knockout library in Saccharomyces cerevisiae, using data-independent acquisition mass spectrometry. We find that global protein expression is driven by a complex interplay of (1) general biological properties, including translation rate, protein turnover, the formation of protein complexes, growth rate, and genome architecture, followed by (2) functional properties, such as the connectivity of a protein in genetic, metabolic, and physical interaction networks. Moreover, we show that functional proteomics complements current gene annotation strategies through the assessment of proteome profile similarity, protein covariation, and reverse proteome profiling. Thus, our study reveals principles that govern protein expression and provides a genome-spanning resource for functional annotation.

Cryptosporidium uses multiple distinct secretory organelles to interact with and modify its host cell.

Cryptosporidium is a leading cause of diarrheal disease in children and an important contributor to early childhood mortality. The parasite invades and extensively remodels intestinal epithelial cells, building an elaborate interface structure. How this occurs at the molecular level and the contributing parasite factors are largely unknown. Here, we generated a whole-cell spatial proteome of the Cryptosporidium sporozoite and used genetic and cell biological experimentation to discover the Cryptosporidium-secreted effector proteome. These findings reveal multiple organelles, including an original secretory organelle, and generate numerous compartment markers by tagging native gene loci. We show that secreted proteins are delivered to the parasite-host interface, where they assemble into different structures including a ring that anchors the parasite into its unique epicellular niche. Cryptosporidium thus uses a complex set of secretion systems during and following invasion that act in concert to subjugate its host cell.

Translation of in vitro-transcribed RNA therapeutics.

In vitro transcribed, modified messenger RNAs (IVTmRNAs) have been used to vaccinate billions of individuals against the SARS-CoV-2 virus, and are currently being developed for many additional therapeutic applications. IVTmRNAs must be translated into proteins with therapeutic activity by the same cellular machinery that also translates native endogenous transcripts. However, different genesis pathways and routes of entry into target cells as well as the presence of modified nucleotides mean that the way in which IVTmRNAs engage with the translational machinery, and the efficiency with which they are being translated, differs from native mRNAs. This review summarises our current knowledge of commonalities and differences in translation between IVTmRNAs and cellular mRNAs, which is key for the development of future design strategies that can generate IVTmRNAs with improved activity in therapeutic applications.

Semi-Supervised Non-Parametric Bayesian Modelling of Spatial Proteomics.

Understanding sub-cellular protein localisation is an essential component in the analysis of context specific protein function. Recent advances in quantitative mass-spectrometry (MS) have led to high resolution mapping of thousands of proteins to sub-cellular locations within the cell. Novel modelling considerations to capture the complex nature of these data are thus necessary. We approach analysis of spatial proteomics data in a non-parametric Bayesian framework, using K-component mixtures of Gaussian process regression models. The Gaussian process regression model accounts for correlation structure within a sub-cellular niche, with each mixture component capturing the distinct correlation structure observed within each niche. The availability of marker proteins (i.e. proteins with a priori known labelled locations) motivates a semi-supervised learning approach to inform the Gaussian process hyperparameters. We moreover provide an efficient Hamiltonian-within-Gibbs sampler for our model. Furthermore, we reduce the computational burden associated with inversion of covariance matrices by exploiting the structure in the covariance matrix. A tensor decomposition of our covariance matrices allows extended Trench and Durbin algorithms to be applied to reduce the computational complexity of inversion and hence accelerate computation. We provide detailed case-studies on Drosophila embryos and mouse pluripotent embryonic stem cells to illustrate the benefit of semi-supervised functional Bayesian modelling of the data.

Inferring differential subcellular localisation in comparative spatial proteomics using BANDLE.

The steady-state localisation of proteins provides vital insight into their function. These localisations are context specific with proteins translocating between different subcellular niches upon perturbation of the subcellular environment. Differential localisation, that is a change in the steady-state subcellular location of a protein, provides a step towards mechanistic insight of subcellular protein dynamics. High-accuracy high-throughput mass spectrometry-based methods now exist to map the steady-state localisation and re-localisation of proteins. Here, we describe a principled Bayesian approach, BANDLE, that uses these data to compute the probability that a protein differentially localises upon cellular perturbation. Extensive simulation studies demonstrate that BANDLE reduces the number of both type I and type II errors compared to existing approaches. Application of BANDLE to several datasets recovers well-studied translocations. In an application to cytomegalovirus infection, we obtain insights into the rewiring of the host proteome. Integration of other high-throughput datasets allows us to provide the functional context of these data.

Operation of a TCA cycle subnetwork in the mammalian nucleus.

Nucleic acid and histone modifications critically depend on the tricarboxylic acid (TCA) cycle for substrates and cofactors. Although a few TCA cycle enzymes have been reported in the nucleus, the corresponding pathways are considered to operate in mitochondria. Here, we show that a part of the TCA cycle is operational also in the nucleus. Using 13C-tracer analysis, we identified activity of glutamine-to-fumarate, citrate-to-succinate, and glutamine-to-aspartate routes in the nuclei of HeLa cells. Proximity labeling mass spectrometry revealed a spatial vicinity of the involved enzymes with core nuclear proteins. We further show nuclear localization of aconitase 2 and 2-oxoglutarate dehydrogenase in mouse embryonic stem cells. Nuclear localization of the latter enzyme, which produces succinyl-CoA, changed from pluripotency to a differentiated state with accompanying changes in the nuclear protein succinylation. Together, our results demonstrate operation of an extended metabolic pathway in the nucleus, warranting a revision of the canonical view on metabolic compartmentalization.

dia-PASEF data analysis using FragPipe and DIA-NN for deep proteomics of low sample amounts.

The dia-PASEF technology uses ion mobility separation to reduce signal interferences and increase sensitivity in proteomic experiments. Here we present a two-dimensional peak-picking algorithm and generation of optimized spectral libraries, as well as take advantage of neural network-based processing of dia-PASEF data. Our computational platform boosts proteomic depth by up to 83% compared to previous work, and is specifically beneficial for fast proteomic experiments and those with low sample amounts. It quantifies over 5300 proteins in single injections recorded at 200 samples per day throughput using Evosep One chromatography system on a timsTOF Pro mass spectrometer and almost 9000 proteins in single injections recorded with a 93-min nanoflow gradient on timsTOF Pro 2, from 200 ng of HeLa peptides. A user-friendly implementation is provided through the incorporation of the algorithms in the DIA-NN software and by the FragPipe workflow for spectral library generation.

CFTR Rescue by Lumacaftor (VX-809) Induces an Extensive Reorganization of Mitochondria in the Cystic Fibrosis Bronchial Epithelium.

BACKGROUND: Cystic Fibrosis (CF) is a genetic disorder affecting around 1 in every 3000 newborns. In the most common mutation, F508del, the defective anion channel, CFTR, is prevented from reaching the plasma membrane (PM) by the quality check control of the cell. Little is known about how CFTR pharmacological rescue impacts the cell proteome. METHODS: We used high-resolution mass spectrometry, differential ultracentrifugation, machine learning and bioinformatics to investigate both changes in the expression and localization of the human bronchial epithelium CF model (F508del-CFTR CFBE41o-) proteome following treatment with VX-809 (Lumacaftor), a drug able to improve the trafficking of CFTR. RESULTS: The data suggested no stark changes in protein expression, yet subtle localization changes of proteins of the mitochondria and peroxisomes were detected. We then used high-content confocal microscopy to further investigate the morphological and compositional changes of peroxisomes and mitochondria under these conditions, as well as in patient-derived primary cells. We profiled several thousand proteins and we determined the subcellular localization data for around 5000 of them using the LOPIT-DC spatial proteomics protocol. CONCLUSIONS: We observed that treatment with VX-809 induces extensive structural and functional remodelling of mitochondria and peroxisomes that resemble the phenotype of healthy cells. Our data suggest additional rescue mechanisms of VX-809 beyond the correction of aberrant folding of F508del-CFTR and subsequent trafficking to the PM.

In Depth Exploration of the Alternative Proteome of Drosophila melanogaster.

Recent studies have shown that hundreds of small proteins were occulted when protein-coding genes were annotated. These proteins, called alternative proteins, have failed to be annotated notably due to the short length of their open reading frame (less than 100 codons) or the enforced rule establishing that messenger RNAs (mRNAs) are monocistronic. Several alternative proteins were shown to be biologically active molecules and seem to be involved in a wide range of biological functions. However, genome-wide exploration of the alternative proteome is still limited to a few species. In the present article, we describe a deep peptidomics workflow which enabled the identification of 401 alternative proteins in Drosophila melanogaster. Subcellular localization, protein domains, and short linear motifs were predicted for 235 of the alternative proteins identified and point toward specific functions of these small proteins. Several alternative proteins had approximated abundances higher than their canonical counterparts, suggesting that these alternative proteins are actually the main products of their corresponding genes. Finally, we observed 14 alternative proteins with developmentally regulated expression patterns and 10 induced upon the heat-shock treatment of embryos, demonstrating stage or stress-specific production of alternative proteins.

Prior Signal Acquisition Software Versions for Orbitrap Underestimate Low Isobaric Mass Tag Intensities, Without Detriment to Differential Abundance Experiments.

Tandem mass tags (TMTs) enable simple and accurate quantitative proteomics for multiplexed samples by relative quantification of tag reporter ions. Orbitrap quantification of reporter ions has been associated with a characteristic notch region in intensity distribution, within which few reporter intensities are recorded. This has been resolved in version 3 of the instrument acquisition software Tune. However, 47% of Orbitrap Fusion, Lumos, or Eclipse submissions to PRIDE were generated using prior software versions. To quantify the impact of the notch on existing quantitative proteomics data, we generated a mixed species benchmark and acquired quantitative data using Tune versions 2 and 3. Intensities below the notch are predominantly underestimated with Tune version 2, leading to overestimation of the true differences in intensities between samples. However, when summarizing reporter ion intensities to higher-level features, such as peptides and proteins, few features are significantly affected. Targeted removal of spectra with reporter ion intensities below the notch is not beneficial for differential peptide or protein testing. Overall, we find that the systematic quantification bias associated with the notch is not detrimental for a typical proteomics experiment.

Drosophila nicotinic acetylcholine receptor subunits and their native interactions with insecticidal peptide toxins.

Drosophila nicotinic acetylcholine receptors (nAChRs) are ligand-gated ion channels that represent a target for insecticides. Peptide neurotoxins are known to block nAChRs by binding to their target subunits, however, a better understanding of this mechanism is needed for effective insecticide design. To facilitate the analysis of nAChRs we used a CRISPR/Cas9 strategy to generate null alleles for all ten nAChR subunit genes in a common genetic background. We studied interactions of nAChR subunits with peptide neurotoxins by larval injections and styrene maleic acid lipid particles (SMALPs) pull-down assays. For the null alleles, we determined the effects of α-Bungarotoxin (α-Btx) and ω-Hexatoxin-Hv1a (Hv1a) administration, identifying potential receptor subunits implicated in the binding of these toxins. We employed pull-down assays to confirm α-Btx interactions with the Drosophila α5 (Dα5), Dα6, Dα7 subunits. Finally, we report the localisation of fluorescent tagged endogenous Dα6 during Drosophila CNS development. Taken together, this study elucidates native Drosophila nAChR subunit interactions with insecticidal peptide toxins and provides a resource for the in vivo analysis of insect nAChRs.