Metabolic labeling with heavy water followed by LC-MS is a high throughput approach to study proteostasis in vivo. Advances in mass spectrometry and sample processing have allowed consistent detection of thousands of proteins at multiple time points. However, freely available automated bioinformatics tools to analyze and extract protein decay rate constants are lacking. Here, we describe d2ome—a robust, automated software solution for in vivo protein turnover analysis. d2ome is highly scalable, uses innovative approaches to nonlinear fitting, implements Grubbs’ outlier detection and removal, uses weighted-averaging of replicates, applies a data dependent elution time windowing, and uses mass accuracy in peak detection. Here, we discuss the application of d2ome in a comparative study of protein turnover in the livers of normal vs Western diet-fed LDLR–/– mice (mouse model of nonalcoholic fatty liver disease), which contained 256 LC-MS experiments. The study revealed reduced stability of 40S ribosomal protein subunits in the Western diet-fed mice.

Few of >150 published cell cycle modeling efforts use significant levels of data for tuning and validation. This reflects the difficultly to generate correlated quantitative data, and it points out a critical uncertainty in modeling efforts. To develop a data-driven model of cell cycle regulation, we used contiguous, dynamic measurements over two time scales (minutes and hours) calculated from static multiparametric cytometry data. The approach provided expression profiles of cyclin A2, cyclin B1, and phospho-S10-histone H3. The model was built by integrating and modifying two previously published models such that the model outputs for cyclins A and B fit cyclin expression measurements and the activation of B cyclin/Cdk1 coincided with phosphorylation of histone H3. The model depends on Cdh1-regulated cyclin degradation during G1, regulation of B cyclin/Cdk1 activity by cyclin A/Cdk via Wee1, and transcriptional control of the mitotic cyclins that reflects some of the current literature. We introduced autocatalytic transcription of E2F, E2F regulated transcription of cyclin B, Cdc20/Cdh1 mediated E2F degradation, enhanced transcription of mitotic cyclins during late S/early G2 phase, and the sustained synthesis of cyclin B during mitosis. These features produced a model with good correlation between state variable output and real measurements. Since the method of data generation is extensible, this model can be continually modified based on new correlated, quantitative data.

BackgroundAn imprecise quantitative sense for the oscillating levels of proteins and their modifications, interactions, and translocations as a function of the cell cycle is fundamentally important for a cartoon/narrative understanding for how the cell cycle works. Mathematical modeling of the same cartoon/narrative models would be greatly enhanced by an open-ended methodology providing precise quantification of many proteins and their modifications, etc. Here we present methodology that fulfills these features.MethodologyMultiparametric flow cytometry was performed on Molt4 cells to measure cyclins A2 and B1, phospho-S10-histone H3, DNA content, and light scatter (cell size). The resulting 5 dimensional data were analyzed as a series of bivariate plots to isolate the data as segments of an N-dimensional “worm” through the data space. Sequential, unidirectional regions of the data were used to assemble expression profiles for each parameter as a function of cell frequency.ResultsAnalysis of synthesized data in which the true values where known validated the approach. Triplicate experiments demonstrated exceptional reproducibility. Comparison of three triplicate experiments stained by two methods (single cyclin or dual cyclin measurements with common DNA and phospho-histone H3 measurements) supported the feasibility of combining an unlimited number of epitopes through this methodology. The sequential degradations of cyclin A2 followed by cyclin B1 followed by de-phosphorylation of histone H3 were precisely mapped. Finally, a two phase expression rate during interphase for each cyclin was robustly identified.ConclusionsVery precise, correlated expression profiles for important cell cycle regulating and regulated proteins and their modifications can be produced, limited only by the number of available high-quality antibodies. These profiles can be assembled into large information libraries for calibration and validation of mathematical models.

Advances in single-cell technologies have highlighted the prevalence and biological significance of cellular heterogeneity. A critical question is how to design experiments that faithfully capture the true range of heterogeneity from samples of cellular populations. Here, we develop a data-driven approach, illustrated in the context of image data, that estimates sampling depth required for prospective investigations of single-cell heterogeneity from an existing collection of samples.