Société de Biomécanique Young Investigator Award 2021: Numerical investigation of the time-dependent stress-strain mechanical behaviour of skeletal muscle tissue in the context of pressure ulcer prevention.

BACKGROUND: Pressure-induced tissue strain is one major pathway for Pressure Ulcer development and, especially, Deep Tissue Injury. Biomechanical investigation of the time-dependent stress-strain mechanical behaviour of skeletal muscle tissue is therefore essential. In the literature, a viscoelastic formulation is generally assumed for the experimental characterization of skeletal muscles, with the limitation that the underlying physical mechanisms that give rise to the time dependent stress-strain behaviour are not known. The objective of this study is to explore the capability of poroelasticity to reproduce the apparent viscoelastic behaviour of passive muscle tissue under confined compression. METHODS: Experimental stress-relaxation response of 31 cylindrical porcine samples tested under fast and slow confined compression by Vaidya and collaborators were used. An axisymmetric Finite Element model was developed in ABAQUS and, for each sample a one-to-one inverse analysis was performed to calibrate the specimen-specific constitutive parameters, namely, the drained Young's modulus, the void ratio, hydraulic permeability, the Poisson's ratio, the solid grain's and fluid's bulk moduli. FINDINGS: The peak stress and consolidation were recovered for most of the samples (N=25) by the poroelastic model (normalised root-mean-square error ≤0.03 for fast and slow confined compression conditions). INTERPRETATION: The strength of the proposed model is its fewer number of variables (N=6 for the proposed poroelastic model versus N=18 for the viscohyperelastic model proposed by Vaidya and collaborators). The incorporation of poroelasticity to clinical models of Pessure Ulcer formation could lead to more precise and mechanistic explorations of soft tissue injury risk factors.

Role of mechanical cues and hypoxia on the growth of tumor cells in strong and weak confinement: A dual in vitro-in silico approach.

Characterization of tumor growth dynamics is of major importance for cancer understanding. By contrast with phenomenological approaches, mechanistic modeling can facilitate disclosing underlying tumor mechanisms and lead to identification of physical factors affecting proliferation and invasive behavior. Current mathematical models are often formulated at the tissue or organ scale with the scope of a direct clinical usefulness. Consequently, these approaches remain empirical and do not allow gaining insight into the tumor properties at the scale of small cell aggregates. Here, experimental and numerical studies of the dynamics of tumor aggregates are performed to propose a physics-based mathematical model as a general framework to investigate tumor microenvironment. The quantitative data extracted from the cellular capsule technology microfluidic experiments allow a thorough quantitative comparison with in silico experiments. This dual approach demonstrates the relative impact of oxygen and external mechanical forces during the time course of tumor model progression.

A tumor growth model with deformable ECM.

Existing tumor growth models based on fluid analogy for the cells do not generally include the extracellular matrix (ECM), or if present, take it as rigid. The three-fluid model originally proposed by the authors and comprising tumor cells (TC), host cells (HC), interstitial fluid (IF) and an ECM, considered up to now only a rigid ECM in the applications. This limitation is here relaxed and the deformability of the ECM is investigated in detail. The ECM is modeled as a porous solid matrix with Green-elastic and elasto-visco-plastic material behavior within a large strain approach. Jauman and Truesdell objective stress measures are adopted together with the deformation rate tensor. Numerical results are first compared with those of a reference experiment of a multicellular tumor spheroid (MTS) growing in vitro, then three different tumor cases are studied: growth of an MTS in a decellularized ECM, growth of a spheroid in the presence of host cells and growth of a melanoma. The influence of the stiffness of the ECM is evidenced and comparison with the case of a rigid ECM is made. The processes in a deformable ECM are more rapid than in a rigid ECM and the obtained growth pattern differs. The reasons for this are due to the changes in porosity induced by the tumor growth. These changes are inhibited in a rigid ECM. This enhanced computational model emphasizes the importance of properly characterizing the biomechanical behavior of the malignant mass in all its components to correctly predict its temporal and spatial pattern evolution.

A two-phase model of plantar tissue: a step toward prediction of diabetic foot ulceration.

