ODE constrained mixture modelling: a method for unraveling subpopulation structures and dynamics.

Area of Science:

• Systems Biology
• Computational Biology
• Cellular Dynamics

Background:

• Cell-to-cell variability is common in biological systems, with genetically identical cells responding differently to stimuli.
• Existing methods for analyzing heterogeneous cell populations have limitations, including inability to handle multiple experimental conditions or incorporate prior biological knowledge.

Purpose of the Study:

• To develop a novel computational method that integrates ordinary differential equation (ODE) models with mixture models to analyze cell-to-cell variability.
• To overcome the limitations of existing methods by enabling simultaneous analysis across experimental conditions and incorporating biological information.

Main Methods:

• Combined ordinary differential equation (ODE) models with mixture models to create ODE constrained mixture models.
• Utilized simulation studies to validate the method's ability to identify subpopulation structures and sources of variability.
• Applied the method to study NGF-induced Erk1/2 phosphorylation in primary sensory neurons.

Main Results:

• ODE constrained mixture models successfully unravel subpopulation structures and identify sources of cell-to-cell variability.
• The method provides reliable estimates for kinetic rates and subpopulation characteristics.
• Demonstrated the model's capability to reconstruct static and dynamic subpopulation characteristics across different experimental conditions for NGF-induced Erk1/2 phosphorylation.

Conclusions:

• ODE constrained mixture models offer a powerful approach to analyze complex biological variability and uncover mechanistic insights.
• The method exhibits high sensitivity and experimental validation confirms its capabilities.
• This approach enhances understanding of cellular processes, such as pain-related signaling pathways.