Quantitative proteome dynamics across embryogenesis in a model chordate

The evolution of gene expression programs underlying the development of vertebrates remains poorly characterized. Here, we present a comprehensive proteome atlas of the model chordate Ciona, covering eight developmental stages and ∼7,000 translated genes, accompanied by a multi-omics analysis of co-evolution with the vertebrate Xenopus. Quantitative proteome comparisons argue against the widely held hourglass model, based solely on transcriptomic profiles, whereby peak conservation is observed during mid-developmental stages. Our analysis reveals maximal divergence at these stages, particularly gastrulation and neurulation. Together, our work provides a valuable resource for evaluating conservation and divergence of multi-omics profiles underlying the diversification of vertebrates.

Mitochondrial ATP generation is more proteome efficient than glycolysis

Metabolic efficiency profoundly influences organismal fitness. Nonphotosynthetic organisms, from yeast to mammals, derive usable energy primarily through glycolysis and respiration. Although respiration is more energy efficient, some cells favor glycolysis even when oxygen is available (aerobic glycolysis, Warburg effect). A leading explanation is that glycolysis is more efficient in terms of ATP production per unit mass of protein (that is, faster). Through quantitative flux analysis and proteomics, we find, however, that mitochondrial respiration is actually more proteome efficient than aerobic glycolysis. This is shown across yeast strains, T cells, cancer cells, and tissues and tumors in vivo. Instead of aerobic glycolysis being valuable for fast ATP production, it correlates with high glycolytic protein expression, which promotes hypoxic growth. Aerobic glycolytic yeasts do not excel at aerobic growth but outgrow respiratory cells during oxygen limitation. We accordingly propose that aerobic glycolysis emerges from cells maintaining a proteome conducive to both aerobic and hypoxic growth.

Evidence for widespread cytoplasmic structuring into mesoscale condensates

Compartmentalization is an essential feature of eukaryotic life and is achieved both via membrane-bound organelles, such as mitochondria, and membrane-less biomolecular condensates, such as the nucleolus. Known biomolecular condensates typically exhibit liquid-like properties and are visualized by microscopy on the scale of ~1 µm (refs. 1,2). They have been studied mostly by microscopy, examining select individual proteins. So far, several dozen biomolecular condensates are known, serving a multitude of functions, for example, in the regulation of transcription3, RNA processing4 or signalling5,6, and their malfunction can cause diseases7,8. However, it remains unclear to what extent biomolecular condensates are utilized in cellular organization and at what length scale they typically form. Here we examine native cytoplasm from Xenopus egg extract on a global scale with quantitative proteomics, filtration, size exclusion and dilution experiments. These assays reveal that at least 18% of the proteome is organized into mesoscale biomolecular condensates at the scale of ~100 nm and appear to be stabilized by RNA or gelation. We confirmed mesoscale sizes via imaging below the diffraction limit by investigating protein permeation into porous substrates with defined pore sizes. Our results show that eukaryotic cytoplasm organizes extensively via biomolecular condensates, but at surprisingly short length scales.

Systematic identification and characterization of genes in the regulation and biogenesis of photosynthetic machinery

Photosynthesis is central to food production and the Earth's biogeochemistry, yet the molecular basis for its regulation remains poorly understood. Here, using high-throughput genetics in the model eukaryotic alga Chlamydomonas reinhardtii, we identify with high confidence (false discovery rate [FDR] < 0.11) 70 poorly characterized genes required for photosynthesis. We then enable the functional characterization of these genes by providing a resource of proteomes of mutant strains, each lacking one of these genes. The data allow assignment of 34 genes to the biogenesis or regulation of one or more specific photosynthetic complexes. Further analysis uncovers biogenesis/regulatory roles for at least seven proteins, including five photosystem I mRNA maturation factors, the chloroplast translation factor MTF1, and the master regulator PMR1, which regulates chloroplast genes via nuclear-expressed factors. Our work provides a rich resource identifying regulatory and functional genes and placing them into pathways, thereby opening the door to a system-level understanding of photosynthesis.

Phase plane dynamics of ERK phosphorylation

The extracellular signal-regulated kinase (ERK) controls multiple critical processes in the cell and is deregulated in human cancers, congenital abnormalities, immune diseases, and neurodevelopmental syndromes. Catalytic activity of ERK requires dual phosphorylation by an upstream kinase, in a mechanism that can be described by two sequential Michaelis-Menten steps. The estimation of individual reaction rate constants from kinetic data in the full mechanism has proved challenging. Here, we present an analytically tractable approach to parameter estimation that is based on the phase plane representation of ERK activation and yields two combinations of six reaction rate constants in the detailed mechanism. These combinations correspond to the ratio of the specificities of two consecutive phosphorylations and the probability that monophosphorylated substrate does not dissociate from the enzyme before the second phosphorylation. The presented approach offers a language for comparing the effects of mutations that disrupt ERK activation and function in vivo. As an illustration, we use phase plane representation to analyze dual phosphorylation under heterozygous conditions, when two enzyme variants compete for the same substrate.

Super-Resolution Mass Spectrometry Enables Rapid, Accurate, and Highly Multiplexed Proteomics at the MS2 Level

