Biomedical subjects

David Cai

Publications and source records attributed to David Cai.

8 recordsLinked to original sources

PubMed2006

Orientation selectivity in visual cortex by fluctuation-controlled criticality.

Within a large-scale neuronal network model of macaque primary visual cortex, we examined how intrinsic dynamic fluctuations in synaptic currents modify the effect of strong recurrent excitation on orientation selectivity. Previously, we showed that, using a strong network inhibition countered by feedforward and recurrent excitation, the cortical model reproduced many observed properties of simple and complex cells. However, that network's complex cells were poorly selective for orientation, and increasing cortical self-excitation led to network instabilities and unrealistically high firing rates. Here, we show that a sparsity of connections in the network produces large, intrinsic fluctuations in the cortico-cortical conductances that can stabilize the network and that there is a critical level of fluctuations (controllable by sparsity) that allows strong cortical gain and the emergence of orientation-selective complex cells. The resultant sparse network also shows near contrast invariance in its selectivity and, in agreement with recent experiments, has extracellular tuning properties that are similar in pinwheel center and iso-orientation regions, whereas intracellular conductances show positional dependencies. Varying the strength of synaptic fluctuations by adjusting the sparsity of network connectivity, we identified a transition between the dynamics of bistability and without bistability. In a network with strong recurrent excitation, this transition is characterized by a near hysteretic behavior and a rapid rise of network firing rates as the synaptic drive or stimulus input is increased. We discuss the connection between this transition and orientation selectivity in our model of primary visual cortex.

Animals↗

PubMed2006

Fast numerical methods for simulating large-scale integrate-and-fire neuronal networks.

We discuss numerical methods for simulating large-scale, integrate-and-fire (I&F) neuronal networks. Important elements in our numerical methods are (i) a neurophysiologically inspired integrating factor which casts the solution as a numerically tractable integral equation, and allows us to obtain stable and accurate individual neuronal trajectories (i.e., voltage and conductance time-courses) even when the I&F neuronal equations are stiff, such as in strongly fluctuating, high-conductance states; (ii) an iterated process of spike-spike corrections within groups of strongly coupled neurons to account for spike-spike interactions within a single large numerical time-step; and (iii) a clustering procedure of firing events in the network to take advantage of localized architectures, such as spatial scales of strong local interactions, which are often present in large-scale computational models-for example, those of the primary visual cortex. (We note that the spike-spike corrections in our methods are more involved than the correction of single neuron spike-time via a polynomial interpolation as in the modified Runge-Kutta methods commonly used in simulations of I&F neuronal networks.) Our methods can evolve networks with relatively strong local interactions in an asymptotically optimal way such that each neuron fires approximately once in [Formula: see text] operations, where N is the number of neurons in the system. We note that quantifications used in computational modeling are often statistical, since measurements in a real experiment to characterize physiological systems are typically statistical, such as firing rate, interspike interval distributions, and spike-triggered voltage distributions. We emphasize that it takes much less computational effort to resolve statistical properties of certain I&F neuronal networks than to fully resolve trajectories of each and every neuron within the system. For networks operating in realistic dynamical regimes, such as strongly fluctuating, high-conductance states, our methods are designed to achieve statistical accuracy when very large time-steps are used. Moreover, our methods can also achieve trajectory-wise accuracy when small time-steps are used.

Algorithms↗

PubMed2006

Maximum-entropy closures for kinetic theories of neuronal network dynamics.

We analyze (1 + 1)D kinetic equations for neuronal network dynamics, which are derived via an intuitive closure from a Boltzmann-like equation governing the evolution of a one-particle (i.e., one-neuron) probability density function. We demonstrate that this intuitive closure is a generalization of moment closures based on the maximum-entropy principle. By invoking maximum-entropy closures, we show how to systematically extend this kinetic theory to obtain higher-order, kinetic equations and to include coupled networks of both excitatory and inhibitory neurons.

Algorithms↗

PubMed2005

Modeling the spatiotemporal cortical activity associated with the line-motion illusion in primary visual cortex.

Our large-scale computational model of the primary visual cortex that incorporates orientation-specific, long-range couplings with slow NMDA conductances operates in a fluctuating dynamic state of intermittent desuppression (IDS), which captures the behavior of coherent spontaneous cortical activity, as revealed by in vivo optical imaging based on voltage-sensitive dyes. Here, we address the functional significance of the IDS cortical operating points by investigating our model cortex response to the Hikosaka line-motion illusion (LMI) stimulus-a cue of a quickly flashed stationary square followed a few milliseconds later by a stationary bar. As revealed by voltage-sensitive dye imaging, there is an intriguing similarity between the cortical spatiotemporal activity in response to (i) the Hikosaka LMI stimulus and (ii) a small moving square. This similarity is believed to be associated with the preattentive illusory motion perception. Our numerical cortex produces similar spatiotemporal patterns in response to the two stimuli above, which are both in very good agreement with experimental results. The essential network mechanisms underpinning the LMI phenomenon in our model are (i) the spatiotemporal structure of the LMI input as sculpted by the lateral geniculate nucleus, (ii) a priming effect of the long-range NMDA-type cortical coupling, and (iii) the NMDA conductance-voltage correlation manifested in the IDS state. This mechanism in our model cortex, in turn, suggests a physiological underpinning for the LMI-associated patterns in the visual cortex of anaesthetized cat.

