Pulsatile signaling in intercellular communication. Periodic stimuli are more efficient than random or chaotic signals in a model based on receptor desensitization

Assessing Different Temporal Scales of Calcium Dynamics in Networks of Beta Cell Populations

Beta cells within the pancreatic islets of Langerhans respond to stimulation with coherent oscillations of membrane potential and intracellular calcium concentration that presumably drive the pulsatile exocytosis of insulin. Their rhythmic activity is multimodal, resulting from networked feedback interactions of various oscillatory subsystems, such as the glycolytic, mitochondrial, and electrical/calcium components. How these oscillatory modules interact and affect the collective cellular activity, which is a prerequisite for proper hormone release, is incompletely understood. In the present work, we combined advanced confocal Ca2+ imaging in fresh mouse pancreas tissue slices with time series analysis and network science approaches to unveil the glucose-dependent characteristics of different oscillatory components on both the intra- and inter-cellular level. Our results reveal an interrelationship between the metabolically driven low-frequency component and the electrically driven high-frequency component, with the latter exhibiting the highest bursting rates around the peaks of the slow component and the lowest around the nadirs. Moreover, the activity, as well as the average synchronicity of the fast component, considerably increased with increasing stimulatory glucose concentration, whereas the stimulation level did not affect any of these parameters in the slow component domain. Remarkably, in both dynamical components, the average correlation decreased similarly with intercellular distance, which implies that intercellular communication affects the synchronicity of both types of oscillations. To explore the intra-islet synchronization patterns in more detail, we constructed functional connectivity maps. The subsequent comparison of network characteristics of different oscillatory components showed more locally clustered and segregated networks of fast oscillatory activity, while the slow oscillations were more global, resulting in several long-range connections and a more cohesive structure. Besides the structural differences, we found a relatively weak relationship between the fast and slow network layer, which suggests that different synchronization mechanisms shape the collective cellular activity in islets, a finding which has to be kept in mind in future studies employing different oscillations for constructing networks.

Tripartite cell networks for glucose homeostasis

Controlling the excess and shortage of energy is a fundamental task for living organisms. Diabetes is a representative metabolic disease caused by the malfunction of energy homeostasis. The islets of Langerhans in the pancreas release long-range messengers, hormones, into the blood to regulate the homeostasis of the primary energy fuel, glucose. The hormone and glucose levels in the blood show rhythmic oscillations with a characteristic period of 5-10 min, and the functional roles of the oscillations are not clear. Each islet has [Formula: see text] and [Formula: see text] cells that secrete glucagon and insulin, respectively. These two counter-regulatory hormones appear sufficient to increase and decrease glucose levels. However, pancreatic islets have a third cell type, [Formula: see text] cells, which secrete somatostatin. The three cell populations have a unique spatial organization in islets, and they interact to perturb their hormone secretions. The mini-organs of islets are scattered throughout the exocrine pancreas. Considering that the human pancreas contains approximately a million islets, the coordination of hormone secretion from the multiple sources of islets and cells within the islets should have a significant effect on human physiology. In this review, we introduce the hierarchical organization of tripartite cell networks, and recent biophysical modeling to systematically understand the oscillations and interactions of [Formula: see text], [Formula: see text], and [Formula: see text] cells. Furthermore, we discuss the functional roles and clinical implications of hormonal oscillations and their phase coordination for the diagnosis of type II diabetes.

Insulin modulates the frequency of Ca2+ oscillations in mouse pancreatic islets

Pancreatic islets can adapt to oscillatory glucose to produce synchronous insulin pulses. Can islets adapt to other oscillatory stimuli, specifically insulin? To answer this question, we stimulated islets with pulses of exogenous insulin and measured their Ca²⁺ oscillations. We observed that sufficiently high insulin (> 500 nM) with an optimal pulse period (~ 4 min) could make islets to produce synchronous Ca²⁺ oscillations. Glucose and insulin, which are key stimulatory factors of islets, modulate islet Ca²⁺ oscillations differently. Glucose increases the active-to-silent ratio of phases, whereas insulin increases the period of the oscillation. To examine the dual modulation, we adopted a phase oscillator model that incorporated the phase and frequency modulations. This mathematical model showed that out-of-phase oscillations of glucose and insulin were more effective at synchronizing islet Ca²⁺ oscillations than in-phase stimuli. This finding suggests that a phase shift in glucose and insulin oscillations can enhance inter-islet synchronization.

A Local Counter-Regulatory Motif Modulates the Global Phase of Hormonal Oscillations

Counter-regulatory elements maintain dynamic equilibrium ubiquitously in living systems. The most prominent example, which is critical to mammalian survival, is that of pancreatic α and β cells producing glucagon and insulin for glucose homeostasis. These cells are not found in a single gland but are dispersed in multiple micro-organs known as the islets of Langerhans. Within an islet, these two reciprocal cell types interact with each other and with an additional cell type: the δ cell. By testing all possible motifs governing the interactions of these three cell types, we found that a unique set of positive/negative intra-islet interactions between different islet cell types functions not only to reduce the superficially wasteful zero-sum action of glucagon and insulin but also to enhance/suppress the synchronization of hormone secretions between islets under high/normal glucose conditions. This anti-symmetric interaction motif confers effective controllability for network (de)synchronization.

