Abstract:Kinetic models of metabolism can be constructed to predict cellular regulation and devise metabolic engineering strategies, and various promising computational workflows have been developed in recent years for this. Due to the uncertainty in the kinetic parameter values required to build kinetic models, these workflows rely on ensemble modeling (EM) principles for sampling and building populations of models describing observed physiologies. Sensitivity coefficients from metabolic control analysis (MCA) of kinetic models can provide important insight about cellular control around a given physiological steady state. However, despite considering populations of kinetic models and their model outputs, current approaches do not provide adequate tools for statistical inference. To derive conclusions from model outputs, such as MCA sensitivity coefficients, it is necessary to rank/compare populations of variables with each other. Currently existing workflows consider confidence intervals (CIs) that are derived independently for each comparable variable. Hence, it is important to derive simultaneous CIs for the variables that we wish to rank/compare. Herein, we used an existing large-scale kinetic model of Escherichia Coli metabolism to present how univariate CIs can lead to incorrect conclusions, and we present a new workflow that applies three different multivariate statistical approaches. We use the Bonferroni and the exact normal methods to build symmetric CIs using the normality assumptions. We then suggest how bootstrapping can compute asymmetric CIs whilst relaxing this normality assumption. We conclude that the Bonferroni and the exact normal methods can provide simple and efficient ways for constructing reliable CIs, with the exact normal method favored over the Bonferroni when the compared variables present dependencies. Bootstrapping, despite its significantly higher computational cost, is recommended when comparing non-normal distributions of variables. Additionally, we show how the Bonferroni method can readily be used to estimate required sample numbers to attain a certain CI size.

Estimating the structural distribution function of cell probabilities

Estimating a structural distribution function by grouping

Approximate Bayesian Computation for Stochastic Single-Cell Time-Lapse Data Using Multivariate Test Statistics

Dynamic Density Estimation in Heterogeneous Cell Populations

Estimation of the lifetime distribution from fluctuations in Bellman-Harris processes

Accounting for randomness in measurement and sampling in studying cancer cell population dynamics

Kernel Estimation of Density Level Sets

Likelihood-based inference, identifiability and prediction using count data from lattice-based random walk models

Accounting for Randomness in Measurement and Sampling in Study of Cancer Cell Population Dynamics

Estimating the Mean of a Non-Linear Function of Conditional Expectation

Probability Density Function Estimation Using Gamma Kernels

Bayesian Computation for Contingency Tables with Incomplete Cell-Counts

Density-based modeling and identification of biochemical networks in cell populations

Estimation of population-level summaries in general semiparametric repeated measures regression models

Estimating Unknown Population Sizes Using the Hypergeometric Distribution

Estimating several survival functions under uniform stochastic ordering

A multi-stage model of cell proliferation and death: tracking cell divisions with Erlang distributions

Bayesian Estimation of Multinomial Cell Probabilities Incorporating Information from Aggregated Observations

Estimation of Structural Reliability for a Nested Multivariate Exponential Distributed Stress

Statistical inference in ensemble modeling of cellular metabolism

Estimations by stable motions and applications