Transforming neural computations and representing time

Motifs of neural circuitry seem surprisingly conserved over different areas of neocortex or of paleocortex, while performing quite different sensory processing tasks. This apparent paradox may be resolved by the fact that seemingly different problems in sensory information processing are related by transformations (changes of variables) that convert one problem into another. The same basic algorithm that is appropriate to the recognition of a known odor quality, independent of the strength of the odor, can be used to recognize a vocalization (e.g., a spoken syllable), independent of whether it is spoken quickly or slowly. To convert one problem into the other, a new representation of time sequences is needed. The time that has elapsed since a recent event must be represented in neural activity. The electrophysiological hallmarks of cells that are involved in generating such a representation of time are discussed. The anatomical relationships between olfactory and auditory pathways suggest relevant experiments. The neurophysiological mechanism for the psychophysical logarithmic encoding of time duration would be of direct use for interconverting olfactory and auditory processing problems. Such reuse of old algorithms in new settings and representations is related to the way that evolution develops new biochemistry.

A logical analysis of T cell activation and anergy

Interaction of the antigen-specific receptor of T lymphocytes with its antigenic ligand can lead either to cell activation or to a state of profound unresponsiveness (anergy). Although subtle changes in the nature of the ligand or of the antigen-presenting cell have been shown to affect the outcome of T cell receptor ligation, the mechanism by which the same receptor can induce alternative cellular responses is not completely understood. A model for explaining both positive (cell proliferation and cytokine production) and negative (anergy induction) signaling of T lymphocytes is described herein. This model relies on the autophosphorylative properties of the tyrosine kinases associated with the T cell receptor. One of its basic assumptions is that the kinase activity of these receptor-associated enzymes remains above background level after ligand removal and is responsible for cellular unresponsiveness. Using a simple Boolean formalism, we show how the timing of the binding and intracellular signal-transduction events can affect the properties of receptor signaling and determine the type of cellular response. The present approach integrates into a common framework a large body of experimental observations and allows specification of conditions leading to cellular activation or to anergy.

Conformation of coenzyme pyrroloquinoline quinone and role of Ca2+ in the catalytic mechanism of quinoprotein methanol dehydrogenase

The ab initio structures of 2,7,9-tricarboxypyrroloquinoline quinone (PQQ), semiquinone (PQQH), and dihydroquinone (PQQH2) have been determined and compared with ab initio structures of the (PQQ)Ca2+, (PQQH)Ca2+, and (PQQH2)Ca2+ complexes as well as the x-ray structure of (PQQ)Ca2+ bound at the active site of the methanol dehydrogenase (MDH) of methyltropic bacteria. Plausible mechanisms for the MDH oxidation of methanol involving the (PQQ)Ca2+ complex are explored via ab initio computations and discussed. Considering the reaction of methanol with PQQ in the absence of Ca2+, nucleophilic addition of methanol to the PQQ C-5 carbonyl followed by a retro-ene elimination is deemed unlikely due to large energy barrier. A much more favorable disposition of the methanol C-5 adduct to provide formaldehyde involves proton ionization of the intermediate followed by elimination of methoxide concerted with hydride transfer to the oxygen of the C-4 carbonyl. Much the same transition state is reached if one searches for the transition state beginning with Asp-303–CO2−general-base removal of the methanol proton of the (PQQ)Ca2+O(H)CH3 complex concerted with hydride transfer to the oxygen at C-4. For such a mechanism the role of the Ca2+ moiety would be to (i) contribute to the formation of the ES complex (ii) provide a modest decrease in the pKa of methanol substrate,; and (iii) polarize the oxygen at C-5.

An empirical explanation of brightness

In this second part of our study on the mechanism of perceived brightness, we explore the effects of manipulating three-dimensional geometry. The additional scenes portrayed here demonstrate that the same luminance profile can elicit different sensations of brightness as a function of how the objects in the scene are arranged in space. This further evidence confirms the implication of the scenes presented in the accompanying paper, namely that sensations of relative brightness—including standard demonstrations of simultaneous brightness contrast—cannot arise by computations of local contrast. The most plausible explanation of the full range of perceptual phenomena we have described is an empirical strategy that links the luminance profile in a visual stimulus with an association (the percept) that represents the profile’s most probable real-world source.

Visual stimuli induce waves of electrical activity in turtle cortex

The computations involved in the processing of a visual scene invariably involve the interactions among neurons throughout all of visual cortex. One hypothesis is that the timing of neuronal activity, as well as the amplitude of activity, provides a means to encode features of objects. The experimental data from studies on cat [Gray, C. M., Konig, P., Engel, A. K. & Singer, W. (1989) Nature (London) 338, 334–337] support a view in which only synchronous (no phase lags) activity carries information about the visual scene. In contrast, theoretical studies suggest, on the one hand, the utility of multiple phases within a population of neurons as a means to encode independent visual features and, on the other hand, the likely existence of timing differences solely on the basis of network dynamics. Here we use widefield imaging in conjunction with voltage-sensitive dyes to record electrical activity from the virtually intact, unanesthetized turtle brain. Our data consist of single-trial measurements. We analyze our data in the frequency domain to isolate coherent events that lie in different frequency bands. Low frequency oscillations (<5 Hz) are seen in both ongoing activity and activity induced by visual stimuli. These oscillations propagate parallel to the afferent input. Higher frequency activity, with spectral peaks near 10 and 20 Hz, is seen solely in response to stimulation. This activity consists of plane waves and spiral-like waves, as well as more complex patterns. The plane waves have an average phase gradient of ≈π/2 radians/mm and propagate orthogonally to the low frequency waves. Our results show that large-scale differences in neuronal timing are present and persistent during visual processing.

