Strategies for calibrating models of biology.

Computational and mathematical modelling has become a valuable tool for investigating biological systems. Modelling enables prediction of how biological components interact to deliver system-level properties and extrapolation of biological system performance to contexts and experimental conditions where this is unknown. A model's value hinges on knowing that it faithfully represents the biology under the contexts of use, or clearly ascertaining otherwise and thus motivating further model refinement. These qualities are evaluated through calibration, typically formulated as identifying model parameter values that align model and biological behaviours as measured through a metric applied to both. Calibration is critical to modelling but is often underappreciated. A failure to appropriately calibrate risks unrepresentative models that generate erroneous insights. Here, we review a suite of strategies to more rigorously challenge a model's representation of a biological system. All are motivated by features of biological systems, and illustrative examples are drawn from the modelling literature. We examine the calibration of a model against distributions of biological behaviours or outcomes, not only average values. We argue for calibration even where model parameter values are experimentally ascertained. We explore how single metrics can be non-distinguishing for complex systems, with multiple-component dynamic and interaction configurations giving rise to the same metric output. Under these conditions, calibration is insufficiently constraining and the model non-identifiable: multiple solutions to the calibration problem exist. We draw an analogy to curve fitting and argue that calibrating a biological model against a single experiment or context is akin to curve fitting against a single data point. Though useful for communicating model results, we explore how metrics that quantify heavily emergent properties may not be suitable for use in calibration. Lastly, we consider the role of sensitivity and uncertainty analysis in calibration and the interpretation of model results. Our goal in this manuscript is to encourage a deeper consideration of calibration, and how to increase its capacity to either deliver faithful models or demonstrate them otherwise.

ASPASIA: A toolkit for evaluating the effects of biological interventions on SBML model behaviour.

A calibrated computational model reflects behaviours that are expected or observed in a complex system, providing a baseline upon which sensitivity analysis techniques can be used to analyse pathways that may impact model responses. However, calibration of a model where a behaviour depends on an intervention introduced after a defined time point is difficult, as model responses may be dependent on the conditions at the time the intervention is applied. We present ASPASIA (Automated Simulation Parameter Alteration and SensItivity Analysis), a cross-platform, open-source Java toolkit that addresses a key deficiency in software tools for understanding the impact an intervention has on system behaviour for models specified in Systems Biology Markup Language (SBML). ASPASIA can generate and modify models using SBML solver output as an initial parameter set, allowing interventions to be applied once a steady state has been reached. Additionally, multiple SBML models can be generated where a subset of parameter values are perturbed using local and global sensitivity analysis techniques, revealing the model's sensitivity to the intervention. To illustrate the capabilities of ASPASIA, we demonstrate how this tool has generated novel hypotheses regarding the mechanisms by which Th17-cell plasticity may be controlled in vivo. By using ASPASIA in conjunction with an SBML model of Th17-cell polarisation, we predict that promotion of the Th1-associated transcription factor T-bet, rather than inhibition of the Th17-associated transcription factor RORγt, is sufficient to drive switching of Th17 cells towards an IFN-γ-producing phenotype. Our approach can be applied to all SBML-encoded models to predict the effect that intervention strategies have on system behaviour. ASPASIA, released under the Artistic License (2.0), can be downloaded from http://www.york.ac.uk/ycil/software.

Leukocyte Motility Models Assessed through Simulation and Multi-objective Optimization-Based Model Selection.

