Intracellular signaling networks receive and process information to control cellular machines. The mitogen-activated protein kinase (MAPK) 1,2/protein kinase C (PKC) system is one such network that regulates many cellular machines, including the cell cycle machinery and autocrine/paracrine factor synthesizing machinery. We used a combination of computational analysis and experiments in mouse NIH-3T3 fibroblasts to understand the design principles of this controller network. We find that the growth factor-stimulated signaling network containing MAPK 1, 2/PKC can operate with one (monostable) or two (bistable) stable states. At low concentrations of MAPK phosphatase, the system exhibits bistable behavior, such that brief stimulus results in sustained MAPK activation. The MAPK-induced increase in the amounts of MAPK phosphatase eliminates the prolonged response capability and moves the network to a monostable state, in which it behaves as a proportional response system responding acutely to stimulus. Thus, the MAPK 1, 2/PKC controller network is flexibly designed, and MAPK phosphatase may be critical for this flexible response.

Many distinct signaling pathways allow the cell to receive, process, and respond to information. Often, components of different pathways interact, resulting in signaling networks. Biochemical signaling networks were constructed with experimentally obtained constants and analyzed by computational methods to understand their role in complex biological processes. These networks exhibit emergent properties such as integration of signals across multiple time scales, generation of distinct outputs depending on input strength and duration, and self-sustaining feedback loops. Feedback can result in bistable behavior with discrete steady-state activities, well-defined input thresholds for transition between states and prolonged signal output, and signal modulation in response to transient stimuli. These properties of signaling networks raise the possibility that information for "learned behavior" of biological systems may be stored within intracellular biochemical reactions that comprise signaling pathways.

Most of the published quantitative models in biology are lost for the community because they are either not made available or they are insufficiently characterized to allow them to be reused. The lack of a standard description format, lack of stringent reviewing and authors' carelessness are the main causes for incomplete model descriptions. With today's increased interest in detailed biochemical models, it is necessary to define a minimum quality standard for the encoding of those models. We propose a set of rules for curating quantitative models of biological systems. These rules define procedures for encoding and annotating models represented in machine-readable form. We believe their application will enable users to (i) have confidence that curated models are an accurate reflection of their associated reference descriptions, (ii) search collections of curated models with precision, (iii) quickly identify the biological phenomena that a given curated model or model constituent represents and (iv) facilitate model reuse and composition into large subcellular models.During the genomic era we have witnessed a vast increase in availability of large amounts of quantitative data. This is motivating a shift in the focus of molecular and cellular research from qualitative descriptions of biochemical interactions towards the quantification of such interactions and their dynamics. One of the tenets of systems biology is the use of quantitative models (see Box 1 for definitions) as a mechanism for capturing precise hypotheses and making predictions 1,2 . Many specialized models exist that attempt to explain aspects of the cellular machinery. However, as has happened with other types of biological information, such as sequences, macromolecular structures or Box 1 GlossarySome terms are used in a very specific way throughout the article. We provide here a precise definition of each one.Quantitative biochemical model. A formal model of a biological system, based on the mathematical description of its molecular and cellular components, and the interactions between those components.Encoded model. A mathematical model written in a formal machine-readable language, such that it can be systematically parsed and employed by simulation and analysis software without further human translation. MIRIAM-compliant model. A model that passes all the tests and fulfills all the conditions listed in MIRIAM.Reference description. A unique document that describes, or references the description of the model, the structure of the model, the numerical values necessary to instantiate a simulation from the model, or to perform a mathematical analysis of the model, and the results one expects from such a simulation or analysis.Curation process. The process by which the compliance of an encoded model with MIRIAM is achieved and/or verified. The curation process may encompass some or all of the following tasks: encoding of the model, verification of the reference correspondence and annotation of the model.Reference correspondence. The fact that the...

Biological signaling pathways interact with one another to form complex networks. Complexity arises from the large number of components, many with isoforms that have partially overlapping functions; from the connections among components; and from the spatial relationship between components. The origins of the complex behavior of signaling networks and analytical approaches to deal with the emergent complexity are discussed here.Signaling in biological systems occurs at multiple levels. In its broad sense, one could use the term "signaling" to describe events ranging from interactions between single molecules to interactions between species in ecological systems. The aim here is to deal with complexity in signaling at a single level: intracellular signaling within a cell. We will outline how current and forthcoming tools in biochemistry, cell and molecular biology, and physiology, as well as theoretical analysis and simulation methods, may be used to study this complex system. In a general sense, the adjective "complex" describes a system or component that by design or function or both is difficult to understand and verify. In the past decade, analysis of complex systems (the field of complexity) has emerged as a distinct facet of mathematical and physical sciences. Understanding of biological systems may be enhanced by analysis of their complex nature. In physical systems, complexity is determined by such factors as the number of components and the intricacy of the interfaces between them, the number and intricacy of conditional branches, the degree of nesting, and the types of data structures. Biological signaling networks possess many of these attributes, as well as dynamic assembly, translocation, degradation, and channeling of chemical reactions. All of these activities occur simultaneously, and each component participates in several different activities.One approach to understanding complexity is to start with a conceptually simple view of signaling and add details that introduce new levels of complexity. As this process unfolds, it becomes clear where experimental data end and how progressively more difficult it becomes to understand the system as a whole in terms of the functional details of individual components. A Signaling WireThe simplest description of signaling may be in terms of elementary chemistry in a homogenous well-stirred cell where all molecules have equal access to each other. Here, the most upstream component of the signaling pathway interacts with an external source and

