VolRoverN: enhancing surface and volumetric reconstruction for realistic dynamical simulation of cellular and subcellular function

Synaptic Information Storage Capacity Measured With Information Theory

Variation in the strength of synapses can be quantified by measuring the anatomical properties of synapses. Quantifying precision of synaptic plasticity is fundamental to understanding information storage and retrieval in neural circuits. Synapses from the same axon onto the same dendrite have a common history of coactivation, making them ideal candidates for determining the precision of synaptic plasticity based on the similarity of their physical dimensions. Here, the precision and amount of information stored in synapse dimensions were quantified with Shannon information theory, expanding prior analysis that used signal detection theory (Bartol et al., 2015). The two methods were compared using dendritic spine head volumes in the middle of the stratum radiatum of hippocampal area CA1 as well-defined measures of synaptic strength. Information theory delineated the number of distinguishable synaptic strengths based on nonoverlapping bins of dendritic spine head volumes. Shannon entropy was applied to measure synaptic information storage capacity (SISC) and resulted in a lower bound of 4.1 bits and upper bound of 4.59 bits of information based on 24 distinguishable sizes. We further compared the distribution of distinguishable sizes and a uniform distribution using Kullback-Leibler divergence and discovered that there was a nearly uniform distribution of spine head volumes across the sizes, suggesting optimal use of the distinguishable values. Thus, SISC provides a new analytical measure that can be generalized to probe synaptic strengths and capacity for plasticity in different brain regions of different species and among animals raised in different conditions or during learning. How brain diseases and disorders affect the precision of synaptic plasticity can also be probed.

Long-Term Potentiation Produces a Sustained Expansion of Synaptic Information Storage Capacity in Adult Rat Hippocampus

Long-term potentiation (LTP) has become a standard model for investigating synaptic mechanisms of learning and memory. Increasingly, it is of interest to understand how LTP affects the synaptic information storage capacity of the targeted population of synapses. Here, structural synaptic plasticity during LTP was explored using three-dimensional reconstruction from serial section electron microscopy. Storage capacity was assessed by applying a new analytical approach, Shannon information theory, to delineate the number of functionally distinguishable synaptic strengths. LTP was induced by delta-burst stimulation of perforant pathway inputs to the middle molecular layer of hippocampal dentate granule cells in adult rats. Spine head volumes were measured as predictors of synaptic strength and compared between LTP and control hemispheres at 30 min and 2 hr after the induction of LTP. Synapses from the same axon onto the same dendrite were used to determine the precision of synaptic plasticity based on the similarity of their physical dimensions. Shannon entropy was measured by exploiting the frequency of spine heads in functionally distinguishable sizes to assess the degree to which LTP altered the number of bits of information storage. Outcomes from these analyses reveal that LTP expanded storage capacity; the distribution of spine head volumes was increased from 2 bits in controls to 3 bits at 30 min and 2.7 bits at 2 hr after the induction of LTP. Furthermore, the distribution of spine head volumes was more uniform across the increased number of functionally distinguishable sizes following LTP, thus achieving more efficient use of coding space across the population of synapses.

Cristae formation is a mechanical buckling event controlled by the inner mitochondrial membrane lipidome

Cristae are high-curvature structures in the inner mitochondrial membrane (IMM) that are crucial for ATP production. While cristae-shaping proteins have been defined, analogous lipid-based mechanisms have yet to be elucidated. Here, we combine experimental lipidome dissection with multi-scale modeling to investigate how lipid interactions dictate IMM morphology and ATP generation. When modulating phospholipid (PL) saturation in engineered yeast strains, we observed a surprisingly abrupt breakpoint in IMM topology driven by a continuous loss of ATP synthase organization at cristae ridges. We found that cardiolipin (CL) specifically buffers the inner mitochondrial membrane against curvature loss, an effect that is independent of ATP synthase dimerization. To explain this interaction, we developed a continuum model for cristae tubule formation that integrates both lipid and protein-mediated curvatures. This model highlighted a snapthrough instability, which drives IMM collapse upon small changes in membrane properties. We also showed that cardiolipin is essential in low-oxygen conditions that promote PL saturation. These results demonstrate that the mechanical function of cardiolipin is dependent on the surrounding lipid and protein components of the IMM.

Cristae formation is a mechanical buckling event controlled by the inner membrane lipidome

