Simulating human sleep spindle MEG and EEG from ion channel and circuit level dynamics.

BACKGROUND: Although they form a unitary phenomenon, the relationship between extracranial M/EEG and transmembrane ion flows is understood only as a general principle rather than as a well-articulated and quantified causal chain. METHOD: We present an integrated multiscale model, consisting of a neural simulation of thalamus and cortex during stage N2 sleep and a biophysical model projecting cortical current densities to M/EEG fields. Sleep spindles were generated through the interactions of local and distant network connections and intrinsic currents within thalamocortical circuits. 32,652 cortical neurons were mapped onto the cortical surface reconstructed from subjects' MRI, interconnected based on geodesic distances, and scaled-up to current dipole densities based on laminar recordings in humans. MRIs were used to generate a quasi-static electromagnetic model enabling simulated cortical activity to be projected to the M/EEG sensors. RESULTS: The simulated M/EEG spindles were similar in amplitude and topography to empirical examples in the same subjects. Simulated spindles with more core-dominant activity were more MEG weighted. COMPARISON WITH EXISTING METHODS: Previous models lacked either spindle-generating thalamic neural dynamics or whole head biophysical modeling; the framework presented here is the first to simultaneously capture these disparate scales. CONCLUSIONS: This multiscale model provides a platform for the principled quantitative integration of existing information relevant to the generation of sleep spindles, and allows the implications of future findings to be explored. It provides a proof of principle for a methodological framework allowing large-scale integrative brain oscillations to be understood in terms of their underlying channels and synapses.

Peptide-oligonucleotide conjugates exhibiting pyrimidine-X cleavage specificity efficiently silence miRNA target acting synergistically with RNase H.

Taking into account the important role of miRNA in carcinogenesis, oncogenic miRNAs are attractive molecules for gene-targeted therapy. Here, we developed a novel series of peptide-oligonucleotide conjugates exhibiting ribonuclease activity targeted to highly oncogenic miRNAs miR-21 and miR-17. When designing the conjugates, we enhanced both nuclease resistance of the targeted oligodeoxyribonucleotide by introducing at its 3'-end mini-hairpin structure displaying high thermostability and robustness against nuclease digestion and the efficiency of its functioning by attachment of the catalytic construction (amide)NH2-Gly(ArgLeu)4-TCAA displaying ribonuclease activity to its 5'-end. Designed miRNases efficiently cleaved miRNA targets, exhibiting Pyr-X specificity, and cleavage specificity had strong dependence on the miRNA sequence in the site of peptide location. In vitro, designed miRNases do not prevent cleavage of miRNA bound with the conjugate by RNase H, and more than an 11-fold enhancement of miRNA cleavage by the conjugate is observed in the presence of RNase H. In murine melanoma cells, miRNase silences mmu-miR-17 with very high efficiency as a result of miR-17 cleavage by miRNase and by recruited RNase H. Thus, miRNases provide a system of double attack of the miRNA molecules, significantly increasing the efficiency of miRNA downregulation in the cells in comparison with antisense oligonucleotide.

Pre-harvest sprouting resistance and haplotype variation of ThVp-1 gene in the collection of wheat-wheatgrass hybrids.

The germplasm collection of 87 wheat-wheatgrass hybrids developed in Tsitisin Main Botanical Garden (Russia, Moscow) was evaluated for resistance to pre-harvest sprouting (PHS) by spike sprouting (SS) and germination index (GI) assays as well as for spike and grain features. The PHS resistance variation and haplotype polymorphism of the wheatgrass ThVp-1 and wheat TaVp-1B genes orthologues of Vp-1 was revealed in the studied collection. Four haplotypes of ThVp-1 were revealed: ThVp-1a (41% of the entries), ThVp-1b (13%), ThVp-1c (29%), and ThVp-1d (15%). The association between the allelic state of ThVp-1 and PHS resistance in the wheat-wheatgrass hybrids was shown: haplotype ThVp-1d of the wheatgrass Vp-1 gene is significantly associated with reduced PHS in the wheat-wheatgrass hybrids (mean SS 0.33, mean GI 0.64). The resistant entries may be perspective as a source of PHS resistance in the development of commercial cultivars of perennial wheat.