The advent of two-photon microscopy now reveals unprecedented, detailed spatio-temporal data on cellular motility and interactions in vivo. Understanding cellular motility patterns is key to gaining insight into the development and possible manipulation of the immune response. Computational simulation has become an established technique for understanding immune processes and evaluating hypotheses in the context of experimental data, and there is clear scope to integrate microscopy-informed motility dynamics. However, determining which motility model best reflects in vivo motility is non-trivial: 3D motility is an intricate process requiring several metrics to characterize. This complicates model selection and parameterization, which must be performed against several metrics simultaneously. Here we evaluate Brownian motion, Lévy walk and several correlated random walks (CRWs) against the motility dynamics of neutrophils and lymph node T cells under inflammatory conditions by simultaneously considering cellular translational and turn speeds, and meandering indices. Heterogeneous cells exhibiting a continuum of inherent translational speeds and directionalities comprise both datasets, a feature significantly improving capture of in vivo motility when simulated as a CRW. Furthermore, translational and turn speeds are inversely correlated, and the corresponding CRW simulation again improves capture of our in vivo data, albeit to a lesser extent. In contrast, Brownian motion poorly reflects our data. Lévy walk is competitive in capturing some aspects of neutrophil motility, but T cell directional persistence only, therein highlighting the importance of evaluating models against several motility metrics simultaneously. This we achieve through novel application of multi-objective optimization, wherein each model is independently implemented and then parameterized to identify optimal trade-offs in performance against each metric. The resultant Pareto fronts of optimal solutions are directly contrasted to identify models best capturing in vivo dynamics, a technique that can aid model selection more generally. Our technique robustly determines our cell populations' motility strategies, and paves the way for simulations that incorporate accurate immune cell motility dynamics.

Imaging Immunity in Lymph Nodes: Past, Present and Future.

Immune responses occur as a result of stochastic interactions between a plethora of different cell types and molecules that regulate the migration and function of innate and adaptive immune cells to drive protection from pathogen infection. The trafficking of immune cells into peripheral tissues during inflammation and then subsequent migration to draining lymphoid tissues has been quantitated using radiolabelled immune cells over 40 years ago. However, how these processes lead to efficient immune responses was unclear. Advances in physics (multi-photon), chemistry (probes) and biology (animal models) have provided immunologists with specialized tools to quantify the molecular and cellular mechanisms driving immune function in lymphoid tissues through directly visualising cellular behaviours in 3-dimensions over time. Through the temporal and spatial resolution of multi-photon confocal microscopy immunologists have developed new insights into normal immune homeostasis, host responses to pathogens, anti-tumour immune responses and processes driving development of autoimmune pathologies, by the quantification of the interactions and cellular migration involved in adaptive immune responses. Advances in deep tissue imaging, including new fluorescent proteins, increased resolution, speed of image acquisition, sensitivity, number of signals and improved data analysis techniques have provided unprecedented capacity to quantify immune responses at the single cell level. This quantitative information has facilitated development of high-fidelity mathematical and computational models of immune function. Together this approach is providing new mechanistic understanding of immune responses and new insights into how immune modulators work. Advances in biophysics have therefore revolutionised our understanding of immune function, directly impacting on the development of next generation immunotherapies and vaccines, and is providing the quantitative basis for emerging technology of simulation-guided experimentation and immunotherapeutic design.

Modelling the impact and control of an infectious disease in a plant nursery with infected plant material inputs.

The ornamental plant trade has been identified as a key introduction pathway for plant pathogens. Establishing effective biosecurity measures to reduce the risk of plant pathogen outbreaks in the live plant trade is therefore important. Management of invasive pathogens has been identified as a weakest link public good, and thus is reliant on the actions of individual private agents. This paper therefore provides an analysis of the impact of the private agents' biosecurity decisions on pathogen prevention and control within the plant trade. We model the impact that an infectious disease has on a plant nursery under a constant pressure of potentially infected input plant materials, like seeds and saplings, where the spread of the disease reduces the value of mature plants. We explore six scenarios to understand the influence of three key bioeconomic parameters; the disease's basic reproductive number, the loss in value of a mature plant from acquiring an infection and the cost-effectiveness of restriction. The results characterise the disease dynamics within the nursery and explore the trade-offs and synergies between the optimal level of efforts on restriction strategies (actions to prevent buying infected inputs), and on removal of infected plants in the nursery. For diseases that can be easily controlled, restriction and removal are substitutable strategies. In contrast, for highly infectious diseases, restriction and removal are often found to be complementary, provided that restriction is cost-effective and the optimal level of removal is non-zero.

