Data-driven integration of hippocampal CA1 synaptic physiology in silico

The anatomy and physiology of monosynaptic connections in rodent hippocampal CA1 have been extensively studied in recent decades. Yet, the resulting knowledge remains disparate and difficult to reconcile. Here, we present a data-driven approach to integrate the current state-of-the-art knowledge on the synaptic anatomy and physiology of rodent hippocampal CA1, including axo-dendritic innervation patterns, number of synapses per connection, quantal conductances, neurotransmitter release probability, and short-term plasticity into a single coherent resource. First, we undertook an extensive literature review of paired recordings of hippocampal neurons and compiled experimental data on their synaptic anatomy and physiology. The data collected in this manner is sparse and inhomogeneous due to the diversity of experimental techniques used by different groups, which necessitates the need for an integrative framework to unify these data. To this end, we extended a previously developed workflow for the neocortex to constrain a unifying in silico reconstruction of the synaptic physiology of CA1 connections. Our work identifies gaps in the existing knowledge and provides a complementary resource toward a more complete quantification of synaptic anatomy and physiology in the rodent hippocampal CA1 region.

Intrinsic morphological diversity of thick-tufted layer 5 pyramidal neurons ensures robust and invariant properties of in silico synaptic connections

Henry Markram, Srikanth Ramaswamy, Felix Schürmann, Yun Wang

The morphology of neocortical pyramidal neurons is not only highly characteristic but also displays an intrinsic diversity that renders each neuron morphologically unique. We investigated the significance of this intrinsic morphological diversity in in silico networks composed of thick-tufted layer 5 (TTL5) pyramidal neurons, by comparing the in silico and in vitro properties of TTL5 synaptic connections. The synaptic locations of in silico connections were determined by placing 3D reconstructed TTL5 neurons randomly in a volume equivalent to that of layer 5 in the juvenile rat somatosensory cortex and using a 'collision-detection' algorithm to identify the incidental loci of axo-dendritic overlap. The activation time of the modelled synapses and their biophysical properties were characterized based on experimental measurements. We found that the anatomical loci of synapses and the physiological properties of the somatically recorded EPSPs closely matched those recorded experimentally without the need for any fine-tuning. Furthermore, perturbations to both the physiological or anatomical parameters of the model did not alter the average physiological properties of the population of modelled synaptic connections. This microcircuit-level robust behaviour was due to the intrinsic diversity of the morphology of pyramidal neurons in the microcircuit. We conclude that synaptic transmission in a network of TTL5 neurons is highly invariant across microcircuits suggesting that intrinsic diversity is a mechanism to ensure the same average synaptic properties in different animals of the same species. Finally, we show that the average physiological properties of the TTL5 microcircuit are surprisingly robust to anatomical and physiological perturbations also partly due to the intrinsic diversity of pyramidal neuron morphology.

2012

Emergent Properties of in silico Synaptic Transmission in a Model of the Rat Neocortical Column

Srikanth Ramaswamy

The cerebral cortex occupies nearly 80% of the entire volume of the mammalian brain and is thought to subserve higher cognitive functions like memory, attention and sensory perception. The neocortex is the newest part in the evolution of the cerebral cortex and is perhaps the most intricate brain region ever studied. The neocortical microcircuit is the smallest Œecosystem‚ of the neocortex that consists of a rich assortment of neurons, which are diverse in both their morphological and electrical properties. In the neocortical microcircuit, neurons are horizontally arranged in 6 distinct sheets called layers. The fundamental operating unit of the neocortical microcircuit is believed to be the Neocortical Column (NCC). Functionally, a single NCC is an arrangement of thousands of neurons in a vertical fashion spanning across all the 6 layers. The structure of the entire neocortex arises from a repeated and stereotypical arrangement of several thousands of such columns, where neurons transmit information to each other through specialized points of information transfer called synapses. The dynamics of synaptic transmission can be as diverse as the neurons defining a connection and are crucial to foster the functional properties of the neocortical microcircuit. The Blue Brain Project (BBP) is the first comprehensive endeavour to build a unifying model of the NCC by systematic data integration and biologically detailed simulations. Through the past 5 years, the BBP has built a facility for a data-constraint driven approach towards modelling and integrating biological information across multiple levels of complexity. Guided by fundamental principles derived from biological experiments, the BBP simulation toolchain has undergone a process of continuous refinement to facilitate the frequent construction of detailed in silico models of the NCC. The focus of this thesis lies in characterizing the functional properties of in silico synaptic transmission by incorporating principles of synaptic communication derived through biological experiments. In order to study in silico synaptic transmission it is crucial to gain an understanding of the key players influencing the manner in which synaptic signals are processed in the neocortical microcircuit - ion channel kinetics and distribution profiles, single neuron models and dynamics of synaptic pathways. First, by means of exhaustive literature survey, I identified ion channel kinetics and their distribution profiles on neocortical neurons to build in silico ion channel models. Thereafter, I developed a prototype framework to analyze the somatic and dendritic features of single neuron models constrained by ion channel kinetics. Finally, within a simulation framework integrating the ion channels, single neuron models and dynamics of synaptic transmission, I replicated in vitro experimental protocols in silico, to characterize the transmission properties of monosynaptic connections. These synaptic connections, arising from the axo-dendritric apposition of neuronal arbours were sampled across many instances of in silico NCC models constructed a priori through the BBP simulation toolchain. In this thesis, I show that when principles of synaptic transmission derived from in vitro experiments are incorporated to model in silico synaptic connections, the resulting anatomy and physiology of synaptic connections modelled from elementary biological rules closely match in vitro data. This thesis work demonstrates that the average synaptic response properties in silico are robust to perturbations in the anatomical and physiological properties of modelled connections in the local neocortical microcircuit. A fundamental discovery through this thesis is an insight into the function of the local neocortical microcircuit by examining the effect of morphological diversity on in silico synaptic transmission. I demonstrate here that intrinsic morphological diversity confers an invariance to the average synaptic response properties in silico in the local neocortical microcircuit, termed "microcircuit level robustness and invariance".

EPFL2011

Uncertainty reduction in biochemical kinetic models: Enforcing desired model properties

Vassily Hatzimanikatis, Ljubisa Miskovic, Michaël Roger Germain Moret

A persistent obstacle for constructing kinetic models of metabolism is uncertainty in the kinetic properties of enzymes. Currently, available methods for building kinetic models can cope indirectly with uncertainties by integrating data from different biological levels and origins into models. In this study, we use the recently proposed computational approach iSCHRUNK (in Silico Approach to Characterization and Reduction of Uncertainty in the Kinetic Models), which combines Monte Carlo parameter sampling methods and machine learning techniques, in the context of Bayesian inference. Monte Carlo parameter sampling methods allow us to exploit synergies between different data sources and generate a population of kinetic models that are consistent with the available data and physicochemical laws. The machine learning allows us to data-mine the a priori generated kinetic parameters together with the integrated datasets and derive posterior distributions of kinetic parameters consistent with the observed physiology. In this work, we used iSCHRUNK to address a design question: can we identify which are the kinetic parameters and what are their values that give rise to a desired metabolic behavior? Such information is important for a wide variety of studies ranging from biotechnology to medicine. To illustrate the proposed methodology, we performed Metabolic Control Analysis, computed the flux control coefficients of the xylose uptake (XTR), and identified parameters that ensure a rate improvement of XTR in a glucose-xylose co-utilizing S. cerevisiae strain. Our results indicate that only three kinetic parameters need to be accurately characterized to describe the studied physiology, and ultimately to design and control the desired responses of the metabolism. This framework paves the way for a new generation of methods that will systematically integrate the wealth of available omics data and efficiently extract the information necessary for metabolic engineering and synthetic biology decisions.

2019