Agent-based modeling of autophagy reveals emergent regulatory behavior of spatio-temporal autophagy dynamics

Autophagy is a vesicle-mediated pathway for lysosomal degradation, essential under basal and stressed conditions. Various cellular components, including specific proteins, protein aggregates, organelles and intracellular pathogens, are targets for autophagic degradation. Thereby, autophagy controls numerous vital physiological and pathophysiological functions, including cell signaling, differentiation, turnover of cellular components and pathogen defense. Moreover, autophagy enables the cell to recycle cellular components to metabolic substrates, thereby permitting prolonged survival under low nutrient conditions. Due to the multi-faceted roles for autophagy in maintaining cellular and organismal homeostasis and responding to diverse stresses, malfunction of autophagy contributes to both chronic and acute pathologies.

Agent-based modeling represents a novel approach to investigate autophagy dynamics, function and dysfunction with high biological realism. Our model accurately recapitulates short-term behavior and cell-to-cell variability under basal and activated conditions of autophagy. Further, this approach also allows investigation of long-term behaviors emerging from biologically-relevant alterations to vesicle trafficking and metabolic state.

We applied a systems biology approach to improve the understanding of this complex cellular process of autophagy. All autophagy pathway vesicle activities, i.e. creation, movement, fusion and degradation, are highly dynamic, temporally and spatially, and under various forms of regulation. We therefore developed an agent-based model (ABM) to represent individual components of the autophagy pathway, subcellular vesicle dynamics and metabolic feedback with the cellular environment, thereby providing a framework to investigate spatio-temporal aspects of autophagy regulation and dynamic behavior. The rules defining our ABM were derived from literature and from high-resolution images of autophagy markers under basal and activated conditions. Key model parameters were fit with an iterative method using a genetic algorithm and a predefined fitness function. From this approach, we found that accurate prediction of spatio-temporal behavior required increasing model complexity by implementing functional integration of autophagy with the cellular nutrient state. The resulting model is able to reproduce short-term autophagic flux measurements (up to 3 hours) under basal and activated autophagy conditions, and to measure the degree of cell-to-cell variability. Moreover, we experimentally confirmed two model predictions, namely (i) peri-nuclear concentration of autophagosomes and (ii) inhibitory lysosomal feedback on mTOR signaling.