In this article we propose a novel mathematical description of biomass growth that combines poroelastic theory of mixtures and cellular population models. The formulation, potentially applicable to general mechanobiological processes, is here used to study the engineered cultivation in bioreactors of articular chondrocytes, a process of Regenerative Medicine characterized by a complex interaction among spatial scales (from nanometers to centimeters), temporal scales (from seconds to weeks) and biophysical phenomena (fluid-controlled nutrient transport, delivery and consumption; mechanical deformation of a multiphase porous medium). The principal contribution of this research is the inclusion of the concept of cellular “force isotropy” as one of the main factors influencing cellular activity. In this description, the induced cytoskeletal tensional states trigger signalling transduction cascades regulating functional cell behavior. This mechanims is modeled by a parameter which estimates the influence of local force isotropy by the norm of the deviatoric part of the total stress tensor. According to the value of the estimator, isotropic mechanical conditions are assumed to be the promoting factor of extracellular matrix production whereas anisotropic conditions are assumed to promote cell proliferation. The resulting mathematical formulation is a coupled system of nonlinear partial differential equations comprising: conservation laws for mass and linear momentum of the growing biomass; advection–diffusion–reaction laws for nutrient (oxygen) transport, delivery and consumption; and kinetic laws for cellular population dynamics. To develop a reliable computational tool for the simulation of the engineered tissue growth process the nonlinear differential problem is numerically solved by: (1) temporal semidiscretization; (2) linearization via a fixed-point map; and (3) finite element spatial approximation. The biophysical accuracy of the mechanobiological model is assessed in the analysis of a simplified 1D geometrical setting. Simulation results show that: (1) isotropic/anisotropic conditions are strongly influenced by both maximum cell specific growth rate and mechanical boundary conditions enforced at the interface between the biomass construct and the interstitial fluid; (2) experimentally measured features of cultivated articular chondrocytes, such as the early proliferation phase and the delayed extracellular matrix production, are well described by the computed spatial and temporal evolutions of cellular populations.

A poroelastic mixture model of mechanobiological processes in biomass growth : theory and application to tissue engineering / R. Sacco, P. Causin, C. Lelli, M..T. Raimondi. - In: MECCANICA. - ISSN 1572-9648. - 52:14(2017 Nov), pp. 3273-3297. (Intervento presentato al 22. convegno Congresso dell’Associazione Italiana di meccanica teorica e applicata (AIMETA): two mini-symposia Soft active materials (Mechanics and Materials 2015, GMA) and Mechanobiology of the cell and morphogenesis of living matter tenutosi a Genova nel 2015) [10.1007/s11012-017-0638-9].

A poroelastic mixture model of mechanobiological processes in biomass growth : theory and application to tissue engineering

P. Causin;
2017

Abstract

In this article we propose a novel mathematical description of biomass growth that combines poroelastic theory of mixtures and cellular population models. The formulation, potentially applicable to general mechanobiological processes, is here used to study the engineered cultivation in bioreactors of articular chondrocytes, a process of Regenerative Medicine characterized by a complex interaction among spatial scales (from nanometers to centimeters), temporal scales (from seconds to weeks) and biophysical phenomena (fluid-controlled nutrient transport, delivery and consumption; mechanical deformation of a multiphase porous medium). The principal contribution of this research is the inclusion of the concept of cellular “force isotropy” as one of the main factors influencing cellular activity. In this description, the induced cytoskeletal tensional states trigger signalling transduction cascades regulating functional cell behavior. This mechanims is modeled by a parameter which estimates the influence of local force isotropy by the norm of the deviatoric part of the total stress tensor. According to the value of the estimator, isotropic mechanical conditions are assumed to be the promoting factor of extracellular matrix production whereas anisotropic conditions are assumed to promote cell proliferation. The resulting mathematical formulation is a coupled system of nonlinear partial differential equations comprising: conservation laws for mass and linear momentum of the growing biomass; advection–diffusion–reaction laws for nutrient (oxygen) transport, delivery and consumption; and kinetic laws for cellular population dynamics. To develop a reliable computational tool for the simulation of the engineered tissue growth process the nonlinear differential problem is numerically solved by: (1) temporal semidiscretization; (2) linearization via a fixed-point map; and (3) finite element spatial approximation. The biophysical accuracy of the mechanobiological model is assessed in the analysis of a simplified 1D geometrical setting. Simulation results show that: (1) isotropic/anisotropic conditions are strongly influenced by both maximum cell specific growth rate and mechanical boundary conditions enforced at the interface between the biomass construct and the interstitial fluid; (2) experimentally measured features of cultivated articular chondrocytes, such as the early proliferation phase and the delayed extracellular matrix production, are well described by the computed spatial and temporal evolutions of cellular populations.
tissue engineering; mechanobiology; mathematical modeling; mixture growth theory; mass transport; continuum mechanics; numerical simulation
Settore MAT/08 - Analisi Numerica
nov-2017
feb-2017
Article (author)
File in questo prodotto:
File Dimensione Formato  
MECC-D-16-00273_revised.pdf

accesso aperto

Tipologia: Post-print, accepted manuscript ecc. (versione accettata dall'editore)
Dimensione 973.85 kB
Formato Adobe PDF
973.85 kB Adobe PDF Visualizza/Apri
10.1007_s11012-017-0638-9.pdf

accesso aperto

Tipologia: Publisher's version/PDF
Dimensione 2.02 MB
Formato Adobe PDF
2.02 MB Adobe PDF Visualizza/Apri
Pubblicazioni consigliate

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/2434/543349
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 20
  • ???jsp.display-item.citation.isi??? 15
social impact