Cristae are high curvature structures in the inner mitochondrial membrane (IMM) that are crucial for ATP production. While cristae-shaping proteins have been defined, analogous mechanisms for lipids have yet to be elucidated. Here we combine experimental lipidome dissection with multi-scale modeling to investigate how lipid interactions dictate IMM morphology and ATP generation. When modulating phospholipid (PL) saturation in engineered yeast strains, we observed a surprisingly abrupt breakpoint in IMM topology driven by a continuous loss of ATP synthase organization at cristae ridges. We found that cardiolipin (CL) specifically buffers the IMM against curvature loss, an effect that is independent of ATP synthase dimerization. To explain this interaction, we developed a continuum model for cristae tubule formation that integrates both lipid and protein-mediated curvatures. The model highlighted a snapthrough instability, which drives IMM collapse upon small changes in membrane properties. We also showed that CL is essential in low oxygen conditions that promote PL saturation. These results demonstrate that the mechanical function of CL is dependent on the surrounding lipid and protein components of the IMM.

CellRemorph: A Toolkit for Transforming, Selecting, and Slicing 3D Cell Structures on the Road to Morphologically Detailed Astrocyte Simulations

Understanding functions of astrocytes can be greatly enhanced by building and simulating computational models that capture their morphological details. Novel computational tools enable utilization of existing morphological data of astrocytes and building models that have appropriate level of details for specific simulation purposes. In addition to analyzing existing computational tools for constructing, transforming, and assessing astrocyte morphologies, we present here the CellRemorph toolkit implemented as an add-on for Blender, a 3D modeling platform increasingly recognized for its utility for manipulating 3D biological data. To our knowledge, CellRemorph is the first toolkit for transforming astrocyte morphologies from polygonal surface meshes into adjustable surface point clouds and vice versa, precisely selecting nanoprocesses, and slicing morphologies into segments with equal surface areas or volumes. CellRemorph is an open-source toolkit under the GNU General Public License and easily accessible via an intuitive graphical user interface. CellRemorph will be a valuable addition to other Blender add-ons, providing novel functionality that facilitates the creation of realistic astrocyte morphologies for different types of morphologically detailed simulations elucidating the role of astrocytes both in health and disease.

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.

Efficient metadata mining of web-accessible neural morphologies

Advancements in neuroscience research have led to steadily accelerating data production and sharing. The online community repository of neural reconstructions NeuroMorpho.Org grew from fewer than 1000 digitally traced neurons in 2006 to more than 140,000 cells today, including glia that now constitute 10.1% of the content. Every reconstruction consists of a detailed 3D representation of branch geometry and connectivity in a standardized format, from which a collection of morphometric features is extracted and stored. Moreover, each entry in the database is accompanied by rich metadata annotation describing the animal subject, anatomy, and experimental details. The rapid expansion of this resource in the past decade was accompanied by a parallel rise in the complexity of the available information, creating both opportunities and challenges for knowledge mining. Here, we introduce a new summary reporting functionality, allowing NeuroMorpho.Org users to efficiently download digests of metadata and morphometrics from multiple groups of similar cells for further analysis. We demonstrate the capabilities of the tool for both glia and neurons and present an illustrative statistical analysis of the resulting data.

Morphological principles of neuronal mitochondria

In the highly dynamic metabolic landscape of a neuron, mitochondrial membrane architectures can provide critical insight into the unique energy balance of the cell. Current theoretical calculations of functional outputs like adenosine triphosphate and heat often represent mitochondria as idealized geometries, and therefore, can miscalculate the metabolic fluxes. To analyze mitochondrial morphology in neurons of mouse cerebellum neuropil, 3D tracings of complete synaptic and axonal mitochondria were constructed using a database of serial transmission electron microscopy (TEM) tomography images and converted to watertight meshes with minimal distortion of the original microscopy volumes with a granularity of 1.64 nanometer isotropic voxels. The resulting in-silico representations were subsequently quantified by differential geometry methods in terms of the mean and Gaussian curvatures, surface areas, volumes, and membrane motifs, all of which can alter the metabolic output of the organelle. Finally, we identify structural motifs present across this population of mitochondria, which may contribute to future modeling studies of mitochondrial physiology and metabolism in neurons.

Predominantly linear summation of metabotropic postsynaptic potentials follows coactivation of neurogliaform interneurons