Utilising a simulation platform to understand the effect of domain model assumptions.

Computational and mathematical modelling approaches are increasingly being adopted in attempts to further our understanding of complex biological systems. This approach can be subjected to strong criticism as substantial aspects of the biological system being captured are not currently known, meaning assumptions need to be made that could have a critical impact on simulation response. We have utilised the CoSMoS process in the development of an agent-based simulation of the formation of Peyer's patches (PP), gut-associated lymphoid organs that have a key role in the initiation of adaptive immune responses to infection. Although the use of genetic tools, imaging technologies and ex vivo culture systems has provided significant insight into the cellular components and associated pathways involved in PP development, interesting questions remain that cannot be addressed using these approaches, and as such well justified assumptions have been introduced into our model to counter this. Here we focus not on the development of the model itself, but instead demonstrate how the resultant simulation can be used to assess how these assumptions impact the simulation response. For example, we consider the impact of our assumption that the migration rate of lymphoid tissue cells into the gut remains constant throughout PP development. We demonstrate that an analysis of the assumptions made in the construction of the domain model may either increase confidence in the model as a representation of the biological system it captures, or may suggest areas where further biological experimentation is required.

A Petri net model of granulomatous inflammation: implications for IL-10 mediated control of Leishmania donovani infection.

Experimental visceral leishmaniasis, caused by infection of mice with the protozoan parasite Leishmania donovani, is characterized by focal accumulation of inflammatory cells in the liver, forming discrete "granulomas" within which the parasite is eventually eliminated. To shed new light on fundamental aspects of granuloma formation and function, we have developed an in silico Petri net model that simulates hepatic granuloma development throughout the course of infection. The model was extensively validated by comparison with data derived from experimental studies in mice, and the model robustness was assessed by a sensitivity analysis. The model recapitulated the progression of disease as seen during experimental infection and also faithfully predicted many of the changes in cellular composition seen within granulomas over time. By conducting in silico experiments, we have identified a previously unappreciated level of inter-granuloma diversity in terms of the development of anti-leishmanial activity. Furthermore, by simulating the impact of IL-10 gene deficiency in a variety of lymphocyte and myeloid cell populations, our data suggest a dominant local regulatory role for IL-10 produced by infected Kupffer cells at the core of the granuloma.

Spartan: a comprehensive tool for understanding uncertainty in simulations of biological systems.

Integrating computer simulation with conventional wet-lab research has proven to have much potential in furthering the understanding of biological systems. Success requires the relationship between simulation and the real-world system to be established: substantial aspects of the biological system are typically unknown, and the abstract nature of simulation can complicate interpretation of in silico results in terms of the biology. Here we present spartan (Simulation Parameter Analysis RToolkit ApplicatioN), a package of statistical techniques specifically designed to help researchers understand this relationship and provide novel biological insight. The tools comprising spartan help identify which simulation results can be attributed to the dynamics of the modelled biological system, rather than artefacts of biological uncertainty or parametrisation, or simulation stochasticity. Statistical analyses reveal the influence that pathways and components have on simulation behaviour, offering valuable biological insight into aspects of the system under study. We demonstrate the power of spartan in providing critical insight into aspects of lymphoid tissue development in the small intestine through simulation. Spartan is released under a GPLv2 license, implemented within the open source R statistical environment, and freely available from both the Comprehensive R Archive Network (CRAN) and http://www.cs.york.ac.uk/spartan. The techniques within the package can be applied to traditional ordinary or partial differential equation simulations as well as agent-based implementations. Manuals, comprehensive tutorials, and example simulation data upon which spartan can be applied are available from the website.

