Search

The human brain operates at multiple levels, from molecules to circuits, and understanding these complex processes requires integrated research efforts. Simulating biophysically-detailed neuron models is a computationally expensive but effective method for studying local neural circuits. Recent innovations have shown that artificial neural networks (ANNs) can accurately predict the behavior of these detailed models in terms of spikes, electrical potentials, and optical readouts. While these methods have the potential to accelerate large network simulations by several orders of magnitude compared to conventional differential equation based modelling, they currently only predict voltage outputs for the soma or a select few neuron compartments. Our novel approach, based on enhanced state-of-the-art architectures for multitask learning (MTL), allows for the simultaneous prediction of membrane potentials in each compartment of a neuron model, at a speed of up to two orders of magnitude faster than classical simulation methods. By predicting all membrane potentials together, our approach not only allows for comparison of model output with a wider range of experimental recordings (patch-electrode, voltage-sensitive dye imaging), it also provides the first stepping stone towards predicting local field potentials (LFPs), electroencephalogram (EEG) signals, and magnetoencephalography (MEG) signals from ANN-based simulations. While LFP and EEG are an important downstream application, the main focus of this paper lies in predicting dendritic voltages within each compartment to capture the entire electrophysiology of a biophysically-detailed neuron model. It further presents a challenging benchmark for MTL architectures due to the large amount of data involved, the presence of correlations between neighbouring compartments, and the non-Gaussian distribution of membrane potentials. Author summary: Our research focuses on cutting-edge techniques in computational neuroscience. We specifically make use of simulations of biophysically detailed neuron models. Traditionally these methods are computationally intensive, but recent advancements using artificial neural networks (ANNs) have shown promise in predicting neural behavior with remarkable accuracy. However, existing ANNs fall short in providing comprehensive predictions across all compartments of a neuron model and only provide information on the activity of a limited number of locations along the extent of a neuron. In our study, we introduce a novel approach leveraging state-of-the-art multitask learning architectures. This approach allows us to simultaneously predict membrane potentials in every compartment of a neuron model. By distilling the underlying electrophysiology into an ANN, we significantly outpace conventional simulation methods. By accurately capturing voltage outputs across the neuron's structure, our method invites comparisons with experimental data and paves the way for predicting complex aggregate signals such as local field potentials and EEG signals. Our findings not only advance our understanding of neural dynamics but also present a significant benchmark for future research in computational neuroscience. [ABSTRACT FROM AUTHOR]

Hippocampal ripple oscillations have been implicated in important cognitive functions such as memory consolidation and planning. Multiple computational models have been proposed to explain the emergence of ripple oscillations, relying either on excitation or inhibition as the main pacemaker. Nevertheless, the generating mechanism of ripples remains unclear. An interesting dynamical feature of experimentally measured ripples, which may advance model selection, is intra-ripple frequency accommodation (IFA): a decay of the instantaneous ripple frequency over the course of a ripple event. So far, only a feedback-based inhibition-first model, which relies on delayed inhibitory synaptic coupling, has been shown to reproduce IFA. Here we use an analytical mean-field approach and numerical simulations of a leaky integrate-and-fire spiking network to explain the mechanism of IFA. We develop a drift-based approximation for the oscillation dynamics of the population rate and the mean membrane potential of interneurons under strong excitatory drive and strong inhibitory coupling. For IFA, the speed at which the excitatory drive changes is critical. We demonstrate that IFA arises due to a speed-dependent hysteresis effect in the dynamics of the mean membrane potential, when the interneurons receive transient, sharp wave-associated excitation. We thus predict that the IFA asymmetry vanishes in the limit of slowly changing drive, but is otherwise a robust feature of the feedback-based inhibition-first ripple model. Author summary: The hippocampus plays a central role in the acquisition and consolidation of explicit memories. During sleep or rest, hippocampal neurons replay recently acquired memories, while the neuronal network exhibits high-frequency oscillations, so-called ripples. To study the potential function of ripples, their generation mechanism needs to be clarified. Here we analyze a model network of inhibitory interneurons as a potential ripple pacemaker. We derive an analytical approximation of the ripple dynamics that explains why the ripple frequency tends to decay over the course of an event (intra-ripple frequency accommodation, IFA). Studying dynamical phenomena such as IFA advances model selection and enables an understanding of how rhythmic brain activity contributes to cognitive functions. [ABSTRACT FROM AUTHOR]

Merkel cells combine with Aβ afferents, producing slowly adapting type 1(SA1) responses to mechanical stimuli. However, how Merkel cells transduce mechanical stimuli into neural signals to Aβ afferents is still unclear. Here we develop a biophysical model of Merkel cells for mechanical transduction by incorporating main ingredients such as Ca2+ and K+ voltage-gated channels, Piezo2 channels, internal Ca2+ stores, neurotransmitters release, and cell deformation. We first validate our model with several experiments. Then we reveal that Ca2+ and K+ channels on the plasma membrane shape the depolarization of membrane potentials, further regulating the Ca2+ transients in the cells. We also show that Ca2+ channels on the plasma membrane mainly inspire the Ca2+ transients, while internal Ca2+ stores mainly maintain the Ca2+ transients. Moreover, we show that though Piezo2 channels are rapidly adapting mechanical-sensitive channels, they are sufficient to inspire sustained Ca2+ transients in Merkel cells, which further induce the release of neurotransmitters for tens of seconds. Thus our work provides a model that captures the membrane potentials and Ca2+ transients features of Merkel cells and partly explains how Merkel cells transduce the mechanical stimuli by Piezo2 channels. Author summary: Touch is an essential way for humans to sense the physical world. It is necessary to figure out how our tactile system works. However, how Merkel cells convey mechanical stimuli into neural signals and deliver them to Aβ afferents is still poorly understood. In this work, we develop a biophysical model of Merkel cells for mechanical transduction. We show that Ca2+ and K+ channels on the plasma membrane control the membrane potentials of Merkel cells and further regulate the Ca2+ transients in the cells. Ca2+ channels on ER and MT directly contribute to the Ca2+ transients. Under indentation, the influx of Ca2+ through Piezo2 channels is sufficient to trigger the Ca2+ transients in Merkel cells and internal Ca2+ stores inherit to maintain the Ca2+ transients, which further result in a sustained release of neurotransmitters that activates the Aβ afferents. Thus our findings partly explain how Merkel cells combine with Aβ afferents to generate an SA1 response to sustained indentation by Piezo2 channels. [ABSTRACT FROM AUTHOR]

Multi-scale mathematical bioelectrical models of organs such as the uterus, stomach or heart present challenges both for accuracy and computational tractability. These multi-scale models are typically founded on models of biological cells derived from the classic Hodkgin-Huxley (HH) formalism. Ion channel behaviour is tracked with dynamical variables representing activation or inactivation of currents that relax to steady-state dependencies on cellular membrane voltage. Timescales for relaxation may be orders of magnitude faster than companion ion channel variables or phenomena of physiological interest for the entire cell (such as bursting sequences of action potentials) or the entire organ (such as electromechanical coordination). Exploiting these time scales with steady-state approximations for relatively fast-acting systems is a well-known but often overlooked approach as evidenced by recent published models. We thus investigate feasibility of an extensive reduction of order for an HH-type cell model with steady-state approximations to the full dynamical activation and inactivation ion channel variables. Our effort utilises a published comprehensive uterine smooth muscle cell model that encompasses 19 ordinary differential equations and 105 formulations overall. The numerous ion channel submodels in the published model exhibit relaxation times ranging from order 10−1 to 105 milliseconds. Substitution of the faster dynamic variables with steady-state formulations demonstrates both an accurate reproduction of the full model and substantial improvements in time-to-solve, for test cases performed. Our demonstration here of an effective and relatively straightforward reduction method underlines the particular importance of considering time scales for model simplification before embarking on large-scale computations or parameter sweeps. As a preliminary complement to more intensive reduction of order methods such as parameter sensitivity and bifurcation analysis, this approach can rapidly and accurately improve computational tractability for challenging multi-scale organ modelling efforts. Author summary: Mathematical modeling of physiological organ systems encompassing intracellular to organ-wide behaviour grapple with intrinsically multi-scaled systems in both space and time. Significant computational challenges arise with their numerical solution, often substantially constraining feasibility of in silico investigations. Naturally, streamlining and reducing the mathematical complexity of these models while maintaining accurate reproduction of experimental data is a persistent concern. We present a straightforward method for exploiting multiple time scales reducing the order of models by way of steady-state approximations at the cellular level. Applied to a uterine smooth muscle cell model, we obtain substantial improvements in time to solve by around a factor of two while accurately reproducing full model results. Our successful demonstration here highlights the importance and relative ease of the method—that is not altogether unknown but often overlooked—with potential application across a wide variety of multi-scale organ modeling efforts. [ABSTRACT FROM AUTHOR]

When a neuron is driven beyond its threshold, it spikes. The fact that it does not communicate its continuous membrane potential is usually seen as a computational liability. Here we show that this spiking mechanism allows neurons to produce an unbiased estimate of their causal influence, and a way of approximating gradient descent-based learning. Importantly, neither activity of upstream neurons, which act as confounders, nor downstream non-linearities bias the results. We show how spiking enables neurons to solve causal estimation problems and that local plasticity can approximate gradient descent using spike discontinuity learning. Author summary: Despite significant research, models of spiking neural networks still lag behind artificial neural networks in terms of performance in machine learning and modeling cognitive tasks. Given this, we may wonder, why do neurons spike? A key problem that must be solved in any neural network is the credit assignment problem. That is, how does a neuron know its effect on downstream computation and rewards, and thus how it should change its synaptic weights to improve? Artificial neural networks solve this problem with the back-propagation algorithm. We are still seeking to understand how biological neural networks effectively solve this problem. In this work we show that the discontinuous, all-or-none spiking response of a neuron can in fact be used to estimate a neuron's causal effect on downstream processes. Inspired by methods from econometrics, we show that the thresholded response of a neuron can be used to get at that neuron's unique contribution to a reward signal, separating it from other neurons whose activity it may be correlated with. This proposal provides insights into a novel function of spiking that we explore in simple networks and learning tasks. [ABSTRACT FROM AUTHOR]

