Collagen-Laponite Nanoclay Hydrogels for Tumor Spheroid Growth.

The extracellular matrix (ECM) plays an important regulatory role in the development and progression of tumoral tissue. Its functions and properties are crucial in determining tumor cell behavior such as invasion, migration, and malignancy development. Our study explores the role of collagen type I in cancer development and spread using engineered tumor models like multicellular spheroids grown in collagen-based hydrogels to simulate early tumor formation. We employ microfluidic techniques to test the hypothesis that (i) adding Laponite nanoclay to collagen hydrogels modifies mechanical and rheological properties and (ii) changing the stiffness of the collagen microenvironment affects tumor spheroid growth. Our findings support our theories and suggest the use of ECM components and engineered tumor models in cancer research, offering a biocompatible and biomimetic method to tailor the mechanical properties of conventional collagen hydrogels.

Force loading explains spatial sensing of ligands by cells.

Cells can sense the density and distribution of extracellular matrix (ECM) molecules by means of individual integrin proteins and larger, integrin-containing adhesion complexes within the cell membrane. This spatial sensing drives cellular activity in a variety of normal and pathological contexts. Previous studies of cells on rigid glass surfaces have shown that spatial sensing of ECM ligands takes place at the nanometre scale, with integrin clustering and subsequent formation of focal adhesions impaired when single integrin-ligand bonds are separated by more than a few tens of nanometres. It has thus been suggested that a crosslinking 'adaptor' protein of this size might connect integrins to the actin cytoskeleton, acting as a molecular ruler that senses ligand spacing directly. Here, we develop gels whose rigidity and nanometre-scale distribution of ECM ligands can be controlled and altered. We find that increasing the spacing between ligands promotes the growth of focal adhesions on low-rigidity substrates, but leads to adhesion collapse on more-rigid substrates. Furthermore, disordering the ligand distribution drastically increases adhesion growth, but reduces the rigidity threshold for adhesion collapse. The growth and collapse of focal adhesions are mirrored by, respectively, the nuclear or cytosolic localization of the transcriptional regulator protein YAP. We explain these findings not through direct sensing of ligand spacing, but by using an expanded computational molecular-clutch model, in which individual integrin-ECM bonds-the molecular clutches-respond to force loading by recruiting extra integrins, up to a maximum value. This generates more clutches, redistributing the overall force among them, and reducing the force loading per clutch. At high rigidity and high ligand spacing, maximum recruitment is reached, preventing further force redistribution and leading to adhesion collapse. Measurements of cellular traction forces and actin flow speeds support our model. Our results provide a general framework for how cells sense spatial and physical information at the nanoscale, precisely tuning the range of conditions at which they form adhesions and activate transcriptional regulation.

A microfluidic-based analysis of 3D macrophage migration after stimulation by Mycobacterium, Salmonella and Escherichia.

Macrophages play an essential role in the process of recognition and containment of microbial infections. These immune cells are recruited to infectious sites to reach and phagocytose pathogens. Specifically, in this article, bacteria from the genus Mycobacterium, Salmonella and Escherichia, were selected to study the directional macrophage movement towards different bacterial fractions. We recreated a three-dimensional environment in a microfluidic device, using a collagen-based hydrogel that simulates the mechanical microarchitecture associated to the Extra Cellular Matrix (ECM). First, we showed that macrophage migration is affected by the collagen concentration of their environment, migrating greater distances at higher velocities with decreasing collagen concentrations. To recreate the infectious microenvironment, macrophages were exposed to lateral gradients of bacterial fractions obtained from the intracellular pathogens M. tuberculosis and S. typhimurium. Our results showed that macrophages migrated directionally, and in a concentration-dependent manner, towards the sites where bacterial fractions are located, suggesting the presence of attractants molecules in all the samples. We confirmed that purified M. tuberculosis antigens, as ESAT-6 and CFP-10, stimulated macrophage recruitment in our device. Finally, we also observed that macrophages migrate towards fractions from non-pathogenic bacteria, such as M. smegmatis and Escherichia coli. In conclusion, our microfluidic device is a useful tool which opens new perspectives to study the recognition of specific antigens by innate immune cells.

