Generating Neuron Geometries for Detailed Three-Dimensional Simulations Using AnaMorph

Synthesis of geometrically realistic and watertight neuronal ultrastructure manifolds for in silico modeling

Understanding the intracellular dynamics of brain cells entails performing three-dimensional molecular simulations incorporating ultrastructural models that can capture cellular membrane geometries at nanometer scales. While there is an abundance of neuronal morphologies available online, e.g. from NeuroMorpho.Org, converting those fairly abstract point-and-diameter representations into geometrically realistic and simulation-ready, i.e. watertight, manifolds is challenging. Many neuronal mesh reconstruction methods have been proposed; however, their resulting meshes are either biologically unplausible or non-watertight. We present an effective and unconditionally robust method capable of generating geometrically realistic and watertight surface manifolds of spiny cortical neurons from their morphological descriptions. The robustness of our method is assessed based on a mixed dataset of cortical neurons with a wide variety of morphological classes. The implementation is seamlessly extended and applied to synthetic astrocytic morphologies that are also plausibly biological in detail. Resulting meshes are ultimately used to create volumetric meshes with tetrahedral domains to perform scalable in silico reaction-diffusion simulations for revealing cellular structure-function relationships. Availability and implementation: Our method is implemented in NeuroMorphoVis, a neuroscience-specific open source Blender add-on, making it freely accessible for neuroscience researchers.

Ultraliser: a framework for creating multiscale, high-fidelity and geometrically realistic 3D models for in silico neuroscience

Ultraliser is a neuroscience-specific software framework capable of creating accurate and biologically realistic 3D models of complex neuroscientific structures at intracellular (e.g. mitochondria and endoplasmic reticula), cellular (e.g. neurons and glia) and even multicellular scales of resolution (e.g. cerebral vasculature and minicolumns). Resulting models are exported as triangulated surface meshes and annotated volumes for multiple applications in in silico neuroscience, allowing scalable supercomputer simulations that can unravel intricate cellular structure-function relationships. Ultraliser implements a high-performance and unconditionally robust voxelization engine adapted to create optimized watertight surface meshes and annotated voxel grids from arbitrary non-watertight triangular soups, digitized morphological skeletons or binary volumetric masks. The framework represents a major leap forward in simulation-based neuroscience, making it possible to employ high-resolution 3D structural models for quantification of surface areas and volumes, which are of the utmost importance for cellular and system simulations. The power of Ultraliser is demonstrated with several use cases in which hundreds of models are created for potential application in diverse types of simulations. Ultraliser is publicly released under the GNU GPL3 license on GitHub (BlueBrain/Ultraliser).

SIGNIFICANCE

There is crystal clear evidence on the impact of cell shape on its signaling mechanisms. Structural models can therefore be insightful to realize the function; the more realistic the structure can be, the further we get insights into the function. Creating realistic structural models from existing ones is challenging, particularly when needed for detailed subcellular simulations. We present Ultraliser, a neuroscience-dedicated framework capable of building these structural models with realistic and detailed cellular geometries that can be used for simulations.

Robust quasi-uniform surface meshing of neuronal morphology using line skeleton-based progressive convolution approximation

Creating high-quality polygonal meshes which represent the membrane surface of neurons for both visualization and numerical simulation purposes is an important yet nontrivial task, due to their irregular and complicated structures. In this paper, we develop a novel approach of constructing a watertight 3D mesh from the abstract point-and-diameter representation of the given neuronal morphology. The membrane shape of the neuron is reconstructed by progressively deforming an initial sphere with the guidance of the neuronal skeleton, which can be regarded as a digital sculpting process. To efficiently deform the surface, a local mapping is adopted to simulate the animation skinning. As a result, only the vertices within the region of influence (ROI) of the current skeletal position need to be updated. The ROI is determined based on the finite-support convolution kernel, which is convolved along the line skeleton of the neuron to generate a potential field that further smooths the overall surface at both unidirectional and bifurcating regions. Meanwhile, the mesh quality during the entire evolution is always guaranteed by a set of quasi-uniform rules, which split excessively long edges, collapse undersized ones, and adjust vertices within the tangent plane to produce regular triangles. Additionally, the local vertices density on the result mesh is decided by the radius and curvature of neurites to achieve adaptiveness.

