Two-dimensional finite element model to study calcium distribution in astrocytes in presence of excess buffer

Abstract

In this paper a finite element model is developed to study cytosolic calcium concentration distribution in astrocytes for a two-dimensional steady-state case in presence of excess buffer. The mathematical model of calcium diffusion in astrocytes leads to a boundary value problem involving elliptical partial differential equation. The model consists of reaction–diffusion phenomena, association and dissociation rates and buffer. A point source of calcium is incorporated in the model. Appropriate boundary conditions have been framed. Finite element method is employed to solve the problem. A MATLAB program has been developed for the entire problem and simulated to compute the numerical results. The numerical results have been used to plot calcium concentration profiles in astrocytes. The effect of EGTA, BAPTA and σCa influx on calcium concentration distribution in astrocytes is studied with the help of numerical results.