Neuronal population models reveal specific linear conductance controllers sufficient to rescue preclinical disease phenotypes.

Preclinical drug candidates are screened for their ability to ameliorate in vitro neuronal electrophysiology, and go/no-go decisions progress drugs to clinical trials based on population means across cells and animals. However, these measures do not mitigate clinical endpoint risk. Population-based modeling captures variability across multiple electrophysiological measures from healthy, disease, and drug phenotypes. We pursued optimizing therapeutic targets by identifying coherent sets of ion channel target modulations for recovering heterogeneous wild-type (WT) population excitability profiles from a heterogeneous Huntington's disease (HD) population. Our approach combines mechanistic simulations with population modeling of striatal neurons using evolutionary optimization algorithms to design 'virtual drugs'. We introduce efficacy metrics to score populations and rank virtual drug candidates. We found virtual drugs using heuristic approaches that performed better than single target modulators and standard classification methods. We compare a real drug to virtual candidates and demonstrate a novel in silico triaging method.

Transient, Consequential Increases in Extracellular Potassium Ions Accompany Channelrhodopsin2 Excitation.

Channelrhodopsin2 (ChR2) optogenetic excitation is widely used to study neurons, astrocytes, and circuits. Using complementary approaches in situ and in vivo, we found that ChR2 stimulation leads to significant transient elevation of extracellular potassium ions by ∼5 mM. Such elevations were detected in ChR2-expressing mice, following local in vivo expression of ChR2(H134R) with adeno-associated viruses (AAVs), in different brain areas and when ChR2 was expressed in neurons or astrocytes. In particular, ChR2-mediated excitation of striatal astrocytes was sufficient to increase medium spiny neuron (MSN) excitability and immediate early gene expression. The effects on MSN excitability were recapitulated in silico with a computational MSN model and detected in vivo as increased action potential firing in awake, behaving mice. We show that transient, physiologically consequential increases in extracellular potassium ions accompany ChR2 optogenetic excitation. This coincidental effect may be important to consider during astrocyte studies employing ChR2 to interrogate neural circuits and animal behavior.

Synaptic Efficacy as a Function of Ionotropic Receptor Distribution: A Computational Study.

Glutamatergic synapses are the most prevalent functional elements of information processing in the brain. Changes in pre-synaptic activity and in the function of various post-synaptic elements contribute to generate a large variety of synaptic responses. Previous studies have explored postsynaptic factors responsible for regulating synaptic strength variations, but have given far less importance to synaptic geometry, and more specifically to the subcellular distribution of ionotropic receptors. We analyzed the functional effects resulting from changing the subsynaptic localization of ionotropic receptors by using a hippocampal synaptic computational framework. The present study was performed using the EONS (Elementary Objects of the Nervous System) synaptic modeling platform, which was specifically developed to explore the roles of subsynaptic elements as well as their interactions, and that of synaptic geometry. More specifically, we determined the effects of changing the localization of ionotropic receptors relative to the presynaptic glutamate release site, on synaptic efficacy and its variations following single pulse and paired-pulse stimulation protocols. The results indicate that changes in synaptic geometry do have consequences on synaptic efficacy and its dynamics.

Influence of ionotropic receptors localization on glutamatergic synaptic response to paired-pulse stimulation protocol.

Paired-pulse protocol is a stimulation pattern that is often used to characterize short-term changes in synaptic potency. Responses to such protocol often yield varying results, going from a depressing response to a facilitated one following the second pulse. Similarly, experimental results have shown that synaptic structures are dynamic and receptors move along the postsynaptic membrane. The present study provides insights on the impact of glutamatergic receptors localization with respect to the neurotransmitters release site on the postsynaptic currents measured; it also proposes an explanation on the diversity of responses observed experimentally. The platform we used is the EONS/RHENOMS modeling platform widely described in the literature, which encompasses a multitude of highly detailed subsynaptic elements to most faithfully replicate synaptic function.

Insights on synaptic paired-pulse response using parametric and non-parametric models.

Paired-pulse protocol is a well-established stimulation pattern used to characterize short-term changes in synaptic potency. Due to the experimental difficulty in accessing and measuring responses and interactions between subsynaptic elements, understanding the mechanisms that shape synaptic response is extremely challenging. We already proposed to address this issue and gain insights on the matter using a complex integrated modeling platform called EONS (Elementary Objects of the Nervous System). The use of this parametric platform provided us with insightful information on the subsynaptic components and how their interactions shape synaptic dynamics. We herein propose to add and combine a non-parametric model to (i) simplify the modeling framework, the number of underlying parameters and the overall computational complexity while faithfully maintaining the desirable synaptic behavior and (ii) provide a clear and concise framework to characterize AMPA and NMDA contributions to the observed paired-pulse responses.

A computational model to investigate astrocytic glutamate uptake influence on synaptic transmission and neuronal spiking.