In tandem mass spectrometry (MS2)-based multiplexed quantitative proteomics, the complement reporter ion approaches (TMTc and TMTproC) were developed to eliminate the ratio-compression problem of conventional MS2-level approaches. Resolving all high m/z complement reporter ions (∼6.32 mDa-spaced) requires mass resolution and scan speeds above the performance levels of OrbitrapTM instruments. Therefore, complement reporter ion quantification with TMT/TMTpro reagents is currently limited to 5 out of 11 (TMT) or 9 out of 18 (TMTpro) channels (∼1 Da spaced). We first demonstrate that a FusionTM LumosTM Orbitrap can resolve 6.32 mDa-spaced complement reporter ions with standard acquisition modes extended with 3 s transients. We then implemented a super-resolution mass spectrometry approach using the least-squares fitting (LSF) method for processing Orbitrap transients to achieve shotgun proteomics-compatible scan rates. The LSF performance resolves the 6.32 mDa doublets for all TMTproC channels in the standard mass range with transients as short as ∼108 ms (Orbitrap resolution setting of 50,000 at m/z 200). However, we observe a slight decrease in measurement precision compared to 1 Da spacing with the 108 ms transients. With 256 ms transients (resolution of 120,000 at m/z 200), coefficients of variation are essentially indistinguishable from 1 Da samples. We thus demonstrate the feasibility of highly multiplexed, accurate, and precise shotgun proteomics at the MS2 level.

Slow TCA flux and ATP production in primary solid tumours but not metastases

Tissues derive ATP from two pathways-glycolysis and the tricarboxylic acid (TCA) cycle coupled to the electron transport chain. Most energy in mammals is produced via TCA metabolism1. In tumours, however, the absolute rates of these pathways remain unclear. Here we optimize tracer infusion approaches to measure the rates of glycolysis and the TCA cycle in healthy mouse tissues, Kras-mutant solid tumours, metastases and leukaemia. Then, given the rates of these two pathways, we calculate total ATP synthesis rates. We find that TCA cycle flux is suppressed in all five primary solid tumour models examined and is increased in lung metastases of breast cancer relative to primary orthotopic tumours. As expected, glycolysis flux is increased in tumours compared with healthy tissues (the Warburg effect2,3), but this increase is insufficient to compensate for low TCA flux in terms of ATP production. Thus, instead of being hypermetabolic, as commonly assumed, solid tumours generally produce ATP at a slower than normal rate. In mouse pancreatic cancer, this is accommodated by the downregulation of protein synthesis, one of this tissue's major energy costs. We propose that, as solid tumours develop, cancer cells shed energetically expensive tissue-specific functions, enabling uncontrolled growth despite a limited ability to produce ATP.

Gut bacterial nutrient preferences quantified in vivo

Great progress has been made in understanding gut microbiomes' products and their effects on health and disease. Less attention, however, has been given to the inputs that gut bacteria consume. Here, we quantitatively examine inputs and outputs of the mouse gut microbiome, using isotope tracing. The main input to microbial carbohydrate fermentation is dietary fiber and to branched-chain fatty acids and aromatic metabolites is dietary protein. In addition, circulating host lactate, 3-hydroxybutyrate, and urea (but not glucose or amino acids) feed the gut microbiome. To determine the nutrient preferences across bacteria, we traced into genus-specific bacterial protein sequences. We found systematic differences in nutrient use: most genera in the phylum Firmicutes prefer dietary protein, Bacteroides dietary fiber, and Akkermansia circulating host lactate. Such preferences correlate with microbiome composition changes in response to dietary modifications. Thus, diet shapes the microbiome by promoting the growth of bacteria that preferentially use the ingested nutrients.

The Shuttling Cascade in Lasso Peptide Benenodin-1 is Controlled by Non-Covalent Interactions

The lasso peptide benenodin-1, a naturally occurring and bacterially produced [1]rotaxane, undergoes a reversible zip tie-like motion under heat activation, in which a peptidic wheel stepwise translates along a molecular thread in a cascade of "tail/loop pulling" equilibria. Conformational and structural analyses of four translational isomers, in solution and in the gas phase, reveal that the equilibrium distribution is controlled by mechanical and non-covalent forces within the lasso peptide. Furthermore, each dynamic pulling step is accompanied by a major restructuring of the intramolecular hydrogen bonding network between wheel and thread, which affects the peptide's physico-chemical properties.

Activity-based RNA-modifying enzyme probing reveals DUS3L-mediated dihydrouridylation

Epitranscriptomic RNA modifications can regulate RNA activity; however, there remains a major gap in our understanding of the RNA chemistry present in biological systems. Here we develop RNA-mediated activity-based protein profiling (RNABPP), a chemoproteomic strategy that relies on metabolic RNA labeling, mRNA interactome capture and quantitative proteomics, to investigate RNA-modifying enzymes in human cells. RNABPP with 5-fluoropyrimidines allowed us to profile 5-methylcytidine (m⁵C) and 5-methyluridine (m⁵U) methyltransferases. Further, we uncover a new mechanism-based crosslink between 5-fluorouridine (5-FUrd)-modified RNA and the dihydrouridine synthase (DUS) homolog DUS3L. We investigate the mechanism of crosslinking and use quantitative nucleoside liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis and 5-FUrd-based crosslinking and immunoprecipitation (CLIP) sequencing to map DUS3L-dependent dihydrouridine (DHU) modifications across the transcriptome. Finally, we show that DUS3L-knockout (KO) cells have compromised protein translation rates and impaired cellular proliferation. Taken together, our work provides a general approach for profiling RNA-modifying enzyme activity in living cells and reveals new pathways for epitranscriptomic RNA regulation.

GCN2 adapts protein synthesis to scavenging-dependent growth

Pancreatic cancer cells with limited access to free amino acids can grow by scavenging extracellular protein. In a murine model of pancreatic cancer, we performed a genome-wide CRISPR screen for genes required for scavenging-dependent growth. The screen identified key mediators of macropinocytosis, peripheral lysosome positioning, endosome-lysosome fusion, lysosomal protein catabolism, and translational control. The top hit was GCN2, a kinase that suppresses translation initiation upon amino acid depletion. Using isotope tracers, we show that GCN2 is not required for protein scavenging. Instead, GCN2 prevents ribosome stalling but without slowing protein synthesis; cells still use all of the limiting amino acids as they emerge from lysosomes. GCN2 also adapts gene expression to the nutrient-poor environment, reorienting protein synthesis away from ribosomes and toward lysosomal hydrolases, such as cathepsin L. GCN2, cathepsin L, and the other genes identified in the screen are potential therapeutic targets in pancreatic cancer.