Summation of ionotropic receptor-mediated responses is critical in neuronal computation by shaping input-output characteristics of neurons. However, arithmetics of summation for metabotropic signals are not known. We characterized the combined ionotropic and metabotropic output of neocortical neurogliaform cells (NGFCs) using electrophysiological and anatomical methods in the rat cerebral cortex. These experiments revealed that GABA receptors are activated outside release sites and confirmed coactivation of putative NGFCs in superficial cortical layers in vivo. Triple recordings from presynaptic NGFCs converging to a postsynaptic neuron revealed sublinear summation of ionotropic GABA_A responses and linear summation of metabotropic GABA_B responses. Based on a model combining properties of volume transmission and distributions of all NGFC axon terminals, we predict that in 83% of cases one or two NGFCs can provide input to a point in the neuropil. We suggest that interactions of metabotropic GABAergic responses remain linear even if most superficial layer interneurons specialized to recruit GABA_B receptors are simultaneously active.

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

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.

Long-term potentiation expands information content of hippocampal dentate gyrus synapses

An approach combining signal detection theory and precise 3D reconstructions from serial section electron microscopy (3DEM) was used to investigate synaptic plasticity and information storage capacity at medial perforant path synapses in adult hippocampal dentate gyrus in vivo. Induction of long-term potentiation (LTP) markedly increased the frequencies of both small and large spines measured 30 minutes later. This bidirectional expansion resulted in heterosynaptic counterbalancing of total synaptic area per unit length of granule cell dendrite. Control hemispheres exhibited 6.5 distinct spine sizes for 2.7 bits of storage capacity while LTP resulted in 12.9 distinct spine sizes (3.7 bits). In contrast, control hippocampal CA1 synapses exhibited 4.7 bits with much greater synaptic precision than either control or potentiated dentate gyrus synapses. Thus, synaptic plasticity altered total capacity, yet hippocampal subregions differed dramatically in their synaptic information storage capacity, reflecting their diverse functions and activation histories.

Heterogeneous anisotropic magnetic susceptibility of the myelin-water layers causes local magnetic field perturbations in axons

One goal of MRI is to determine the myelin water fraction in neural tissue. One approach is to measure the reduction in T₂ * arising from microscopic perturbations in the magnetic field caused by heterogeneities in the magnetic susceptibility of myelin. In this paper, analytic expressions for the induced magnetic field distribution are derived within and around an axon, assuming that the myelin susceptibility is anisotropic. Previous models considered the susceptibility to be piecewise continuous, whereas this model considers a sinusoidally varying susceptibility. Many conclusions are common in both models. When the magnetic field is applied perpendicular to the axon, the magnetic field in the intraaxonal space is uniformly perturbed, the magnetic field in the myelin sheath oscillates between the lipid and water layers, and the magnetic field in the extracellular space just outside the myelin sheath is heterogeneous. These field heterogeneities cause the spins to dephase, shortening T₂ *. When the magnetic field is applied along the axon, the field is homogeneous within water-filled regions, including between lipid layers. Therefore the spins do not dephase and the magnetic susceptibility has no effect on T₂ *. Generally, the response of an axon is given as the superposition of these two contributions. The sinusoidal model uses a different set of approximations compared with the piecewise model, so their common predictions indicate that the models are not too sensitive to the details of the myelin-water distribution. Other predictions, such as the sensitivity to water diffusion between myelin and water layers, may highlight differences between the two approaches. Copyright © 2016 John Wiley & Sons, Ltd.

Nanoconnectomic upper bound on the variability of synaptic plasticity

Information in a computer is quantified by the number of bits that can be stored and recovered. An important question about the brain is how much information can be stored at a synapse through synaptic plasticity, which depends on the history of probabilistic synaptic activity. The strong correlation between size and efficacy of a synapse allowed us to estimate the variability of synaptic plasticity. In an EM reconstruction of hippocampal neuropil we found single axons making two or more synaptic contacts onto the same dendrites, having shared histories of presynaptic and postsynaptic activity. The spine heads and neck diameters, but not neck lengths, of these pairs were nearly identical in size. We found that there is a minimum of 26 distinguishable synaptic strengths, corresponding to storing 4.7 bits of information at each synapse. Because of stochastic variability of synaptic activation the observed precision requires averaging activity over several minutes.

DOI:http://dx.doi.org/10.7554/eLife.10778.001

What is the memory capacity of a human brain? The storage capacity in a computer memory is measured in bits, each of which can have a value of 0 or 1. In the brain, information is stored in the form of synaptic strength, a measure of how strongly activity in one neuron influences another neuron to which it is connected. The number of different strengths can be measured in bits. The total storage capacity of the brain therefore depends on both the number of synapses and the number of distinguishable synaptic strengths.

