An in-silico approach to meniscus tissue regeneration: Modeling, numerical simulation, and experimental analysis

Elise Grosjean,Alex Keilmann,Henry Jäger,Shimi Mohanan,Claudia Redenbach,Bernd Simeon,Christina Surulescu,Luisa de Roy,Andreas Seitz,Graciosa Teixeira,Martin Dauner,Carsten Linti,Günter Schmidt

2024-03-09

Abstract:We develop a model the dynamics of human mesenchymal stem cells (hMSCs) and chondrocytes evolving in a nonwoven polyethylene terephtalate (PET) scaffold impregnated with hyaluron and supplied with a differentiation medium. The scaffold and the cells are assumed to be contained in a bioreactor with fluid perfusion. The differentiation of hMSCs into chondrocytes favors the production of extracellular matrix (ECM) and is influenced by fluid stress. The model takes deformations of ECM and PET scaffold into account. The scaffold structure is explicitly included by statistical assessment of the fibre distribution from CT images. The effective macroscopic equations are obtained by appropriate upscaling from dynamics on lower (microscopic and mesoscopic) scales and feature in the motility terms an explicit cell diffusion tensor encoding the assessed anisotropic scaffold structure. Numerical simulations show its influence on the overall cell and tissue dynamics.

Tissues and Organs,Analysis of PDEs,Numerical Analysis,Cell Behavior,Quantitative Methods

What problem does this paper attempt to address?

### Problems the paper attempts to solve

This paper aims to solve the problem of meniscus tissue regeneration through a computational method (i.e., the **in - silico** method). Specifically, the researchers developed a mathematical model to describe the dynamic behaviors of human mesenchymal stem cells (hMSCs) and chondrocytes in non - woven polyethylene terephthalate (PET) scaffolds. This scaffold was impregnated with hyaluronic acid and provided a differentiation medium. The model took into account the influence of fluid stress on cell differentiation, as well as the deformation of the extracellular matrix (ECM) and the PET scaffold. In addition, the model also explicitly included the scaffold structure by statistically evaluating the fiber distribution in CT images. Finally, numerical simulations showed the influence of these factors on the overall cell and tissue dynamics.

### Main research contents

1. **Model development**:
- A multi - scale model was developed to describe the dynamic behaviors of hMSCs and chondrocytes in PET scaffolds from the micro - scale, meso - scale to the macro - scale.
- The influence of fluid stress on cell differentiation and the deformation of ECM and scaffolds were considered.
- The scaffold structure was explicitly included by statistically evaluating the fiber distribution in CT images.

2. **Experimental verification**:
- Experiments were carried out. In vitro studies were conducted using non - woven PET scaffolds in perfusion chambers to verify the validity of the model.
- The experimental data provided detailed information on the key stimulatory factors for cell proliferation, migration and differentiation.

3. **Numerical simulation**:
- The predictive ability of the model was demonstrated through numerical simulations, especially the cell and tissue dynamics under different environmental parameters.

### Key features of the model

- **Multi - scale modeling**: The model starts from the receptor - binding kinetics of a single cell, passes through the cell distribution function at the meso - scale, and finally obtains the population density equation at the macro - scale.
- **Mechanical effects**: The mechanical stresses caused by fluid flow in the bioreactor were considered, and these stresses have an important impact on cell differentiation and tissue regeneration.
- **Statistical evaluation**: The fiber distribution was statistically evaluated by CT images and explicitly included in the model to more accurately describe the influence of the scaffold structure on cell behavior.

### Mathematical model

The core of the model is a set of reaction - diffusion - chemotaxis equations (RDTEs), which describe the dynamic behaviors of hMSCs and chondrocytes. The specific equations are as follows:

\[
\frac{\partial c_1}{\partial t}-\nabla\cdot(\nabla\cdot(D_1 c_1))+\nabla\cdot\left(\frac{k'\lambda_{11}}{B(h,\tau)^2(B(h,\tau)+\lambda_{10})}D_1\nabla B(h,\tau)c_1\right)=-\alpha_1(\chi,S)c_1+\alpha_2(\chi,S)\frac{\omega_1}{\omega_2}c_2+\beta c_1\left(1 - \frac{c_1}{C_1^*}-\frac{c_2}{C_2^*}\right)
\]

\[
\frac{\partial c_2}{\partial t}-\nabla\cdot(\nabla\cdot(D_2 c_2))=\alpha_1(\chi,S)\frac{\omega_2}{\omega_1}c_1-\alpha_2(\chi,S)c_2
\]

where:
- \(c_1\) and \(c_2\) represent the macroscopic densities of hMSCs and chondrocytes respectively.
- \(D_1\) and \(D_2\) are diffusion tensors.
- \(\alpha_1(\chi,S)\) and \(\alpha_2(\chi,S)\) represent the rates of hMSCs differentiating into chondrocytes and chondrocytes dedifferentiating respectively.
- \(\beta\) is the growth rate of hMSCs.
- \(C_1\)