Binocular vision in insects: How mantids solve the correspondence problem

Praying mantids use binocular cues to judge whether their prey is in striking distance. When there are several moving targets within their binocular visual field, mantids need to solve the correspondence problem. They must select between the possible pairings of retinal images in the two eyes so that they can strike at a single real target. In this study, mantids were presented with two targets in various configurations, and the resulting fixating saccades that precede the strike were analyzed. The distributions of saccades show that mantids consistently prefer one out of several possible matches. Selection is in part guided by the position and the spatiotemporal features of the target image in each eye. Selection also depends upon the binocular disparity of the images, suggesting that insects can perform local binocular computations. The pairing rules ensure that mantids tend to aim at real targets and not at “ghost” targets arising from false matches.

Computation, prediction, and experimental tests of fitness for bacteriophage T7 mutants with permuted genomes

We created a simulation based on experimental data from bacteriophage T7 that computes the developmental cycle of the wild-type phage and also of mutants that have an altered genome order. We used the simulation to compute the fitness of more than 105 mutants. We tested these computations by constructing and experimentally characterizing T7 mutants in which we repositioned gene 1, coding for T7 RNA polymerase. Computed protein synthesis rates for ectopic gene 1 strains were in moderate agreement with observed rates. Computed phage-doubling rates were close to observations for two of four strains, but significantly overestimated those of the other two. Computations indicate that the genome organization of wild-type T7 is nearly optimal for growth: only 2.8% of random genome permutations were computed to grow faster, the highest 31% faster, than wild type. Specific discrepancies between computations and observations suggest that a better understanding of the translation efficiency of individual mRNAs and the functions of qualitatively “nonessential” genes will be needed to improve the T7 simulation. In silico representations of biological systems can serve to assess and advance our understanding of the underlying biology. Iteration between computation, prediction, and observation should increase the rate at which biological hypotheses are formulated and tested.

Gamma rhythms and beta rhythms have different synchronization properties

Experimental and modeling efforts suggest that rhythms in the CA1 region of the hippocampus that are in the beta range (12–29 Hz) have a different dynamical structure than that of gamma (30–70 Hz). We use a simplified model to show that the different rhythms employ different dynamical mechanisms to synchronize, based on different ionic currents. The beta frequency is able to synchronize over long conduction delays (corresponding to signals traveling a significant distance in the brain) that apparently cannot be tolerated by gamma rhythms. The synchronization properties are consistent with data suggesting that gamma rhythms are used for relatively local computations whereas beta rhythms are used for higher level interactions involving more distant structures.

Anatomy of word and sentence meaning

Reading and listening involve complex psychological processes that recruit many brain areas. The anatomy of processing English words has been studied by a variety of imaging methods. Although there is widespread agreement on the general anatomical areas involved in comprehending words, there are still disputes about the computations that go on in these areas. Examination of the time relations (circuitry) among these anatomical areas can aid in understanding their computations. In this paper, we concentrate on tasks that involve obtaining the meaning of a word in isolation or in relation to a sentence. Our current data support a finding in the literature that frontal semantic areas are active well before posterior areas. We use the subject’s attention to amplify relevant brain areas involved either in semantic classification or in judging the relation of the word to a sentence to test the hypothesis that frontal areas are concerned with lexical semantics and posterior areas are more involved in comprehension of propositions that involve several words.

Image recognition: Visual grouping, recognition, and learning

Vision extracts useful information from images. Reconstructing the three-dimensional structure of our environment and recognizing the objects that populate it are among the most important functions of our visual system. Computer vision researchers study the computational principles of vision and aim at designing algorithms that reproduce these functions. Vision is difficult: the same scene may give rise to very different images depending on illumination and viewpoint. Typically, an astronomical number of hypotheses exist that in principle have to be analyzed to infer a correct scene description. Moreover, image information might be extracted at different levels of spatial and logical resolution dependent on the image processing task. Knowledge of the world allows the visual system to limit the amount of ambiguity and to greatly simplify visual computations. We discuss how simple properties of the world are captured by the Gestalt rules of grouping, how the visual system may learn and organize models of objects for recognition, and how one may control the complexity of the description that the visual system computes.

Neural codes: Firing rates and beyond

