Algorithmic statistical analysis of electrophysiological data for the investigation of structure-activity relationship in single neurons

We are developing computational tools supporting the detailed analysis of the dependence of neural electrophysiological response on dendritic morphology. We approach this problem by combining simulations of faithful models of neurons (experimental real life morphological data with known models of channel kinetics) with algorithmic extraction of morphological and physiological parameters and statistical analysis. In this paper, we present the novel method for an automatic recognition of spike trains in voltage traces, which eliminates the need for human intervention. This enables classification of waveforms with consistent criteria across all the analyzed traces and so it amounts to reduction of the noise in the data. This method allows for an automatic extraction of relevant physiological parameters necessary for further statistical analysis. In order to illustrate the usefulness of this procedure to analyze voltage traces, we characterized the influence of the somatic current injection level on several electrophysiological parameters in a set of modeled neurons. This application suggests that such an algorithmic processing of physiological data extracts parameters in a suitable form for further investigation of structure-activity relationship in single neurons.

Projeto multicritério de filtros harmônicos passivos para instalações industriais utilizando técnicas de inteligência computacional

Devido ao auge do crescimento industrial na Região Norte e, em especial, o Pólo Industrial de Manaus (PIM), são necessários obter ferramentas matemáticas que facilitem ao especialista tomar decisões sobre a seleção e dimensionamento dos filtros harmônicos que proporcionam neutralizar os efeitos prejudiciais dos harmônicos gerados pelas cargas não lineares da indústria e alcançar conformidade com os padrões das normas de qualidade de energia correspondentes. Além disso, como os filtros harmônicos passivos têm a capacidade de gerar potência reativa à rede, estes meios são eficazes compensadores de potência reativa e, portanto, podem conseguir uma economia significativa no faturamento de energia elétrica consumida por essas instalações industriais. Esta tese tem como objetivo geral desenvolver um método matemático e uma ferramenta computacional para a seleção da configuração e parâmetros do projeto de um conjunto de filtros harmônicos passivos para sistemas elétricos industriais. Nesta ótica, o problema de otimização da compensação de harmônicos por meio de filtros passivos foi formulado como um problema multiobjetivo que considera tanto os objetivos da redução da distorção harmônica como da efetividade econômica do projeto considerando as características das tarifas brasileiras. Todavia, a formulação apresentada considera as restrições relevantes impostas pelas normas brasileiras e estrangeiras. A solução computacional para este problema foi conseguida, usando o algoritmo genético NSGA-II que determina um conjunto de soluções ótimas de Pareto (Fronteira) que permitem ao projetista escolher as soluções mais adequadas para o problema. Por conseguinte, a ferramenta computacional desenvolvida tem várias novidades como: não só calcula os parâmetros que caracterizam os filtros, como também seleciona o tipo de configuração e o número de ramos do filtro em cada barra candidata de acordo com um conjunto de configurações pré-estabelecidas; têm implementada duas normas para a avaliação das restrições de qualidade de energia (Prodist-Módulo 8 e IEEE 519-92) que podem ser selecionadas pelo usuário; determina soluções com bons indicadores de desempenho para vários cenários característicos e não característicos do sistema que permitem a representação das as variações diárias da carga; das variações dos parâmetros do sistema e dos filtros; avalia o custo das contas de energia numa rede elétrica industrial que tem diferentes condições de operação (cenários característicos); e avalia o efeito econômico de filtros de harmônicos como compensadores de potência reativa. Para desenvolver a ferramenta computacional adequada desta tese, foi empregado um modelo trifásico em coordenadas de fase para redes de energia elétrica industriais e de serviços onde foram feitos vários programas utilizando várias ferramentas computacionais adicionais. Estas ferramentas compreendem um programa de varredura de freqüência, um programa do fluxo de harmônicos por injeção de correntes e um programa de fluxo de potência à freqüência fundamental. Os resultados positivos desta tese, a partir da análise de vários exemplos práticos, mostram as vantagens do método desenvolvido.

Fluxo de potência continuado baseado em equações de injeção de corrente expressas em coordenadas polares

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Resonance properties of different neuronal populations in the immature mouse neocortex

