From affinity selection to kinetic selection in Germinal Centre modelling.

Affinity maturation is an evolutionary process by which the affinity of antibodies (Abs) against specific antigens (Ags) increases through rounds of B-cell proliferation, somatic hypermutation, and positive selection in germinal centres (GC). The positive selection of B cells depends on affinity, but the underlying mechanisms of affinity discrimination and affinity-based selection are not well understood. It has been suggested that selection in GC depends on both rapid binding of B-cell receptors (BcRs) to Ags which is kinetically favourable and tight binding of BcRs to Ags, which is thermodynamically favourable; however, it has not been shown whether a selection bias for kinetic properties is present in the GC. To investigate the GC selection bias towards rapid and tight binding, we developed an agent-based model of GC and compared the evolution of founder B cells with initially identical low affinities but with different association/dissociation rates for Ag presented by follicular dendritic cells in three Ag collection mechanisms. We compared an Ag collection mechanism based on association/dissociation rates of B-cell interaction with presented Ag, which includes a probabilistic rupture of bonds between the B-cell and Ag (Scenario-1) with a reference scenario based on an affinity-based Ag collection mechanism (Scenario-0). Simulations showed that the mechanism of Ag collection affects the GC dynamics and the GC outputs concerning fast/slow (un)binding of B cells to FDC-presented Ags. In particular, clones with lower dissociation rates outcompete clones with higher association rates in Scenario-1, while remaining B cells from clones with higher association rates reach higher affinities. Accordingly, plasma cell and memory B cell populations were biased towards B-cell clones with lower dissociation rates. Without such probabilistic ruptures during the Ag extraction process (Scenario-2), the selective advantage for clones with very low dissociation rates diminished, and the affinity maturation level of all clones decreased to the reference level.

Human Blood Lipoprotein Predictions from 1H NMR Spectra: Protocol, Model Performances, and Cage of Covariance.

Lipoprotein subfractions are biomarkers for the early diagnosis of cardiovascular diseases. The reference method, ultracentrifugation, for measuring lipoproteins is time-consuming, and there is a need to develop a rapid method for cohort screenings. This study presents partial least-squares regression models developed using 1H nuclear magnetic resonance (NMR) spectra and concentrations of lipoproteins as measured by ultracentrifugation on 316 healthy Danes. This study explores, for the first time, different regions of the 1H NMR spectrum representing signals of molecules in lipoprotein particles and different lipid species to develop parsimonious, reliable, and optimal prediction models. A total of 65 lipoprotein main and subfractions were predictable with high accuracy, Q2 of >0.6, using an optimal spectral region (1.4-0.6 ppm) containing methylene and methyl signals from lipids. The models were subsequently tested on an independent cohort of 290 healthy Swedes with predicted and reference values matching by up to 85-95%. In addition, an open software tool was developed to predict lipoproteins concentrations in human blood from standardized 1H NMR spectral recordings.

Training for translation between disciplines: a philosophy for life and data sciences curricula.

Motivation: Our society has become data-rich to the extent that research in many areas has become impossible without computational approaches. Educational programmes seem to be lagging behind this development. At the same time, there is a growing need not only for strong data science skills, but foremost for the ability to both translate between tools and methods on the one hand, and application and problems on the other. Results: Here we present our experiences with shaping and running a masters' programme in bioinformatics and systems biology in Amsterdam. From this, we have developed a comprehensive philosophy on how translation in training may be achieved in a dynamic and multidisciplinary research area, which is described here. We furthermore describe two requirements that enable translation, which we have found to be crucial: sufficient depth and focus on multidisciplinary topic areas, coupled with a balanced breadth from adjacent disciplines. Finally, we present concrete suggestions on how this may be implemented in practice, which may be relevant for the effectiveness of life science and data science curricula in general, and of particular interest to those who are in the process of setting up such curricula. Supplementary information: Supplementary data are available at Bioinformatics online.

Repeatability and reproducibility of lipoprotein particle profile measurements in plasma samples by ultracentrifugation.

