Understanding molecular signaling cascades in neural disease using multi-resolution models.

If the genome defines the program for the operations of a cell, signaling networks execute it. These cascades of chemical, cell-biological, structural, and trafficking events span milliseconds (e.g., synaptic release) to potentially a lifetime (e.g., stabilization of dendritic spines). In principle almost every aspect of neuronal function, particularly at the synapse, depends on signaling. Thus dysfunction of these cascades, whether through mutations, local dysregulation, or infection, leads to disease. The sheer complexity of these pathways is matched by the range of diseases and the diversity of their phenotypes. In this review, we discuss how to build computational models, how these models are essential to tackle this complexity, and the benefits of using families of models at different levels of detail to understand signaling in health and disease.

Combining hypothesis- and data-driven neuroscience modeling in FAIR workflows.

Modeling in neuroscience occurs at the intersection of different points of view and approaches. Typically, hypothesis-driven modeling brings a question into focus so that a model is constructed to investigate a specific hypothesis about how the system works or why certain phenomena are observed. Data-driven modeling, on the other hand, follows a more unbiased approach, with model construction informed by the computationally intensive use of data. At the same time, researchers employ models at different biological scales and at different levels of abstraction. Combining these models while validating them against experimental data increases understanding of the multiscale brain. However, a lack of interoperability, transparency, and reusability of both models and the workflows used to construct them creates barriers for the integration of models representing different biological scales and built using different modeling philosophies. We argue that the same imperatives that drive resources and policy for data - such as the FAIR (Findable, Accessible, Interoperable, Reusable) principles - also support the integration of different modeling approaches. The FAIR principles require that data be shared in formats that are Findable, Accessible, Interoperable, and Reusable. Applying these principles to models and modeling workflows, as well as the data used to constrain and validate them, would allow researchers to find, reuse, question, validate, and extend published models, regardless of whether they are implemented phenomenologically or mechanistically, as a few equations or as a multiscale, hierarchical system. To illustrate these ideas, we use a classical synaptic plasticity model, the Bienenstock-Cooper-Munro rule, as an example due to its long history, different levels of abstraction, and implementation at many scales.

Adult brain neurons require continual expression of the schizophrenia-risk gene Tcf4 for structural and functional integrity.

The schizophrenia-risk gene Tcf4 has been widely studied in the context of brain development using mouse models of haploinsufficiency, in utero knockdown and embryonic deletion. However, Tcf4 continues to be abundantly expressed in adult brain neurons where its functions remain unknown. Given the importance of Tcf4 in psychiatric diseases, we investigated its role in adult neurons using cell-specific deletion and genetic tracing in adult animals. Acute loss of Tcf4 in adult excitatory neurons in vivo caused hyperexcitability and increased dendritic complexity of neurons, effects that were distinct from previously observed effects in embryonic-deficiency models. Interestingly, transcriptomic analysis of genetically traced adult-deleted FACS-sorted Tcf4-knockout neurons revealed that Tcf4 targets in adult neurons are distinct from those in the embryonic brain. Meta-analysis of the adult-deleted neuronal transcriptome from our study with the existing datasets of embryonic Tcf4 deficiencies revealed plasma membrane and ciliary genes to underlie Tcf4-mediated structure-function regulation specifically in adult neurons. The profound changes both in the structure and excitability of adult neurons upon acute loss of Tcf4 indicates that proactive regulation of membrane-related processes underlies the functional and structural integrity of adult neurons. These findings not only provide insights for the functional relevance of continual expression of a psychiatric disease-risk gene in the adult brain but also identify previously unappreciated gene networks underpinning mature neuronal regulation during the adult lifespan.

Precise excitation-inhibition balance controls gain and timing in the hippocampus.

Excitation-inhibition (EI) balance controls excitability, dynamic range, and input gating in many brain circuits. Subsets of synaptic input can be selected or 'gated' by precise modulation of finely tuned EI balance, but assessing the granularity of EI balance requires combinatorial analysis of excitatory and inhibitory inputs. Using patterned optogenetic stimulation of mouse hippocampal CA3 neurons, we show that hundreds of unique CA3 input combinations recruit excitation and inhibition with a nearly identical ratio, demonstrating precise EI balance at the hippocampus. Crucially, the delay between excitation and inhibition decreases as excitatory input increases from a few synapses to tens of synapses. This creates a dynamic millisecond-range window for postsynaptic excitation, controlling membrane depolarization amplitude and timing via subthreshold divisive normalization. We suggest that this combination of precise EI balance and dynamic EI delays forms a general mechanism for millisecond-range input gating and subthreshold gain control in feedforward networks.

Subunit exchange enhances information retention by CaMKII in dendritic spines.

Molecular bistables are strong candidates for long-term information storage, for example, in synaptic plasticity. Calcium/calmodulin-dependent protein Kinase II (CaMKII) is a highly expressed synaptic protein which has been proposed to form a molecular bistable switch capable of maintaining its state for years despite protein turnover and stochastic noise. It has recently been shown that CaMKII holoenzymes exchange subunits among themselves. Here, we used computational methods to analyze the effect of subunit exchange on the CaMKII pathway in the presence of diffusion in two different micro-environments, the post synaptic density (PSD) and spine cytosol. We show that CaMKII exhibits multiple timescales of activity due to subunit exchange. Further, subunit exchange enhances information retention by CaMKII both by improving the stability of its switching in the PSD, and by slowing the decay of its activity in the spine cytosol. The existence of diverse timescales in the synapse has important theoretical implications for memory storage in networks.