Subthreshold resonance is a characteristic membrane property of different neuronal classes, is critically involved in the generation of network oscillations, and tunes the integration of synaptic inputs to particular frequency ranges. In order to investigate whether resonance properties of distinct neuronal populations in the immature neocortex contribute to these network oscillations, I performed whole-cell patch-clamp recordings from visually identified neurons in tangential and coronal neocortical slices from postnatal day (P) P0-P7 C57Bl/6 and P6-P13 GAD67-GFP knock-in mice. Subthreshold resonance was analyzed by sinusoidal current injection of varying frequency. All Cajal-Retzius cells showed subthreshold resonance with an average frequency of 2.6 ± 0.1 Hz (n=60), which was massively reduced by ZD7288, a blocker of hyperpolarization-activated cation currents. About 65.6% (n=61) of the supragranular pyramidal neurons showed subthreshold resonance with an average frequency of 1.4 ± 0.1 Hz (n=40). Application of 1 mM Ni2+ suppressed subthreshold resonance, suggesting that low-threshold Ca2+ currents contribute to resonance in these neurons. About 63.6% (n=77) of the layer V pyramidal neurons showed subthreshold resonance with an average frequency of 1.4 ± 0.2 Hz (n=49), which was abolished by ZD7288. Only 44.1% (n=59) of the subplate neurons showed subthreshold resonance with an average frequency of 1.3 ± 0.2 Hz (n=26) and a small resonance strength. Finally, 50% of the investigated GABAergic interneurons showed subthreshold resonance with an average frequency of 2.0 ± 0.2 Hz (n=42). Membrane hyperpolarization to –86 mV attenuated the frequency and strength of subthreshold resonance. Subthreshold resonance was virtually abolished in the presence of 1 mM Ni2+, suggesting that t-type Ca2+ currents are critically involved in the generation of resonance, while ZD7288 had no effect. Application of 0.4 µM TTX suppressed subthreshold resonance at depolarized, but not hyperpolarized membrane potential, suggesting that persistent Na+ current contribute to the amplification of membrane resonance. rnIn summary, these results demonstrate that all investigated neuronal subpopulations reveal resonance behavior, with either hyperpolarization-activated cation or low-threshold Ca2+ currents contributing to the subthreshold resonance. GABAergic interneurons also express subthreshold resonance at low frequencies, with t-type Ca2+ and persistent Na+ currents underlying the generation of membrane resonance. The membrane resonance of immature neurons may contribute to the generation of slow oscillatory activity pattern in the immature neocortex and enhance the temporal precision of synaptic integration in developing cortical neurons.rn

Predicting spike timing of neocortical pyramidal neurons by simple threshold models

Neurons generate spikes reliably with millisecond precision if driven by a fluctuating current--is it then possible to predict the spike timing knowing the input? We determined parameters of an adapting threshold model using data recorded in vitro from 24 layer 5 pyramidal neurons from rat somatosensory cortex, stimulated intracellularly by a fluctuating current simulating synaptic bombardment in vivo. The model generates output spikes whenever the membrane voltage (a filtered version of the input current) reaches a dynamic threshold. We find that for input currents with large fluctuation amplitude, up to 75% of the spike times can be predicted with a precision of +/-2 ms. Some of the intrinsic neuronal unreliability can be accounted for by a noisy threshold mechanism. Our results suggest that, under random current injection into the soma, (i) neuronal behavior in the subthreshold regime can be well approximated by a simple linear filter; and (ii) most of the nonlinearities are captured by a simple threshold process.

Autoinhibition of serotonin cells: An intrinsic regulatory mechanism sensitive to the pattern of usage of the cells

After periods of high-frequency firing, the normal rhythmically active serotonin (5HT)-containing neurosecretory neurons of the lobster ventral nerve cord display a period of suppressed spike generation and reduced synaptic input that we refer to as “autoinhibition.” The duration of this autoinhibition is directly related to the magnitude and duration of the current injection triggering the high-frequency firing. More interesting, however, is that the autoinhibition is inversely related to the initial firing frequency of these cells within their normal range of firing (0.5–3 Hz). This allows more active 5HT neurons to resume firing after shorter durations of inhibition than cells that initially fired at slower rates. Although superfused 5HT inhibits the spontaneous firing of these cells, the persistence of autoinhibition in saline with no added calcium, in cadmium-containing saline, and in lobsters depleted of serotonin suggests that intrinsic membrane properties account for the autoinhibition. A similar autoinhibition is seen in spontaneously active octopamine neurons but is absent from spontaneously active γ-aminobutyric acid cells. Thus, this might be a characteristic feature of amine-containing neurosecretory neurons. The 5HT cells of vertebrate brain nuclei share similarities in firing frequencies, spike shapes, and inhibition by 5HT with the lobster cells that were the focus of this study. However, the mechanism suggested to underlie autoinhibition in vertebrate neurons is that 5HT released from activated or neighboring cells acts back on inhibitory autoreceptors that are found on the dendrites and cell bodies of these neurons.

Democratization in a passive dendritic tree: an analytical investigation