A theoretical analysis of the scale separation in a model to predict solid tumour growth.

Solid tumour growth depends on a host of factors which affect the cell life cycle and extracellular matrix vascularization that leads to a favourable environment. The whole solid tumour can either grow or wither in response to the action of the immune system and therapeutics. A personalised mathematical model of such behaviour must consider both the intra- and inter-cellular dynamics and the mechanics of the solid tumour and its microenvironment. However, such wide range of spatial and temporal scales can hardly be modelled in a single model, and require the so-called multiscale models, defined as orchestrations of single-scale component models, connected by relation models that transform the data for one scale to another. While multiscale models are becoming common, there is a well-established engineering approach to the definition of the scale separation, e.g., how the spatiotemporal continuum is split in the various component models. In most studies scale separation is defined as natural, linked to anatomical concepts such as organ, tissue, or cell; but these do not provide reliable definition of scales: for examples skeletal organs can be as large as 500 mm (femur), or as small as 3 mm (stapes). Here we apply a recently proposed scale-separation approach based on the actual experimental and computational limitations to a patient-specific model of the growth of neuroblastoma. The resulting multiscale model can be properly informed with the available experimental data and solved in a reasonable timeframe with the available computational resources.

Extracellular matrix density regulates the formation of tumour spheroids through cell migration.

In this work, we show how the mechanical properties of the cellular microenvironment modulate the growth of tumour spheroids. Based on the composition of the extracellular matrix, its stiffness and architecture can significantly vary, subsequently influencing cell movement and tumour growth. However, it is still unclear exactly how both of these processes are regulated by the matrix composition. Here, we present a centre-based computational model that describes how collagen density, which modulates the steric hindrance properties of the matrix, governs individual cell migration and, consequently, leads to the formation of multicellular clusters of varying size. The model was calibrated using previously published experimental data, replicating a set of experiments in which cells were seeded in collagen matrices of different collagen densities, hence producing distinct mechanical properties. At an initial stage, we tracked individual cell trajectories and speeds. Subsequently, the formation of multicellular clusters was also analysed by quantifying their size. Overall, the results showed that our model could accurately replicate what was previously seen experimentally. Specifically, we showed that cells seeded in matrices with low collagen density tended to migrate more. Accordingly, cells strayed away from their original cluster and thus promoted the formation of small structures. In contrast, we also showed that high collagen densities hindered cell migration and produced multicellular clusters with increased volume. In conclusion, this model not only establishes a relation between matrix density and individual cell migration but also showcases how migration, or its inhibition, modulates tumour growth.

Balance of mechanical forces drives endothelial gap formation and may facilitate cancer and immune-cell extravasation.

The formation of gaps in the endothelium is a crucial process underlying both cancer and immune cell extravasation, contributing to the functioning of the immune system during infection, the unfavorable development of chronic inflammation and tumor metastasis. Here, we present a stochastic-mechanical multiscale model of an endothelial cell monolayer and show that the dynamic nature of the endothelium leads to spontaneous gap formation, even without intervention from the transmigrating cells. These gaps preferentially appear at the vertices between three endothelial cells, as opposed to the border between two cells. We quantify the frequency and lifetime of these gaps, and validate our predictions experimentally. Interestingly, we find experimentally that cancer cells also preferentially extravasate at vertices, even when they first arrest on borders. This suggests that extravasating cells, rather than initially signaling to the endothelium, might exploit the autonomously forming gaps in the endothelium to initiate transmigration.

Image-based Characterization of 3D Collagen Networks and the Effect of Embedded Cells.