Functional forms of biophysically-realistic neuron models are constrained by neurobiological and anatomical considerations, such as cell morphologies and the presence of known ion channels. Despite these constraints, neuron models still contain unknown static parameters which must be inferred from experiment. This inference task is most readily cast into the framework of state-space models, which systematically takes into account partial observability and measurement noise. Inferring only dynamical state variables such as membrane voltages is a well-studied problem, and has been approached with a wide range of techniques beginning with the well-known Kalman filter. Inferring both states and fixed parameters, on the other hand, is less straightforward. Here, we develop a method for joint parameter and state inference that combines traditional state space modeling with chaotic synchronization and optimal control. Our methods are tailored particularly to situations with considerable measurement noise, sparse observability, very nonlinear or chaotic dynamics, and highly uninformed priors. We illustrate our approach both in a canonical chaotic model and in a phenomenological neuron model, showing that many unknown parameters can be uncovered reliably and accurately from short and noisy observed time traces. Our method holds promise for estimation in larger-scale systems, given ongoing improvements in calcium reporters and genetically-encoded voltage indicators. Author Summary: Systems neuroscience aims to understand how individual neurons and neural networks process external stimuli into behavioral responses. Underlying this characterization are mathematical models intimately shaped by experimental observations. But neural systems are high-dimensional and contain highly nonlinear interactions, so developing accurate models remains a challenge given current experimental capabilities. In practice, this means that the dynamical equations characterizing neural activity have many unknown parameters, and these parameters must be inferred from data. This inference problem is nontrivial owing to model nonlinearity, system and measurement noise, and the sparsity of observations from electrode recordings. Here, we present a novel method for inferring model parameters of neural systems. Our technique combines ideas from control theory and optimization, and amounts to using data to "control" estimates toward the best fit. Our method compares well in accuracy against other state-of-the-art inference methods, both in phenomenological chaotic systems and biophysical neuron models. Our work shows that many unknown model parameters of interest can be inferred from voltage measurements, despite signaling noise, instrument noise, and low observability. [ABSTRACT FROM AUTHOR]

Identifying plausible mechanisms for compartmentalization and accumulation of the organic intermediates of early metabolic cycles in primitive cells has been a major challenge in theories of life's origins. Here, we propose a mechanism, where positive membrane potentials elevate the concentration of the organic intermediates. Positive membrane potentials are generated by positively charged surfaces of protocell membranes due to accumulation of transition metals. We find that (i) positive membrane potentials comparable in magnitude to those of modern cells can increase the concentration of the organic intermediates by several orders of magnitude; (ii) generation of large membrane potentials destabilize ion distributions; (iii) violation of electroneutrality is necessary to induce nonzero membrane potentials; and (iv) violation of electroneutrality enhances osmotic pressure and diminishes reaction efficiency, resulting in an evolutionary driving force for the formation of lipid membranes, specialized ion channels, and active transport systems. Author summary: The building blocks of life are constantly synthesized and broken down through concurrent cycles of chemical transformations. Tracing these reactions back 4 billion years to their origins has been a long-standing goal of evolutionary biology. The first metabolic cycles at the origin of life must have overcome several obstacles to spontaneously start and sustain their nonequilibrium states. Notably, maintaining the concentration of organic intermediates at high levels to support their continued operation and evolution would have been challenging in abiotic cells lacking evolutionarily tuned lipid membranes and enzymes. We propose a mechanism, by which the concentrations of organic intermediates are elevated to drive early metabolic cycles in abiotic cells and demonstrate its feasibility in a protocell model from first principles. [ABSTRACT FROM AUTHOR]

Resonance is defined as maximal response of a system to periodic inputs in a limited frequency band. Resonance may serve to optimize inter-neuronal communication, and has been observed at multiple levels of neuronal organization. However, it is unknown how neuronal resonance observed at the network level is generated and how network resonance depends on the properties of the network building blocks. Here, we first develop a metric for quantifying spike timing resonance in the presence of background noise, extending the notion of spiking resonance for in vivo experiments. Using conductance-based models, we find that network resonance can be inherited from resonances at other levels of organization, or be intrinsically generated by combining mechanisms across distinct levels. Resonance of membrane potential fluctuations, postsynaptic potentials, and single neuron spiking can each be generated independently of resonance at any other level and be propagated to the network level. At all levels of organization, interactions between processes that give rise to low- and high-pass filters generate the observed resonance. Intrinsic network resonance can be generated by the combination of filters belonging to different levels of organization. Inhibition-induced network resonance can emerge by inheritance from resonance of membrane potential fluctuations, and be sharpened by presynaptic high-pass filtering. Our results demonstrate a multiplicity of qualitatively different mechanisms that can generate resonance in neuronal systems, and provide analysis tools and a conceptual framework for the mechanistic investigation of network resonance in terms of circuit components, across levels of neuronal organization. Author summary: How one part of the brain responds to periodic input from another part depends on resonant circuit properties. Resonance is a basic property of physical systems, and has been experimentally observed at various levels of neuronal organization both in vitro and in vivo. Yet how resonance is generated in neuronal networks is largely unknown. In particular, whether resonance can be generated directly at the level of a network of spiking neurons remains to be determined. Using detailed biophysical modeling, we develop a conceptual framework according to which resonance at a given level of organization is generated by the interplay of low- and high-pass filters, implemented at either the same or across levels of neuronal organization. We tease apart representative, biophysically-plausible generative mechanisms of resonance at four different levels of organization: membrane potential fluctuations, single neuron spiking, synaptic transmission, and neuronal networks. We identify conditions under which resonance at one level can be inherited to another level of organization, provide conclusive evidence that resonance at each level can be generated without resonance at any other level, and describe a number of representative routes to network resonance. The proposed framework facilitates the investigation of resonance in neuronal systems. [ABSTRACT FROM AUTHOR]

Rhythmic neuronal network activity underlies brain oscillations. To investigate how connected neuronal networks contribute to the emergence of the α-band and to the regulation of Up and Down states, we study a model based on synaptic short-term depression-facilitation with afterhyperpolarization (AHP). We found that the α-band is generated by the network behavior near the attractor of the Up-state. Coupling inhibitory and excitatory networks by reciprocal connections leads to the emergence of a stable α-band during the Up states, as reflected in the spectrogram. To better characterize the emergence and stability of thalamocortical oscillations containing α and δ rhythms during anesthesia, we model the interaction of two excitatory networks with one inhibitory network, showing that this minimal topology underlies the generation of a persistent α-band in the neuronal voltage characterized by dominant Up over Down states. Finally, we show that the emergence of the α-band appears when external inputs are suppressed, while fragmentation occurs at small synaptic noise or with increasing inhibitory inputs. To conclude, α-oscillations could result from the synaptic dynamics of interacting excitatory neuronal networks with and without AHP, a principle that could apply to other rhythms. Author summary: Brain oscillations, recorded from electroencephalograms characterize behaviors such as sleep, wakefulness, brain evoked responses, coma or anesthesia. The underlying rhythms for these oscillations are associated at a neuronal population level to fluctuations of the membrane potential between Up (depolarized) and Down (hyperpolarized) states. During anesthesia with propofol, a dominant α-band (8–12Hz) can emerge or disappear, but the underlying mechanism remains unclear. Using modeling, we report that the α-band appears during Up states in neuronal populations driven by short-term synaptic plasticity and synaptic noise. Moreover, we show that three connected neuronal networks representing the thalamocortical loop reproduce the dynamics of the α-band, which emerges following the arrest of excitatory stimulations, but that can disappear by increasing inhibitory inputs. To conclude, short-term plasticity in well connected neuronal networks can explain the emergence and fragmentation of the α-band. [ABSTRACT FROM AUTHOR]

Mutations are known to cause perturbations in essential functional features of integral membrane proteins, including ion channels. Even restricted or point mutations can result in substantially changed properties of ion currents. The additive effect of these alterations for a specific ion channel can result in significantly changed properties of the action potential (AP). Both AP shortening and AP prolongation can result from known mutations, and the consequences can be life-threatening. Here, we present a computational method for identifying new drugs utilizing combinations of existing drugs. Based on the knowledge of theoretical effects of existing drugs on individual ion currents, our aim is to compute optimal combinations that can 'repair' the mutant AP waveforms so that the baseline AP-properties are restored. More specifically, we compute optimal, combined, drug concentrations such that the waveforms of the transmembrane potential and the cytosolic calcium concentration of the mutant cardiomyocytes (CMs) becomes as similar as possible to their wild type counterparts after the drug has been applied. In order to demonstrate the utility of this method, we address the question of computing an optimal drug for the short QT syndrome type 1 (SQT1). For the SQT1 mutation N588K, there are available data sets that describe the effect of various drugs on the mutated K+ channel. These published findings are the basis for our computational analysis which can identify optimal compounds in the sense that the AP of the mutant CMs resembles essential biomarkers of the wild type CMs. Using recently developed insights regarding electrophysiological properties among myocytes from different species, we compute optimal drug combinations for hiPSC-CMs, rabbit ventricular CMs and adult human ventricular CMs with the SQT1 mutation. Since the 'composition' of ion channels that form the AP is different for the three types of myocytes under consideration, so is the composition of the optimal drug. Author summary: Poly-pharmacology (using multiple drugs to treat disease) has been proposed for improving cardiac anti-arrhythmic therapy for at least two decades. However, the specific arrhythmia contexts in which polytherapy is likely to be both safe and effective have remained elusive. Type 1 short QT syndrome (SQT1) is a rare form of cardiac arrhythmia that results from mutations to the human Ether-á-go-go Related Gene (hERG) potassium channel. Functionally, these mutations are remarkably consistent in that they permit the channel to open earlier during each heart beat. While hundreds of compounds are known to inhibit hERG channels, the specific effect of SQT1 mutations that allows for early channel opening also limits the ability of most of those compounds to correct SQT1 dysfunction. Here, we have applied a suite of ventricular cardiomyocyte computational models to ask whether polytherapy may offer a more effective therapeutic strategy in SQT1, and if so, what the likely characteristics of that strategy are. Our analyses suggest that simultaneous induction of late sodium current and partial hERG blockade offers a promising strategy. While no activators of late sodium current have been clinically approved, several experimental compounds are available and may provide a basis for interrogating this strategy. The method presented here can be used to compute optimal drug combinations provided that the effect of each drug on every relevant ion channel is known. [ABSTRACT FROM AUTHOR]

Lung cancer is still a leading cause of death worldwide. In recent years, knowledge has been obtained of the mechanisms modulating ion channel kinetics and thus of cell bioelectric properties, which is promising for oncological biomarkers and targets. The complex interplay of channel expression and its consequences on malignant processes, however, is still insufficiently understood. We here introduce the first approach of an in-silico whole-cell ion current model of a cancer cell, in particular of the A549 human lung adenocarcinoma, including the main functionally expressed ion channels in the plasma membrane as so far known. This hidden Markov-based model represents the electrophysiology behind proliferation of the A549 cell, describing its rhythmic oscillation of the membrane potential able to trigger the transition between cell cycle phases, and it predicts membrane potential changes over the cell cycle provoked by targeted ion channel modulation. This first A549 in-silico cell model opens up a deeper insight and understanding of possible ion channel interactions in tumor development and progression, and is a valuable tool for simulating altered ion channel function in lung cancer electrophysiology. Author summary: Advances in the understanding of functional alterations at genetic, epigenetic or protein expression and the expanding knowledge in mechanisms modulating ion channel kinetics and thus the cells' bioelectric properties have arisen as promising cancer biomarkers and oncological targets. Our hidden Markov-based in-silico cell model represents the electrophysiology behind proliferation of the A549 cell line, explaining the cell's rhythmic oscillation from hyperpolarized to depolarized states of the membrane potential, able to trigger the transition between cell cycle phases. The model enables the prediction of membrane potential changes over the cell cycle provoked by targeted modulation of specific ion channels, leading to cell cycle promotion or interruption. We are encouraged that the availability of this first cancer cell model will provide profound insight into possible roles and interactions of ion channels in tumor development and progression, and may aid in the testing of research hypotheses in lung cancer electrophysiology. [ABSTRACT FROM AUTHOR]