One way to achieve amplification of distal synaptic inputs on a dendritic tree is to scale the amplitude and/or duration of the synaptic conductance with its distance from the soma. This is an example of what is often referred to as “dendritic democracy”. Although well studied experimentally, to date this phenomenon has not been thoroughly explored from a mathematical perspective. In this paper we adopt a passive model of a dendritic tree with distributed excitatory synaptic conductances and analyze a number of key measures of democracy. In particular, via moment methods we derive laws for the transport, from synapse to soma, of strength, characteristic time, and dispersion. These laws lead immediately to synaptic scalings that overcome attenuation with distance. We follow this with a Neumann approximation of Green’s representation that readily produces the synaptic scaling that democratizes the peak somatic voltage response. Results are obtained for both idealized geometries and for the more realistic geometry of a rat CA1 pyramidal cell. For each measure of democratization we produce and contrast the synaptic scaling associated with treating the synapse as either a conductance change or a current injection. We find that our respective scalings agree up to a critical distance from the soma and we reveal how this critical distance decreases with decreasing branch radius.

Neuronal networks with gap junctions: A study of piece-wise linear planar neuron models

The presence of gap junction coupling among neurons of the central nervous systems has been appreciated for some time now. In recent years there has been an upsurge of interest from the mathematical community in understanding the contribution of these direct electrical connections between cells to large-scale brain rhythms. Here we analyze a class of exactly soluble single neuron models, capable of producing realistic action potential shapes, that can be used as the basis for understanding dynamics at the network level. This work focuses on planar piece-wise linear models that can mimic the firing response of several different cell types. Under constant current injection the periodic response and phase response curve (PRC) is calculated in closed form. A simple formula for the stability of a periodic orbit is found using Floquet theory. From the calculated PRC and the periodic orbit a phase interaction function is constructed that allows the investigation of phase-locked network states using the theory of weakly coupled oscillators. For large networks with global gap junction connectivity we develop a theory of strong coupling instabilities of the homogeneous, synchronous and splay state. For a piece-wise linear caricature of the Morris-Lecar model, with oscillations arising from a homoclinic bifurcation, we show that large amplitude oscillations in the mean membrane potential are organized around such unstable orbits.

Mode locking in a spatially extended neuron model: active soma and compartmental tree

Understanding the mode-locked response of excitable systems to periodic forcing has important applications in neuroscience. For example it is known that spatially extended place cells in the hippocampus are driven by the theta rhythm to generate a code conveying information about spatial location. Thus it is important to explore the role of neuronal dendrites in generating the response to periodic current injection. In this paper we pursue this using a compartmental model, with linear dynamics for each compartment, coupled to an active soma model that generates action potentials. By working with the piece-wise linear McKean model for the soma we show how the response of the whole neuron model (soma and dendrites) can be written in closed form. We exploit this to construct a stroboscopic map describing the response of the spatially extended model to periodic forcing. A linear stability analysis of this map, together with a careful treatment of the non-differentiability of the soma model, allows us to construct the Arnol'd tongue structure for 1:q states (one action potential for q cycles of forcing). Importantly we show how the presence of quasi-active membrane in the dendrites can influence the shape of tongues. Direct numerical simulations confirm our theory and further indicate that resonant dendritic membrane can enlarge the windows in parameter space for chaotic behavior. These simulations also show that the spatially extended neuron model responds differently to global as opposed to point forcing. In the former case spatio-temporal patterns of activity within an Arnol'd tongue are standing waves, whilst in the latter they are traveling waves.

Regulation of electrical transmission by protein kinase C in Aplysia bag cell neurons

Electrical synapses are composed of gap junctions, made from paired hemi-channels that allow for the transfer of current from one neuron to another. Gap junctions mediate electrical transmission in neurons, where they synchronize spiking and promote rapid transmission, thereby influencing the coordination, pattern, and frequency of firing. In the marine snail, Aplysia calfornica, two clusters of neuroendocrine bag cell neurons use electrical synapses to synchronize a 30-min burst of action potentials, known as the afterdischarge, which releases egg-laying hormone and induces reproduction. In culture, paired bag cell neurons present a junctional conductance that is non-rectifying and largely voltage-independent. During the afterdischarge, PKC is activated, which is known to increase voltage-gated Ca2+ current; yet, little is understood as to how this pathway impacts electrical transmission. The transfer of presynaptic spike-like waveforms (generated in voltage-clamp) to the postsynaptic cell (measured in current-clamp) was monitored with or without PKC activation. It was found that pretreatment with the PKC activator, phorbol-12-myristate-13-acetate (PMA), enhanced junctional conductance between bag cell neurons. Furthermore, in control, presynaptic action potential waveforms mainly evoked postsynaptic electrotonic potentials at both -60 and -40 mV. However, with PKC activation the presynaptic stimulus consistently elicited postsynaptic action potentials from resting potentials of -40 mV, and would occasionally result in firing from repetitive input at -60 mV. Moreover, to assess whether this enhanced electrical transmission genuinely reflects a greater junctional conductance or a change in postsynaptic responsiveness, a fast-phase junctional-like current was applied to single bag cell neurons. Neurons in PMA always fired action potentials in response to current injection as opposed to control, which were less likely to spike. This outcome did not change when the junctional-like current was artificially enhanced in control conditions. Also, in response to fast- and slow-phase electrotonic potential (ETP) waveforms, Ca2+ current was markedly larger in single PMA-treated neurons. These findings suggest that PKC activation may contribute to afterdischarge fidelity by recruiting postsynaptic Ca2+ current to promote synchronous network firing. Finally, Aplysia gap junction genes (innexins) were transfected into mouse N2A cells and characterized. This revealed a biophysical and pharmacological profile similar to native gap junctions.