Phase modulation of insulin pulses enhances glucose regulation and enables inter-islet synchronization

Insulin is secreted in a pulsatile manner from multiple micro-organs called the islets of Langerhans. The amplitude and phase (shape) of insulin secretion are modulated by numerous factors including glucose. The role of phase modulation in glucose homeostasis is not well understood compared to the obvious contribution of amplitude modulation. In the present study, we measured Ca²⁺ oscillations in islets as a proxy for insulin pulses, and we observed their frequency and shape changes under constant/alternating glucose stimuli. Here we asked how the phase modulation of insulin pulses contributes to glucose regulation. To directly answer this question, we developed a phenomenological oscillator model that drastically simplifies insulin secretion, but precisely incorporates the observed phase modulation of insulin pulses in response to glucose stimuli. Then, we mathematically modeled how insulin pulses regulate the glucose concentration in the body. The model of insulin oscillation and glucose regulation describes the glucose-insulin feedback loop. The data-based model demonstrates that the existence of phase modulation narrows the range within which the glucose concentration is maintained through the suppression/enhancement of insulin secretion in conjunction with the amplitude modulation of this secretion. The phase modulation is the response of islets to glucose perturbations. When multiple islets are exposed to the same glucose stimuli, they can be entrained to generate synchronous insulin pulses. Thus, we conclude that the phase modulation of insulin pulses is essential for glucose regulation and inter-islet synchronization.

On the translocation of bacteria and their lipopolysaccharides between blood and peripheral locations in chronic, inflammatory diseases: the central roles of LPS and LPS-induced cell death

We have recently highlighted (and added to) the considerable evidence that blood can contain dormant bacteria. By definition, such bacteria may be resuscitated (and thus proliferate). This may occur under conditions that lead to or exacerbate chronic, inflammatory diseases that are normally considered to lack a microbial component. Bacterial cell wall components, such as the endotoxin lipopolysaccharide (LPS) of Gram-negative strains, are well known as potent inflammatory agents, but should normally be cleared. Thus, their continuing production and replenishment from dormant bacterial reservoirs provides an easy explanation for the continuing, low-grade inflammation (and inflammatory cytokine production) that is characteristic of many such diseases. Although experimental conditions and determinants have varied considerably between investigators, we summarise the evidence that in a great many circumstances LPS can play a central role in all of these processes, including in particular cell death processes that permit translocation between the gut, blood and other tissues. Such localised cell death processes might also contribute strongly to the specific diseases of interest. The bacterial requirement for free iron explains the strong co-existence in these diseases of iron dysregulation, LPS production, and inflammation. Overall this analysis provides an integrative picture, with significant predictive power, that is able to link these processes via the centrality of a dormant blood microbiome that can resuscitate and shed cell wall components.

Band-pass processing in a GPCR signaling pathway selects for NFAT transcription factor activation

Many biological processes are rhythmic and proper timing is increasingly appreciated as being critical for development and maintenance of physiological functions. To understand how temporal modulation of an input signal influences downstream responses, we employ microfluidic pulsatile stimulation of a G-protein coupled receptor, the muscarinic M3 receptor, in single cells with simultaneous real-time imaging of both intracellular calcium and NFAT nuclear localization. Interestingly, we find that reduced stimulation with pulses of ligand can give more efficient transcription factor activation, if stimuli are timed appropriately. Our experiments and computational analyses show that M3 receptor-induced calcium oscillations form a low pass filter while calcium-induced NFAT translocation forms a high pass filter. The combination acts as a band-pass filter optimized for intermediate frequencies of stimulation. We demonstrate that receptor desensitization and NFAT translocation rates determine critical features of the band-pass filter and that the band-pass may be shifted for different receptors or NFAT dynamics. As an example, we show that the two NFAT isoforms (NFAT4 and NFAT1) have shifted band-pass windows for the same receptor. While we focus specifically on the M3 muscarinic receptor and NFAT translocation, band-pass processing is expected to be a general theme that applies to multiple signaling pathways.

Accurate encoding and decoding by single cells: amplitude versus frequency modulation

Cells sense external concentrations and, via biochemical signaling, respond by regulating the expression of target proteins. Both in signaling networks and gene regulation there are two main mechanisms by which the concentration can be encoded internally: amplitude modulation (AM), where the absolute concentration of an internal signaling molecule encodes the stimulus, and frequency modulation (FM), where the period between successive bursts represents the stimulus. Although both mechanisms have been observed in biological systems, the question of when it is beneficial for cells to use either AM or FM is largely unanswered. Here, we first consider a simple model for a single receptor (or ion channel), which can either signal continuously whenever a ligand is bound, or produce a burst in signaling molecule upon receptor binding. We find that bursty signaling is more accurate than continuous signaling only for sufficiently fast dynamics. This suggests that modulation based on bursts may be more common in signaling networks than in gene regulation. We then extend our model to multiple receptors, where continuous and bursty signaling are equivalent to AM and FM respectively, finding that AM is always more accurate. This implies that the reason some cells use FM is related to factors other than accuracy, such as the ability to coordinate expression of multiple genes or to implement threshold crossing mechanisms.