The anatomical and functional organization of neurons and astrocytes at 'tripartite synapses' is essential for reliable neurotransmission, which critically depends on ATP. In low energy conditions, synaptic transmission fails, accompanied by a breakdown of ion gradients, changes in membrane potentials and cell swelling. The resulting cellular damage and cell death are causal to the often devastating consequences of an ischemic stroke. The severity of ischemic damage depends on the age and the brain region in which a stroke occurs, but the reasons for this differential vulnerability are far from understood. In the present study, we address this question by developing a comprehensive biophysical model of a glutamatergic synapse to identify key determinants of synaptic failure during energy deprivation. Our model is based on fundamental biophysical principles, includes dynamics of the most relevant ions, i.e., Na+, K+, Ca2+, Cl− and glutamate, and is calibrated with experimental data. It confirms the critical role of the Na+/K+-ATPase in maintaining ion gradients, membrane potentials and cell volumes. Our simulations demonstrate that the system exhibits two stable states, one physiological and one pathological. During energy deprivation, the physiological state may disappear, forcing a transit to the pathological state, which can be reverted when blocking voltage-gated Na+ and K+ channels. Our model predicts that the transition to the pathological state is favoured if the extracellular space fraction is small. A reduction in the extracellular space volume fraction, as, e.g. observed with ageing, will thus promote the brain's susceptibility to ischemic damage. Our work provides new insights into the brain's ability to recover from energy deprivation, with translational relevance for diagnosis and treatment of ischemic strokes. Author summary: The brain consumes energy to keep ion concentrations at normal working conditions. In the case of energy deprivation (ED), e.g., during a stroke, synaptic communication fails first. Inspired by our recent experimental work on ED, we formulated a novel computational model to explore initial events during ED. Our model reproduces time courses for several ions from different experimental data. In some cases, the system returns to baseline upon restoring energy supply. In others, we observe that neurons and astrocytes cannot recover accompanied by cell swelling. There is a threshold depending on the depth and duration of ATP depletion differentiating these cases. Also, smaller extracellular spaces hamper recovery more. This result may explain clinical observations of increased vulnerability to stroke as the size of the extracellular space shrinks with ageing. [ABSTRACT FROM AUTHOR]

The electroencephalogram (EEG) is a major tool for non-invasively studying brain function and dysfunction. Comparing experimentally recorded EEGs with neural network models is important to better interpret EEGs in terms of neural mechanisms. Most current neural network models use networks of simple point neurons. They capture important properties of cortical dynamics, and are numerically or analytically tractable. However, point neurons cannot generate an EEG, as EEG generation requires spatially separated transmembrane currents. Here, we explored how to compute an accurate approximation of a rodent's EEG with quantities defined in point-neuron network models. We constructed different approximations (or proxies) of the EEG signal that can be computed from networks of leaky integrate-and-fire (LIF) point neurons, such as firing rates, membrane potentials, and combinations of synaptic currents. We then evaluated how well each proxy reconstructed a ground-truth EEG obtained when the synaptic currents of the LIF model network were fed into a three-dimensional network model of multicompartmental neurons with realistic morphologies. Proxies based on linear combinations of AMPA and GABA currents performed better than proxies based on firing rates or membrane potentials. A new class of proxies, based on an optimized linear combination of time-shifted AMPA and GABA currents, provided the most accurate estimate of the EEG over a wide range of network states. The new linear proxies explained 85–95% of the variance of the ground-truth EEG for a wide range of network configurations including different cell morphologies, distributions of presynaptic inputs, positions of the recording electrode, and spatial extensions of the network. Non-linear EEG proxies using a convolutional neural network (CNN) on synaptic currents increased proxy performance by a further 2–8%. Our proxies can be used to easily calculate a biologically realistic EEG signal directly from point-neuron simulations thus facilitating a quantitative comparison between computational models and experimental EEG recordings. Author summary: Networks of point neurons are widely used to model neural dynamics. Their output, however, cannot be directly compared to the electroencephalogram (EEG), which is one of the most used tools to non-invasively measure brain activity. To allow a direct integration between neural network theory and empirical EEG data, here we derived a new mathematical expression, termed EEG proxy, which estimates with high accuracy the EEG based simply on the variables available from simulations of point-neuron network models. To compare and validate these EEG proxies, we computed a realistic ground-truth EEG produced by a network of simulated neurons with realistic 3D morphologies that receive the same synaptic input of the simpler network of point neurons. The new obtained EEG proxies outperformed previous approaches and worked well under a wide range of network configurations with different cell morphologies, distribution of presynaptic inputs, position of the recording electrode and spatial extension of the network. The new proxies approximated well both EEG spectra and EEG evoked potentials. Our work provides important mathematical tools that allow a better interpretation of experimentally measured EEGs in terms of neural models of brain function. [ABSTRACT FROM AUTHOR]

The voltage-dependent anion channel (VDAC) is a critical β-barrel membrane protein of the mitochondrial outer membrane, which regulates the transport of ions and ATP between mitochondria and the cytoplasm. In addition, VDAC plays a central role in the control of apoptosis and is therefore of great interest in both cancer and neurodegenerative diseases. Although not fully understood, it is presumed that the gating mechanism of VDAC is governed by its N-terminal region which, in the open state of the channel, exhibits an α-helical structure positioned midway inside the pore and strongly interacting with the β-barrel wall. In the present work, we performed molecular simulations with a recently developed force field for disordered systems to shed new light on known experimental results, showing that the N-terminus of VDAC is an intrinsically disordered region (IDR). First, simulation of the N-terminal segment as a free peptide highlighted its disordered nature and the importance of using an IDR-specific force field to properly sample its conformational landscape. Secondly, accelerated dynamics simulation of a double cysteine VDAC mutant under applied voltage revealed metastable low conducting states of the channel representative of closed states observed experimentally. Related structures were characterized by partial unfolding and rearrangement of the N-terminal tail, that led to steric hindrance of the pore. Our results indicate that the disordered properties of the N-terminus are crucial to properly account for the gating mechanism of VDAC. Author summary: The voltage-dependent anion channel (VDAC) is a membrane protein playing a pivotal role in the transport of ions or ATP across the mitochondrial outer membrane as well as in the induction of apoptosis. At high enough membrane potential, VDAC is known to transition from an open state to multiple closed states, reducing the flow of ions through the channel and blocking the passage of large metabolites. While the structure of the open state was resolved more than a decade ago, a molecular description of the gating mechanism of the channel is still missing. Here we show that the N-terminus of VDAC is an intrinsically disordered region and that such a property has a profound impact on its dynamics either as a free peptide or as part of the channel. By taking disordered properties of the N-terminus into account, we managed to generate long-lived closed conformations of the channel at experimental values of the membrane potential. Our results provide new insights into the molecular mechanism driving the gating of VDAC. [ABSTRACT FROM AUTHOR]

A synaptic theory of Working Memory (WM) has been developed in the last decade as a possible alternative to the persistent spiking paradigm. In this context, we have developed a neural mass model able to reproduce exactly the dynamics of heterogeneous spiking neural networks encompassing realistic cellular mechanisms for short-term synaptic plasticity. This population model reproduces the macroscopic dynamics of the network in terms of the firing rate and the mean membrane potential. The latter quantity allows us to gain insight of the Local Field Potential and electroencephalographic signals measured during WM tasks to characterize the brain activity. More specifically synaptic facilitation and depression integrate each other to efficiently mimic WM operations via either synaptic reactivation or persistent activity. Memory access and loading are related to stimulus-locked transient oscillations followed by a steady-state activity in the β-γ band, thus resembling what is observed in the cortex during vibrotactile stimuli in humans and object recognition in monkeys. Memory juggling and competition emerge already by loading only two items. However more items can be stored in WM by considering neural architectures composed of multiple excitatory populations and a common inhibitory pool. Memory capacity depends strongly on the presentation rate of the items and it maximizes for an optimal frequency range. In particular we provide an analytic expression for the maximal memory capacity. Furthermore, the mean membrane potential turns out to be a suitable proxy to measure the memory load, analogously to event driven potentials in experiments on humans. Finally we show that the γ power increases with the number of loaded items, as reported in many experiments, while θ and β power reveal non monotonic behaviours. In particular, β and γ rhythms are crucially sustained by the inhibitory activity, while the θ rhythm is controlled by excitatory synapses. Author summary: Working Memory (WM) is the ability to temporarily store and manipulate stimuli representations that are no longer available to the senses. We have developed an innovative coarse-grained population model able to mimic several operations associated to WM. The novelty of the model consists in reproducing exactly the dynamics of spiking neural networks with realistic synaptic plasticity composed of hundreds of thousands of neurons in terms of a few macroscopic variables. These variables give access to experimentally measurable quantities such as local field potentials and electroencephalographic signals. Memory operations are joined to sustained or transient oscillations emerging in different frequency bands, in accordance with experimental results for primate and humans performing WM tasks. We have designed an architecture composed of many excitatory populations and a common inhibitory pool able to store and retain several memory items. The capacity of our multi-item architecture is around 3–5 items, a value similar to the WM capacities measured in many experiments. Furthermore, the maximal capacity is achievable only for presentation rates within an optimal frequency range. Finally, we have defined a measure of the memory load analogous to the event-related potentials employed to test humans' WM capacity during visual memory tasks. [ABSTRACT FROM AUTHOR]