Evaluating the Arrhenius equation for developmental processes

The famous Arrhenius equation is well suited to describing the temperature dependence of chemical reactions but has also been used for complicated biological processes. Here, we evaluate how well the simple Arrhenius equation predicts complex multi-step biological processes, using frog and fruit fly embryogenesis as two canonical models. We find that the Arrhenius equation provides a good approximation for the temperature dependence of embryogenesis, even though individual developmental intervals scale differently with temperature. At low and high temperatures, however, we observed significant departures from idealized Arrhenius Law behavior. When we model multi-step reactions of idealized chemical networks, we are unable to generate comparable deviations from linearity. In contrast, we find the two enzymes GAPDH and β-galactosidase show non-linearity in the Arrhenius plot similar to our observations of embryonic development. Thus, we find that complex embryonic development can be well approximated by the simple Arrhenius equation regardless of non-uniform developmental scaling and propose that the observed departure from this law likely results more from non-idealized individual steps rather than from the complexity of the system.

TMTpro Complementary Ion Quantification Increases Plexing and Sensitivity for Accurate Multiplexed Proteomics at the MS2 Level

Multiplexed proteomics is a powerful tool to assay cell states in health and disease, but accurate quantification of relative protein changes is impaired by interference from co-isolated peptides. Interference can be reduced by using MS3-based quantification, but this reduces sensitivity and requires specialized instrumentation. An alternative approach is quantification by complementary ions, the balancer group-peptide conjugates, which allows accurate and precise multiplexed quantification at the MS2 level and is compatible with most proteomics instruments. However, complementary ions of the popular TMT-tag form inefficiently and multiplexing is limited to five channels. Here, we evaluate and optimize complementary ion quantification for the recently released TMTpro-tag, which increases complementary ion plexing capacity to eight channels (TMTproC). Furthermore, the beneficial fragmentation properties of TMTpro increase sensitivity for TMTproC, resulting in ∼65% more proteins quantified compared to TMTpro-MS3 and ∼18% more when compared to real-time-search TMTpro-MS3 (RTS-SPS-MS3). TMTproC quantification is more accurate than TMTpro-MS2 and even superior to RTS-SPS-MS3. We provide the software for quantifying TMTproC data as an executable that is compatible with the MaxQuant analysis pipeline. Thus, TMTproC advances multiplexed proteomics data quality and widens access to accurate multiplexed proteomics beyond laboratories with MS3-capable instrumentation.

Precise Temporal Regulation of Post-transcriptional Repressors Is Required for an Orderly Drosophila Maternal-to-Zygotic Transition

In animal embryos, the maternal-to-zygotic transition (MZT) hands developmental control from maternal to zygotic gene products. We show that the maternal proteome represents more than half of the protein-coding capacity of Drosophila melanogaster’s genome, and that 2% of this proteome is rapidly degraded during the MZT. Cleared proteins include the post-transcriptional repressors Cup, Trailer hitch (TRAL), Maternal expression at 31B (ME31B), and Smaug (SMG). Although the ubiquitin-proteasome system is necessary for clearance of these repressors, distinct E3 ligase complexes target them: the C-terminal to Lis1 Homology (CTLH) complex targets Cup, TRAL, and ME31B for degradation early in the MZT and the Skp/Cullin/F-box-containing (SCF) complex targets SMG at the end of the MZT. Deleting the C-terminal 233 amino acids of SMG abrogates F-box protein interaction and confers immunity to degradation. Persistent SMG downregulates zygotic re-expression of mRNAs whose maternal contribution is degraded by SMG. Thus, clearance of SMG permits an orderly MZT.

Cao et al. show that 2% of the proteome is degraded in early Drosophila embryos, including a repressive ribonucleoprotein complex. Two E3 ubiquitin ligases separately act on distinct components of this complex to phase their clearance. Failure to degrade a key component, the Smaug RNA-binding protein, disrupts an orderly maternal-to-zygotic transition.

Inference of Multisite Phosphorylation Rate Constants and Their Modulation by Pathogenic Mutations

Multisite protein phosphorylation plays a critical role in cell regulation [1-3]. It is widely appreciated that the functional capabilities of multisite phosphorylation depend on the order and kinetics of phosphorylation steps, but kinetic aspects of multisite phosphorylation remain poorly understood [4-6]. Here, we focus on what appears to be the simplest scenario, when a protein is phosphorylated on only two sites in a strict, well-defined order. This scenario describes the activation of ERK, a highly conserved cell-signaling enzyme. We use Bayesian parameter inference in a structurally identifiable kinetic model to dissect dual phosphorylation of ERK by MEK, a kinase that is mutated in a large number of human diseases [7-12]. Our results reveal how enzyme processivity and efficiencies of individual phosphorylation steps are altered by pathogenic mutations. The presented approach, which connects specific mutations to kinetic parameters of multisite phosphorylation mechanisms, provides a systematic framework for closing the gap between studies with purified enzymes and their effects in the living organism.

A Click-Chemistry-Based Enrichable Crosslinker for Structural and Protein Interaction Analysis by Mass Spectrometry