A new computational model, based on the thermodynamically constrained averaging theory, has been recently proposed to predict tumor initiation and proliferation. A similar mathematical approach is proposed here as an aid in diabetic ulcer prevention. The common aspects at the continuum level are the macroscopic balance equations governing the flow of the fluid phase, diffusion of chemical species, tissue mechanics, and some of the constitutive equations. The soft plantar tissue is modeled as a two-phase system: a solid phase consisting of the tissue cells and their extracellular matrix, and a fluid one (interstitial fluid and dissolved chemical species). The solid phase may become necrotic depending on the stress level and on the oxygen availability in the tissue. Actually, in diabetic patients, peripheral vascular disease impacts tissue necrosis; this is considered in the model via the introduction of an effective diffusion coefficient that governs transport of nutrients within the microvasculature. The governing equations of the mathematical model are discretized in space by the finite element method and in time domain using the θ-Wilson Method. While the full mathematical model is developed in this paper, the example is limited to the simulation of several gait cycles of a healthy foot.

The TNF-family cytokine TL1A promotes allergic immunopathology through group 2 innate lymphoid cells.

The tumor necrosis factor (TNF)-family cytokine TL1A (TNFSF15) costimulates T cells and promotes diverse T cell-dependent models of autoimmune disease through its receptor DR3. TL1A polymorphisms also confer susceptibility to inflammatory bowel disease. Here, we find that allergic pathology driven by constitutive TL1A expression depends on interleukin-13 (IL-13), but not on T, NKT, mast cells, or commensal intestinal flora. Group 2 innate lymphoid cells (ILC2) express surface DR3 and produce IL-13 and other type 2 cytokines in response to TL1A. DR3 is required for ILC2 expansion and function in the setting of T cell-dependent and -independent models of allergic disease. By contrast, DR3-deficient ILC2 can still differentiate, expand, and produce IL-13 when stimulated by IL-25 or IL-33, and mediate expulsion of intestinal helminths. These data identify costimulation of ILC2 as a novel function of TL1A important for allergic lung disease, and suggest that TL1A may be a therapeutic target in these settings.

A multiphase model for three-dimensional tumor growth.

Several mathematical formulations have analyzed the time-dependent behaviour of a tumor mass. However, most of these propose simplifications that compromise the physical soundness of the model. Here, multiphase porous media mechanics is extended to model tumor evolution, using governing equations obtained via the Thermodynamically Constrained Averaging Theory (TCAT). A tumor mass is treated as a multiphase medium composed of an extracellular matrix (ECM); tumor cells (TC), which may become necrotic depending on the nutrient concentration and tumor phase pressure; healthy cells (HC); and an interstitial fluid (IF) for the transport of nutrients. The equations are solved by a Finite Element method to predict the growth rate of the tumor mass as a function of the initial tumor-to-healthy cell density ratio, nutrient concentration, mechanical strain, cell adhesion and geometry. Results are shown for three cases of practical biological interest such as multicellular tumor spheroids (MTS) and tumor cords. First, the model is validated by experimental data for time-dependent growth of an MTS in a culture medium. The tumor growth pattern follows a biphasic behaviour: initially, the rapidly growing tumor cells tend to saturate the volume available without any significant increase in overall tumor size; then, a classical Gompertzian pattern is observed for the MTS radius variation with time. A core with necrotic cells appears for tumor sizes larger than 150 µm, surrounded by a shell of viable tumor cells whose thickness stays almost constant with time. A formula to estimate the size of the necrotic core is proposed. In the second case, the MTS is confined within a healthy tissue. The growth rate is reduced, as compared to the first case - mostly due to the relative adhesion of the tumor and healthy cells to the ECM, and the less favourable transport of nutrients. In particular, for tumor cells adhering less avidly to the ECM, the healthy tissue is progressively displaced as the malignant mass grows, whereas tumor cell infiltration is predicted for the opposite condition. Interestingly, the infiltration potential of the tumor mass is mostly driven by the relative cell adhesion to the ECM. In the third case, a tumor cord model is analyzed where the malignant cells grow around microvessels in a 3D geometry. It is shown that tumor cells tend to migrate among adjacent vessels seeking new oxygen and nutrient. This model can predict and optimize the efficacy of anticancer therapeutic strategies. It can be further developed to answer questions on tumor biophysics, related to the effects of ECM stiffness and cell adhesion on tumor cell proliferation.

Tumor growth modeling from the perspective of multiphase porous media mechanics.

Multiphase porous media mechanics is used for modeling tumor growth, using governing equations obtained via the thermodynamically constrained averaging theory (TCAT). This approach incorporates the interaction of more phases than legacy tumor growth models. The tumor is treated as a multiphase system composed of an extracellular matrix, tumor cells which may become necrotic depending on nutrient level and pressure, healthy cells and an interstitial fluid which transports nutrients. The governing equations are numerically solved within a Finite Element framework for predicting the growth rate of the tumor mass, and of its individual components, as a function of the initial tumor-to-healthy cell ratio, nutrient concentration, and mechanical strain. Preliminary results are shown.