Cerebellar stellate cells (CSCs) are spontaneously active, tonically firing (5-30 Hz), inhibitory interneurons that synapse onto Purkinje cells. We previously analyzed the excitability properties of CSCs, focusing on four key features: type I excitability, non-monotonic first-spike latency, switching in responsiveness and runup (i.e., temporal increase in excitability during whole-cell configuration). In this study, we extend this analysis by using whole-cell configuration to show that these neurons can also burst when treated with certain pharmacological agents separately or jointly. Indeed, treatment with (4-Aminopyridine) 4-AP, a partial blocker of delayed rectifier and A-type K+ channels, at low doses induces a bursting profile in CSCs significantly different than that produced at high doses or when it is applied at low doses but with cadmium (Cd2+), a blocker of high voltage-activated (HVA) Ca2+ channels. By expanding a previously revised Hodgkin–Huxley type model, through the inclusion of Ca2+-activated K+ (K(Ca)) and HVA currents, we explain how these bursts are generated and what their underlying dynamics are. Specifically, we demonstrate that the expanded model preserves the four excitability features of CSCs, as well as captures their bursting patterns induced by 4-AP and Cd2+. Model investigation reveals that 4-AP is potentiating HVA, inducing square-wave bursting at low doses and pseudo-plateau bursting at high doses, whereas Cd2+ is potentiating K(Ca), inducing pseudo-plateau bursting when applied in combination with low doses of 4-AP. Using bifurcation analysis, we show that spike adding in square-wave bursts is non-sequential when gradually changing HVA and K(Ca) maximum conductances, delayed Hopf is responsible for generating the plateau segment within the active phase of pseudo-plateau bursts, and bursting can become "chaotic" when HVA and K(Ca) maximum conductances are made low and high, respectively. These results highlight the secondary effects of the drugs applied and suggest that CSCs have all the ingredients needed for bursting. Author summary: Excitable cells, including neurons, fire action potentials (APs) in their membrane voltage that allow them to communicate with each other and to serve certain physiological purposes. They do so either tonically by firing APs periodically, or episodically by repeatedly firing clusters of APs (called bursts) separated by quiescent periods. Each one of those firing patterns can be neuron-specific and dependent on synaptic inputs and/or their physiological environment. Cerebellar stellate cells (CSCs) that synapse onto Purkinje cells, the sole output of the cerebellum responsible for motor control, are spontaneously active inhibitory interneurons that fire APs tonically. We previously studied the excitability properties of these neurons and showed that they possess several important key features, including type I excitability, runup, non-monotonic first spike latency and switching in responsiveness. In this study, we show that CSCs can also exhibit two modes of burst firing, called square-wave and pseudo-plateau, when treated with certain pharmacological agents. Using bifurcation theory, we demonstrate that spike adding in the square-wave burst is non-sequential, changing by several spikes when certain conductances are altered gradually. This study thus sheds lights onto the overall effects of the pharmacological agents and highlights the ability of CSCs to burst in certain biological conditions. [ABSTRACT FROM AUTHOR]

An important problem in computational neuroscience is to understand how networks of spiking neurons can carry out various computations underlying behavior. Balanced spiking networks (BSNs) provide a powerful framework for implementing arbitrary linear dynamical systems in networks of integrate-and-fire neurons. However, the classic BSN model requires near-instantaneous transmission of spikes between neurons, which is biologically implausible. Introducing realistic synaptic delays leads to an pathological regime known as "ping-ponging", in which different populations spike maximally in alternating time bins, causing network output to overshoot the target solution. Here we document this phenomenon and provide a novel solution: we show that a network can have realistic synaptic delays while maintaining accuracy and stability if neurons are endowed with conditionally Poisson firing. Formally, we propose two alternate formulations of Poisson balanced spiking networks: (1) a "local" framework, which replaces the hard integrate-and-fire spiking rule within each neuron by a "soft" threshold function, such that firing probability grows as a smooth nonlinear function of membrane potential; and (2) a "population" framework, which reformulates the BSN objective function in terms of expected spike counts over the entire population. We show that both approaches offer improved robustness, allowing for accurate implementation of network dynamics with realistic synaptic delays between neurons. Both Poisson frameworks preserve the coding accuracy and robustness to neuron loss of the original model and, moreover, produce positive correlations between similarly tuned neurons, a feature of real neural populations that is not found in the deterministic BSN. This work unifies balanced spiking networks with Poisson generalized linear models and suggests several promising avenues for future research. Author summary: A central idea in neuroscience is that populations of neurons work together to efficiently perform computations, although just how they do that remains unclear. Boerlin et al (2013) proposed a powerful framework for embedding linear dynamical systems into populations of spiking neurons, which they called balanced spiking networks (BSNs). Their approach starts from the principle that neurons greedily fire spikes to reduce error in the network output. Here we focus on a key limitation of this framework, which is that the network may become unbalanced in the presence of physiologically plausible communication delays. To overcome this shortcoming, propose two different extensions of the BSN framework that rely on probabilistic spiking. In our first model, we replace deterministic spiking of the original BSN with a Poisson spiking rule. In the second, we re-formulate the BSN objective so that Poisson spiking emerges as a way to reduce the expected network error. Our work brings the BSN framework closer to biological realism by increasing the stability and, most importantly, allowing communication delays between neurons without sacrificing accuracy. Furthermore, both probabilistic approaches reproduce key experimentally observed spiking behaviors of neural populations. [ABSTRACT FROM AUTHOR]

A biophysically detailed description of the mechanisms of the primary vision is still being developed. We have incorporated a simplified, filter-based description of retino-thalamic visual signal processing into the detailed, conductance-based refractory density description of the neuronal population activity of the primary visual cortex. We compared four mechanisms of the direction selectivity (DS), three of them being based on asymmetrical projections of different types of thalamic neurons to the cortex, distinguishing between (i) lagged and nonlagged, (ii) transient and sustained, and (iii) On and Off neurons. The fourth mechanism implies a lack of subcortical bias and is an epiphenomenon of intracortical interactions between orientation columns. The simulations of the cortical response to moving gratings have verified that first three mechanisms provide DS to an extent compared with experimental data and that the biophysical model realistically reproduces characteristics of the visual cortex activity, such as membrane potential, firing rate, and synaptic conductances. The proposed model reveals the difference between the mechanisms of both the intact and the silenced cortex, favoring the second mechanism. In the fourth case, DS is weaker but significant; it completely vanishes in the silenced cortex.DS in the On-Off mechanism derives from the nonlinear interactions within the orientation map. Results of simulations can help to identify a prevailing mechanism of DS in V1. This is a step towards a comprehensive biophysical modeling of the primary visual system in the frameworks of the population rate coding concept. Author summary: A major mechanism that underlies tuning of cortical neurons to the direction of a moving stimulus is still debated. Considering the visual cortex structured with orientation-selective columns, we have realized and compared in our biophysically detailed mathematical model four hypothetical mechanisms of the direction selectivity (DS) known from experiments. The present model accomplishes our previous model that was tuned to experimental data on excitability in slices and reproduces orientation tuning effects in vivo. In simulations, we have found that the convergence of inputs from so-called transient and sustained (or lagged and nonlagged) thalamic neurons in the cortex provides an initial bias for DS, whereas cortical interactions amplify the tuning. In the absence of any bias, DS emerges as an epiphenomenon of the orientation map. In the case of a biased convergence of On- and Off- thalamic inputs, DS emerges with the help of the intracortical interactions on the orientation map, also. Thus, we have proposed a comprehensive description of the primary vision and revealed characteristic features of different mechanisms of DS in the visual cortex with columnar structure. [ABSTRACT FROM AUTHOR]

Biological excitable media, such as cardiac or neural cells and tissue, exhibit memory in which a change in the present excitation may affect the behaviors in the next excitation. For example, a change in calcium (Ca2+) concentration in a cell in the present excitation may affect the Ca2+ dynamics in the next excitation via bi-directional coupling between voltage and Ca2+, forming a delayed feedback loop. Since the Ca2+ dynamics inside the excitable cells are spatiotemporal while the membrane voltage is a global signal, the feedback loop is then a delayed global feedback (DGF) loop. In this study, we investigate the roles of DGF in the genesis and stability of spatiotemporal excitation patterns in periodically-paced excitable media using mathematical models with different levels of complexity: a model composed of coupled FitzHugh-Nagumo units, a 3-dimensional physiologically-detailed ventricular myocyte model, and a coupled map lattice model. We investigate the dynamics of excitation patterns that are temporal period-2 (P2) and spatially concordant or discordant, such as subcellular concordant or discordant Ca2+alternans in cardiac myocytes or spatially concordant or discordant Ca2+ and repolarization alternans in cardiac tissue. Our modeling approach allows both computer simulations and rigorous analytical treatments, which lead to the following results and conclusions. When DGF is absent, concordant and discordant P2 patterns occur depending on initial conditions with the discordant P2 patterns being spatially random. When the DGF is negative, only concordant P2 patterns exist. When the DGF is positive, both concordant and discordant P2 patterns can occur. The discordant P2 patterns are still spatially random, but they satisfy that the global signal exhibits a temporal period-1 behavior. The theoretical analyses of the coupled map lattice model reveal the underlying instabilities and bifurcations caused by the DGF for the genesis, selection, and stability of spatiotemporal excitation patterns. Author summary: Understanding the mechanisms of pattern formation in biological systems is of great importance. Here we investigate the dynamical mechanisms by which delayed global feedback affects excitation pattern formation and stability in periodically-paced biological excitable media, such as cardiac or neural cells and tissue. We focus on the formation and stability of the temporal period-2 and spatially in-phase and out-of-phase excitation patterns. Using models of different levels of complexity, we show that when the delayed global feedback is negative, only the spatially in-phase patterns are stable. When the feedback is positive, both spatially in-phase and out-of-phase patterns are stable, and the out-of-phase patterns are spatially random but satisfy that the global signals are temporal period-1 solutions. [ABSTRACT FROM AUTHOR]

A basic—yet nontrivial—function which neocortical circuitry must satisfy is the ability to maintain stable spiking activity over time. Stable neocortical activity is asynchronous, critical, and low rate, and these features of spiking dynamics contribute to efficient computation and optimal information propagation. However, it remains unclear how neocortex maintains this asynchronous spiking regime. Here we algorithmically construct spiking neural network models, each composed of 5000 neurons. Network construction synthesized topological statistics from neocortex with a set of objective functions identifying naturalistic low-rate, asynchronous, and critical activity. We find that simulations run on the same topology exhibit sustained asynchronous activity under certain sets of initial membrane voltages but truncated activity under others. Synchrony, rate, and criticality do not provide a full explanation of this dichotomy. Consequently, in order to achieve mechanistic understanding of sustained asynchronous activity, we summarized activity as functional graphs where edges between units are defined by pairwise spike dependencies. We then analyzed the intersection between functional edges and synaptic connectivity- i.e. recruitment networks. Higher-order patterns, such as triplet or triangle motifs, have been tied to cooperativity and integration. We find, over time in each sustained simulation, low-variance periodic transitions between isomorphic triangle motifs in the recruitment networks. We quantify the phenomenon as a Markov process and discover that if the network fails to engage this stereotyped regime of motif dominance "cycling", spiking activity truncates early. Cycling of motif dominance generalized across manipulations of synaptic weights and topologies, demonstrating the robustness of this regime for maintenance of network activity. Our results point to the crucial role of excitatory higher-order patterns in sustaining asynchronous activity in sparse recurrent networks. They also provide a possible explanation why such connectivity and activity patterns have been prominently reported in neocortex. Author summary: Neocortical spiking activity tends to be low-rate, non-rhythmic, and operates near the critical point of a phase transition. It remains unclear how this kind of spiking activity can be maintained within a neuronal network. Neurons are leaky and individual synaptic connections are sparse and weak, making the maintenance of an asynchronous regime a nontrivial problem. Higher order patterns involving more than two units abound in neocortex, and several lines of evidence suggest that they may be instrumental for brain function. For example, stable activity in vivo displays elevated clustering dominated by specific three-node (triplet) motifs. In this study we demonstrate a link between the maintenance of asynchronous activity and triplet motifs. We algorithmically build spiking neural network models to mimic the topology of neocortex and the spiking statistics that characterize wakefulness. We show that higher order coordination of synapses is always present during sustained asynchronous activity. Coordination takes the form of transitions in time between specific triangle motifs. These motifs summarize the way spikes traverse the underlying synaptic topology. The results of our model are consistent with numerous experimental observations, and their generalizability to other weakly and sparsely connected networks is predicted. [ABSTRACT FROM AUTHOR]

