Excitatory and inhibitory effects of HCN channel modulation on excitability of layer V pyramidal cells

Dendrites of cortical pyramidal cells are densely populated by hyperpolarization-activated cyclic nucleotide-gated (HCN) channels, a.k.a. Ih channels. Ih channels are targeted by multiple neuromodulatory pathways, and thus are one of the key ion-channel populations regulating the pyramidal cell activity. Previous observations and theories attribute opposing effects of the Ih channels on neuronal excitability due to their mildly hyperpolarized reversal potential. These effects are difficult to measure experimentally due to the fine spatiotemporal landscape of the Ih activity in the dendrites, but computational models provide an efficient tool for studying this question in a reduced but generalizable setting. In this work, we build upon existing biophysically detailed models of thick-tufted layer V pyramidal cells and model the effects of over- and under-expression of Ih channels as well as their neuromodulation. We show that Ih channels facilitate the action potentials of layer V pyramidal cells in response to proximal dendritic stimulus while they hinder the action potentials in response to distal dendritic stimulus at the apical dendrite. We also show that the inhibitory action of the Ih channels in layer V pyramidal cells is due to the interactions between Ih channels and a hot zone of low voltage-activated Ca2+ channels at the apical dendrite. Our simulations suggest that a combination of Ih-enhancing neuromodulation at the proximal part of the apical dendrite and Ih-inhibiting modulation at the distal part of the apical dendrite can increase the layer V pyramidal excitability more than either of the two alone. Our analyses uncover the effects of Ih-channel neuromodulation of layer V pyramidal cells at a single-cell level and shed light on how these neurons integrate information and enable higher-order functions of the brain.publishedVersionPeer reviewe

Modelling Structure and Dynamics of Complex Systems: Applications to Neuronal Networks

Complex systems theory is a mathematical framework for studying interconnected dynamical objects. Usually these objects themselves are by construction simple, and their temporal behavior in isolation is easily predictable, but the way they are interconnected into a network allows emergence of complex, non-obvious phenomena. The emergent phenomena and their stability are dependent on both the intrinsic dynamics of the objects, the types of interactions between the objects, and the connectivity patterns between the objects. This work focuses on the third aspect, i.e., the structure of the network, although the other two aspects are inherently present in the study as well. Tools from graph theory are applied to generate and analyze the network structure, and the effect of the structure on the network dynamics is analyzed by various methods. The objects of interest are biological and physical systems, and special attention is given to spiking neuronal networks, i.e., networks of nerve cells that communicate by transmitting and receiving action potentials. In this thesis, methods for modelling spiking neuronal networks are introduced. Different point neuron models, including the integrate-and-fire model, are presented and applied to study the collective behaviour of the neurons. Special focus is placed on the emergence of network bursts, i.e., short periods of network-wide high-frequency firing. The occurrence of this behaviour is stable in certain regimes of connection strengths. This work shows that the network bursting is found to be more frequent in locally connected networks than in non-local networks, such as randomly connected networks. To gain a deeper insight, the aspects of structure that promote the bursting behaviour are analyzed by graph-theoretic means. The clustering coefficient and the maximal eigenvalue of the connectivity matrix are found the most important measures of structure in this matter, both expressing their relevance under different structural conditions. A range of different network structures are applied to confirm this result. A special class of connectivity is studied in more detail, namely, the connectivity patterns produced by simulations of growing and interconnecting neurons placed on a 2-dimensional array. Two simulators of growth are applied for this purpose. In addition, a more abstract class of dynamical systems, the Boolean networks, are considered. These systems were originally introduced as a model for genetic regulatory networks, but have thereafter been extensively used for more general studies of complex systems. In this work, measures of information diversity and complexity are applied to several types of systems that obey Boolean dynamics. The random Boolean networks are shown to possess high temporal complexity prior to reaching an attractor. Similarly, high values of complexity are found at a transition stage of another dynamical system, the lattice gas automaton, which can be formulated using the Boolean network framework as well. The temporal maximization of the complexity near the transitions between different dynamical regimes could therefore be a more general phenomenon in complex networks. The applicability of the information-theoretic framework is also confirmed in a study of bursting neuronal networks, where different types of networks are shown to be separable by the intrinsic information distance distributions they produce. The connectivities of the networks studied in this thesis are analyzed using graph-theoretic tools. The graph theory provides a mathematical framework for studying the structure of complex systems and how it affects the system dynamics. In the studies of the nervous system, detailed maps on the connections between neurons have been collected, although such data are yet scarce and laborious to obtain experimentally. This work shows which aspects of the structure are relevant for the dynamics of spontaneously bursting neuronal networks. Such information could be useful in directing the experiments to measure only the relevant aspects of the structure instead of assessing the whole connectome. In addition, the framework of generating the network structure by animating the growth of the neurons, as presented in this thesis, could serve in simulations of the nervous system as a reliable alternative to importing the experimentally obtained connectome

The effect of alterations of schizophrenia-associated genes on gamma band oscillations

