Integration of biochemical and electrical signaling-multiscale model of the medium spiny neuron of the striatum

Comparison of rule- and ordinary differential equation-based dynamic model of DARPP-32 signalling network

Dynamic modelling has considerably improved our understanding of complex molecular mechanisms. Ordinary differential equations (ODEs) are the most detailed and popular approach to modelling the dynamics of molecular systems. However, their application in signalling networks, characterised by multi-state molecular complexes, can be prohibitive. Contemporary modelling methods, such as rule- based (RB) modelling, have addressed these issues. The advantages of RB modelling over ODEs have been presented and discussed in numerous reviews. In this study, we conduct a direct comparison of the time courses of a molecular system founded on the same reaction network but encoded in the two frameworks. To make such a comparison, a set of reactions that underlie an ODE model was manually encoded in the Kappa language, one of the RB implementations. A comparison of the models was performed at the level of model specification and dynamics, acquired through model simulations. In line with previous reports, we confirm that the Kappa model recapitulates the general dynamics of its ODE counterpart with minor differences. These occur when molecules have multiple sites binding the same interactor. Furthermore, activation of these molecules in the RB model is slower than in the ODE one. As reported for other molecular systems, we find that, also for the DARPP-32 reaction network, the RB representation offers a more expressive and flexible syntax that facilitates access to fine details of the model, easing model reuse. In parallel with these analyses, we report a refactored model of the DARPP-32 interaction network that can serve as a canvas for the development of more complex dynamic models to study this important molecular system.

A computational model of dopaminergic modulation of hippocampal Schaffer collateral-CA1 long-term plasticity

Dopamine plays a critical role in modulating the long-term synaptic plasticity of the hippocampal Schaffer collateral-CA1 pyramidal neuron synapses (SC-CA1), a widely accepted cellular model of learning and memory. Limited results from hippocampal slice experiments over the last four decades have shown that the timing of the activation of dopamine D1/D5 receptors relative to a high/low-frequency stimulation (HFS/LFS) in SC-CA1 synapses regulates the modulation of HFS/LFS-induced long-term potentiation/depression (LTP/LTD) in these synapses. However, the existing literature lacks a complete picture of how various concentrations of D1/D5 agonists and the relative timing between the activation of D1/D5 receptors and LTP/LTD induction by HFS/LFS, affect the spatiotemporal modulation of SC-CA1 synaptic dynamics. In this paper, we have developed a computational model, a first of its kind, to make quantitative predictions of the temporal dose-dependent modulation of the HFS/LFS induced LTP/LTD in SC-CA1 synapses by various D1/D5 agonists. Our model combines the biochemical effects with the electrical effects at the electrophysiological level. We have estimated the model parameters from the published electrophysiological data, available from diverse hippocampal CA1 slice experiments, in a Bayesian framework. Our modeling results demonstrate the capability of our model in making quantitative predictions of the available experimental results under diverse HFS/LFS protocols. The predictions from our model show a strong nonlinear dependency of the modulated LTP/LTD by D1/D5 agonists on the relative timing between the activated D1/D5 receptors and the HFS/LFS protocol and the applied concentration of D1/D5 agonists.

SBML Level 3: an extensible format for the exchange and reuse of biological models

Systems biology has experienced dramatic growth in the number, size, and complexity of computational models. To reproduce simulation results and reuse models, researchers must exchange unambiguous model descriptions. We review the latest edition of the Systems Biology Markup Language (SBML), a format designed for this purpose. A community of modelers and software authors developed SBML Level 3 over the past decade. Its modular form consists of a core suited to representing reaction-based models and packages that extend the core with features suited to other model types including constraint-based models, reaction-diffusion models, logical network models, and rule-based models. The format leverages two decades of SBML and a rich software ecosystem that transformed how systems biologists build and interact with models. More recently, the rise of multiscale models of whole cells and organs, and new data sources such as single-cell measurements and live imaging, has precipitated new ways of integrating data with models. We provide our perspectives on the challenges presented by these developments and how SBML Level 3 provides the foundation needed to support this evolution.

