Simple, Fast and Accurate Implementation of the Diffusion Approximation Algorithm for Stochastic Ion Channels with Multiple States

The phenomena that emerge from the interaction of the stochastic opening and closing of ion channels (channel noise) with the non-linear neural dynamics are essential to our understanding of the operation of the nervous system. The effects that channel noise can have on neural dynamics are generally studied using numerical simulations of stochastic models. Algorithms based on discrete Markov Chains (MC) seem to be the most reliable and trustworthy, but even optimized algorithms come with a non-negligible computational cost. Diffusion Approximation (DA) methods use Stochastic Differential Equations (SDE) to approximate the behavior of a number of MCs, considerably speeding up simulation times. However, model comparisons have suggested that DA methods did not lead to the same results as in MC modeling in terms of channel noise statistics and effects on excitability. Recently, it was shown that the difference arose because MCs were modeled with coupled activation subunits, while the DA was modeled using uncoupled activation subunits. Implementations of DA with coupled subunits, in the context of a specific kinetic scheme, yielded similar results to MC. However, it remained unclear how to generalize these implementations to different kinetic schemes, or whether they were faster than MC algorithms. Additionally, a steady state approximation was used for the stochastic terms, which, as we show here, can introduce significant inaccuracies. We derived the SDE explicitly for any given ion channel kinetic scheme. The resulting generic equations were surprisingly simple and interpretable - allowing an easy and efficient DA implementation. The algorithm was tested in a voltage clamp simulation and in two different current clamp simulations, yielding the same results as MC modeling. Also, the simulation efficiency of this DA method demonstrated considerable superiority over MC methods.Comment: 32 text pages, 10 figures, 1 supplementary text + figur

Brain energy metabolism: A roadmap for future research

Although we have learned much about how the brain fuels its functions over the last decades, there remains much still to discover in an organ that is so complex. This article lays out major gaps in our knowledge of interrelationships between brain metabolism and brain function, including biochemical, cellular, and subcellular aspects of functional metabolism and its imaging in adult brain, as well as during development, aging, and disease. The focus is on unknowns in metabolism of major brain substrates and associated transporters, the roles of insulin and of lipid droplets, the emerging role of metabolism in microglia, mysteries about the major brain cofactor and signaling molecule NAD+, as well as unsolved problems underlying brain metabolism in pathologies such as traumatic brain injury, epilepsy, and metabolic downregulation during hibernation. It describes our current level of understanding of these facets of brain energy metabolism as well as a roadmap for future research

Detection of the PAX3-FKHR fusion gene in paediatric rhabdomyosarcoma: a reproducible predictor of outcome?

Rhabdomyosarcoma has 2 major histological subtypes, embryonal and alveolar. Alveolar histology is associated with the fusion genes PAX3-FKHR and PAX7-FKHR. Definition of alveolar has been complicated by changes in terminology and subjectivity. It is currently unclear whether adverse clinical behaviour is better predicted by the presence of these fusion genes or by alveolar histology. We have determined the presence of the PAX3/7-FKHR fusion genes in 91 primary rhabdomyosarcoma tumours using a combination of classical cytogenetics, FISH and RT-PCR, with a view to determining the clinical characteristics of tumours with and without the characteristic translocations. There were 37 patients with t(2;13)/PAX3-FKHR, 8 with t(1;13) PAX7-FKHR and 46 with neither translocation. One or other of the characteristic translocations was found in 31/38 (82%) of alveolar cases. Univariate survival analysis revealed the presence of the translocation t(2;13)/PAX3-FKHR to be an adverse prognostic factor. With the difficulties in morphological diagnosis of alveolar rhabdomyosarcoma on increasingly used small needle biopsy specimens, these data suggest that molecular analysis for PAX3-FKHR will be a clinically useful tool in treatment stratification in the future. This hypothesis requires testing in a prospective study. Variant t(1;13)/PAX7-FKHR appears biologically different, occurring in younger patients with more localised disease. © 2001 Cancer Research Campaignhttp://www.bjcancer.co