Signals, and hence information, can generally be transmitted either by amplitude (AM) or frequency (FM) modulation, as used, for example, in the transmission of radio waves since the 1930s. Both types of modulation are known to play a role in biology with AM conventionally associated with signaling and gene expression, and FM used to reliably transmit electrical signals over large distances between neurons. Surprisingly, FM was recently also observed in gene regulation, making their roles less distinct than previously thought. Although the engineering advantages and disadvantages of AM and FM are well understood, the equivalent question in biological systems is still largely unsolved. Here, we propose a simple model of signaling by receptors (or ion channels) with subsequent gene regulation, thus implementing both AM and FM in different types of biological pathways. We then compare the accuracy in the production of target proteins. We find that FM can be more accurate than AM only for a single receptor with fast signaling, whereas AM is more accurate in slow gene regulation and with signaling by multiple receptors. Finally, we propose possible reasons that cells use FM despite the potential decrease in accuracy.

Leaps and lulls in the developmental transcriptome of Dictyostelium discoideum

BACKGROUND

Development of the soil amoeba Dictyostelium discoideum is triggered by starvation. When placed on a solid substrate, the starving solitary amoebae cease growth, communicate via extracellular cAMP, aggregate by tens of thousands and develop into multicellular organisms. Early phases of the developmental program are often studied in cells starved in suspension while cAMP is provided exogenously. Previous studies revealed massive shifts in the transcriptome under both developmental conditions and a close relationship between gene expression and morphogenesis, but were limited by the sampling frequency and the resolution of the methods.

RESULTS

Here, we combine the superior depth and specificity of RNA-seq-based analysis of mRNA abundance with high frequency sampling during filter development and cAMP pulsing in suspension. We found that the developmental transcriptome exhibits mostly gradual changes interspersed by a few instances of large shifts. For each time point we treated the entire transcriptome as single phenotype, and were able to characterize development as groups of similar time points separated by gaps. The grouped time points represented gradual changes in mRNA abundance, or molecular phenotype, and the gaps represented times during which many genes are differentially expressed rapidly, and thus the phenotype changes dramatically. Comparing developmental experiments revealed that gene expression in filter developed cells lagged behind those treated with exogenous cAMP in suspension. The high sampling frequency revealed many genes whose regulation is reproducibly more complex than indicated by previous studies. Gene Ontology enrichment analysis suggested that the transition to multicellularity coincided with rapid accumulation of transcripts associated with DNA processes and mitosis. Later development included the up-regulation of organic signaling molecules and co-factor biosynthesis. Our analysis also demonstrated a high level of synchrony among the developing structures throughout development.

CONCLUSIONS

Our data describe D. discoideum development as a series of coordinated cellular and multicellular activities. Coordination occurred within fields of aggregating cells and among multicellular bodies, such as mounds or migratory slugs that experience both cell-cell contact and various soluble signaling regimes. These time courses, sampled at the highest temporal resolution to date in this system, provide a comprehensive resource for studies of developmental gene expression.

Motivations and methods for analyzing pulsatile hormone secretion

Endocrine glands communicate with remote target cells via a mixture of continuous and intermittent signal exchange. Continuous signaling allows slowly varying control, whereas intermittency permits large rapid adjustments. The control systems that mediate such homeostatic corrections operate in a species-, gender-, age-, and context-selective fashion. Significant progress has been made in understanding mechanisms of adaptive interglandular signaling in vivo. Principal goals are to understand the physiological origins, significance, and mechanisms of pulsatile hormone secretion. Key analytical issues are: 1) to quantify the number, size, shape, and uniformity of pulses, nonpulsatile (basal) secretion, and elimination kinetics; 2) to evaluate regulation of the axis as a whole; and 3) to reconstruct dose-response interactions without disrupting hormone connections. This review will focus on the motivations driving and the methodologies used for such analyses.

Dual chemotaxis signalling regulates Dictyostelium development: intercellular cyclic AMP pulses and intracellular F-actin disassembly waves induce each other

Aggregating Dictyostelium discoideum amoebae periodically emit and relay cAMP, which regulates their chemotaxis and morphogenesis into a multicellular, differentiated organism. Cyclic AMP also stimulates F-actin assembly and chemotactic pseudopodium extension. We used actin-GFP expression to visualise for the first time intracellular F-actin assembly as a spatio-temporal indicator of cell reactions to cAMP, and thus the kinematics of cell communication, in aggregating streams. Every natural cAMP signal pulse induces an autowave of F-actin disassembly, which propagates from each cell's leading end to its trailing end at a linear rate, much slower than the calculated and measured velocities of cAMP diffusion in aggregating Dictyostelium. A sequence of transient reactions follows behind the wave, including anterior F-actin assembly, chemotactic pseudopodium extension and cell advance at the cell front and, at the back, F-actin assembly, extension of a small retrograde pseudopodium (forcing a brief cell retreat) and chemotactic stimulation of the following cell, yielding a 20s cAMP relay delay. These dynamics indicate that stream cell behaviour is mediated by a dual signalling system: a short-range cAMP pulse directed from one cell tail to an immediately following cell front and a slower, long-range wave of intracellular F-actin disassembly, each inducing the other.