Computational neuroscience has contributed significantly to our understanding of higher brain function by combining experimental neurobiology, psychophysics, modeling, and mathematical analysis. This article reviews recent advances in a key area: neural coding and information processing. It is shown that synapses are capable of supporting computations based on highly structured temporal codes. Such codes could provide a substrate for unambiguous representations of complex stimuli and be used to solve difficult cognitive tasks, such as the binding problem. Unsupervised learning rules could generate the circuitry required for precise temporal codes. Together, these results indicate that neural systems perform a rich repertoire of computations based on action potential timing.

The generation of oscillations in networks of electrically coupled cells

In several biological systems, the electrical coupling of nonoscillating cells generates synchronized membrane potential oscillations. Because the isolated cell is nonoscillating and electrical coupling tends to equalize the membrane potentials of the coupled cells, the mechanism underlying these oscillations is unclear. Here we present a dynamic mechanism by which the electrical coupling of identical nonoscillating cells can generate synchronous membrane potential oscillations. We demonstrate this mechanism by constructing a biologically feasible model of electrically coupled cells, characterized by an excitable membrane and calcium dynamics. We show that strong electrical coupling in this network generates multiple oscillatory states with different spatio-temporal patterns and discuss their possible role in the cooperative computations performed by the system.

Temporal fluctuations in coherence of brain waves.

As a measure of dynamical structure, short-term fluctuations of coherence between 0.3 and 100 Hz in the electroencephalogram (EEG) of humans were studied from recordings made by chronic subdural macroelectrodes 5-10 mm apart, on temporal, frontal, and parietal lobes, and from intracranial probes deep in the temporal lobe, including the hippocampus, during sleep, alert, and seizure states. The time series of coherence between adjacent sites calculated every second or less often varies widely in stability over time; sometimes it is stable for half a minute or more. Within 2-min samples, coherence commonly fluctuates by a factor up to 2-3, in all bands, within the time scale of seconds to tens of seconds. The power spectrum of the time series of these fluctuations is broad, extending to 0.02 Hz or slower, and is weighted toward the slower frequencies; little power is faster than 0.5 Hz. Some records show conspicuous swings with a preferred duration of 5-15s, either irregularly or quasirhythmically with a broad peak around 0.1 Hz. Periodicity is not statistically significant in most records. In our sampling, we have not found a consistent difference between lobes of the brain, subdural and depth electrodes, or sleeping and waking states. Seizures generally raise the mean coherence in all frequencies and may reduce the fluctuations by a ceiling effect. The coherence time series of different bands is positively correlated (0.45 overall); significant nonindependence extends for at least two octaves. Coherence fluctuations are quite local; the time series of adjacent electrodes is correlated with that of the nearest neighbor pairs (10 mm) to a coefficient averaging approximately 0.4, falling to approximately 0.2 for neighbors-but-one (20 mm) and to < 0.1 for neighbors-but-two (30 mm). The evidence indicates fine structure in time and space, a dynamic and local determination of this measure of cooperativity. Widely separated frequencies tending to fluctuate together exclude independent oscillators as the general or usual basis of the EEG, although a few rhythms are well known under special conditions. Broad-band events may be the more usual generators. Loci only a few millimeters apart can fluctuate widely in seconds, either in parallel or independently. Scalp EEG coherence cannot be predicted from subdural or deep recordings, or vice versa, and intracortical microelectrodes show still greater coherence fluctuation in space and time. Widely used computations of chaos and dimensionality made upon data from scalp or even subdural or depth electrodes, even when reproducible in successive samples, cannot be considered representative of the brain or the given structure or brain state but only of the scale or view (receptive field) of the electrodes used. Relevant to the evolution of more complex brains, which is an outstanding fact of animal evolution, we believe that measures of cooperativity are likely to be among the dynamic features by which major evolutionary grades of brains differ.

Calicheamicin-DNA complexes: warhead alignment and saccharide recognition of the minor groove.

The solution structures of calicheamicin gamma 1I, its cycloaromatized analog (calicheamicin epsilon), and its aryl tetrasaccharide complexed to a common DNA hairpin duplex have been determined by NMR and distance-refined molecular dynamics computations. Sequence specificity is associated with carbohydrate-DNA recognition that places the aryl tetrasaccharide component of all three ligands in similar orientations in the minor groove at the d(T-C-C-T).d(A-G-G-A) segment. The complementary fit of the ligands and the DNA minor groove binding site creates numerous van der Waals contacts as well as hydrogen bonding interactions. Notable are the iodine and sulfur atoms of calicheamicin that hydrogen bond with the exposed amino proton of the 5'- and 3'-guanines, respectively, of the d(A-G-G-A) segment. The sequence-specific carbohydrate binding orients the enediyne aglycone of calicheamicin gamma 1I such that its C3 and C6 proradical centers are adjacent to the cleavage sites. While the enediyne aglycone of calicheamicin gamma 1I is tilted relative to the helix axis and spans the minor groove, the cycloaromatized aglycone is aligned approximately parallel to the helix axis in the respective complexes. Specific localized conformational perturbations in the DNA have been identified from imino proton complexation shifts and changes in specific sugar pucker patterns on complex formation. The helical parameters for the carbohydrate binding site are comparable with corresponding values in B-DNA fibers while a widening of the groove is observed at the adjacent aglycone binding site.