FEM-based oxygen consumption and cell viability models for avascular pancreatic islets

<p>Abstract</p> <p>Background</p> <p>The function and viability of cultured, transplanted, or encapsulated pancreatic islets is often limited by hypoxia because these islets have lost their vasculature during the isolation process and have to rely on gradient-driven passive diffusion, which cannot provide adequate oxygen transport. Pancreatic islets (islets of Langerhans) are particularly susceptible due to their relatively large size, large metabolic demand, and increased sensitivity to hypoxia. Here, finite element method (FEM) based multiphysics models are explored to describe oxygen transport and cell viability in avascular islets both in static and in moving culture media.</p> <p>Methods</p> <p>Two- and three-dimensional models were built in COMSOL Multiphysics using the convection and diffusion as well as the incompressible Navier-Stokes fluid dynamics application modes. Oxygen consumption was assumed to follow Michaelis-Menten-type kinetics and to cease when local concentrations fell below a critical threshold; in a dynamic model, it was also allowed to increase with increasing glucose concentration.</p> <p>Results</p> <p>Partial differential equation (PDE) based exploratory cellular-level oxygen consumption and cell viability models incorporating physiologically realistic assumptions have been implemented for fully scaled cell culture geometries with 100, 150, and 200 <it>μ</it>m diameter islets as representative. Calculated oxygen concentrations and intra-islet regions likely to suffer from hypoxia-related necrosis obtained for traditional flask-type cultures, oxygen-permeable silicone-rubber membrane bottom cultures, and perifusion chambers with flowing media and varying incoming glucose levels are presented in detail illustrated with corresponding colour-coded figures and animations.</p> <p>Conclusion</p> <p>Results of the computational models are, as a first estimate, in good quantitative agreement with existing experimental evidence, and they confirm that during culture, hypoxia is often a problem for non-vascularised islet and can lead to considerable cell death (necrosis), especially in the core region of larger islets. Such models are of considerable interest to improve the function and viability of cultured, transplanted, or encapsulated islets. The present implementation allows convenient extension to true multiphysics applications that solve coupled physics phenomena such as diffusion and consumption with convection due to flowing or moving media.</p

A local glucose-and oxygen concentration-based insulin secretion model for pancreatic islets

<p>Abstract</p> <p>Background</p> <p>Because insulin is the main regulator of glucose homeostasis, quantitative models describing the dynamics of glucose-induced insulin secretion are of obvious interest. Here, a computational model is introduced that focuses not on organism-level concentrations, but on the quantitative modeling of local, cellular-level glucose-insulin dynamics by incorporating the detailed spatial distribution of the concentrations of interest within isolated avascular pancreatic islets.</p> <p>Methods</p> <p>All nutrient consumption and hormone release rates were assumed to follow Hill-type sigmoid dependences on local concentrations. Insulin secretion rates depend on both the glucose concentration and its time-gradient, resulting in second-and first-phase responses, respectively. Since hypoxia may also be an important limiting factor in avascular islets, oxygen and cell viability considerations were also built in by incorporating and extending our previous islet cell oxygen consumption model. A finite element method (FEM) framework is used to combine reactive rates with mass transport by convection and diffusion as well as fluid-mechanics.</p> <p>Results</p> <p>The model was calibrated using experimental results from dynamic glucose-stimulated insulin release (GSIR) perifusion studies with isolated islets. Further optimization is still needed, but calculated insulin responses to stepwise increments in the incoming glucose concentration are in good agreement with existing experimental insulin release data characterizing glucose and oxygen dependence. The model makes possible the detailed description of the intraislet spatial distributions of insulin, glucose, and oxygen levels. In agreement with recent observations, modeling also suggests that smaller islets perform better when transplanted and/or encapsulated.</p> <p>Conclusions</p> <p>An insulin secretion model was implemented by coupling local consumption and release rates to calculations of the spatial distributions of all species of interest. The resulting glucose-insulin control system fits in the general framework of a sigmoid proportional-integral-derivative controller, a generalized PID controller, more suitable for biological systems, which are always nonlinear due to the maximum response being limited. Because of the general framework of the implementation, simulations can be carried out for arbitrary geometries including cultured, perifused, transplanted, and encapsulated islets.</p

The influence of inertia and contact angle on the instability of partially wetting liquid strips : a numerical analysis study

The stability of a thread of fluid deposited on a flat solid substrate is studied numerically by means of the Finite Element Method in combination with an Arbitrary Lagrangian-Eulerian technique. A good agreement is observed when our results are compared with predictions of linear stability analysis obtained by other authors. Moreover, we also analysed the influence of inertia for different contact angles and found that inertia strongly affects the growth rate of the instability when contact angles are large. By contrast, the wave number of the fastest growing mode does not show important variations with inertia. The numerical technique allows us to follow the evolution of the free surface instability until comparatively late stages, where the filament begins to break into droplets. The rupture pattern observed for several cases shows that the number of principal droplets agrees reasonably well with an estimation based on the fastest growing modes