Fin-Tail Coordination during Escape and Predatory Behavior in Larval Zebrafish

Larval zebrafish innately perform a suite of behaviors that are tightly linked to their evolutionary past, notably escape from threatening stimuli and pursuit and capture of prey. These behaviors have been carefully examined in the past, but mostly with regard to the movements of the trunk and tail of the larvae. Here, we employ kinematics analyses to describe the movements of the pectoral fins during escape and predatory behavior. In accord with previous studies, we find roles for the pectoral fins in slow swimming and immediately after striking prey. We find novel roles for the pectoral fins in long-latency, but not in short-latency C-bends. We also observe fin movements that occur during orienting J-turns and S-starts that drive high-velocity predatory strikes. Finally, we find that the use of pectoral fins following a predatory strike is scaled to the velocity of the strike, supporting a role for the fins in braking. The implications of these results for central control of coordinated movements are discussed, and we hope that these results will provide baselines for future analyses of cross-body coordination using mutants, morphants, and transgenic approaches

Stochastic Ion Channel Gating in Dendritic Neurons: Morphology Dependence and Probabilistic Synaptic Activation of Dendritic Spikes

Neuronal activity is mediated through changes in the probability of stochastic transitions between open and closed states of ion channels. While differences in morphology define neuronal cell types and may underlie neurological disorders, very little is known about influences of stochastic ion channel gating in neurons with complex morphology. We introduce and validate new computational tools that enable efficient generation and simulation of models containing stochastic ion channels distributed across dendritic and axonal membranes. Comparison of five morphologically distinct neuronal cell types reveals that when all simulated neurons contain identical densities of stochastic ion channels, the amplitude of stochastic membrane potential fluctuations differs between cell types and depends on sub-cellular location. For typical neurons, the amplitude of membrane potential fluctuations depends on channel kinetics as well as open probability. Using a detailed model of a hippocampal CA1 pyramidal neuron, we show that when intrinsic ion channels gate stochastically, the probability of initiation of dendritic or somatic spikes by dendritic synaptic input varies continuously between zero and one, whereas when ion channels gate deterministically, the probability is either zero or one. At physiological firing rates, stochastic gating of dendritic ion channels almost completely accounts for probabilistic somatic and dendritic spikes generated by the fully stochastic model. These results suggest that the consequences of stochastic ion channel gating differ globally between neuronal cell-types and locally between neuronal compartments. Whereas dendritic neurons are often assumed to behave deterministically, our simulations suggest that a direct consequence of stochastic gating of intrinsic ion channels is that spike output may instead be a probabilistic function of patterns of synaptic input to dendrites

Opto-Current-Clamp Actuation of Cortical Neurons Using a Strategically Designed Channelrhodopsin

BACKGROUND: Optogenetic manipulation of a neuronal network enables one to reveal how high-order functions emerge in the central nervous system. One of the Chlamydomonas rhodopsins, channelrhodopsin-1 (ChR1), has several advantages over channelrhodopsin-2 (ChR2) in terms of the photocurrent kinetics. Improved temporal resolution would be expected by the optogenetics using the ChR1 variants with enhanced photocurrents. METHODOLOGY/PRINCIPAL FINDINGS: The photocurrent retardation of ChR1 was overcome by exchanging the sixth helix domain with its counterpart in ChR2 producing Channelrhodopsin-green receiver (ChRGR) with further reform of the molecule. When the ChRGR photocurrent was measured from the expressing HEK293 cells under whole-cell patch clamp, it was preferentially activated by green light and has fast kinetics with minimal desensitization. With its kinetic advantages the use of ChRGR would enable one to inject a current into a neuron by the time course as predicted by the intensity of the shedding light (opto-current clamp). The ChRGR was also expressed in the motor cortical neurons of a mouse using Sindbis pseudovirion vectors. When an oscillatory LED light signal was applied sweeping through frequencies, it robustly evoked action potentials synchronized to the oscillatory light at 5-10 Hz in layer 5 pyramidal cells in the cortical slice. The ChRGR-expressing neurons were also driven in vivo with monitoring local field potentials (LFPs) and the time-frequency energy distribution of the light-evoked response was investigated using wavelet analysis. The oscillatory light enhanced both the in-phase and out-phase responses of LFP at the preferential frequencies of 5-10 Hz. The spread of activity was evidenced by the fact that there were many c-Fos-immunoreactive neurons that were negative for ChRGR in a region of the motor cortex. CONCLUSIONS/SIGNIFICANCE: The opto-current-clamp study suggests that the depolarization of a small number of neurons wakes up the motor cortical network over some critical point to the activated state