Calcium oscillations and waves induce depolarization in cardiac cells which are believed to cause life-threathening arrhythimas. In this work, we study the conditions for the appearance of calcium oscillations in both a detailed subcellular model of calcium dynamics and a minimal model that takes into account just the minimal ingredients of the calcium toolkit. To avoid the effects of homeostatic changes and the interaction with the action potential we consider the somewhat artificial condition of a cell without pacing and with no calcium exchange with the extracellular medium. Both the full subcellular model and the minimal model present the same scenarios depending on the calcium load: two stationary states, one with closed ryanodine receptors (RyR) and most calcium in the cell stored in the sarcoplasmic reticulum (SR), and another, with open RyRs and a depleted SR. In between, calcium oscillations may appear. The robustness of these oscillations is determined by the amount of calsequestrin (CSQ). The lack of this buffer in the SR enhances the appearance of oscillations. The minimal model allows us to relate the stability of the oscillating state to the nullcline structure of the system, and find that its range of existence is bounded by a homoclinic and a Hopf bifurcation, resulting in a sudden transition to the oscillatory regime as the cell calcium load is increased. Adding a small amount of noise to the RyR behavior increases the parameter region where oscillations appear and provides a gradual transition from the resting state to the oscillatory regime, as observed in the subcellular model and experimentally. Author summary: In cardiac cells, calcium plays a very important role. An increase in calcium levels is the trigger used by the cell to initiate contraction. Besides, calcium modulates several transmembrane currents, affecting the cell transmembrane potential. Thus, dysregulations in calcium handling have been associated with the appearance of arrhythmias. Often, this dysregulation results in the appearance of periodic calcium waves or global oscillations, providing a pro-arrhythmic substrate. In this paper, we study the onset of calcium oscillations in cardiac cells using both a detailed subcellular model of calcium dynamics and a minimal model that takes into account the essential ingredients of the calcium toolkit. Both reproduce the main experimental results and link this behavior with the presence of different steady-state solutions and bifurcations that depend on the total amount of calcium in the cell and in the level of buffering present. We expect that this work will help to clarify the conditions under which calcium oscillations appear in cardiac myocytes and, therefore, will represent a step further in the understanding of the origin of cardiac arrhythmias. [ABSTRACT FROM AUTHOR]

The estimation of parameters controlling the electrical properties of biological neurons is essential to determine their complement of ion channels and to understand the function of biological circuits. By synchronizing conductance models to time series observations of the membrane voltage, one may construct models capable of predicting neuronal dynamics. However, identifying the actual set of parameters of biological ion channels remains a formidable theoretical challenge. Here, we present a regularization method that improves convergence towards this optimal solution when data are noisy and the model is unknown. Our method relies on the existence of an offset in parameter space arising from the interplay between model nonlinearity and experimental error. By tuning this offset, we induce saddle-node bifurcations from sub-optimal to optimal solutions. This regularization method increases the probability of finding the optimal set of parameters from 67% to 94.3%. We also reduce parameter correlations by implementing adaptive sampling and stimulation protocols compatible with parameter identifiability requirements. Our results show that the optimal model parameters may be inferred from imperfect observations provided the conditions of observability and identifiability are fulfilled. Author summary: The accurate estimation of neuronal parameters inaccessible to experiment is essential to our understanding of intracellular dynamics and to predicting the behaviour of biocircuits. However, this program is met with challenges including our lack of knowledge of the precise equations of biological neurons, their highly nonlinear response to stimulation and error introduced by the measurement apparatus. The imprecise knowledge of model and data introduces uncertainty in the parameter field. Our work describes a regularization method that arrives at the optimal parameter solution with a probability of 94%. The uncertainty on parameter estimates is further reduced with the help of an adaptive sampling method that maximises the duration of the assimilation window while keeping the size of the problem constant. Our work shows that the true configuration of a neuronal system may be inferred from time series observations provided external stimuli are calibrated to drive the system over its entire dynamic range. [ABSTRACT FROM AUTHOR]

Several decades of research suggest that weak electric fields may influence neural processing, including those induced by neuronal activity and proposed as a substrate for a potential new cellular communication system, i.e., ephaptic transmission. Here we aim to model mesoscopic ephaptic activity in the human brain and explore its trajectory during aging by characterizing the electric field generated by cortical dipoles using realistic finite element modeling. Extrapolating from electrophysiological measurements, we first observe that modeled endogenous field magnitudes are comparable to those in measurements of weak but functionally relevant self-generated fields and to those produced by noninvasive transcranial brain stimulation, and therefore possibly able to modulate neuronal activity. Then, to evaluate the role of these fields in the human cortex in large MRI databases, we adapt an interaction approximation that considers the relative orientation of neuron and field to estimate the membrane potential perturbation in pyramidal cells. We use this approximation to define a simplified metric (EMOD1) that weights dipole coupling as a function of distance and relative orientation between emitter and receiver and evaluate it in a sample of 401 realistic human brain models from healthy subjects aged 16–83. Results reveal that ephaptic coupling, in the simplified mesoscopic modeling approach used here, significantly decreases with age, with higher involvement of sensorimotor regions and medial brain structures. This study suggests that by providing the means for fast and direct interaction between neurons, ephaptic modulation may contribute to the complexity of human function for cognition and behavior, and its modification across the lifespan and in response to pathology. Author summary: We study the potential role of a less understood type of communication between neurons in the brain. While the principal mechanism for neuron communication is synaptic, active neurons generate electric fields. Whether this is physiologically relevant at the systems level or merely an epiphenomenon is uncertain because these fields are rather weak at the mesoscopic scale (i.e., between single neuron and entire brain scales). We first review converging evidence from in-vitro and in-vivo studies that suggest the former, and then, using realistic finite element modeling, we show that the electric fields generated by transcranial electrical current stimulation are of the same magnitude as endogenous ones in this scale. We then develop a method to estimate the amount of potential ephaptic or electric-field interaction in an individual brain using finite element modeling and show its decrease with age in a large cohort of 401 subjects. Ephaptic interaction may be important for complex processing in biological neural networks, because it travels at very fast speeds and provides a potential communication link across distant neurons in the cortex. Assessing the physiological relevance of this mechanism may be key in understanding some brain disorders and to design improved tES protocols. [ABSTRACT FROM AUTHOR]

In most neuronal models, ion concentrations are assumed to be constant, and effects of concentration variations on ionic reversal potentials, or of ionic diffusion on electrical potentials are not accounted for. Here, we present the electrodiffusive Pinsky-Rinzel (edPR) model, which we believe is the first multicompartmental neuron model that accounts for electrodiffusive ion concentration dynamics in a way that ensures a biophysically consistent relationship between ion concentrations, electrical charge, and electrical potentials in both the intra- and extracellular space. The edPR model is an expanded version of the two-compartment Pinsky-Rinzel (PR) model of a hippocampal CA3 neuron. Unlike the PR model, the edPR model includes homeostatic mechanisms and ion-specific leakage currents, and keeps track of all ion concentrations (Na+, K+, Ca2+, and Cl−), electrical potentials, and electrical conductivities in the intra- and extracellular space. The edPR model reproduces the membrane potential dynamics of the PR model for moderate firing activity. For higher activity levels, or when homeostatic mechanisms are impaired, the homeostatic mechanisms fail in maintaining ion concentrations close to baseline, and the edPR model diverges from the PR model as it accounts for effects of concentration changes on neuronal firing. We envision that the edPR model will be useful for the field in three main ways. Firstly, as it relaxes commonly made modeling assumptions, the edPR model can be used to test the validity of these assumptions under various firing conditions, as we show here for a few selected cases. Secondly, the edPR model should supplement the PR model when simulating scenarios where ion concentrations are expected to vary over time. Thirdly, being applicable to conditions with failed homeostasis, the edPR model opens up for simulating a range of pathological conditions, such as spreading depression or epilepsy. Author summary: Neurons generate their electrical signals by letting ions pass through their membranes. Despite this fact, most models of neurons apply the simplifying assumption that ion concentrations remain effectively constant during neural activity. This assumption is often quite good, as neurons contain a set of homeostatic mechanisms that make sure that ion concentrations vary quite little under normal circumstances. However, under some conditions, these mechanisms can fail, and ion concentrations can vary quite dramatically. Standard models are thus not able to simulate such conditions. Here, we present what to our knowledge is the first multicompartmental neuron model that accounts for ion concentration variations in a way that ensures complete and consistent ion concentration and charge conservation. In this work, we use the model to explore under which activity conditions the ion concentration variations become important for predicting the neurodynamics. We expect the model to be of great value for the field of neuroscience, as it can be used to simulate a range of pathological conditions, such as spreading depression or epilepsy, which are associated with large changes in extracellular ion concentrations. [ABSTRACT FROM AUTHOR]

This work reports a dynamical Markov state model of CLC-2 "fast" (pore) gating, based on 600 microseconds of molecular dynamics (MD) simulation. In the starting conformation of our CLC-2 model, both outer and inner channel gates are closed. The first conformational change in our dataset involves rotation of the inner-gate backbone along residues S168-G169-I170. This change is strikingly similar to that observed in the cryo-EM structure of the bovine CLC-K channel, though the volume of the intracellular (inner) region of the ion conduction pathway is further expanded in our model. From this state (inner gate open and outer gate closed), two additional states are observed, each involving a unique rotameric flip of the outer-gate residue GLUex. Both additional states involve conformational changes that orient GLUex away from the extracellular (outer) region of the ion conduction pathway. In the first additional state, the rotameric flip of GLUex results in an open, or near-open, channel pore. The equilibrium population of this state is low (∼1%), consistent with the low open probability of CLC-2 observed experimentally in the absence of a membrane potential stimulus (0 mV). In the second additional state, GLUex rotates to occlude the channel pore. This state, which has a low equilibrium population (∼1%), is only accessible when GLUex is protonated. Together, these pathways model the opening of both an inner and outer gate within the CLC-2 selectivity filter, as a function of GLUex protonation. Collectively, our findings are consistent with published experimental analyses of CLC-2 gating and provide a high-resolution structural model to guide future investigations. Author summary: In contrast to sodium-, potassium-, and calcium-selective ion channels, the roles of chloride-selective ion channels in mammalian physiology have been studied much less and are not sufficiently understood, despite known associations of chloride-channel defects with a variety of pathological conditions. CLC-2 is a voltage-activated chloride channel (one of 9 human CLC homologs) with broad tissue and organ distribution. In this work, we use simulations to model the conformational dynamics of the CLC-2 chloride ion channel selectivity filter (SF), which is the part of the protein that controls whether the channel is in an ion-conducting or non-conducting state. Our analysis identifies four primary conformational states and a specific progression through these states. Our results are consistent with structural and functional data in the literature and provide a high-resolution model for guiding further studies of CLC-2. These results will inform our understanding of how CLC-2 regulates electrical activity and ion homeostasis in the many tissues where it is expressed. [ABSTRACT FROM AUTHOR]

