Disruptions of energy supply to the brain are associated with many neurodegenerative pathologies and are difficult to study due to numerous interlinked metabolic pathways. We explored the effects of diminished energy supply on brain metabolism using a computational model of the neuro-glia-vasculature ensemble, in the form of a neuron, an astrocyte and local blood supply. As a case study, we investigated the glucose transporter type-1 deficiency syndrome (GLUT1-DS), a childhood affliction characterized by impaired glucose utilization and associated with phenotypes including seizures. Compared to neurons, astrocytes exhibited markedly higher metabolite concentration variabilities for all but a few redox species. This effect could signal a role for astrocytes in absorbing the shock of blood nutrient fluctuations. Redox balances were disrupted in GLUT1-DS with lower levels of reducing equivalent carriers NADH and ATP. The best non-glucose nutrient or pharmacotherapies for re-establishing redox normalcy involved lactate, the keto-diet (β-hydroxybutyrate), NAD and Q10 supplementation, suggesting a possible glucose sparing mechanism. GLUT1-DS seizures resulted from after-discharge neuronal firing caused by post-stimulus ATP reductions and impaired Na+/K+-ATPase, which can be rescued by restoring either normal glucose or by relatively small increases in neuronal ATP.

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 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.

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.

MOTIVATION: Accurate morphological models of brain vasculature are key to modeling and simulating cerebral blood flow in realistic vascular networks. This in silico approach is fundamental to revealing the principles of neurovascular coupling. Validating those vascular morphologies entails performing certain visual analysis tasks that cannot be accomplished with generic visualization frameworks. This limitation has a substantial impact on the accuracy of the vascular models employed in the simulation. RESULTS: We present VessMorphoVis, an integrated suite of toolboxes for interactive visualization and analysis of vast brain vascular networks represented by morphological graphs segmented originally from imaging or microscopy stacks. Our workflow leverages the outstanding potentials of Blender, aiming to establish an integrated, extensible and domain-specific framework capable of interactive visualization, analysis, repair, high-fidelity meshing and high-quality rendering of vascular morphologies. Based on the initial feedback of the users, we anticipate that our framework will be an essential component in vascular modeling and simulation in the future, filling a gap that is at present largely unfulfilled. AVAILABILITY AND IMPLEMENTATION: VessMorphoVis is freely available under the GNU public license on Github at https://github.com/BlueBrain/VessMorphoVis. The morphology analysis, visualization, meshing and rendering modules are implemented as an add-on for Blender 2.8 based on its Python API (application programming interface). The add-on functionality is made available to users through an intuitive graphical user interface, as well as through exhaustive configuration files calling the API via a feature-rich command line interface running Blender in background mode. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Glucose is the main energy substrate in the adult brain under normal conditions. Accumulating evidence, however, indicates that lactate produced in astrocytes (a type of glial cell) can also fuel neuronal activity. The quantitative aspects of this so-called astrocyte-neuron lactate shuttle (ANLS) are still debated. To address this question, we developed a detailed biophysical model of the brain's metabolic interactions. Our model integrates three modeling approaches, the Buxton-Wang model of vascular dynamics, the Hodgkin-Huxley formulation of neuronal membrane excitability and a biophysical model of metabolic pathways. This approach provides a template for large-scale simulations of the neuron-glia-vasculature (NGV) ensemble, and for the first time integrates the respective timescales at which energy metabolism and neuronal excitability occur. The model is constrained by relative neuronal and astrocytic oxygen and glucose utilization, by the concentration of metabolites at rest and by the temporal dynamics of NADH upon activation. These constraints produced four observations. First, a transfer of lactate from astrocytes to neurons emerged in response to activity. Second, constrained by activity-dependent NADH transients, neuronal oxidative metabolism increased first upon activation with a subsequent delayed astrocytic glycolysis increase. Third, the model correctly predicted the dynamics of extracellular lactate and oxygen as observed in vivo in rats. Fourth, the model correctly predicted the temporal dynamics of tissue lactate, of tissue glucose and oxygen consumption, and of the BOLD signal as reported in human studies. These findings not only support the ANLS hypothesis but also provide a quantitative mathematical description of the metabolic activation in neurons and glial cells, as well as of the macroscopic measurements obtained during brain imaging.