© The Author(s) 2022. This article is licensed under a Creative Commons Attribution 4.0 International License. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.Abnormalities in the synchronized oscillatory activity of neurons in general and, specifically in the gamma band, might play a crucial role in the pathophysiology of schizophrenia. While these changes in oscillatory activity have traditionally been linked to alterations at the synaptic level, we demonstrate here, using computational modeling, that common genetic variants of ion channels can contribute strongly to this effect. Our model of primary auditory cortex highlights multiple schizophrenia-associated genetic variants that reduce gamma power in an auditory steady-state response task. Furthermore, we show that combinations of several of these schizophrenia-associated variants can produce similar effects as the more traditionally considered synaptic changes. Overall, our study provides a mechanistic link between schizophrenia-associated common genetic variants, as identified by genome-wide association studies, and one of the most robust neurophysiological endophenotypes of schizophrenia.Peer reviewedFinal Published versio

Functional effects of schizophrenia-linked genetic variants on intrinsic single-neuron excitability: A modeling study

Background: Recent genome-wide association studies (GWAS) have identified a large number of genetic risk factors for schizophrenia (SCZ) featuring ion channels and calcium transporters. For some of these risk factors, independent prior investigations have examined the effects of genetic alterations on the cellular electrical excitability and calcium homeostasis. In the present proof-of-concept study, we harnessed these experimental results for modeling of computational properties on layer V cortical pyramidal cell and identify possible common alterations in behavior across SCZ-related genes. Methods: We applied a biophysically detailed multi-compartmental model to study the excitability of a layer V pyramidal cell. We reviewed the literature on functional genomics for variants of genes associated with SCZ, and used changes in neuron model parameters to represent the effects of these variants. Results: We present and apply a framework for examining the effects of subtle single nucleotide polymorphisms in ion channel and Ca2+ transporter-encoding genes on neuron excitability. Our analysis indicates that most of the considered SCZ- related genetic variants affect the spiking behavior and intracellular calcium dynamics resulting from summation of inputs across the dendritic tree. Conclusions: Our results suggest that alteration in the ability of a single neuron to integrate the inputs and scale its excitability may constitute a fundamental mechanistic contributor to mental disease, alongside with the previously proposed deficits in synaptic communication and network behavior

A stepwise neuron model fitting procedure designed for recordings with high spatial resolution : Application to layer 5 pyramidal cells

© 2017 The Author(s). Published by Elsevier B. V. This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution, and reproduction in any medium, provided the original work is properly cited.Background Recent progress in electrophysiological and optical methods for neuronal recordings provides vast amounts of high-resolution data. In parallel, the development of computer technology has allowed simulation of ever-larger neuronal circuits. A challenge in taking advantage of these developments is the construction of single-cell and network models in a way that faithfully reproduces neuronal biophysics with subcellular level of details while keeping the simulation costs at an acceptable level. New method In this work, we develop and apply an automated, stepwise method for fitting a neuron model to data with fine spatial resolution, such as that achievable with voltage sensitive dyes (VSDs) and Ca2+ imaging. Result We apply our method to simulated data from layer 5 pyramidal cells (L5PCs) and construct a model with reduced neuronal morphology. We connect the reduced-morphology neurons into a network and validate against simulated data from a high-resolution L5PC network model. Comparison with existing methods Our approach combines features from several previously applied model-fitting strategies. The reduced-morphology neuron model obtained using our approach reliably reproduces the membrane-potential dynamics across the dendrites as predicted by the full-morphology model. Conclusions The network models produced using our method are cost-efficient and predict that interconnected L5PCs are able to amplify delta-range oscillatory inputs across a large range of network sizes and topologies, largely due to the medium after hyperpolarization mediated by the Ca2+-activated SK current.Peer reviewedFinal Published versio

Genetic mechanisms for impaired synaptic plasticity in schizophrenia revealed by computational modeling

Schizophrenia phenotypes are suggestive of impaired cortical plasticity in the disease, but the mechanisms of these deficits are unknown. Genomic association studies have implicated a large number of genes that regulate neuromodulation and plasticity, indicating that the plasticity deficits have a genetic origin. Here, we used biochemically detailed computational modeling of postsynaptic plasticity to investigate how schizophrenia-associated genes regulate long-term potentiation (LTP) and depression (LTD). We combined our model with data from postmortem RNA expression studies (CommonMind gene-expression datasets) to assess the consequences of altered expression of plasticity-regulating genes for the amplitude of LTP and LTD. Our results show that the expression alterations observed post mortem, especially those in the anterior cingulate cortex, lead to impaired protein kinase A (PKA)-pathway-mediated LTP in synapses containing GluR1 receptors. We validated these findings using a genotyped electroencephalogram (EEG) dataset where polygenic risk scores for synaptic and ion channel-encoding genes as well as modulation of visual evoked potentials were determined for 286 healthy controls. Our results provide a possible genetic mechanism for plasticity impairments in schizophrenia, which can lead to improved understanding and, ultimately, treatment of the disorder.Peer reviewe