In silico investigation into dendritic cell regulation of CD8Treg mediated killing of Th1 cells in murine experimental autoimmune encephalomyelitis.

BACKGROUND: Experimental autoimmune encephalomyelitis has been used extensively as an animal model of T cell mediated autoimmunity. A down-regulatory pathway through which encephalitogenic CD4Th1 cells are killed by CD8 regulatory T cells (Treg) has recently been proposed. With the CD8Treg cells being primed by dendritic cells, regulation of recovery may be occuring around these antigen presenting cells. CD4Treg cells provide critical help within this process, by licensing dendritic cells to prime CD8Treg cells, however the spatial and temporal aspects of this help in the CTL response is currently unclear. RESULTS: We have previously developed a simulator of experimental autoimmune encephalomyelitis (ARTIMMUS). We use ARTIMMUS to perform novel in silico experimentation regarding the priming of CD8Treg cells by dendritic cells, and the resulting CD8Treg mediated killing of encephalitogenic CD4Th1 cells. Simulations using dendritic cells that present antigenic peptides in a mutually exclusive manner (either MBP or TCR-derived, but not both) suggest that there is no significant reliance on dendritic cells that can prime both encephalitogenic CD4Th1 and Treg cells. Further, in silico experimentation suggests that dynamics of CD8Treg priming are significantly influenced through their spatial competition with CD4Treg cells and through the timing of Qa-1 expression by dendritic cells. CONCLUSION: There is no requirement for the encephalitogenic CD4Th1 cells and cytotoxic CD8Treg cells to be primed by the same dendritic cells. We conjecture that no significant portion of CD4Th1 regulation by Qa-1 restricted CD8Treg cells occurs around individual dendritic cells, and as such, that CD8Treg mediated killing of CD4Th1 cells occurring around dendritic cells is not critical for recovery from the murine autoimmune disease. Furthermore, the timing of the CD4Treg licensing of dendritic cells and the spatial competition between CD4Treg and CD8Treg cells around the dendritic cell is critical for the size of the cytotoxic T lymphocyte response, because dendritic cells have a limited lifespan. If treatments can be found to either speed up the licensing process, or increase the spatial competitiveness of CD8Treg cells, the magnitude of the cytotoxic T lymphocyte response can be increased.

Practical hardware for evolvable robots.

The evolutionary robotics field offers the possibility of autonomously generating robots that are adapted to desired tasks by iteratively optimising across successive generations of robots with varying configurations until a high-performing candidate is found. The prohibitive time and cost of actually building this many robots means that most evolutionary robotics work is conducted in simulation, but to apply evolved robots to real-world problems, they must be implemented in hardware, which brings new challenges. This paper explores in detail the design of an example system for realising diverse evolved robot bodies, and specifically how this interacts with the evolutionary process. We discover that every aspect of the hardware implementation introduces constraints that change the evolutionary space, and exploring this interplay between hardware constraints and evolution is the key contribution of this paper. In simulation, any robot that can be defined by a suitable genetic representation can be implemented and evaluated, but in hardware, real-world limitations like manufacturing/assembly constraints and electrical power delivery mean that many of these robots cannot be built, or will malfunction in operation. This presents the novel challenge of how to constrain an evolutionary process within the space of evolvable phenotypes to only those regions that are practically feasible: the viable phenotype space. Methods of phenotype filtering and repair were introduced to address this, and found to degrade the diversity of the robot population and impede traversal of the exploration space. Furthermore, the degrees of freedom permitted by the hardware constraints were found to be poorly matched to the types of morphological variation that would be the most useful in the target environment. Consequently, the ability of the evolutionary process to generate robots with effective adaptations was greatly reduced. The conclusions from this are twofold. 1) Designing a hardware platform for evolving robots requires different thinking, in which all design decisions should be made with reference to their impact on the viable phenotype space. 2) It is insufficient to just evolve robots in simulation without detailed consideration of how they will be implemented in hardware, because the hardware constraints have a profound impact on the evolutionary space.