Efficient Simulation of 3D Reaction-Diffusion in Models of Neurons and Networks

Neuronal activity is the result of both the electrophysiology and chemophysiology. A neuron can be well-represented for the purposes of electrophysiological simulation as a tree composed of connected cylinders. This representation is also apt for 1D simulations of their chemophysiology, provided the spatial scale is larger than the diameter of the cylinders and there is radial symmetry. Higher dimensional simulation is necessary to accurately capture the dynamics when these criteria are not met, such as with wave curvature, spines, or diffusion near the soma. We have developed a solution to enable efficient finite volume method simulation of reaction-diffusion kinetics in intracellular 3D regions in neuron and network models and provide an implementation within the NEURON simulator. An accelerated version of the CTNG 3D reconstruction algorithm transforms morphologies suitable for ion-channel based simulations into consistent 3D voxelized regions. Kinetics are then solved using a parallel algorithm based on Douglas-Gunn that handles the irregular 3D geometry of a neuron; these kinetics are coupled to NEURON's 1D mechanisms for ion channels, synapses, pumps, and so forth. The 3D domain may cover the entire cell or selected regions of interest. Simulations with dendritic spines and of the soma reveal details of dynamics that would be missed in a pure 1D simulation. We describe and validate the methods and discuss their performance.

Multi-scale modeling toolbox for single neuron and subcellular activity under Transcranial Magnetic Stimulation

BACKGROUND

Transcranial Magnetic Stimulation (TMS) is a widely used non-invasive brain stimulation method. However, its mechanism of action and the neural response to TMS are still poorly understood. Multi-scale modeling can complement experimental research to study the subcellular neural effects of TMS. At the macroscopic level, sophisticated numerical models exist to estimate the induced electric fields. However, multi-scale computational modeling approaches to predict TMS cellular and subcellular responses, crucial to understanding TMS plasticity inducing protocols, are not available so far.

OBJECTIVE

We develop an open-source multi-scale toolbox Neuron Modeling for TMS (NeMo-TMS) to address this problem.

METHODS

NeMo-TMS generates accurate neuron models from morphological reconstructions, couples them to the external electric fields induced by TMS, and simulates the cellular and subcellular responses of single-pulse and repetitive TMS.

RESULTS

We provide examples showing some of the capabilities of the toolbox.

CONCLUSION

NeMo-TMS toolbox allows researchers a previously not available level of detail and precision in realistically modeling the physical and physiological effects of TMS.

Metaball skinning of synthetic astroglial morphologies into realistic mesh models for visual analytics and in silico simulations

Motivation

Astrocytes, the most abundant glial cells in the mammalian brain, have an instrumental role in developing neuronal circuits. They contribute to the physical structuring of the brain, modulating synaptic activity and maintaining the blood–brain barrier in addition to other significant aspects that impact brain function. Biophysically, detailed astrocytic models are key to unraveling their functional mechanisms via molecular simulations at microscopic scales. Detailed, and complete, biological reconstructions of astrocytic cells are sparse. Nonetheless, data-driven digital reconstruction of astroglial morphologies that are statistically identical to biological counterparts are becoming available. We use those synthetic morphologies to generate astrocytic meshes with realistic geometries, making it possible to perform these simulations.

Results

We present an unconditionally robust method capable of reconstructing high fidelity polygonal meshes of astroglial cells from algorithmically-synthesized morphologies. Our method uses implicit surfaces, or metaballs, to skin the different structural components of astrocytes and then blend them in a seamless fashion. We also provide an end-to-end pipeline to produce optimized two- and three-dimensional meshes for visual analytics and simulations, respectively. The performance of our pipeline has been assessed with a group of 5000 astroglial morphologies and the geometric metrics of the resulting meshes are evaluated. The usability of the meshes is then demonstrated with different use cases.