[Distal reconstructions under critical lower limb ischemia in elderly patients].

Performong reconstructive operations on arteriae is the only way to avoid limb amputation in patients with critical limb ischemia. Whereas reconstructive manipulations on distal arteriae, especially in patients older then 75 years, are often entailed by thrombotic complications. High peripheric vessel resistance is considered to be the crucial factor of the bypass thrombosis. The conducted study demonstrated the efficacy of the use of vena saphena magna collateral tributaries as bypass unload collectors. Femoral distal autovenous in situ reconstruction decreases vessel resistance. Due to the arterial reconstruction with bypass unload using unligated autovenous tributaries the cumulative bypass patency was 75,1% a year after the procedure. 5-year cumulative bypass patency was 53,8%. 5 - year lethality was 28,6%.

Effect of synaptic connectivity on long-range synchronization of fast cortical oscillations.

Cortical gamma oscillations in the 20- to 80-Hz range are associated with attentiveness and sensory perception and have strong connections to both cognitive processing and temporal binding of sensory stimuli. These gamma oscillations become synchronized within a few milliseconds over distances spanning a few millimeters in spite of synaptic delays. In this study using in vivo recordings and large-scale cortical network models, we reveal a critical role played by the network geometry in achieving precise long-range synchronization in the gamma frequency band. Our results indicate that the presence of many independent synaptic pathways in a two-dimensional network facilitate precise phase synchronization of fast gamma band oscillations with nearly zero phase delays between remote network sites. These findings predict a common mechanism of precise oscillatory synchronization in neuronal networks.

Oscillations in large-scale cortical networks: map-based model.

We develop a new computationally efficient approach for the analysis of complex large-scale neurobiological networks. Its key element is the use of a new phenomenological model of a neuron capable of replicating important spike pattern characteristics and designed in the form of a system of difference equations (a map). We developed a set of map-based models that replicate spiking activity of cortical fast spiking, regular spiking and intrinsically bursting neurons. Interconnected with synaptic currents these model neurons demonstrated responses very similar to those found with Hodgkin-Huxley models and in experiments. We illustrate the efficacy of this approach in simulations of one- and two-dimensional cortical network models consisting of regular spiking neurons and fast spiking interneurons to model sleep and activated states of the thalamocortical system. Our study suggests that map-based models can be widely used for large-scale simulations and that such models are especially useful for tasks where the modeling of specific firing patterns of different cell classes is important.

Potassium model for slow (2-3 Hz) in vivo neocortical paroxysmal oscillations.

In slow neocortical paroxysmal oscillations, the de- and hyperpolarizing envelopes in neocortical neurons are large compared with slow sleep oscillations. Increased local synchrony of membrane potential oscillations during seizure is reflected in larger electroencephalographic oscillations and the appearance of spike- or polyspike-wave complex recruitment at 2- to 3-Hz frequencies. The oscillatory mechanisms underlying this paroxysmal activity were investigated in computational models of cortical networks. The extracellular K(+) concentration ([K(+)](o)) was continuously computed based on neuronal K(+) currents and K(+) pumps as well as glial buffering. An increase of [K(+)](o) triggered a transition from normal awake-like oscillations to 2- to 3-Hz seizure-like activity. In this mode, the cells fired periodic bursts and nearby neurons oscillated highly synchronously; in some cells depolarization led to spike inactivation lasting 50-100 ms. A [K(+)](o) increase, sufficient to produce oscillations could result from excessive firing (e.g., induced by external stimulation) or inability of K(+) regulatory system (e.g., when glial buffering was blocked). A combination of currents including high-threshold Ca(2+), persistent Na(+) and hyperpolarization-activated depolarizing (I(h)) currents was sufficient to maintain 2- to 3-Hz activity. In a network model that included lateral K(+) diffusion between cells, increase of [K(+)](o) in a small region was generally sufficient to maintain paroxysmal oscillations in the whole network. Slow changes of [K(+)](o) modulated the frequency of bursting and, in some case, led to fast oscillations in the 10- to 15-Hz frequency range, similar to the fast runs observed during seizures in vivo. These results suggest that modifications of the intrinsic currents mediated by increase of [K(+)](o) can explain the range of neocortical paroxysmal oscillations in vivo.