Decoding of calcium oscillations by phosphorylation cycles: analytic results

Experimental studies have demonstrated that Ca(2+)-regulated proteins are sensitive to the frequency of Ca(2+) oscillations, and several mathematical models for specific proteins have provided insight into the mechanisms involved. Because of the large number of Ca(2+)-regulated proteins in signal transduction, metabolism and gene expression, it is desirable to establish in general terms which molecular properties shape the response to oscillatory Ca(2+) signals. Here we address this question by analyzing in detail a model of a prototypical Ca(2+)-decoding module, consisting of a target protein whose activity is controlled by a Ca(2+)-activated kinase and the counteracting phosphatase. We show that this module can decode the frequency of Ca(2+) oscillations, at constant average Ca(2+) signal, provided that the Ca(2+) spikes are narrow and the oscillation frequency is sufficiently low--of the order of the phosphatase rate constant or below. Moreover, Ca(2+) oscillations activate the target more efficiently than a constant signal when Ca(2+) is bound cooperatively and with low affinity. Thus, the rate constants and the Ca(2+) affinities of the target-modifying enzymes can be tuned in such a way that the module responds optimally to Ca(2+) spikes of a certain amplitude and frequency. Frequency sensitivity is further enhanced when the limited duration of the external stimulus driving Ca(2+) signaling is accounted for. Thus, our study identifies molecular parameters that may be involved in establishing the specificity of cellular responses downstream of Ca(2+) oscillations.

Dynamic model of hormonal systems coupled by negative feedback

Most hormone concentrations in the body are regulated by negative feedback mechanisms in which the production and release of hormones are regulated according to the concentration of related species. Also, it has been observed that several hormones are released in a variety of pulsatile patterns. In most cases, the mechanism driving these complex patterns is not well understood. Our model of two cells coupled through negative feedback to their external products demonstrates periodic, aperiodic and chaotic oscillations. The coupling between the cells seems to be responsible for these dynamic behaviors. The variety of dynamic behaviors observed in the model demonstrates that a simple physiological feedback loop mimicking the coupling between circulatory hormones and production centers could be the source of complex hormone release patterns observed in vivo.

Oscillations and waves of cyclic AMP in Dictyostelium: a prototype for spatio-temporal organization and pulsatile intercellular communication

The amoebae Dictyostelium discoideum aggregate after starvation in a wavelike manner in response to periodic pulses of cyclic AMP (cAMP) secreted by cells which behave as aggregation centers. In addition to autonomous oscillations, the cAMP signaling system that controls aggregation is also capable of excitable behavior, which consists in the transient amplification of suprathreshold pulses of extracellular cAMP. Since the first theoretical model for slime mold aggregation proposed by Keller and Segel in 1970, many theoretical studies have addressed various aspects of the mechanism and function of cAMP signaling in Dictyostelium. This paper presents a brief overview of these developments as well as some reminiscences of the author's collaboration with Lee Segel in modeling the dynamics of cAMP relay and oscillations. Considered in turn are models for cAMP signaling in Dictyostelium, the developmental path followed by the cAMP signaling system after starvation, the frequency encoding of cAMP signals, and the origin of concentric or spiral waves of cAMP.

Role of cascades in converting oscillatory signals into stationary step-like responses

In biological signal transduction pathways intermediates are often oscillatory and need to be converted into smooth output signals at the end. We show by mathematical modelling that protein kinase cascades enable converting oscillatory signals into sharp stationary step-like outputs. The importance of this result is demonstrated for the switch-like protein activation by calcium oscillations, which is of biological importance for regulating different cellular processes. In addition, we found that protein kinase cascades cause memory effects in the protein activation, which might be of a physiological advantage since a smaller amount of calcium transported in the cell is required for an effective activation of cellular processes.

A minimal model for decoding of time-limited Ca2+ oscillations

Calcium oscillations regulate several cellular processes by activating particular proteins. Most theoretical studies focused on the idealized situation of infinitely long oscillations. Here we analyze information transfer by time-limited calcium spike trains. We show that proteins can be selectively activated in a resonance-like manner by time-limited spike trains of different frequencies, while infinitely long oscillations do not show this resonance phenomenon. We found that proteins are activated more specifically by shorter oscillatory signals with narrower spikes.

Selective regulation of cellular processes via protein cascades acting as band-pass filters for time-limited oscillations

We show by mathematical modelling that a two-level protein cascade can act as a band-pass filter for time-limited oscillations. The band-pass filters are then combined into a network of three-level signalling cascades that by filtering the frequency of time-limited oscillations selectively switches cellular processes on and off. The physiological relevance for the selective regulation of cellular processes is demonstrated for the case of regulation by time-limited calcium oscillations.

