Search

Many biological, physical, and social interactions have a particular dependence on where they take place. In living cells, protein movement between the nucleus and cytoplasm affects cellular response (i.e., proteins must be present in the nucleus to regulate their target genes). Here we use recent developments from dynamical systems and chemical reaction network theory to identify and characterize the key-role of the spatial organization of eukaryotic cells in cellular information processing. In particular the existence of distinct compartments plays a pivotal role in whether a system is capable of multistationarity (multiple response states), and is thus directly linked to the amount of information that the signaling molecules can represent in the nucleus. Multistationarity provides a mechanism for switching between different response states in cell signaling systems and enables multiple outcomes for cellular-decision making. We find that introducing species localization can alter the capacity for multistationarity and mathematically demonstrate that shuttling confers flexibility for and greater control of the emergence of an all-or-none response.

We present determinant criteria for the preclusion of non-degenerate multiple steady states in networks of interacting species. A network is modeled as a system of ordinary differential equations in which the form of the species formation rate function is restricted by the reactions of the network and how the species influence each reaction. We characterize families of so-called power-law kinetics for which the associated species formation rate function is injective within each stoichiometric class and thus the network cannot exhibit multistationarity. The criterion for power-law kinetics is derived from the determinant of the Jacobian of the species formation rate function. Using this characterization we further derive similar determinant criteria applicable to general sets of kinetics. The criteria are conceptually simple, computationally tractable and easily implemented. Our approach embraces and extends previous work on multistationarity, such as work in relation to chemical reaction networks with dynamics defined by mass-action or non-catalytic kinetics, and also work based on graphical analysis of the interaction graph associated to the system. Further, we interpret the criteria in terms of circuits in the so-called DSR-graph, Comment: To appear in SIAM Journal on Applied Dynamical Systems

We provide a Jacobian criterion that applies to arbitrary chemical reaction networks taken with mass-action kinetics to preclude the existence of multiple positive steady states within any stoichiometric class for any choice of rate constants. We are concerned with the characterization of injective networks, that is, networks for which the species formation rate function is injective in the interior of the positive orthant within each stoichiometric class. We show that a network is injective if and only if the determinant of the Jacobian of a certain function does not vanish. The function consists of components of the species formation rate function and a maximal set of independent conservation laws. The determinant of the function is a polynomial in the species concentrations and the rate constants (linear in the latter) and its coefficients are fully determined. The criterion also precludes the existence of degenerate steady states. Further, we relate injectivity of a chemical reaction network to that of the chemical reaction network obtained by adding outflow, or degradation, reactions for all species., Comment: Version updated, including change of title. Previous title: "Injectivity of chemical reaction networks with mass-action kinetics revisited"

Phosphorelays are a class of signaling mechanisms used by cells to respond to changes in their environment. Phosphorelays (of which two-component systems constitute a special case) are particularly abundant in prokaryotes and have been shown to be involved in many fundamental processes such as stress response, osmotic regulation, virulence, and chemotaxis. We develop a general model of phosphorelays extending existing models of phosphorelays and two-component systems. We analyze the model analytically under the assumption of mass-action kinetics and prove that a phosphorelay has a unique stable steady-state. Furthermore, we derive explicit functions relating stimulus to the response in any layer of a phosphorelay and show that a limited degree of ultrasensitivity (the ability to respond to changes in stimulus in a switch-like manner) in the bottom layer of a phosphorelay is an intrinsic feature which does not depend on any reaction rates or substrate amounts. On the other hand, we show how adjusting reaction rates and substrate amounts may lead to higher degrees of ultrasensitivity in intermediate layers. The explicit formulas also enable us to prove how the response changes with alterations in stimulus, kinetic parameters, and substrate amounts. Aside from providing biological insight, the formulas may also be used to avoid time-consuming simulations in numerical analyses and simulations.

We consider chemical reaction networks taken with mass action kinetics. The steady states of such a system are solutions to a system of polynomial equations. Even for small systems the task of finding the solutions is daunting. We develop an algebraic framework and procedure for linear elimination of variables. The procedure reduces the variables in the system to a set of "core" variables by eliminating variables corresponding to a set of non-interacting species. The steady states are parameterized algebraically by the core variables, and a graphical condition is given for when a steady state with positive core variables necessarily have all variables positive. Further, we characterize graphically the sets of eliminated variables that are constrained by a conservation law and show that this conservation law takes a specific form.

We define a subclass of Chemical Reaction Networks called Post-Translational Modification systems. Important biological examples of such systems include MAPK cascades and two-component systems which are well-studied experimentally as well as theoretically. The steady states of such a system are solutions to a system of polynomial equations with as many variables as equations. Even for small systems the task of finding the solutions is daunting. We develop a mathematical framework based on the notion of a cut, which provides a linear elimination procedure to reduce the number of variables in the system. The steady states are parameterized algebraically by a set of "core" variables, and the non-negative steady states correspond to non-negative values of the core variables. Further, minimal cuts are the connected components in the species graph and provide conservation laws. A criterion for when a set of independent conservation laws can be derived from cuts is given.