Model of transient oscillatory synchronization in the locust antennal lobe.

Transient pairwise synchronization of locust antennal lobe (AL) projection neurons (PNs) occurs during odor responses. In a Hodgkin-Huxley-type model of the AL, interactions between excitatory PNs and inhibitory local neurons (LNs) created coherent network oscillations during odor stimulation. GABAergic interconnections between LNs led to competition among them such that different groups of LNs oscillated with periodic Ca(2+) spikes during different 50-250 ms temporal epochs, similar to those recorded in vivo. During these epochs, LN-evoked IPSPs caused phase-locked, population oscillations in sets of postsynaptic PNs. The model shows how alternations of the inhibitory drive can temporally encode sensory information in networks of neurons without precisely tuned intrinsic oscillatory properties.

Model of cellular and network mechanisms for odor-evoked temporal patterning in the locust antennal lobe.

Locust antennal lobe (AL) projection neurons (PNs) respond to olfactory stimuli with sequences of depolarizing and hyperpolarizing epochs, each lasting hundreds of milliseconds. A computer simulation of an AL network was used to test the hypothesis that slow inhibitory connections between local neurons (LNs) and PNs are responsible for temporal patterning. Activation of slow inhibitory receptors on PNs by the same GABAergic synapses that underlie fast oscillatory synchronization of PNs was sufficient to shape slow response modulations. This slow stimulus- and neuron-specific patterning of AL activity was resistant to blockade of fast inhibition. Fast and slow inhibitory mechanisms at synapses between LNs and PNs can thus form dynamical PN assemblies whose elements synchronize transiently and oscillate collectively, as observed not only in the locust AL, but also in the vertebrate olfactory bulb.

Origin of slow cortical oscillations in deafferented cortical slabs.

An in vivo preparation has been developed to study the mechanisms underlying spontaneous sleep oscillations. Dual and triple simultaneous intracellular recordings were made from neurons in small isolated cortical slabs (10 mm x 6 mm) in anesthetized cats. Spontaneously occurring slow sleep oscillations, present in the adjacent intact cortex, were absent in small slabs. However, the isolated slabs displayed brief active periods separated by long periods of silence, up to 60 s in duration. During these silent periods, 60% of neurons showed non-linear amplification of low-amplitude depolarizing activity. Nearly 40% of the cells, twice as many as in intact cortex, were classified as intrinsically bursting. In cortical network models based on Hodgkin-Huxley-like neurons, the summation of simulated spontaneous miniature excitatory postsynaptic potentials was sufficient to activate a persistent sodium current, initiating action potentials in single neurons that then spread through the network. Consistent with this model, enlarging the isolated cortical territory to an isolated gyrus (30 mm x 20 mm) increased the probability of initiating large-scale activity. In these larger territories, both the frequency and regularity of the slow oscillation approached that generated in intact cortex. The frequency of active periods in an analytical model of the cortical network accurately predicted the scaling observed in simulations and from recordings in cortical slabs of increasing size.

Spiking-bursting activity in the thalamic reticular nucleus initiates sequences of spindle oscillations in thalamic networks.

Recent intracellular and local field potential recordings from thalamic reticular (RE) neurons in vivo as well as computational modeling of the isolated RE nucleus suggest that, at relatively hyperpolarized levels of membrane potentials, the inhibitory postsynaptic potentials (IPSPs) between RE cells can be reversed and gamma-aminobutyric acid-A (GABA(A)) -mediated depolarization can generate persistent spatio-temporal patterns in the RE nucleus. Here we investigate how this activity affects the spatio-temporal properties of spindle oscillations with computer models of interacting RE and thalamocortical (TC) cells. In a one-dimensional network of RE and TC cells, sequences of spindle oscillations alternated with localized patterns of spike-burst activity propagating inside the RE network. New sequences of spindle oscillations were initiated after removal of I(h)-mediated depolarization of the TC cells. The length of the interspindle lulls depended on the intrinsic and synaptic properties of RE and TC cells and was in the range of 3-20 s. In a two-dimensional model, GABA(A)-mediated 2-3 Hz oscillations persisted in the RE nucleus during interspindle lulls and initiated spindle sequences at many foci within the RE-TC network simultaneously. This model predicts that the intrinsic properties of the reticular thalamus may contribute to the synchrony of spindle oscillations observed in vivo.

