Dataset Information

Reconstruction and Simulation of a Scaffold Model of the Cerebellar Network.

ABSTRACT: Reconstructing neuronal microcircuits through computational models is fundamental to simulate local neuronal dynamics. Here a scaffold model of the cerebellum has been developed in order to flexibly place neurons in space, connect them synaptically, and endow neurons and synapses with biologically-grounded mechanisms. The scaffold model can keep neuronal morphology separated from network connectivity, which can in turn be obtained from convergence/divergence ratios and axonal/dendritic field 3D geometries. We first tested the scaffold on the cerebellar microcircuit, which presents a challenging 3D organization, at the same time providing appropriate datasets to validate emerging network behaviors. The scaffold was designed to integrate the cerebellar cortex with deep cerebellar nuclei (DCN), including different neuronal types: Golgi cells, granule cells, Purkinje cells, stellate cells, basket cells, and DCN principal cells. Mossy fiber inputs were conveyed through the glomeruli. An anisotropic volume (0.077 mm³) of mouse cerebellum was reconstructed, in which point-neuron models were tuned toward the specific discharge properties of neurons and were connected by exponentially decaying excitatory and inhibitory synapses. Simulations using both pyNEST and pyNEURON showed the emergence of organized spatio-temporal patterns of neuronal activity similar to those revealed experimentally in response to background noise and burst stimulation of mossy fiber bundles. Different configurations of granular and molecular layer connectivity consistently modified neuronal activation patterns, revealing the importance of structural constraints for cerebellar network functioning. The scaffold provided thus an effective workflow accounting for the complex architecture of the cerebellar network. In principle, the scaffold can incorporate cellular mechanisms at multiple levels of detail and be tuned to test different structural and functional hypotheses. A future implementation using detailed 3D multi-compartment neuron models and dynamic synapses will be needed to investigate the impact of single neuron properties on network computation.

SUBMITTER: Casali S

PROVIDER: S-EPMC6530631 | biostudies-literature | 2019

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Reconstruction and Simulation of a Scaffold Model of the Cerebellar Network.

Casali Stefano S Marenzi Elisa E Medini Chaitanya C Casellato Claudia C D'Angelo Egidio E

Frontiers in neuroinformatics 20190515

Reconstructing neuronal microcircuits through computational models is fundamental to simulate local neuronal dynamics. Here a <i>scaffold model</i> of the cerebellum has been developed in order to flexibly place neurons in space, connect them synaptically, and endow neurons and synapses with biologically-grounded mechanisms. The scaffold model can keep neuronal morphology separated from network connectivity, which can in turn be obtained from convergence/divergence ratios and axonal/dendritic fi ...[more]

PMID: 31156416

Similar Datasets

Project description:Large-scale simulation of detailed computational models of neuronal microcircuits plays a prominent role in reproducing and predicting the dynamics of the microcircuits. To reconstruct a microcircuit, one must choose neuron and synapse models, placements, connectivity, and numerical simulation methods according to anatomical and physiological constraints. For reconstruction and refinement, it is useful to be able to replace one module easily while leaving the others as they are. One way to achieve this is via a scaffolding approach, in which a simulation code is built on independent modules for placements, connections, and network simulations. Owing to the modularity of functions, this approach enables researchers to improve the performance of the entire simulation by simply replacing a problematic module with an improved one. Casali et al. (2019) developed a spiking network model of the cerebellar microcircuit using this approach, and while it reproduces electrophysiological properties of cerebellar neurons, it takes too much computational time. Here, we followed this scaffolding approach and replaced the simulation module with an accelerated version on graphics processing units (GPUs). Our cerebellar scaffold model ran roughly 100 times faster than the original version. In fact, our model is able to run faster than real time, with good weak and strong scaling properties. To demonstrate an application of real-time simulation, we implemented synaptic plasticity mechanisms at parallel fiber-Purkinje cell synapses, and carried out simulation of behavioral experiments known as gain adaptation of optokinetic response. We confirmed that the computer simulation reproduced experimental findings while being completed in real time. Actually, a computer simulation for 2 s of the biological time completed within 750 ms. These results suggest that the scaffolding approach is a promising concept for gradual development and refactoring of simulation codes for large-scale elaborate microcircuits. Moreover, a real-time version of the cerebellar scaffold model, which is enabled by parallel computing technology owing to GPUs, may be useful for large-scale simulations and engineering applications that require real-time signal processing and motor control.

