Dataset Information

1D-3D hybrid modeling-from multi-compartment models to full resolution models in space and time.

ABSTRACT: Investigation of cellular and network dynamics in the brain by means of modeling and simulation has evolved into a highly interdisciplinary field, that uses sophisticated modeling and simulation approaches to understand distinct areas of brain function. Depending on the underlying complexity, these models vary in their level of detail, in order to cope with the attached computational cost. Hence for large network simulations, single neurons are typically reduced to time-dependent signal processors, dismissing the spatial aspect of each cell. For single cell or networks with relatively small numbers of neurons, general purpose simulators allow for space and time-dependent simulations of electrical signal processing, based on the cable equation theory. An emerging field in Computational Neuroscience encompasses a new level of detail by incorporating the full three-dimensional morphology of cells and organelles into three-dimensional, space and time-dependent, simulations. While every approach has its advantages and limitations, such as computational cost, integrated and methods-spanning simulation approaches, depending on the network size could establish new ways to investigate the brain. In this paper we present a hybrid simulation approach, that makes use of reduced 1D-models using e.g., the NEURON simulator-which couples to fully resolved models for simulating cellular and sub-cellular dynamics, including the detailed three-dimensional morphology of neurons and organelles. In order to couple 1D- and 3D-simulations, we present a geometry-, membrane potential- and intracellular concentration mapping framework, with which graph- based morphologies, e.g., in the swc- or hoc-format, are mapped to full surface and volume representations of the neuron and computational data from 1D-simulations can be used as boundary conditions for full 3D simulations and vice versa. Thus, established models and data, based on general purpose 1D-simulators, can be directly coupled to the emerging field of fully resolved, highly detailed 3D-modeling approaches. We present the developed general framework for 1D/3D hybrid modeling and apply it to investigate electrically active neurons and their intracellular spatio-temporal calcium dynamics.

SUBMITTER: Grein S

PROVIDER: S-EPMC4114301 | biostudies-literature | 2014

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Kinetic modeling is widely used to analyze dynamic imaging data, estimating kinetic parameters that quantify functional or physiologic processes in vivo. Typical kinetic models give rise to nonlinear solution equations in multiple dimensions, presenting a complex fitting environment. This work generalizes previously described separable nonlinear least-squares techniques for fitting serial compartment models with up to three tissue compartments and five rate parameters.The approach maximally separates the linear and nonlinear aspects of the modeling equations, using a formulation modified from previous basis function methods to avoid a potential mathematical degeneracy. A fast and robust algorithm for solving the linear subproblem with full user-defined constraints is also presented. The generalized separable parameter space technique effectively reduces the dimensionality of the nonlinear fitting problem to one dimension for 2K-3K compartment models, and to two dimensions for 4K-5K models.Exhaustive search fits, which guarantee identification of the true global minimum fit, required approximately 10 ms for 2K-3K and 1.1 s for 4K-5K models, respectively. The technique is also amenable to fast gradient-descent iterative fitting algorithms, where the reduced dimensionality offers improved convergence properties. The objective function for the separable parameter space nonlinear subproblem was characterized and found to be generally well-behaved with a well-defined global minimum. Separable parameter space fits with the Levenberg-Marquardt algorithm required fewer iterations than comparable fits for conventional model formulations, averaging 1 and 7 ms for 2K-3K and 4K-5K models, respectively. Sensitivity to initial conditions was likewise reduced.The separable parameter space techniques described herein generalize previously described techniques to encompass 1K-5K compartment models, enable robust solution of the linear subproblem with full user-defined constraints, and are amenable to rapid and robust fitting using iterative gradient-descent type algorithms.

Project description:OBJECTIVE:Focal cortical dysplasia (FCD) lesion detection and subtyping remain challenging on conventional MRI. New diffusion models such as the spherical mean technique (SMT) and neurite orientation dispersion and density imaging (NODDI) provide measurements that potentially produce more specific maps of abnormal tissue microstructure. This study aims to assess the SMT and NODDI maps for computational and radiological lesion characterization compared to standard fractional anisotropy (FA) and mean diffusivity (MD). METHODS:SMT, NODDI, FA, and MD maps were calculated for 33 pediatric patients with suspected FCD (18 histologically confirmed). Two neuroradiologists scored lesion visibility on clinical images and diffusion maps. Signal profile changes within lesions and homologous regions were quantified using a surface-based approach. Diffusion parameter changes at multiple cortical depths were statistically compared between FCD type IIa and type IIb. RESULTS:Compared to fluid-attenuated inversion recovery (FLAIR) or T1-weighted imaging, lesions conspicuity on NODDI intracellular volume fraction (ICVF) maps was better/equal/worse in 5/14/14 patients, respectively, while on SMT intra-neurite volume fraction (INVF) in 3/3/27. Compared to FA or MD, lesion conspicuity on the ICVF was better/equal/worse in 27/4/2, while on the INVF in 20/7/6. Quantitative signal profiling demonstrated significant ICVF and INVF reductions in the lesions, whereas SMT microscopic mean, radial, and axial diffusivities were significantly increased. FCD type IIb exhibited greater changes than FCD type IIa. No changes were detected on FA or MD profiles. SIGNIFICANCE:FCD lesion-specific signal changes were found in ICVF and INVF but not in FA and MD maps. ICVF and INVF showed greater contrast than FLAIR in some cases and had consistent signal changes specific to FCD, suggesting that they could improve current presurgical pediatric epilepsy imaging protocols and can provide features useful for automated lesion detection.

Project description:Background: Multi-compartment cable models of auditory nerve fibers have been developed to assist in the improvement of cochlear implants. With the advancement of computational technology and the results obtained from in vivo and in vitro experiments, these models have evolved to incorporate a considerable degree of morphological and physiological details. They have also been combined with three-dimensional volume conduction models of the cochlea to simulate neural responses to electrical stimulation. However, no specific rules have been provided on choosing the appropriate cable model, and most models adopted in recent studies were chosen without a specific reason or by inheritance. Methods: Three of the most cited biophysical multi-compartment cable models of the human auditory nerve, i.e., Rattay et al. (2001b), Briaire and Frijns (2005), and Smit et al. (2010), were implemented in this study. Several properties of single fibers were compared among the three models, including threshold, conduction velocity, action potential shape, latency, refractory properties, as well as stochastic and temporal behaviors. Experimental results regarding these properties were also included as a reference for comparison. Results: For monophasic single-pulse stimulation, the ratio of anodic vs. cathodic thresholds in all models was within the experimental range despite a much larger ratio in the model by Briaire and Frijns. For biphasic pulse-train stimulation, thresholds as a function of both pulse rate and pulse duration differed between the models, but none matched the experimental observations even coarsely. Similarly, for all other properties including the conduction velocity, action potential shape, and latency, the models presented different outcomes and not all of them fell within the range observed in experiments. Conclusions: While all three models presented similar values in certain single fiber properties to those obtained in experiments, none matched all experimental observations satisfactorily. In particular, the adaptation and temporal integration behaviors were completely missing in all models. Further extensions and analyses are required to explain and simulate realistic auditory nerve fiber responses to electrical stimulation.

Project description:Tuberculosis is a worldwide health problem with 2 billion people infected with Mycobacterium tuberculosis (Mtb, the bacteria causing TB). The hallmark of infection is the emergence of organized structures of immune cells forming primarily in the lung in response to infection. Granulomas physically contain and immunologically restrain bacteria that cannot be cleared. We have developed several models that spatially characterize the dynamics of the host-mycobacterial interaction, and identified mechanisms that control granuloma formation and development. In particular, we published several agent-based models (ABMs) of granuloma formation in TB that include many subtypes of T cell populations, macrophages as well as key cytokine and chemokine effector molecules. These ABM studies emphasize the important role of T-cell related mechanisms in infection progression, such as magnitude and timing of T cell recruitment, and macrophage activation. In these models, the priming and recruitment of T cells from the lung draining lymph node (LN) was captured phenomenologically. In addition to these ABM studies, we have also developed several multi-organ models using ODEs to examine trafficking of cells between, for example, the lung and LN. While we can predict temporal dynamic behaviors, those models are not coupled to the spatial aspects of granuloma. To this end, we have developed a multi-organ model that is hybrid: an ABM for the lung compartment and a non-linear system of ODE representing the lymph node compartment. This hybrid multi-organ approach to study TB granuloma formation in the lung and immune priming in the LN allows us to dissect protective mechanisms that cannot be achieved using the single compartment or multi-compartment ODE system. The main finding of this work is that trafficking of important cells known as antigen presenting cells from the lung to the lymph node is a key control mechanism for protective immunity: the entire spectrum of infection outcomes can be regulated by key immune cell migration rates. Our hybrid multi-organ implementation suggests that effector CD4+ T cells can rescue the system from a persistent infection and lead to clearance once a granuloma is fully formed. This could be effective as an immunotherapy strategy for latently infected individuals.

Dataset Information

1D-3D hybrid modeling-from multi-compartment models to full resolution models in space and time.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets