Control of clustered action potential firing in a mathematical model of entorhinal cortex stellate cells.

This is the author accepted manuscript. The final version is available from Elsevier via the DOI in this record.The entorhinal cortex is a crucial component of our memory and spatial navigation systems and is one of the first areas to be affected in dementias featuring tau pathology, such as Alzheimer's disease and frontotemporal dementia. Electrophysiological recordings from principle cells of medial entorhinal cortex (layer II stellate cells, mEC-SCs) demonstrate a number of key identifying properties including subthreshold oscillations in the theta (4-12 Hz) range and clustered action potential firing. These single cell properties are correlated with network activity such as grid firing and coupling between theta and gamma rhythms, suggesting they are important for spatial memory. As such, experimental models of dementia have revealed disruption of organised dorsoventral gradients in clustered action potential firing. To better understand the mechanisms underpinning these different dynamics, we study a conductance based model of mEC-SCs. We demonstrate that the model, driven by extrinsic noise, can capture quantitative differences in clustered action potential firing patterns recorded from experimental models of tau pathology and healthy animals. The differential equation formulation of our model allows us to perform numerical bifurcation analyses in order to uncover the dynamic mechanisms underlying these patterns. We show that clustered dynamics can be understood as subcritical Hopf/homoclinic bursting in a fast-slow system where the slow sub-system is governed by activation of the persistent sodium current and inactivation of the slow A-type potassium current. In the full system, we demonstrate that clustered firing arises via flip bifurcations as conductance parameters are varied. Our model analyses confirm the experimentally suggested hypothesis that the breakdown of clustered dynamics in disease occurs via increases in AHP conductance.The contribution of MG, KTR and JB was generously supported by a Wellcome Trust Institutional Strategic Support Award (WT105618MA). MG and KT gratefully acknowledge the financial support of the EPSRC via grant EP/N014391/1. LT’s doctoral studentship is supported by the Alzheimer’s Society in partnership with the Garfield Weston Foundation (grant reference 231). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript

Multi-Scale Mathematical Modelling of Brain Networks in Alzheimer's Disease

Perturbations to brain network dynamics on a range of spatial and temporal scales are believed to underpin neurological disorders such as Alzheimer’s disease (AD). This thesis combines quantitative data analysis with tools such as dynamical systems and graph theory to understand how the network dynamics of the brain are altered in AD and experimental models of related pathologies. Firstly, we use a biophysical neuron model to elucidate ionic mechanisms underpinning alterations to the dynamics of principal neurons in the brain’s spatial navigation systems in an animal model of tauopathy. To uncover how synaptic deficits result in alterations to brain dynamics, we subsequently study an animal model featuring local and long-range synaptic degeneration. Synchronous activity (functional connectivity; FC) between neurons within a region of the cortex is analysed using two-photon calcium imaging data. Long-range FC between regions of the brain is analysed using EEG data. Furthermore, a computational model is used to study relationships between networks on these different spatial scales. The latter half of this thesis studies EEG to characterize alterations to macro-scale brain dynamics in clinical AD. Spectral and FC measures are correlated with cognitive test scores to study the hypothesis that impaired integration of the brain’s processing systems underpin cognitive impairment in AD. Whole brain computational modelling is used to gain insight into the role of spectral slowing on FC, and elucidate potential synaptic mechanisms of FC differences in AD. On a finer temporal scale, microstate analyses are used to identify changes to the rapid transitioning behaviour of the brain’s resting state in AD. Finally, the electrophysiological signatures of AD identified throughout the thesis are combined into a predictive model which can accurately separate people with AD and healthy controls based on their EEG, results which are validated on an independent patient cohort. Furthermore, we demonstrate in a small preliminary cohort that this model is a promising tool for predicting future conversion to AD in patients with mild cognitive impairment

Multistability in Bursting Patterns in a Model of a Multifunctional Central Pattern Generator.

A multifunctional central pattern generator (CPG) can produce bursting polyrhythms that determine locomotive activity in an animal: for example, swimming and crawling in a leech. Each rhythm corresponds to a specific attractor of the CPG. We employ a Hodgkin-Huxley type model of a bursting leech heart interneuron, and connect three such neurons by fast inhibitory synapses to form a ring. This network motif exhibits multistable co-existing bursting rhythms. The problem of determining rhythmic outcomes is reduced to an analysis of fixed points of Poincare mappings and their attractor basins, in a phase plane defined by the interneurons\u27 phase differences along bursting orbits. Using computer assisted analysis, we examine stability, bifurcations of attractors, and transformations of their basins in the phase plane. These structures determine the global bursting rhythms emitted by the CPG. By varying the coupling synaptic strength, we examine the dynamics and patterns produced by inhibitory networks

Modelling human choices: MADeM and decision‑making

Research supported by FAPESP 2015/50122-0 and DFG-GRTK 1740/2. RP and AR are also part of the Research, Innovation and Dissemination Center for Neuromathematics FAPESP grant (2013/07699-0). RP is supported by a FAPESP scholarship (2013/25667-8). ACR is partially supported by a CNPq fellowship (grant 306251/2014-0)

29th Annual Computational Neuroscience Meeting: CNS*2020

Meeting abstracts This publication was funded by OCNS. The Supplement Editors declare that they have no competing interests. Virtual | 18-22 July 202

DEVELOPMENT OF A CEREBELLAR MEAN FIELD MODEL: THE THEORETICAL FRAMEWORK, THE IMPLEMENTATION AND THE FIRST APPLICATION

Brain modeling constantly evolves to improve the accuracy of the simulated brain dynamics with the ambitious aim to build a digital twin of the brain. Specific models tuned on brain regions specific features empower the brain simulations introducing bottom-up physiology properties into data-driven simulators. Despite the cerebellum contains 80 % of the neurons and is deeply involved in a wide range of functions, from sensorimotor to cognitive ones, a specific cerebellar model is still missing. Furthermore, its quasi-crystalline multi-layer circuitry deeply differs from the cerebral cortical one, therefore is hard to imagine a unique general model suitable for the realistic simulation of both cerebellar and cerebral cortex. The present thesis tackles the challenge of developing a specific model for the cerebellum. Specifically, multi-neuron multi-layer mean field (MF) model of the cerebellar network, including Granule Cells, Golgi Cells, Molecular Layer Interneurons, and Purkinje Cells, was implemented, and validated against experimental data and the corresponding spiking neural network microcircuit model. The cerebellar MF model was built using a system of interdependent equations, where the single neuronal populations and topological parameters were captured by neuron-specific inter- dependent Transfer Functions. The model time resolution was optimized using Local Field Potentials recorded experimentally with high-density multielectrode array from acute mouse cerebellar slices. The present MF model satisfactorily captured the average discharge of different microcircuit neuronal populations in response to various input patterns and was able to predict the changes in Purkinje Cells firing patterns occurring in specific behavioral conditions: cortical plasticity mapping, which drives learning in associative tasks, and Molecular Layer Interneurons feed-forward inhibition, which controls Purkinje Cells activity patterns. The cerebellar multi-layer MF model thus provides a computationally efficient tool that will allow to investigate the causal relationship between microscopic neuronal properties and ensemble brain activity in health and pathological conditions. Furthermore, preliminary attempts to simulate a pathological cerebellum were done in the perspective of introducing our multi-layer cerebellar MF model in whole-brain simulators to realize patient-specific treatments, moving ahead towards personalized medicine. Two preliminary works assessed the relevant impact of the cerebellum on whole-brain dynamics and its role in modulating complex responses in causal connected cerebral regions, confirming that a specific model is required to further investigate the cerebellum-on- cerebrum influence. The framework presented in this thesis allows to develop a multi-layer MF model depicting the features of a specific brain region (e.g., cerebellum, basal ganglia), in order to define a general strategy to build up a pool of biology grounded MF models for computationally feasible simulations. Interconnected bottom-up MF models integrated in large-scale simulators would capture specific features of different brain regions, while the applications of a virtual brain would have a substantial impact on the reality ranging from the characterization of neurobiological processes, subject-specific preoperative plans, and development of neuro-prosthetic devices