Differential regulation of proteins by bursting calcium oscillations--a theoretical study

Calcium in ionic form is a second messenger connecting several input signals to several target processes in the cell. The question arises how one second messenger can transmit more than one signal simultaneously (bow-tie structure of signalling). Experimental data on calcium dynamics often show patterns of successive low-peak and high-peak oscillatory phases, known as bursting. Here, we propose that bursting calcium oscillations can perform the function of simultaneous transmission of two signals at physiological calcium concentrations, for example, by selective activation of two calcium-binding proteins. This differential regulation by periodic bursting is investigated in a theoretical model. The two proteins are assumed to be activated by calcium, and one of them is assumed to be subject to biphasic regulation due to additional inhibitory binding sites. To explore which characteristics of the complex signal could be responsible for independent regulation of low-peak activated and spike activated targets, different bursting patterns of simplified square pulses are applied. Depending on the change in the bursting pattern, one protein can be gradually activated at a constant level of the other protein's activity, or the two proteins can be activated simultaneously, or one protein can be activated while the other one is deactivated simultaneously. Thus, the two proteins can be regulated virtually independently.

Estimation of parameters for a mathematical model of growth hormone secretion

Here, we describe partial calibration of a parsimonious mathematical model of growth hormone (GH) secretion. From first principles, we derived a model of the effects on GH secretion from pituitary somatotrophs of stimulation by GH-releasing factor (GRF) or GH secretagogue, and of inhibition by somatostatin. We obtained a concise model by collapsing the many processes of the signal transduction cascade into a single step broadly reflecting the initial binding of GRF to its receptors. In the model, GH secretion is proportional to the rate of binding of GRF to activatable receptors. Desensitization occurs because of reduction of free receptors/available effector units, and resensitization occurs as those lost are replaced. This replacement is speeded up in the presence of somatostatin, which also inhibits GH secretion by reducing the constant of proportionality between the rate of GH secretion and the rate of GRF binding. We derived simple mathematical equations for the rate of GH secretion and cumulative secretion. Using these, we tested the model against data obtained from experiments performed in vitro, and made it quantitative using rigorous statistical approaches to optimize parameter estimates. The behaviour of the calibrated model matches experimental observations closely.

Signal processing times in neutrophil activation: dependence on ligand concentration and the relative phase of metabolic oscillations

Intracellular NAD(P)H oscillations exhibited by polarized neutrophils display congruent with 20 s periods, which are halved to congruent with 10 s upon stimulation with chemotactic peptides such as FNLPNTL (N-formyl-nle-leu-phe-nle-tyr-lys). By monitoring this frequency change, we have measured accurately the time interval between stimulus and metabolic frequency changes. A microscope flow chamber was designed to allow rapid delivery of FNLPNTL to adherent cells. Using fluorescein as a marker, we found delivery to be complete and stable throughout the chamber within approximately 400 ms. Peptides were injected into the chamber at concentrations ranging from 10(-6) to 10(-9) M. Injections also varied with respect to the relative phase of a cell's NAD(P)H oscillations. The time interval between injection of 10(-6) M FNLPNTL and the acquisition of congruent with 10 s period metabolic oscillations was found to be 12.2+/-3.3 s when injections occurred at the NAD(P)H oscillation peak whereas the lag time was 22.5+/-4.8 s when coinciding with a trough. At 10(-8) M FNLPNTL, lag times were found to be 26.1+/-5.2 and 30.5+/-7.3 s for injections at NAD(P)H peaks and troughs, respectively. FNLPNTL at 10(-9) M had no effect on metabolic oscillations, consistent with previous studies. Our experiments show that the kinetics of transmembrane signal processing, in contrast to a simple transmembrane chemical reaction, can depend upon both ligand dose and its temporal relationship with intracellular metabolic oscillations.

Structure of the GnRH receptor-stimulated signaling network: insights from genomics

The GnRH receptor influences gene expression in the gonadotrope through activating signaling cascades that modulate transcription factor expression and activity. A longstanding question in neuroendocrinology is how instructions received at the membrane in the form of the pattern of receptor stimulation are processed into specific biosynthetic changes at each gonadotropin promoter. Signal transduction from the membrane to preformed transcription factors relies on recognition of altered conformations. Signal transduction through the layers of the gene network also requires the biosynthesis of new transcription factors. The signal processing of this system depends on its molecular connectivity map and its feedback and feed-forward loops. Review of signal transduction, gene control, and genomic studies provide evidence of key loops that cross between cellular and nuclear compartments. Genomic studies suggest that the signal transduction and gene network form a continuum. We propose that information transfer in the gonadotrope depends on robust signaling modules that serve to integrate events at different time scales across cytoplasmic and nuclear compartments.

Time-patterned drug administration: insights from a modeling approach