Over the past decades, our view of astrocytes has switched from passive support cells to active processing elements in the brain. The current view is that astrocytes shape neuronal communication and also play an important role in many neurodegenerative diseases. Despite the growing awareness of the importance of astrocytes, the exact mechanisms underlying neuron-astrocyte communication and the physiological consequences of astrocytic-neuronal interactions remain largely unclear. In this work, we define a modeling framework that will permit to address unanswered questions regarding the role of astrocytes. Our computational model of a detailed glutamatergic synapse facilitates the analysis of neural system responses to various stimuli and conditions that are otherwise difficult to obtain experimentally, in particular the readouts at the sub-cellular level. In this paper, we extend a detailed glutamatergic synaptic model, to include astrocytic glutamate transporters. We demonstrate how these glial transporters, responsible for the majority of glutamate uptake, modulate synaptic transmission mediated by ionotropic AMPA and NMDA receptors at glutamatergic synapses. Furthermore, we investigate how these local signaling effects at the synaptic level are translated into varying spatio-temporal patterns of neuron firing. Paired pulse stimulation results reveal that the effect of astrocytic glutamate uptake is more apparent when the input inter-spike interval is sufficiently long to allow the receptors to recover from desensitization. These results suggest an important functional role of astrocytes in spike timing dependent processes and demand further investigation of the molecular basis of certain neurological diseases specifically related to alterations in astrocytic glutamate uptake, such as epilepsy.

Influence of ionotropic receptor location on their dynamics at glutamatergic synapses.

In this paper we study the effects of the location of ionotropic receptors, especially AMPA and NMDA receptors, on their function at excitatory glutamatergic synapses. As few computational models only allow to evaluate the influence of receptor location on state transition and receptor dynamics, we present an elaborate computational model of a glutamatergic synapse that takes into account detailed parametric models of ionotropic receptors along with glutamate diffusion within the synaptic cleft. Our simulation results underscore the importance of the wide spread distribution of AMPA receptors which is required to avoid massive desensitization of these receptors following a single glutamate release event while NMDA receptor location is potentially optimal relative to the glutamate release site thus, emphasizing the contribution of location dependent effects of the two major ionotropic receptors to synaptic efficacy.

Modeling of the nervous system: from modulation of glutamatergic and gabaergic molecular dynamics to neuron spiking activity.

One of the fundamental characteristics of the brain is its hierarchical and temporal organization: scales in both space and time must be considered to fully grasp the system's underlying mechanisms and their impact on brain function. Complex interactions taking place at the molecular level regulate neuronal activity that further modifies the function of millions of neurons connected by trillions of synapses, ultimately giving rise to complex function and behavior at the system level. Likewise, the spatial complexity is accompanied by a complex temporal integration of events taking place at the microsecond scale leading to slower changes occurring at the second, minute and hour scales. These integrations across hierarchies of the nervous system are sufficiently complex to have impeded the development of routine multi-level modeling methodologies. The present study describes an example of our multiscale efforts to rise from the biomolecular level to the neuron level. We more specifically describe how we integrate biomolecular mechanisms taking place at glutamatergic and gabaergic synapses and integrate them to study the impact of these modifications on spiking activity of a CA1 pyramidal cell in the hippocampus.

Simulation of postsynaptic glutamate receptors reveals critical features of glutamatergic transmission.

Activation of several subtypes of glutamate receptors contributes to changes in postsynaptic calcium concentration at hippocampal synapses, resulting in various types of changes in synaptic strength. Thus, while activation of NMDA receptors has been shown to be critical for long-term potentiation (LTP) and long term depression (LTD) of synaptic transmission, activation of metabotropic glutamate receptors (mGluRs) has been linked to either LTP or LTD. While it is generally admitted that dynamic changes in postsynaptic calcium concentration represent the critical elements to determine the direction and amplitude of the changes in synaptic strength, it has been difficult to quantitatively estimate the relative contribution of the different types of glutamate receptors to these changes under different experimental conditions. Here we present a detailed model of a postsynaptic glutamatergic synapse that incorporates ionotropic and mGluR type I receptors, and we use this model to determine the role of the different receptors to the dynamics of postsynaptic calcium with different patterns of presynaptic activation. Our modeling framework includes glutamate vesicular release and diffusion in the cleft and a glutamate transporter that modulates extracellular glutamate concentration. Our results indicate that the contribution of mGluRs to changes in postsynaptic calcium concentration is minimal under basal stimulation conditions and becomes apparent only at high frequency of stimulation. Furthermore, the location of mGluRs in the postsynaptic membrane is also a critical factor, as activation of distant receptors contributes significantly less to calcium dynamics than more centrally located ones. These results confirm the important role of glutamate transporters and of the localization of mGluRs in postsynaptic sites in their signaling properties, and further strengthen the notion that mGluR activation significantly contributes to postsynaptic calcium dynamics only following high-frequency stimulation. They also provide a new tool to analyze the interactions between metabotropic and ionotropic glutamate receptors.

Integrated multiscale modeling of the nervous system: predicting changes in hippocampal network activity by a positive AMPA receptor modulator.

