skip to main content

Explore Research Products in the NSF-PAR

  • NSF Public Access
  • Search Results
  • Mathematical and Computational Modeling of Poroelastic Cell Scaffolds Used in the Design of an Implantable Bioartificial Pancreas
Title: Mathematical and Computational Modeling of Poroelastic Cell Scaffolds Used in the Design of an Implantable Bioartificial Pancreas
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
Award ID(s):
2011319
NSF-PAR ID:
10429493
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Fluids
Volume:
7
Issue:
7
ISSN:
2311-5521
Page Range / eLocation ID:
222
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ( , Current Protocols)
    Abstract

    Stem cell−based therapies for various ailments have attracted significant attention for over a decade. However, low retention of transplanted cells at the damaged site has hindered their potential for use in therapy. Tissue engineered grafts with fibrillar structures mimicking the extracellular matrix (ECM) can be potentially used to increase the retention and engraftment of stem cells at the damaged site. Moreover, these grafts may also provide mechanical stability at the damaged site to enhance function and regeneration. Among all the methods to produce fibrillar structures developed in recent years, electrospinning is a simple and versatile method to produce fibrous structures ranging from a few nanometers to micrometers. Coaxial electrospinning enables production of a mechanically stable core with a cell‐binding sheath for enhanced cell adhesion and proliferation. Furthermore, this process provides an alternative to functionalized engineered scaffolds with specific compositions. The present article describes the protocol for developing a polycaprolactone (PCL) core and gelatin/gelatin methacrylate (GelMA) sheath laden with stem cells for various regenerative engineering applications. © 2021 Wiley Periodicals LLC.

    This article was corrected on 18 July 2022. See the end of the full text for details.

    Basic Protocol 1: Uniaxial PCL electrospinning

    Basic Protocol 2: Coaxial electrospinning

    Support Protocol 1: Scaffold characterization for Basic Protocols 1 and 2

    Basic Protocol 3: Cell seeding on uniaxial and coaxial electrospun scaffolds and MTS assay

    Support Protocol 2: Preparation of scaffold with cells for scanning electron microscopy

     
    more » « less
  2. ; ; ; ( , Proceedings of the Institution of Mechanical Engineers, Part B: Journal of Engineering Manufacture)
    Pore size, external shape, and internal complexity of additively manufactured porous titanium scaffolds are three primary determinants of cell viability and structural strength of scaffolds in bone tissue engineering. To obtain an optimal design with the combination of all three determinants, four scaffolds each with a unique topology (external geometry and internal structure) were designed and varied the pore sizes of each scaffold 3 times. For each topology, scaffolds with pore sizes of 300, 400, and 500 µm were designed. All designed scaffolds were additively manufactured in material Ti6Al4V by the direct metal laser melting machine. Compression test was conducted on the scaffolds to assure meeting minimum compressive strength of human bone. The effects of pore size and topology on the cell viability of the scaffolds were analyzed. The 12 scaffolds were ultrasonically cleaned and seeded with NIH3T3 cells. Each scaffold was seeded with 1 million cells. After 32 days of culturing, the cells were fixed for their three-dimensional architecture preservation and to obtain scanning electron microscope images. 
    more » « less
  3. ( , PLOS computational biology)
    Recent advances in modeling oxygen supply to cortical brain tissue have begun to elucidate the functional mechanisms of neurovascular coupling. While the principal mechanisms of blood flow regulation after neuronal firing are generally known, mechanistic hemodynamic simulations cannot yet pinpoint the exact spatial and temporal coordination between the network of arteries, arterioles, capillaries and veins for the entire brain. Because of the potential significance of blood flow and oxygen supply simulations for illuminating spatiotemporal regulation inside the cortical microanatomy, there is a need to create mathematical models of the entire cerebral circulation with realistic anatomical detail. Our hypothesis is that an anatomically accurate reconstruction of the cerebrocirculatory architecture will inform about possible regulatory mechanisms of the neurovascular interface. In this article, we introduce large-scale networks of the murine cerebral circulation spanning the Circle of Willis, main cerebral arteries connected to the pial network down to the microcirculation in the capillary bed. Several multiscale models were generated from state-of-the-art neuroimaging data. Using a vascular network construction algorithm, the entire circulation of the middle cerebral artery was synthesized. Blood flow simulations indicate a consistent trend of higher hematocrit in deeper cortical layers, while surface layers with shorter vascular path lengths seem to carry comparatively lower red blood cell (RBC) concentrations. Moreover, the variability of RBC flux decreases with cortical depth. These results support the notion that plasma skimming serves a self-regulating function for maintaining uniform oxygen perfusion to neurons irrespective of their location in the blood supply hierarchy. Our computations also demonstrate the practicality of simulating blood flow for large portions of the mouse brain with existing computer resources. The efficient simulation of blood flow throughout the entire middle cerebral artery (MCA) territory is a promising milestone towards the final aim of predicting blood flow patterns for the entire brain. 
    more » « less
  4. ; ; ; ; ( , 2021 Design of Medical Devices Conference)
    Abstract

    Metastatic cancer in bones is incurable, which causes significant mobility and mortality to the patients. In this work, we investigate the role of interstitial fluid flow on cancer cells' growth within the interconnected pores of human bone. In-vitro experiments were carried out in a bio-reactor which includes bone-like scaffold specimens. A pump is used to maintain a laminar flow condition inside the bioreactor to resemble fluid flow in bones. The scaffold specimens are harvested after 23 days in the bioreactor. The scaffold specimen is scanned with Micro-CT under the resolution of 70 micrometers. We created a full-scale 3D computational model of the scaffold based on the micro-CT data using the open-source software Seg3D and Meshmixer. Based on the geometrical models, we generated the computational grids using the commercial software Gridgen. We performed Computational Fluid Dynamics (CFD) simulations with the immersed boundary method (Gilmanov, Le, Sotiropoulos, JCP 300, 1, 2015) to investigate the flow patterns inside the pores of the scaffolds. The results reveal a non-uniform flow distribution in the vicinity of the scaffold. The flow velocity and the shear stress distributions inside the scaffold are shown to be convoluted and very sensitive to the pore sizes. Our future work will further quantify these distributions and correlate them to cancer cells' growth observed in the experiments.

     
    more » « less
  5. ; ; ; ; ( , Proceedings of the 2022 Design of Medical Devices Conference)
    Our study aims to identify the role of fluid flow in the growth of human bone cancer cells during metastasis. In our experiments, the cancer cells are seeded on the surface of cylindrical scaffolds in a bioreactor. The flow is laminar flow, which mimics the physiological conditions of the human body. A full-scale 3D high-resolution computational mesh of scaffold was created based on the physical scaffold's Micro-CT scans using open-source imaging software Slicer3D and Meshmixer. To investigate the influences of the flow on the seeded cells, we performed Computational Fluid Dynamics (CFD) simulations with the immersed boundary method (Gilmanov, Le, Sotiropoulos, JCP 300, 1, 2015). The computational domain was generated using the commercial software Gridgen. Our results show that the fluid flow velocity is highly dependent on the shape and pore sizes. In addition, the magnitude of the velocity on the surface where the cells are seeded is in between [0-0.05] μm/sallowing the cells to grow without being detached from the surface of the scaffold. Our future work will focus on (i) investigating the role of the shear stress on the distribution and orientation of the cancer cells. (ii) Simulating multiple scaffolds within the bioreactor to further quantify the impact of the gap on the flow velocity and shear. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email