Muscle Regeneration Agent-Based Model Predicts Enhanced Recovery Outcomes with Altered Cytokine Dynamics
Megan Haase
Department of Biomedical Engineering, University of Virginia, Charlottesville, VA, USA�Email: mh2uk@virginia.edu
Twitter: @biomeganics
IMAG/MSM Working Group
June 8th, 2023
2
Muscle injuries account for more than 30% of all injuries and are one of the most common complaints to orthopedics
Exercise-induced muscle damage
Skeletal muscle has remarkable regenerative capabilities
Background
Muscle regeneration is essential for everyday life
Valle et al. (2011) Br J Sports Med
Muscle regeneration is a highly coordinated cellular process dependent on various signaling factors
3
Healthy muscle fiber
Extracellular matrix (ECM)
Inflammatory phase
Muscle damage
Neutrophil
Damaged fiber
Remodeling phase
Differentiated SSC
Anti-inflammatory macrophage
Regeneration and muscle growth phase
Fusing SSC
Pro-inflammatory macrophage
Activated satellite stem cell (SSC)
Fibroblast
Background
Skeletal muscle regeneration requires numerous cell types and microenvironmental signals
Yang & Hu (2018) J Orthop Translat
4
new experiments and treatment
insight into cell and cytokine cross-talk
predict time-varying spatial information
test altered cytokine dynamics
new
experiments
targets
for cytokine
therapies
decades of excellent muscle regeneration research
computational models
+ physiologically based rules
Computational modeling allows us to study complex dynamics and gain new insight into muscle regeneration
Background
5
ABMs allow simulation of intricate interactions between cells
cell division
fusion
activation
differentiation
Example model agents
Satellite stem cell
Neutrophil
phagocytosis of damage
apoptosis
Macrophage
phagocytosis of damage and apoptotic cells
HGF
MCP
VEGF
TNF
TGF
MMP
IL-10
HGF
TNF
TGF
Division signal: TNF-α + VEGF - TGF-β
Differentiation signal: 3*IL10 - HGF - TNF - TGF-β
Background
6
Limited diffusion and simplified cell dynamics
Previous regeneration models weren’t built for studying cytokine dynamics
Martin et al. (2015) J Appl Physiol
Focuses on spatial cell dynamics and diffusion
No spatial diffusion or spatial dynamics of inflammatory cells
Virgilio et al. (2018) J Appl Physiol
Westman et al. (2021) PLOS Comp Bio
Background
7
Goal: Build ABM to predict regeneration outcomes with altered cytokine dynamics
8
Complex cell behaviors and cytokines were modeled spatially
quiescent
SSC
activated
SSC
satellite cell lineage
myoblast
myocyte
myotube
fibroblast lineage
fibroblast
myofibroblast
M2
Inflammatory cells
monocyte
neutrophil
M1
ECM
fiber
capillary
Microstructure
Example behaviors:
Cytokines
MMP
HGF
VEGF
MCP
TNF
TGF
IL-10
Model Development
9
Spatial characteristics are defined from histology
User defined damage
Damaged fibers secrete HGF and TGF
Assign microstructure
elements
Randomize placement
of baseline cells
Import muscle
histology
lymphatic
vessel
muscle
fiber
ECM
capillary
necrosis
macrophage
fibroblast
SSC
Model Development
10
Microvascular growth and remodeling plays a key role in muscle regeneration
Umek et al. (2019) Histochem and Cell Bio
Capillary
New capillary
formation
Model Development
11
Dynamic ECM properties were integrated to more accurately represent cytokine diffusion with altered collagen density
ECM
Altered
diffusivity
Collagen changes
Filion & Popel (2005) Am J Physiol - Hear Circ Physiol
rs = cytokine radius
rf = fiber radius
ϕ = fiber volume fraction
D∞ = cytokine diffusivity in free solution
D = diffusivity of cytokine in matrix
Model Development
12
Cellular behaviors are governed by literature derived rules
Model Development
13
Cellular behaviors are governed by literature derived rules
Model Development
14
HGF
MMP
TGF-β
TNF-α
VEGF
Muscle Cross-Section
Cellular Interactions
IL-10
MCP
15
CaliPro can identify a robust parameter space for stochastic biological models
Joslyn et al. (2020) Cellular and Molecular Bioengineering
Model Calibration
16
Rules that have not been experimentally quantified create unknown parameters
Model Calibration
See supplemental slides for 38 rule references
17
Parameter density estimation and partial rank correlation coefficient to calibrate the model
Model Calibration
18
CSA recovery, SSC and fibroblast count outputs were fit to literature data
Ochoa et al. (2007) Am J Physiol - Regul Integr Comp Physiol
Murphy et al. (2011) Development
Model Calibration
19
Model replicates inflammatory cell count and capillary count per fiber area
Hardy et al. (2016) PLoS One
Wang et al. (2018) J Neuroimmunol
Nguyen et al. (2011) Sci World J
Ochoa et al. (2007) Am J Physiol - Regul Integr Comp Physiol
Model Validation
20
1Teixeira et al. (2003) Muscle & Nerve
2Raimondo & Mooney (2018) PNAS
3Deng et al. (2012) J Immunol
4Hardy et al. (2019) Skelet Muscle
5Arsic et al. (2004) Mol Ther
6Lu et al. (2011) FASEB J.
7Chen et al. (2005) Am J Physiol Cell Physiol.
8Liu et al. (2016) Cell Bio Int.
Model Validation
Model inputs altered to simulate various experimental conditions
21
Model Predictions
ABM provides deeper understanding of response to altered angiogenesis during muscle regeneration
22
Model Predictions
Cytokine knockout perturbations revealed crosstalk and temporal interplay between cytokines
23
| CSA
| SSC
| Fibroblasts | Non-perfused capillaries | Myoblasts | Myocytes | Neutrophils | M1
| M2 |
Day | 16.7 | 6.3 | 10.5 | 8.4 | 6.3 | 8.4 | 8.4 | 4.2 | 6.3 |
HGF decay | - | - | - | + | - | - | + |
| + |
TGF-β decay | + | + | + | - | + |
|
|
| - |
MMP decay | + | + | + | - | + | + |
|
|
|
TNF-α decay |
|
| + |
|
|
|
|
| - |
VEGF decay |
|
|
| + |
|
|
|
|
|
MCP decay |
|
|
|
|
|
|
| + | + |
MCP diffusion |
| + |
| - |
|
|
| + |
|
LHS-PRCC relationships between cytokine parameters and key regeneration metrics
Model Predictions
24
Combined alterations of cytokine dynamics enhance muscle regeneration
Model Predictions
25
Agent-based modeling provides new hypotheses for in vivo experiments
Conclusion
Acknowledgements
26
M3 Lab Members:
Silvia Blemker, PhD
Xiao Hu, PhD
Emily McCain, PhD
Matthew DiSalvo
Ridhi Sahani
Allison McCrady
Mario Garcia
Jacob Dunn
Undergraduates:
Autumn Routt
Anne Felipe
Brendan Shea
Keerthana Vijayaragavan
Josiah Calhoun
CC3D Collaborators:
James Glazier
T.J. Sego, PhD
Tien Comlekoglu
Alexa Petrucciani
Funding:
Figures created with BioRender.com
27
Thank You!
28
Thank You!