Self-sustained rhythmic activity in the thalamic reticular nucleus mediated by depolarizing GABAA receptor potentials.

Intracellular recordings from reticular thalamic (RE) neurons in vivo revealed inhibitory postsynaptic potentials (IPSPs) between RE cells that reversed and became depolarizing at the hyperpolarized membrane potentials that occur during sleep. These excitatory IPSPs can directly trigger low-threshold spikes (LTSs). The oscillatory mechanisms underlying IPSP-triggered LTSs crowned by spike bursts were investigated in models of isolated RE networks. In a one-dimensional network model, external stimulation evoked waves of excitation propagating at a constant velocity of 25-150 cells per second. In a large-scale, two-dimensional model of the reticular nucleus, the network showed transient or self-sustained oscillations controlled by the maximum conductance of the low-threshold calcium current and the membrane potential. This model predicts that the isolated reticular nucleus could initiate sequences of spindle oscillations in thalamocortical networks in vivo.

Computational models of thalamocortical augmenting responses.

Repetitive stimulation of the dorsal thalamus at 7-14 Hz produces an increasing number of spikes at an increasing frequency in neocortical neurons during the first few stimuli. Possible mechanisms underlying these cortical augmenting responses were analyzed with a computer model that included populations of thalamocortical cells, thalamic reticular neurons, up to two layers of cortical pyramidal cells, and cortical inhibitory interneurons. Repetitive thalamic stimulation produced a low-threshold intrathalamic augmentation in the model based on the deinactivation of the low-threshold Ca2+ current in thalamocortical cells, which in turn induced cortical augmenting responses. In the cortical model, augmenting responses were more powerful in the "input" layer compared with those in the "output" layer. Cortical stimulation of the network model produced augmenting responses in cortical neurons in distant cortical areas through corticothalamocortical loops and low-threshold intrathalamic augmentation. Thalamic stimulation was more effective in eliciting augmenting responses than cortical stimulation. Intracortical inhibition had an important influence on the genesis of augmenting responses in cortical neurons: A shift in the balance between intracortical excitation and inhibition toward excitation transformed an augmenting responses to long-lasting paroxysmal discharge. The predictions of the model were compared with in vivo recordings from neurons in cortical area 4 and thalamic ventrolateral nucleus of anesthetized cats. The known intrinsic properties of thalamic cells and thalamocortical interconnections can account for the basic properties of cortical augmenting responses.

Cellular and network models for intrathalamic augmenting responses during 10-Hz stimulation.

Repetitive stimulation of the thalamus at 7-14 Hz evokes responses of increasing amplitude in the thalamus and the areas of the neocortex to which the stimulated foci project. Possible mechanisms underlying the thalamic augmenting responses during repetitive stimulation were investigated with computer models of interacting thalamocortical (TC) and thalamic reticular (RE) cells. The ionic currents in these cells were modeled with Hodgkin-Huxley type of kinetics, and the results of the model were compared with in vivo thalamic recordings from decorticated cats. The simplest network model demonstrating an augmenting response was a single pair of coupled RE and TC cells, in which RE-induced inhibitory postsynaptic potentials (IPSPs) in the TC cell led to progressive deinactivation of a low-threshold Ca2+ current. The augmenting responses in two reciprocally interacting chains of RE and TC cells depended also on gamma-aminobutyric acid-B (GABAB) IPSPs. Lateral GABAA inhibition between identical RE cells, which weakened bursts in these cells, diminished GABAB IPSPs and delayed the augmenting response in TC cells. The results of these simulations show that the interplay between existing mechanisms in the thalamus explains the basic properties of the intrathalamic augmenting responses.