Mass spectrometry is the method of choice for the characterisation of proteomes. Most proteins operate in protein complexes, in which their close association modulates their function. However, with standard MS analysis, information on protein–protein interactions is lost and no structural information is retained. To gain structural and interactome data, new crosslinking reagents are needed that freeze inter‐ and intramolecular interactions. Herein, the development of a new reagent, which has several features that enable highly sensitive crosslinking MS, is reported. The reagent enables enrichment of crosslinked peptides from the majority of background peptides to facilitate efficient detection of low‐abundant crosslinked peptides. Due to the special cleavable properties, the reagent can be used for MS² and potentially for MS³ experiments. Thus, the new crosslinking reagent, in combination with high‐end MS, should enable sensitive analysis of interactomes, which will help researchers to obtain important insights into cellular states in health and diseases.

Bayesian Confidence Intervals for Multiplexed Proteomics Integrate Ion-statistics with Peptide Quantification Concordance

Multiplexed proteomics has emerged as a powerful tool to measure relative protein expression levels across multiple conditions. The relative protein abundances are inferred by comparing the signals generated by isobaric tags, which encode the samples' origins. Intuitively, the trust associated with a protein measurement depends on the similarity of ratios from the protein's peptides and the signal-strength of these measurements. However, typically the average peptide ratio is reported as the estimate of relative protein abundance, which is only the most likely ratio with a very naive model. Moreover, there is no sense on the confidence in these measurements. Here, we present a mathematically rigorous approach that integrates peptide signal strengths and peptide-measurement agreement into an estimation of the true protein ratio and the associated confidence (BACIQ). The main advantages of BACIQ are: (1) It removes the need to threshold reported peptide signal based on an arbitrary cut-off, thereby reporting more measurements from a given experiment; (2) Confidence can be assigned without replicates; (3) For repeated experiments BACIQ provides confidence intervals for the union, not the intersection, of quantified proteins; (4) For repeated experiments, BACIQ confidence intervals are more predictive than confidence intervals based on protein measurement agreement. To demonstrate the power of BACIQ we reanalyzed previously published data on subcellular protein movement on treatment with an Exportin-1 inhibiting drug. We detect ∼2× more highly significant movers, down to subcellular localization changes of ∼1%. Thus, our method drastically increases the value obtainable from quantitative proteomics experiments, helping researchers to interpret their data and prioritize resources. To make our approach easily accessible we distribute it via a Python/Stan package.

The Synthetic Phenotype of ΔbamB ΔbamE Double Mutants Results from a Lethal Jamming of the Bam Complex by the Lipoprotein RcsF

The selective permeability of the Gram-negative outer membrane (OM) is maintained by integral β-barrel outer membrane proteins (OMPs). The heteropentomeric β-barrel assembly machine (Bam) folds and inserts OMPs into the OM. Coordination of the essential proteins BamA and BamD is critical for OMP assembly and therefore the viability of the cell. The role of the nonessential lipoproteins BamBCE has yet to be characterized; however, genetic evidence suggests that they have nonoverlapping roles in OMP assembly. In this work, we quantify changes of the proteome in the conditional lethal ΔbamB ΔbamE double mutant. We show that cells lacking BamB and BamE have a global OMP defect that is a result of a lethal obstruction of an assembly-competent Bam complex by the lipoprotein RcsF. RcsF is a stress-sensing lipoprotein that is threaded through the lumen of abundant β-barrel OMPs by the Bam complex to expose the amino terminus on the cell surface. We demonstrate that simply removing this lipoprotein corrects the severe OMP assembly defect of the double mutant nearly as efficiently as a previously isolated suppressor mutation in bamA We propose that BamB and BamE play crucial, nonoverlapping roles to coordinate the activities of BamA and BamD during OMP biogenesis.IMPORTANCE Protein assembly into lipid bilayers is an essential process that ensures the viability of diverse organisms. In Gram-negative bacteria, the heteropentomeric β-barrel assembly machine (Bam) folds and inserts proteins into the outer membrane. Due to its essentiality, outer membrane protein (OMP) assembly by the Bam complex is an attractive target for antibiotic development. Here, we show that the conditional lethal phenotype of a mutant lacking two of the three nonessential lipoproteins, BamB and BamE, is caused by lethal jamming of the stripped-down Bam complex by a normally surface-exposed lipoprotein, RcsF. The heterotrimeric Bam complex (BamA, BamD, BamC) is nearly as efficient as the wild-type complex in OMP assembly if RcsF is removed. Our study highlights the importance of BamB and BamE in regulating the interaction between BamA and BamD and expands our understanding of the role of the Bam complex in outer membrane biogenesis.

A Review on Quantitative Multiplexed Proteomics

Over the last few decades, mass spectrometry-based proteomics has become an increasingly powerful tool that is now able to routinely detect and quantify thousands of proteins. A major advance for global protein quantification was the introduction of isobaric tags, which, in a single experiment, enabled the global quantification of proteins across multiple samples. Herein, these methods are referred to as multiplexed proteomics. The principles, advantages, and drawbacks of various multiplexed proteomics techniques are discussed and compared with alternative approaches. We also discuss how the emerging combination of multiplexing with targeted proteomics might enable the reliable and high-quality quantification of very low abundance proteins across multiple conditions. Lastly, we suggest that fusing multiplexed proteomics with data-independent acquisition approaches might enable the comparison of hundreds of different samples without missing values, while maintaining the superb measurement precision and accuracy obtainable with isobaric tag quantification.

Proteotoxicity from aberrant ribosome biogenesis compromises cell fitness

To achieve maximal growth, cells must manage a massive economy of ribosomal proteins (r-proteins) and RNAs (rRNAs) to produce thousands of ribosomes every minute. Although ribosomes are essential in all cells, natural disruptions to ribosome biogenesis lead to heterogeneous phenotypes. Here, we model these perturbations in Saccharomyces cerevisiae and show that challenges to ribosome biogenesis result in acute loss of proteostasis. Imbalances in the synthesis of r-proteins and rRNAs lead to the rapid aggregation of newly synthesized orphan r-proteins and compromise essential cellular processes, which cells alleviate by activating proteostasis genes. Exogenously bolstering the proteostasis network increases cellular fitness in the face of challenges to ribosome assembly, demonstrating the direct contribution of orphan r-proteins to cellular phenotypes. We propose that ribosome assembly is a key vulnerability of proteostasis maintenance in proliferating cells that may be compromised by diverse genetic, environmental, and xenobiotic perturbations that generate orphan r-proteins.