Collagen microstructure is closely related to the mechanical properties of tissues and affects cell migration through the extracellular matrix. To study these structures, three-dimensional (3D) in vitro collagen-based gels are often used, attempting to mimic the natural environment of cells. Some key parameters of the microstructure of these gels are fiber orientation, fiber length, or pore size, which define the mechanical properties of the network and therefore condition cell behavior. In the present study, an automated tool to reconstruct 3D collagen networks is used to extract the aforementioned parameters of gels of different collagen concentration and determine how their microstructure is affected by the presence of cells. Two different experiments are presented to test the functionality of the method: first, collagen gels are embedded within a microfluidic device and collagen fibers are imaged by using confocal fluorescence microscopy; second, collagen gels are directly polymerized in a cell culture dish and collagen fibers are imaged by confocal reflection microscopy. Finally, we investigate and compare the collagen microstructure far from and in the vicinities of MDA-MB 23 cells, finding that cell activity during migration was able to strongly modify the orientation of the collagen fibers and the porosity-related values.

The role of nuclear mechanics in cell deformation under creeping flows.

Despite the relevant regulatory role that nuclear deformation plays in cell behaviour, a thorough understanding of how fluid flow modulates the deformation of the cell nucleus in non-confined environments is lacking. In this work, we investigated the dynamics of cell deformation under different creeping flows as a general simulation tool for predicting nuclear stresses and strains. Using this solid-fluid modelling interaction framework, we assessed the stress and strain levels that the cell nucleus experiences as a function of different microenvironmental conditions, such as physical constraints, fluid flows, cytosol properties, and nucleus properties and size. Therefore, the simulation methodology proposed here allows the design of deformability-based experiments involving fluid flow, such as real-time deformability cytometry and dynamic cell culture in bioreactors or microfluidic devices.

A coupled mechano-biochemical model for bone adaptation.

Bone remodelling is a fundamental biological process that controls bone microrepair, adaptation to environmental loads and calcium regulation among other important processes. It is not surprising that bone remodelling has been subject of intensive both experimental and theoretical research. In particular, many mathematical models have been developed in the last decades focusing in particular aspects of this complicated phenomenon where mechanics, biochemistry and cell processes strongly interact. In this paper, we present a new model that combines most of these essential aspects in bone remodelling with especial focus on the effect of the mechanical environment into the biochemical control of bone adaptation mainly associated to the well known RANKL-RANK-OPG pathway. The predicted results show a good correspondence with experimental and clinical findings. For example, our results indicate that trabecular bone is more severely affected both in disuse and disease than cortical bone what has been observed in osteoporotic bones. In future, the methodology proposed would help to new therapeutic strategies following the evolution of bone tissue distribution in osteoporotic patients.

An agent-based method to estimate 3D cell migration trajectories from 2D measurements: Quantifying and comparing T vs CAR-T 3D cell migration.

BACKGROUND AND OBJECTIVE: Immune cell migration is one of the key features that enable immune cells to find invading pathogens, control tissue damage, and eliminate primary developing tumors. Chimeric antigen receptor (CAR) T-cell therapy is a novel strategy in the battle against various cancers. It has been successful in treating hematological tumors, yet it still faces many challenges in the case of solid tumors. In this work, we evaluate the three-dimensional (3D) migration capacity of T and CAR-T cells within dense collagen-based hydrogels. Quantifying three-dimensional (3D) cell migration requires microscopy techniques that may not be readily accessible. Thus, we introduce a straightforward mathematical model designed to infer 3D trajectories of cells from two-dimensional (2D) cell trajectories. METHODS: We develop a 3D agent-based model (ABM) that simulates the temporal changes in the direction of migration with an inverse transform sampling method. Then, we propose an optimization procedure to accurately orient cell migration over time to reproduce cell migration from 2D experimental cell trajectories. With this model, we simulate cell migration assays of T and CAR-T cells in microfluidic devices conducted under hydrogels with different concentrations of type I collagen and validate our 3D cell migration predictions with light-sheet microscopy. RESULTS: Our findings indicate that CAR-T cell migration is more sensitive to collagen concentration increases than T cells, resulting in a more pronounced reduction in their invasiveness. Moreover, our computational model reveals significant differences in 3D movement patterns between T and CAR-T cells. T cells exhibit migratory behavior in 3D whereas that CAR-T cells predominantly move within the XY plane, with limited movement in the Z direction. However, upon the introduction of a CXCL12 chemical gradient, CAR-T cells present migration patterns that closely resemble those of T cells. CONCLUSIONS: This framework demonstrates that 2D projections of 3D trajectories may not accurately represent real migration patterns. Moreover, it offers a tool to estimate 3D migration patterns from 2D experimental data, which can be easily obtained with automatic quantification algorithms. This approach helps reduce the need for sophisticated and expensive microscopy equipment required in laboratories, as well as the computational burden involved in producing and analyzing 3D experimental data.