Background Characterization of lipoprotein particle profiles (LPPs) (including main classes and subclasses) by means of ultracentrifugation (UC) is highly requested given its clinical potential. However, rapid methods are required to replace the very labor-intensive UC method and one solution is to calibrate rapid nuclear magnetic resonance (NMR)-based prediction models, but the reliability of the UC-response method required for the NMR calibration has been largely overlooked. Methods This study provides a comprehensive repeatability and reproducibility study of various UC-based lipid measurements (cholesterol, triglycerides [TGs], free cholesterol, phospholipids, apolipoprotein [apo]A1 and apoB) in different main classes and subclasses of 25 duplicated fresh plasma samples and of 42 quality control (QC) frozen pooled plasma samples of healthy individuals. Results Cholesterol, apoA1 and apoB measurements were very repeatable in all classes (intraclass correlation coefficient [ICC]: 92.93%-99.54%). Free cholesterol and phospholipid concentrations in main classes and subclasses and TG concentrations in high-density lipoproteins (HDL), HDL subclasses and low-density lipoproteins (LDL) subclasses, showed worse repeatability (ICC: 19.21%-99.08%) attributable to low concentrations, variability introduced during UC and assay limitations. On frozen QC samples, the reproducibility of cholesterol, apoA1 and apoB concentrations was found to be better than for the free cholesterol, phospholipids and TGs concentrations. Conclusions This study shows that for LPPs measurements near or below the limit of detection (LOD) in some of the subclasses, as well as the use of frozen samples, results in worsened repeatability and reproducibility. Furthermore, we show that the analytical assay coupled to UC for free cholesterol and phospholipids have different repeatability and reproducibility. All of this needs to be taken into account when calibrating future NMR-based models.

Stoichiometry, Absolute Abundance, and Localization of Proteins in the Bacillus cereus Spore Coat Insoluble Fraction Determined Using a QconCAT Approach.

Spores of Bacillus cereus pose a threat to food safety due to their high resistance to the heat or acid treatments commonly used to make food microbiologically safe. Spores may survive these treatments and later resume growth either on foodstuffs or, after ingestion, upon entering the gut they are capable of producing toxins, which cause either vomiting or diarrhea. The outer layers of the spore, the spore coat and exosporium, consist primarily of proteins that may serve as potential biomarkers for detection. The major morphogenetic protein CotE is important for correct assembly and attachment of the outermost layer, the exosporium, and by extension retention of many proteins. However, characterization of the proteins affected by deletion of CotE has been limited to electrophoretic patterns. Here we report the effect of CotE deletion on the insoluble fraction of the spore proteome through liquid chromatography-Fourier transform tandem mass spectrometry (LC-FTMS/MS) analysis. A total of 560 proteins have been identified in both mutant and wild-type spore coat isolates. A further 163 proteins were identified exclusively in wild-type spore isolates indicating that they are dependent on CotE for their association with the spore. Several of these are newly confirmed as associated with the exosporium, namely BC_2569 (BclF), BC_3345, BC_2427, BC_2878, BC_0666, BC_2984, BC_3481, and BC_2570. A total of 153 proteins were only identified in ΔCotE spore isolates. This was observed for proteins that are known or likely to be interacting with or are encased by CotE. Crucial spore proteins were quantified using a QconCAT reference standard, the first time this was used in a biochemically heterogeneous system. This allowed us to determine the absolute abundance of 21 proteins, which spanned across three orders of magnitude and together covered 5.66% ± 0.51 of the total spore weight. Applying the QconCAT methodology to the ΔCotE mutant allowed us to quantify 4.13% ± 0.14 of the spore total weight and revealed a reduction in abundance for most known exosporium associated proteins upon CotE deletion. In contrast, several proteins, either known or likely to be interacting with or encased by CotE (i.e., GerQ), were more abundant. The results obtained provide deeper insight into the layered spore structure such as which proteins are exposed on the outside of the spore. This information is important for developing detection methods for targeting spores in a food safety setting. Furthermore, protein stoichiometry and determination of the abundance of germination mediating enzymes provides useful information for germination and outgrowth model development.

Global testing of shifts in metabolic phenotype.