Cells are made up of thousands of different proteins that perform unique roles required for life. To create all of these proteins, cells use machines called ribosomes that are partly formed of elements known as r-proteins. When cells grow and divide, the ribosomes have to make copies of themselves through a process called ribosome biogenesis.

Although all cells need ribosomes, certain types of cells are especially sensitive to events that interfere with ribosome biogenesis. For example, patients that have mutations in genes needed for ribosome biogenesis produce fewer red blood cells, but their other cells and tissues are mostly healthy. It is not clear why some cells are more sensitive than others.

Ribosome biogenesis is very similar between different organisms, so researchers often use budding yeast as a model to study the process. Here, Tye et al. used genetic and chemical tools to interfere with ribosome biogenesis on short time scales, which made it possible to detect early on what was going wrong in the cells.

The experiments found that when ribosome biogenesis was disrupted, r-proteins that were waiting to be assembled into ribosomes quickly stuck to one another and formed clumps that reduced the ability of the yeast cells to grow. The cells responded by switching on a protein called Hsf1, which restored their ability to grow. Yeast cells that were growing quickly, and therefore making more ribosomes, were more sensitive to abnormal ribosome biogenesis than slow-growing cells.

These results indicate that how actively a cell is growing, and its ability to cope with r-proteins sticking together, may in part explain why certain cells are more vulnerable to events that interfere with ribosome biogenesis. Since human cells also have Hsf1, future experiments could investigate whether turning it on might also protect fast-growing human cells from such events.

Proteomics of nucleocytoplasmic partitioning

The partitioning of the proteome between nucleus and cytoplasm affects nearly every aspect of eukaryotic biology. Despite this central role, we still have a poor understanding of which proteins localize in the nucleus and how this varies in different cell types and conditions. Recent advances in quantitative proteomics and high-throughput imaging are starting to close this knowledge gap. Studies on protein interaction are beginning to reveal the spectrum of cargos of nuclear import and export receptors. We anticipate that it will soon be possible to predict each protein's nucleocytoplasmic localization based on its importin/exportin interactions and its estimated diffusion rate through the nuclear pore. This insight is likely to provide us with a fundamental understanding of how cells use nucleocytoplasmic partitioning to encode and relay information.

Accurate, Sensitive, and Precise Multiplexed Proteomics Using the Complement Reporter Ion Cluster

Quantitative analysis of proteomes across multiple time points, organelles, and perturbations is essential for understanding both fundamental biology and disease states. The development of isobaric tags (e.g., TMT) has enabled the simultaneous measurement of peptide abundances across several different conditions. These multiplexed approaches are promising in principle because of advantages in throughput and measurement quality. However, in practice, existing multiplexing approaches suffer from key limitations. In its simple implementation (TMT-MS2), measurements are distorted by chemical noise leading to poor measurement accuracy. The current state-of-the-art (TMT-MS3) addresses this but requires specialized quadrupole-iontrap-Orbitrap instrumentation. The complement reporter ion approach (TMTc) produces high accuracy measurements and is compatible with many more instruments, like quadrupole-Orbitraps. However, the required deconvolution of the TMTc cluster leads to poor measurement precision. Here, we introduce TMTc+, which adds the modeling of the MS2-isolation step into the deconvolution algorithm. The resulting measurements are comparable in precision to TMT-MS3/MS2. The improved duty cycle and lower filtering requirements make TMTc+ more sensitive than TMT-MS3 and comparable with TMT-MS2. At the same time, unlike TMT-MS2, TMTc+ is exquisitely able to distinguish signal from chemical noise even outperforming TMT-MS3. Lastly, we compare TMTc+ to quantitative label-free proteomics of total HeLa lysate and find that TMTc+ quantifies 7.8k versus 3.9k proteins in a 5-plex sample. At the same time, the median coefficient of variation improves from 13% to 4%. Thus, TMTc+ advances quantitative proteomics by enabling accurate, sensitive, and precise multiplexed experiments on more commonly used instruments.

A Sulfoxide-Based Isobaric Labelling Reagent for Accurate Quantitative Mass Spectrometry

Modern proteomics requires reagents for exact quantification of peptides in complex mixtures. Peptide labelling is most typically achieved with isobaric tags that consist of a balancer and a reporter part that separate in the gas phase. An ingenious distribution of stable isotopes provides multiple reagents with identical molecular weight but a different mass of the reporter groups, allowing relative quantification of multiple samples in one measurement. Here we report a new isobaric labelling reagent, where the balancer and the reporter are linked by a sulfoxide group, which, based on the sulfoxide pyrolysis, leads to easy and asymmetric cleavage at low fragmentation energy. The fragmentation of our new design is significantly improved, yielding more intense complementary ion signals, allowing complementary ion cluster analysis as well.

Quantitative Proteomics for Xenopus Embryos II, Data Analysis

The oocytes, embryos, and cell-free lysates of the frog Xenopus laevis have emerged as powerful models for quantitative proteomic experiments. In the accompanying paper (Chapter 13) we describe how to prepare samples and acquire multiplexed proteomics spectra from those. As an illustrative example we use a 10-stage developmental time series from the egg to stage 35 (just before hatching). Here, we outline how to convert the ~700,000 acquired mass spectra from this time series into protein expression dynamics for ~9000 proteins. We first outline a preliminary quality-control analysis to discover any errors that occurred during sample preparation. We discuss how peptide and protein identification error rates are controlled, and how peptide and protein species are quantified. Our analysis relies on the freely available MaxQuant proteomics pipeline. Finally, we demonstrate how to start interpreting this large dataset by clustering and gene-set enrichment analysis.