The physiological effects of a drug depend not only on its molecular structure but also on the time-pattern of its administration. One of the main reasons for the importance of temporal patterns in drug action is biological rhythms--particularly those of circadian period. These rhythms affect most physiological functions as well as drug metabolism, clearance, and dynamic processes that may alter drug availability and target cell responsiveness with reference to biological time. We present an overview of the importance of time-patterned signals in physiology focused on the insights provided by a modeling approach. We first discuss examples of pulsatile intercellular communication by hormones such as gonadotropin-releasing hormone, and by cyclic adenosine monophosphate (cAMP) signals in Dictyostelium amoebae. Models based on reversible receptor desensitization account in both cases for the existence of optimal patterns of pulsatile signaling. Turning to circadian rhythms, we examine how models can be used to account for the response of 24h patterns to external stimuli such as light pulses or gene expression, and to predict how to restore the physiological characteristics of altered rhythms. Time-patterned treatments of cancer involve two distinct lines of research. The first, currently evaluated in clinical trials, relies on circadian chronomodulation of anticancer drugs, while the second, mostly based on theoretical studies, involves a resonance phenomenon with the cell-cycle length. We discuss the implications of modeling studies to improve the temporal patterning of drug administration.

Nonequilibrium / nonlinear chemical oscillation in the virtual absence of gravity

The Belousov-Zhabotinsky (BZ) reactions were used as typical examples of a nonlinear system far from equilibrium in connection with biological evolution. The virtual absence of gravity in the present work was given from the free-fall facility of Japan Microgravity Center (JAMIC) in Hokkaido. The reaction solution of BZ reaction was composed of bromate in sulfuric acid, 1,4-cyclohexanedione and ferroin to visualize the time development of patterns of chemical oscillations in the reaction-diffusion system. It is a bubble-free constitution in the aging of the reaction. Therefore, the setup constructed to collect image data via CCD cameras was simplified. The operation sequences of necessary devices were comprised of simple solid state relays which were started by a command from the operation room of JAMIC. The propagation profile of chemical patterns under microgravity of 10(-5) g was collected as image data for 9.8 s, and processed by a software of STM-STS2. In the aqueous solutions, propagation velocity of chemical patterns under microgravity was decreased to 80.9 % of that under normal gravity, owing to suppression of convection. On the other hand, in gel matrix, gravity did not influence the propagation velocity.

Simulation of prokaryotic genetic circuits

Biochemical and genetic approaches have identified the molecular mechanisms of many genetic reactions, particularly in bacteria. Now a comparably detailed understanding is needed of how groupings of genes and related protein reactions interact to orchestrate cellular functions over the cell cycle, to implement preprogrammed cellular development, or to dynamically change a cell's processes and structures in response to environmental signals. Simulations using realistic, molecular-level models of genetic mechanisms and of signal transduction networks are needed to analyze dynamic behavior of multigene systems, to predict behavior of mutant circuits, and to identify the design principles applicable to design of genetic regulatory circuits. When the underlying design rules for regulatory circuits are understood, it will be far easier to recognize common circuit motifs, to identify functions of individual proteins in regulation, and to redesign circuits for altered functions.

Cellular memory: neutrophil orientation reverses during temporally decreasing chemoattractant concentrations

Cell directional orientation or shape polarization is the first cellular step in neutrophil locomotion. To better understand how chemoattractants interact with cells, we studied neutrophil polarization (or shape changes) during exposure to a temporally decreasing chemoattractant signal of N-formyl-methionyl-leucyl-phenylalanine (FMLP) in the absence of a spatial concentration gradient. To accomplish this objective, we used a manifold of differing FMLP concentrations attached to a stopped-flow microscope chamber. Spatial gradients of a fluorescent chemotactic peptide could not be detected in the chamber by using microfluorometry. When FMLP was injected at continually increasing concentrations at 10-s intervals, the shape and relative direction of the neutrophil persisted. However, when temporally decreasing FMLP concentrations were injected, approximately 80% of the cells changed their direction with 44% of the total cells swinging about to 180 degrees +/- 15 degrees. Most of these directional changes involved dissolution of both the lamellipodium and uropod and reformation of these structures 180 degrees from their original positions. This research suggests that neutrophils reverse their morphological polarity when exposed to temporally decreasing ligand concentrations by "remembering" their ligand exposure history and relative direction.

The organization and function of endogenous antinociceptive systems

Much progress has been made the understanding of endogenous pain-controlling systems. Recently, new concepts and ideas which are derived from neurobiology, chaos research and from research on learning and memory have been introduced into pain research and shed further light on the organization and function of endogenous antinociception. These most recent developments will be reviewed here. Three principles of endogenous antinociception have been identified, as follows. (1) Supraspinal descending inhibition: the patterns of neuronal activity in diencephalon, brainstem and spinal cord during antinociceptive stimulation in midbrain periaqueductal gray (PAG) or medullary nucleus raphe magnus have now been mapped on the cellular level, using the c-Fos technique. Results demonstrate that characteristic activity patterns result within and outside the PAG when stimulating at its various subdivisions. The descending systems may not only depress mean discharge rates of nociceptive spinal dorsal horn neurons, but also may modify harmonic oscillations and nonlinear dynamics (dimensionality) of discharges. (2) Propriospinal, heterosegmental inhibition: antinociceptive, heterosegmental interneurons exist which may be activated by noxious stimulation or by supraspinal descending pathways. (3) Segmental spinal inhibition: a robust long-term depression of primary afferent neurotransmission in A delta fibers has been identified in superficial spinal dorsal horn which may underlie long-lasting antinociception by afferent stimulation, e.g. by physical therapy or acupuncture.

