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.

Mechanical principles governing the shapes of dendritic spines

“…We do not explicitly consider the remodeling of the actin network or the dynamics of the associated proteins, but use force as a lumped parameter. Additionally, the use of axisymmetric coordinates restricts our ability to obtain realistic spine shapes [124]. We note that the technology required to simulate these processes is quite challenging [125] and is ongoing research in our group.…”

Deciphering the postsynaptic calcium-mediated energy homeostasis through mitochondria-endoplasmic reticulum contact sites using systems modeling

“…First, all our geometries are simplified, allowing us to construct a spatial model of multiple organelles. We note that current efforts include using realistic geometries of spines and organelles derived from microscopy to investigate how these play a role in Ca 2+ and ATP dynamics [52,58,59]. Second, we note that we assume that Ca 2+ and ATP are present in large enough amounts to justify a deterministic approach.…”

Transfer function for YAP/TAZ nuclear translocation revealed through spatial systems modeling

“…Cell-cell and cell matrix adhesions can have synergistic and antagonistic effects (4) , and regulation of YAP/TAZ by cell-cell communication through LATs and Merlin are additional key features to explore in future work. Incorporating cell type specific protein expression levels, realistic cellular geometries and physiologically relevant activation areas such as number, size, and localization of focal adhesions (83)(84)(85)(86)(87)(88) will also offer further insight into how YAP/TAZ translocation is mediated by substrate stiffness and cell shape. Likewise, the timescale of YAP/TAZ nuclear localization could be important (89) .…”