Transition from stochastic to deterministic behavior in calcium oscillations

Simulation and modeling is becoming more and more important when studying complex biochemical systems. Most often, ordinary differential equations are employed for this purpose. However, these are only applicable when the numbers of participating molecules in the biochemical systems are large enough to be treated as concentrations. For smaller systems, stochastic simulations on discrete particle basis are more accurate. Unfortunately, there are no general rules for determining which method should be employed for exactly which problem to get the most realistic result. Therefore, we study the transition from stochastic to deterministic behavior in a widely studied system, namely the signal transduction via calcium, especially calcium oscillations. We observe that the transition occurs within a range of particle numbers, which roughly corresponds to the number of receptors and channels in the cell, and depends heavily on the attractive properties of the phase space of the respective systems dynamics. We conclude that the attractive properties of a system, expressed, e.g., by the divergence of the system, are a good measure for determining which simulation algorithm is appropriate in terms of speed and realism

On the correct determination of rotational angles for twisted-profiled sweep objects

In the context of the paper of V. Akman and A. Arslan, "Sweeping with All Graphical Ingredients in a Topological Picturebook" [Computers & Graphics, Vol. 16(3), pp. 273-281, 1992], a new construction of the rotational matrix is presented. This fixes a bug discovered by the first author. © 1994

Dose-Response Aligned Circuits in Signaling Systems

Cells use biological signal transduction pathways to respond to environmental stimuli and the behavior of many cell types depends on precise sensing and transmission of external information. A notable property of signal transduction that was characterized in the Saccharomyces cerevisiae yeast cell and many mammalian cells is the alignment of dose-response curves. It was found that the dose response of the receptor matches closely the dose responses of the downstream. This dose-response alignment (DoRA) renders equal sensitivities and concordant responses in different parts of signaling system and guarantees a faithful information transmission. The experimental observations raise interesting questions about the nature of the information transmission through DoRA signaling networks and design principles of signaling systems with this function. Here, we performed an exhaustive computational analysis on network architectures that underlie the DoRA function in simple regulatory networks composed of two and three enzymes. The minimal circuits capable of DoRA were examined with Michaelis-Menten kinetics. Several motifs that are essential for the dynamical function of DoRA were identified. Systematic analysis of the topology space of robust DoRA circuits revealed that, rather than fine-tuning the network's parameters, the function is primarily realized by enzymatic regulations on the controlled node that are constrained in limiting regions of saturation or linearity

Mechanisms for the Intracellular Manipulation of Organelles by Conventional Electroporation

Conventional electroporation (EP) changes both the conductance and molecular permeability of the plasma membrane (PM) of cells and is a standard method for delivering both biologically active and probe molecules of a wide range of sizes into cells. However, the underlying mechanisms at the molecular and cellular levels remain controversial. Here we introduce a mathematical cell model that contains representative organelles (nucleus, endoplasmic reticulum, mitochondria) and includes a dynamic EP model, which describes formation, expansion, contraction, and destruction for the plasma and all organelle membranes. We show that conventional EP provides transient electrical pathways into the cell, sufficient to create significant intracellular fields. This emerging intracellular electrical field is a secondary effect due to EP and can cause transmembrane voltages at the organelles, which are large enough and long enough to gate organelle channels, and even sufficient, at some field strengths, for the poration of organelle membranes. This suggests an alternative to nanosecond pulsed electric fields for intracellular manipulations.National Science Foundation (U.S.) (NSF Graduate Research Fellowship)National Institutes of Health (U.S.) (grant No. R01-GM63857)Aegis Industries, Inc

Mathematical modeling of intracellular signaling pathways

Dynamic modeling and simulation of signal transduction pathways is an important topic in systems biology and is obtaining growing attention from researchers with experimental or theoretical background. Here we review attempts to analyze and model specific signaling systems. We review the structure of recurrent building blocks of signaling pathways and their integration into more comprehensive models, which enables the understanding of complex cellular processes. The variety of mechanisms found and modeling techniques used are illustrated with models of different signaling pathways. Focusing on the close interplay between experimental investigation of pathways and the mathematical representations of cellular dynamics, we discuss challenges and perspectives that emerge in studies of signaling systems