Intracellular calcium ([Ca2+]i) is a basic and ubiquitous cellular signal controlling a wide variety of biological processes. A remarkable example is the steering of sea urchin spermatozoa towards the conspecific egg by a spatially and temporally orchestrated series of [Ca2+]i spikes. Although this process has been an experimental paradigm for reproduction and sperm chemotaxis studies, the composition and regulation of the signalling network underlying the cytosolic calcium fluctuations are hitherto not fully understood. Here, we used a differential equations model of the signalling network to assess which set of channels can explain the characteristic envelope and temporal organisation of the [Ca2+]i-spike trains. The signalling network comprises an initial membrane hyperpolarisation produced by an Upstream module triggered by the egg-released chemoattractant peptide, via receptor activation, cGMP synthesis and decay. Followed by downstream modules leading to intraflagellar pH (pHi), voltage and [Ca2+]i fluctuations. The Upstream module outputs were fitted to kinetic data on cGMP activity and early membrane potential changes measured in bulk cell populations. Two candidate modules featuring voltage-dependent Ca2+-channels link these outputs to the downstream dynamics and can independently explain the typical decaying envelope and the progressive spacing of the spikes. In the first module, [Ca2+]i-spike trains require the concerted action of a classical CaV-like channel and a potassium channel, BK (Slo1), whereas the second module relies on pHi-dependent CatSper dynamics articulated with voltage-dependent neutral sodium-proton exchanger (NHE). We analysed the dynamics of these two modules alone and in mixed scenarios. We show that the [Ca2+]i dynamics observed experimentally after sustained alkalinisation can be reproduced by a model featuring the CatSper and NHE module but not by those including the pH-independent CaV and BK module or proportionate mixed scenarios. We conclude in favour of the module containing CatSper and NHE and highlight experimentally testable predictions that would corroborate this conclusion. Author summary: Fertilisation in marine invertebrates, such as the sea urchin, occurs during broadcast spawning events in which males and females of co-localised species ejaculate sperm and spawn eggs synchronously. During these events, spermatozoa have to find and fertilise conspecific eggs in the midst of all the other ones, which is a remarkable navigation and mating choice achievement. Sperm cells do this by navigating towards the source of species-specific peptides released by the egg, steered by spatial and temporally orchestrated peaks in intracellular calcium concentration that trigger sudden reorientations. How these calcium spikes are regulated and timed remains elusive. Different calcium channels have been implicated by indirect experimental evidence giving rise to a complex network of putative interacting components. We gained insight into the structure and function of this network by modelling it as a set of candidate modules that could be studied separately. By using this 'divide and conquer' approach to the complexity of the network, we could characterise the potential dynamics of each module and confront these dynamics with specific quantitative data. Our results indicate that the channel mediating calcium signals in sea urchin sperm is likely CatSper, a calcium channel necessary for human male fertility. [ABSTRACT FROM AUTHOR]

Message passing between components of a distributed physical system is non-instantaneous and contributes to determine the time scales of the emerging collective dynamics. In biological neuron networks this is due in part to local synaptic filtering of exchanged spikes, and in part to the distribution of the axonal transmission delays. How differently these two kinds of communication protocols affect the network dynamics is still an open issue due to the difficulties in dealing with the non-Markovian nature of synaptic transmission. Here, we develop a mean-field dimensional reduction yielding to an effective Markovian dynamics of the population density of the neuronal membrane potential, valid under the hypothesis of small fluctuations of the synaptic current. Within this limit, the resulting theory allows us to prove the formal equivalence between the two transmission mechanisms, holding for any synaptic time scale, integrate-and-fire neuron model, spike emission regimes and for different network states even when the neuron number is finite. The equivalence holds even for larger fluctuations of the synaptic input, if white noise currents are incorporated to model other possible biological features such as ionic channel stochasticity. [ABSTRACT FROM AUTHOR]

Dopaminergic neurons (DAs) of the rodent substantia nigra pars compacta (SNc) display varied electrophysiological properties in vitro. Despite this, projection patterns and functional inputs from DAs to other structures are conserved, so in vivo delivery of consistent, well-timed dopamine modulation to downstream circuits must be coordinated. Here we show robust coordination by linear parameter controllers, discovered through powerful mathematical analyses of data and models, and from which consistent control of DA subthreshold oscillations (STOs) and spontaneous firing emerges. These units of control represent coordinated intracellular variables, sufficient to regulate complex cellular properties with radical simplicity. Using an evolutionary algorithm and dimensionality reduction, we discovered metaparameters, which when regressed against STO features, revealed a 2-dimensional control plane for the neuron’s 22-dimensional parameter space that fully maps the natural range of DA subthreshold electrophysiology. This plane provided a basis for spiking currents to reproduce a large range of the naturally occurring spontaneous firing characteristics of SNc DAs. From it we easily produced a unique population of models, derived using unbiased parameter search, that show good generalization to channel blockade and compensatory intracellular mechanisms. From this population of models, we then discovered low-dimensional controllers for regulating spontaneous firing properties, and gain insight into how currents active in different voltage regimes interact to produce the emergent activity of SNc DAs. Our methods therefore reveal simple regulators of neuronal function lurking in the complexity of combined ion channel dynamics. [ABSTRACT FROM AUTHOR]

In many brain areas, sensory responses are heavily modulated by factors including attentional state, context, reward history, motor preparation, learned associations, and other cognitive variables. Modelling the effect of these modulatory factors on sensory responses has proven challenging, mostly due to the time-varying and nonlinear nature of the underlying computations. Here we present a computational model capable of capturing and dissociating multiple time-varying modulatory effects on neuronal responses on the order of milliseconds. The model’s performance is tested on extrastriate perisaccadic visual responses in nonhuman primates. Visual neurons respond to stimuli presented around the time of saccades differently than during fixation. These perisaccadic changes include sensitivity to the stimuli presented at locations outside the neuron’s receptive field, which suggests a contribution of multiple sources to perisaccadic response generation. Current computational approaches cannot quantitatively characterize the contribution of each modulatory source in response generation, mainly due to the very short timescale on which the saccade takes place. In this study, we use a high spatiotemporal resolution experimental paradigm along with a novel extension of the generalized linear model framework (GLM), termed the sparse-variable GLM, to allow for time-varying model parameters representing the temporal evolution of the system with a resolution on the order of milliseconds. We used this model framework to precisely map the temporal evolution of the spatiotemporal receptive field of visual neurons in the middle temporal area during the execution of a saccade. Moreover, an extended model based on a factorization of the sparse-variable GLM allowed us to disassociate and quantify the contribution of individual sources to the perisaccadic response. Our results show that our novel framework can precisely capture the changes in sensitivity of neurons around the time of saccades, and provide a general framework to quantitatively track the role of multiple modulatory sources over time. [ABSTRACT FROM AUTHOR]

The mechanism(s) of action of most commonly used pharmacological blockers of voltage-gated ion channels are well understood; however, this knowledge is rarely considered when interpreting experimental data. Effects of blockade are often assumed to be equivalent, regardless of the mechanism of the blocker involved. Using computer simulations, we demonstrate that this assumption may not always be correct. We simulate the blockade of a persistent sodium current (INaP), proposed to underlie rhythm generation in pre-Bötzinger complex (pre-BötC) respiratory neurons, via two distinct pharmacological mechanisms: (1) pore obstruction mediated by tetrodotoxin and (2) altered inactivation dynamics mediated by riluzole. The reported effects of experimental application of tetrodotoxin and riluzole in respiratory circuits are diverse and seemingly contradictory and have led to considerable debate within the field as to the specific role of INaP in respiratory circuits. The results of our simulations match a wide array of experimental data spanning from the level of isolated pre-BötC neurons to the level of the intact respiratory network and also generate a series of experimentally testable predictions. Specifically, in this study we: (1) provide a mechanistic explanation for seemingly contradictory experimental results from in vitro studies of INaP block, (2) show that the effects of INaP block in in vitro preparations are not necessarily equivalent to those in more intact preparations, (3) demonstrate and explain why riluzole application may fail to effectively block INaP in the intact respiratory network, and (4) derive the prediction that effective block of INaP by low concentration tetrodotoxin will stop respiratory rhythm generation in the intact respiratory network. These simulations support a critical role for INaP in respiratory rhythmogenesis in vivo and illustrate the importance of considering mechanism when interpreting and simulating data relating to pharmacological blockade. [ABSTRACT FROM AUTHOR]

Membrane-protein design is an exciting and increasingly successful research area which has led to landmarks including the design of stable and accurate membrane-integral proteins based on coiled-coil motifs. Design of topologically more complex proteins, such as most receptors, channels, and transporters, however, demands an energy function that balances contributions from intra-protein contacts and protein-membrane interactions. Recent advances in water-soluble all-atom energy functions have increased the accuracy in structure-prediction benchmarks. The plasma membrane, however, imposes different physical constraints on protein solvation. To understand these constraints, we recently developed a high-throughput experimental screen, called dsTβL, and inferred apparent insertion energies for each amino acid at dozens of positions across the bacterial plasma membrane. Here, we express these profiles as lipophilicity energy terms in Rosetta and demonstrate that the new energy function outperforms previous ones in modelling and design benchmarks. Rosetta ab initio simulations starting from an extended chain recapitulate two-thirds of the experimentally determined structures of membrane-spanning homo-oligomers with <2.5Å root-mean-square deviation within the top-predicted five models (available online: ). Furthermore, in two sequence-design benchmarks, the energy function improves discrimination of stabilizing point mutations and recapitulates natural membrane-protein sequences of known structure, thereby recommending this new energy function for membrane-protein modelling and design. [ABSTRACT FROM AUTHOR]

This study presents a modelling framework in which information on muscle fiber direction and orientation during contraction is derived from diffusion tensor imaging (DTI) and incorporated in a computational model of the surface electromyographic (EMG) signal. The proposed model makes use of the principle of reciprocity to simultaneously calculate the electric potentials produced at the recording electrode by charges distributed along an arbitrary number of muscle fibers within the muscle, allowing for a computationally efficient evaluation of extracellular motor unit action potentials. The approach is applied to the complex architecture of the first dorsal interosseous (FDI) muscle of the hand to simulate EMG during index finger flexion and abduction. Using diffusion tensor imaging methods, the results show how muscle fiber orientation and curvature in this intrinsic hand muscle change during flexion and abduction. Incorporation of anatomically accurate muscle architecture and other hand tissue morphologies enable the model to capture variations in extracellular action potential waveform shape across the motor unit population and to predict experimentally observed differences in EMG signal features when switching from index finger abduction to flexion. The simulation results illustrate how structural and electrical properties of the tissues comprising the volume conductor, in combination with fiber direction and curvature, shape the detected action potentials. Using the model, the relative contribution of motor units of different sizes located throughout the muscle under both conditions is examined, yielding a prediction of the detection profile of the surface EMG electrode array over the muscle cross-section. [ABSTRACT FROM AUTHOR]