Availability and implementation

Our metaball skinning algorithm is implemented in Blender 2.82 relying on its Python API (Application Programming Interface). To make it accessible to computational biologists and neuroscientists, the implementation has been integrated into NeuroMorphoVis, an open source and domain specific package that is primarily designed for neuronal morphology visualization and meshing.

Supplementary information

Supplementary data are available at Bioinformatics online.

NeuroBox: Computational Mathematics in Multiscale Neuroscience

The brain is a complex organ operating on multiple scales. From molecular events that inform electrical and biochemical cellular responses, the brain interconnects processes all the way up to the massive network size of billions of brain cells. This strongly coupled, nonlinear, system has been subject to research that has turned increasingly multidisciplinary. The seminal work of Hodgkin and Huxley in the 1950s made use of experimental data to derive a coherent physical model of electrical signaling in neurons, which can be solved using mathematical and computational methods, thus bringing together neuroscience, physics, mathematics, and computer science. Over the last decades numerous projects have been dedicated to modeling and simulation of specific parts of molecular dynamics, neuronal signaling, and neural network behavior. Simulators have been developed around a specific objective and scale, in order to cope with the underlying computational complexity. Often times a dimension reduction approach allows larger scale simulations, this however has the inherent drawback of losing insight into structure-function interplay at the cellular level. This paper gives an overview of the project NeuroBox that has the objective of integrating multiple brain scales and associated physical models into one unified framework. NeuroBox hosts geometry and anatomical reconstruction methods, such that detailed three-dimensional domains can be integrated into numerical simulations of models based on partial differential equations. The project further focusses on deriving numerical methods for handling complex computational domains, and to couple multiple spatial dimensions. The latter allows the user to specify in which parts of the biological problem high-dimensional representations are necessary and where low-dimensional approximations are acceptable. NeuroBox offers workflow user interfaces that are automatically generated with VRL-Studio and can be controlled by non-experts. The project further uses uG4 as the numerical backend, and therefore accesses highly advanced discretization methods as well as hierarchical and scalable numerical solvers for very large neurobiological problems.

How does the presence of neural probes affect extracellular potentials?

OBJECTIVE

Mechanistic modeling of neurons is an essential component of computational neuroscience that enables scientists to simulate, explain, and explore neural activity. The conventional approach to simulation of extracellular neural recordings first computes transmembrane currents using the cable equation and then sums their contribution to model the extracellular potential. This two-step approach relies on the assumption that the extracellular space is an infinite and homogeneous conductive medium, while measurements are performed using neural probes. The main purpose of this paper is to assess to what extent the presence of the neural probes of varying shape and size impacts the extracellular field and how to correct for them.

APPROACH

We apply a detailed modeling framework allowing explicit representation of the neuron and the probe to study the effect of the probes and thereby estimate the effect of ignoring it. We use meshes with simplified neurons and different types of probe and compare the extracellular action potentials with and without the probe in the extracellular space. We then compare various solutions to account for the probes' presence and introduce an efficient probe correction method to include the probe effect in modeling of extracellular potentials.

MAIN RESULTS

Our computations show that microwires hardly influence the extracellular electric field and their effect can therefore be ignored. In contrast, multi-electrode arrays (MEAs) significantly affect the extracellular field by magnifying the recorded potential. While MEAs behave similarly to infinite insulated planes, we find that their effect strongly depends on the neuron-probe alignment and probe orientation.

SIGNIFICANCE

Ignoring the probe effect might be deleterious in some applications, such as neural localization and parameterization of neural models from extracellular recordings. Moreover, the presence of the probe can improve the interpretation of extracellular recordings, by providing a more accurate estimation of the extracellular potential generated by neuronal models.