Multistationarity in biological systems is a mechanism of cellular decision making. In particular, signaling pathways regulated by protein phosphorylation display features that facilitate a variety of responses to different biological inputs. The features that lead to multistationarity are of particular interest to determine as well as the stability properties of the steady states. In this paper we determine conditions for the emergence of multistationarity in small motifs without feedback that repeatedly occur in signaling pathways. We derive an explicit mathematical relationship between the concentration of a chemical species at steady state and a conserved quantity of the system such as the total amount of substrate available. We show that the relation determines the number of steady states and provides a necessary condition for a steady state to be stable, that is, to be biologically attainable. Further, we identify characteristics of the motifs that lead to multistationarity, and extend the view that multistationarity in signaling pathways arises from multisite phosphorylation. Our approach relies on mass-action kinetics and the conclusions are drawn in full generality without resorting to simulations or random generation of parameters. The approach is extensible to other systems., Comment: 67 pages, 8 figures

Let $n\geq 2$ be an integer, and let $i\in\{0,...,n-1\}$. An $i$-th dimension edge in the $n$-dimensional hypercube $Q_n$ is an edge ${v_1}{v_2}$ such that $v_1,v_2$ differ just at their $i$-th entries. The parity of an $i$-th dimension edge $\edg{v_1}{v_2}$ is the number of 1's modulus 2 of any of its vertex ignoring the $i$-th entry. We prove that the number of $i$-th dimension edges appearing in a given Hamiltonian cycle of $Q_n$ with parity zero coincides with the number of edges with parity one. As an application of this result it is introduced and explored the conjecture of the inscribed squares in Hamiltonian cycles of the hypercube: Any Hamiltonian cycle in $Q_n$ contains two opposite edges in a 4-cycle. We prove this conjecture for $n \le 7$, and for any Hamiltonian cycle containing more than $2^{n-2}$ edges in the same dimension. This bound is finally improved considering the equi-independence number of $Q_{n-1}$, which is a concept introduced in this paper for bipartite graphs.

Posttranslational modification of proteins is key in transmission of signals in cells. Many signaling pathways contain several layers of modification cycles that mediate and change the signal through the pathway. Here, we study a simple signaling cascade consisting of n layers of modification cycles, such that the modified protein of one layer acts as modifier in the next layer. Assuming mass-action kinetics and taking the formation of intermediate complexes into account, we show that the steady states are solutions to a polynomial in one variable, and in fact that there is exactly one steady state for any given total amounts of substrates and enzymes. We demonstrate that many steady state concentrations are related through rational functions, which can be found recursively. For example, stimulus-response curves arise as inverse functions to explicit rational functions. We show that the stimulus-response curves of the modified substrates are shifted to the left as we move down the cascade. Further, our approach allows us to study enzyme competition, sequestration and how the steady state changes in response to changes in the total amount of substrates. Our approach is essentially algebraic and follows recent trends in the study of posttranslational modification systems.

We define a general mathematical framework for studying post-translational modification processes under the assumption of mass action kinetics., Comment: This paper has been withdrawn by the author due to a few inaccuracies in the text. An improved version of the results of this paper can be found in arXiv:1107.3531 [q-bio.MN]

In this paper we show that the regulator defined by Goncharov from higher algebraic Chow groups to Deligne-Beilinson cohomology agrees with Beilinson's regulator. We give a direct comparison of Goncharov's regulator to the construction given by Burgos and Feliu. As a consequence, we show that the higher arithmetic Chow groups defined by Goncharov agree, for all projective arithmetic varieties over an arithmetic field, with the ones defined by Burgos and Feliu.

We give a new construction of higher arithmetic Chow groups for quasi-projective arithmetic varieties over a field. Our definition agrees with the higher arithmetic Chow groups defined by Goncharov for projective arithmetic varieties over a field. These groups are the analogue, in the Arakelov context, of the higher algebraic Chow groups defined by Bloch. The degree zero group agrees with the arithmetic Chow groups of Burgos. Our new construction is shown to be a contravariant functor and is endowed with a product structure, which is commutative and associative.

We construct Adams operations on the rational higher arithmetic K-groups of a proper arithmetic variety. The definition applies to the higher arithmetic K-groups given by Takeda as well as to the groups suggested by Deligne and Soule, by means of the homotopy groups of the homotopy fiber of the regulator map. They are compatible with the Adams operations on algebraic K-theory. The definition relies on the chain morphism representing Adams operations in higher algebraic K-theory given previously by the author. In this paper it is shown that a slight modification of this chain morphism commutes strictly with the representative of the Beilinson regulator given by Burgos and Wang.