INTRODUCTION: Current metabolomics approaches to unravel impact of diet- or lifestyle induced phenotype variation and shifts predominantly deploy univariate or multivariate approaches, with a posteriori interpretation at pathway level. This however often provides only a fragmented view on the involved metabolic pathways. OBJECTIVES: To demonstrate the feasibility of using Goeman's global test (GGT) for assessment of variation and shifts in metabolic phenotype at the level of a priori defined pathways. METHODS: Two intervention studies with identified phenotype variations and shifts were examined. In a weight loss (WL) intervention study obese subjects received a mixed meal challenge before and after WL. In a polyphenol (PP) intervention study obese subjects received a high fat mixed meal challenge (61E% fat) before and after a PP intervention. Plasma samples were obtained at fasting and during the postprandial response. Besides WL- and PP-induced phenotype shifts, also correlation of plasma metabolome with phenotype descriptors was assessed at pathway level. The plasma metabolome covered organic acids, amino acids, biogenic amines, acylcarnitines and oxylipins. RESULTS: For the population of the WL study, GGT revealed that HOMA correlated with the fasting levels of the TCA cycle, BCAA catabolism, the lactate, arginine-proline and phenylalanine-tyrosine pathways. For the population of the PP study, HOMA correlated with fasting metabolite levels of TCA cycle, fatty acid oxidation and phenylalanine-tyrosine pathways. These correlations were more pronounced for metabolic pathways in the fasting state, than during the postprandial response. The effect of the WL and PP intervention on a priori defined metabolic pathways, and correlation of pathways with insulin sensitivity as described by HOMA was in line with previous studies. CONCLUSION: GGT confirmed earlier biological findings in a hypothesis led approach. A main advantage of GGT is that it provides a direct view on involvement of a priori defined pathways in phenotype shifts.

Toward Reliable Lipoprotein Particle Predictions from NMR Spectra of Human Blood: An Interlaboratory Ring Test.

Lipoprotein profiling of human blood by 1H nuclear magnetic resonance (NMR) spectroscopy is a rapid and promising approach to monitor health and disease states in medicine and nutrition. However, lack of standardization of measurement protocols has prevented the use of NMR-based lipoprotein profiling in metastudies. In this study, a standardized NMR measurement protocol was applied in a ring test performed across three different laboratories in Europe on plasma and serum samples from 28 individuals. Data was evaluated in terms of (i) spectral differences, (ii) differences in LPD predictions obtained using an existing prediction model, and (iii) agreement of predictions with cholesterol concentrations in high- and low-density lipoproteins (HDL and LDL) particles measured by standardized clinical assays. ANOVA-simultaneous component analysis (ASCA) of the ring test spectral ensemble that contains methylene and methyl peaks (1.4-0.6 ppm) showed that 97.99% of the variance in the data is related to subject, 1.62% to sample type (serum or plasma), and 0.39% to laboratory. This interlaboratory variation is in fact smaller than the maximum acceptable intralaboratory variation on quality control samples. It is also shown that the reproducibility between laboratories is good enough for the LPD predictions to be exchangeable when the standardized NMR measurement protocol is followed. With the successful implementation of this protocol, which results in reproducible prediction of lipoprotein distributions across laboratories, a step is taken toward bringing NMR more into scope of prognostic and diagnostic biomarkers, reducing the need for less efficient methods such as ultracentrifugation or high-performance liquid chromatography (HPLC).

Time-series analysis of the transcriptome and proteome of Escherichia coli upon glucose repression.

Time-series transcript- and protein-profiles were measured upon initiation of carbon catabolite repression in Escherichia coli, in order to investigate the extent of post-transcriptional control in this prototypical response. A glucose-limited chemostat culture was used as the CCR-free reference condition. Stopping the pump and simultaneously adding a pulse of glucose, that saturated the cells for at least 1h, was used to initiate the glucose response. Samples were collected and subjected to quantitative time-series analysis of both the transcriptome (using microarray analysis) and the proteome (through a combination of 15N-metabolic labeling and mass spectrometry). Changes in the transcriptome and corresponding proteome were analyzed using statistical procedures designed specifically for time-series data. By comparison of the two sets of data, a total of 96 genes were identified that are post-transcriptionally regulated. This gene list provides candidates for future in-depth investigation of the molecular mechanisms involved in post-transcriptional regulation during carbon catabolite repression in E. coli, like the involvement of small RNAs.

Bioinformatics and systems biology: bridging the gap between heterogeneous student backgrounds.

Teaching students with very diverse backgrounds can be extremely challenging. This article uses the Bioinformatics and Systems Biology MSc in Amsterdam as a case study to describe how the knowledge gap for students with heterogeneous backgrounds can be bridged. We show that a mix in backgrounds can be turned into an advantage by creating a stimulating learning environment for the students. In the MSc Programme, conversion classes help to bridge differences between students, by mending initial knowledge and skill gaps. Mixing students from different backgrounds in a group to solve a complex task creates an opportunity for the students to reflect on their own abilities. We explain how a truly interdisciplinary approach to teaching helps students of all backgrounds to achieve the MSc end terms. Moreover, transferable skills obtained by the students in such a mixed study environment are invaluable for their later careers.