Structurally, neurons consist of a cell body that influences other neurons through a cable-like axon. The cell body bears numerous short branches called dendrites, which are covered in tiny protrusions, or “spines”. Most excitatory synapses are formed between the axon of one neuron and a dendritic spine on another. When two neurons on either side of a synapse are active simultaneously, that synapse becomes stronger, a form of memory. The dendritic spine also becomes larger to accommodate the extra molecular machinery needed to support a stronger synapse.

Some axons form two or more synapses with the same dendrite, but on different dendritic spines. These synapses should be the same strength because they will have experienced the same history of neural activity. Bartol et al. used a technique called serial section electron microscopy to create a 3D reconstruction of part of the brain that allowed the sizes of the dendritic spines these synapses form on to be compared. This revealed that the synaptic areas and volumes of the spine heads were nearly identical. This remarkable similarity can be used to estimate the number of bits of information that a single synapse can store, since the size of dendritic spines and their synapses can be used as proxies for synaptic strength.

Measurements in a small cube of brain tissue revealed 26 different dendritic spine sizes, each associated with a distinct synaptic strength. This number translates into a storage capacity of roughly 4.7 bits of information per synapse. This estimate is markedly higher than previous suggestions. It implies that the total memory capacity of the brain – with its many trillions of synapses – may have been underestimated by an order of magnitude. Additional measurements in the same and other brain regions are needed to confirm this possibility.

DOI:http://dx.doi.org/10.7554/eLife.10778.002

Computational reconstitution of spine calcium transients from individual proteins

We have built a stochastic model in the program MCell that simulates Ca(2+) transients in spines from the principal molecular components believed to control Ca(2+) entry and exit. Proteins, with their kinetic models, are located within two segments of dendrites containing 88 intact spines, centered in a fully reconstructed 6 × 6 × 5 μm(3) cube of hippocampal neuropil. Protein components include AMPA- and NMDA-type glutamate receptors, L- and R-type voltage-dependent Ca(2+) channels, Na(+)/Ca(2+) exchangers, plasma membrane Ca(2+) ATPases, smooth endoplasmic reticulum Ca(2+) ATPases, immobile Ca(2+) buffers, and calbindin. Kinetic models for each protein were taken from published studies of the isolated proteins in vitro. For simulation of electrical stimuli, the time course of voltage changes in the dendritic spine was generated with the desired stimulus in the program NEURON. Voltage-dependent parameters were then continuously re-adjusted during simulations in MCell to reproduce the effects of the stimulus. Nine parameters of the model were optimized within realistic experimental limits by a process that compared results of simulations to published data. We find that simulations in the optimized model reproduce the timing and amplitude of Ca(2+) transients measured experimentally in intact neurons. Thus, we demonstrate that the characteristics of individual isolated proteins determined in vitro can accurately reproduce the dynamics of experimentally measured Ca(2+) transients in spines. The model will provide a test bed for exploring the roles of additional proteins that regulate Ca(2+) influx into spines and for studying the behavior of protein targets in the spine that are regulated by Ca(2+) influx.

A resource from 3D electron microscopy of hippocampal neuropil for user training and tool development

Resurgent interest in synaptic circuitry and plasticity has emphasized the importance of 3D reconstruction from serial section electron microscopy (3DEM). Three volumes of hippocampal CA1 neuropil from adult rat were imaged at X-Y resolution of ~2 nm on serial sections of ~50-60 nm thickness. These are the first densely reconstructed hippocampal volumes. All axons, dendrites, glia, and synapses were reconstructed in a cube (~10 μm(3)) surrounding a large dendritic spine, a cylinder (~43 μm(3)) surrounding an oblique dendritic segment (3.4 μm long), and a parallelepiped (~178 μm(3)) surrounding an apical dendritic segment (4.9 μm long). The data provide standards for identifying ultrastructural objects in 3DEM, realistic reconstructions for modeling biophysical properties of synaptic transmission, and a test bed for enhancing reconstruction tools. Representative synapses are quantified from varying section planes, and microtubules, polyribosomes, smooth endoplasmic reticulum, and endosomes are identified and reconstructed in a subset of dendrites. The original images, traces, and Reconstruct software and files are freely available and visualized at the Open Connectome Project (Data Citation 1).