Animals↗

PubMed2005

Renormalized waves and discrete breathers in Beta-Fermi-Pasta-Ulam chains.

We demonstrate via numerical simulation that in the strongly nonlinear limit the Beta-Fermi-Pasta-Ulam (Beta-FPU) system in thermal equilibrium behaves surprisingly like weakly nonlinear waves in properly renormalized normal variables. This arises because the collective effect of strongly nonlinear interactions effectively renormalizes linear dispersion frequency and leads to effectively weak interaction among these renormalized waves. Furthermore, we show that the dynamical scenario for thermalized Beta-FPU chains is spatially highly localized discrete breathers riding chaotically on spatially extended, renormalized waves.

Journal Article↗

PubMed2005

Architectural and synaptic mechanisms underlying coherent spontaneous activity in V1.

To investigate the existence and the characteristics of possible cortical operating points of the primary visual cortex, as manifested by the coherent spontaneous ongoing activity revealed by real-time optical imaging based on voltage-sensitive dyes, we studied numerically a very large-scale ( approximately 5 x 10(5)) conductance-based, integrate-and-fire neuronal network model of an approximately 16-mm(2) patch of 64 orientation hypercolumns, which incorporates both isotropic local couplings and lateral orientation-specific long-range connections with a slow NMDA component. A dynamic scenario of an intermittent desuppressed state (IDS) is identified in the computational model, which is a dynamic state of (i) high conductance, (ii) strong inhibition, and (iii) large fluctuations that arise from intermittent spiking events that are strongly correlated in time as well as in orientation domains, with the correlation time of the fluctuations controlled by the NMDA decay time scale. Our simulation results demonstrate that the IDS state captures numerically many aspects of experimental observation related to spontaneous ongoing activity, and the specific network mechanism of the IDS may suggest cortical mechanisms and the cortical operating point underlying observed spontaneous activity.

Animals↗

PubMed2004

An embedded network approach for scale-up of fluctuation-driven systems with preservation of spike information.

To address computational "scale-up" issues in modeling large regions of the cortex, many coarse-graining procedures have been invoked to obtain effective descriptions of neuronal network dynamics. However, because of local averaging in space and time, these methods do not contain detailed spike information and, thus, cannot be used to investigate, e.g., cortical mechanisms that are encoded through detailed spike-timing statistics. To retain high-order statistical information of spikes, we develop a hybrid theoretical framework that embeds a subnetwork of point neurons within, and fully interacting with, a coarse-grained network of dynamical background. We use a newly developed kinetic theory for the description of the coarse-grained background, in combination with a Poisson spike reconstruction procedure to ensure that our method applies to the fluctuation-driven regime as well as to the mean-driven regime. This embedded-network approach is verified to be dynamically accurate and numerically efficient. As an example, we use this embedded representation to construct "reverse-time correlations" as spiked-triggered averages in a ring model of orientation-tuning dynamics.

Cerebral Cortex↗

PubMed2004

An effective kinetic representation of fluctuation-driven neuronal networks with application to simple and complex cells in visual cortex.

A coarse-grained representation of neuronal network dynamics is developed in terms of kinetic equations, which are derived by a moment closure, directly from the original large-scale integrate-and-fire (I&F) network. This powerful kinetic theory captures the full dynamic range of neuronal networks, from the mean-driven limit (a limit such as the number of neurons N --> infinity, in which the fluctuations vanish) to the fluctuation-dominated limit (such as in small N networks). Comparison with full numerical simulations of the original I&F network establishes that the reduced dynamics is very accurate and numerically efficient over all dynamic ranges. Both analytical insights and scale-up of numerical representation can be achieved by this kinetic approach. Here, the theory is illustrated by a study of the dynamical properties of networks of various architectures, including excitatory and inhibitory neurons of both simple and complex type, which exhibit rich dynamic phenomena, such as, transitions to bistability and hysteresis, even in the presence of large fluctuations. The implication for possible connections between the structure of the bifurcations and the behavior of complex cells is discussed. Finally, I&F networks and kinetic theory are used to discuss orientation selectivity of complex cells for "ring-model" architectures that characterize changes in the response of neurons located from near "orientation pinwheel centers" to far from them.

Computer Simulation↗