Capturing Arabidopsis root architecture dynamics with ROOT-FIT reveals diversity in responses to salinity.

The plant root is the first organ to encounter salinity stress, but the effect of salinity on root system architecture (RSA) remains elusive. Both the reduction in main root (MR) elongation and the redistribution of the root mass between MRs and lateral roots (LRs) are likely to play crucial roles in water extraction efficiency and ion exclusion. To establish which RSA parameters are responsive to salt stress, we performed a detailed time course experiment in which Arabidopsis (Arabidopsis thaliana) seedlings were grown on agar plates under different salt stress conditions. We captured RSA dynamics with quadratic growth functions (root-fit) and summarized the salt-induced differences in RSA dynamics in three growth parameters: MR elongation, average LR elongation, and increase in number of LRs. In the ecotype Columbia-0 accession of Arabidopsis, salt stress affected MR elongation more severely than LR elongation and an increase in LRs, leading to a significantly altered RSA. By quantifying RSA dynamics of 31 different Arabidopsis accessions in control and mild salt stress conditions, different strategies for regulation of MR and LR meristems and root branching were revealed. Different RSA strategies partially correlated with natural variation in abscisic acid sensitivity and different Na(+)/K(+) ratios in shoots of seedlings grown under mild salt stress. Applying root-fit to describe the dynamics of RSA allowed us to uncover the natural diversity in root morphology and cluster it into four response types that otherwise would have been overlooked.

A critical assessment of feature selection methods for biomarker discovery in clinical proteomics.

In this paper, we compare the performance of six different feature selection methods for LC-MS-based proteomics and metabolomics biomarker discovery-t test, the Mann-Whitney-Wilcoxon test (mww test), nearest shrunken centroid (NSC), linear support vector machine-recursive features elimination (SVM-RFE), principal component discriminant analysis (PCDA), and partial least squares discriminant analysis (PLSDA)-using human urine and porcine cerebrospinal fluid samples that were spiked with a range of peptides at different concentration levels. The ideal feature selection method should select the complete list of discriminating features that are related to the spiked peptides without selecting unrelated features. Whereas many studies have to rely on classification error to judge the reliability of the selected biomarker candidates, we assessed the accuracy of selection directly from the list of spiked peptides. The feature selection methods were applied to data sets with different sample sizes and extents of sample class separation determined by the concentration level of spiked compounds. For each feature selection method and data set, the performance for selecting a set of features related to spiked compounds was assessed using the harmonic mean of the recall and the precision (f-score) and the geometric mean of the recall and the true negative rate (g-score). We conclude that the univariate t test and the mww test with multiple testing corrections are not applicable to data sets with small sample sizes (n = 6), but their performance improves markedly with increasing sample size up to a point (n > 12) at which they outperform the other methods. PCDA and PLSDA select small feature sets with high precision but miss many true positive features related to the spiked peptides. NSC strikes a reasonable compromise between recall and precision for all data sets independent of spiking level and number of samples. Linear SVM-RFE performs poorly for selecting features related to the spiked compounds, even though the classification error is relatively low.

Generic framework for high-dimensional fixed-effects ANOVA.

In functional genomics it is more rule than exception that experimental designs are used to generate the data. The samples of the resulting data sets are thus organized according to this design and for each sample many biochemical compounds are measured, e.g. typically thousands of gene-expressions or hundreds of metabolites. This results in high-dimensional data sets with an underlying experimental design. Several methods have recently become available for analyzing such data while utilizing the underlying design. We review these methods by putting them in a unifying and general framework to facilitate understanding the (dis-)similarities between the methods. The biological question dictates which method to use and the framework allows for building new methods to accommodate a range of such biological questions. The framework is built on well known fixed-effect ANOVA models and subsequent dimension reduction. We present the framework both in matrix algebra as well as in more insightful geometrical terms. We show the workings of the different special cases of our framework with a real-life metabolomics example from nutritional research and a gene-expression example from the field of virology.

A critical assessment of SELDI-TOF-MS for biomarker discovery in serum and tissue of patients with an ovarian mass.