Despite its critical importance in experimental and clinical neuroscience, at present there is no systematic method to predict which neural elements will be activated by a given stimulation regime. Here we develop a novel approach to model the effect of cortical stimulation on spiking probability of neurons in a volume of tissue, by applying an analytical estimate of stimulation-induced activation of different cell types across cortical layers. We utilize the morphology and properties of axonal arborization profiles obtained from publicly available anatomical reconstructions of the twelve main categories of neocortical neurons to derive the dependence of activation probability on cell type, layer and distance from the source. We then propagate this activity through the local network incorporating connectivity, synaptic and cellular properties. Our work predicts that (a) intracranial cortical stimulation induces selective activation across cell types and layers; (b) superficial anodal stimulation is more effective than cathodal at cell activation; (c) cortical surface stimulation focally activates layer I axons, and (d) there is an optimal stimulation intensity capable of eliciting cell activation lasting beyond the end of stimulation. We conclude that selective effects of cortical electrical stimulation across cell types and layers are largely driven by their different axonal arborization and myelination profiles. [ABSTRACT FROM AUTHOR]

Multiple cellular organelles tightly orchestrate intracellular calcium (Ca2+) dynamics to regulate cellular activities and maintain homeostasis. The interplay between the endoplasmic reticulum (ER), a major store of intracellular Ca2+, and mitochondria, an important source of adenosine triphosphate (ATP), has been the subject of much research, as their dysfunction has been linked with metabolic diseases. Interestingly, throughout the cell’s cytosolic domain, these two organelles share common microdomains called mitochondria-associated ER membranes (MAMs), where their membranes are in close apposition. The role of MAMs is critical for intracellular Ca2+ dynamics as they provide hubs for direct Ca2+ exchange between the organelles. A recent experimental study reported correlation between obesity and MAM formation in mouse liver cells, and obesity-related cellular changes that are closely associated with the regulation of Ca2+ dynamics. We constructed a mathematical model to study the effects of MAM Ca2+ dynamics on global Ca2+ activities. Through a series of model simulations, we investigated cellular mechanisms underlying the altered Ca2+ dynamics in the cells under obesity. We predict that, as the dosage of stimulus gradually increases, liver cells from obese mice will reach the state of saturated cytosolic Ca2+ concentration at a lower stimulus concentration, compared to cells from healthy mice. [ABSTRACT FROM AUTHOR]

Pituitary endocrine cells fire action potentials (APs) to regulate their cytosolic Ca2+ concentration and hormone secretion rate. Depending on animal species, cell type, and biological conditions, pituitary APs are generated either by TTX-sensitive Na+ currents (INa), high-voltage activated Ca2+ currents (ICa), or by a combination of the two. Previous computational models of pituitary cells have mainly been based on data from rats, where INa is largely inactivated at the resting potential, and spontaneous APs are predominantly mediated by ICa. Unlike in rats, spontaneous INa-mediated APs are consistently seen in pituitary cells of several other animal species, including several species of fish. In the current work we develop a computational model of gonadotropin releasing cells in the teleost fish medaka (Oryzias latipes). The model stands out from previous modeling efforts by being (1) the first model of a pituitary cell in teleosts, (2) the first pituitary cell model that fires sponateous APs that are predominantly mediated by INa, and (3) the first pituitary cell model where the kinetics of the depolarizing currents, INa and ICa, are directly fitted to voltage-clamp data. We explore the firing properties of the model, and compare it to the properties of previous models that fire ICa-based APs. We put a particular focus on how the big conductance K+ current (IBK) modulates the AP shape. Interestingly, we find that IBK can prolong AP duration in models that fire ICa-based APs, while it consistently shortens the duration of the predominantly INa-mediated APs in the medaka gonadotroph model. Although the model is constrained to experimental data from gonadotroph cells in medaka, it may likely provide insights also into other pituitary cell types that fire INa-mediated APs. [ABSTRACT FROM AUTHOR]

The formation of a membrane-enveloped virus starts with the assembly of a curved layer of capsid proteins lining the interior of the plasma membrane (PM) of the host cell. This layer develops into a spherical shell (capsid) enveloped by a lipid-rich membrane. In many cases, the budding process stalls prior to the release of the virus. Recently, Brownian dynamics simulations of a coarse-grained model system reproduced protracted pausing and stalling, which suggests that the origin of pausing/stalling is to be found in the physics of the budding process. Here, we propose that the pausing/stalling observed in the simulations can be understood as a purely kinetic phenomenon associated with the neck geometry. A geometrical potential energy barrier develops during the budding that must be overcome by capsid proteins diffusing along the membrane prior to incorporation into the capsid. The barrier is generated by a conflict between the positive Gauss curvature of the assembling capsid and the negative Gauss curvature of the neck region. A continuum theory description is proposed and is compared with the Brownian simulations of the budding of enveloped viruses. [ABSTRACT FROM AUTHOR]

Spontaneous sub-cellular calcium release events (SCRE) are conjectured to promote rapid arrhythmias associated with conditions such as heart failure and atrial fibrillation: they can underlie the emergence of spontaneous action potentials in single cells which can lead to arrhythmogenic triggers in tissue. The multi-scale mechanisms of the development of SCRE into arrhythmia triggers, and their dynamic interaction with the tissue substrate, remain elusive; rigorous and simultaneous study of dynamics from the nanometre to the centimetre scale is a major challenge. The aim of this study was to develop a computational approach to overcome this challenge and study potential bi-directional coupling between sub-cellular and tissue-scale arrhythmia phenomena. A framework comprising a hierarchy of computational models was developed, which includes detailed single-cell models describing spatio-temporal calcium dynamics in 3D, efficient non-spatial cell models, and both idealised and realistic tissue models. A phenomenological approach was implemented to reproduce SCRE morphology and variability in the efficient cell models, comprising the definition of analytical Spontaneous Release Functions (SRF) whose parameters may be randomly sampled from appropriate distributions in order to match either the 3D cell models or experimental data. Pro-arrhythmogenic pacing protocols were applied to initiate re-entry and promote calcium overload, leading to the emergence of SCRE. The SRF accurately reproduced the dynamics of SCRE and its dependence on environment variables under multiple different conditions. Sustained re-entrant excitation promoted calcium overload, and led to the emergence of focal excitations after termination. A purely functional mechanism of re-entry and focal activity localisation was demonstrated, related to the unexcited spiral wave core. In conclusion, a novel approach has been developed to dynamically model SCRE at the tissue scale, which facilitates novel, detailed multi-scale mechanistic analysis. It was revealed that complex re-entrant excitation patterns and SCRE may be bi-directionally coupled, promoting novel mechanisms of arrhythmia perpetuation. [ABSTRACT FROM AUTHOR]

Deep brain stimulation (DBS) is known to be an effective treatment for a variety of neurological disorders, including Parkinson’s disease and essential tremor (ET). At present, it involves administering a train of pulses with constant frequency via electrodes implanted into the brain. New ‘closed-loop’ approaches involve delivering stimulation according to the ongoing symptoms or brain activity and have the potential to provide improvements in terms of efficiency, efficacy and reduction of side effects. The success of closed-loop DBS depends on being able to devise a stimulation strategy that minimizes oscillations in neural activity associated with symptoms of motor disorders. A useful stepping stone towards this is to construct a mathematical model, which can describe how the brain oscillations should change when stimulation is applied at a particular state of the system. Our work focuses on the use of coupled oscillators to represent neurons in areas generating pathological oscillations. Using a reduced form of the Kuramoto model, we analyse how a patient should respond to stimulation when neural oscillations have a given phase and amplitude, provided a number of conditions are satisfied. For such patients, we predict that the best stimulation strategy should be phase specific but also that stimulation should have a greater effect if applied when the amplitude of brain oscillations is lower. We compare this surprising prediction with data obtained from ET patients. In light of our predictions, we also propose a new hybrid strategy which effectively combines two of the closed-loop strategies found in the literature, namely phase-locked and adaptive DBS. [ABSTRACT FROM AUTHOR]

We have previously shown that the physiological size of postsynaptic currents maximises energy efficiency rather than information transfer across the retinothalamic relay synapse. Here, we investigate information transmission and postsynaptic energy use at the next synapse along the visual pathway: from relay neurons in the thalamus to spiny stellate cells in layer 4 of the primary visual cortex (L4SS). Using both multicompartment Hodgkin-Huxley-type simulations and electrophysiological recordings in rodent brain slices, we find that increasing or decreasing the postsynaptic conductance of the set of thalamocortical inputs to one L4SS cell decreases the energy efficiency of information transmission from a single thalamocortical input. This result is obtained in the presence of random background input to the L4SS cell from excitatory and inhibitory corticocortical connections, which were simulated (both excitatory and inhibitory) or injected experimentally using dynamic-clamp (excitatory only). Thus, energy efficiency is not a unique property of strong relay synapses: even at the relatively weak thalamocortical synapse, each of which contributes minimally to the output firing of the L4SS cell, evolutionarily-selected postsynaptic properties appear to maximise the information transmitted per energy used. [ABSTRACT FROM AUTHOR]

A central issue in the study of the neural generation of respiratory rhythms is the role of the intrinsic pacemaking capabilities that some respiratory neurons exhibit. The debate on this issue has occurred in parallel to investigations of interactions among respiratory network neurons and how these contribute to respiratory behavior. In this computational study, we demonstrate how these two issues are inextricably linked. We use simulations and dynamical systems analysis to show that once a conditional respiratory pacemaker, which can be tuned across oscillatory and non-oscillatory dynamic regimes in isolation, is embedded into a respiratory network, its dynamics become masked: the network exhibits similar dynamical properties regardless of the conditional pacemaker node’s tuning, and that node’s outputs are dominated by network influences. Furthermore, the outputs of the respiratory central pattern generator as a whole are invariant to these changes of dynamical properties, which ensures flexible and robust performance over a wide dynamic range. [ABSTRACT FROM AUTHOR]

Orientation selectivity is a key property of primary visual cortex that contributes, downstream, to object recognition. The origin of orientation selectivity, however, has been debated for decades. It is known that on- and off-centre subcortical pathways converge onto single neurons in primary visual cortex, and that the spatial offset between these pathways gives rise to orientation selectivity. On- and off-centre pathways are intermingled, however, so it is unclear how their inputs to cortex come to be spatially segregated. We here describe a model in which the segregation occurs through Hebbian strengthening and weakening of geniculocortical synapses during the development of the visual system. Our findings include the following. 1. Neighbouring on- and off-inputs to cortex largely cancelled each other at the start of development. At each receptive field location, the Hebbian process increased the strength of one input sign at the expense of the other sign, producing a spatial segregation of on- and off-inputs. 2. The resulting orientation selectivity was precise in that the bandwidths of the orientation tuning functions fell within empirical estimates. 3. The model produced maps of preferred orientation–complete with iso-orientation domains and pinwheels–similar to those found in real cortex. 4. These maps did not originate in cortical processes, but from clustering of off-centre subcortical pathways and the relative location of neighbouring on-centre clusters. We conclude that a model with intermingled on- and off-pathways shaped by Hebbian synaptic plasticity can explain both the origin and development of orientation selectivity. [ABSTRACT FROM AUTHOR]