Current Sheet Regulation of Solar Near-relativistic Electron Injection Histories

We present a sample of three large near-relativistic (>50 keV) electron events observed in 2001 by both the ACE and the Ulysses spacecraft, when Ulysses was at high-northern latitudes (>60◦) and close to 2 AU. Despite the large latitudinal distance between the two spacecraft, electrons injected near the Sun reached both heliospheric locations. All three events were associated with large solar flares, strong decametric type II radio bursts and accompanied by wide (>212◦) and fast (>1400 km s−1) coronal mass ejections (CMEs). We use advanced interplanetary transport simulations and make use of the directional intensities observed in situ by the spacecraft to infer the electron injection profile close to the Sun and the interplanetary transport conditions at both low and high latitudes. For the three selected events, we find similar interplanetary transport conditions at different heliolatitudes for a given event, with values of the mean free path ranging from 0.04 AU to 0.27 AU. We find differences in the injection profiles inferred for each spacecraft. We investigate the role that sector boundaries of the heliospheric current sheet (HCS) have on determining the characteristics of the electron injection profiles. Extended injection profiles, associated with coronal shocks, are found if the magnetic footpoints of the spacecraft lay in the same magnetic sector as the associated flare, while intermittent sparse injection episodes appear when the spacecraft footpoints are in the opposite sector or a wrap in the HCS bounded the CME structure.

Current Sheet Regulation of Solar Near-Relativistic Electron Injection Histories

We present a sample of three large near-relativistic (>50 keV) electron events observed in 2001 by both the ACE and the Ulysses spacecraft, when Ulysses was at high-northern latitudes (>60°) and close to 2 AU. Despite the large latitudinal distance between the two spacecraft, electrons injected near the Sun reached both heliospheric locations. All three events were associated with large solar flares, strong decametric type II radio bursts and accompanied by wide (>212°) and fast (>1400 km s-1) coronal mass ejections (CMEs). We use advanced interplanetary transport simulations and make use of the directional intensities observed in situ by the spacecraft to infer the electron injection profile close to the Sun and the interplanetary transport conditions at both low and high latitudes. For the three selected events, we find similar interplanetary transport conditions at different heliolatitudes for a given event, with values of the mean free path ranging from 0.04 AU to 0.27 AU. We find differences in the injection profiles inferred for each spacecraft. We investigate the role that sector boundaries of the heliospheric current sheet (HCS) have on determining the characteristics of the electron injection profiles. Extended injection profiles, associated with coronal shocks, are found if the magnetic footpoints of the spacecraft lay in the same magnetic sector as the associated flare, while intermittent sparse injection episodes appear when the spacecraft footpoints are in the opposite sector or a wrap in the HCS bounded the CME structure.

High current-density anodic electro-dissolution in flow-injection systems for the determination of aluminium, copper and zinc in non-ferroalloys by flame atomic absorption spectrometry

An automatic Procedure with a high current-density anodic electrodissolution unit (HDAE) is proposed for the determination of aluminium, copper and zinc in non-ferroalloys by flame atonic absorption spectrometry, based on the direct solid analysis. It consists of solenoid valve-based commutation in a flow-injection system for on-line sample electro-dissolution and calibration with one multi-element standard, an electrolytic cell equipped with two electrodes (a silver needle acts as cathode, and sample as anode), and an intelligent unit. The latter is assembled in a PC-compatible microcomputer for instrument control, and far data acquisition and processing. General management of the process is achieved by use of software written in Pascal. Electrolyte compositions, flow rates, commutation times, applied current and electrolysis time mere investigated. A 0.5 mol l(-1) HNO3 solution was elected as electrolyte and 300 A/cm(2) as the continuous current pulse. The performance of the proposed system was evaluated by analysing aluminium in Al-allay samples, and copper/zinc in brass and bronze samples, respectively. The system handles about 50 samples per hour. Results are precise (R.S.D < 2%) and in agreement with those obtained by ICP-AES and spectrophotometry at a 95% confidence level.

Evaluation of current underground injection of industrial waste in Illinois : final draft report /

"HWRIC RR 008." -- Cover.