Linking brain behavior to underlying cellular mechanisms via large-scale brain modeling and simulation

To understand brain behaviors, it is important to directly associate the network level activities to the underlying biophysical mechanisms, which require large-scale simulations with biophysically realistic neural models like Hodgkin-Huxley models. However, when simulations are conducted on models with sufficient biophysical details, great challenges arise from limited computer power, thereby restricting most existing computational works with biophysical models only to small-scale networks. On the other hand, with the emergence of powerful computing platforms, many recent works are geared to performing large-scale simulations with simple spiking models. However, the applicability of those works is limited by the nature of the underlying phenomenological model. To bridge the gap, an intermediate step is taken to construct a scalable brain model with sufficient biophysical details. In this work, great efforts are devoted to taking into account not only local cortical microcircuits but also the global brain architecture, and efficient techniques are proposed and adopted to address the associated computational challenges in simulation of networks of such complexity. With the customized simulator developed, we are able to simulate the brain model to generate not only sleep spindle and delta waves but also the spike-and-wave pattern of absence seizures, and directly link those behaviors to underlying biophysical mechanism. Those initial results are interesting because they show the possibility to determine underlying causes of diseases by simulating the biologically realistic brain model. With further development, the work is geared to assisting the clinicians in selecting the optimal treatment on an individual basis in the future.

[1]  Boyuan Yan,et al.  An Integrative View of Mechanisms Underlying Generalized Spike-and-Wave Epileptic Seizures and Its Implication on Optimal Therapeutic Treatments , 2011, PloS one.

[2]  V. A. Saletore A Distributed and Adaptive Dynamic Load Balancing Scheme for Parallel Processing of Medium-Grain Tasks , 1990, Proceedings of the Fifth Distributed Memory Computing Conference, 1990..

[3]  A. Destexhe Spike-and-Wave Oscillations Based on the Properties of GABAB Receptors , 1998, The Journal of Neuroscience.

[4]  G. Deco,et al.  Emerging concepts for the dynamical organization of resting-state activity in the brain , 2010, Nature Reviews Neuroscience.

[5]  Marc-Oliver Gewaltig,et al.  Efficient Parallel Simulation of Large-Scale Neuronal Networks on Clusters of Multiprocessor Computers , 2007, Euro-Par.

[6]  A. Depaulis,et al.  Involvement of intrathalamic GABA b neurotransmission in the control of absence seizures in the rat , 1992, Neuroscience.

[7]  Dharmendra S. Modha,et al.  The cat is out of the bag: cortical simulations with 109 neurons, 1013 synapses , 2009, Proceedings of the Conference on High Performance Computing Networking, Storage and Analysis.

[8]  Donald O. Walter,et al.  Mass action in the nervous system , 1975 .

[9]  Eugene M. Izhikevich,et al.  Which model to use for cortical spiking neurons? , 2004, IEEE Transactions on Neural Networks.

[10]  Stiliyan Kalitzin,et al.  Dynamical diseases of brain systems: different routes to epileptic seizures , 2003, IEEE Transactions on Biomedical Engineering.

[11]  The selective GABAB antagonist CGP-35348 blocks spike-wave bursts in the cholesterol synthesis rat absence epilepsy model , 1996, Brain Research.

[12]  James M. Bower,et al.  The Book of GENESIS , 1994, Springer New York.

[13]  John R. Terry,et al.  On the genesis of spike-wave oscillations in a mean-field model of human thalamic and corticothalamic dynamics , 2006 .

[14]  Henry Markram,et al.  Minimal Hodgkin–Huxley type models for different classes of cortical and thalamic neurons , 2008, Biological Cybernetics.

[15]  Stephen Coombes,et al.  Large-scale neural dynamics: Simple and complex , 2010, NeuroImage.

[16]  M Steriade,et al.  Interneuronal epileptic discharges related to spike-and-wave cortical seizures in behaving monkeys. , 1974, Electroencephalography and clinical neurophysiology.

[17]  S. Amari,et al.  Characteristics of Random Nets of Analog Neuron-Like Elements , 1972, IEEE Trans. Syst. Man Cybern..

[18]  J. Martinerie,et al.  The brainweb: Phase synchronization and large-scale integration , 2001, Nature Reviews Neuroscience.

[19]  Maxim Bazhenov,et al.  Thalamocortical oscillations , 2006, Scholarpedia.

[20]  J. Parra,et al.  Epileptic Transitions: Model Predictions and Experimental Validation , 2005, Journal of clinical neurophysiology : official publication of the American Electroencephalographic Society.

[21]  N. Spruston Pyramidal neurons: dendritic structure and synaptic integration , 2008, Nature Reviews Neuroscience.

[22]  D. McCormick,et al.  Noradrenergic and serotonergic modulation of a hyperpolarization‐activated cation current in thalamic relay neurones. , 1990, The Journal of physiology.

[23]  Fabrice Wendling,et al.  Computational modeling of high-frequency oscillations at the onset of neocortical partial seizures: From ‘altered structure’ to ‘dysfunction’ , 2010, NeuroImage.

[24]  Peng Li,et al.  Simulation of large neuronal networks with biophysically accurate models on graphics processors , 2011, The 2011 International Joint Conference on Neural Networks.

[25]  P. Gloor,et al.  Transition from spindles to generalized spike and wave discharges in the cat: Simultaneous single-cell recordings in cortex and thalamus , 1984, Experimental Neurology.

[26]  Sartaj Sahni,et al.  Programming a hypercube multicomputer , 1988, IEEE Software.

[27]  Arndt Bode Load balancing in distributed memory multiprocessors , 1991, [1991] Proceedings, Advanced Computer Technology, Reliable Systems and Applications.

[28]  Karl J. Friston,et al.  The Dynamic Brain: From Spiking Neurons to Neural Masses and Cortical Fields , 2008, PLoS Comput. Biol..

[29]  Andrew Zalesky,et al.  A DTI-Derived Measure of Cortico-Cortical Connectivity , 2009, IEEE Transactions on Medical Imaging.

[30]  B. Knight The Relationship between the Firing Rate of a Single Neuron and the Level of Activity in a Population of Neurons , 1972, The Journal of general physiology.

[31]  W. A. Wilson,et al.  The role of GABAB receptor activation in absence seizures of lethargic (lh/lh) mice. , 1992, Science.

[32]  Duane Q. Nykamp,et al.  A Population Density Approach That Facilitates Large-Scale Modeling of Neural Networks: Analysis and an Application to Orientation Tuning , 2004, Journal of Computational Neuroscience.

[33]  T J Sejnowski,et al.  In vivo, in vitro, and computational analysis of dendritic calcium currents in thalamic reticular neurons , 1996, The Journal of neuroscience : the official journal of the Society for Neuroscience.

[34]  P. Robinson,et al.  Dynamics of large-scale brain activity in normal arousal states and epileptic seizures. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[35]  Ioannis G. Kevrekidis,et al.  Projective Methods for Stiff Differential Equations: Problems with Gaps in Their Eigenvalue Spectrum , 2002, SIAM J. Sci. Comput..

[36]  Germán Mato,et al.  Synchrony in Excitatory Neural Networks , 1995, Neural Computation.

[37]  Richard M. Leahy,et al.  Brainstorm: A User-Friendly Application for MEG/EEG Analysis , 2011, Comput. Intell. Neurosci..

[38]  O. Snead Evidence for GABAB-mediated mechanisms in experimental generalized absence seizures. , 1992, European journal of pharmacology.

[39]  A. Hodgkin,et al.  A quantitative description of membrane current and its application to conduction and excitation in nerve , 1990 .

[40]  George Cybenko,et al.  Dynamic Load Balancing for Distributed Memory Multiprocessors , 1989, J. Parallel Distributed Comput..

[41]  A. Peters,et al.  Numerical relationships between geniculocortical afferents and pyramidal cell modules in cat primary visual cortex. , 1993, Cerebral cortex.

[42]  C. Cepeda,et al.  Morphological and electrophysiological characterization of abnormal cell types in pediatric cortical dysplasia , 2003, Journal of neuroscience research.

[43]  M L Hines,et al.  Neuron: A Tool for Neuroscientists , 2001, The Neuroscientist : a review journal bringing neurobiology, neurology and psychiatry.

[45]  Bruce W. Knight,et al.  Dynamics of Encoding in a Population of Neurons , 1972, The Journal of general physiology.

[46]  John R. Terry,et al.  Onset of polyspike complexes in a mean-field model of human electroencephalography and its application to absence epilepsy , 2009, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[47]  John R. Terry,et al.  A unifying explanation of primary generalized seizures through nonlinear brain modeling and bifurcation analysis. , 2006, Cerebral cortex.

[48]  J. Cowan,et al.  Excitatory and inhibitory interactions in localized populations of model neurons. , 1972, Biophysical journal.

[49]  Frédérique Amor,et al.  Cortical local and long-range synchronization interplay in human absence seizure initiation , 2009, NeuroImage.

[50]  Ioannis G. Kevrekidis,et al.  Telescopic projective methods for parabolic differential equations , 2003 .

[51]  J. Hablitz,et al.  Abnormal pyramidal cell morphology and HCN channel expression in cortical dysplasia , 2010, Epilepsia.

[52]  Christophe Bernard,et al.  Synaptic integration of NMDA and non-NMDA receptors in large neuronal network models solved by means of differential equations , 2004, Biological Cybernetics.

[53]  M. Deschenes,et al.  Abolition of spindle oscillations in thalamic neurons disconnected from nucleus reticularis thalami. , 1985, Journal of neurophysiology.

[54]  C. J. Honeya,et al.  Predicting human resting-state functional connectivity from structural connectivity , 2009 .

[55]  B. Connors,et al.  Intrinsic firing patterns of diverse neocortical neurons , 1990, Trends in Neurosciences.

[56]  M. Dichter,et al.  Cellular mechanisms of epilepsy: a status report. , 1987, Science.

[57]  Jacques Duysens,et al.  Thalamic multiple-unit activity underlying spike-wave discharges in anesthetized rats , 1993, Brain Research.

[58]  Fabrice Wendling,et al.  Relevance of nonlinear lumped-parameter models in the analysis of depth-EEG epileptic signals , 2000, Biological Cybernetics.

[59]  T. Sejnowski,et al.  G protein activation kinetics and spillover of gamma-aminobutyric acid may account for differences between inhibitory responses in the hippocampus and thalamus. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[60]  J Gotman,et al.  An analysis of penicillin-induced generalized spike and wave discharges using simultaneous recordings of cortical and thalamic single neurons. , 1983, Journal of neurophysiology.

[61]  G. Buzsáki,et al.  Petit mal epilepsy and parkinsonian tremor: Hypothesis of a common pacemaker , 1990, Neuroscience.

[62]  G. Edelman,et al.  Large-scale model of mammalian thalamocortical systems , 2008, Proceedings of the National Academy of Sciences.

[63]  S. Amari Homogeneous nets of neuron-like elements , 1975, Biological Cybernetics.

[64]  Gunnar Rätsch,et al.  Support Vector Machines and Kernels for Computational Biology , 2008, PLoS Comput. Biol..

[65]  D. Modha,et al.  Network architecture of the long-distance pathways in the macaque brain , 2010, Proceedings of the National Academy of Sciences.

[66]  D. Liley,et al.  Understanding the Transition to Seizure by Modeling the Epileptiform Activity of General Anesthetic Agents , 2005, Journal of clinical neurophysiology : official publication of the American Electroencephalographic Society.

[67]  A. Daffertshofer,et al.  Multivariate Ornstein-Uhlenbeck processes with mean-field dependent coefficients: application to postural sway. , 2000, Physical review. E, Statistical, nonlinear, and soft matter physics.

[68]  J. Bellanger,et al.  Interictal to Ictal Transition in Human Temporal Lobe Epilepsy: Insights From a Computational Model of Intracerebral EEG , 2005, Journal of clinical neurophysiology : official publication of the American Electroencephalographic Society.

[69]  D. Pollen INTRACELLULAR STUDIES OF CORTICAL NEURONS DURING THALAMIC INDUCED WAVE AND SPIKE. , 1964, Electroencephalography and clinical neurophysiology.

[70]  R. Douglas,et al.  A Quantitative Map of the Circuit of Cat Primary Visual Cortex , 2004, The Journal of Neuroscience.

[71]  G. Edelman,et al.  Spike-timing dynamics of neuronal groups. , 2004, Cerebral cortex.

[72]  Ingo Bojak,et al.  Realistic mean field forward predictions for the integration of co-registered EEG/fMRI , 2009, BMC Neuroscience.

[73]  G. Kohring,et al.  COMMENT: Convergence time and finite size effects in neural networks , 1990 .

[74]  Karl J. Friston,et al.  Population dynamics under the Laplace assumption , 2009, NeuroImage.

[75]  Olaf Sporns,et al.  Network structure of cerebral cortex shapes functional connectivity on multiple time scales , 2007, Proceedings of the National Academy of Sciences.

[76]  Karl J. Friston,et al.  Dynamic causal modelling of distributed electromagnetic responses , 2009, NeuroImage.

[77]  Murray Shanahan,et al.  Accelerated simulation of spiking neural networks using GPUs , 2010, The 2010 International Joint Conference on Neural Networks (IJCNN).

[78]  G. Blasdel,et al.  Intrinsic connections of macaque striate cortex: afferent and efferent connections of lamina 4C , 1985, The Journal of neuroscience : the official journal of the Society for Neuroscience.

[79]  P. Gloor,et al.  Pathophysiology of generalized penicillin epilepsy in the cat: the role of cortical and subcortical structures. I. Systemic application of penicillin. , 1977, Electroencephalography and clinical neurophysiology.

[80]  M. Deschenes,et al.  The thalamus as a neuronal oscillator , 1984, Brain Research Reviews.

[81]  D. McCormick,et al.  Cellular mechanisms of a synchronized oscillation in the thalamus. , 1993, Science.

[82]  T. Sejnowski,et al.  Ionic mechanisms underlying synchronized oscillations and propagating waves in a model of ferret thalamic slices. , 1996, Journal of neurophysiology.

[83]  J. Bellanger,et al.  Epileptic fast activity can be explained by a model of impaired GABAergic dendritic inhibition , 2002, The European journal of neuroscience.

[84]  Arkady Pikovsky,et al.  Finite-size effects in a population of interacting oscillators , 1999 .

[85]  L F Quesney,et al.  Pathophysiology of generalized penicillin epilepsy in the cat: the role of cortical and subcortical structures. II. Topical application of penicillin to the cerebral cortex and to subcortical structures. , 1977, Electroencephalography and clinical neurophysiology.

[86]  A. Louisa,et al.  コロイド混合体における有効力 空乏引力から集積斥力へ | 文献情報 | J-GLOBAL 科学技術総合リンクセンター , 2002 .

[87]  H. Sompolinsky,et al.  Relaxational dynamics of the Edwards-Anderson model and the mean-field theory of spin-glasses , 1982 .

[88]  J. Winer,et al.  Evolution of GABAergic circuitry in the mammalian medial geniculate body. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[89]  D. Tranchina,et al.  Population density methods for large-scale modelling of neuronal networks with realistic synaptic kinetics: cutting the dimension down to size. , 2001, Network.

[90]  Moritz Grosse-Wentrup,et al.  Multisubject Learning for Common Spatial Patterns in Motor-Imagery BCI , 2011, Comput. Intell. Neurosci..

[91]  Christian Windischberger,et al.  Toward discovery science of human brain function , 2010, Proceedings of the National Academy of Sciences.

[92]  R. Traub,et al.  Neuronal Networks of the Hippocampus , 1991 .

[93]  D. Prince,et al.  A novel T-type current underlies prolonged Ca(2+)-dependent burst firing in GABAergic neurons of rat thalamic reticular nucleus , 1992, The Journal of neuroscience : the official journal of the Society for Neuroscience.

[94]  P. Bressloff Spatiotemporal dynamics of continuum neural fields , 2012 .

[95]  S. Sherman,et al.  Relative distribution of synapses in the A‐laminae of the lateral geniculate nucleus of the cat , 2000, The Journal of comparative neurology.

[96]  W. Singer,et al.  Neural Synchrony in Brain Disorders: Relevance for Cognitive Dysfunctions and Pathophysiology , 2006, Neuron.

[97]  D. McCormick,et al.  Neurotransmitter control of neocortical neuronal activity and excitability. , 1993, Cerebral cortex.

[98]  Seyed Hossein Hosseini,et al.  Analysis of a Graph Coloring Based Distributed Load Balancing Algorithm , 1990, J. Parallel Distributed Comput..

[99]  F. L. D. Silva,et al.  Dynamics of non-convulsive epileptic phenomena modeled by a bistable neuronal network , 2004, Neuroscience.

[100]  M. Kramer,et al.  Pathological pattern formation and cortical propagation of epileptic seizures , 2005, Journal of The Royal Society Interface.

[101]  Karl J. Friston,et al.  Dynamic causal modeling of evoked responses in EEG and MEG , 2006, NeuroImage.

[102]  Lawrence Sirovich,et al.  On the Simulation of Large Populations of Neurons , 2004, Journal of Computational Neuroscience.

[103]  Thom F. Oostendorp,et al.  Towards a model-based integration of co-registered electroencephalography/functional magnetic resonance imaging data with realistic neural population meshes , 2011, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[104]  G. Blasdel,et al.  Intrinsic connections of macaque striate cortex: axonal projections of cells outside lamina 4C , 1985, The Journal of neuroscience : the official journal of the Society for Neuroscience.

[105]  Viktor K. Jirsa,et al.  Noise during Rest Enables the Exploration of the Brain's Dynamic Repertoire , 2008, PLoS Comput. Biol..

[106]  D. Contreras,et al.  Mechanisms underlying the synchronizing action of corticothalamic feedback through inhibition of thalamic relay cells. , 1998, Journal of neurophysiology.

[107]  M. Steriade,et al.  Contribution of intrinsic neuronal factors in the generation of cortically driven electrographic seizures. , 2004, Journal of neurophysiology.

[108]  F. H. Lopes da Silva,et al.  Models of neuronal populations: the basic mechanisms of rhythmicity. , 1976, Progress in brain research.

[109]  R. Kötter,et al.  Connecting Mean Field Models of Neural Activity to EEG and fMRI Data , 2010, Brain Topography.