Computational Analysis of Fibroblast Subpopulation Dynamics as a Driver of Fibrotic Foci Formation

Fibroblasts maintain the extracellular matrix (ECM) to support tissue homeostasis and wound healing. In fibrotic diseases, fibroblasts are a primary driver of disease progression through excess collagen secretion and enhanced contractility. Replacing native tissue with a collagen rich fibrotic scar leads to a decline in tissue function. Recent research into idiopathic pulmonary fibrosis (IPF) has identified fibroblast subpopulations that may be primed for the hyper-activation that leads to increased progression of fibrotic disease. Understanding the contribution of these subpopulations to disease progression requires integrating experimental and computational techniques to understand their dynamic contributions to tissue phenotype. Herein, we introduce a framework for modeling subpopulations using a multiscale mechanistic computational model to understand differences within subpopulations, at the intracellular level and how these differences contribute to cell-and tissue-level pathology. We build and validate this framework using two well-defined subpopulations of fibroblasts in IPF. The subpopulations are defined by the presence or absence of Thy-1, a cell-surface protein that regulates fibroblast mechanosensing. We first developed a logic-based network model of a fibroblast. We then applied this model to identify sub-networks that regulate myofibroblast marker expression in the two subpopulations. Coupling this with an agent-based model (ABM) of the lung microenvironment, we observed how different rules regulating cell fate decisions in each subpopulation affected collagen content. Computational image outputs were analyzed with the open-source biological image analysis software QuPath to quantify how changes in subpopulation dynamics change model-predicted foci characteristics such as size and collagen density. We find that the ability for Thy-1 ⁺ fibroblasts to transition to Thy-1 ^- fibroblasts significantly increases total collagen content, as well as influences fibrotic foci characteristics. Overall, we present a combined experimental and computational framework for studying how dynamic changes in fibroblast subpopulations lead to tissue-level disease phenotypes.

In wound healing fibroblasts are responsible for rebuilding the scaffolding, called the extracellular matrix (ECM), of the damaged tissue to aid in regeneration. In fibrosis, the normal wound healing processes are hijacked leading to overproduction of ECM proteins, such as collagen, by fibroblasts leading to fibrosis. In fibrotic diseases with no known cause, such as idiopathic pulmonary fibrosis (IPF), subpopulations of fibroblasts have been identified as possible drivers of disease. Herein, we use a multiscale computational model that represents intracellular, cellular, and tissue level signaling to study how the presence or absence of a single protein on a fibroblast's surface, Thy-1, affects the formation of fibrotic scar in IPF. Thy-1 regulates how a fibroblast senses the stiffness of the microenvironment. The absence of Thy-1 impacted intracellular signaling and when combined across many cells led to fibrosis at the tissue level. Additionally, if normal Thy-1 ⁺ fibroblasts could dynamically lose Thy-1 expression the amount of collagen they produced correlated with that in end stage IPF lungs. This model presents a framework for studying different fibroblast subpopulations using multiscale computational modeling to understand how changes in a single protein in one cell can, over an entire population, affect the dynamics of disease progression.