Effects of location and extent of spine clustering on synaptic integration in striatal medium spiny neurons-a computational study

The nucleus accumbens (NAc) is known widely for its role in the reward circuit, which is dysregulated in a number of psychological disorders. Recent evidence also suggests the contribution of this structure in spatial and gustatory memories. Because of its role in different types of memories, similar to the hippocampus, we assumed the formation of spine clusters, which are engrams of memory, to be present on dendrites of medium spiny neurons (MSNs). We found that the activation of clustered inputs resulted in sublinear summation when clusters were present on the same branch and also when inputs were distributed on different branches. The size, as well as the location of clusters, was found to affect the summation. With an increase in cluster size and distance from soma, the summation was increasingly sublinear. When the temporal integration window was measured for clustered spines, it was found to be narrower as compared to that for a single spine. Also, distally located clusters resulted in a wider temporal window, as compared to proximal clusters. Our results suggest that depending on the location of clusters, the modes of integration will differ in MSNs possessing clustered spines.

Respectful Modeling: Addressing Uncertainty in Dynamic System Models for Molecular Biology

Although there is still some skepticism in the biological community regarding the value and significance of quantitative computational modeling, important steps are continually being taken to enhance its accessibility and predictive power. We view these developments as essential components of an emerging 'respectful modeling' framework which has two key aims: (i) respecting the models themselves and facilitating the reproduction and update of modeling results by other scientists, and (ii) respecting the predictions of the models and rigorously quantifying the confidence associated with the modeling results. This respectful attitude will guide the design of higher-quality models and facilitate the use of models in modern applications such as engineering and manipulating microbial metabolism by synthetic biology.

Reproducibility in Computational Neuroscience Models and Simulations

OBJECTIVE

Like all scientific research, computational neuroscience research must be reproducible. Big data science, including simulation research, cannot depend exclusively on journal articles as the method to provide the sharing and transparency required for reproducibility.

METHODS

Ensuring model reproducibility requires the use of multiple standard software practices and tools, including version control, strong commenting and documentation, and code modularity.

RESULTS

Building on these standard practices, model-sharing sites and tools have been developed that fit into several categories: 1) standardized neural simulators; 2) shared computational resources; 3) declarative model descriptors, ontologies, and standardized annotations; and 4) model-sharing repositories and sharing standards.

CONCLUSION

A number of complementary innovations have been proposed to enhance sharing, transparency, and reproducibility. The individual user can be encouraged to make use of version control, commenting, documentation, and modularity in development of models. The community can help by requiring model sharing as a condition of publication and funding.

SIGNIFICANCE

Model management will become increasingly important as multiscale models become larger, more detailed, and correspondingly more difficult to manage by any single investigator or single laboratory. Additional big data management complexity will come as the models become more useful in interpreting experiments, thus increasing the need to ensure clear alignment between modeling data, both parameters and results, and experiment.

Efficient Integration of Coupled Electrical-Chemical Systems in Multiscale Neuronal Simulations

Multiscale modeling and simulations in neuroscience is gaining scientific attention due to its growing importance and unexplored capabilities. For instance, it can help to acquire better understanding of biological phenomena that have important features at multiple scales of time and space. This includes synaptic plasticity, memory formation and modulation, homeostasis. There are several ways to organize multiscale simulations depending on the scientific problem and the system to be modeled. One of the possibilities is to simulate different components of a multiscale system simultaneously and exchange data when required. The latter may become a challenging task for several reasons. First, the components of a multiscale system usually span different spatial and temporal scales, such that rigorous analysis of possible coupling solutions is required. Then, the components can be defined by different mathematical formalisms. For certain classes of problems a number of coupling mechanisms have been proposed and successfully used. However, a strict mathematical theory is missing in many cases. Recent work in the field has not so far investigated artifacts that may arise during coupled integration of different approximation methods. Moreover, in neuroscience, the coupling of widely used numerical fixed step size solvers may lead to unexpected inefficiency. In this paper we address the question of possible numerical artifacts that can arise during the integration of a coupled system. We develop an efficient strategy to couple the components comprising a multiscale test problem in neuroscience. We introduce an efficient coupling method based on the second-order backward differentiation formula (BDF2) numerical approximation. The method uses an adaptive step size integration with an error estimation proposed by Skelboe (2000). The method shows a significant advantage over conventional fixed step size solvers used in neuroscience for similar problems. We explore different coupling strategies that define the organization of computations between system components. We study the importance of an appropriate approximation of exchanged variables during the simulation. The analysis shows a substantial impact of these aspects on the solution accuracy in the application to our multiscale neuroscientific test problem. We believe that the ideas presented in the paper may essentially contribute to the development of a robust and efficient framework for multiscale brain modeling and simulations in neuroscience.