BACKGROUND: Less than 25% of patients with a pelvic mass who are presented to a gynecologist will eventually be diagnosed with epithelial ovarian cancer. Since there is no reliable test to differentiate between different ovarian tumors, accurate classification could facilitate adequate referral to a gynecological oncologist, improving survival. The goal of our study was to assess the potential value of a SELDI-TOF-MS based classifier for discriminating between patients with a pelvic mass. METHODS: Our study design included a well-defined patient population, stringent protocols and an independent validation cohort. We compared serum samples of 53 ovarian cancer patients, 18 patients with tumors of low malignant potential, and 57 patients with a benign ovarian tumor on different ProteinChip arrays. In addition, from a subset of 84 patients, tumor tissues were collected and microdissection was used to isolate a pure and homogenous cell population. RESULTS: Diagonal Linear Discriminant Analysis (DLDA) and Support Vector Machine (SVM) classification on serum samples comparing cancer versus benign tumors, yielded models with a classification accuracy of 71-81% (cross-validation), and 73-81% on the independent validation set. Cancer and benign tissues could be classified with 95-99% accuracy using cross-validation. Tumors of low malignant potential showed protein expression patterns different from both benign and cancer tissues. Remarkably, none of the peaks differentially expressed in serum samples were found to be differentially expressed in the tissue lysates of those same groups. CONCLUSION: Although SELDI-TOF-MS can produce reliable classification results in serum samples of ovarian cancer patients, it will not be applicable in routine patient care. On the other hand, protein profiling of microdissected tumor tissue may lead to a better understanding of oncogenesis and could still be a source of new serum biomarkers leading to novel methods for differentiating between different histological subtypes.

To aggregate or not to aggregate high-dimensional classifiers.

BACKGROUND: High-throughput functional genomics technologies generate large amount of data with hundreds or thousands of measurements per sample. The number of sample is usually much smaller in the order of ten or hundred. This poses statistical challenges and calls for appropriate solutions for the analysis of this kind of data. RESULTS: Principal component discriminant analysis (PCDA), an adaptation of classical linear discriminant analysis (LDA) for high-dimensional data, has been selected as an example of a base learner. The multiple versions of PCDA models from repeated double cross-validation were aggregated, and the final classification was performed by majority voting. The performance of this approach was evaluated by simulation, genomics, proteomics and metabolomics data sets. CONCLUSIONS: The aggregating PCDA learner can improve the prediction performance, provide more stable result, and help to know the variability of the models. The disadvantage and limitations of aggregating were also discussed.

Time alignment algorithms based on selected mass traces for complex LC-MS data.

Time alignment of complex LC-MS data remains a challenge in proteomics and metabolomics studies. This work describes modifications of the Dynamic Time Warping (DTW) and the Parametric Time Warping (PTW) algorithms that improve the alignment quality for complex, highly variable LC-MS data sets. Regular DTW or PTW use one-dimensional profiles such as the Total Ion Chromatogram (TIC) or Base Peak Chromatogram (BPC) resulting in correct alignment if the signals have a relatively simple structure. However, when aligning the TICs of chromatograms from complex mixtures with large concentration variability such as serum or urine, both algorithms often lead to misalignment of peaks and thus incorrect comparisons in the subsequent statistical analysis. This is mainly due to the fact that compounds with different m/z values but similar retention times are not considered separately but confounded in the benefit function of the algorithms using only one-dimensional information. Thus, it is necessary to treat the information of different mass traces separately in the warping function to ensure that compounds having the same m/z value and retention time are aligned to each other. The Component Detection Algorithm (CODA) is widely used to calculate the quality of an LC-MS mass trace. By combining CODA with the warping algorithms of DTW or PTW (DTW-CODA or PTW-CODA), we include only high quality mass traces measured by CODA in the benefit function. Our results show that using several CODA selected high quality mass traces in DTW-CODA and PTW-CODA significantly improves the alignment quality of three different, highly complex LC-MS data sets. Moreover, DTW-CODA leads to better preservation of peak shape as compared to the original DTW-TIC algorithm, which often suffers from a substantial peak shape distortion. Our results show that combination of CODA selected mass traces with different time alignment algorithm is a general principle that provide accurate alignment for highly complex samples with large concentration variability.

Endocrine pulse identification using penalized methods and a minimum set of assumptions.