Neuroscience models come in a wide range of scales and specificity, from mean-field rate models to large-scale networks of spiking neurons. There are potential trade-offs between simplicity and realism, versatility and computational speed. This paper is about large-scale cortical network models, and the question we address is one of scalability: would scaling down cell density impact a network’s ability to reproduce cortical dynamics and function? We investigated this problem using a previously constructed realistic model of the monkey visual cortex that is true to size. Reducing cell density gradually up to 50-fold, we studied changes in model behavior. Size reduction without parameter adjustment was catastrophic. Surprisingly, relatively minor compensation in synaptic weights guided by a theoretical algorithm restored mean firing rates and basic function such as orientation selectivity to models 10-20 times smaller than the real cortex. Not all was normal in the reduced model cortices: intracellular dynamics acquired a character different from that of real neurons, and while the ability to relay feedforward inputs remained intact, reduced models showed signs of deficiency in functions that required dynamical interaction among cortical neurons. These findings are not confined to models of the visual cortex, and modelers should be aware of potential issues that accompany size reduction. Broader implications of this study include the importance of homeostatic maintenance of firing rates, and the functional consequences of feedforward versus recurrent dynamics, ideas that may shed light on other species and on systems suffering cell loss. [ABSTRACT FROM AUTHOR]

Neurons in different layers of sensory cortex generally have different functional properties. But what determines firing rates and tuning properties of neurons in different layers? Orientation selectivity in primary visual cortex (V1) is an interesting case to study these questions. Thalamic projections essentially determine the preferred orientation of neurons that receive direct input. But how is this tuning propagated though layers, and how can selective responses emerge in layers that do not have direct access to the thalamus? Here we combine numerical simulations with mathematical analyses to address this problem. We find that a large-scale network, which just accounts for experimentally measured layer and cell-type specific connection probabilities, yields firing rates and orientation selectivities matching electrophysiological recordings in rodent V1 surprisingly well. Further analysis, however, is complicated by the fact that neuronal responses emerge in a dynamic fashion and cannot be directly inferred from static neuroanatomy, as some connections tend to have unintuitive effects due to recurrent interactions and strong feedback loops. These emergent phenomena can be understood by linearizing and coarse-graining. In fact, we were able to derive a low-dimensional linear dynamical system effectively describing stimulus-driven activity layer by layer. This low-dimensional system explains layer-specific firing rates and orientation tuning by accounting for the different gain factors of the aggregate system. Our theory can also be used to design novel optogenetic stimulation experiments, thus facilitating further exploration of the interplay between connectivity and function. [ABSTRACT FROM AUTHOR]

The dimensionality of a network’s collective activity is of increasing interest in neuroscience. This is because dimensionality provides a compact measure of how coordinated network-wide activity is, in terms of the number of modes (or degrees of freedom) that it can independently explore. A low number of modes suggests a compressed low dimensional neural code and reveals interpretable dynamics [], while findings of high dimension may suggest flexible computations [, ]. Here, we address the fundamental question of how dimensionality is related to connectivity, in both autonomous and stimulus-driven networks. Working with a simple spiking network model, we derive three main findings. First, the dimensionality of global activity patterns can be strongly, and systematically, regulated by local connectivity structures. Second, the dimensionality is a better indicator than average correlations in determining how constrained neural activity is. Third, stimulus evoked neural activity interacts systematically with neural connectivity patterns, leading to network responses of either greater or lesser dimensionality than the stimulus. [ABSTRACT FROM AUTHOR]

Neurons utilize bursts of action potentials as an efficient and reliable way to encode information. It is likely that the intrinsic membrane properties of neurons involved in burst generation may also participate in preserving its temporal features. Here we examined the contribution of the persistent and resurgent components of voltage-gated Na+ currents in modulating the burst discharge in sensory neurons. Using mathematical modeling, theory and dynamic-clamp electrophysiology, we show that, distinct from the persistent Na+ component which is important for membrane resonance and burst generation, the resurgent Na+ can help stabilize burst timing features including the duration and intervals. Moreover, such a physiological role for the resurgent Na+ offered noise tolerance and preserved the regularity of burst patterns. Model analysis further predicted a negative feedback loop between the persistent and resurgent gating variables which mediate such gain in burst stability. These results highlight a novel role for the voltage-gated resurgent Na+ component in moderating the entropy of burst-encoded neural information. [ABSTRACT FROM AUTHOR]

In order to further our understanding of how gene expression contributes to key functional properties of neurons, we combined publicly accessible gene expression, electrophysiology, and morphology measurements to identify cross-cell type correlations between these data modalities. Building on our previous work using a similar approach, we distinguished between correlations which were “class-driven,” meaning those that could be explained by differences between excitatory and inhibitory cell classes, and those that reflected graded phenotypic differences within classes. Taking cell class identity into account increased the degree to which our results replicated in an independent dataset as well as their correspondence with known modes of ion channel function based on the literature. We also found a smaller set of genes whose relationships to electrophysiological or morphological properties appear to be specific to either excitatory or inhibitory cell types. Next, using data from PatchSeq experiments, allowing simultaneous single-cell characterization of gene expression and electrophysiology, we found that some of the gene-property correlations observed across cell types were further predictive of within-cell type heterogeneity. In summary, we have identified a number of relationships between gene expression, electrophysiology, and morphology that provide testable hypotheses for future studies. [ABSTRACT FROM AUTHOR]

The ventromedial nucleus of the hypothalamus (VMN) has an important role in diverse behaviours. The common involvement in these of sex steroids, nutritionally-related signals, and emotional inputs from other brain areas, suggests that, at any given time, its output is in one of a discrete number of possible states corresponding to discrete motivational drives. Here we explored how networks of VMN neurons might generate such a decision-making architecture. We began with minimalist assumptions about the intrinsic properties of VMN neurons inferred from electrophysiological recordings of these neurons in rats in vivo, using an integrate-and-fire based model modified to simulate activity-dependent post-spike changes in neuronal excitability. We used a genetic algorithm based method to fit model parameters to the statistical features of spike patterning in each cell. The spike patterns in both recorded cells and model cells were assessed by analysis of interspike interval distributions and of the index of dispersion of firing rate over different binwidths. Simpler patterned cells could be closely matched by single neuron models incorporating a hyperpolarising afterpotential and either a slow afterhyperpolarisation or a depolarising afterpotential, but many others could not. We then constructed network models with the challenge of explaining the more complex patterns. We assumed that neurons of a given type (with heterogeneity introduced by independently random patterns of external input) were mutually interconnected at random by excitatory synaptic connections (with a variable delay and a random chance of failure). Simple network models of one or two cell types were able to explain the more complex patterns. We then explored the information processing features of such networks that might be relevant for a decision-making network. We concluded that rhythm generation (in the slow theta range) and bistability arise as emergent properties of networks of heterogeneous VMN neurons. [ABSTRACT FROM AUTHOR]

The conduction of electrical signals through cardiac tissue is essential for maintaining the function of the heart, and conduction abnormalities are known to potentially lead to life-threatening arrhythmias. The properties of cardiac conduction have therefore been the topic of intense study for decades, but a number of questions related to the mechanisms of conduction still remain unresolved. In this paper, we demonstrate how the so-called EMI model may be used to study some of these open questions. In the EMI model, the extracellular space, the cell membrane, the intracellular space and the cell connections are all represented as separate parts of the computational domain, and the model therefore allows for study of local properties that are hard to represent in the classical homogenized bidomain or monodomain models commonly used to study cardiac conduction. We conclude that a non-uniform sodium channel distribution increases the conduction velocity and decreases the time delays over gap junctions of reduced coupling in the EMI model simulations. We also present a theoretical optimal cell length with respect to conduction velocity and consider the possibility of ephaptic coupling (i.e. cell-to-cell coupling through the extracellular potential) acting as an alternative or supporting mechanism to gap junction coupling. We conclude that for a non-uniform distribution of sodium channels and a sufficiently small intercellular distance, ephaptic coupling can influence the dynamics of the sodium channels and potentially provide cell-to-cell coupling when the gap junction connection is absent. [ABSTRACT FROM AUTHOR]

Several recent studies have shown that neural activity in vivo tends to be constrained to a low-dimensional manifold. Such activity does not arise in simulated neural networks with homogeneous connectivity and it has been suggested that it is indicative of some other connectivity pattern in neuronal networks. In particular, this connectivity pattern appears to be constraining learning so that only neural activity patterns falling within the intrinsic manifold can be learned and elicited. Here, we use three different models of spiking neural networks (FORCE learning, the Neural Engineering Framework and Efficient Coding) to demonstrate how the intrinsic manifold can be made a direct consequence of the circuit connectivity. Using this relationship between the circuit connectivity and the intrinsic manifold, we show that learning patterns outside the intrinsic manifold corresponds to much larger changes in synaptic weights than learning patterns within the intrinsic manifold. Assuming larger changes to synaptic weights requires extensive learning, this observation provides an explanation of why learning is easier when it does not require the neural activity to leave its intrinsic manifold. [ABSTRACT FROM AUTHOR]

Electroencephalography (EEG) provides a non-invasive measure of brain electrical activity. Neural population models, where large numbers of interacting neurons are considered collectively as a macroscopic system, have long been used to understand features in EEG signals. By tuning dozens of input parameters describing the excitatory and inhibitory neuron populations, these models can reproduce prominent features of the EEG such as the alpha-rhythm. However, the inverse problem, of directly estimating the parameters from fits to EEG data, remains unsolved. Solving this multi-parameter non-linear fitting problem will potentially provide a real-time method for characterizing average neuronal properties in human subjects. Here we perform unbiased fits of a 22-parameter neural population model to EEG data from 82 individuals, using both particle swarm optimization and Markov chain Monte Carlo sampling. We estimate how much is learned about individual parameters by computing Kullback-Leibler divergences between posterior and prior distributions for each parameter. Results indicate that only a single parameter, that determining the dynamics of inhibitory synaptic activity, is directly identifiable, while other parameters have large, though correlated, uncertainties. We show that the eigenvalues of the Fisher information matrix are roughly uniformly spaced over a log scale, indicating that the model is sloppy, like many of the regulatory network models in systems biology. These eigenvalues indicate that the system can be modeled with a low effective dimensionality, with inhibitory synaptic activity being prominent in driving system behavior. [ABSTRACT FROM AUTHOR]