An official website of the United States government
We present a multi-scale mathematical model and a novel numerical solver to study blood plasma flow and oxygen concentration in a prototype model of an implantable Bioartificial Pancreas (iBAP) that operates under arteriovenous pressure differential without the need for immunosuppressive therapy. The iBAP design consists of a poroelastic cell scaffold containing the healthy transplanted cells, encapsulated between two semi-permeable nano-pore size membranes to prevent the patient’s own immune cells from attacking the transplant. The device is connected to the patient’s vascular system via an anastomosis graft bringing oxygen and nutrients to the transplanted cells of which oxygen is the limiting factor for long-term viability. Mathematically, we propose a (nolinear) fluid–poroelastic structure interaction model to describe the flow of blood plasma through the scaffold containing the cells, and a set of (nonlinear) advection–reaction–diffusion equations defined on moving domains to study oxygen supply to the cells. These macro-scale models are solved using finite element method based solvers. One of the novelties of this work is the design of a novel second-order accurate fluid–poroelastic structure interaction solver, for which we prove that it is unconditionally stable. At the micro/nano-scale, Smoothed Particle Hydrodynamics (SPH) simulations are used to capture the micro/nano-structure (architecture) of cell scaffolds and obtain macro-scale parameters, such as hydraulic conductivity/permeability, from the micro-scale scaffold-specific architecture. To avoid expensive micro-scale simulations based on SPH simulations for every new scaffold architecture, we use Encoder–Decoder Convolution Neural Networks. Based on our numerical simulations, we propose improvements in the current prototype design. For example, we show that highly elastic scaffolds have a higher capacity for oxygen transfer, which is an important finding considering that scaffold elasticity can be controlled during their fabrication, and that elastic scaffolds improve cell viability. The mathematical and computational approaches developed in this work provide a benchmark tool for computational analysis of not only iBAP, but also, more generally, of cell encapsulation strategies used in the design of devices for cell therapy and bio-artificial organs.
more »
« less
This content will become publicly available on November 11, 2025
A computational algorithm for optimal design of a bioartificial organ scaffold architecture
We develop a computational algorithm based on a diffuse interface approach to study the design of bioartificial organ scaffold architectures. These scaffolds, composed of poroelastic hydrogels housing transplanted cells, are linked to the patient’s blood circulation via an anastomosis graft. Before entering the scaffold, the blood flow passes through a filter, and the resulting filtered blood plasma transports oxygen and nutrients to sustain the viability of transplanted cells over the long term. A key issue in maintaining cell viability is the design of ultrafiltrate channels within the hydrogel scaffold to facilitate advection-enhanced oxygen supply ensuring oxygen levels remain above a critical threshold to prevent hypoxia. In this manuscript, we develop a computational algorithm to analyze the plasma flow and oxygen concentration within hydrogels featuring various channel geometries. Our objective is to identify the optimal hydrogel channel architecture that sustains oxygen concentration throughout the scaffold above the critical hypoxic threshold. The computational algorithm we introduce here employs a diffuse interface approach to solve a multi-physics problem. The corresponding model couples the time-dependent Stokes equations, governing blood plasma flow through the channel network, with the time-dependent Biot equations, characterizing Darcy velocity, pressure, and displacement within the poroelastic hydrogel containing the transplanted cells. Subsequently, the calculated plasma velocity is utilized to determine oxygen concentration within the scaffold using a diffuse interface advection-reaction-diffusion model. Our investigation yields a scaffold architecture featuring a hexagonal network geometry that meets the desired oxygen concentration criteria. Unlike classical sharp interface approaches, the diffuse interface approach we employ is particularly adept at addressing problems with intricate interface geometries, such as those encountered in bioartificial organ scaffold design. This study is significant because recent developments in hydrogel fabrication make it now possible to control hydrogel rheology and utilize computational results to generate optimized scaffold architectures.
more »
« less
- PAR ID:
- 10583554
- Editor(s):
- Haugh, Jason M
- Publisher / Repository:
- PLoS
- Date Published:
- Journal Name:
- PLOS Computational Biology
- Volume:
- 20
- Issue:
- 11
- ISSN:
- 1553-7358
- Page Range / eLocation ID:
- e1012079
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Thanks to their softness, biocompatibility, porosity, and ready availability, hydrogels are commonly used in microfluidic assays and organ-on-chip devices as a matrix for cells. They not only provide a supporting scaffold for the differentiating cells and the developing organoids, but also serve as the medium for transmitting oxygen, nutrients, various chemical factors, and mechanical stimuli to the cells. From a bioengineering viewpoint, the transmission of forces from fluid perfusion to the cells through the hydrogel is critical to the proper function and development of the cell colony. In this paper, we develop a poroelastic model to represent the fluid flow through a hydrogel containing a biological cell modeled as a hyperelastic inclusion. In geometries representing shear and normal flows that occur frequently in microfluidic experiments, we use finite-element simulations to examine how the perfusion engenders interstitial flow in the gel and displaces and deforms the embedded cell. The results show that pressure is the most important stress component in moving and deforming the cell, and the model predicts the velocity in the gel and stress transmitted to the cell that is comparable to in vitro and in vivo data. This work provides a computational tool to design the geometry and flow conditions to achieve optimal flow and stress fields inside the hydrogels and around the cell.more » « less
-
Abstract Inadequate oxygenation is a major challenge in cell encapsulation, a therapy which holds potential to treat many diseases including type I diabetes. In such systems, cellular oxygen (O 2 ) delivery is limited to slow passive diffusion from transplantation sites through the poorly O 2 -soluble encapsulating matrix, usually a hydrogel. This constrains the maximum permitted distance between the encapsulated cells and host site to within a few hundred micrometers to ensure cellular function. Inspired by the natural gas-phase tracheal O 2 delivery system of insects, we present herein the design of a biomimetic scaffold featuring internal continuous air channels endowed with 10,000-fold higher O 2 diffusivity than hydrogels. We incorporate the scaffold into a bulk hydrogel containing cells, which facilitates rapid O 2 transport through the whole system to cells several millimeters away from the device-host boundary. A computational model, validated by in vitro analysis, predicts that cells and islets maintain high viability even in a thick (6.6 mm) device. Finally, the therapeutic potential of the device is demonstrated through the correction of diabetes in immunocompetent mice using rat islets for over 6 months.more » « less
-
Abstract 3D bioprinting improves orientation ofin vitrotumor models by offering layer by layer positioning of cancer cells and cancer associated fibroblasts (CAFs) which can replicate tumor microenvironment. Aim of this study was to develop a sodium alginate -gelatin (SA-GL) hydrogel by optimizing rheological parameters to print non-small cell lung cancer (NSCLC) patient derived xenograft (PDX) cells and lung CAFs co-cultures. SA-GL hydrogels were prepared, and rheological properties were evaluated. Both the cells were mixed with the hydrogel and printed using INKREDIBLE bioprinter. Hydrogels prepared with 3.25% and 3.5% (w/v) SA and 4% (w/v) GL showed higher printability and cell viability. A significant decline in viscosity with shear rate was observed in these hydrogels suggesting the shear thinning property of hydrogels. Spheroid size distribution after 15 days was in the diameter range of 50–1100 µm. Up-regulation of vimentin, α-SMA and loss of E-cadherin in co-culture spheroids confirmed cellular crosstalk. This study demonstrates that rheological optimization of SA-GL hydrogel enhances printability and viability of NSCLC PDX and CAF co-culture which allows 3D co-culture spheroid formation within the printed scaffold. Therefore, this model can be used for studying high throughput drug screening and other pre-clinical applications.more » « less
-
Abstract The human placenta is a critical organ, mediating the exchange of nutrients, oxygen, and waste products between fetus and mother. Placental malaria (PM) resulted from Plasmodium falciparum infections causes up to 200 thousand newborn deaths annually, mainly due to low birth weight, as well as 10 thousand mother deaths. In this work, a placenta-on-a-chip model is developed to mimic the nutrient exchange between the fetus and mother under the influence of PM. In this model, trophoblasts cells (facing infected or uninfected blood simulating maternal blood and termed “trophoblast side”) and human umbilical vein endothelial cells (facing uninfected blood simulating fetal blood and termed “endothelial” side) are cultured on the opposite sides of an extracellular matrix gel in a compartmental microfluidic system, forming a physiological barrier between the co-flow tubular structure to mimic a simplified maternal–fetal interface in placental villi. The influences of infected erythrocytes (IEs) sequestration through cytoadhesion to chondroitin sulfate A (CSA) expressed on the surface of trophoblast cells, a critical feature of PM, on glucose transfer efficiency across the placental barrier was studied. To create glucose gradients across the barrier, uninfected erythrocyte or IE suspension with a higher glucose concentration was introduced into the “trophoblast side” and a culture medium with lower glucose concentration was introduced into the “endothelial side”. The glucose levels in the endothelial channel in response to CSA-adherent erythrocytes infected with CS2 line of parasites in trophoblast channel under flow conditions was monitored. Uninfected erythrocytes served as a negative control. The results demonstrated that CSA-binding IEs added resistance to the simulated placental barrier for glucose perfusion and decreased the glucose transfer across this barrier. The results of this study can be used for better understanding of PM pathology and development of models useful in studying potential treatment of PM.more » « less
