A data-informed mathematical model of microglial cell dynamics during ischemic stroke in the middle cerebral artery

Sara Amato,Andrea Arnold

2024-03-22

Abstract:Neuroinflammation immediately follows the onset of ischemic stroke in the middle cerebral artery. During this process, microglial cells are activated in and recruited to the penumbra. Microglial cells can be activated into two different phenotypes: M1, which can worsen brain injury; or M2, which can aid in long-term recovery. In this study, we contribute a summary of experimental data on microglial cell counts in the penumbra following ischemic stroke induced by middle cerebral artery occlusion (MCAO) in mice and compile available data sets into a single set suitable for time series analysis. Further, we formulate a mathematical model of microglial cells in the penumbra during ischemic stroke due to MCAO. Through use of global sensitivity analysis and Markov Chain Monte Carlo (MCMC)-based parameter estimation, we analyze the effects of the model parameters on the number of M1 and M2 cells in the penumbra and fit identifiable parameters to the compiled experimental data set. We utilize results from MCMC parameter estimation to ascertain uncertainty bounds and forward predictions for the number of M1 and M2 microglial cells over time. Results demonstrate the significance of parameters related to M1 and M2 activation on the number of M1 and M2 microglial cells. Simulations further suggest that potential outliers in the observed data may be omitted and forecast predictions suggest a lingering inflammatory response.

Cell Behavior,Quantitative Methods,Applications

What problem does this paper attempt to address?

The paper attempts to address the dynamic changes of microglia in the penumbra during ischemic stroke induced by middle cerebral artery occlusion (MCAO). Specifically, the study aims to:

1. **Summarize Experimental Data**: Collect and organize experimental data on the number of microglia in the penumbra after ischemic stroke in MCAO-induced mouse models, forming a dataset suitable for time series analysis.
2. **Establish a Mathematical Model**: Construct a mathematical model based on experimental data to describe the dynamic behavior of microglia in the penumbra during ischemic stroke induced by MCAO.
3. **Parameter Estimation and Sensitivity Analysis**: Estimate model parameters through global sensitivity analysis and Markov Chain Monte Carlo (MCMC) methods, and analyze the impact of these parameters on the number of M1 and M2 microglia.
4. **Prediction and Uncertainty Quantification**: Use MCMC parameter estimation results to determine the uncertainty range of parameters and perform forward predictions to assess the changes in the number of M1 and M2 microglia over time.

### Main Contributions
- **Data Organization**: Provided a summary of experimental data on the number of microglia in the penumbra in MCAO-induced mouse models and integrated data from multiple similar studies into a time series dataset suitable for parameter estimation.
- **Mathematical Model**: Proposed a new data-driven mathematical model to describe the dynamic behavior of microglia in the penumbra during ischemic stroke induced by MCAO.
- **Parameter Estimation and Prediction**: Conducted data-driven parameter estimation and forward prediction by combining sensitivity analysis, parameter estimation, and uncertainty quantification methods.

### Research Background
- **Ischemic Stroke**: Stroke is one of the leading causes of death and disability in the United States, with 87% of strokes caused by ischemic stroke due to cerebral vascular occlusion.
- **Neuroinflammation**: Neuroinflammation occurs immediately after ischemic stroke, with microglia being activated and recruited to the penumbra. Microglia can be activated into two phenotypes: M1 (which exacerbates the inflammatory response, leading to further brain damage) and M2 (which aids in long-term recovery).
- **Dynamic Interactions**: M1 and M2 microglia interact during the neuroinflammatory process and can switch phenotypes.

### Methods
- **Data Collection**: Collected data from multiple experimental studies to form a dataset of microglia numbers at multiple time points.
- **Model Construction**: Proposed a set of nonlinear ordinary differential equations (ODEs) to describe the changes in the number of M1 and M2 microglia.
- **Parameter Estimation**: Used global sensitivity analysis and MCMC methods to estimate model parameters and analyze their impact on cell numbers.
- **Prediction and Validation**: Used the estimated parameters for forward prediction and evaluated the model's predictive ability.

### Results
- **Parameter Sensitivity**: Identified key parameters affecting the number of M1 and M2 microglia.
- **Model Prediction**: Simulation results showed that M2 microglia dominate in the early stages of inflammation, followed by a gradual increase and dominance of M1 microglia.
- **Uncertainty Quantification**: Determined the uncertainty range of parameters through MCMC methods and performed forward predictions to assess the model's predictive accuracy.

### Future Work
- **Further Validation**: Validate the model's predictive ability with more experimental data.
- **Model Expansion**: Consider the impact of other factors (such as cytokines) to further refine the model.

This paper establishes a data-driven mathematical model to describe the dynamics of microglia after ischemic stroke and provides new insights into the neuroinflammatory process through parameter estimation and prediction.