Biologically detailed single neuron and network models are important for understanding how ion channels, synapses and anatomical connectivity underlie the complex electrical behavior of the brain. While neuronal simulators such as NEURON, GENESIS, MOOSE, NEST, and PSICS facilitate the development of these data-driven neuronal models, the specialized languages they employ are generally not interoperable, limiting model accessibility and preventing reuse of model components and cross-simulator validation. To overcome these problems we have used an Open Source software approach to develop NeuroML, a neuronal model description language based on XML (Extensible Markup Language). This enables these detailed models and their components to be defined in a standalone form, allowing them to be used across multiple simulators and archived in a standardized format. Here we describe the structure of NeuroML and demonstrate its scope by converting into NeuroML models of a number of different voltage- and ligand-gated conductances, models of electrical coupling, synaptic transmission and short-term plasticity, together with morphologically detailed models of individual neurons. We have also used these NeuroML-based components to develop an highly detailed cortical network model. NeuroML-based model descriptions were validated by demonstrating similar model behavior across five independently developed simulators. Although our results confirm that simulations run on different simulators converge, they reveal limits to model interoperability, by showing that for some models convergence only occurs at high levels of spatial and temporal discretisation, when the computational overhead is high. Our development of NeuroML as a common description language for biophysically detailed neuronal and network models enables interoperability across multiple simulation environments, thereby improving model transparency, accessibility and reuse in computational neuroscience.

Sensory inputs frequently converge on the brain in a spatially organized manner, often with overlapping inputs to multiple target neurons. Whether the responses of target neurons with common inputs become decorrelated depends on the contribution of local circuit interactions. We addressed this issue in the olfactory system using newly generated transgenic mice expressing channelrhodopsin-2 in all olfactory sensory neurons. By selectively stimulating individual glomeruli with light, we identified mitral/tufted (M/T) cells that receive common input (sister cells). Sister M/T cells had highly correlated responses to odors as measured by average spike rates, but their spike timing in relation to respiration was differentially altered. In contrast, non-sister M/T cells correlated poorly on both these measures. We suggest that sister M/T cells carry two different channels of information: average activity representing shared glomerular input, and phase-specific information that refines odor representations and is substantially independent for sister M/T cells.

It has been hypothesized that rats and other mammals can use stereo cues to localize odor sources, but there is limited behavioral evidence to support this hypothesis. We found that rats trained on an odor-localization task can localize odors accurately in one or two sniffs. Bilateral sampling was essential for accurate odor localization, with internasal intensity and timing differences as directional cues. If the stimulus arrived at the correct point of the respiration cycle, internasal timing differences as short as 50 milliseconds sufficed. Neuronal recordings show that bulbar neurons responded differentially to stimuli from the left and stimuli from the right.

1. Detailed compartmental computer simulations of single mitral and granule cells of the vertebrate olfactory bulb were constructed using previously published geometric data. Electrophysiological properties were determined by comparing model output to previously published experimental data, mainly current-clamp recordings. 2. The passive electrical properties of each model were explored by comparing model output with intracellular potential data from hyperpolarizing current injection experiments. The results suggest that membrane resistivity in both cells is nonuniform, with somatas having a substantially lower resistivity than the dendrites. 3. The active properties of these cells were explored by incorporating active ion channels into modeled compartments. On the basis of evidence from the literature, the mitral cell model included six channel types: fast sodium, fast delayed rectifier (Kfast), slow delayed rectifier (K), transient outward potassium current (KA), voltage- and calcium-dependent potassium current (KCa), and L-type calcium current. The granule cell model included four channel types: rat brain sodium, K, KA, and the non-inactivating muscarinic potassium current (KM). Modeled channels were based on the Hodgkin-Huxley formalism. 4. Representative kinetics for each of the channel classes above were obtained from the literature. The experimentally unknown spatial distributions of each included channel were obtained by systematic parameter searches. These were conducted in two ways: large-scale simulation series, in which each parameter was varied in turn, and an adaptation of a multidimensional conjugate gradient method. In each case, the simulated results were compared wtih experimental data using a curve-matching function evaluating mean squared differences of several aspects of the simulated and experimental voltage waveforms. 5. Systematic parameter variations revealed a single distinct region of parameter space in which the mitral cell model best fit the data. This region of parameter space was also very robust to parameter variations. Specifically, optimum performance was obtained when calcium and slow K channels were concentrated in the glomeruli, with a lower density in the soma and proximal secondary dendrites. The distribution of sodium and fast potassium channels, on the other hand, was highest at the soma and axon, with a much lighter distribution throughout the secondary dendrites. The KA and KCa channels were also concentrated near the soma. 6. The parameter search of the granule cell model was much less restrained by experimental data. Several parameter regimes were found that gave a good match to the data.(ABSTRACT TRUNCATED AT 400 WORDS)