A study of error diversity in robotic swarms for task partitioning in foraging tasks.

Often in swarm robotics, an assumption is made that all robots in the swarm behave the same and will have a similar (if not the same) error model. However, in reality, this is not the case, and this lack of uniformity in the error model, and other operations, can lead to various emergent behaviors. This paper considers the impact of the error model and compares robots in a swarm that operate using the same error model (uniform error) against each robot in the swarm having a different error model (thus introducing error diversity). Experiments are presented in the context of a foraging task. Simulation and physical experimental results show the importance of the error model and diversity in achieving the expected swarm behavior.

Simulating CXCR5 Dynamics in Complex Tissue Microenvironments.

To effectively navigate complex tissue microenvironments, immune cells sense molecular concentration gradients using G-protein coupled receptors. However, due to the complexity of receptor activity, and the multimodal nature of chemokine gradients in vivo, chemokine receptor activity in situ is poorly understood. To address this issue, we apply a modelling and simulation approach that permits analysis of the spatiotemporal dynamics of CXCR5 expression within an in silico B-follicle with single-cell resolution. Using this approach, we show that that in silico B-cell scanning is robust to changes in receptor numbers and changes in individual kinetic rates of receptor activity, but sensitive to global perturbations where multiple parameters are altered simultaneously. Through multi-objective optimization analysis we find that the rapid modulation of CXCR5 activity through receptor binding, desensitization and recycling is required for optimal antigen scanning rates. From these analyses we predict that chemokine receptor signaling dynamics regulate migration in complex tissue microenvironments to a greater extent than the total numbers of receptors on the cell surface.

Elucidation of T cell signalling models.

A potential mechanism that allows T cells to reliably discriminate pMHC ligands involves an interplay between kinetic proofreading, negative feedback and a destruction of this negative feedback. We analyse a detailed model of these mechanisms which involves the TCR, SHP1 and ERK. We discover that the behaviour of pSHP1 negative feedback is of primary importance, and particularly the influence of a kinetic proofreading base negative feedback state on pSHP1 dynamics. The CD8 co-receptor is shown to benefit from a kinetic proofreading locking mechanism and is able to overcome pSHP1 negative influences to sensitise a T cell.

Using Emulation to Engineer and Understand Simulations of Biological Systems.

Modeling and simulation techniques have demonstrated success in studying biological systems. As the drive to better capture biological complexity leads to more sophisticated simulators, it becomes challenging to perform statistical analyses that help translate predictions into increased understanding. These analyses may require repeated executions and extensive sampling of high-dimensional parameter spaces: analyses that may become intractable due to time and resource limitations. Significant reduction in these requirements can be obtained using surrogate models, or emulators, that can rapidly and accurately predict the output of an existing simulator. We apply emulation to evaluate and enrich understanding of a previously published agent-based simulator of lymphoid tissue organogenesis, showing an ensemble of machine learning techniques can reproduce results obtained using a suite of statistical analyses within seconds. This performance improvement permits incorporation of previously intractable analyses, including multi-objective optimization to obtain parameter sets that yield a desired response, and Approximate Bayesian Computation to assess parametric uncertainty. To facilitate exploitation of emulation in simulation-focused studies, we extend our open source statistical package, spartan, to provide a suite of tools for emulator development, validation, and application. Overcoming resource limitations permits enriched evaluation and refinement, easing translation of simulator insights into increased biological understanding.

B cell zone reticular cell microenvironments shape CXCL13 gradient formation.