Characterization of Membrane Potential Dependency of Mitochondrial Ca2+ Uptake by an Improved Biophysical Model of Mitochondrial Ca2+ Uniporter

Mitochondrial Ca2+ uniporter is the primary influx pathway for Ca2+ into respiring mitochondria, and hence plays a key role in mitochondrial Ca2+ homeostasis. Though the mechanism of extra-matrix Ca2+ dependency of mitochondrial Ca2+ uptake has been well characterized both experimentally and mathematically, the mechanism of membrane potential (ΔΨ) dependency of mitochondrial Ca2+ uptake has not been completely characterized. In this paper, we perform a quantitative reevaluation of a previous biophysical model of mitochondrial Ca2+ uniporter that characterized the possible mechanism of ΔΨ dependency of mitochondrial Ca2+ uptake. Based on a model simulation analysis, we show that model predictions with a variant assumption (Case 2: external and internal Ca2+ binding constants for the uniporter are distinct), that provides the best possible description of the ΔΨ dependency, are highly sensitive to variation in matrix [Ca2+], indicating limitations in the variant assumption (Case 2) in providing physiologically plausible description of the observed ΔΨ dependency. This sensitivity is attributed to negative estimate of a biophysical parameter that characterizes binding of internal Ca2+ to the uniporter. Reparameterization of the model with additional nonnengativity constraints on the biophysical parameters showed that the two variant assumptions (Case 1 and Case 2) are indistinguishable, indicating that the external and internal Ca2+ binding constants for the uniporter may be equal (Case 1). The model predictions in this case are insensitive to variation in matrix [Ca2+] but do not match the ΔΨ dependent data in the domain ΔΨ≤120 mV. To effectively characterize this ΔΨ dependency, we reformulate the ΔΨ dependencies of the rate constants of Ca2+ translocation via the uniporter by exclusively redefining the biophysical parameters associated with the free-energy barrier of Ca2+ translocation based on a generalized, non-linear Goldman-Hodgkin-Katz formulation. This alternate uniporter model has all the characteristics of the previous uniporter model and is also able to characterize the possible mechanisms of both the extra-matrix Ca2+ and ΔΨ dependencies of mitochondrial Ca2+ uptake. In addition, the model is insensitive to variation in matrix [Ca2+], predicting relatively stable physiological operation. The model is critical in developing mechanistic, integrated models of mitochondrial bioenergetics and Ca2+ handling

Fundamental properties of Ca²⁺ signals

Background Ca²⁺ is a ubiquitous and versatile second messenger that transmits information through changes of the cytosolic Ca²⁺ concentration. Recent investigations changed basic ideas on the dynamic character of Ca²⁺ signals and challenge traditional ideas on information transmission. Scope of review We present recent findings on key characteristics of the cytosolic Ca²⁺ dynamics and theoretical concepts that explain the wide range of experimentally observed Ca²⁺ signals. Further, we relate properties of the dynamical regulation of the cytosolic Ca²⁺ concentration to ideas about information transmission by stochastic signals. Major conclusions We demonstrate the importance of the hierarchal arrangement of Ca²⁺ release sites on the emergence of cellular Ca²⁺ spikes. Stochastic Ca²⁺ signals are functionally robust and adaptive to changing environmental conditions. Fluctuations of interspike intervals (ISIs) and the moment relation derived from ISI distributions contain information on the channel cluster open probability and on pathway properties. General significance Robust and reliable signal transduction pathways that entail Ca²⁺ dynamics are essential for eukaryotic organisms. Moreover, we expect that the design of a stochastic mechanism which provides robustness and adaptivity will be found also in other biological systems. Ca2 + dynamics demonstrate that the fluctuations of cellular signals contain information on molecular behavior. This article is part of a Special Issue entitled Biochemical, biophysical and genetic approaches to intracellular calcium signaling. Highlights ► We review recent findings on key characteristics of cytosolic Ca²⁺ dynamics. ► We demonstrate the importance of the hierarchal arrangement of Ca²⁺ release sites. ► New theoretical concepts exploit emergent behavior of cellular Ca²⁺ spikes. ► We relate the dynamical regulation of [Ca²⁺] to information transmission. ► Stochastic Ca²⁺ signals are functionally robust and adaptive to changing conditions