One of the fundamental characteristics of the brain is its hierarchical organization. Scales in both space and time that must be considered when integrating across hierarchies of the nervous system are sufficiently great as to have impeded the development of routine multilevel modeling methodologies. Complex molecular interactions at the level of receptors and channels regulate activity at the level of neurons; interactions between multiple populations of neurons ultimately give rise to complex neural systems function and behavior. This spatial complexity takes place in the context of a composite temporal integration of multiple, different events unfolding at the millisecond, second, minute, hour, and longer time scales. In this study, we present a multiscale modeling methodology that integrates synaptic models into single neuron, and multineuron, network models. We have applied this approach to the specific problem of how changes at the level of kinetic parameters of a receptor-channel model are translated into changes in the temporal firing pattern of a single neuron, and ultimately, changes in the spatiotemporal activity of a network of neurons. These results demonstrate how this powerful methodology can be applied to understand the effects of a given local process within multiple hierarchical levels of the nervous system.

Modeling of the nervous system: from molecular dynamics and synaptic modulation to neuron spiking activity.

The brain is a perfect example of an integrated multi-scale system, as the complex interactions taking place at the molecular level regulate neuronal activity that further modifies the function of millions of neurons connected by trillions of synapses, ultimately giving rise to complex function and behavior at the system level. Likewise, the spatial complexity is accompanied by a complex temporal integration of events taking place at the microsecond scale leading to slower changes occurring at the second, minute and hour scales. In the present study we illustrate our approach to model and simulate the spatio-temporal complexity of the nervous system by developing a multi-scale model integrating synaptic models into the neuronal and ultimately network levels. We apply this approach to a concrete example and demonstrate how changes at the level of kinetic parameters of a receptor model are translated into significant changes in the firing of a pyramidal neuron. These results illustrate the abilities of our modeling approach and support its direct application to the evaluation of the effects of drugs, from functional target to integrated system.

Modeling neuron-glia interactions: from parametric model to neuromorphic hardware.

Recent experimental evidence suggests that glial cells are more than just supporting cells to neurons - they play an active role in signal transmission in the brain. We herein propose to investigate the importance of these mechanisms and model neuron-glia interactions at synapses using three approaches: A parametric model that takes into account the underlying mechanisms of the physiological system, a non-parametric model that extracts its input-output properties, and an ultra-low power, fast processing, neuromorphic hardware model. We use the EONS (Elementary Objects of the Nervous System) platform, a highly elaborate synaptic modeling platform to investigate the influence of astrocytic glutamate transporters on postsynaptic responses in the detailed micro-environment of a tri-partite synapse. The simulation results obtained using EONS are then used to build a non-parametric model that captures the essential features of glutamate dynamics. The structure of the non-parametric model we use is specifically designed for efficient hardware implementation using ultra-low power subthreshold CMOS building blocks. The utilization of the approach described allows us to build large-scale models of neuron/glial interaction and consequently provide useful insights on glial modulation during normal and pathological neural function.

Paired-pulse stimulation at glutamatergic synapses - pre- and postsynaptic components.

Paired-pulse stimulation is a standard protocol that has been used for decades to characterize dynamic systems: the differences in responses to two sequential identical stimuli as a function of inter-stimulus interval provide quantitative information on the dynamics of the system. In neuroscience, the paired-pulse protocol is also widely used at multiple levels of analysis, from behavioral conditioning to synaptic plasticity, and in particular to define the biomolecular mechanism of learning and memory. In a system as small and complex as synapses, it is extremely challenging - if not impossible - to experimentally gain access to the multitude of possible readouts. In the present study, we first introduce a computational synaptic modeling platform that we developed and calibrated based on experimental data from both our laboratories and a variety of publications. We then show how this platform allows not only to replicate experimental data, but also to go beyond technological boundaries and investigate the main parameters responsible for regulation of synaptic transmission and plasticity. The results provide critical information regarding the respective role of various subsynaptic processes and of their interactions. Additionally, this approach can strengthen our understanding of potential dysfunctions (pathologies) and suggest potential approaches to re-establish normal function.

Modeling glutamatergic synapses: insights into mechanisms regulating synaptic efficacy.

The hippocampal formation is critically involved for the long-term storage of various forms of information, and it is widely believed that the phenomenon of long-term potentiation (LTP) of synaptic transmission is a molecular/cellular mechanism participating in memory formation. Although several high level models of hippocampal function have been developed, they do not incorporate detailed molecular information of the type necessary to understand the contribution of individual molecular events to the mechanisms underlying LTP and learning and memory. We are therefore developing new technological tools based on mathematical modeling and computer simulation of the molecular processes taking place in realistic biological networks to reach such an understanding. This article briefly summarizes the approach we are using and illustrates it by presenting data regarding the effects of changing the number of AMPA receptors on various features of glutamatergic transmission, including NMDA receptor-mediated responses and paired-pulse facilitation. We conclude by discussing the significance of these results and providing some ideas for future directions with this approach.