Low dimensional attractors in discharges of sensory neurons of the rat spinal dorsal horn are maintained by supraspinal descending systems

Background activity was recorded from sensory neurons in laminae I-V of the lumbar spinal dorsal horn of the rat prior to and during cold block spinalization at the cervical cord. To graphically and quantitatively describe the complexity of the discharge patterns, phase space portraits were plotted and the correlation dimension D2 was calculated by using the Grassberger-Procaccia algorithm adopted for point processes, i.e. for series of interspike intervals. The algorithm is validated both for the Baker transformation, which is a simple point process, and for the Lorenz model, whereby a transformation from continuous to point process variables is achieved. A method of surrogate data is provided in order to show the difference between original neuronal patterns and their surrogate stochastic data. Therefore, this method shows that neuronal discharge patterns cannot be fully described in terms of interspike interval histograms. However, in the intact cord most (73%) of the neurons displayed background activity with low (0.28-4.3) D2 values. During spinalization, D2 values significantly increased in 68% of the neurons showing previously low D2 values, irrespective of classification and laminar location of neurons, thus proving that tonic descending systems may maintain a high order in the discharge of sensory dorsal horn neurons.

Modulation of caffeine contractures in mammalian skeletal muscles by variation of extracellular potassium

Caffeine contractures were induced after K(+)-conditioning of skeletal muscles from pigs and mice. K(+)-conditioning is defined as the partial depolarization caused by increasing external potassium (K+0) with [K+]x[Cl-] constant. Conditioning depolarizations that rendered muscles refractory to brief electrical stimulation still enhanced the contracture tension elicited by subsequent direct caffeine stimulation of sarcoplasmic reticulum (SR) calcium release. The effects of K(+)-conditioning on caffeine-induced contractures of intact cell bundles reached a maximum at 15-30 mM K+0 and then progressively declined at higher [K+]0. Conditioning with 30 mM K+ for 5 min, which inactivates excitation-contraction (EC) coupling in response to action potentials, both increased the magnitude of caffeine contractures 2-10-fold and shifted the contracture threshold toward lower caffeine concentrations. Enhanced sensitivity to caffeine was inhibited by dantrolene (20 microM) and its watersoluble analogue azumolene (150 microM). These drugs decreased caffeine-induced contractures following depolarization with 4-15 mM K+ to 25-50% of control tension. The inorganic anion perchlorate (CIO-4), which like caffeine potentiates twitches, increased caffeine-induced contractures approximately twofold after K(+)-conditioning (> 4 mM). The results suggest that CIO-4 and dantrolene, in addition to caffeine, also influence SR calcium release either directly or by mechanism(s) subsequent to depolarization of the sarcolemma. Moreover, since CIO-4 is known to shift the voltage-dependence of intramembrane charge movement, CIO-4 may exert effects on the transverse-tubule voltage sensors as well as the SR.

The effects of extrasynaptic substance P on nociceptive neurons in laminae I and II in rat lumbar spinal dorsal horn

Inflammation of the skin induces release and extrasynaptic spread of neuropeptides such as substance P mainly in spinal laminae I and II and causes changes in discharge properties of nociceptive neurons in spinal dorsal horn. To evaluate the role of extrasynaptic substance P we have superfused the spinal cord at recording segment with artificial cerebrospinal fluid or with substance P. A total of 102 multireceptive neurons responding to both noxious and innocuous skin stimulation were recorded in laminae I or II of lumbar spinal dorsal horn in pentobarbital anaesthetized rats. During superfusion with substance P (10 or 100 microM) significant increases of background activities (from 2.2 +/- 0.6 to 8.4 +/- 3.2 imp./s, mean +/- S.E.M.), enlargement of cutaneous receptive fields (from 359.9 +/- 60.4 to 465.5 +/- 77.3 mm2) and enhanced responses to mechanical (from 89.1 +/- 22.7 to 147.0 +/- 27.5 imp./5 s) but not thermal noxious skin stimuli were observed in the 22 neurons tested. Noxious heat-evoked responses and C-fibre-evoked responses were changed in both directions. In 50 other neurons, the coefficients of dispersion of interspike intervals, which is an indicator of burst-like discharges, were significantly reduced (from 60.4 +/- 5.5 to 52.7 +/- 5.3) after application of substance P. Substance P induced oscillations in background activities in 13 of 40 non-rhythmic neurons and depressed oscillations in 2 of 11 neurons. Cross-correlations of discharges of pairs simultaneously recorded neurons were flat (n = 4), or had a central peak (n = 19) or a central trough (n = 2) and were not changed qualitatively by extrasynaptic substance P. Thus, extrasynaptic substance P can modify not only discharge patterns in the spinal dorsal horn.