Physico-chemical characterization of the tumour microenvironment of pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive lethal malignancy that accounts for more than 90% of pancreatic cancer diagnoses. Our research is focused on the physico-chemical properties of the tumour microenvironment (TME), including its tumoural extracellular matrix (tECM), as they may have an important impact on the success of cancer therapies. PDAC xenografts and their decellularized tECM offer a great material source for research in terms of biomimicry with the original human tumour. Our aim was to evaluate and quantify the physico-chemical properties of the PDAC TME. Both cellularized (native TME) and decellularized (tECM) patient-derived PDAC xenografts were analyzed. A factorial design of experiments identified an optimal combination of factors for effective xenograft decellularization. Our results provide a complete advance in our understanding of the PDAC TME and its corresponding stroma, showing that it presents an interconnected porous architecture with very low permeability and small pores due to the contractility of the cellular components. This fact provides a potential therapeutic strategy based on the therapeutic agent size.

Culture of human bone marrow-derived mesenchymal stem cells on of poly(L-lactic acid) scaffolds: potential application for the tissue engineering of cartilage.

PURPOSE: Due to the attractive properties of poly(L-lactic acid) (PLLA) for tissue engineering, the aim was to determine the growth and differentiation capacity of mesenchymal stromal cells (MSCs) in PLLA scaffolds and their potential use in the treatment of cartilage diseases. METHODS: MSCs were cultured in PLLA films and thin porous membranes to study adherence and proliferation. Permeability and porosity were determined for the different scaffolds employed. The optimal conditions for cell seeding were first determined, as well as cell density and distribution inside the PLLA. Scaffolds were then maintained in expansion or chondrogenic differentiation media for 21 days. Apoptosis, proliferation and chondrogenic differentiation was assessed after 21 days in culture by immunohistochemistry. Mechanical characteristics of scaffolds were determined before and after cell seeding. RESULTS: MSCs uniformly adhered to PLLA films as well as to porous membranes. Proliferation was detected only in monolayers of pure PLLA, but was no longer detected after 10 days. Mechanical characterization of PLLA scaffolds showed differences in the apparent compression elastic modulus for the two sizes used. After determining high efficiencies of seeding, the production of extracellular matrix (ECM) was determined and contained aggrecan and collagens type I and X. ECM produced by the cells induced a twofold increase in the apparent elastic modulus of the composite. CONCLUSIONS: Biocompatible PLLA scaffolds have been developed that can be efficiently loaded with MSCs. The scaffold supports chondrogenic differentiation and ECM deposition that improves the mechanics of the scaffold. Although this improvement does not met the expectations of a hyaline-like cartilage ECM, in part due to the lack of a mechanical stimulation, their potential use in the treatment of cartilage pathologies encourages to improve the mechanical component.

Hybrid computational models of multicellular tumour growth considering glucose metabolism.

Cancer cells metabolize glucose through metabolic pathways that differ from those used by healthy and differentiated cells. In particular, tumours have been shown to consume more glucose than their healthy counterparts and to use anaerobic metabolic pathways, even under aerobic conditions. Nevertheless, scientists have still not been able to explain why cancer cells evolved to present an altered metabolism and what evolutionary advantage this might provide them. Experimental and computational models have been increasingly used in recent years to understand some of these biological questions. Multicellular tumour spheroids are effective experimental models as they replicate the initial stages of avascular solid tumour growth. Furthermore, these experiments generate data which can be used to calibrate and validate computational studies that aim to simulate tumour growth. Hybrid models are of particular relevance in this field of research because they model cells as individual agents while also incorporating continuum representations of the substances present in the surrounding microenvironment that may participate in intracellular metabolic networks as concentration or density distributions. Henceforth, in this review, we explore the potential of computational modelling to reveal the role of metabolic reprogramming in tumour growth.

