3D mesh processing using GAMer 2 to enable reaction-diffusion simulations in realistic cellular geometries

Abstract

Recent advances in electron microscopy have enabled the imaging of single cells in 3D at nanometer length scale resolutions. An uncharted frontier for in silico biology is the ability to simulate cellular processes using these observed geometries. Enabling such simulations requires watertight meshing of electron micrograph images into 3D volume meshes, which can then form the basis of computer simulations of such processes using numerical techniques such as the Finite Element Method. In this paper, we describe the use of our recently rewritten mesh processing software, GAMer 2, to bridge the gap between poorly conditioned meshes generated from segmented micrographs and boundary marked tetrahedral meshes which are compatible with simulation. We demonstrate the application of a workflow using GAMer 2 to a series of electron micrographs of neuronal dendrite morphology explored at three different length scales and show that the resulting meshes are suitable for finite element simulations. This work is an important step towards making physical simulations of biological processes in realistic geometries routine. Innovations in algorithms to reconstruct and simulate cellular length scale phenomena based on emerging structural data will enable realistic physical models and advance discovery at the interface of geometry and cellular processes. We posit that a new frontier at the intersection of computational technologies and single cell biology is now open.

39 pages, 14 figures. High resolution figures and supplemental movies available upon request

References in corpus (7)

• A roadmap for the computation of persistent homology
• Finite element exterior calculus: from Hodge theory to numerical stability
• Geometric variational crimes: Hilbert complexes, finite element exterior calculus, and problems on hypersurfaces
• Accurate Reaction-Diffusion Operator Splitting on Tetrahedral Meshes for Parallel Stochastic Molecular Simulations
• Semilinear mixed problems on Hilbert complexes and their numerical approximation
• Finite Element Exterior Calculus for Evolution Problems
• Finite Element Exterior Calculus for Parabolic Evolution Problems On Riemannian Hypersurfaces

Cited by in corpus (3)

• Biomembranes undergo complex, non-axisymmetric deformations governed by Kirchhoff-Love kinematics and revealed by a three dimensional computational framework
• Can biophysical models give insight into the synaptic changes associated with addiction?
• Describing the movement of molecules in reduced-dimension models