Identification and characterization of rhythmic nociceptive and non-nociceptive spinal dorsal horn neurons in the rat

The properties of rhythmic low-threshold and multireceptive spinal dorsal horn neurons were determined. Multiple neuron recordings were made via a single electrode in the lumbar spinal cord of pentobarbital-anesthetized or decerebrate, unanesthetized, spinalized rats. The background activity of a total of 223 neurons was analysed: 21.0% of 176 fully characterized neurons were low threshold, 73.3% multireceptive and 5.7% nociceptive-specific neurons. Twenty of 100 neurons tested were driven by antidromic stimulation at the upper cervical cord. To identify and evaluate rhythmic harmonic oscillations in the discharges of spinal dorsal horn neurons during background activity and steady-state noxious heat-evoked responses, interspike interval, autocorrelation and autospectral analysis were performed. The background activity of 99 of the 223 neurons (44.4%) of our sample was rhythmic. The distribution of the fundamental spectral frequencies has a bimodal shape, the first band between 0.5 and 2 Hz and the second between 6 and 13 Hz. Low-threshold and multireceptive neurons had a similar incidence of rhythmicity (54.1 and 43.4%, respectively). Only one of 20 neurons with long ascending projections presented rhythmic background activity. Activation of heat-sensitive nociceptors within the cutaneous receptive fields of the neurons had a strong anti-rhythmic effect in nine of 15 (60%) neurons. No change was observed in the pattern of autospectra of non-rhythmic neurons or low-threshold neurons during noxious stimulation. Twenty-four of 37 (66.6%) rhythmic neurons retained their rhythmic background discharges during reversible cold-block spinalization at the upper thoracic cord. The incidence of neurons with burst-like discharges was highest among multireceptive neurons (98/129, 75.9%) and non-rhythmic neurons (89/124, 71.8%). Thus, rhythmicity exists in sensory neurons of the spinal dorsal horn probably generated within its local neuronal network and partially modulated by supraspinal descending systems. Rhythmicity is depressed by activity in primary afferent nociceptors. The role of rhythmicity for information transfer and neuronal plasticity is discussed.

The regulation of chemotaxis and chemokinesis in Dictyostelium amoebae by temporal signals and spatial gradients of cyclic AMP

The tactic and kinetic locomotion of Dictyostelium discoideum amoebae were examined in cyclic AMP (cAMP) spatial gradient and temporal signal fields. The distributions of migrating cells were examined within 150 microns-thick micropore filters after incubation with different cAMP concentrations, [cAMP], applied in three ways across the fields: as positively or negatively developing gradients, generated either by increasing or decreasing the [cAMP] on one side of the filter, respectively, or as static, linear gradients after negative development. Chemotaxis was only induced by oriented, temporally increasing [cAMP]. Pulses propagated by molecular diffusion or mechanical flow were equally effective. Negatively developing cAMP gradients had no initial effect on cell accumulation. However, if the subsequent static spatial gradient was maintained by an infusion system, some gradients also induced cell accumulation, whose degree and direction depended on the gradient [cAMP]. The basis of this new effect was examined by tracking individual cells by computer-assisted videomicroscopy during locomotion in different [cAMP]. Cells produced a triphasic [cAMP]-dependent response, with optimal cell motility induced by 10–30 nM. The results demonstrate that cell accumulation either up-field or down-field in spatial gradients is governed by the field locations of the attractant concentrations that induce the relative locomotory maxima and minima in the gradient field. Cells perceive the ambient [cAMP], but cannot read the spatial gradient orientation in static or yet steeper regions of developing gradients. Accumulation in static spatial gradients is a function of klino- and orthokinesis, but chemotaxis requires an oriented cAMP pulse or impulse.(ABSTRACT TRUNCATED AT 250 WORDS)

Immunofluorescence localization of the Na-Ca exchanger in heart cells

We investigated the localization of the Na-Ca exchanger in fixed, isolated heart cells from rat and guinea pig using immunocytochemical methods with epifluorescence and confocal microscopy. We found that the Na-Ca exchanger is distributed throughout all membranes in contact with the extracellular space, including the sarcolemma, the transverse tubules (T-tubules), and the intercalated disks. Microscopic nonuniformities in the fluorescent labeling appear to reflect varying views of the membranes containing Na-Ca exchanger protein. Confocal thin-section imaging reveals a regular grid of discrete foci of fluorescence, which represent Na-Ca exchanger in T-tubules viewed en face. These foci are 1.80 +/- 0.01 microns apart from sarcomere to sarcomere and are aligned with the Z-line. Along each Z-line, these foci are spaced at 1.22 +/- 0.11-microns intervals. Longitudinal sections of the sarcolemma-T-tubule junction show a comblike appearance, with T-tubules extending inward from the heavily labeled sarcolemma. Our finding that the Na-Ca exchanger is widely distributed over the cell surface may provide further insight into the role of Na-Ca exchange in the heart.