Through the formation of concentration gradients, morphogens drive graded responses to extracellular signals, thereby fine-tuning cell behaviors in complex tissues. Here we show that the chemokine CXCL13 forms both soluble and immobilized gradients. Specifically, CXCL13+ follicular reticular cells form a small-world network of guidance structures, with computer simulations and optimization analysis predicting that immobilized gradients created by this network promote B cell trafficking. Consistent with this prediction, imaging analysis show that CXCL13 binds to extracellular matrix components in situ, constraining its diffusion. CXCL13 solubilization requires the protease cathepsin B that cleaves CXCL13 into a stable product. Mice lacking cathepsin B display aberrant follicular architecture, a phenotype associated with effective B cell homing to but not within lymph nodes. Our data thus suggest that reticular cells of the B cell zone generate microenvironments that shape both immobilized and soluble CXCL13 gradients.

Optimizing amino acid groupings for GPCR classification.

MOTIVATION: There is much interest in reducing the complexity inherent in the representation of the 20 standard amino acids within bioinformatics algorithms by developing a so-called reduced alphabet. Although there is no universally applicable residue grouping, there are numerous physiochemical criteria upon which one can base groupings. Local descriptors are a form of alignment-free analysis, the efficiency of which is dependent upon the correct selection of amino acid groupings. RESULTS: Within the context of G-protein coupled receptor (GPCR) classification, an optimization algorithm was developed, which was able to identify the most efficient grouping when used to generate local descriptors. The algorithm was inspired by the relatively new computational intelligence paradigm of artificial immune systems. A number of amino acid groupings produced by this algorithm were evaluated with respect to their ability to generate local descriptors capable of providing an accurate classification algorithm for GPCRs.

Application of Modeling Approaches to Explore Vaccine Adjuvant Mode-of-Action.

Novel adjuvant technologies have a key role in the development of next-generation vaccines, due to their capacity to modulate the duration, strength and quality of the immune response. The AS01 adjuvant is used in the malaria vaccine RTS,S/AS01 and in the licensed herpes-zoster vaccine (Shingrix) where the vaccine has proven its ability to generate protective responses with both robust humoral and T-cell responses. For many years, animal models have provided insights into adjuvant mode-of-action (MoA), generally through investigating individual genes or proteins. Furthermore, modeling and simulation techniques can be utilized to integrate a variety of different data types; ranging from serum biomarkers to large scale "omics" datasets. In this perspective we present a framework to create a holistic integration of pre-clinical datasets and immunological literature in order to develop an evidence-based hypothesis of AS01 adjuvant MoA, creating a unified view of multiple experiments. Furthermore, we highlight how holistic systems-knowledge can serve as a basis for the construction of models and simulations supporting exploration of key questions surrounding adjuvant MoA. Using the Systems-Biology-Graphical-Notation, a tool for graphical representation of biological processes, we have captured high-level cellular behaviors and interactions, and cytokine dynamics during the early immune response, which are substantiated by a series of diagrams detailing cellular dynamics. Through explicitly describing AS01 MoA we have built a consensus of understanding across multiple experiments, and so we present a framework to integrate modeling approaches into exploring adjuvant MoA, in order to guide experimental design, interpret results and inform rational design of vaccines.

Tissue-specific transcriptomic changes associated with AmBisome® treatment of BALB/c mice with experimental visceral leishmaniasis.