Computational modelling of the mechanical behaviour of protein-based hydrogels.

Protein-based hydrogels have been extensively studied in the field of biomaterials given their ability to mimic living tissues and their special resemblance to the extracellular matrix. Despite this, the methods used for the control of mechanical properties of hydrogels are very limited, focusing mainly on their elasticity, with an often unrealistic characterization of mechanical properties such as extensibility, stiffness and viscoelasticity. Being able to control these properties is essential for the development of new biomaterials, since it has been demonstrated that mechanical properties affect cell behaviour and biological processes. To better understand the mechanical behaviour of these biopolymers, a computational model is here developed to characterize the mechanical behaviour of two different protein-based hydrogels. Strain-stress tests and stress-relaxation tests are evaluated computationally and compared to the results obtained experimentally in a previous work. To achieve this goal the Finite Element Method is used, combining hyperelastic and viscoelastic models. Different hyperelastic constitutive models (Mooney-Rivlin, Neo-Hookean, first and third order Ogden, and Yeoh) are proposed to estimate the mechanical properties of the protein-based hydrogels by least-square fitting of the in-vitro uniaxial test results. Among these models, the first order Ogden model with a viscoelastic model defined in Prony parameters better reproduces the strain-stress response and the change of stiffness with strain observed in the in-vitro tests.

Fluid flow to mimic organ function in 3D in vitro models.

Many different strategies can be found in the literature to model organ physiology, tissue functionality, and disease in vitro; however, most of these models lack the physiological fluid dynamics present in vivo. Here, we highlight the importance of fluid flow for tissue homeostasis, specifically in vessels, other lumen structures, and interstitium, to point out the need of perfusion in current 3D in vitro models. Importantly, the advantages and limitations of the different current experimental fluid-flow setups are discussed. Finally, we shed light on current challenges and future focus of fluid flow models applied to the newest bioengineering state-of-the-art platforms, such as organoids and organ-on-a-chip, as the most sophisticated and physiological preclinical platforms.

Effects of tumour heterogeneous properties on modelling the transport of radiative particles.

Dose calculation plays a critical role in radiotherapy (RT) treatment planning, and there is a growing need to develop accurate dose deposition models that incorporate heterogeneous tumour properties. Deterministic models have demonstrated their capability in this regard, making them the focus of recent treatment planning studies as they serve as a basis for simplified models in RT treatment planning. In this study, we present a simplified deterministic model for photon transport based on the Boltzmann transport equation (BTE) as a proof-of-concept to illustrate the impact of heterogeneous tumour properties on RT treatment planning. We employ the finite element method (FEM) to simulate the photon flux and dose deposition in real cases of diffuse intrinsic pontine glioma (DIPG) and neuroblastoma (NB) tumours. Importantly, in light of the availability of pipelines capable of extracting tumour properties from magnetic resonance imaging (MRI) data, we highlight the significance of such data. Specifically, we utilise cellularity data extracted from DIPG and NB MRI images to demonstrate the importance of heterogeneity in dose calculation. Our model simplifies the process of simulating a RT treatment system and can serve as a useful starting point for further research. To simulate a full RT treatment system, one would need a comprehensive model that couples the transport of electrons and photons.

PhysiCOOL: A generalized framework for model Calibration and Optimization Of modeLing projects.

In silico models of biological systems are usually very complex and rely on a large number of parameters describing physical and biological properties that require validation. As such, parameter space exploration is an essential component of computational model development to fully characterize and validate simulation results. Experimental data may also be used to constrain parameter space (or enable model calibration) to enhance the biological relevance of model parameters. One widely used computational platform in the mathematical biology community is PhysiCell, which provides a standardized approach to agent-based models of biological phenomena at different time and spatial scales. Nonetheless, one limitation of PhysiCell is the lack of a generalized approach for parameter space exploration and calibration that can be run without high-performance computing access. Here, we present PhysiCOOL, an open-source Python library tailored to create standardized calibration and optimization routines for PhysiCell models.

