Hybrid finite element method for describing the electrical response of biological cells to applied fields.
A novel hybrid finite element method (FEM) for modeling the response of passive and active biological membranes to external stimuli is presented. The method is based on the differential equations that describe the conservation of electric flux and membrane currents. By introducing the electric flux through the cell membrane as an additional variable, the algorithm decouples the linear partial differential equation part from the nonlinear ordinary differential equation part that defines the membrane dynamics of interest. This conveniently results in two subproblems: a linear interface problem and a nonlinear initial value problem. The linear interface problem is solved with a hybrid FEM. The initial value problem is integrated by a standard ordinary differential equation solver such as the Euler and Runge-Kutta methods. During time integration, these two subproblems are solved alternatively. The algorithm can be used to model the interaction of stimuli with multiple cells of almost arbitrary geometries and complex ion-channel gating at the plasma membrane. Numerical experiments are presented demonstrating the uses of the method for modeling field stimulation and action potential propagation.
Ying
Wenjun
W
Department of Biomedical Engineering, Duke University, Durham, NC 27708-0281, USA. ying@cel-mail.bme.duke.edu
Henriquez
Craig S
CS
D000818:Animals; D002462:Cell Membrane / Q000502:physiology* / Q000528:radiation effects*; D002468:Cell Physiological Phenomena / Q000528:radiation effects; D003198:Computer Simulation; D004307:Dose-Response Relationship, Radiation; D004558:Electric Stimulation / Q000379:methods*; D004574:Electromagnetic Fields; D020342:Finite Element Analysis; D006801:Humans; D008564:Membrane Potentials / Q000502:physiology* / Q000528:radiation effects*; D008954:Models, Biological; D011829:Radiation Dosage