The Melbourne epidemic thunderstorm asthma event 2016: an investigation of environmental triggers, effect on health services, and patient risk factors.

BACKGROUND: A multidisciplinary collaboration investigated the world's largest, most catastrophic epidemic thunderstorm asthma event that took place in Melbourne, Australia, on Nov 21, 2016, to inform mechanisms and preventive strategies. METHODS: Meteorological and airborne pollen data, satellite-derived vegetation index, ambulance callouts, emergency department presentations, and data on hospital admissions for Nov 21, 2016, as well as leading up to and following the event were collected between Nov 21, 2016, and March 31, 2017, and analysed. We contacted patients who presented during the epidemic thunderstorm asthma event at eight metropolitan health services (each including up to three hospitals) via telephone questionnaire to determine patient characteristics, and investigated outcomes of intensive care unit (ICU) admissions. FINDINGS: Grass pollen concentrations on Nov 21, 2016, were extremely high (>100 grains/m3). At 1800 AEDT, a gust front crossed Melbourne, plunging temperatures 10°C, raising humidity above 70%, and concentrating particulate matter. Within 30 h, there were 3365 (672%) excess respiratory-related presentations to emergency departments, and 476 (992%) excess asthma-related admissions to hospital, especially individuals of Indian or Sri Lankan birth (10% vs 1%, p<0·0001) and south-east Asian birth (8% vs 1%, p<0·0001) compared with previous 3 years. Questionnaire data from 1435 (64%) of 2248 emergency department presentations showed a mean age of 32·0 years (SD 18·6), 56% of whom were male. Only 28% had current doctor-diagnosed asthma. 39% of the presentations were of Asian or Indian ethnicity (25% of the Melbourne population were of this ethnicity according to the 2016 census, relative risk [RR] 1·93, 95% CI 1·74-2·15, p <0·0001). Of ten individuals who died, six were Asian or Indian (RR 4·54, 95% CI 1·28-16·09; p=0·01). 35 individuals were admitted to an intensive care unit, all had asthma, 12 took inhaled preventers, and five died. INTERPRETATION: Convergent environmental factors triggered a thunderstorm asthma epidemic of unprecedented magnitude, tempo, and geographical range and severity on Nov 21, 2016, creating a new benchmark for emergency and health service escalation. Asian or Indian ethnicity and current doctor-diagnosed asthma portended life-threatening exacerbations such as those requiring admission to an ICU. Overall, the findings provide important public health lessons applicable to future event forecasting, health care response coordination, protection of at-risk populations, and medical management of epidemic thunderstorm asthma. FUNDING: None

The Immunological Synapse: a Dynamic Platform for Local Signaling

The immunological synapse (IS) as a concept has evolved from a static view of the junction between T cells and their antigen-presenting cell partners. The entire process of IS formation and extinction is now known to entail a dynamic reorganization of membrane domains and proteins within and adjacent to those domains. Discussion The entire process is also intricately tied to the motility machinery—both as that machinery directs “scanning” prior to T-cell receptor engagement and as it is appropriated during the ongoing developments at the IS. While the synapse often remains dynamic in order to encourage surveillance of new antigen-presenting surfaces, cytoskeletal forces also regulate the development of signals, likely including the assembly of ion channels. In both neuronal and immunological synapses, localized Ca 2+ signals and accumulation or depletion of ions in microdomains accompany the concentration of signaling molecules in the synapse. Such spatiotemporal signaling in the synapse greatly accelerates kinetics and provides essential checkpoints to validate effective cell–cell communication