Volterra representation enables modeling of complex synaptic nonlinear dynamics in large-scale simulations

Chemical synapses are comprised of a wide collection of intricate signaling pathways involving complex dynamics. These mechanisms are often reduced to simple spikes or exponential representations in order to enable computer simulations at higher spatial levels of complexity. However, these representations cannot capture important nonlinear dynamics found in synaptic transmission. Here, we propose an input-output (IO) synapse model capable of generating complex nonlinear dynamics while maintaining low computational complexity. This IO synapse model is an extension of a detailed mechanistic glutamatergic synapse model capable of capturing the input-output relationships of the mechanistic model using the Volterra functional power series. We demonstrate that the IO synapse model is able to successfully track the nonlinear dynamics of the synapse up to the third order with high accuracy. We also evaluate the accuracy of the IO synapse model at different input frequencies and compared its performance with that of kinetic models in compartmental neuron models. Our results demonstrate that the IO synapse model is capable of efficiently replicating complex nonlinear dynamics that were represented in the original mechanistic model and provide a method to replicate complex and diverse synaptic transmission within neuron network simulations.

Quantitative and logic modelling of molecular and gene networks

Behaviours of complex biomolecular systems are often irreducible to the elementary properties of their individual components. Explanatory and predictive mathematical models are therefore useful for fully understanding and precisely engineering cellular functions. The development and analyses of these models require their adaptation to the problems that need to be solved and the type and amount of available genetic or molecular data. Quantitative and logic modelling are among the main methods currently used to model molecular and gene networks. Each approach comes with inherent advantages and weaknesses. Recent developments show that hybrid approaches will become essential for further progress in synthetic biology and in the development of virtual organisms.

Revealing neuronal function through microelectrode array recordings

Microelectrode arrays and microprobes have been widely utilized to measure neuronal activity, both in vitro and in vivo. The key advantage is the capability to record and stimulate neurons at multiple sites simultaneously. However, unlike the single-cell or single-channel resolution of intracellular recording, microelectrodes detect signals from all possible sources around every sensor. Here, we review the current understanding of microelectrode signals and the techniques for analyzing them. We introduce the ongoing advancements in microelectrode technology, with focus on achieving higher resolution and quality of recordings by means of monolithic integration with on-chip circuitry. We show how recent advanced microelectrode array measurement methods facilitate the understanding of single neurons as well as network function.

COMBINE archive and OMEX format: one file to share all information to reproduce a modeling project

BACKGROUND

With the ever increasing use of computational models in the biosciences, the need to share models and reproduce the results of published studies efficiently and easily is becoming more important. To this end, various standards have been proposed that can be used to describe models, simulations, data or other essential information in a consistent fashion. These constitute various separate components required to reproduce a given published scientific result.

RESULTS

We describe the Open Modeling EXchange format (OMEX). Together with the use of other standard formats from the Computational Modeling in Biology Network (COMBINE), OMEX is the basis of the COMBINE Archive, a single file that supports the exchange of all the information necessary for a modeling and simulation experiment in biology. An OMEX file is a ZIP container that includes a manifest file, listing the content of the archive, an optional metadata file adding information about the archive and its content, and the files describing the model. The content of a COMBINE Archive consists of files encoded in COMBINE standards whenever possible, but may include additional files defined by an Internet Media Type. Several tools that support the COMBINE Archive are available, either as independent libraries or embedded in modeling software.

CONCLUSIONS