Project description:BackgroundPositron-emission tomography (PET) simulators are frequently used for development and performance evaluation of segmentation methods or quantitative uptake metrics. To date, most PET simulation tools are based on Monte Carlo simulations, which are computationally demanding. Other analytical simulation tools lack the implementation of time of flight (TOF) or resolution modelling (RM). In this study, a fast and easy-to-use PET simulation-reconstruction package, SiMulAtion and ReconsTruction (SMART)-PET, is developed and validated, which includes both TOF and RM. SMART-PET, its documentation and instructions to calibrate the tool to a specific PET/CT system are available on Zenodo. SMART-PET allows the fast generation of 3D PET images. As input, it requires one image representing the activity distribution and one representing the corresponding CT image/attenuation map. It allows the user to adjust different parameters, such as reconstruction settings (TOF/RM), noise level or scan duration. Furthermore, a random spatial shift can be included, representing patient repositioning. To evaluate the tool, simulated images were compared with real scan data of the NEMA NU 2 image quality phantom. The scan was acquired as a 60-min list-mode scan and reconstructed with and without TOF and/or RM. For every reconstruction setting, ten statistically equivalent images, representing 30, 60, 120 and 300 s scan duration, were generated. Simulated and real-scan data were compared regarding coefficient of variation in the phantom background and activity recovery coefficients (RCs) of the spheres. Furthermore, standard deviation images of each of the ten statistically equivalent images were compared.ResultsSMART-PET produces images comparable to actual phantom data. The image characteristics of simulated and real PET images varied in similar ways as function of reconstruction protocols and noise levels. The change in image noise with variation of simulated TOF settings followed the theoretically expected behaviour. RC as function of sphere size agreed within 0.3-11% between simulated and actual phantom data.ConclusionsSMART-PET allows for rapid and easy simulation of PET data. The user can change various acquisition and reconstruction settings (including RM and TOF) and noise levels. The images obtained show similar image characteristics as those seen in actual phantom data.

Project description:The distribution of multiscale pores and fractures in coal and rock is an important basis for reflecting the capacity of fluid flow in coal seam seepage passages. Accurate extraction and qualitative and quantitative analysis of pore-fracture structures are helpful in revealing the flow characteristics of fluid in seepage channels. The relationship between pore and fracture connectivity can provide a scientific reference for optimizing coal seam water injection parameters. Therefore, to analyse the change in permeability caused by the variability in the coal pore-fracture network structure, a CT scanning technique was used to scan coal samples from the Leijia District, Fuxin. A total of 720 sets of original images were collected, a median filter was used to filter out the noise in the obtained images, and to form the basis of a model, the reconstruction and analysis of the three-dimensional pore-fracture morphology of coal samples were carried out. A pore-fracture network model of the coal body was extracted at different scales. Using the maximum sphere algorithm combined with the coordination number, the effect of different quantitative relationships between pore size and pore throat channel permeability was studied. Avizo software was used to simulate the flow path of fluid in the seepage channels. The change trend of the fluid velocity between different seepage channels was discussed. The results of the pore-fracture network models at different scales show that the pore-fracture structure is nonuniform and vertically connected, and the pores are connected at connecting points. The pore size distribution ranges from 104 μm to 9425 μm. The pore throat channel length distribution ranges from 4206 μm to 48073 μm. The size of the coordination number determines the connectivity and thus the porosity of the coal seam. The more connected pore channels there are, the larger the pore diameters and the stronger the percolation ability. During flow in the seepage channels of the coal, the velocity range is divided into a low-speed region, medium-speed region and high-speed region. The fluid seepage in the coal seam is driven by the following factors: pore connectivity > pore and pore throat dimensions > pore and pore throat structure distribution. Ultimately, the pore radius and pore connectivity directly affect the permeability of the coal seam.

Project description:The empirical harmonic potential function of elastic network models (ENMs) is augmented by three- and four-body interactions as well as by a parameter-free connection rule. In the new bend-twist-stretch (BTS) model the complexity of the parametrization is shifted from the spatial level of detail to the potential function, enabling an arbitrary coarse graining of the network. Compared to distance cutoff-based Hookean springs, the approach yields a more stable parametrization of coarse-grained ENMs for biomolecular dynamics. Traditional ENMs give rise to unbounded zero-frequency vibrations when (pseudo)atoms are connected to fewer than three neighbors. A large cutoff is therefore chosen in an ENM (about twice the average nearest-neighbor distance), resulting in many false-positive connections that reduce the spatial detail that can be resolved. More importantly, the required three-neighbor connectedness also limits the coarse graining, i.e., the network must be dense, even in the case of low-resolution structures that exhibit few spatial features. The new BTS model achieves such coarse graining by extending the ENM potential to include three-and four-atom interactions (bending and twisting, respectively) in addition to the traditional two-atom stretching. Thus, the BTS model enables reliable modeling of any three-dimensional graph irrespective of the atom connectedness. The additional potential terms were parametrized using continuum elastic theory of elastic rods, and the distance cutoff was replaced by a competitive Hebb connection rule, setting all free parameters in the model. We validate the approach on a carbon-alpha representation of adenylate kinase and illustrate its use with electron microscopy maps of E. coli RNA polymerase, E. coli ribosome, and eukaryotic chaperonin containing T-complex polypeptide 1, which were difficult to model with traditional ENMs. For adenylate kinase, we find excellent reproduction (>90% overlap) of the ENM modes and B factors when BTS is applied to the carbon-alpha representation as well as to coarser descriptions. For the volumetric maps, coarse BTS yields similar motions (70%-90% overlap) to those obtained from significantly denser representations with ENM. Our Python-based algorithms of ENM and BTS implementations are freely available.

Dataset Information

Reconstruction and Simulation of a Scaffold Model of the Cerebellar Network.

Publications

Reconstruction and Simulation of a Scaffold Model of the Cerebellar Network.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets