Showing papers in "arXiv: Neurons and Cognition in 2006"

Ising models for networks of real neurons

Abstract: Ising models with pairwise interactions are the least structured, or maximum-entropy, probability distributions that exactly reproduce measured pairwise correlations between spins. Here we use this equivalence to construct Ising models that describe the correlated spiking activity of populations of 40 neurons in the retina, and show that pairwise interactions account for observed higher-order correlations. By first finding a representative ensemble for observed networks we can create synthetic networks of 120 neurons, and find that with increasing size the networks operate closer to a critical point and start exhibiting collective behaviors reminiscent of spin glasses.

The astrocyte as a gatekeeper of synaptic information transfer

Abstract: We present a simple biophysical model for the coupling between synaptic transmission and the local calcium concentration on an enveloping astrocytic domain. This interaction enables the astrocyte to modulate the information flow from presynaptic to postsynaptic cells in a manner dependent on previous activity at this and other nearby synapses. Our model suggests a novel, testable hypothesis for the spike timing statistics measured for rapidly-firing cells in culture experiments.

Conservation rules, their breakdown, and optimality in Caenorhabditis sinusoidal locomotion

Abstract: Undulatory locomotion is common to nematodes as well as to limbless vertebrates, but its control is not understood in spite of the identification of hundred of genes involved in Caenorhabditis elegans locomotion. To reveal the mechanisms of nematode undulatory locomotion, we quantitatively analyzed the movement of C. elegans with genetic perturbations to neurons, muscles, and skeleton (cuticle). We also compared locomotion of different Caenorhabditis species. We constructed a theoretical model that combines mechanics and biophysics, and that is constrained by the observations of propulsion and muscular velocities, as well as wavelength and amplitude of undulations. We find that normalized wavelength is a conserved quantity among wild-type C. elegans individuals, across mutants, and across different species. The velocity of forward propulsion scales linearly with the velocity of the muscular wave and the corresponding slope is also a conserved quantity and almost optimal; the exceptions are in some mutants affecting cuticle structure. In theoretical terms, the optimality of the slope is equivalent to the exact balance between muscular and visco-elastic body reaction bending moments. We find that the amplitude and frequency of undulations are inversely correlated and provide a theoretical explanation for this fact. These experimental results are valid both for young adults and for all larval stages of wild type C. elegans. In particular, during development, the amplitude scales linearly with the wavelength, consistent with our theory. We also investigated the influence of substrate firmness on motion parameters, and found that it does not affect the above invariants. In general, our biomechanical model can explain the observed robustness of the mechanisms controlling nematode undulatory locomotion.

Glial activation in white matter following ischemia in the neonatal P7 rat brain

Abstract: This study examines cell death and proliferation in the white matter after neonatal stroke. In post-natal day 7 injured rat, there was a marked reduction in myelin basic protein (MBP) immunostaining mainly corresponding to numerous pyknotic immature oligodendrocytes and TUNEL-positive astrocytes in the ipsilateral external capsule. In contrast, a substantial restoration of MBP, as indicated by the MBP ratio of left-toright, occurred in the cingulum at 48 (1.27 +- 0.12) and 72 (1.30 +- 0.18, p<0.05) hours of recovery as compared to age-matched controls (1.03 +- 0.14). Ki-67 immunostaining revealed a first peak of newly-generated cells in the dorsolateral hippocampal subventricular zone and cingulum at 72 hours after reperfusion. Double immunofluorescence revealed that most of the Ki-67-positive cells were astrocytes at 48 hours and NG2 pre-oligodendrocytes at 72 hours of recovery. Microglia infiltration occurs over several days in the cingulum and a huge quantity of macrophages reached the subcortical white matter where they engulfed immature oligodendrocytes. The overall results suggest that the persistent activation of microglia involves a chronic component of immunoinflammation, which overwhelms repair processes and contributes to cystic growth in the developing brain.

Coding static natural images using spiking event times: do neurons cooperate?

Abstract: To understand possible strategies of temporal spike coding in the central nervous system, we study functional neuromimetic models of visual processing for static images. We will first present the retinal model which was introduced by Van Rullen and Thorpe and which represents the multiscale contrast values of the image using an orthonormal wavelet transform. These analog values activate a set of spiking neurons which each fire once to produce an asynchronous wave of spikes. According to this model, the image may be progressively reconstructed from this spike wave thanks to regularities in the statistics of the coefficients determined with natural images. Here, we study mathematically how the quality of information transmission carried by this temporal representation varies over time. In particular, we study how these regularities can be used to optimize information transmission by using a form of temporal cooperation of neurons to code analog values. The original model used wavelet transforms that are close to orthogonal. However, the selectivity of realistic neurons overlap, and we propose an extension of the previous model by adding a spatial cooperation between filters. This model extends the previous scheme for arbitrary--and possibly nonorthogonal--representations of features in the images. In particular, we compared the performance of increasingly over-complete representations in the retina. Results show that this algorithm provides an efficient spike coding strategy for low-level visual processing which may adapt to the complexity of the visual input.

Simulation of networks of spiking neurons: A review of tools and strategies

Abstract: We review different aspects of the simulation of spiking neural networks. We start by reviewing the different types of simulation strategies and algorithms that are currently implemented. We next review the precision of those simulation strategies, in particular in cases where plasticity depends on the exact timing of the spikes. We overview different simulators and simulation environments presently available (restricted to those freely available, open source and documented). For each simulation tool, its advantages and pitfalls are reviewed, with an aim to allow the reader to identify which simulator is appropriate for a given task. Finally, we provide a series of benchmark simulations of different types of networks of spiking neurons, including Hodgkin-Huxley type, integrate-and-fire models, interacting with current-based or conductance-based synapses, using clock-driven or event-driven integration strategies. The same set of models are implemented on the different simulators, and the codes are made available. The ultimate goal of this review is to provide a resource to facilitate identifying the appropriate integration strategy and simulation tool to use for a given modeling problem related to spiking neural networks.

Single neuron computation: from dynamical system to feature detector

Abstract: White noise methods are a powerful tool for characterizing the computation performed by neural systems. These methods allow one to identify the feature or features that a neural system extracts from a complex input, and to determine how these features are combined to drive the system's spiking response. These methods have also been applied to characterize the input/output relations of single neurons driven by synaptic inputs, simulated by direct current injection. To interpret the results of white noise analysis of single neurons, we would like to understand how the obtained feature space of a single neuron maps onto the biophysical properties of the membrane, in particular the dynamics of ion channels. Here, through analysis of a simple dynamical model neuron, we draw explicit connections between the output of a white noise analysis and the underlying dynamical system. We find that under certain assumptions, the form of the relevant features is well defined by the parameters of the dynamical system. Further, we show that under some conditions, the feature space is spanned by the spike-triggered average and its successive order time derivatives.

Competition between synaptic depression and facilitation in attractor neural networks

Abstract: We study the effect of competition between short-term synaptic depression and facilitation on the dynamical properties of attractor neural networks, using Monte Carlo simulation and a mean field analysis. Depending on the balance between depression, facilitation and the noise, the network displays different behaviours, including associative memory and switching of the activity between different attractors. We conclude that synaptic facilitation enhances the attractor instability in a way that (i) intensifies the system adaptability to external stimuli, which is in agreement with experiments, and (ii) favours the retrieval of information with less error during short--time intervals.

Ecosystems in the Mind: Fuzzy Cognitive Maps of the Kizilirmak Delta Wetlands in Turkey

Abstract: Sustainability of ecosystems & ecosystem management are increasingly accepted societal goals Can conservation programs be improved by incorporating local people's understanding of ecosystems? The Kizilirmak Delta is one of Turkey's most important wetland complexes It is also one of the most productive agricultural deltas in Turkey We obtained 31 cognitive models of the social & ecological system These models were converted to adjacency matrices, analyzed using graph theoretical methods, & augmented into social cognitive maps Causal "What-if" scenarios were run to determine the trajectory of the ecosystem based on models defined by stakeholders Villagers had significantly larger numbers of variables, more complex maps, a broader understanding of all the variables that affect the Kizilirmak Delta, & mentioned more variables that control the ecosystem than did NGO and government officials Villagers adapting to changing ecological & social conditions actively changed & challenged conditions through the political process Villagers were faced with many important forcing functions that they could not control Most of the variables defined by villagers were related to agriculture and animal husbandry Conservation policies & ecosystem management must encompass larger environmental issues & villagers' cognitive maps must be reconciled with that of NGOs & government officials Cognitive maps can serve as a basis for discussion when policies & management options are formulated A villager-centered cognitive mapping approach is not only necessary because villagers resist conservation projects, or because top down projects that do not take local knowledge systems into account fail, but because it is the ethical and responsible way of doing ecosystem management

Searching for memories, Sudoku, implicit check-bits, and the iterative use of not-always-correct rapid neural computation

Abstract: The algorithms that simple feedback neural circuits representing a brain area can rapidly carry out are often adequate to solve only easy problems, and for more difficult problems can return incorrect answers. A new excitatory-inhibitory circuit model of associative memory displays the common human problem of failing to rapidly find a memory when only a small clue is present. The memory model and a related computational network for solving Sudoku puzzles produce answers that contain implicit check-bits in the representation of information across neurons, allowing a rapid evaluation of whether the putative answer is correct or incorrect through a computation related to visual 'pop-out'. This fact may account for our strong psychological feeling of right or wrong when we retrieve a nominal memory from a minimal clue. This information allows more difficult computations or memory retrievals to be done in a serial fashion by using the fast but limited capabilities of a computational module multiple times. The mathematics of the excitatory-inhibitory circuits for associative memory and for Sudoku, both of which are understood in terms of 'energy' or Lyapunov functions, is described in detail.

Discovering Functional Communities in Dynamical Networks

Abstract: Many networks are important because they are substrates for dynamical systems, and their pattern of functional connectivity can itself be dynamic -- they can functionally reorganize, even if their underlying anatomical structure remains fixed. However, the recent rapid progress in discovering the community structure of networks has overwhelmingly focused on that constant anatomical connectivity. In this paper, we lay out the problem of discovering_functional communities_, and describe an approach to doing so. This method combines recent work on measuring information sharing across stochastic networks with an existing and successful community-discovery algorithm for weighted networks. We illustrate it with an application to a large biophysical model of the transition from beta to gamma rhythms in the hippocampus.

Calculating event-triggered average synaptic conductances from the membrane potential

Abstract: The optimal patterns of synaptic conductances for spike generation in central neurons is a subject of considerable interest. Ideally, such conductance time courses should be extracted from membrane potential (Vm) activity, but this is difficult because the nonlinear contribution of conductances to the Vm renders their estimation from the membrane equation extremely sensitive. We outline here a solution to this problem based on a discretization of the time axis. This procedure can extract the time course of excitatory and inhibitory conductances solely from the analysis of Vm activity. We test this method by calculating spike-triggered averages of synaptic conductances using numerical simulations of the integrate-and-fire model subject to colored conductance noise. The procedure was also tested successfully in biological cortical neurons using conductance noise injected with dynamic-clamp. This method should allow the extraction of synaptic conductances from Vm recordings in vivo.

Capacitance fluctuations causing channel noise reduction in stochastic Hodgkin-Huxley systems

Abstract: Voltage-dependent ion channels determine the electric properties of axonal cell membranes. They not only allow the passage of ions through the cell membrane but also contribute to an additional charging of the cell membrane resulting in the so-called capacitance loading. The switching of the channel gates between an open and a closed configuration is intrinsically related to the movement of gating charge within the cell membrane. At the beginning of an action potential the transient gating current is opposite to the direction of the current of sodium ions through the membrane. Therefore, the excitability is expected to become reduced due to the influence of a gating current. Our stochastic Hodgkin-Huxley like modeling takes into account both the channel noise -- i.e. the fluctuations of the number of open ion channels -- and the capacitance fluctuations that result from the dynamics of the gating charge. We investigate the spiking dynamics of membrane patches of variable size and analyze the statistics of the spontaneous spiking. As a main result, we find that the gating currents yield a drastic reduction of the spontaneous spiking rate for sufficiently large ion channel clusters. Consequently, this demonstrates a prominent mechanism for channel noise reduction.

The brain near the edge

Abstract: When viewed at a certain coarse grain, the brain seems a relatively small dynamical system composed by a few dozen interacting areas, performing a number of stereotypical behaviors. It is known that, even relatively small dynamical systems can reliably generate robust and flexible behavior if they are possed near a second order phase transition, because of the abundance of metastable states at the critical point. The approach pursued here assumes that some of the most fundamental properties of the functioning brain are possible because it is spontaneously possed at the border of such instability. In this notes we review the motivation, the arguments and recent results as well as the implications of this view of the functioning brain.

Toward an adequate mathematical model of mental space: conscious/unconscious dynamics on $m$-adic trees

Abstract: We try to perform geometrization of cognitive science and psychology by representing information states of cognitive systems by points of {\it mental space} given by a hierarchic $m$-adic tree. Associations are represented by balls and ideas by collections of balls. We consider dynamics of ideas based on lifting of dynamics of mental points. We apply our dynamical model for modeling of flows of unconscious and conscious information in the human brain. In series of models, Models 1-3, we consider cognitive systems with increasing complexity of psychological behavior determined by structure of flows of associations and ideas.

Synaptic plasticity of Inhibitory synapse promote synchrony in inhibitory network in presence of heterogeneity and noise

Abstract: Recently spike timing dependent plasticity was observed in inhibitory synapse in the layer II of entorhinal cortex. The rule provides an interesting zero in the region of $\Delta t=t_{post}-t_{pre}=0$ and in addition the dynamic range of the rule lie in gamma frequency band. We propose a robust mechanism based on this observed synaptic plasticity rule for inhibitory synapses for two mutually coupled interneurons to phase lock in synchrony in the presence of intrisic heterogeneity in firing. We study the stability of the phase locked solution by defining a map for spike times dependent on the phase response curve for the coupled neurons. Finally we present results on robustness of synchronization in the presence of noise.

Frequency and phase synchronization of two coupled neurons with channel noise

Abstract: We study the frequency and phase synchronization in two coupled identical and nonidentical neurons with channel noise. The occupation number method is used to model the neurons in the context of stochastic Hodgkin-Huxley model in which the strength of of channel noise is represented by ion channel cluster size of the initiation region of neuron. It is shown that frequency synchronization only was achieved at arbitrary value of couple strength as long as two neurons' channel cluster sizes are the same. We also show that the relative phase of neurons can display profuse dynamic behavior under the combined action of coupling and channel noise. Both qualitative and quantitative descriptions are applied to describe the transitions between those behaviors. Relevance of our findings to controlling neural synchronization experimentally is discussed.

From systems biology to dynamical neuropharmacology: proposal for a new methodology

Abstract: The concepts and methods of Systems Biology are being extended to neuropharmacology, to test and design drugs against neurological and psychiatric disorders. Computational modeling by integrating compartmental neural modeling technique and detailed kinetic description of pharmacological modulation of transmitter-receptor interaction is offered as a method to test the electrophysiological and behavioral effects of putative drugs. Even more, an inverse method is suggested as a method for controlling a neural system to realize a prescribed temporal pattern. In particular, as an application of the proposed new methodology a computational platform is offered to analyze the generation and pharmacological modulation of theta rhythm related to anxiety is analyzed here in more detail.

On the use of analytic expressions for the voltage distribution to analyze intracellular recordings

Abstract: Different analytic expressions for the membrane potential distribution of membranes subject to synaptic noise have been proposed, and can be very helpful to analyze experimental data. However, all of these expressions are either approximations or limit cases, and it is not clear how they compare, and which expression should be used in a given situation. In this note, we provide a comparison of the different approximations available, with an aim to delineate which expression is most suitable for analyzing experimental data.

Neuronal networks and controlled symmetries, a generic framework

Abstract: The extraordinary computational power of the brain may be related in part to the fact that each of the smaller neural networks that compose it can behave transiently in many different ways, depending on its inputs. Mathematically, input continuity helps to show how a large network, constructed recursively from smaller blocks, can exhibit robust specific properties according to its input. By extending earlier work on synchrony and symmetry, we exploit input continuity of contracting systems to ensure robust control of diverse spatial and spatio-temporal symmetries of the output signal in such a network.

Synergy from Silence in a Combinatorial Neural Code

Abstract: The manner in which groups of neurons represent events in the external world is fundamental to neuroscience. Here, we analyze the population code of the retina during naturalistic stimulation and show that the information conveyed by specific multi-neuronal firing patterns can be very different from the sum of the pattern's parts. Synchronous spikes convey more information than either of the participating cells, but almost always less than their sum, making them redundant coding symbols. Surprisingly, patterns of spiking and silence are mostly synergistic - carrying information that is accessible only by observing the whole pattern of activity, rather than its components, and signifying unique features in the stimulus. These results demonstrate that the retina uses a combinatorial code and that the brain can benefit significantly from recognizing multi-neuronal firing patterns.

Cerebellar Purkinje Cell Loss in Heterozygous Rora+/-Mice: A Longitudinal Study

Abstract: The staggerer (sg/sg) mutation is a spontaneous deletion in the Rora gene that prevents the translation of the ligand-binding domain (LBD), leading to the loss of ROR\alpha activity. The homozygous Rorasg/sg mutant mouse, whose most obvious phenotype is ataxia associated with cerebellar degeneration, also displays a variety of other phenotypes. The heterozygous Rora+/sg is able to develop a cerebellum which is qualitatively normal but with advancing age suffers a significant loss of cerebellar neuronal cells. A truncated protein synthesized by the mutated allele may play a role, both in Rorasg/sg and Rora+/sg. To determine the effects during life span of true haplo-insufficiency of the ROR\alpha protein, derived from the invalidation of the gene, we compared the evolution of Purkinje cell numbers in heterozygous Rora knock-out males (Rora+/-) and in their wildtype counterparts from 1 to 24 months of age. We also compared the evolution of Purkinje cell numbers in Rora+/- and Rora+/sg males from 1 to 9 months. The main finding is that in Rora+/- mice, when only a half dose of protein is synthesized, the deficit was already established at 1 month and did not change during life span....

Frequency decomposition of conditional Granger causality and application to multivariate neural field potential data

Abstract: It is often useful in multivariate time series analysis to determine statistical causal relations between different time series. Granger causality is a fundamental measure for this purpose. Yet the traditional pairwise approach to Granger causality analysis may not clearly distinguish between direct causal influences from one time series to another and indirect ones acting through a third time series. In order to differentiate direct from indirect Granger causality, a conditional Granger causality measure in the frequency domain is derived based on a partition matrix technique. Simulations and an application to neural field potential time series are demonstrated to validate the method.

Triangular lattice neurons may implement an advanced numeral system to precisely encode rat position over large ranges

Abstract: We argue by observation of the neural data that neurons in area dMEC of rats, which fire whenever the rat is on any vertex of a regular triangular lattice that tiles 2-d space, may be using an advanced numeral system to reversibly encode rat position. We interpret measured dMEC properties within the framework of a residue number system (RNS), and describe how RNS encoding -- which breaks the non-periodic variable of rat position into a set of narrowly distributed periodic variables -- allows a small set of cells to compactly represent and efficiently update rat position with high resolution over a large range. We show that the uniquely useful properties of RNS encoding still hold when the encoded and encoding quantities are relaxed to be real numbers with built-in uncertainties, and provide a numerical and functional estimate of the range and resolution of rat positions that can be uniquely encoded in dMEC. The use of a compact, `arithmetic-friendly' numeral system to encode a metric variable, as we propose is happening in dMEC, is qualitatively different from all previously identified examples of coding in the brain. We discuss the numerous neurobiological implications and predictions of our hypothesis.

Neural integrator - a sandpile model

Abstract: We investigated a model for the neural integrator based on hysteretic units connected by positive feedback. Hysteresis is assumed to emerge from the intrinsic properties of the cells. We consider the recurrent networks containing either bistable or multistable neurons. We apply our analysis to the oculomotor velocity-to-position neural integrator that calculates the eye positions from the inputs that carry information about eye angular velocity. Using the analysis of the system in the parameter space we show the following. The direction of hysteresis in the neuronal response may be reversed for the system with recurrent connections compared to the case of unconnected neurons. Thus, for the NMDA receptor based bistability the firing rates after ON saccades may be higher than after OFF saccades for the same eye position. We suggest that this is an emergent property due to the presence of global recurrent feedback. The reversal of hysteresis occurs only when the size of hysteresis differs from neuron to neuron. We also relate the macroscopic leak time-constant of the integrator to the rate of microscopic spontaneous noise-driven transitions in the hysteretic units. Finally, we argue that the presence of neurons with small hysteresis may remove the threshold for integration.

Neural coding of a natural stimulus ensemble: Uncovering information at sub-millisecond resolution

Abstract: Our knowledge of the sensory world is encoded by neurons in sequences of discrete, identical pulses termed action potentials or spikes. There is persistent controversy about the extent to which the precise timing of these spikes is relevant to the function of the brain. We revisit this issue, using the motion-sensitive neurons of the fly visual system as a test case. New experimental methods allow us to deliver more nearly natural visual stimuli, comparable to those which flies encounter in free, acrobatic flight, and new mathematical methods allow us to draw more reliable conclusions about the information content of neural responses even when the set of possible responses is very large. We find that significant amounts of visual information are represented by details of the spike train at millisecond and sub-millisecond precision, even though the sensory input has a correlation time of ~60ms; different patterns of spike timing represent distinct motion trajectories, and the absolute timing of spikes points to particular features of these trajectories with high precision. Under these naturalistic conditions, the system continues to transmit more information at higher photon flux, even though individual photoreceptors are counting more than one million photons per second, and removes redundancy in the stimulus to generate a more efficient neural code.

The role of synaptic facilitation in coincidence spike detection

Abstract: Using a realistic model of activity dependent dynamical synapses and a standard integrate and fire neuron model we study, both analytically and numerically, the conditions in which a postsynaptic neuron efficiently detects temporal coincidences of spikes arriving at certain frequency from N different afferents. We extend a previous work that only considers synaptic depression as the most important mechanism in the transmission of information through synapses, to a more general situation including also synaptic facilitation. Our study shows that: 1) facilitation enhances the detection of correlated signals arriving from a subset of presynaptic excitatory neurons, with different degrees of correlation among this subset, and 2) the presence of facilitation allows for a better detection of firing rate changes. Finally, we also observed that facilitation determines the existence of an optimal input frequency which allows the best performance for a wide (maximum) range of the neuron firing threshold. This optimal frequency can be controlled by means of facilitation parameters.

Spatial embedding of neuronal trees modeled by diffusive growth

Abstract: The relative importance of the intrinsic and extrinsic factors determining the variety of geometric shapes exhibited by dendritic trees remains unclear. This question was addressed by developing a model of the growth of dendritic trees based on diffusion-limited aggregation process. The model reproduces diverse neuronal shapes (i.e., granule cells, Purkinje cells, the basal and apical dendrites of pyramidal cells, and the axonal trees of interneurons) by changing only the size of the growth area, the time span of pruning, and the spatial concentration of 'neurotrophic particles'. Moreover, the presented model shows how competition between neurons can affect the shape of the dendritic trees. The model reveals that the creation of complex (but reproducible) dendrite-like trees does not require precise guidance or an intrinsic plan of the dendrite geometry. Instead, basic environmental factors and the simple rules of diffusive growth adequately account for the spatial embedding of different types of dendrites observed in the cortex. An example demonstrating the broad applicability of the algorithm to model diverse types of tree structures is also presented. Key words: Diffusion-limited aggregation; Neuronal morphology; Dendrites; Growth model; tree; dendritic geometry.

Using Adaptive Dynamic Programming to Understand and Replicate Brain Intelligence: the Next Level Design

Abstract: Since the 1960s I proposed that we could understand and replicate the highest level of intelligence seen in the brain, by building ever more capable and general systems for adaptive dynamic programming (ADP), which is like reinforcement learning but based on approximating the Bellman equation and allowing the controller to know its utility function. Growing empirical evidence on the brain supports this approach. Adaptive critic systems now meet tough engineering challenges and provide a kind of first-generation model of the brain. Lewis, Prokhorov and myself have early second-generation work. Mammal brains possess three core capabilities, creativity/imagination and ways to manage spatial and temporal complexity, even beyond the second generation. This paper reviews previous progress, and describes new tools and approaches to overcome the spatial complexity gap.

Structured thalamocortical connectivity revealed by random walks on complex networks

Abstract: The segregated regions of the mammalian cerebral cortex and thalamus form an extensive and complex network, whose structure and function are still only incompletely understood. The present article describes an application of the concepts of complex networks and random walks that allows the identification of non-random, highly structured features of thalamocortical connections, and their potential effects on dynamic interactions between cortical areas in the cat brain. Utilizing large-scale anatomical data sets of this thalamocortical system, we investigate uniform random walks in such a network by considering the steady state eigenvector of the respective stochastic matrix. It is shown that thalamocortical connections are organized in such a way as to guarantee strong correlation between the outdegree and occupancy rate (a stochastic measure potentially related to activation) of each cortical area. Possible organizational principles underlying this effect are identified and discussed. The relationship between the topology of structural connection patterns and cortical dynamics currently represents a significant challenge to brain theory. Numerous neuroanatomical studies have revealed that the pathways of the mammalian thalamocortical system exhibit specific patterns ranging in scale from interconnections linking whole brain regions to intra-areal patterns of connections between cell populations or individual cortical neurons. Comprehensive descriptions of large-scale anatomical patterns of cortical connectivity have been collated for several mammalian species (e.g. [1, 2]). Quantitative analysis has revealed that these patterns are neither completely regular nor completely random [3], but exhibit specific attributes, such as short path lengths combined with high clustering [4, 5], short wiring and distinctive hierarchical features [6]. The current article proposes a novel macroscopic approach to studying neuronal dynamics in cortical maps based on concepts from random walks (e.g. [7, 8]) and complex networks (e.g. [9, 10, 11]). Specifically, we examine the interactions between cortex and thalamus in the mammalian central nervous system. The cerebral cortex consists of a network of interconnected functionally specialized regions, each of which maintains reciprocal connections with a specific set of thalamic nuclei. These nuclei relay cortical inputs and outputs, and are thought to play a role in coordinating cortical interactions. The thalamocortical architecture is represented in terms of a complex network, obtained by assigning a node to each cortical region and thalamic nucleus (nodes) and