Quantitative Proteomics of Xenopus Embryos I, Sample Preparation

Xenopus oocytes and embryos are model systems optimally suited for quantitative proteomics. This is due to the availability of large amount of protein material and the ease of physical manipulation. Furthermore, facile in vitro fertilization provides superbly synchronized embryos for cell cycle and developmental stages. Here, we detail protocols developed over the last few years for sample preparation of multiplexed proteomics with TMT-tags followed by quantitative mass spectrometry analysis using the MultiNotch MS3 approach. In this approach, each condition is barcoded with an isobaric tag at the peptide level. After barcoding, samples are combined and the relative abundance of ~100,000 peptides is quantified on a mass spectrometer. High reproducibility of the sample preparation process prior to peptides being tagged and combined is of upmost importance for obtaining unbiased data. Otherwise, differences in sample handling can inadvertently appear as biological changes. We detail and exemplify the application of our sample workflow on an embryonic time-series of ten developmental stages of Xenopus laevis embryos ranging from the egg to stage 35 (just before hatching). Our accompanying paper (Chapter 14 ) details a bioinformatics pipeline to analyze the quality of the given sample preparation and strategies to convert spectra of X. laevis peptides into biologically interpretable data.

Degradation of the BAF Complex Factor BRD9 by Heterobifunctional Ligands

The bromodomain-containing protein BRD9, a subunit of the human BAF (SWI/SNF) nucleosome remodeling complex, has emerged as an attractive therapeutic target in cancer. Despite the development of chemical probes targeting the BRD9 bromodomain, there is a limited understanding of BRD9 function beyond acetyl-lysine recognition. We have therefore created the first BRD9-directed chemical degraders, through iterative design and testing of heterobifunctional ligands that bridge the BRD9 bromodomain and the cereblon E3 ubiquitin ligase complex. Degraders of BRD9 exhibit markedly enhanced potency compared to parental ligands (10- to 100-fold). Parallel study of degraders with divergent BRD9-binding chemotypes in models of acute myeloid leukemia resolves bromodomain polypharmacology in this emerging drug class. Together, these findings reveal the tractability of non-BET bromodomain containing proteins to chemical degradation, and highlight lead compound dBRD9 as a tool for the study of BRD9.

A Strategy to Combine Sample Multiplexing with Targeted Proteomics Assays for High-Throughput Protein Signature Characterization

Targeted mass spectrometry assays for protein quantitation monitor peptide surrogates, which are easily multiplexed to target many peptides in a single assay. However, these assays have generally not taken advantage of sample multiplexing, which allows up to ten analyses to occur in parallel. We present a two-dimensional multiplexing workflow that utilizes synthetic peptides for each protein to prompt the simultaneous quantification of >100 peptides from up to ten mixed sample conditions. We demonstrate that targeted analysis of unfractionated lysates (2 hr) accurately reproduces the quantification of fractionated lysates (72 hr analysis) while obviating the need for peptide detection prior to quantification. We targeted 131 peptides corresponding to 69 proteins across all 60 National Cancer Institute cell lines in biological triplicate, analyzing 180 samples in only 48 hr (the equivalent of 16 min/sample). These data further elucidated a correlation between the expression of key proteins and their cellular response to drug treatment.

Vertebrate Embryonic Cleavage Pattern Determination

The pattern of the earliest cell divisions in a vertebrate embryo lays the groundwork for later developmental events such as gastrulation, organogenesis, and overall body plan establishment. Understanding these early cleavage patterns and the mechanisms that create them is thus crucial for the study of vertebrate development. This chapter describes the early cleavage stages for species representing ray-finned fish, amphibians, birds, reptiles, mammals, and proto-vertebrate ascidians and summarizes current understanding of the mechanisms that govern these patterns. The nearly universal influence of cell shape on orientation and positioning of spindles and cleavage furrows and the mechanisms that mediate this influence are discussed. We discuss in particular models of aster and spindle centering and orientation in large embryonic blastomeres that rely on asymmetric internal pulling forces generated by the cleavage furrow for the previous cell cycle. Also explored are mechanisms that integrate cell division given the limited supply of cellular building blocks in the egg and several-fold changes of cell size during early development, as well as cytoskeletal specializations specific to early blastomeres including processes leading to blastomere cohesion. Finally, we discuss evolutionary conclusions beginning to emerge from the contemporary analysis of the phylogenetic distributions of cleavage patterns. In sum, this chapter seeks to summarize our current understanding of vertebrate early embryonic cleavage patterns and their control and evolution.

Asymmetries in Cell Division, Cell Size, and Furrowing in the Xenopus laevis Embryo

Asymmetric cell divisions produce two daughter cells with distinct fate. During embryogenesis, this mechanism is fundamental to build tissues and organs because it generates cell diversity. In adults, it remains crucial to maintain stem cells. The enthusiasm for asymmetric cell division is not only motivated by the beauty of the mechanism and the fundamental questions it raises, but has also very pragmatic reasons. Indeed, misregulation of asymmetric cell divisions is believed to have dramatic consequences potentially leading to pathogenesis such as cancers. In diverse model organisms, asymmetric cell divisions result in two daughter cells, which differ not only by their fate but also in size. This is the case for the early Xenopus laevis embryo, in which the two first embryonic divisions are perpendicular to each other and generate two pairs of blastomeres, which usually differ in size: one pair of blastomeres is smaller than the other. Small blastomeres will produce embryonic dorsal structures, whereas the larger pair will evolve into ventral structures. Here, we present a speculative model on the origin of the asymmetry of this cell division in the Xenopus embryo. We also discuss the apparently coincident asymmetric distribution of cell fate determinants and cell-size asymmetry of the 4-cell stage embryo. Finally, we discuss the asymmetric furrowing during epithelial cell cytokinesis occurring later during Xenopus laevis embryo development.

Generic Theoretical Models to Predict Division Patterns of Cleaving Embryos

Life for all animals starts with a precise 3D choreography of reductive divisions of the fertilized egg, known as cleavage patterns. These patterns exhibit conserved geometrical features and striking interspecies invariance within certain animal classes. To identify the generic rules that may govern these morphogenetic events, we developed a 3D-modeling framework that iteratively infers blastomere division positions and orientations, and consequent multicellular arrangements. From a minimal set of parameters, our model predicts detailed features of cleavage patterns in the embryos of fishes, amphibians, echinoderms, and ascidians, as well as the genetic and physical perturbations that alter these patterns. This framework demonstrates that a geometrical system based on length-dependent microtubule forces that probe blastomere shape and yolk gradients, biased by cortical polarity domains, may dictate division patterns and overall embryo morphogenesis. These studies thus unravel the default self-organization rules governing early embryogenesis and how they are altered by deterministic regulatory layers.

Amyloid-like Self-Assembly of a Cellular Compartment

Most vertebrate oocytes contain a Balbiani body, a large, non-membrane-bound compartment packed with RNA, mitochondria, and other organelles. Little is known about this compartment, though it specifies germline identity in many non-mammalian vertebrates. We show Xvelo, a disordered protein with an N-terminal prion-like domain, is an abundant constituent of Xenopus Balbiani bodies. Disruption of the prion-like domain of Xvelo, or substitution with a prion-like domain from an unrelated protein, interferes with its incorporation into Balbiani bodies in vivo. Recombinant Xvelo forms amyloid-like networks in vitro. Amyloid-like assemblies of Xvelo recruit both RNA and mitochondria in binding assays. We propose that Xenopus Balbiani bodies form by amyloid-like assembly of Xvelo, accompanied by co-recruitment of mitochondria and RNA. Prion-like domains are found in germ plasm organizing proteins in other species, suggesting that Balbiani body formation by amyloid-like assembly could be a conserved mechanism that helps oocytes function as long-lived germ cells.

The BioPlex Network: A Systematic Exploration of the Human Interactome

Protein interactions form a network whose structure drives cellular function and whose organization informs biological inquiry. Using high-throughput affinity-purification mass spectrometry, we identify interacting partners for 2,594 human proteins in HEK293T cells. The resulting network (BioPlex) contains 23,744 interactions among 7,668 proteins with 86% previously undocumented. BioPlex accurately depicts known complexes, attaining 80%-100% coverage for most CORUM complexes. The network readily subdivides into communities that correspond to complexes or clusters of functionally related proteins. More generally, network architecture reflects cellular localization, biological process, and molecular function, enabling functional characterization of thousands of proteins. Network structure also reveals associations among thousands of protein domains, suggesting a basis for examining structurally related proteins. Finally, BioPlex, in combination with other approaches, can be used to reveal interactions of biological or clinical significance. For example, mutations in the membrane protein VAPB implicated in familial amyotrophic lateral sclerosis perturb a defined community of interactors.

The Nuclear Proteome of a Vertebrate

The composition of the nucleoplasm determines the behavior of key processes such as transcription, yet there is still no reliable and quantitative resource of nuclear proteins. Furthermore, it is still unclear how the distinct nuclear and cytoplasmic compositions are maintained. To describe the nuclear proteome quantitatively, we isolated the large nuclei of frog oocytes via microdissection and measured the nucleocytoplasmic partitioning of ∼9,000 proteins by mass spectrometry. Most proteins localize entirely to either nucleus or cytoplasm; only ∼17% partition equally. A protein's native size in a complex, but not polypeptide molecular weight, is predictive of localization: partitioned proteins exhibit native sizes larger than ∼100 kDa, whereas natively smaller proteins are equidistributed. To evaluate the role of nuclear export in maintaining localization, we inhibited Exportin 1. This resulted in the expected re-localization of proteins toward the nucleus, but only 3% of the proteome was affected. Thus, complex assembly and passive retention, rather than continuous active transport, is the dominant mechanism for the maintenance of nuclear and cytoplasmic proteomes.

Size Scaling of Microtubule Assemblies in Early Xenopus Embryos

The first 12 cleavage divisions in Xenopus embryos provide a natural experiment in size scaling, as cell radius decreases ∼16-fold with little change in biochemistry. Analyzing both natural cleavage and egg extract partitioned into droplets revealed that mitotic spindle size scales with cell size, with an upper limit in very large cells. We discuss spindle-size scaling in the small- and large-cell regimes with a focus on the "limiting-component" hypotheses. Zygotes and early blastomeres show a scaling mismatch between spindle and cell size. This problem is solved, we argue, by interphase asters that act to position the spindle and transport chromosomes to the center of daughter cells. These tasks are executed by the spindle in smaller cells. We end by discussing possible mechanisms that limit mitotic aster size and promote interphase aster growth to cell-spanning dimensions.

Organization of early frog embryos by chemical waves emanating from centrosomes

The large cells in early vertebrate development face an extreme physical challenge in organizing their cytoplasm. For example, amphibian embryos have to divide cytoplasm that spans hundreds of micrometres every 30 min according to a precise geometry, a remarkable accomplishment given the extreme difference between molecular and cellular scales in this system. How do the biochemical reactions occurring at the molecular scale lead to this emergent behaviour of the cell as a whole? Based on recent findings, we propose that the centrosome plays a crucial role by initiating two autocatalytic reactions that travel across the large cytoplasm as chemical waves. Waves of mitotic entry and exit propagate out from centrosomes using the Cdk1 oscillator to coordinate the timing of cell division. Waves of microtubule-stimulated microtubule nucleation propagate out to assemble large asters that position spindles for the following mitosis and establish cleavage plane geometry. By initiating these chemical waves, the centrosome rapidly organizes the large cytoplasm during the short embryonic cell cycle, which would be impossible using more conventional mechanisms such as diffusion or nucleation by structural templating. Large embryo cells provide valuable insights to how cells control chemical waves, which may be a general principle for cytoplasmic organization.

Spatial organization of cytokinesis signaling reconstituted in a cell-free system

During animal cell division, the cleavage furrow is positioned by microtubules that signal to the actin cortex at the cell midplane. We developed a cell-free system to recapitulate cytokinesis signaling using cytoplasmic extract from Xenopus eggs. Microtubules grew out as asters from artificial centrosomes and met to organize antiparallel overlap zones. These zones blocked the interpenetration of neighboring asters and recruited cytokinesis midzone proteins, including the chromosomal passenger complex (CPC) and centralspindlin. The CPC was transported to overlap zones, which required two motor proteins, Kif4A and a Kif20A paralog. Using supported lipid bilayers to mimic the plasma membrane, we observed the recruitment of cleavage furrow markers, including an active RhoA reporter, at microtubule overlaps. This system opens further approaches to understanding the biophysics of cytokinesis signaling.

MultiNotch MS3 enables accurate, sensitive, and multiplexed detection of differential expression across cancer cell line proteomes

Multiplexed quantitation via isobaric chemical tags (e.g., tandem mass tags (TMT) and isobaric tags for relative and absolute quantitation (iTRAQ)) has the potential to revolutionize quantitative proteomics. However, until recently the utility of these tags was questionable due to reporter ion ratio distortion resulting from fragmentation of coisolated interfering species. These interfering signals can be negated through additional gas-phase manipulations (e.g., MS/MS/MS (MS3) and proton-transfer reactions (PTR)). These methods, however, have a significant sensitivity penalty. Using isolation waveforms with multiple frequency notches (i.e., synchronous precursor selection, SPS), we coisolated and cofragmented multiple MS2 fragment ions, thereby increasing the number of reporter ions in the MS3 spectrum 10-fold over the standard MS3 method (i.e., MultiNotch MS3). By increasing the reporter ion signals, this method improves the dynamic range of reporter ion quantitation, reduces reporter ion signal variance, and ultimately produces more high-quality quantitative measurements. To demonstrate utility, we analyzed biological triplicates of eight colon cancer cell lines using the MultiNotch MS3 method. Across all the replicates we quantified 8 378 proteins in union and 6 168 proteins in common. Taking into account that each of these quantified proteins contains eight distinct cell-line measurements, this data set encompasses 174 704 quantitative ratios each measured in triplicate across the biological replicates. Herein, we demonstrate that the MultiNotch MS3 method uniquely combines multiplexing capacity with quantitative sensitivity and accuracy, drastically increasing the informational value obtainable from proteomic experiments.

Accurate multiplexed proteomics at the MS2 level using the complement reporter ion cluster

Isobaric labeling strategies, such as isobaric tags for relative and absolute quantitation (iTRAQ) or tandem mass tags (TMT), have promised to dramatically increase the power of quantitative proteomics. However, when applied to complex mixtures, both the accuracy and precision are undermined by interfering peptide ions that coisolate and cofragment with the target peptide. Additional gas-phase isolation steps, such as proton-transfer ion-ion reactions (PTR) or higher-order MS3 scans, can almost completely eliminate this problem. Unfortunately, these methods come at the expense of decreased acquisition speed and sensitivity. Here we present a method that allows accurate quantification of TMT-labeled peptides at the MS2 level without additional ion purification. Quantification is based on the fragment ion cluster that carries most of the TMT mass balance. In contrast to the use of low m/z reporter ions, the localization of these complement TMT (TMT(C)) ions in the spectrum is precursor-specific; coeluting peptides do not generally affect the measurement of the TMT(C) ion cluster of interest. Unlike the PTR or MS3 strategies, this method can be implemented on a wide range of high-resolution mass spectrometers like the quadrupole Orbitrap instruments (QExactive). A current limitation of the method is that the efficiency of TMT(C) ion formation is affected by both peptide sequence and peptide ion charge state; we discuss potential routes to overcome this problem. Finally, we show that the complement reporter ion approach allows parallelization of multiplexed quantification and therefore holds the potential to multiply the number of distinct peptides that can be quantified in a given time frame.

Growth, interaction, and positioning of microtubule asters in extremely large vertebrate embryo cells

Ray Rappaport spent many years studying microtubule asters, and how they induce cleavage furrows. Here, we review recent progress on aster structure and dynamics in zygotes and early blastomeres of Xenopus laevis and Zebrafish, where cells are extremely large. Mitotic and interphase asters differ markedly in size, and only interphase asters span the cell. Growth of interphase asters occurs by a mechanism that allows microtubule density at the aster periphery to remain approximately constant as radius increases. We discuss models for aster growth, and favor a branching nucleation process. Neighboring asters that grow into each other interact to block further growth at the shared boundary. We compare the morphology of interaction zones formed between pairs of asters that grow out from the poles of the same mitotic spindle (sister asters) and between pairs not related by mitosis (non-sister asters) that meet following polyspermic fertilization. We argue growing asters recognize each other by interaction between antiparallel microtubules at the mutual boundary, and discuss models for molecular organization of interaction zones. Finally, we discuss models for how asters, and the centrosomes within them, are positioned by dynein-mediated pulling forces so as to generate stereotyped cleavage patterns. Studying these problems in extremely large cells is starting to reveal how general principles of cell organization scale with cell size.