Abstract:Background: Multi-cellular organisms as well as bacteria are faced with diverse, ever changing environments. To ensure survival and evolutionary success, bacterial populations exploit cell-to-cell variability, which increase the robustness against environmental changes. More complex organisms, such as mammals, evolved strategies to actively detect and respond to environmental changes. The building blocks for the necessary structures and functional units are cell types with distinct properties. Several statistical methods and modeling approaches have been introduced to analyze such structured populations, they are however either incapable of simultaneously analyzing different experimental conditions, computationally demanding and difficult to apply, or neglect information available in the literature. Methods: To improve upon existing methods, we introduced ordinary differential equation (ODE) constrained mixture models. This hybrid modeling approach integrates ideas from statistics, namely mixture modeling, and mechanistic dynamical modeling. The mechanistic model enables the simultaneous analysis of multiple experimental conditions while the mixture components can exploit distributional information contained in the data. The fitting of different ODE constrained mixture models, representing different model hypotheses, to the available data allows for the assessment of the origin of cell-to-cell variability as well as subpopulation structures and dynamics. To assess the proposed method we performed different simulation studies. Furthermore, ODE constrained mixture modeling has been used to study NGF-induced Erk1/2 phosphorylation in primary sensory neurons, a process relevant in inflammatory and neuropathic pain. We developed a simple pathway model for NGF-induced Erk1/2 phosphorylation and compared different hypothesis regarding the subpopulation structure and dynamics using the available time and dose response data. ODE constrained mixture modeling suggested that the population consists of two subpopulation with different NGF receptor abundances, which has been validated using additional experiments. Discussion and conclusions: The ODE constrained mixture modeling introduced in this work bears great potential for the analysis of heterogeneous cell populations. The proposed method allows for a simple but pervasive analysis of heterogeneous cell systems and more profound, mechanistic insights. ODE constrained mixture models facilitate the reconstruction of subpopulation sizes and dynamics even in situations, where the subpopulations are hardly distinguishable. This is shown for simulation examples as well as for the process of NGF-induced Erk1/2 phosphorylation in primary sensory neurons. ODE constrained mixture modeling is also highly flexible as the subpopulation dynamics can either be deterministic or stochastic. This only requires alterations in the ODE part, namely reaction rate equations or moment equations are used to model the dynamics.

Nonlinear Mixed‐Effects Model Development and Simulation Using nlmixr and Related R Open‐Source Packages

Linear mixed effects models for non‐Gaussian continuous repeated measurement data

NLMEModeling: A Wolfram Mathematica Package for Nonlinear Mixed Effects Modeling of Dynamical Systems

Multi-loop nonlinear internal model controller design under nonlinear dynamic PLS framework using ARX-neural network model

Extending the code in the open-source saemix package to fit joint models of longitudinal and time-to-event data

Low‐dimensional neural ordinary differential equations accounting for inter‐individual variability implemented in Monolix and NONMEM

Marginalization in nonlinear mixed-effects models with an application to dose-response analysis

Embedding R inside the PhysPK Bio-simulation Software for Pharmacokinetics Population Analysis

Approximate Inferences for Nonlinear Mixed Effects Models with Scale Mixtures of Skew-Normal Distributions

Linear Mixed-Effects Models for Non-Gaussian Repeated Measurement Data

Estimation of nonlinear mixed-effects continuous-time models using the continuous-discrete extended Kalman filter

A fully automatic tool for development of population pharmacokinetic models

Fast selection of nonlinear mixed effect models using penalized likelihood

Unraveling Subpopulation Structures and Dynamics with ODE Constrained Mixture Modeling

A Finite Mixture of Nonlinear Random Coefficient Models for Continuous Repeated Measures Data

An amortized approach to non-linear mixed-effects modeling based on neural posterior estimation

Use of nonlinear mixed effect modeling for the meta-analysis of preclinical pharmacokinetic data: application to S 20342 in the rat

Development and comparison of model‐integrated evidence approaches for bioequivalence studies with pharmacokinetic end points

CPhaMAS: An Online Platform for Pharmacokinetic Data Analysis Based on Optimized Parameter Fitting Algorithm

Using Fisher Information Matrix to predict uncertainty in covariate effects and power to detect their relevance in Non-Linear Mixed Effect Models in pharmacometrics

PC-based differential model fitting as a support for clinical research