A hybrid physics-based and data-driven framework for cellular biological systems: Application to the morphogenesis of organoids.

How cells orchestrate their cellular functions remains a crucial question to unravel how they organize in different patterns. We present a framework based on artificial intelligence to advance the understanding of how cell functions are coordinated spatially and temporally in biological systems. It consists of a hybrid physics-based model that integrates both mechanical interactions and cell functions with a data-driven model that regulates the cellular decision-making process through a deep learning algorithm trained on image data metrics. To illustrate our approach, we used data from 3D cultures of murine pancreatic ductal adenocarcinoma cells (PDAC) grown in Matrigel as tumor organoids. Our approach allowed us to find the underlying principles through which cells activate different cell processes to self-organize in different patterns according to the specific microenvironmental conditions. The framework proposed here expands the tools for simulating biological systems at the cellular level, providing a novel perspective to unravel morphogenetic patterns.

Hydrocolloids of Egg White and Gelatin as a Platform for Hydrogel-Based Tissue Engineering.

Innovative materials are needed to produce scaffolds for various tissue engineering and regenerative medicine (TERM) applications, including tissue models. Materials derived from natural sources that offer low production costs, easy availability, and high bioactivity are highly preferred. Chicken egg white (EW) is an overlooked protein-based material. Whilst its combination with the biopolymer gelatin has been investigated in the food technology industry, mixed hydrocolloids of EW and gelatin have not been reported in TERM. This paper investigates these hydrocolloids as a suitable platform for hydrogel-based tissue engineering, including 2D coating films, miniaturized 3D hydrogels in microfluidic devices, and 3D hydrogel scaffolds. Rheological assessment of the hydrocolloid solutions suggested that temperature and EW concentration can be used to fine-tune the viscosity of the ensuing gels. Fabricated thin 2D hydrocolloid films presented globular nano-topography and in vitro cell work showed that the mixed hydrocolloids had increased cell growth compared with EW films. Results showed that hydrocolloids of EW and gelatin can be used for creating a 3D hydrogel environment for cell studies inside microfluidic devices. Finally, 3D hydrogel scaffolds were fabricated by sequential temperature-dependent gelation followed by chemical cross-linking of the polymeric network of the hydrogel for added mechanical strength and stability. These 3D hydrogel scaffolds displayed pores, lamellae, globular nano-topography, tunable mechanical properties, high affinity for water, and cell proliferation and penetration properties. In conclusion, the large range of properties and characteristics of these materials provide a strong potential for a large variety of TERM applications, including cancer models, organoid growth, compatibility with bioprinting, or implantable devices.

Tumour growth: An approach to calibrate parameters of a multiphase porous media model based on in vitro observations of Neuroblastoma spheroid growth in a hydrogel microenvironment.

To unravel processes that lead to the growth of solid tumours, it is necessary to link knowledge of cancer biology with the physical properties of the tumour and its interaction with the surrounding microenvironment. Our understanding of the underlying mechanisms is however still imprecise. We therefore developed computational physics-based models, which incorporate the interaction of the tumour with its surroundings based on the theory of porous media. However, the experimental validation of such models represents a challenge to its clinical use as a prognostic tool. This study combines a physics-based model with in vitro experiments based on microfluidic devices used to mimic a three-dimensional tumour microenvironment. By conducting a global sensitivity analysis, we identify the most influential input parameters and infer their posterior distribution based on Bayesian calibration. The resulting probability density is in agreement with the scattering of the experimental data and thus validates the proposed workflow. This study demonstrates the huge challenges associated with determining precise parameters with usually only limited data for such complex processes and models, but also demonstrates in general how to indirectly characterise the mechanical properties of neuroblastoma spheroids that cannot feasibly be measured experimentally.