Background: Liposomal amphotericin B (AmBisome®) as a treatment modality for visceral leishmaniasis (VL) has had significant impact on patient care in some but not all regions where VL is endemic.  As the mode of action of AmBisome® in vivo is poorly understood, we compared the tissue-specific transcriptome in drug-treated vs untreated mice with experimental VL.    Methods:  BALB/c mice infected with L. donovani were treated with 8mg/kg AmBisome®, resulting in parasite elimination from liver and spleen over a 7-day period. At day 1 and day 7 post treatment (R x+1 and R x+7), transcriptomic profiling was performed on spleen and liver tissue from treated and untreated mice and uninfected mice.  BALB/c mice infected with M. bovis BCG (an organism resistant to amphotericin B) were analysed to distinguish between direct effects of AmBisome® and those secondary to parasite death.   Results: AmBisome® treatment lead to rapid parasitological clearance.  At R x+1, spleen and liver displayed only 46 and 88 differentially expressed (DE) genes (P<0.05; 2-fold change) respectively. In liver, significant enrichment was seen for pathways associated with TNF, fatty acids and sterol biosynthesis.  At R x+7, the number of DE genes was increased (spleen, 113; liver 400).  In spleen, these included many immune related genes known to be involved in anti-leishmanial immunity. In liver, changes in transcriptome were largely accounted for by loss of granulomas.   PCA analysis indicated that treatment only partially restored homeostasis.  Analysis of BCG-infected mice treated with AmBisome® revealed a pattern of immune modulation mainly targeting macrophage function.   Conclusions: Our data indicate that the tissue response to AmBisome® treatment varies between target organs and that full restoration of homeostasis is not achieved at parasitological cure.  The pathways required to restore homeostasis deserve fuller attention, to understand mechanisms associated with treatment failure and relapse and to promote more rapid restoration of immune competence.

Exploring Self-Repair in a Coupled Spiking Astrocyte Neural Network.

It is now known that astrocytes modulate the activity at the tripartite synapses where indirect signaling via the retrograde messengers, endocannabinoids, leads to a localized self-repairing capability. In this paper, a self-repairing spiking astrocyte neural network (SANN) is proposed to demonstrate a distributed self-repairing capability at the network level. The SANN uses a novel learning rule that combines the spike-timing-dependent plasticity (STDP) and Bienenstock, Cooper, and Munro (BCM) learning rules (hereafter referred to as the BSTDP rule). In this learning rule, the synaptic weight potentiation is not only driven by the temporal difference between the presynaptic and postsynaptic neuron firing times but also by the postsynaptic neuron activity. We will show in this paper that the BSTDP modulates the height of the plasticity window to establish an input-output mapping (in the learning phase) and also maintains this mapping (via self-repair) if synaptic pathways become dysfunctional. It is the functional dependence of postsynaptic neuron firing activity on the height of the plasticity window that underpins how the proposed SANN self-repairs on the fly. The SANN also uses the coupling between the tripartite synapses and γ -GABAergic interneurons. This interaction gives rise to a presynaptic neuron frequency filtering capability that serves to route information, represented as spike trains, to different neurons in the subsequent layers of the SANN. The proposed SANN follows a feedforward architecture with multiple interneuron pathways and astrocytes modulate synaptic activity at the hidden and output neuronal layers. The self-repairing capability will be demonstrated in a robotic obstacle avoidance application, and the simulation results will show that the SANN can maintain learned maneuvers at synaptic fault densities of up to 80% regardless of the fault locations.

GABA Regulation of Burst Firing in Hippocampal Astrocyte Neural Circuit: A Biophysical Model.

It is now widely accepted that glia cells and gamma-aminobutyric acidergic (GABA) interneurons dynamically regulate synaptic transmission and neuronal activity in time and space. This paper presents a biophysical model that captures the interaction between an astrocyte cell, a GABA interneuron and pre/postsynaptic neurons. Specifically, GABA released from a GABA interneuron triggers in astrocytes the release of calcium (Ca 2+) from the endoplasmic reticulum via the inositol 1, 4, 5-trisphosphate (IP 3) pathway. This results in gliotransmission which elevates the presynaptic transmission probability rate (PR) causing weight potentiation and a gradual increase in postsynaptic neuronal firing, that eventually stabilizes. However, by capturing the complex interactions between IP 3, generated from both GABA and the 2-arachidonyl glycerol (2-AG) pathway, and PR, this paper shows that this interaction not only gives rise to an initial weight potentiation phase but also this phase is followed by postsynaptic bursting behavior. Moreover, the model will show that there is a presynaptic frequency range over which burst firing can occur. The proposed model offers a novel cellular level mechanism that may underpin both seizure-like activity and neuronal synchrony across different brain regions.