The detection of hormone secretion episodes is important for understanding normal and abnormal endocrine functioning, but pulse identification from hormones measured with short interval sampling is challenging. Furthermore, to obtain useable results, the model underlying hormone secretion and clearance must be augmented with restrictions based on biologically acceptable assumptions. Here, using the assumption that there are only a few time points at which a hormone is secreted, we used a modern penalized nonlinear least-squares setup to select the number of secretion events. We did not assume a particular shape or frequency distribution for the secretion pulses. Our pulse identfication method, VisPulse, worked well with luteinizing hormone (LH), cortisol, growth hormone, or testosterone. In particular, applying our modeling strategy to previous LH data revealed a good correlation between the modeled and measured LH hormone concentrations, the estimated secretion pattern was sparse, and the small and structureless residuals indicated a proper model with a good fit. We benchmarked our method to AutoDecon, a commonly used hormone secretion model, and performed releasing hormone infusion experiments. The results of these experiments confirmed that our method is accurate and outperforms AutoDecon, especially for detecting silent periods and small secretion events, suggesting a high-secretion event resolution. Method validation using (releasing hormone) infusion data revealed sensitivities and selectivities of 0.88 and 0.95 and of 0.69 and 0.91 for VisPulse and AutoDecon, respectively.

The photographer and the greenhouse: how to analyse plant metabolomics data.

INTRODUCTION: Plant metabolomics experiments yield large amounts of data, too much to be interpretable by eye. Multivariate data analyses are therefore essential to extract and visualise the information of interest. OBJECTIVE: Because multivariate statistical methods may be remote from the expertise of many scientists working in the metabolomics field, this overview provides a step-by-step description of a multivariate data analysis, starting from the experiment and ending with the figures appearing in scientific journals. METHODOLOGY: We developed a thought experiment that explores the relationship between the differences in nutrient levels and three plant developmental descriptors through photography of the greenhouse they grow in. Through this, multivariate data analysis, data preprocessing and model validation are illustrated. Finally some of the presented methods are illustrated by the analysis of a plant metabolomics dataset. CONCLUSION: This paper will familiarize non-specialised researchers with the main concepts in multivariate data analysis and allow them to develop and evaluate metabolomic data analyses more critically.

Crossfit analysis: a novel method to characterize the dynamics of induced plant responses.

BACKGROUND: Many plant species show induced responses that protect them against exogenous attacks. These responses involve the production of many different bioactive compounds. Plant species belonging to the Brassicaceae family produce defensive glucosinolates, which may greatly influence their favorable nutritional properties for humans. Each responding compound may have its own dynamic profile and metabolic relationships with other compounds. The chemical background of the induced response is therefore highly complex and may therefore not reveal all the properties of the response in any single model. RESULTS: This study therefore aims to describe the dynamics of the glucosinolate response, measured at three time points after induction in a feral Brassica, by a three-faceted approach, based on Principal Component Analysis. First the large-scale aspects of the response are described in a 'global model' and then each time-point in the experiment is individually described in 'local models' that focus on phenomena that occur at specific moments in time. Although each local model describes the variation among the plants at one time-point as well as possible, the response dynamics are lost. Therefore a novel method called the 'Crossfit' is described that links the local models of different time-points to each other. CONCLUSIONS: Each element of the described analysis approach reveals different aspects of the response. The crossfit shows that smaller dynamic changes may occur in the response that are overlooked by global models, as illustrated by the analysis of a metabolic profiling dataset of the same samples.

Improved cholesterol phenotype analysis by a model relating lipoprotein life cycle processes to particle size.

Increased plasma cholesterol is a known risk factor for cardiovascular disease. Lipoprotein particles transport both cholesterol and triglycerides through the blood. It is thought that the size distribution of these particles codetermines cardiovascular disease risk. New types of measurements can determine the concentration of many lipoprotein size-classes but exactly how each small class relates to disease risk is difficult to clear up. Because relating physiological process status to disease risk seems promising, we propose investigating how lipoprotein production, lipolysis, and uptake processes depend on particle size. To do this, we introduced a novel model framework (Particle Profiler) and evaluated its feasibility. The framework was tested using existing stable isotope flux data. The model framework implementation we present here reproduced the flux data and derived lipoprotein size pattern changes that corresponded to measured changes. It also sensitively indicated changes in lipoprotein metabolism between patient groups that are biologically plausible. Finally, the model was able to reproduce the cholesterol and triglyceride phenotype of known genetic diseases like familial hypercholesterolemia and familial hyperchylomicronemia. In the future, Particle Profiler can be applied for analyzing detailed lipoprotein size profile data and deriving rates of various lipolysis and uptake processes if an independent production estimate is given.