Synaptic input sequence discrimination on behavioral timescales mediated by reaction-diffusion chemistry in dendrites.

Sequences of events are ubiquitous in sensory, motor, and cognitive function. Key computational operations, including pattern recognition, event prediction, and plasticity, involve neural discrimination of spatio-temporal sequences. Here, we show that synaptically-driven reaction-diffusion pathways on dendrites can perform sequence discrimination on behaviorally relevant time-scales. We used abstract signaling models to show that selectivity arises when inputs at successive locations are aligned with, and amplified by, propagating chemical waves triggered by previous inputs. We incorporated biological detail using sequential synaptic input onto spines in morphologically, electrically, and chemically detailed pyramidal neuronal models based on rat data. Again, sequences were recognized, and local channel modulation downstream of putative sequence-triggered signaling could elicit changes in neuronal firing. We predict that dendritic sequence-recognition zones occupy 5 to 30 microns and recognize time-intervals of 0.2 to 5 s. We suggest that this mechanism provides highly parallel and selective neural computation in a functionally important time range.

Robust and Rapid Air-Borne Odor Tracking without Casting.

Casting behavior (zigzagging across an odor stream) is common in air/liquid-borne odor tracking in open fields; however, terrestrial odor localization often involves path selection in a familiar environment. To study this, we trained rats to run toward an odor source in a multi-choice olfactory arena with near-laminar airflow. We find that rather than casting, rats run directly toward an odor port, and if this is incorrect, they serially sample other sources. This behavior is consistent and accurate in the presence of perturbations, such as novel odors, background odor, unilateral nostril stitching, and turbulence. We developed a model that predicts that this run-and-scan tracking of air-borne odors is faster than casting, provided there are a small number of targets at known locations. Thus, the combination of best-guess target selection with fallback serial sampling provides a rapid and robust strategy for finding odor sources in familiar surroundings.

Rats track odour trails accurately using a multi-layered strategy with near-optimal sampling.

Tracking odour trails is a crucial behaviour for many animals, often leading to food, mates or away from danger. It is an excellent example of active sampling, where the animal itself controls how to sense the environment. Here we show that rats can track odour trails accurately with near-optimal sampling. We trained rats to follow odour trails drawn on paper spooled through a treadmill. By recording local field potentials (LFPs) from the olfactory bulb, and sniffing rates, we find that sniffing but not LFPs differ between tracking and non-tracking conditions. Rats can track odours within ~1 cm, and this accuracy is degraded when one nostril is closed. Moreover, they show path prediction on encountering a fork, wide 'casting' sweeps on encountering a gap and detection of reappearance of the trail in 1-2 sniffs. We suggest that rats use a multi-layered strategy, and achieve efficient sampling and high accuracy in this complex task.

Trafficking motifs as the basis for two-compartment signaling systems to form multiple stable states.

Transport of molecules in cells is a central part of cell biology. Frequently such trafficking is not just for material transport, but also for information propagation, and serves to couple signaling circuits across cellular compartments. Here, I show that trafficking transforms simple local signaling pathways into self-organizing systems that span compartments and confer distinct states and identities to these compartments. I find that three motifs encapsulate the responses of most single-compartment signaling pathways in the context of trafficking. These motifs combine with different trafficking reactions to generate a diverse set of cellular functions. For example, trafficked bistable switches can oscillate or become quad- or tristable, depending on trafficking mechanisms and rates. Furthermore, the analysis shows how compartments participating in traffic can settle to distinct molecular compositions characteristic of distinct organelle identities. This general framework shows how the interplay between molecular movement and local reactions can generate many system functions, and give distinct identities to different parts of the cell.

Odor representations in the mammalian olfactory bulb.

A first key step in studying a sensory modality is to define how the brain represents the features of the sensory stimulus. This has proven to be a challenge in olfaction, where even the stimulus features have been a matter of considerable debate. In this review, we focus on olfactory representations in the first stage of the olfactory pathway, the olfactory bulb (OB). We examine the diverging viewpoints on spatially organized versus distributed representations. We then consider how odor sampling through respiration is a key part of the odorant code. Finally, we ask how the bulb handles the challenging task of representing mixtures. We suggest that current evidence points toward a representation that is spatially organized at the inputs but later distributed, with the spatial organization not being used for much computation. Nevertheless, this is a simple representation that effectively represents multiple individual odorants, as well as odor mixtures.

Developing complex signaling models using GENESIS/Kinetikit.

The development of biologically realistic models of signaling pathways is a demanding process, involving computational challenges as well as those arising from the complexity of detailed pathway models. We have developed the General Neural Simulation System (GENESIS) and Kinetikit (GENESIS/Kinetikit), a graphical simulation environment for modeling biochemical signaling pathways using deterministic and stochastic methods. A library of models of several common signaling pathways complements the software. This combination of numerical computation engines, graphical modeling tools, and library of models is designed to build on the cumulative development of models and techniques from many sources. The complete simulation environment and demonstration models are available from (http://stke.sciencemag.org/cgi/content/full/sigtrans;2004/219/pl4/DC1; also at http://www.ncbs.res.in/~bhalla/kkit/download.html). The associated library of signaling pathways is based on published experimental and simulation studies and is curated to ensure that the simulation outcomes match published results. Models in the library are maintained in a database (http://doqcs.ncbs.res.in). Individual pathway models can be combined to build complex signaling network simulations. The overall goal of this process is to attain sufficient biological realism in models to directly compare their outcomes with experiments and to improve our understanding of complex signaling.