1
|
Tangella N, Cess CG, Ildefonso GV, Finley SD. Integrating mechanism-based T cell phenotypes into a model of tumor-immune cell interactions. APL Bioeng 2024; 8:036111. [PMID: 39175956 PMCID: PMC11341129 DOI: 10.1063/5.0205996] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2024] [Revised: 08/21/2024] [Accepted: 07/25/2024] [Indexed: 08/24/2024] Open
Abstract
Interactions between cancer cells and immune cells in the tumor microenvironment influence tumor growth and can contribute to the response to cancer immunotherapies. It is difficult to gain mechanistic insights into the effects of cell-cell interactions in tumors using a purely experimental approach. However, computational modeling enables quantitative investigation of the tumor microenvironment, and agent-based modeling, in particular, provides relevant biological insights into the spatial and temporal evolution of tumors. Here, we develop a novel agent-based model (ABM) to predict the consequences of intercellular interactions. Furthermore, we leverage our prior work that predicts the transitions of CD8+ T cells from a naïve state to a terminally differentiated state using Boolean modeling. Given the details incorporated to predict T cell state, we apply the integrated Boolean-ABM framework to study how the properties of CD8+ T cells influence the composition and spatial organization of tumors and the efficacy of an immune checkpoint blockade. Overall, we present a mechanistic understanding of tumor evolution that can be leveraged to study targeted immunotherapeutic strategies.
Collapse
Affiliation(s)
- Neel Tangella
- Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, California 90089, USA
| | - Colin G. Cess
- Alfred E. Mann Department of Biomedical Engineering, University of Southern California, Los Angeles, California 90089, USA
| | - Geena V. Ildefonso
- Alfred E. Mann Department of Biomedical Engineering, University of Southern California, Los Angeles, California 90089, USA
| | | |
Collapse
|
2
|
Flanary SM, Peak KE, Barocas VH. A Graphical Approach to Visualize and Interpret Biochemically Coupled Biomechanical Models. J Biomech Eng 2024; 146:054504. [PMID: 38421368 PMCID: PMC11005857 DOI: 10.1115/1.4064970] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2023] [Revised: 02/21/2024] [Accepted: 02/22/2024] [Indexed: 03/02/2024]
Abstract
The last decade has seen the emergence of progressively more complex mechanobiological models, often coupling biochemical and biomechanical components. The complexity of these models makes interpretation difficult, and although computational tools can solve model equations, there is considerable potential value in a simple method to explore the interplay between different model components. Pump and system performance curves, long utilized in centrifugal pump selection and design, inspire the development of a graphical technique to depict visually the performance of biochemically-coupled mechanical models. Our approach is based on a biochemical performance curve (analogous to the classical pump curve) and a biomechanical performance curve (analogous to the system curve). Upon construction of the two curves, their intersection, or lack thereof, describes the coupled model's equilibrium state(s). One can also observe graphically how an applied perturbation shifts one or both curves, and thus how the other component will respond, without rerunning the full model. While the upfront cost of generating the performance curve graphic varies with the efficiency of the model components, the easily interpretable visual depiction of what would otherwise be nonintuitive model behavior is valuable. Herein, we outline how performance curves can be constructed and interpreted for biochemically-coupled biomechanical models and apply the technique to two independent models in the cardiovascular space. The performance curve approach can illustrate and help identify weaknesses in model construction, inform user-applied perturbations and fitting procedures to generate intended behaviors, and improve the efficiency of the model generation and application process.
Collapse
Affiliation(s)
- Shannon M. Flanary
- Department of Chemical Engineering & Materials Science, University of Minnesota, Nils Hasselmo Hall, Room 7-115, 312 Church St SE, Minneapolis, MN 55455
| | - Kara E. Peak
- Department of Biomedical Engineering, University of Minnesota, Nils Hasselmo Hall, Room 7-115, 312 Church St SE, Minneapolis, MN 55455
- University of Minnesota
| | - Victor H. Barocas
- Department of Biomedical Engineering, University of Minnesota, Nils Hasselmo Hall, Room 7-115, 312 Church St SE, Minneapolis, MN 55455
| |
Collapse
|
3
|
Chen P, Ye C, Huang Y, Xu B, Wu T, Dong Y, Jin Y, Zhao L, Hu C, Mao J, Wu R. Glutaminolysis regulates endometrial fibrosis in intrauterine adhesion via modulating mitochondrial function. Biol Res 2024; 57:13. [PMID: 38561846 PMCID: PMC10983700 DOI: 10.1186/s40659-024-00492-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2023] [Accepted: 03/22/2024] [Indexed: 04/04/2024] Open
Abstract
BACKGROUND Endometrial fibrosis, a significant characteristic of intrauterine adhesion (IUA), is caused by the excessive differentiation and activation of endometrial stromal cells (ESCs). Glutaminolysis is the metabolic process of glutamine (Gln), which has been implicated in multiple types of organ fibrosis. So far, little is known about whether glutaminolysis plays a role in endometrial fibrosis. METHODS The activation model of ESCs was constructed by TGF-β1, followed by RNA-sequencing analysis. Changes in glutaminase1 (GLS1) expression at RNA and protein levels in activated ESCs were verified experimentally. Human IUA samples were collected to verify GLS1 expression in endometrial fibrosis. GLS1 inhibitor and glutamine deprivation were applied to ESCs models to investigate the biological functions and mechanisms of glutaminolysis in ESCs activation. The IUA mice model was established to explore the effect of glutaminolysis inhibition on endometrial fibrosis. RESULTS We found that GLS1 expression was significantly increased in activated ESCs models and fibrotic endometrium. Glutaminolysis inhibition by GLS1 inhibitor bis-2-(5-phenylacetamido-1,2,4-thiadiazol-2-yl) ethyl sulfide (BPTES or glutamine deprivation treatment suppressed the expression of two fibrotic markers, α-SMA and collagen I, as well as the mitochondrial function and mTORC1 signaling in ESCs. Furthermore, inhibition of the mTORC1 signaling pathway by rapamycin suppressed ESCs activation. In IUA mice models, BPTES treatment significantly ameliorated endometrial fibrosis and improved pregnancy outcomes. CONCLUSION Glutaminolysis and glutaminolysis-associated mTOR signaling play a role in the activation of ESCs and the pathogenesis of endometrial fibrosis through regulating mitochondrial function. Glutaminolysis inhibition suppresses the activation of ESCs, which might be a novel therapeutic strategy for IUA.
Collapse
Affiliation(s)
- Pei Chen
- Department of Obstetrics and Gynecology, Women's Hospital, Zhejiang University School of Medicine, Hangzhou, China
- Key Laboratory of Women's Reproductive Health of Zhejiang Province, Hangzhou, China
| | - Chaoshuang Ye
- Department of Obstetrics and Gynecology, Women's Hospital, Zhejiang University School of Medicine, Hangzhou, China
- Key Laboratory of Women's Reproductive Health of Zhejiang Province, Hangzhou, China
| | - Yunke Huang
- Department of Obstetrics and Gynecology, Women's Hospital, Zhejiang University School of Medicine, Hangzhou, China
- Key Laboratory of Women's Reproductive Health of Zhejiang Province, Hangzhou, China
| | - Bingning Xu
- Department of Obstetrics and Gynecology, Women's Hospital, Zhejiang University School of Medicine, Hangzhou, China
- Key Laboratory of Women's Reproductive Health of Zhejiang Province, Hangzhou, China
| | - Tianyu Wu
- Department of Obstetrics and Gynecology, Women's Hospital, Zhejiang University School of Medicine, Hangzhou, China
- Key Laboratory of Women's Reproductive Health of Zhejiang Province, Hangzhou, China
| | - Yuanhang Dong
- Department of Obstetrics and Gynecology, Women's Hospital, Zhejiang University School of Medicine, Hangzhou, China
- Key Laboratory of Women's Reproductive Health of Zhejiang Province, Hangzhou, China
| | - Yang Jin
- Department of Obstetrics and Gynecology, Women's Hospital, Zhejiang University School of Medicine, Hangzhou, China
- Key Laboratory of Women's Reproductive Health of Zhejiang Province, Hangzhou, China
| | - Li Zhao
- Department of Obstetrics and Gynecology, Women's Hospital, Zhejiang University School of Medicine, Hangzhou, China
- Key Laboratory of Women's Reproductive Health of Zhejiang Province, Hangzhou, China
| | - Changchang Hu
- Department of Obstetrics and Gynecology, Women's Hospital, Zhejiang University School of Medicine, Hangzhou, China
- Key Laboratory of Women's Reproductive Health of Zhejiang Province, Hangzhou, China
| | - Jingxia Mao
- Department of Obstetrics and Gynecology, Women's Hospital, Zhejiang University School of Medicine, Hangzhou, China
- Key Laboratory of Women's Reproductive Health of Zhejiang Province, Hangzhou, China
| | - Ruijin Wu
- Department of Obstetrics and Gynecology, Women's Hospital, Zhejiang University School of Medicine, Hangzhou, China.
- Key Laboratory of Women's Reproductive Health of Zhejiang Province, Hangzhou, China.
| |
Collapse
|
4
|
Leonard-Duke J, Agro SMJ, Csordas DJ, Bruce AC, Eggertsen TG, Tavakol TN, Barker TH, Bonham CA, Saucerman JJ, Taite LJ, Peirce SM. Multiscale computational model predicts how environmental changes and drug treatments affect microvascular remodeling in fibrotic disease. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.03.15.585249. [PMID: 38559112 PMCID: PMC10979947 DOI: 10.1101/2024.03.15.585249] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
Investigating the molecular, cellular, and tissue-level changes caused by disease, and the effects of pharmacological treatments across these biological scales, necessitates the use of multiscale computational modeling in combination with experimentation. Many diseases dynamically alter the tissue microenvironment in ways that trigger microvascular network remodeling, which leads to the expansion or regression of microvessel networks. When microvessels undergo remodeling in idiopathic pulmonary fibrosis (IPF), functional gas exchange is impaired due to loss of alveolar structures and lung function declines. Here, we integrated a multiscale computational model with independent experiments to investigate how combinations of biomechanical and biochemical cues in IPF alter cell fate decisions leading to microvascular remodeling. Our computational model predicted that extracellular matrix (ECM) stiffening reduced microvessel area, which was accompanied by physical uncoupling of endothelial cell (ECs) and pericytes, the cells that comprise microvessels. Nintedanib, an FDA-approved drug for treating IPF, was predicted to further potentiate microvessel regression by decreasing the percentage of quiescent pericytes while increasing the percentage of pericytes undergoing pericyte-myofibroblast transition (PMT) in high ECM stiffnesses. Importantly, the model suggested that YAP/TAZ inhibition may overcome the deleterious effects of nintedanib by promoting EC-pericyte coupling and maintaining microvessel homeostasis. Overall, our combination of computational and experimental modeling can explain how cell decisions affect tissue changes during disease and in response to treatments.
Collapse
Affiliation(s)
- Julie Leonard-Duke
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, USA
- Robert M. Berne Cardiovascular Research Center, University of Virginia, Charlottesville, Virginia, USA
| | - Samuel M. J. Agro
- Department of Chemical Engineering, University of Virginia, Charlottesville, Virginia, USA
| | - David J. Csordas
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, USA
- Robert M. Berne Cardiovascular Research Center, University of Virginia, Charlottesville, Virginia, USA
| | - Anthony C. Bruce
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, USA
| | - Taylor G. Eggertsen
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, USA
- Robert M. Berne Cardiovascular Research Center, University of Virginia, Charlottesville, Virginia, USA
| | - Tara N. Tavakol
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, USA
| | - Thomas H. Barker
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, USA
- Robert M. Berne Cardiovascular Research Center, University of Virginia, Charlottesville, Virginia, USA
| | - Catherine A. Bonham
- Department of Pulmonary and Critical Care Medicine, University of Virginia, Charlottesville, Virginia, USA
| | - Jeffery J. Saucerman
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, USA
- Robert M. Berne Cardiovascular Research Center, University of Virginia, Charlottesville, Virginia, USA
| | - Lakeshia J. Taite
- Department of Chemical Engineering, University of Virginia, Charlottesville, Virginia, USA
| | - Shayn M. Peirce
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, USA
- Robert M. Berne Cardiovascular Research Center, University of Virginia, Charlottesville, Virginia, USA
| |
Collapse
|
5
|
Bergman DR, Norton KA, Jain HV, Jackson T. Connecting Agent-Based Models with High-Dimensional Parameter Spaces to Multidimensional Data Using SMoRe ParS: A Surrogate Modeling Approach. Bull Math Biol 2023; 86:11. [PMID: 38159216 PMCID: PMC10757706 DOI: 10.1007/s11538-023-01240-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2023] [Accepted: 11/22/2023] [Indexed: 01/03/2024]
Abstract
Across a broad range of disciplines, agent-based models (ABMs) are increasingly utilized for replicating, predicting, and understanding complex systems and their emergent behavior. In the biological and biomedical sciences, researchers employ ABMs to elucidate complex cellular and molecular interactions across multiple scales under varying conditions. Data generated at these multiple scales, however, presents a computational challenge for robust analysis with ABMs. Indeed, calibrating ABMs remains an open topic of research due to their own high-dimensional parameter spaces. In response to these challenges, we extend and validate our novel methodology, Surrogate Modeling for Reconstructing Parameter Surfaces (SMoRe ParS), arriving at a computationally efficient framework for connecting high dimensional ABM parameter spaces with multidimensional data. Specifically, we modify SMoRe ParS to initially confine high dimensional ABM parameter spaces using unidimensional data, namely, single time-course information of in vitro cancer cell growth assays. Subsequently, we broaden the scope of our approach to encompass more complex ABMs and constrain parameter spaces using multidimensional data. We explore this extension with in vitro cancer cell inhibition assays involving the chemotherapeutic agent oxaliplatin. For each scenario, we validate and evaluate the effectiveness of our approach by comparing how well ABM simulations match the experimental data when using SMoRe ParS-inferred parameters versus parameters inferred by a commonly used direct method. In so doing, we show that our approach of using an explicitly formulated surrogate model as an interlocutor between the ABM and the experimental data effectively calibrates the ABM parameter space to multidimensional data. Our method thus provides a robust and scalable strategy for leveraging multidimensional data to inform multiscale ABMs and explore the uncertainty in their parameters.
Collapse
Affiliation(s)
- Daniel R Bergman
- Department of Mathematics, University of Michigan, 530 Church Street, Ann Arbor, MI, 48109, USA
| | - Kerri-Ann Norton
- Computational Biology Laboratory, Computer Science Program, Bard College, 30 Campus Road, Annandale-on-Hudson, NY, 12504, USA
| | - Harsh Vardhan Jain
- Department of Mathematics & Statistics, University of Minnesota Duluth, 1117 University Drive, Duluth, MN, 55812, USA
| | - Trachette Jackson
- Department of Mathematics, University of Michigan, 530 Church Street, Ann Arbor, MI, 48109, USA.
| |
Collapse
|
6
|
Islam MA, Getz M, Macklin P, Ford Versypt AN. An agent-based modeling approach for lung fibrosis in response to COVID-19. PLoS Comput Biol 2023; 19:e1011741. [PMID: 38127835 PMCID: PMC10769079 DOI: 10.1371/journal.pcbi.1011741] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2023] [Revised: 01/05/2024] [Accepted: 12/05/2023] [Indexed: 12/23/2023] Open
Abstract
The severity of the COVID-19 pandemic has created an emerging need to investigate the long-term effects of infection on patients. Many individuals are at risk of suffering pulmonary fibrosis due to the pathogenesis of lung injury and impairment in the healing mechanism. Fibroblasts are the central mediators of extracellular matrix (ECM) deposition during tissue regeneration, regulated by anti-inflammatory cytokines including transforming growth factor beta (TGF-β). The TGF-β-dependent accumulation of fibroblasts at the damaged site and excess fibrillar collagen deposition lead to fibrosis. We developed an open-source, multiscale tissue simulator to investigate the role of TGF-β sources in the progression of lung fibrosis after SARS-CoV-2 exposure, intracellular viral replication, infection of epithelial cells, and host immune response. Using the model, we predicted the dynamics of fibroblasts, TGF-β, and collagen deposition for 15 days post-infection in virtual lung tissue. Our results showed variation in collagen area fractions between 2% and 40% depending on the spatial behavior of the sources (stationary or mobile), the rate of activation of TGF-β, and the duration of TGF-β sources. We identified M2 macrophages as primary contributors to higher collagen area fraction. Our simulation results also predicted fibrotic outcomes even with lower collagen area fraction when spatially-localized latent TGF-β sources were active for longer times. We validated our model by comparing simulated dynamics for TGF-β, collagen area fraction, and macrophage cell population with independent experimental data from mouse models. Our results showed that partial removal of TGF-β sources changed the fibrotic patterns; in the presence of persistent TGF-β sources, partial removal of TGF-β from the ECM significantly increased collagen area fraction due to maintenance of chemotactic gradients driving fibroblast movement. The computational findings are consistent with independent experimental and clinical observations of collagen area fractions and cell population dynamics not used in developing the model. These critical insights into the activity of TGF-β sources may find applications in the current clinical trials targeting TGF-β for the resolution of lung fibrosis.
Collapse
Affiliation(s)
- Mohammad Aminul Islam
- Department of Chemical and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York, United States of America
| | - Michael Getz
- Department of Intelligent Systems Engineering, Indiana University, Bloomington, Indiana, United States of America
| | - Paul Macklin
- Department of Intelligent Systems Engineering, Indiana University, Bloomington, Indiana, United States of America
| | - Ashlee N. Ford Versypt
- Department of Chemical and Biological Engineering, University at Buffalo, The State University of New York, Buffalo, New York, United States of America
- Department of Biomedical Engineering, University at Buffalo, The State University of New York, Buffalo, New York, United States of America
- Institute for Artificial Intelligence and Data Science, University at Buffalo, The State University of New York, Buffalo, New York, United States of America
| |
Collapse
|
7
|
Islam MA, Getz M, Macklin P, Versypt ANF. An agent-based modeling approach for lung fibrosis in response to COVID-19. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2022.10.03.510677. [PMID: 36238719 PMCID: PMC9558432 DOI: 10.1101/2022.10.03.510677] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
The severity of the COVID-19 pandemic has created an emerging need to investigate the long-term effects of infection on patients. Many individuals are at risk of suffering pulmonary fibrosis due to the pathogenesis of lung injury and impairment in the healing mechanism. Fibroblasts are the central mediators of extracellular matrix (ECM) deposition during tissue regeneration, regulated by anti-inflammatory cytokines including transforming growth factor beta (TGF-β). The TGF-β-dependent accumulation of fibroblasts at the damaged site and excess fibrillar collagen deposition lead to fibrosis. We developed an open-source, multiscale tissue simulator to investigate the role of TGF-β sources in the progression of lung fibrosis after SARS-CoV-2 exposure, intracellular viral replication, infection of epithelial cells, and host immune response. Using the model, we predicted the dynamics of fibroblasts, TGF-β, and collagen deposition for 15 days post-infection in virtual lung tissue. Our results showed variation in collagen area fractions between 2% and 40% depending on the spatial behavior of the sources (stationary or mobile), the rate of activation of TGF-β, and the duration of TGF-β sources. We identified M2 macrophages as primary contributors to higher collagen area fraction. Our simulation results also predicted fibrotic outcomes even with lower collagen area fraction when spatially-localized latent TGF-β sources were active for longer times. We validated our model by comparing simulated dynamics for TGF-β, collagen area fraction, and macrophage cell population with independent experimental data from mouse models. Our results showed that partial removal of TGF-β sources changed the fibrotic patterns; in the presence of persistent TGF-β sources, partial removal of TGF-β from the ECM significantly increased collagen area fraction due to maintenance of chemotactic gradients driving fibroblast movement. The computational findings are consistent with independent experimental and clinical observations of collagen area fractions and cell population dynamics not used in developing the model. These critical insights into the activity of TGF-β sources may find applications in the current clinical trials targeting TGF-β for the resolution of lung fibrosis. Author summary COVID-19 survivors are at risk of lung fibrosis as a long-term effect. Lung fibrosis is the excess deposition of tissue materials in the lung that hinder gas exchange and can collapse the whole organ. We identified TGF-β as a critical regulator of fibrosis. We built a model to investigate the mechanisms of TGF-β sources in the process of fibrosis. Our results showed spatial behavior of sources (stationary or mobile) and their activity (activation rate of TGF-β, longer activation of sources) could lead to lung fibrosis. Current clinical trials for fibrosis that target TGF-β need to consider TGF-β sources' spatial properties and activity to develop better treatment strategies.
Collapse
|
8
|
Chowkwale M, Lindsey ML, Saucerman JJ. Intercellular model predicts mechanisms of inflammation-fibrosis coupling after myocardial infarction. J Physiol 2023; 601:2635-2654. [PMID: 35862254 PMCID: PMC9859968 DOI: 10.1113/jp283346] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2022] [Accepted: 07/18/2022] [Indexed: 01/25/2023] Open
Abstract
After myocardial infarction (MI), cardiac cells work together to regulate wound healing of the infarct. The pathological response to MI yields cardiac remodelling comprising inflammatory and fibrosis phases, and the interplay of cellular dynamics that underlies these phases has not been elucidated. This study developed a computational model to identify cytokine and cellular dynamics post-MI to predict mechanisms driving post-MI inflammation, resolution of inflammation, and scar formation. Additionally, this study evaluated the interdependence between inflammation and fibrosis. Our model bypassed limitations of in vivo approaches in achieving cellular specificity and performing specific perturbations such as global knockouts of chemical factors. The model predicted that inflammation is a graded response to initial infarct size that is amplified by a positive feedback loop between neutrophils and interleukin 1β (IL-1β). Resolution of inflammation was driven by degradation of IL-1β, matrix metalloproteinase 9, and transforming growth factor β (TGF-β), as well as apoptosis of neutrophils. Inflammation regulated TGFβ secretion directly through immune cell recruitment and indirectly through upregulation of macrophage phagocytosis. Lastly, we found that mature collagen deposition was an ultrasensitive switch in response to inflammation, which was amplified primarily by cardiac fibroblast proliferation. These findings describe the relationship between inflammation and fibrosis and highlight how the two responses work together post-MI. This model revealed that post-MI inflammation and fibrosis are dynamically coupled, which provides rationale for designing novel anti-inflammatory, pro-resolving or anti-fibrotic therapies that may improve the response to MI. KEY POINTS: Inflammation and matrix remodelling are two processes involved in wound healing after a heart attack. Cardiac cells work together to facilitate these processes; this is done by secreting cytokines that then regulate the cells themselves or other cells surrounding them. This study developed a computational model of the dynamics of cardiac cells and cytokines to predict mechanisms through which inflammation and matrix remodelling is regulated. We show the roles of various cytokines and signalling motifs in driving inflammation, resolution of inflammation and fibrosis. The novel concept of inflammation-fibrosis coupling, based on the model prediction that inflammation and fibrosis are dynamically coupled, provides rationale for future studies and for designing therapeutics to improve the response after a heart attack.
Collapse
Affiliation(s)
- Mukti Chowkwale
- Department of Biomedical Engineering, University of Virginia, Charlottesville, VA
| | - Merry L. Lindsey
- School of Graduate Studies and Research, Meharry Medical College, Nashville, TN
- Research Service, Nashville VA Medical Center, Nashville, TN
| | - Jeffrey J. Saucerman
- Department of Biomedical Engineering, University of Virginia, Charlottesville, VA
- Robert M. Berne Cardiovascular Research Center, University of Virginia, Charlottesville, VA
| |
Collapse
|
9
|
Cogno N, Bauer R, Durante M. An Agent-Based Model of Radiation-Induced Lung Fibrosis. Int J Mol Sci 2022; 23:ijms232213920. [PMID: 36430398 PMCID: PMC9693125 DOI: 10.3390/ijms232213920] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2022] [Revised: 11/03/2022] [Accepted: 11/05/2022] [Indexed: 11/16/2022] Open
Abstract
Early- and late-phase radiation-induced lung injuries, namely pneumonitis and lung fibrosis (RILF), severely constrain the maximum dose and irradiated volume in thoracic radiotherapy. As the most radiosensitive targets, epithelial cells respond to radiation either by undergoing apoptosis or switching to a senescent phenotype that triggers the immune system and damages surrounding healthy cells. Unresolved inflammation stimulates mesenchymal cells' proliferation and extracellular matrix (ECM) secretion, which irreversibly stiffens the alveolar walls and leads to respiratory failure. Although a thorough understanding is lacking, RILF and idiopathic pulmonary fibrosis share multiple pathways and would mutually benefit from further insights into disease progression. Furthermore, current normal tissue complication probability (NTCP) models rely on clinical experience to set tolerance doses for organs at risk and leave aside mechanistic interpretations of the undergoing processes. To these aims, we implemented a 3D agent-based model (ABM) of an alveolar duct that simulates cell dynamics and substance diffusion following radiation injury. Emphasis was placed on cell repopulation, senescent clearance, and intra/inter-alveolar bystander senescence while tracking ECM deposition. Our ABM successfully replicates early and late fibrotic response patterns reported in the literature along with the ECM sigmoidal dose-response curve. Moreover, surrogate measures of RILF severity via a custom indicator show qualitative agreement with published fibrosis indices. Finally, our ABM provides a fully mechanistic alveolar survival curve highlighting the need to include bystander damage in lung NTCP models.
Collapse
Affiliation(s)
- Nicolò Cogno
- Biophysics Department, GSI Helmholtzzentrum für Schwerionenforschung GmbH, 64291 Darmstadt, Germany
- Institute for Condensed Matter Physics, Technische Universität Darmstadt, 64289 Darmstadt, Germany
| | - Roman Bauer
- Department of Computer Science, University of Surrey, Guildford GU2 7XH, UK
| | - Marco Durante
- Biophysics Department, GSI Helmholtzzentrum für Schwerionenforschung GmbH, 64291 Darmstadt, Germany
- Institute for Condensed Matter Physics, Technische Universität Darmstadt, 64289 Darmstadt, Germany
- Correspondence: or
| |
Collapse
|
10
|
Sivakumar N, Warner HV, Peirce SM, Lazzara MJ. A computational modeling approach for predicting multicell spheroid patterns based on signaling-induced differential adhesion. PLoS Comput Biol 2022; 18:e1010701. [PMID: 36441822 PMCID: PMC9747056 DOI: 10.1371/journal.pcbi.1010701] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2021] [Revised: 12/13/2022] [Accepted: 11/01/2022] [Indexed: 11/29/2022] Open
Abstract
Physiological and pathological processes including embryogenesis and tumorigenesis rely on the ability of individual cells to work collectively to form multicell patterns. In these heterogeneous multicell systems, cell-cell signaling induces differential adhesion between cells that leads to tissue-level patterning. However, the sensitivity of pattern formation to changes in the strengths of signaling or cell adhesion processes is not well understood. Prior work has explored these issues using synthetically engineered heterogeneous multicell spheroid systems, in which cell subpopulations engage in bidirectional intercellular signaling to regulate the expression of different cadherins. While engineered cell systems provide excellent experimental tools to observe pattern formation in cell populations, computational models of these systems may be leveraged to explore more systematically how specific combinations of signaling and adhesion parameters can drive the emergence of unique patterns. We developed and validated two- and three-dimensional agent-based models (ABMs) of spheroid patterning for previously described cells engineered with a bidirectional signaling circuit that regulates N- and P-cadherin expression. Systematic exploration of model predictions, some of which were experimentally validated, revealed how cell seeding parameters, the order of signaling events, probabilities of induced cadherin expression, and homotypic adhesion strengths affect pattern formation. Unsupervised clustering was also used to map combinations of signaling and adhesion parameters to these unique spheroid patterns predicted by the ABM. Finally, we demonstrated how the model may be deployed to design new synthetic cell signaling circuits based on a desired final multicell pattern.
Collapse
Affiliation(s)
- Nikita Sivakumar
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, United States of America
| | - Helen V. Warner
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, United States of America
| | - Shayn M. Peirce
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, United States of America
| | - Matthew J. Lazzara
- Department of Biomedical Engineering, University of Virginia, Charlottesville, Virginia, United States of America
- Department of Chemical Engineering, University of Virginia, Charlottesville, Virginia, United States of America
- * E-mail:
| |
Collapse
|
11
|
Aminian-Dehkordi J, Valiei A, Mofrad MRK. Emerging computational paradigms to address the complex role of gut microbial metabolism in cardiovascular diseases. Front Cardiovasc Med 2022; 9:987104. [PMID: 36299869 PMCID: PMC9589059 DOI: 10.3389/fcvm.2022.987104] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2022] [Accepted: 09/20/2022] [Indexed: 11/13/2022] Open
Abstract
The human gut microbiota and its associated perturbations are implicated in a variety of cardiovascular diseases (CVDs). There is evidence that the structure and metabolic composition of the gut microbiome and some of its metabolites have mechanistic associations with several CVDs. Nevertheless, there is a need to unravel metabolic behavior and underlying mechanisms of microbiome-host interactions. This need is even more highlighted when considering that microbiome-secreted metabolites contributing to CVDs are the subject of intensive research to develop new prevention and therapeutic techniques. In addition to the application of high-throughput data used in microbiome-related studies, advanced computational tools enable us to integrate omics into different mathematical models, including constraint-based models, dynamic models, agent-based models, and machine learning tools, to build a holistic picture of metabolic pathological mechanisms. In this article, we aim to review and introduce state-of-the-art mathematical models and computational approaches addressing the link between the microbiome and CVDs.
Collapse
Affiliation(s)
| | | | - Mohammad R. K. Mofrad
- Department of Bioengineering and Mechanical Engineering, University of California, Berkeley, Berkeley, CA, United States
| |
Collapse
|
12
|
Guo Y, Mofrad MRK, Tepole AB. On modeling the multiscale mechanobiology of soft tissues: Challenges and progress. BIOPHYSICS REVIEWS 2022; 3:031303. [PMID: 38505274 PMCID: PMC10903412 DOI: 10.1063/5.0085025] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/12/2022] [Accepted: 07/12/2022] [Indexed: 03/21/2024]
Abstract
Tissues grow and remodel in response to mechanical cues, extracellular and intracellular signals experienced through various biological events, from the developing embryo to disease and aging. The macroscale response of soft tissues is typically nonlinear, viscoelastic anisotropic, and often emerges from the hierarchical structure of tissues, primarily their biopolymer fiber networks at the microscale. The adaptation to mechanical cues is likewise a multiscale phenomenon. Cell mechanobiology, the ability of cells to transform mechanical inputs into chemical signaling inside the cell, and subsequent regulation of cellular behavior through intra- and inter-cellular signaling networks, is the key coupling at the microscale between the mechanical cues and the mechanical adaptation seen macroscopically. To fully understand mechanics of tissues in growth and remodeling as observed at the tissue level, multiscale models of tissue mechanobiology are essential. In this review, we summarize the state-of-the art modeling tools of soft tissues at both scales, the tissue level response, and the cell scale mechanobiology models. To help the interested reader become more familiar with these modeling frameworks, we also show representative examples. Our aim here is to bring together scientists from different disciplines and enable the future leap in multiscale modeling of tissue mechanobiology.
Collapse
Affiliation(s)
- Yifan Guo
- School of Mechanical Engineering, Purdue University, West Lafayette, Indiana 47907, USA
| | - Mohammad R. K. Mofrad
- Departments of Bioengineering and Mechanical Engineering, University of California Berkeley, Berkeley, California 94720, USA
| | - Adrian Buganza Tepole
- School of Mechanical Engineering, Purdue University, West Lafayette, Indiana 47907, USA
| |
Collapse
|
13
|
Bermudez-Lekerika P, Crump KB, Tseranidou S, Nüesch A, Kanelis E, Alminnawi A, Baumgartner L, Muñoz-Moya E, Compte R, Gualdi F, Alexopoulos LG, Geris L, Wuertz-Kozak K, Le Maitre CL, Noailly J, Gantenbein B. Immuno-Modulatory Effects of Intervertebral Disc Cells. Front Cell Dev Biol 2022; 10:924692. [PMID: 35846355 PMCID: PMC9277224 DOI: 10.3389/fcell.2022.924692] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2022] [Accepted: 05/20/2022] [Indexed: 11/29/2022] Open
Abstract
Low back pain is a highly prevalent, chronic, and costly medical condition predominantly triggered by intervertebral disc degeneration (IDD). IDD is often caused by structural and biochemical changes in intervertebral discs (IVD) that prompt a pathologic shift from an anabolic to catabolic state, affecting extracellular matrix (ECM) production, enzyme generation, cytokine and chemokine production, neurotrophic and angiogenic factor production. The IVD is an immune-privileged organ. However, during degeneration immune cells and inflammatory factors can infiltrate through defects in the cartilage endplate and annulus fibrosus fissures, further accelerating the catabolic environment. Remarkably, though, catabolic ECM disruption also occurs in the absence of immune cell infiltration, largely due to native disc cell production of catabolic enzymes and cytokines. An unbalanced metabolism could be induced by many different factors, including a harsh microenvironment, biomechanical cues, genetics, and infection. The complex, multifactorial nature of IDD brings the challenge of identifying key factors which initiate the degenerative cascade, eventually leading to back pain. These factors are often investigated through methods including animal models, 3D cell culture, bioreactors, and computational models. However, the crosstalk between the IVD, immune system, and shifted metabolism is frequently misconstrued, often with the assumption that the presence of cytokines and chemokines is synonymous to inflammation or an immune response, which is not true for the intact disc. Therefore, this review will tackle immunomodulatory and IVD cell roles in IDD, clarifying the differences between cellular involvements and implications for therapeutic development and assessing models used to explore inflammatory or catabolic IVD environments.
Collapse
Affiliation(s)
- Paola Bermudez-Lekerika
- Tissue Engineering for Orthopaedics and Mechanobiology, Bone and Joint Program, Department for BioMedical Research (DBMR), Faculty of Medicine, University of Bern, Bern, Switzerland.,Department of Orthopaedic Surgery and Traumatology, Inselspital, Bern University Hospital, Medical Faculty, University of Bern, Bern, Switzerland
| | - Katherine B Crump
- Tissue Engineering for Orthopaedics and Mechanobiology, Bone and Joint Program, Department for BioMedical Research (DBMR), Faculty of Medicine, University of Bern, Bern, Switzerland.,Department of Orthopaedic Surgery and Traumatology, Inselspital, Bern University Hospital, Medical Faculty, University of Bern, Bern, Switzerland
| | | | - Andrea Nüesch
- Biomolecular Sciences Research Centre, Sheffield Hallam University, Sheffield, United Kingdom
| | - Exarchos Kanelis
- ProtATonce Ltd., Athens, Greece.,School of Mechanical Engineering, National Technical University of Athens, Zografou, Greece
| | - Ahmad Alminnawi
- GIGA In Silico Medicine, University of Liège, Liège, Belgium.,Skeletal Biology and Engineering Research Center, KU Leuven, Leuven, Belgium
| | | | | | - Roger Compte
- Twin Research and Genetic Epidemiology, St Thomas' Hospital, King's College London, London, United Kingdom
| | - Francesco Gualdi
- Institut Hospital Del Mar D'Investigacions Mèdiques (IMIM), Barcelona, Spain
| | - Leonidas G Alexopoulos
- ProtATonce Ltd., Athens, Greece.,School of Mechanical Engineering, National Technical University of Athens, Zografou, Greece
| | - Liesbet Geris
- GIGA In Silico Medicine, University of Liège, Liège, Belgium.,Skeletal Biology and Engineering Research Center, KU Leuven, Leuven, Belgium.,Biomechanics Research Unit, KU Leuven, Leuven, Belgium
| | - Karin Wuertz-Kozak
- Department of Biomedical Engineering, Rochester Institute of Technology, Rochester, NY, United States.,Spine Center, Schön Klinik München Harlaching Academic Teaching Hospital and Spine Research Institute of the Paracelsus Private Medical University Salzburg (Austria), Munich, Germany
| | - Christine L Le Maitre
- Biomolecular Sciences Research Centre, Sheffield Hallam University, Sheffield, United Kingdom
| | | | - Benjamin Gantenbein
- Tissue Engineering for Orthopaedics and Mechanobiology, Bone and Joint Program, Department for BioMedical Research (DBMR), Faculty of Medicine, University of Bern, Bern, Switzerland.,Department of Orthopaedic Surgery and Traumatology, Inselspital, Bern University Hospital, Medical Faculty, University of Bern, Bern, Switzerland
| |
Collapse
|
14
|
In silico optimization of heparin microislands in microporous annealed particle (MAP) hydrogel for endothelial cell migration. Acta Biomater 2022; 148:171-180. [PMID: 35660016 DOI: 10.1016/j.actbio.2022.05.049] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2022] [Revised: 05/03/2022] [Accepted: 05/27/2022] [Indexed: 11/24/2022]
Abstract
Biomaterials capable of generating growth factor gradients have shown success in guiding tissue regeneration, as growth factor gradients are a physiologic driver of cell migration. Of particular importance, a focus on promoting endothelial cell migration is vital to angiogenesis and new tissue formation. Microporous Annealed Particle (MAP) scaffolds represent a unique niche in the field of regenerative biomaterials research as an injectable biomaterial with an open porosity that allows cells to freely migrate independent of material degradation. Recently, we have used the MAP platform to heterogeneously include spatially isolated heparin-modified microgels (heparin microislands) which can sequester growth factors and guide cell migration. In in vitro sprouting angiogenesis assays, we observed a parabolic relationship between the percentage of heparin microislands and cell migration, where 10% heparin microislands had more endothelial cell migration compared to 1% and 100%. Due to the low number of heparin microisland ratios tested, we hypothesize the spacing between microgels can be further optimized. Rather than use purely empirical methods, which are both expensive and time intensive, we believe this challenge represents an opportunity to use computational modeling. Here we present the first agent-based model of a MAP scaffold to optimize the ratio of heparin microislands. Specifically, we develop a two-dimensional model in Hybrid Automata Library (HAL) of endothelial cell migration within the unique MAP scaffold geometry. Finally, we present how our model can accurately predict cell migration trends in vitro, and these studies provide insight on how computational modeling can be used to design particle-based biomaterials. STATEMENT OF SIGNIFICANCE: : While the combination of experimental and computational approaches is increasingly being used to gain a better understanding of cellular processes, their combination in biomaterials development has been relatively limited. Heparin microislands are spatially isolated heparin microgels; when located within a microporous annealed particle (MAP) scaffold, they can sequester and release growth factors. Importantly, we present the first agent-based model of MAP scaffolds to optimize the ratio of heparin microislands within the scaffold to promote endothelial cell migration. We demonstrate this model can accurately predict trends in vitro, thus opening a new avenue of research to aid in the design of MAP scaffolds.
Collapse
|
15
|
Shou Y, Johnson SC, Quek YJ, Li X, Tay A. Integrative lymph node-mimicking models created with biomaterials and computational tools to study the immune system. Mater Today Bio 2022; 14:100269. [PMID: 35514433 PMCID: PMC9062348 DOI: 10.1016/j.mtbio.2022.100269] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2022] [Revised: 04/16/2022] [Accepted: 04/18/2022] [Indexed: 11/17/2022]
Abstract
The lymph node (LN) is a vital organ of the lymphatic and immune system that enables timely detection, response, and clearance of harmful substances from the body. Each LN comprises of distinct substructures, which host a plethora of immune cell types working in tandem to coordinate complex innate and adaptive immune responses. An improved understanding of LN biology could facilitate treatment in LN-associated pathologies and immunotherapeutic interventions, yet at present, animal models, which often have poor physiological relevance, are the most popular experimental platforms. Emerging biomaterial engineering offers powerful alternatives, with the potential to circumvent limitations of animal models, for in-depth characterization and engineering of the lymphatic and adaptive immune system. In addition, mathematical and computational approaches, particularly in the current age of big data research, are reliable tools to verify and complement biomaterial works. In this review, we first discuss the importance of lymph node in immunity protection followed by recent advances using biomaterials to create in vitro/vivo LN-mimicking models to recreate the lymphoid tissue microstructure and microenvironment, as well as to describe the related immuno-functionality for biological investigation. We also explore the great potential of mathematical and computational models to serve as in silico supports. Furthermore, we suggest how both in vitro/vivo and in silico approaches can be integrated to strengthen basic patho-biological research, translational drug screening and clinical personalized therapies. We hope that this review will promote synergistic collaborations to accelerate progress of LN-mimicking systems to enhance understanding of immuno-complexity.
Collapse
|
16
|
Abstract
Understanding the pathophysiology of lung fibrosis is of paramount importance to elaborate targeted and effective therapies. As it onsets, the randomly accumulating extracellular matrix (ECM) breaks the symmetry of the branching lung structure. Interestingly, similar pathways have been reported for both idiopathic pulmonary fibrosis and radiation-induced lung fibrosis (RILF). Individuals suffering from the disease, the worldwide incidence of which is growing, have poor prognosis and a short mean survival time. In this context, mathematical and computational models have the potential to shed light on key underlying pathological mechanisms, shorten the time needed for clinical trials, parallelize hypotheses testing, and improve personalized drug development. Agent-based modeling (ABM) has proven to be a reliable and versatile simulation tool, whose features make it a good candidate for recapitulating emergent behaviors in heterogeneous systems, such as those found at multiple scales in the human body. In this paper, we detail the implementation of a 3D agent-based model of lung fibrosis using a novel simulation platform, namely, BioDynaMo, and prove that it can qualitatively and quantitatively reproduce published results. Furthermore, we provide additional insights on late-fibrosis patterns through ECM density distribution histograms. The model recapitulates key intercellular mechanisms, while cell numbers and types are embodied by alveolar segments that act as agents and are spatially arranged by a custom algorithm. Finally, our model may hold potential for future applications in the context of lung disorders, ranging from RILF (by implementing radiation-induced cell damage mechanisms) to COVID-19 and inflammatory diseases (such as asthma or chronic obstructive pulmonary disease).
Collapse
|
17
|
Myers PJ, Lee SH, Lazzara MJ. MECHANISTIC AND DATA-DRIVEN MODELS OF CELL SIGNALING: TOOLS FOR FUNDAMENTAL DISCOVERY AND RATIONAL DESIGN OF THERAPY. CURRENT OPINION IN SYSTEMS BIOLOGY 2021; 28:100349. [PMID: 35935921 PMCID: PMC9348571 DOI: 10.1016/j.coisb.2021.05.010] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
A full understanding of cell signaling processes requires knowledge of protein structure/function relationships, protein-protein interactions, and the abilities of pathways to control phenotypes. Computational models offer a valuable framework for integrating that knowledge to predict the effects of system perturbations and interventions in health and disease. Whereas mechanistic models are well suited for understanding the biophysical basis for signal transduction and principles of therapeutic design, data-driven models are particularly suited to distill complex signaling relationships among samples and between multivariate signaling changes and phenotypes. Both approaches have limitations and provide incomplete representations of signaling biology, but their careful implementation and integration can provide new understanding for how manipulating system variables impacts cellular decisions.
Collapse
Affiliation(s)
- Paul J. Myers
- Department of Chemical Engineering, Charlottesville, VA 22904
| | - Sung Hyun Lee
- Department of Chemical Engineering, Charlottesville, VA 22904
| | - Matthew J. Lazzara
- Department of Chemical Engineering, Charlottesville, VA 22904
- Department of Biomedical Engineering University of Virginia, Charlottesville, VA 22904
| |
Collapse
|
18
|
Corti A, Colombo M, Migliavacca F, Rodriguez Matas JF, Casarin S, Chiastra C. Multiscale Computational Modeling of Vascular Adaptation: A Systems Biology Approach Using Agent-Based Models. Front Bioeng Biotechnol 2021; 9:744560. [PMID: 34796166 PMCID: PMC8593007 DOI: 10.3389/fbioe.2021.744560] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2021] [Accepted: 10/04/2021] [Indexed: 12/20/2022] Open
Abstract
The widespread incidence of cardiovascular diseases and associated mortality and morbidity, along with the advent of powerful computational resources, have fostered an extensive research in computational modeling of vascular pathophysiology field and promoted in-silico models as a support for biomedical research. Given the multiscale nature of biological systems, the integration of phenomena at different spatial and temporal scales has emerged to be essential in capturing mechanobiological mechanisms underlying vascular adaptation processes. In this regard, agent-based models have demonstrated to successfully embed the systems biology principles and capture the emergent behavior of cellular systems under different pathophysiological conditions. Furthermore, through their modular structure, agent-based models are suitable to be integrated with continuum-based models within a multiscale framework that can link the molecular pathways to the cell and tissue levels. This can allow improving existing therapies and/or developing new therapeutic strategies. The present review examines the multiscale computational frameworks of vascular adaptation with an emphasis on the integration of agent-based approaches with continuum models to describe vascular pathophysiology in a systems biology perspective. The state-of-the-art highlights the current gaps and limitations in the field, thus shedding light on new areas to be explored that may become the future research focus. The inclusion of molecular intracellular pathways (e.g., genomics or proteomics) within the multiscale agent-based modeling frameworks will certainly provide a great contribution to the promising personalized medicine. Efforts will be also needed to address the challenges encountered for the verification, uncertainty quantification, calibration and validation of these multiscale frameworks.
Collapse
Affiliation(s)
- Anna Corti
- Laboratory of Biological Structure Mechanics (LaBS), Department of Chemistry, Materials and Chemical Engineering "Giulio Natta", Politecnico di Milano, Milan, Italy
| | - Monika Colombo
- Laboratory of Biological Structure Mechanics (LaBS), Department of Chemistry, Materials and Chemical Engineering "Giulio Natta", Politecnico di Milano, Milan, Italy.,Institute for Chemical and Bioengineering, Department of Chemistry and Applied Biosciences, ETH Zürich, Switzerland
| | - Francesco Migliavacca
- Laboratory of Biological Structure Mechanics (LaBS), Department of Chemistry, Materials and Chemical Engineering "Giulio Natta", Politecnico di Milano, Milan, Italy
| | - Jose Felix Rodriguez Matas
- Laboratory of Biological Structure Mechanics (LaBS), Department of Chemistry, Materials and Chemical Engineering "Giulio Natta", Politecnico di Milano, Milan, Italy
| | - Stefano Casarin
- Department of Surgery, Houston Methodist Hospital, Houston, TX, United States.,Center for Computational Surgery, Houston Methodist Research Institute, Houston, TX, United States.,Houston Methodist Academic Institute, Houston, TX, United States
| | - Claudio Chiastra
- Laboratory of Biological Structure Mechanics (LaBS), Department of Chemistry, Materials and Chemical Engineering "Giulio Natta", Politecnico di Milano, Milan, Italy.,PoliToMed Lab, Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Turin, Italy
| |
Collapse
|
19
|
Abstract
Multiscale computational modeling aims to connect the complex networks of effects at different length and/or time scales. For example, these networks often include intracellular molecular signaling, crosstalk, and other interactions between neighboring cell populations, and higher levels of emergent phenomena across different regions of tissues and among collections of tissues or organs interacting with each other in the whole body. Recent applications of multiscale modeling across intracellular, cellular, and/or tissue levels are highlighted here. These models incorporated the roles of biochemical and biomechanical modulation in processes that are implicated in the mechanisms of several diseases including fibrosis, joint and bone diseases, respiratory infectious diseases, and cancers.
Collapse
|
20
|
Lafci Büyükkahraman M, Sabine GK, Kojouharov HV, Chen-Charpentier BM, McMahan SR, Liao J. Using models to advance medicine: mathematical modeling of post-myocardial infarction left ventricular remodeling. Comput Methods Biomech Biomed Engin 2021; 25:298-307. [PMID: 34266318 DOI: 10.1080/10255842.2021.1953487] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Abstract
The heart is an organ with limited capacity for regeneration and repair. The irreversible cell death and corresponding diminished ability of the heart to repair after myocardial infarction (MI), is a leading cause of morbidity and mortality worldwide. In this paper, a new mathematical model is presented to study the left ventricular (LV) remodeling and associated events after MI. The model accurately describes and predicts the interactions among heart cells and the immune system post-MI in the absence of medical interventions. The resulting system of nonlinear ordinary differential equations is studied both analytically and numerically in order to demonstrate the functionality and performance of the new model. To the best of our knowledge, this model is the only one of its kind to consider and correctly apply all of the known factors in diseased heart LV modeling. This model has the potential to provide researchers with a predictive computational tool to better understand the MI pathology and develop various cell-based treatment options, with benefits of lowering the cost and reducing the development time.
Collapse
Affiliation(s)
- Mehtap Lafci Büyükkahraman
- Department of Mathematics, The University of Texas at Arlington, Arlington, TX, USA.,Department of Mathematics, Uşak University, Uşak, Turkey
| | - Gavin K Sabine
- Department of Mathematics, The University of Texas at Arlington, Arlington, TX, USA
| | - Hristo V Kojouharov
- Department of Mathematics, The University of Texas at Arlington, Arlington, TX, USA
| | | | - Sara R McMahan
- Department of Bioengineering, The University of Texas at Arlington, Arlington, TX, USA
| | - Jun Liao
- Department of Bioengineering, The University of Texas at Arlington, Arlington, TX, USA
| |
Collapse
|
21
|
Gregg RW, Shabnam F, Shoemaker JE. Agent-based modeling reveals benefits of heterogeneous and stochastic cell populations during cGAS-mediated IFNβ production. Bioinformatics 2021; 37:1428-1434. [PMID: 33196784 DOI: 10.1093/bioinformatics/btaa969] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2020] [Revised: 10/13/2020] [Accepted: 11/04/2020] [Indexed: 01/25/2023] Open
Abstract
MOTIVATION The cGAS pathway is a component of the innate immune system responsible for the detection of pathogenic DNA and upregulation of interferon beta (IFNβ). Experimental evidence shows that IFNβ signaling occurs in highly heterogeneous cells and is stochastic in nature; however, the benefits of these attributes remain unclear. To investigate how stochasticity and heterogeneity affect IFNβ production, an agent-based model is developed to simulate both DNA transfection and viral infection. RESULTS We show that heterogeneity can enhance IFNβ responses during infection. Furthermore, by varying the degree of IFNβ stochasticity, we find that only a percentage of cells (20-30%) need to respond during infection. Going beyond this range provides no additional protection against cell death or reduction of viral load. Overall, these simulations suggest that heterogeneity and stochasticity are important for moderating immune potency while minimizing cell death during infection. AVAILABILITY AND IMPLEMENTATION Model repository is available at: https://github.com/ImmuSystems-Lab/AgentBasedModel-cGASPathway. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Robert W Gregg
- Department of Chemical and Petroleum Engineering, 15260, Pittsburgh, PA 15260, USA
| | - Fathima Shabnam
- Department of Chemical and Petroleum Engineering, 15260, Pittsburgh, PA 15260, USA
| | - Jason E Shoemaker
- Department of Chemical and Petroleum Engineering, 15260, Pittsburgh, PA 15260, USA.,McGowan Institute for Regenerative Medicine, 15219, Pittsburgh, PA 15260, USA.,Department of Computational and Systems Biology, University of Pittsburgh, Pittsburgh, PA 15260, USA
| |
Collapse
|
22
|
Westman AM, Peirce SM, Christ GJ, Blemker SS. Agent-based model provides insight into the mechanisms behind failed regeneration following volumetric muscle loss injury. PLoS Comput Biol 2021; 17:e1008937. [PMID: 33970905 PMCID: PMC8110270 DOI: 10.1371/journal.pcbi.1008937] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2020] [Accepted: 04/01/2021] [Indexed: 12/22/2022] Open
Abstract
Skeletal muscle possesses a remarkable capacity for repair and regeneration following a variety of injuries. When successful, this highly orchestrated regenerative process requires the contribution of several muscle resident cell populations including satellite stem cells (SSCs), fibroblasts, macrophages and vascular cells. However, volumetric muscle loss injuries (VML) involve simultaneous destruction of multiple tissue components (e.g., as a result of battlefield injuries or vehicular accidents) and are so extensive that they exceed the intrinsic capability for scarless wound healing and result in permanent cosmetic and functional deficits. In this scenario, the regenerative process fails and is dominated by an unproductive inflammatory response and accompanying fibrosis. The failure of current regenerative therapeutics to completely restore functional muscle tissue is not surprising considering the incomplete understanding of the cellular mechanisms that drive the regeneration response in the setting of VML injury. To begin to address this profound knowledge gap, we developed an agent-based model to predict the tissue remodeling response following surgical creation of a VML injury. Once the model was able to recapitulate key aspects of the tissue remodeling response in the absence of repair, we validated the model by simulating the tissue remodeling response to VML injury following implantation of either a decellularized extracellular matrix scaffold or a minced muscle graft. The model suggested that the SSC microenvironment and absence of pro-differentiation SSC signals were the most important aspects of failed muscle regeneration in VML injuries. The major implication of this work is that agent-based models may provide a much-needed predictive tool to optimize the design of new therapies, and thereby, accelerate the clinical translation of regenerative therapeutics for VML injuries.
Collapse
Affiliation(s)
- Amanda M. Westman
- Biomedical Engineering, University of Virginia, Charlottesville, Virginia, United States of America
| | - Shayn M. Peirce
- Biomedical Engineering, University of Virginia, Charlottesville, Virginia, United States of America
- Ophthalmology, University of Virginia, Charlottesville, Virginia, United States of America
| | - George J. Christ
- Biomedical Engineering, University of Virginia, Charlottesville, Virginia, United States of America
- Orthopaedic Surgery, University of Virginia, Charlottesville, Virginia, United States of America
- * E-mail: (GJC); (SSB)
| | - Silvia S. Blemker
- Biomedical Engineering, University of Virginia, Charlottesville, Virginia, United States of America
- Ophthalmology, University of Virginia, Charlottesville, Virginia, United States of America
- Orthopaedic Surgery, University of Virginia, Charlottesville, Virginia, United States of America
- Mechanical and Aerospace Engineering, University of Virginia, Charlottesville, Virginia, United States of America
- * E-mail: (GJC); (SSB)
| |
Collapse
|
23
|
Bloomingdale P, Karelina T, Cirit M, Muldoon SF, Baker J, McCarty WJ, Geerts H, Macha S. Quantitative systems pharmacology in neuroscience: Novel methodologies and technologies. CPT-PHARMACOMETRICS & SYSTEMS PHARMACOLOGY 2021; 10:412-419. [PMID: 33719204 PMCID: PMC8129713 DOI: 10.1002/psp4.12607] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/01/2020] [Revised: 12/02/2020] [Accepted: 02/03/2021] [Indexed: 11/12/2022]
Abstract
The development and application of quantitative systems pharmacology models in neuroscience have been modest relative to other fields, such as oncology and immunology, which may reflect the complexity of the brain. Technological and methodological advancements have enhanced the quantitative understanding of brain physiology and pathophysiology and the effects of pharmacological interventions. To maximize the knowledge gained from these novel data types, pharmacometrics modelers may need to expand their toolbox to include additional mathematical and statistical frameworks. A session was held at the 10th annual American Conference on Pharmacometrics (ACoP10) to highlight several recent advancements in quantitative and systems neuroscience. In this mini‐review, we provide a brief overview of technological and methodological advancements in the neuroscience therapeutic area that were discussed during the session and how these can be leveraged with quantitative systems pharmacology modeling to enhance our understanding of neurological diseases. Microphysiological systems using human induced pluripotent stem cells (IPSCs), digital biomarkers, and large‐scale imaging offer more clinically relevant experimental datasets, enhanced granularity, and a plethora of data to potentially improve the preclinical‐to‐clinical translation of therapeutics. Network neuroscience methodologies combined with quantitative systems models of neurodegenerative disease could help bridge the gap between cellular and molecular alterations and clinical end points through the integration of information on neural connectomics. Additional topics, such as the neuroimmune system, microbiome, single‐cell transcriptomic technologies, and digital device biomarkers, are discussed in brief.
Collapse
Affiliation(s)
- Peter Bloomingdale
- Quantitative Pharmacology and Pharmacometrics, Merck & Co. Inc, Kenilworth, New Jersey, USA
| | | | - Murat Cirit
- Javelin Biotech, Inc, Woburn, Massachusetts, USA
| | - Sarah F Muldoon
- Mathematics Department, CDSE Program, Neuroscience Program, University at Buffalo, SUNY, Buffalo, New York, USA
| | - Justin Baker
- Harvard Medical School, Boston, Massachusetts, USA
| | | | | | - Sreeraj Macha
- Quantitative Pharmacology, Sanofi, Bridgewater, New Jersey, USA
| |
Collapse
|
24
|
Hall C, Gehmlich K, Denning C, Pavlovic D. Complex Relationship Between Cardiac Fibroblasts and Cardiomyocytes in Health and Disease. J Am Heart Assoc 2021; 10:e019338. [PMID: 33586463 PMCID: PMC8174279 DOI: 10.1161/jaha.120.019338] [Citation(s) in RCA: 77] [Impact Index Per Article: 25.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Cardiac fibroblasts are the primary cell type responsible for deposition of extracellular matrix in the heart, providing support to the contracting myocardium and contributing to a myriad of physiological signaling processes. Despite the importance of fibrosis in processes of wound healing, excessive fibroblast proliferation and activation can lead to pathological remodeling, driving heart failure and the onset of arrhythmias. Our understanding of the mechanisms driving the cardiac fibroblast activation and proliferation is expanding, and evidence for their direct and indirect effects on cardiac myocyte function is accumulating. In this review, we focus on the importance of the fibroblast-to-myofibroblast transition and the cross talk of cardiac fibroblasts with cardiac myocytes. We also consider the current use of models used to explore these questions.
Collapse
Affiliation(s)
- Caitlin Hall
- Institute of Cardiovascular Sciences University of Birmingham United Kingdom
| | - Katja Gehmlich
- Institute of Cardiovascular Sciences University of Birmingham United Kingdom.,Division of Cardiovascular Medicine Radcliffe Department of Medicine and British Heart Foundation Centre of Research Excellence Oxford University of Oxford United Kingdom
| | - Chris Denning
- Biodiscovery Institute University of Nottingham United Kingdom
| | - Davor Pavlovic
- Institute of Cardiovascular Sciences University of Birmingham United Kingdom
| |
Collapse
|
25
|
Irons L, Latorre M, Humphrey JD. From Transcript to Tissue: Multiscale Modeling from Cell Signaling to Matrix Remodeling. Ann Biomed Eng 2021; 49:1701-1715. [PMID: 33415527 DOI: 10.1007/s10439-020-02713-8] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2020] [Accepted: 12/15/2020] [Indexed: 02/06/2023]
Abstract
Tissue-level biomechanical properties and function derive from underlying cell signaling, which regulates mass deposition, organization, and removal. Here, we couple two existing modeling frameworks to capture associated multiscale interactions-one for vessel-level growth and remodeling and one for cell-level signaling-and illustrate utility by simulating aortic remodeling. At the vessel level, we employ a constrained mixture model describing turnover of individual wall constituents (elastin, intramural cells, and collagen), which has proven useful in predicting diverse adaptations as well as disease progression using phenomenological constitutive relations. Nevertheless, we now seek an improved mechanistic understanding of these processes; we replace phenomenological relations in the mixture model with a logic-based signaling model, which yields a system of ordinary differential equations predicting changes in collagen synthesis, matrix metalloproteinases, and cell proliferation in response to altered intramural stress, wall shear stress, and exogenous angiotensin II. This coupled approach promises improved understanding of the role of cell signaling in achieving tissue homeostasis and allows us to model feedback between vessel mechanics and cell signaling. We verify our model predictions against data from the hypertensive murine infrarenal abdominal aorta as well as results from validated phenomenological models, and consider effects of noisy signaling and heterogeneous cell populations.
Collapse
Affiliation(s)
- Linda Irons
- Department of Biomedical Engineering, Yale University, New Haven, CT, USA.
| | - Marcos Latorre
- Department of Biomedical Engineering, Yale University, New Haven, CT, USA
| | - Jay D Humphrey
- Department of Biomedical Engineering, Yale University, New Haven, CT, USA
| |
Collapse
|
26
|
Joslyn LR, Kirschner DE, Linderman JJ. CaliPro: A Calibration Protocol That Utilizes Parameter Density Estimation to Explore Parameter Space and Calibrate Complex Biological Models. Cell Mol Bioeng 2020; 14:31-47. [PMID: 33643465 DOI: 10.1007/s12195-020-00650-z] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2020] [Accepted: 09/02/2020] [Indexed: 12/15/2022] Open
Abstract
Introduction Mathematical and computational modeling have a long history of uncovering mechanisms and making predictions for biological systems. However, to create a model that can provide relevant quantitative predictions, models must first be calibrated by recapitulating existing biological datasets from that system. Current calibration approaches may not be appropriate for complex biological models because: 1) many attempt to recapitulate only a single aspect of the experimental data (such as a median trend) or 2) Bayesian techniques require specification of parameter priors and likelihoods to experimental data that cannot always be confidently assigned. A new calibration protocol is needed to calibrate complex models when current approaches fall short. Methods Herein, we develop CaliPro, an iterative, model-agnostic calibration protocol that utilizes parameter density estimation to refine parameter space and calibrate to temporal biological datasets. An important aspect of CaliPro is the user-defined pass set definition, which specifies how the model might successfully recapitulate experimental data. We define the appropriate settings to use CaliPro. Results We illustrate the usefulness of CaliPro through four examples including predator-prey, infectious disease transmission, and immune response models. We show that CaliPro works well for both deterministic, continuous model structures as well as stochastic, discrete models and illustrate that CaliPro can work across diverse calibration goals. Conclusions We present CaliPro, a new method for calibrating complex biological models to a range of experimental outcomes. In addition to expediting calibration, CaliPro may be useful in already calibrated parameter spaces to target and isolate specific model behavior for further analysis.
Collapse
Affiliation(s)
- Louis R Joslyn
- Department of Chemical Engineering, University of Michigan, G045W NCRC B28, 2800 Plymouth Rd, Ann Arbor, MI 48109-2136 USA.,Department of Microbiology and Immunology, University of Michigan Medical School, 1150 W Medical Center Drive, 5641 Medical Science II, Ann Arbor, MI 48109-5620 USA
| | - Denise E Kirschner
- Department of Microbiology and Immunology, University of Michigan Medical School, 1150 W Medical Center Drive, 5641 Medical Science II, Ann Arbor, MI 48109-5620 USA
| | - Jennifer J Linderman
- Department of Chemical Engineering, University of Michigan, G045W NCRC B28, 2800 Plymouth Rd, Ann Arbor, MI 48109-2136 USA
| |
Collapse
|
27
|
Leonard-Duke J, Evans S, Hannan RT, Barker TH, Bates JHT, Bonham CA, Moore BB, Kirschner DE, Peirce SM. Multi-scale models of lung fibrosis. Matrix Biol 2020; 91-92:35-50. [PMID: 32438056 DOI: 10.1016/j.matbio.2020.04.003] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2020] [Revised: 03/13/2020] [Accepted: 04/15/2020] [Indexed: 02/08/2023]
Abstract
The architectural complexity of the lung is crucial to its ability to function as an organ of gas exchange; the branching tree structure of the airways transforms the tracheal cross-section of only a few square centimeters to a blood-gas barrier with a surface area of tens of square meters and a thickness on the order of a micron or less. Connective tissue comprised largely of collagen and elastic fibers provides structural integrity for this intricate and delicate system. Homeostatic maintenance of this connective tissue, via a balance between catabolic and anabolic enzyme-driven processes, is crucial to life. Accordingly, when homeostasis is disrupted by the excessive production of connective tissue, lung function deteriorates rapidly with grave consequences leading to chronic lung conditions such as pulmonary fibrosis. Understanding how pulmonary fibrosis develops and alters the link between lung structure and function is crucial for diagnosis, prognosis, and therapy. Further information gained could help elaborate how the healing process breaks down leading to chronic disease. Our understanding of fibrotic disease is greatly aided by the intersection of wet lab studies and mathematical and computational modeling. In the present review we will discuss how multi-scale modeling has facilitated our understanding of pulmonary fibrotic disease as well as identified opportunities that remain open and have produced techniques that can be incorporated into this field by borrowing approaches from multi-scale models of fibrosis beyond the lung.
Collapse
Affiliation(s)
- Julie Leonard-Duke
- Department of Biomedical Engineering, University of Virginia, Charlottesville, VA 22908, USA
| | - Stephanie Evans
- Department of Microbiology and Immunology, University of Michigan Medical School, Ann Arbor, MI 48109, USA
| | - Riley T Hannan
- Department of Pathology, University of Virginia, Charlottesville, VA 22908, USA
| | - Thomas H Barker
- Department of Biomedical Engineering, University of Virginia, Charlottesville, VA 22908, USA
| | - Jason H T Bates
- Department of Medicine, Vermont Lung Center, University of Vermont College of Medicine, Burlington, VT 05405, USA
| | - Catherine A Bonham
- Division of Pulmonary and Critical Care Medicine, University of Virginia, Charlottesville VA 22908, USA
| | - Bethany B Moore
- Division of Pulmonary and Critical Care Medicine, Department of Internal Medicine, and Department of Microbiology and Immunology, University of Michigan Medical Center, Ann Arbor, MI, 48109, USA
| | - Denise E Kirschner
- Department of Microbiology and Immunology, University of Michigan Medical School, Ann Arbor, MI 48109, USA
| | - Shayn M Peirce
- Department of Biomedical Engineering, University of Virginia, Charlottesville, VA 22908, USA; Robert M. Berne Cardiovascular Research Center, University of Virginia, Charlottesville, VA 22908, USA.
| |
Collapse
|