We present a biochemical model of the wall shear stress-induced activation of endothelial nitric oxide synthase (eNOS) in an endothelial cell. The model includes three key mechanotransducers: mechanosensing ion channels, integrins, and G protein-coupled receptors. The reaction cascade consists of two interconnected parts. The first is rapid activation of calcium, which results in formation of calcium-calmodulin complexes, followed by recruitment of eNOS from caveolae. The second is phosphorylation of eNOS by protein kinases PKC and AKT. The model also includes a negative feedback loop due to inhibition of calcium influx into the cell by cyclic guanosine monophosphate (cGMP). In this feedback, increased nitric oxide (NO) levels cause an increase in cGMP levels, so that cGMP inhibition of calcium influx can limit NO production. The model was used to predict the dynamics of NO production by an endothelial cell subjected to a step increase of wall shear stress from zero to a finite physiologically relevant value. Among several experimentally observed features, the model predicts a highly nonlinear, biphasic transient behavior of eNOS activation and NO production: a rapid initial activation due to the very rapid influx of calcium into the cytosol (occurring within 1-5 min) is followed by a sustained period of activation due to protein kinases.

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.

Advances in imaging methods such as electron microscopy, tomography, and other modalities are enabling high-resolution reconstructions of cellular and organelle geometries. Such advances pave the way for using these geometries for biophysical and mathematical modeling once these data can be represented as a geometric mesh, which, when carefully conditioned, enables the discretization and solution of partial differential equations. In this work, we outline the steps for a naïve user to approach the Geometry-preserving Adaptive MeshER software version 2, a mesh generation code written in C++ designed to convert structural data sets to realistic geometric meshes while preserving the underlying shapes. We present two example cases: 1) mesh generation at the subcellular scale as informed by electron tomography and 2) meshing a protein with a structure from x-ray crystallography. We further demonstrate that the meshes generated by the Geometry-preserving Adaptive MeshER software are suitable for use with numerical methods. Together, this collection of libraries and tools simplifies the process of constructing realistic geometric meshes from structural biology data.

Recent advances in electron microscopy have, for the first time, enabled imaging of single cells in 3D at a nanometer length scale resolution. An uncharted frontier for in silico biology is the ability to simulate cellular processes using these observed geometries. Enabling such simulations will require a system for going from electron micrographs to 3D volume meshes, which can then form the basis of computer simulations of such processes using numerical techniques such as the Finite Element Method (FEM). In this paper, we develop an end-to-end pipeline for this task by adapting and extending computer graphics mesh processing and smoothing algorithms. Our workflow makes use of our recently rewritten mesh processing software, GAMer 2, which implements several mesh conditioning algorithms and serves as a platform to connect different pipeline steps. We apply this pipeline to a series of electron micrographs of neuronal dendrite morphology explored at three different length scales and demonstrate that the resultant meshes are suitable for finite element simulations. Our pipeline, which consists of free and open-source community driven tools, is a step towards routine physical simulations of biological processes in realistic geometries. We posit that a new frontier at the intersection of computational technologies and single cell biology is now open. Innovations in algorithms to reconstruct and simulate cellular length scale phenomena based on emerging structural data will enable realistic physical models and advance discovery. Author summary3D imaging of cellular components and associated reconstruction methods have made great strides in the past decade, opening windows into the complex intraceullar organization. These advances also mean that computational tools need to be developed to work with these images not just for purposes of visualization but also for biophysical simulations. Here, we describe a pipeline that takes images from electron microscopy as input and produces smooth surface and volume meshes as output. These meshes are suitable for building high-quality finite element simulations of cellular processes modeled by ordinary and partial differential equations, bringing us closer to realizing the goal of generating high-resolution simulations of such phenomena in realistic geometries. We demonstrate the utility of this pipeline by meshing 3D reconstructions of dendritic spines, calculating 1 the curvatures of the different component membranes, and conducting finite-element simulations of reaction-diffusion equations using the generated meshes. The software tools employed in our pipeline are community driven, open source, and free. We believe that technologies such as those presented will enable a new frontier in biophysical simulations in realistic geometries.

Advances in imaging methods such as electron microscopy, tomography and other modalities are enabling high-resolution reconstructions of cellular and organelle geometries. Such advances pave the way for using these geometries for biophysical and mathematical modeling once these data can be represented as a geometric mesh, which, when carefully conditioned, enables the discretization and solution of partial differential equations. In this study, we outline the steps for a naïve user to approach GAMer 2, a mesh generation code written in C++ designed to convert structural datasets to realistic geometric meshes, while preserving the underlying shapes. We present two example cases, 1) mesh generation at the subcellular scale as informed by electron tomography, and 2) meshing a protein with structure from x-ray crystallography. We further demonstrate that the meshes generated by GAMer are suitable for use with numerical methods. Together, this collection of libraries and tools simplifies the process of constructing realistic geometric meshes from structural biology data. SIGNIFICANCE As biophysical structure determination methods improve, the rate of new structural data is increasing. New methods that allow the interpretation, analysis, and reuse of such structural information will thus take on commensurate importance. In particular, geometric meshes, such as those commonly used in graphics and mathematics, can enable a myriad of mathematical analysis. In this work, we describe GAMer 2, a mesh generation library designed for biological datasets. Using GAMer 2 and associated tools PyGAMer and BlendGAMer, biologists can robustly generate computer and algorithm friendly geometric mesh representations informed by structural biology data. We expect that GAMer 2 will be a valuable tool to bring realistic geometries to biophysical models.