Diffusion on a curved surface coupled to diffusion in the volume: Application to cell biology

  • Authors:
  • Igor L. Novak;Fei Gao;Yung-Sze Choi;Diana Resasco;James C. Schaff;Boris M. Slepchenko

  • Affiliations:
  • Center for Cell Analysis and Modeling, Department of Cell Biology, University of Connecticut Health Center, Farmington, CT 06030, United States;Center for Cell Analysis and Modeling, Department of Cell Biology, University of Connecticut Health Center, Farmington, CT 06030, United States;Department of Mathematics, University of Connecticut, Storrs, CT 06269, United States;Department of Computer Science, Yale University, New Haven, CT 06520, United States;Center for Cell Analysis and Modeling, Department of Cell Biology, University of Connecticut Health Center, Farmington, CT 06030, United States;Center for Cell Analysis and Modeling, Department of Cell Biology, University of Connecticut Health Center, Farmington, CT 06030, United States

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2007

Quantified Score

Hi-index 31.46

Visualization

Abstract

An algorithm is presented for solving a diffusion equation on a curved surface coupled to diffusion in the volume, a problem often arising in cell biology. It applies to pixilated surfaces obtained from experimental images and performs at low computational cost. In the method, the Laplace-Beltrami operator is approximated locally by the Laplacian on the tangential plane and then a finite volume discretization scheme based on a Voronoi decomposition is applied. Convergence studies show that mass conservation built in the discretization scheme and cancellation of sampling error ensure convergence of the solution in space with an order between 1 and 2. The method is applied to a cell-biological problem where a signaling molecule, G-protein Rac, cycles between the cytoplasm and cell membrane thus coupling its diffusion in the membrane to that in the cell interior. Simulations on realistic cell geometry are performed to validate, and determine the accuracy of, a recently proposed simplified quantitative analysis of fluorescence loss in photobleaching. The method is implemented within the Virtual Cell computational framework freely accessible at www.vcell.org.