The COMBINE Archive facilitates the reproduction of modeling and simulation experiments in biology by embedding all the relevant information in one file. Having all the information stored and exchanged at once also helps in building activity logs and audit trails. We anticipate that the COMBINE Archive will become a significant help for modellers, as the domain moves to larger, more complex experiments such as multi-scale models of organs, digital organisms, and bioengineering.

Parametric Anatomical Modeling: a method for modeling the anatomical layout of neurons and their projections

Computational models of neural networks can be based on a variety of different parameters. These parameters include, for example, the 3d shape of neuron layers, the neurons' spatial projection patterns, spiking dynamics and neurotransmitter systems. While many well-developed approaches are available to model, for example, the spiking dynamics, there is a lack of approaches for modeling the anatomical layout of neurons and their projections. We present a new method, called Parametric Anatomical Modeling (PAM), to fill this gap. PAM can be used to derive network connectivities and conduction delays from anatomical data, such as the position and shape of the neuronal layers and the dendritic and axonal projection patterns. Within the PAM framework, several mapping techniques between layers can account for a large variety of connection properties between pre- and post-synaptic neuron layers. PAM is implemented as a Python tool and integrated in the 3d modeling software Blender. We demonstrate on a 3d model of the hippocampal formation how PAM can help reveal complex properties of the synaptic connectivity and conduction delays, properties that might be relevant to uncover the function of the hippocampus. Based on these analyses, two experimentally testable predictions arose: (i) the number of neurons and the spread of connections is heterogeneously distributed across the main anatomical axes, (ii) the distribution of connection lengths in CA3-CA1 differ qualitatively from those between DG-CA3 and CA3-CA3. Models created by PAM can also serve as an educational tool to visualize the 3d connectivity of brain regions. The low-dimensional, but yet biologically plausible, parameter space renders PAM suitable to analyse allometric and evolutionary factors in networks and to model the complexity of real networks with comparatively little effort.

Multiscale modeling and synaptic plasticity

Synaptic plasticity is a major convergence point for theory and computation, and the process of plasticity engages physiology, cell, and molecular biology. In its many manifestations, plasticity is at the hub of basic neuroscience questions about memory and development, as well as more medically themed questions of neural damage and recovery. As an important cellular locus of memory, synaptic plasticity has received a huge amount of experimental and theoretical attention. If computational models have tended to pick specific aspects of plasticity, such as STDP, and reduce them to an equation, some experimental studies are equally guilty of oversimplification each time they identify a new molecule and declare it to be the last word in plasticity and learning. Multiscale modeling begins with the acknowledgment that synaptic function spans many levels of signaling, and these are so tightly coupled that we risk losing essential features of plasticity if we focus exclusively on any one level. Despite the technical challenges and gaps in data for model specification, an increasing number of multiscale modeling studies have taken on key questions in plasticity. These have provided new insights, but importantly, they have opened new avenues for questioning. This review discusses a wide range of multiscale models in plasticity, including their technical landscape and their implications.

Modeling intracellular signaling underlying striatal function in health and disease

Striatum, which is the input nucleus of the basal ganglia, integrates cortical and thalamic glutamatergic inputs with dopaminergic afferents from the substantia nigra pars compacta. The combination of dopamine and glutamate strongly modulates molecular and cellular properties of striatal neurons and the strength of corticostriatal synapses. These actions are performed via intracellular signaling networks, containing several intertwined feedback loops. Understanding the role of dopamine and other neuromodulators requires the development of quantitative dynamical models for describing the intracellular signaling, in order to provide precise unambiguous descriptions and quantitative predictions. Building such models requires integration of data from multiple data sources containing information regarding the molecular interactions, the strength of these interactions, and the subcellular localization of the molecules. Due to the uncertainty, variability, and sparseness of these data, parameter estimation techniques are critical for inferring or constraining the unknown parameters, and sensitivity analysis evaluates which parameters are most critical for a given observed macroscopic behavior. Here, we briefly review the modeling approaches and tools that have been used to investigate biochemical signaling in the striatum, along with some of the models built around striatum. We also suggest a future direction for the development of such models from the, now becoming abundant, high-throughput data.