Research article

Fixed point theorems in R-metric spaces with applications

  • Received: 24 October 2019 Accepted: 10 March 2020 Published: 24 March 2020
  • MSC : 54H25, 47H10

  • The purpose of this paper is to introduce the notion of R-metric spaces and give a real generalization of Banach fixed point theorem. Also, we give some conditions to construct the Brouwer fixed point. As an application, we find the existence of solution for a fractional integral equation.

    Citation: Siamak Khalehoghli, Hamidreza Rahimi, Madjid Eshaghi Gordji. Fixed point theorems in R-metric spaces with applications[J]. AIMS Mathematics, 2020, 5(4): 3125-3137. doi: 10.3934/math.2020201

    Related Papers:

    [1] Nathaniel J. Merrill, Zheming An, Sean T. McQuade, Federica Garin, Karim Azer, Ruth E. Abrams, Benedetto Piccoli . Stability of metabolic networks via Linear-in-Flux-Expressions. Networks and Heterogeneous Media, 2019, 14(1): 101-130. doi: 10.3934/nhm.2019006
    [2] Maya Briani, Emiliano Cristiani . An easy-to-use algorithm for simulating traffic flow on networks: Theoretical study. Networks and Heterogeneous Media, 2014, 9(3): 519-552. doi: 10.3934/nhm.2014.9.519
    [3] Gildas Besançon, Didier Georges, Zohra Benayache . Towards nonlinear delay-based control for convection-like distributed systems: The example of water flow control in open channel systems. Networks and Heterogeneous Media, 2009, 4(2): 211-221. doi: 10.3934/nhm.2009.4.211
    [4] Thibault Liard, Raphael Stern, Maria Laura Delle Monache . A PDE-ODE model for traffic control with autonomous vehicles. Networks and Heterogeneous Media, 2023, 18(3): 1190-1206. doi: 10.3934/nhm.2023051
    [5] Edward S. Canepa, Alexandre M. Bayen, Christian G. Claudel . Spoofing cyber attack detection in probe-based traffic monitoring systems using mixed integer linear programming. Networks and Heterogeneous Media, 2013, 8(3): 783-802. doi: 10.3934/nhm.2013.8.783
    [6] Mahmoud Saleh, Endre Kovács, Nagaraja Kallur . Adaptive step size controllers based on Runge-Kutta and linear-neighbor methods for solving the non-stationary heat conduction equation. Networks and Heterogeneous Media, 2023, 18(3): 1059-1082. doi: 10.3934/nhm.2023046
    [7] Paola Goatin, Chiara Daini, Maria Laura Delle Monache, Antonella Ferrara . Interacting moving bottlenecks in traffic flow. Networks and Heterogeneous Media, 2023, 18(2): 930-945. doi: 10.3934/nhm.2023040
    [8] Ginestra Bianconi, Riccardo Zecchina . Viable flux distribution in metabolic networks. Networks and Heterogeneous Media, 2008, 3(2): 361-369. doi: 10.3934/nhm.2008.3.361
    [9] Riccardo Bonetto, Hildeberto Jardón Kojakhmetov . Nonlinear diffusion on networks: Perturbations and consensus dynamics. Networks and Heterogeneous Media, 2024, 19(3): 1344-1380. doi: 10.3934/nhm.2024058
    [10] Didier Georges . Infinite-dimensional nonlinear predictive control design for open-channel hydraulic systems. Networks and Heterogeneous Media, 2009, 4(2): 267-285. doi: 10.3934/nhm.2009.4.267
  • The purpose of this paper is to introduce the notion of R-metric spaces and give a real generalization of Banach fixed point theorem. Also, we give some conditions to construct the Brouwer fixed point. As an application, we find the existence of solution for a fractional integral equation.


    Quantitative Systems Pharmacology (QSP) aims to gain more information about a potential drug treatment on a human patient before the more expensive stages of development begin [22]. QSP models allow us to perform in silico experiments on a simulated metabolic system that predicts the response of perturbing a flux. A drug may be metabolized differently by various patients, and modelers working in pharmacology must anticipate these differences. Building a profile of how the drug affects different classes of simulated patients will help the developers of new drugs understand the viability of a treatment and acquire insight into the mechanisms by which the drug acts.

    A recent advancement in QSP modeling called Linear-in-Flux-Expressions (LIFE) is a method of analyzing systems of Ordinary Differential Equations (ODEs) [19,18]. Originally, LIFE was designed to analyze metabolic systems, which are composed of Fluxes and Metabolites. Fluxes in the metabolic system are the rates of chemical reactions in the human body, and they determine the dynamics on the metabolites, which are the various chemical compounds involved in metabolism. Modeling these systems depends on choosing fluxes, which are difficult to measure directly, so that the system effectively simulates human metabolism.

    To implement LIFE on a metabolic network, the network must be written as a directed graph [18]. The edges of the graph represent the reaction rates (fluxes), and the vertices represent quantities of chemical compounds (metabolites). From the graph we construct the stoichiometric matrix of the system. This stoichiometric matrix is not the classical one mentioned by [21,24]. The LIFE method is also different from QSP models whose dynamics traditionally depend on a matrix containing information about the flux of the system. In these classical QSP models the dynamics of the metabolites are linear with respect to the metabolites. By contrast, systems using the LIFE method are linear in fluxes and have a stoichiometric matrix that is dependent on the metabolites.

    Initially, the LIFE method was developed using the human cholesterol metabolism network [19]. LIFE enables us to simply describe the correlations among the fluxes of the model at steady state. There are generally many correlations among fluxes, and maintaining these correlations leads to a more consistent response to perturbing the fluxes in the system. This was advantageous to QSP modelers, who previously analyzed flux perturbations with little to no consideration to relationships among fluxes [1]. Now, we expand our study of these systems, showing that with few assumptions, systems that are linear in the flux is stable.

    The LIFE method evolved from methods in systems biology [21]. Systems biology, in conjunction with network flows [12,9], Markov chains [6], laplacian dynamics [13], control theory [3], and compartmental systems [4,14] allow us to better understand biological networks on which pharmacology models are based. The field of compartmental systems focuses on models based on directed graphs. Vertices of the graph represent quantities whose dynamics are determined by the edges of the graph, which represent fluxes among compartments. Markov chains study dynamics on directed graphs as well, but by contrast, this field focuses on stochastic processes. Control theory studies the way an external agent can alter the natural evolution of a system, given a set of admissible controls. In pharmacology, metabolism follows its natural evolution, and drugs serve as our controls. These fields have much to contribute to systems pharmacology, and we summarize useful results. We identify assumptions which are usually satisfied by real metabolic networks, guaranteeing stability of the metabolic system at a unique equilibrium.

    The paper is organized as follows. In section 2, we describe the model system for the LIFE approach in the form of a system of ODEs associated to a metabolic networks, then show existence of positive solutions and provide results of equilibria under general assumptions. Also special classes of LIFE systems are introduced. Section 3 investigates the flow vectors for which a given metabolite vector x is an equilibrium of the network, including the extreme pathways approach. On the other side, Section 4 studies the opposite problem: find the metabolite vectors which are equilibria of the network for a fixed flow vector. This is done first investigating the relationships between linear LIFE systems and Markov chains, Laplacian dynamics and linear compartmental systems. Then we deal with special classes of nonlinear LIFE systems. Finally, a comparison between zero-deficiency theory is discussed. The paper ends with conclusions in Section 5 and an Appendix containing examples.

    We indicate by R+=[0,+) the set of positive real numbers, by Rn the Euclidean real space of dimension n and by Mn×m the set of n×m matrices with real entries. Given a matrix S, we indicate by ST its transpose. Given d1,,dnR, diag(d1,,dn) is the diagonal matrix with entries di on the diagonal. We denote by 1 a column vector with all entries equal to 1, of size clear from the context.

    We introduce some terminology commonly used in graph theory. A directed graph is a couple G=(V,E), with V={v1,,vn} the set of vertices and EV×V the set of edges. For a graph with n vertices and m edges, ordering the edges lexicographically, the incidence matrix is a matrix, ΓMn×m such that Γij=1 if the jth edge is (vk,vi) for some vertex vk, Γij=1 if the jth edge is (vi,vk) for some vertex vk, and Γij=0 otherwise. A path is a sequence of distinct vertices vi1vik, with (vij,vij+1)E for j=1,,k1. A graph is strongly connected if there exists a path between every pair of vertices. A strongly connected component of a directed graph is a maximal strongly connected subgraph.

    A terminal component of a directed graph G=(V,E) is a strongly connected component G=(V,E), with VV, EE, such that there exists no edge e=(v,v), with vV and vVV. An undirected path is a sequence of distinct vertices vi1,,vik, with either (vij,vij+1)E or (vij+1,vij)E for j=1,,k1. A directed graph is weakly connected if there exists an undirected path between every pair of vertices. A weakly connected component of a directed graph is a maximal weakly connected subgraph. A directed graph G=(V,E) is weakly reversible if every weakly connected component is also strongly connected.

    We indicate by xRn the metabolite variables and by fRm the flux variables. A general system of ODEs which governs the quantities of x and f is written as

    dxdt=F(x,f), (1)
    dfdt=G(x,f), (2)

    where F:Rn×RmRn and G:Rn×RmRm. In [11,15], the authors show that the dynamics described by (1) evolve over a much smaller time-scale than (2). This is referred to as "time-scale separation". Based on time-scale separation arguments of metabolic systems, we approximate the dynamics of the fluxes with G0, and our work focuses on the dynamics of the metabolites (1), with the fluxes playing the role of constant parameters.

    The dynamics (1) is very general and we restrict to special system, which are linear in the fluxes, thus can be written as:

    ˙x=S(x)f (3)

    where S:RnMn×m is called the stoichiometric matrix. One constructs the stoichiometric matrix from the metabolites and the reactions that comprise a biochemical system. Each reaction corresponds to a flux f that connects two distinct metabolites or represents an intake or an excretion from the network. Each row of S corresponds to a metabolite and each column of S corresponds to a flux.

    We now illustrate how to construct a directed graph from the metabolic network for the system (3). We represent metabolites with vertices V={v1,,vn}. We construct a set of edges EV×V to represent reactions; each edge is associated to a flux from one metabolite to another, notice that we do not have loops. To represent intakes and excretions, we introduce two virtual vertices, v0 and vn+1, not associated with any metabolite but rather representing the external environment. We denote by I,X the set of vertices attached to v0,vn+1, the vertices in I and X are called intake vertices and excretion vertices, respectively. We also introduce edges (v0,w) with wIV representing intakes, and (w,vn+1), wXV representing excretions. We use the extended graph ˜G=(˜V,˜E) defined by ˜V=V{v0,vn+1}={v0,v1,,vn,vn+1} and ˜E collecting edges in E together with intake and excretion edges. The rows of the matrix S can be indexed by vertices in V and the columns by edges in ˜E, thus we write Sve for the entry corresponding to vertex v and edge e. Moreover we denote by xv the metabolite corresponding to vertex v and by fe the flux corresponding to edge e.

    Example 2.1. To illustrate the concepts of graph with virtual vertices and stoichiometric matrix related to a metabolic network, we provide a toy example with linear dynamics. Consider the system given by the following stoichiometric matrix and fluxes vector:

    S(x)=(1x1x1000x40x10x200000x1x2x3x300000x30x4),f=(f(v0,v1)f(v1,v2)f(v1,v3)f(v2,v3)f(v3,v4)f(v3,v5)f(v4,v1)).

    Then the corresponding graph is represented in Figure 2.

    Figure 2. 

    A directed graph ˜G=(˜V,˜E) illustrating Proposition 2. Vertices v3 and v4 form a terminal component. There exists a path from v0 to v4 yet there is no path from v4 to v5

    .

    Remark 1. It is worth a reminder that our stoichiometric matrix is different from the traditional one defined by [21,24], in which entries are stoichiometric coefficients, i.e. do not depend on metabolites.

    To correctly represent the reactions corresponding to fluxes (which take always strictly positive values), we assume:

    (A) For x(R+)n, it holds

    Sve(x)={He(x)>0if e=(w,v), wV and xv>0 or e=(v0,v), vIHe(x)<0if e=(v,w), wV and xv>0 or e=(v,vn+1), vX and xv>00otherwise,
    Figure 1. 

    A directed graph ˜G=(˜V,˜E) representing a biochemical system. The rectangles indicate virtual vertices and the subgraph of circular vertices and edges connecting them is G=(V,E)

    .

    where He:RnR is a positive continuous function.

    Notice that Assumption (A) implies that, for each vV,

    vVSve(x)={He(x)e=(v0,ˉv), ˉvI,He(x)e=(ˉv,vn+1), ˉvX,0otherwise, (4)

    namely all columns of S have zero sum, except those corresponding to intakes and excretions, which have positive and negative sum, respectively. Under Assumption (A), the dynamics (3) can be interpreted as mass conservation law. Indeed, re-writing (3) entrywise and using (A), we have

    ˙xv=e˜ESve(x)fe=w:(w,v)˜EH(w,v)(x)f(w,v)w:(v,w)˜EH(v,w)(x)f(v,w), (5)

    which is the mass balance for metabolite xv: its variation is given by the sum of the incoming flows, minus the sum of the outgoing flows. This is the analogous Kirchhoff's current law for electrical networks, with the difference that currents are allowed to take negative values as well, while here metabolite variables are non-negative.

    The total mass in the system is m=vVxv. From (5) we have

    ˙m=vIH(v0,v)(x)f(v0,v)vXH(v,vn+1)(x)f(v,vn+1).

    Clearly, in the case without intakes nor excretions, ˙m=0, i.e. the total mass of a closed system is constant in time.

    Another remark which is useful is that, under Assumption (A), S(x) can be re-written as S(x)=ΓD(x), where D(x) is a diagonal matrix of size m×m, with diagonal entries given by He(x)'s, and Γ is obtained from the incidence matrix of ˜G by removing the first and last rows (corrseponding to v0 and vn+1). In the particular case without intakes nor excretions, Γ is the incidence matrix of G.

    In the remainder of this section we study the dynamics (3) under the very general Assumption (A), while later in the paper we add other assumptions, restricting our attention to systems for which stronger statements can be obtained. A first important general property is that positivity of solution is guaranteed:

    Proposition 1. Consider a system (3) satisfying (A) and the Cauchy problem with initial datum xv(0)=xv0. Assume that S is locally Lipschitz. If fe>0 for every eE and xv00 for every vV, then there exists a local solution xv() defined on [0,T], T>0, and xv(t)0 for every t[0,T].

    Proof. Existence follows from Lipschitz condition, while positivity of solution follows from the invariance of the set {x:xv0}.

    In the next Proposition we show that existence of nontrivial equilibria implies some structure on the network: every vertex v for which there is a directed path from some wI to v, must also have a directed path from v to some yX. This result refines the space of networks with which we are concerned. More precisely:

    Proposition 2. Consider a system (3) satisfying (A). Assume there exists an equilibrium ˉx(R+)n for a flux vector f such that fe>0 for every e˜E. Then for every vertex vV for which there exists a path from I to v, there exists a path from v to X.

    Proof. Assume there exists an equilibrium ˉx as in the statement and, by contradiction, a vertex v for which there exists a path from wI to v, but there exists no yX to which v is connected. Since there is no path from v to X, either v belongs to a terminal component with no excretion, or there is a path from v to a terminal component with no excretion. Denote by GT=(VT,ET) such a terminal component. Since there are a path from wI to v and a (possibly trivial) path from v to VT, then there is also a path from v0 to VT. Denote by v0,v1=w,,v1,vVT one such a path, such that v1VT (possibly the path is a single edge, in case wVT).

    It is easy to show that xvi>0 for all i=1,,, as follows. Considering e=(v0,v1), by (A) we have Sv0,e(ˉx)=He(ˉx)>0. Similarly, for every e=(v1,w), xv1=0 implies Sv1,e(x)=0. This means we must have ˉxv1>0, otherwise we would have ˙xv1feHe(ˉx)>0 (where e=(v0,v1)), contradicting ˉx being an equilibrium. Having proved that ˉxv1>0, we can proceed by induction: for i=1,,1, ˉxvi>0 implies ˉxvi+1>0. The argument is the same as above, with a slight modification: looking at e=(vi,vi+1), Svi,e(ˉx)=He(ˉx)>0 thanks to Assumption (A) together with ˉxi>0, while above we were in the case of an intake.

    Finally we have a terminal component GT=(VT,ET) with no excretion, and an edge (v1,v) with vVT and v1VT, such that either v1=v0 or ˉxv1>0. In either case, considering e=(v1,v), by (A) we have Sv1,e(ˉx)=He(ˉx)>0. Now consider the variation of mass in the component GT: since there is no outgoing edge from GT, and there is at least the incoming edge e, we have ddtvVTxv=vVT˙xvHefe>0, contradicting the fact that ˉx is an equilibrium.

    The system associated to the graph in Figure 2 provides an explicit example of the contradiction argument of Proposition 2. The vertices v3 and v4 violate the condition of Proposition 2, since there is a path from the intake vertex v1 to v3 and v4, and there is no path from v3 and v4 to the excretion vertex v2; hence, the system does not admit any equilibrium. The proof argument, specialized to this example, is to notice that v3 and v4 form a terminal component with no excretion, and to look at the path v0,v1,v2,v4. Assuming that there is an equilibrium ˉx, one shows first that ˉxv1>0 due to the intake, and then that ˉxv1>0 implies also ˉxv2>0. Notice that ˙xv3+˙xv4Sv2,(v2,v4)fv2,v4 (in this particular example, equality is actually true), and the latter is >0 since ˉxv2>0, thus contradicting the fact that ˉx is an equilibrium: the mass xv3+xv4 grows unbounded.

    The same system of Figure 2 shows the necessity of the assumption that He(x)=0 when xv=0, e=(v,w) for Proposition 2 to hold. Indeed assume that H(v2,v5)(x)=1 for every x such that xv2=0 and that the other functions He, e=(v,w), satisfy He(x)=xv. Let f be a flux vector with strictly positive components such that f(v1,v2)=f(v2,v5). Then ˉx=(f(v0,v1)f(v1,v2),0,f(v4,v3)f(v3,v4),1) is an equilibrium.

    Another example not satisfying the assumptions of Proposition 2 is shown in Figure 3. This system admits an equilibrium under Assumption (A) even though the vertex v3 does not have a path to v5. This is because there does not exist a path from v0 to v3.

    Figure 3. 

    A directed graph where vertices v3 and v4 do not have a path from v0 and also have no path to v5. For an equilibrium, ˉx, of this system under Assumption (A), ˉxv4=0 and ˉxv3xv3(0)

    .

    For general systems (3) with Assumption (A) it is not possible to prove other general conclusions about equilibria, beside Propositions 1 and 2. Indeed consider the simple metabolic network given in Figure 4. Under Assumption (A), the dynamics is written as:

    Figure 4. 

    A directed cycle graph G=(V,E) with n vertices and no intakes nor excretions. On such a LIFE system one can prescribe any desired dynamics

    .
    {˙xv1=f(v1,v2)H(v1,v2)(x)+f(vn,v1)H(vn,v1)(x)˙xv2=f(v2,v3)H(v2,v3)(x)+f(v1,v2)H(v1,v2)(x)˙xv3=f(v3,v4)H(v3,v4)(x)+f(v2,v3)H(v2,v3)(x)˙xvn=f(vn,v1)H(vn,v1)(x)+f(vn1,vn)H(vn1,vn)(x). (6)

    We want to show that for any dynamical system in Rn on a compact set there exists an equivalent dynamics defined on the cycle graph of Figure 4. In other words, Asumption (A) is so general that we give Proposition 3 to show arbitrary dynamics can be defined, and we focus on more specialized cases (Assumption (B), and (C)). More precisely for every general dynamics

    {˙xv1=F1(xv1,xv2,,xvn)˙xv2=F2(xv1,xv2,,xvn)˙xvn=Fn(xv1,xv2,,xvn). (7)

    we look for a choice of the functions He and the fluxes fe realizing such equivalence. Note that from (4), we have xvn=Cx1xn1 which implies that

    ˙xvn=˙xv1˙xvn1=n1i=1Fi. (8)

    Define the set T=ni=1{x:xj0 for all j, xi=0 and x1++xn1} then we have the following:

    Proposition 3. Consider system (7) and assume F1==Fn=0 on the set T. Then there exist functions He and fluxes fe such that the dynamics (6) is equivalent to (7) on the bounded set delimited by T.

    Proof. Assign the functions H(vi,vj)(x) according to the following rule:

    H(vi,vi+1)(x)=ki[Fk]+>i[F]+if  i<n,H(vn,v1)(x)=nk=1[Fk]+, (9)

    where [F]+=max{F,0} and [F]=max{F,0}. It is easy to verify that ˙xvi=Fi(x) for i<n and ˙xvn=ni=1Fi=Fn(x) using (8). Moreover, because of the assumption on Fi's, the functions He satisfy (A), thus we are done.

    Here we introduce a special class of systems of type (3) with simplified dynamics. We consider Assumption (A), and we impose further restriction on the functions He(x): for an edge e=(v,w), we assume He to depend on xv only, and moreover we impose the scalar function He(xv) to be strictly increasing. More precisely, Assumption (B) is the following.

    (B) It holds

    Sve(x)={He(xv)e=(v,w), vV,wV{vn+1}He(xw)e=(w,v), wV1e=(v0,v) vI0otherwise, (10)

    and each He is a positive, differentiable, and strictly increasing function He:RR, with He(0)=0.

    A typical example of a system verifying (B) is given by metabolic networks with Hill functions representing reactions, i.e. He(xv)=xpevK+xpev with peN, where K>0 is the dissociation constant.

    A further simplification occurs in the case where He(xv) is the same function Hv(xv) for all edges e having v as an initial vertex, i.e. such that e=(v,w) for some w. This gives Assumption (C), as follows.

    (C) It holds

    Sve(x)={Hv(xv)e=(v,w), vV,wV{vn+1}Hw(xw)e=(w,v), wV1e=(v0,v) vI0otherwise,

    and each Hv is a positive, differentiable, and strictly increasing function Hv:RR, with Hv(0)=0.

    Under Assumption (C), the system ˙x=S(x)f can be equivalently re-written as

    ˙x=J(f)h(x)+ϕ, (11)

    where J(f)Mn×n is defined by

    Jij(f)={f(vj,vi)if (vj,vi)Ew:(vi,w)˜Ef(vi,w)ifj=i0otherwise,

    where h(x) is a vector of size n given by hi(x)=Hvi(xvi) and ϕ is a vector of size n given by ϕi=f(v0,vi) if (v0,vi)˜E, and ϕi=0 otherwise.

    Finally, the simplest class of LIFE models we consider are linear systems, namely systems satisfying Assumption (C), where each Hv(xv) is the identity function Hv(xv)=xv.

    Example 2.1 (continued). This example is a linear LIFE system. We can equivalently re-write its dynamics ˙x=S(x)f as ˙x=J(f)x+ϕ, with

    J(f)=(f(v1,v2)f(v1,v3)00f(v4,v1)f(v1,v2)f(v2,v3)00f(v1,v3)f(v2,v3)f(v3,v4)f(v3,v5)000f(v3,v4)f(v4,v1)),ϕ=(f(v0,v1)000).

    In this section we consider equilibrium solutions of the system (3) satisfying Assumption (A). In general one is interested in conditions guaranteeing existence of an equilibrium and also in conditions necessary for uniqueness and stability of such an equilibrium. There are two problems: for fixed metabolite concentrations x find all flux vectors f for which x is an equilibrium, i.e. ˙x=S(x)f=0, and, vice versa, for a fixed flux vector f, find all x that are equilibria. In this section we focus on the first, while the latter is investigated in Section 4.

    The set of flux vectors for which x is an equilibrium formed by all vectors f that solve the equation S(x)f=0, i.e. the nullspace N(S(x)) of S(x). First, we discuss the dimension of N(S(x)). Then, since fluxes must be positive (to have a correct biological meaning), we focus on describing the cone N(S(x))(R+)m. Recall that, under Assumption (A), we can rewrite S(x) as S(x)=ΓD(x) where D(x)Mm×m is a diagonal matrix with He(x)'s as entries and Γ is obtained from the incidence matrix of ˜G by removing the first and last rows; in case there are no intakes nor excretions, then Γ is the incidence matrix of G. Assuming that x has strictly positive entries, (A) implies that all diagonal elements of D(x) are strictly positive, so that D(x) is invertible. Hence, for Γ(D(x)f)=0 to have non trivial solutions, the nullspace of Γ must have dimension greater than zero.

    We extend one of the results of [18] to get:

    Proposition 4. Consider the system (3) with no intakes nor excretions satisfying Assumption (A) and let G be the the associated graph. Let xRn have strictly positive entries and S(x)Mn×m be the stoichiometric matrix, then

    rank(S(x))=n,

    where is the number of weakly connected components of G.

    Proof. We have S(x)=ΓD(x). As shown in Proposition 4.3 of [2] the rank of an incidence matrix Γ is rank(Γ)=n. Since D(x) is full rank we have rank(S(x))=n.

    Next we consider systems that contain both intakes and excretions. It was shown in [18] that if intakes and excretions are added to a graph satisfying the conditions of Proposition 4 then rank(S(x))=n. We extend this result to get:

    Proposition 5. Consider the system (3) satisfying Assumption (A) and let G be the the associated graph. Let xRn have strictly positive entries and S(x)Mn×m be the stoichiometric matrix, then rank(S(x))=nk, where k is the number of weakly connected components containing neither intake nor excretion vertices.

    Proof. It was shown in [18] that if G is weakly connected and contains an intake vertex then rank(S(x))=n. The same argument from [18] can be also used for an excretion vertex and so if G is weakly connected and contains an excretion vertex we also have rank(S(x))=n. In the case where G is not weakly connected S(x) can be rewritten as a block diagonal matrix

    S(x)=(S1(x)00S(x))

    where each diagonal element Si(x) is the stoichiometric matrix of a weakly connected component. If Si(x) contains an intake or excretion vertex then rank(Si(x))=ni, where ni is the number of metabolites in Si(x). If Si(x) contains neither intake nor excretion vertices then rank(Si(x))=ni1. This implies that S(x) has ni linearly independent rows for each Si(x) with intake or excretion vertices and ni1 linearly independent rows for each Si(x) with neither intake nor excretion vertices. Let k represents the number of Si(x) with neither intake nor excretion vertices, the total number of linearly independent rows in S is now expressed (i=1(ni))k=nk. Thus rank(S(x))=nk.

    Notice that the rank of S(x) depends on the number of weakly connected components of the graph, which is the same irrespective of the orientation of edges. However, when we focus on the existence of non-trivial positive flows admitting an equilibrium with strictly positive entries, the orientation of edges does matter, as we can see e.g. from Proposition 2, where a necessary condition is given for existence of equilibria, in terms of existence of suitable paths.

    The problem of finding positive flows admitting an equilibrium with strictly positive entries has been extensively studied in the operations research literature under the name of 'network flow' problems [12]. In network flow problems one considers a directed graph where edges represent flows between the vertices. The maximum flow which an edge can support is called the capacity of the edge. In addition to capacity each edge may also have a cost associated to it. Network flow problems assume that the system is at equilibrium with respect to the vertices i.e. that the flow entering and leaving a vertex must be the same. A flow that satisfies this assumption is called a feasible flow. Finding feasible flows is the same as finding flows in the cone N(S(x))(R+)m with the further constraint that each flux must be less than its capacity.

    A flow is a mapping from f:E(R+)m that satisfies 0f(vi,vj)c(vi,vj) (where c(vi,vj) is the capacity of the edge (vi,vj)) and (vi,vj)Ef(vi,vj)(vj,vk)Ef(vj,vk)=0 [17]. One of the most common network flow problems is to find the maximum flow of the network, i.e. the largest amount of total flow from a source to a sink. In [9] the authors consider a network which contains exactly one source (v0) and one sink (vn+1). They then prove the following result, known as max-flow min-cut theorem, which characterizes the maximum flow as the minimum cut cacapity, where a cut set is a set of edges whose removal disconnects the source from the sink, and the capacity of a cut is the sum of the capacities of its edges.

    Proposition 6 (Max-Flow Min-Cut Theorem). The maximum flow value obtainable in a network is the minimum capacity of all cut sets which disconnect v0 and vn+1.

    An implication of the max-flow min-cut theorem is that a feasible flow exists if there is a path from v0 and vn+1. The following proposition extends this relationship:

    Proposition 7. Given a LIFE system with graph ˜G and stoichiometric matrix S(x) satisfying (A), fix an x(R+)n with strictly positive entries and for every vI fix ˉfv0,v>0, i.e. fix arbitrary values for the intake flows. There exists fN(S(x))(R+)m such that fv0,v=ˉfv0,v for all vI if and only if for each vI there exists a path to X.

    Proof. Consider the maximum flow problem on the graph ˜G, where edges (v0,v) have capacity ˉfv0,v, and all other edges have infinite capacity. The feasible flows φ for this network are in one-to-one correspondence with the equilibrium flows fN(S(x))(R+)m such that fv0,vˉfv0,v for all vI; the correspondence is simply given by fv0,v=φv0,v/Hv0,v(x) for all vI and fv,w=φv,w/Hv,w(x) for all wV. If for all vI there is a path from v to X (and hence a path from v to vn+1), the minimum cut is the removal of all edges (v0,v). The maximum flow φ then has φv0,v=ˉfv0,v, thus also ensuring the existence of an equilibrium flow f satisfying the same.

    If for some viI there is no path from v to X, then all feasible flows φ satisfy φv0,v=0, and hence all equilibrium flows f satisfy fv0,v=0 which contradicts the assumption.

    There are many standard methods for computing a basis of the nullspace of a matrix [20], and hence to describe N(S(x)). Since we are interested in positive fluxes, we are rather interested in the cone N(S(x))(R+)m. One method to describe this cone is to look for a positive basis, which is a minimal set of vectors generating the whole cone via linear combinations with positive scalars, called positive combinations. The use of positive combinations ensures that all generated vectors belong to the cone. Minimality is equivalent to ask for the vectors to be positively linearly independent, i.e. no vector can be expressed as a positive combination of the others. Since we allow only positive combinations (and not every linear combination), typically a positive basis has more vectors than a traditional basis. In [23] the authors prove that a flux cone has a positive basis, and that the vectors of the positive basis are unique up to multiplication by a positive scalar.

    Finding a positive basis for a cone is not as easy as finding a traditional basis. In [21], Palsson describes a method to find a positive basis by using extreme pathways of a metabolic network. The extreme pathways are a unique set of positively independent flux vectors that represent the edges of the cone N(S(x))(R+)m [21,23].

    Here we summarize a method to build extreme pathways described by [23]. We start considering the equation V0S(x)T=C0, and the solution given by V0=Im, the m×m identity matrix, S(x) the stoichiometric matrix and C0=S(x)T. At each iteration new matrices ViMni×m, CiMni×n are defined which satisfy ViS(x)T=Ci.

    Algorithm 1. The algorithm consists of the following steps:

    Step 1. Select from Ci the first non-zero column (lexicographically), whose corresponding metabolite is neither an intake vertex nor an excretion vertex, say ji. If no such column exists jump to Step 5.

    Step 2. Add to Ci all possible new rows, obtained using positive linear combinations of two other rows, which have a zero on the column ji. Then define Ci+1 by removing the old rows used to generate the new ones. Notice that Ci+1 has all zeros in column ji.

    Step 3. Define Vi+1 by adding to Vi the rows generated by the same combination as those of Step 2 (to keep the validity of the equation VS(x)T=C) and removing the old rows.

    Step 4. Remove positively linearly dependent rows from Vi+1. Remove the corresponding rows from Ci+1.

    Step 5. Repeat the steps 1-4 now considering first excretion vertices and then intake vertices until all columns are 0.

    If the algorithm stops at step m, then Cm=0 (the matrix with all zero entries). The rows of Vm are the vectors of a positive basis. Notice also that the number of rows ni may increase or decrease during the various steps.

    Originally in [21] the extreme pathways method was proposed for stoichiometric matrices not dependent on x. We apply it here to the most general case of systems (3) with Assumption (A), by fixing a desired metabolite variables vector x(R+)n and using the extreme pathways method to characterize all fluxes vectors with positive entries such that the corresponding LIFE system (3) admits x as an equilibrium.

    To illustrate the process we report an example using the Reverse Cholesterol Transport Network (RCT) from [19]. The RCT network is shown in Figure 5.

    Figure 5. 

    Reverse Cholesterol Transport Network from [19]. This network contains 6 vertices which represent metabolites, 10 edges which represent fluxes and 2 virtual vertices v0,vn+1. There are three intake vertices v1,v2,v3 and 1 excretion vertex v6

    .

    Example 3.1. Finding Extreme Pathways of RCT network

    We use Palsson's algorithm for finding extreme pathways on the stoichiometric matrix from the RCT network. The stoichiometric matrix for RCT is given by:

    S(x)=(100x10000000100x20000000100x30000000x1x2x3x4x400000000x40x500000000x4x5x6)

    For space, the complete calculations are continued in the Appendix.

    Let us first recall Farkas' Lemma [10]. Here we use the notation x0 to indicate that every entry in the vector x is positive.

    Lemma 3.1 (Farkas' Lemma). Let AMn×m and bRn. For the equation Ax=b exactly one is true:

    1. There exists xRm such that Ax=b and x0.

    2. There exists yRn such that ATy0 and bTy<0.

    Proposition 8. The extreme pathways of S(x) form a positive basis of N(S(x))(R+)m.

    Proof. There are three statements which must be shown to prove this Proposition.

    Statement 1: The extreme pathway vectors (vi) are positively independent (systematically independent in terminology used in [21]).

    Statement 2: The span of the extreme pathways vectors is contained in (N(S(x))(R+)m) i.e. λi0λivi(N(S(x))(R+)m).

    Statement 3: The extreme pathways vectors span (N(S(x))(R+)m) i.e. (N(S(x))(R+)m)λi0λivi.

    Statement 1 follows from Step 4 of the algorithm.

    Let us now prove Statement 2. Since the vectors corresponding to extreme pathways are positive combinations of rows of matrices with positive entries, then vi(R+)m for all i. We have VnS(x)T=Cn=0, equivalently S(x)VTn=CTn=0. This implies that the columns of VTn in the nullspace of S(x), thus the rows of Vn in the nullspace of S(x). Since both N(S(x) and (R+)m are closed for positive linear combinations we are done.

    Finally it remains to prove statement 3, in order to show (N(S(x))(R+)m)λi0λivi we proceed as follows. First we show that after one iteration of the algorithm the rows of V1 form a positive basis for vectors w such that the j1-th element of S(x)w equals zero. After n iterations, we show that Vn is a positive basis for vectors wN(S(x)).

    For simplicity of notation, we assume j1=1. Define W1 to be the set of vectors w such that the first entry of S(x)w equals zero. It is clear from the construction of V1 that its rows belong to W1. It remains to be shown that the rows of V1 span W1. Since V1S(x)T=C1, thus we can alternatively show that the columns of CT1 span the subset of the range of S(x) with first entry equal to zero. In other words, we must show that if λi0λici,1=0, where ci,1 represents the first entry in the ci-th row of C0, that the sum λi0λici=0 can also be represented using rows from C1. We assume an ordering such that ci,10 for i and ci,1<0 for i>, each row in C1 can be represented by ci+αi,jcj where i,j> and αi,j=ci,1cj,10. We need to find μi,j such that,

    μi,j(ci+αi,jcj)=ni=1λici. (12)

    We can split the sum on the right hand side by considering the positive and negative entries separately,

    μi,j(ci+αi,jcj)=i=1λici+ni=+1λici.

    For λi, i{+1,,n} we can write:

    k=1μk,iαk,i=λi

    which gives:

    k=1μk,i (ck,1)=ci,1λi,i{+1,,n}. (13)

    For λi, i{1,,} we can write:

    nk=+1μi,k=λi

    which gives:

    nk=+1μi,kci,1=ci,1λi,i{1,,}. (14)

    The equations (13), (14) can be written in the following way,

    00 rows0{n rows{([c1,1]000[c2,1]000[c,1]diag(c1,1)diag(c2,1)diag(c,1))(μ1,+1μ1,+2μ1,nμ2,+1μ2,nμ,n)=(c1,1λ1c2,1λ2cn,1λn) (15)

    where [ci,1] represents a 1×(n) vector with all entries equal to ci,1 and diag(ci,1)M(n)×(n) a diagonal matrix with ci,1 on the diagonal. Let A, b be the matrix on the left-hand side and the vector on the right-hand of equation (15), respectively. We can then apply Farkas' Lemma. Consider a vector yRn such that ATy0, then for the first row we have c1,1y1c1,1y+10. Further for any i=1,, and j=+1,,n we have ci,1yici,1yj0 which implies yiyj. Setting ˉy=max(y+1,,yn) we have yiˉy for i=1,,. Since all cj,1<0 for j=+1,,n and all λj0 we have that yjcj,1λjˉycj,1λj. Next we define ˉy=(ˉy,,ˉy)Rn and note that bTybTˉy. Since λi0λici,1=0 we have bTˉy=0. Therefore condition 2 of Farkas' Lemma fails. This implies that condition 1 is true and thus there exists a set of positive μ's that solve the system (15). Since a set of μ's can be found we have that after the first iteration of the algorithm the rows of V1 span W1. By similar argument the rows of the Vi in subsequent iterations span the set of vectors w such that the first i entries of S(x)w equal zero. After n iterations, the rows of Vn span N(S(x))(R+)m, thus we are done.

    Example 3.2. Comparison of the standard and positive bases

    Using the RCT network, an example is presented which compares a standard basis of the nullspace and positive basis for the nullspace. Here we provide the calculations needed to show Statement 2 and Statement 3 of Proposition 8 for the RCT. The full details for this example are contained in the Appendix.

    In this section we characterize equilibria and the stability of LIFE systems with a fixed flux vector. We start with linear systems, and then we consider special LIFE systems, with Assumptions (B) and (C). LIFE systems satisfying at least Assumption (A) are known as compartmental systems in the automatic control community [14,4]. We build upon the well-established results on linear compartmental systems to get a full understanding of linear LIFE systems, as well as for special LIFE systems satisfying Assumption (C). Notice that weakly connected components of ˜G correspond to subsystems having no interaction with each other, so that we can study each weakly connected component separately from the others. Hence, throughout this section, we assume ˜G is weakly connected, without loss of generality.

    In this section we study the properties of linear systems, i.e. special LIFE system (3) satisfying Assumption (C) with fluxes Hv(xv)=xv, or, equivalently, a system (11), with h(x)=x. We first focus on metabolic networks with no intakes nor excretions, recalling results from a rich literature from different communities. For the case with intakes and excretions, we make use of results for Compartmental Systems [4].

    In the case with no intakes nor excretions (also known as free closed system), the linear LIFE dynamics become ˙x=J(f)x. Notice that J(f) is a Metzler matrix (i.e. has non-negative off-diagonal entries), and all its columns sum to zero (1TJ(f)=0T).

    This system is known as Laplacian dynamics, because L=J(f)T is a weighted Laplacian of the graph G defined as follows. Given the graph G and given strictly positive weights fe associated with its edges (for us, the weights are the fluxes), the weighted adjacency matrix A is defined by Aij=f(vi,vj) if (vi,vj)E, and Aij=0 otherwise. The corresponding weighted Laplacian is L=DA, where D=diag(A1) is the diagonal matrix containing the weighted out-degrees, i.e. row-sums of A.

    The eigenvalues of the Laplacian, and in particular its eigenvalue 0 and the corresponding eigenspace which is actually the set of equilibria of ˙x=J(f)x, have been studied in graph theory (see e.g. [5]). The study of the spectrum of the Laplacian has received extensive attention also in the automatic control community, on one hand because of the recent interest in the so-called consensus dynamics (see [4,Chapter 7]) ˙x=Lx (which are different from the dynamics ˙x=J(f)x, since L=J(f)T, but the eigenvalues of L and J(f) are the same), and on the other hand because of the interest in the linear compartmental system ˙x=J(f)x (see [4,Chapter 9]).

    Laplacian dynamics have been introduced in the mathematical biology literature by [13], together with a complete study of their equilibria and convergence properties. Notice that in [13] the matrix J(f) itself is called a Laplacian, while in the graph theory and control theory communities the name Laplacian refers to L=J(f)T.

    The dynamics ˙x=J(f)x is also very related to continuous-time Markov chains (see e.g. the textbook [6]). Thus we elaborate on this, a homogeneous continuous-time Markov chain over a finite state V={1,,n} is a stochastic process X(t) taking values in V, that satisfies the Markov property:

    for all times t0t1thth+1,

    Pr(Xth+1=ih+1Xt0=i0,Xt1=i1,,Xth=ih)=Pr(Xth+1=ih+1Xth=ih),

    and which is homogeneous in time:

    Pr(Xth+1=ih+1Xth=ih)=Pr(Xth+1th=ih+1Xt0=ih)=Pih,ih+1(th+1th).

    P(t) is the matrix whose (i,j)-th entry represents the probability of transition from state i to state j in an interval of time of length t. P(t) is assumed to be right-differentiable, and with time derivative defined by ˙P(t)=limh0+(P(t+h)P(t))/h. The Markov chain is fully described by its generator matrix (or transition rate matrix) defined by Q=˙P(0). Indeed, P(t) is related to Q by the Kolmogorov forward equation ˙P(t)=P(t)Q. The unique solution of the Cauchy problem ˙P(t)=P(t)Q with P(0)=I is P(t)=eQt. Denoting by π(t) a column vector whose ith entry is πi(t)=Pr(X(t)=i), we have that πT(t)=πT(0)P(t)=πT(0)eQt and hence is a solution of ˙πT=πTQ with given initial condition π(0).

    The generator matrix Q is a Metzler matrix whose rows sum to 0 (Q1=0), and hence there is an immediate equivalence between the dynamics ˙πT=πTQ and the linear LIFE system ˙x=J(f)x, simply by taking Q=J(f)T. In Markov chains, one looks at π being a probability vector, namely having positive entries and iπi=1. In LIFE systems, we are also interested in a vector x with positive entries, but the total mass could be arbitrary, so that one should set πi(0)=xvi(0)/m0, with m0=ixvi(0).

    Based on all this rich literature, we recall here the results about equilibria and their stability. The spectrum of J(f) is characterized as follows (see e.g. [13,Proposition. 1]):

    Proposition 9. Assume there are no intakes nor excretions, then:

    All eigenvalues of J(f) are either 0 or have strictly negative real part.

    The dimension of the nullspace of J(f) is equal to the algebraic multiplicity of the 0 eigenvalue 1, and is equal to the number of terminal components in G.

    1The algebraic multiplicity of an eigenvalue is its multiplicity as a root of the characteristic polynomial, while its geometric multiplicity is the dimension of the corresponding eigenspace. In the case of the zero eigenvalue, the eigenspace is the nullspace of the matrix.

    Moreover, denoting by G1,,Gk the terminal components of G, there exists a basis of the nullspace of J(f) composed of vectors π1,,πk, such that vector πi has strictly positive entries in correspondence of vertices in Gi, and has vanishing entries otherwise.

    It is customary to 'normalize' the vectors π1,,πk so that their entries sum to 1, so that they can be interpreted as probability vectors. Under this choice, the restriction of vector πi to the terminal component Gi is the stationary distribution of the Markov chain restricted to such component, i.e. its unique equilibrium with mass 1 (uniqueness is obtained from the fact that the terminal component is a strongly connected component).

    Notice that the dimension of the nullspace of an incidence matrix Γ (and hence of a stoichiometrix matrix S(x)) is related to the number of weakly connected components, while the dimension of the nullspace of a Laplacian matrix is related to strongly connected components, and not all components matter, but only the terminal ones.

    In the case where G has a unique terminal component, the nullspace of J(f) has dimension one: the equilibrium is unique, up to a multiplicative factor which is the initial total mass m0=ixvi(0). In the case where G is strongly connected, the whole graph is a single strongly connected component. In this case, not only the nullspace of J(f) has dimension one, but also it is generated by a vector whose entries are all strictly positive.

    The spectral properties of J(f) given in Proposition 9 fully characterize the asymptotic behavior of the linear system ˙x=J(f)x, by standard theory of linear systems, giving the following:

    Proposition 10. Consider the linear LIFE system with no intakes nor excretions and denote by G1,,Gk the terminal components of the associated graph G. The following properties hold:

    The total mass of the system m=vVxv=m0 is constant in time.

    From any positive initial condition x(0), the system converges to an equilibrium, having strictly positive entries in correspondence of vertices of terminal components, and 0 elsewhere.

    Moreover, if there is a unique terminal component, then the equilibrium is uniquely determined by the initial total mass m0.

    We now focus on linear systems with intakes and/or excretions, using results from linear compartmental systems as summarized in [4,Chapter 9]. A first remark is that, having introduced a single v0 from which all intake edges are originated, and a single vn+1 to which all excretion edges are headed, we have a single weakly connected component containing intakes and/or excretions. We restrict our attention to such a weakly connected component, while other possible components with no intakes nor excretions have a behavior described in Sect. 4.1.1.

    The dynamics are given by ˙x=J(f)x+ϕ, where vector ϕ represents the intakes. Due to excretions J(f) does not have all column-sums equal to 0. This means that J(f)T is not any more a Laplacian, but it is a grounded Laplacian. The term grounded Laplacian refers to a matrix obtained from a larger Laplacian matrix by deleting the row and column corresponding to a given vertex; the name 'grounded' has an interpretation for electrical networks, where this corresponds to connecting the given vertex to the ground. Consider the subgraph of ˜G where we have removed v0 but not vn+1, and consider the weighted Laplacian LM(n+1)×(n+1) of such a graph (with weights equal to the fluxes), then define Lg by deleting the last row and last column (associated with vn+1). The resulting grounded Laplacian is Lg is such that J(f)=LTg.

    The spectral properties of J(f) are summarized as follows (see [4,Theorem 9.5 and Lemma 9.12]):

    Proposition 11. Consider a linear system with intakes and/or excretions, then:

    All eigenvalues of J(f) are either 0 or have strictly negative real part.

    The dimension of the nullspace of J(f) is equal to the algebraic multiplicity of the 0 eigenvalue, and is equal to the number of terminal components not containing any excretion.

    In particular, the following are equivalent:

    (a) For every vV there is a path from v to X.

    (b) J(f) is Hurwitz stable (i.e. all its eigenvalues have strictly negative real part).

    (c) J(f) is invertible.

    Moreover, when J(f) is invertible, all entries of J(f)1 are positive; if G is strongly connected, they are strictly positive.

    Equilibria of the dynamics ˙x=J(f)x+ϕ are the solutions of the linear system of equations J(f)x=ϕ. By Proposition 11, if all vertices v have a path to some excretions, then there is a unique equilibrium ˉx=J(f)1ϕ, and moreover all entries of ˉx are positive. To study the general case, where vertices might or might not have a path to some excretions, it is convenient to partition the system into two subsystems, as follows. Partition the vertex set as V=V1V2, with V1 the set of vV such that there is a path from v to vn+1. Without loss of generality, we can re-label vertices in V so as to have vertices v1,,vrV1 and vr+1,,vnV2. According to this decomposition, partition the vector x into two blocks x1 corresponding to V1 and x2 corresponding to V2, and similarly partition ϕ as ϕ1, ϕ2. Notice that there is no edge from V2 to V1, since an edge (w,v) with vV1 implies that there is a path from v to vn+1 and also a path from w to vn+1. Hence, we have:

    J(f)=[J10J21J2],

    and

    {˙x1=J1x1+ϕ1˙x2=J2x2+J21x1+ϕ2. (16)

    The first subsystem is called the reduced system, and its evolution is not affected by the second subsystem. The matrix J1 is equal to the matrix J(f) of the graph ˜H obtained from ˜G as follows: for any edge (v,w) with vV1, wV2, remove the edge (v,w) and add the edge (v,vn+1) with flow fv,vn+1=fv,w; then remove all vertices in V2 and all corresponding edges. By definition of V1, all vertices of ˜H have a path to vn+1, and hence, by Proposition 11, J1 is invertible and Hurwitz stable, and J11 has positive entries.

    For the second subsystem, it is easy to see that J2 is the matrix J(f) of the subgraph K of G corresponding to vertices in V2. Hence, L=JT2 is a Laplacian matrix, and by Proposition 9 its nullspace is generated by vectors π1,,πk, where G1,,Gk are the terminal components in K (same as the terminal components with no excretions in G), and each vector πi has strictly positive entries corresponding to vertices in Gi and is 0 elsewhere. Choosing each πi so that its entries sum to 1, the restriction of πi to Gi is the stationary distribution of the corresponding Markov chain restricted to Gi, namely the Markov chain whose generator matrix is the transpose of the submatrix of J(f) corresponding to Gi.

    These remarks, together with standard tools of analysis of linear dynamical systems, lead to the following proposition:

    Proposition 12 ([4], Theorem 9.13). Consider a weakly connected linear LIFE system with positive flows on all edges of ˜G. From any positive initial condition, the reduced system (namely the subsystem connected to X) converges to its unique equilibrium with positive entries ˉx1=J11ϕ1.

    If there are some vertices not connected to X, then:

    If there is a terminal component GT with no excretions such that there is a path from I to GT, then the mass in GT grows unbounded, and hence also limtx2(t)=+ and limtx(t)=+.

    If for all terminal components with no excretions G1,,Gk there is no path from I to Gi, then the mass of the system remains bounded, and moreover limtx2(t) is some equilibrium point ˉx2 (depending on the initial condition), such that all entries of ˉx2 corresponding to non-terminal components are 0, while the restriction of ˉx2 to a terminal component Gi is proportional to the stationary distribution of the Markov chain with generator matrix Q equal to the transpose of the submatrix of J(f) corresponding to Gi.

    Example 4.1 In this example a simulation is used to verify the calculated equilibrium found using Proposition 12. Here we use the RCT network shown in 5 and calculate the equilibrium. Because every vertex of the RCT is connected to X we have that ˙x=J(f)x+ϕ with,

    J(f)=(f(v1,v4)000000f(v2v4)000000f(v3v4)000f(v1v4)f(v2v4)f(v3v4)f(v4v5)f(v4v6)00000f(v4v5)f(v5v6)0000f(v4v6)f(v5v6)f(v6vn+1)),ϕ=(f(v0,v1)f(v0,v2)f(v0,v3)000). (17)

    The vector of fluxes and initial metabolite values were randomized to obtain

    f=(f(v0,v1)f(v0,v2)f(v0,v3)f(v1,v4)f(v2,v4)f(v3,v4)f(v4,v5)f(v5,v6)f(v5,v6)f(v6,vn+1))=(0.27290.03720.67330.42960.45170.60990.05940.31580.77270.6964),x0=(0.12530.13020.09240.00780.42310.6556). (18)

    Using these values the equilibrium was calculated to be ˉx=(0.6354,0.0824,1.1040,2.6211,0.2015,1.4122). The RCT network with the randomized initial values was simulated for 50 hours and the simulation results closely matched ˉx. The simulation results for the first 25 hours are shown in figure 6.

    Figure 6. 

    The trajectories of the values of metabolites over 25 hours

    .

    In this section, we focus on special LIFE systems. We first recall some interesting results from [16] valid under Assumption (B), and then we exploit them together with the spectral properties of J(f) described in Sect. 4.1 in order to fully characterize equilibria and convergence of special LIFE systems under Assumption (C). We refer the reader to [14] for some results valid under great generality, with assumptions less restrictive than Assumption (B). As for linear systems, we first notice that weakly connected components correspond to subsystems which have no influence on each other, and hence can be studied separately; if there are any weakly connected components with no intake nor excretion, they are subsystems with constant total mass.

    Proposition 13. ([16,Theorem 6]) Consider the special LIFE system under Assumption (B) with no intakes nor excretions. The following properties hold:

    The total mass of the system m=vVxv is constant in time.

    From any positive initial condition x(0), the system tends to the equilibrium set.

    Moreover, if there is a unique terminal component, then there exists a unique equilibrium with positive entries with the same mass as the initial mass, and the system converges to it.

    For the general case with intakes and/or excretions, the following result holds on the asymptotic behavior of the dynamics.

    Proposition 14. ([16,Theorems 2 and 3]) For special LIFE systems under Assumption (B), with a positive initial condition,

    Trajectories are bounded if and only if there exists an equilibrium with positive entries;

    If trajectories are bounded, then they approach an equilibrium set for t, and if moreover the equilibrium set consists of isolated points, then they converge to some equilibrium.

    A caveat reported in [16] is that the second item in Proposition 14 does not rule out the possibility to have non-periodic oscillatory trajectories that approach the equilibrium set lying outside it, and approaching it rotating infinitely many times; this can happen in the case of a connected compact equilibrium set.

    Remark 2. Recall that, by Proposition 2, if there exists a vertex vV such that there is a path from I to v and no path from v to X, then there exists no equilibrium. By the first item of Proposition 14, this further implies that all trajectories are unbounded.

    The following results concern the existence and uniqueness of equilibria.

    Proposition 15. ([16,Theorems 4 and 5]) For special LIFE systems under Assumption (B), the following holds.

    There exists an equilibrium with positive entries for arbitrary constant intakes if and only if for all vV there is a path to X such that all edges in the path have limxvHe(xv)=+;

    If there exists an equilibrium with positive entries, and if there exists a path from all vV to X, then the equilibrium is unique.

    Notice that in the first item, the 'only if' part is true only when we require existence of equilibria with positive entries for completely arbitrary intakes: arbitrary intake set I and arbitrary positive values of the corresponding fluxes; this strong condition is not necessary to have an equilibrium with positive entries for a given set I and a given value of the intake fluxes. Also notice that this statement considers fixed fluxes fe, differently from Proposition 5, where x is fixed and fluxes fe are allowed to vary (except for intake fluxes); the crucial difference is that the product He(x)fe can be made arbitrarily large in the context of Proposition 5, even in the case where He(x) is bounded.

    A remark about the second item is that, by Proposition 14, the existence and uniqueness of the equilibrium with positive entries further implies that all trajectories converge to such equilibrium.

    To apply the first item of Proposition 14 or the second item of Proposition 15, one needs to already have the knowledge about existence of an equilibrium with positive entries. This happens for example in the case where one starts by fixing a desired equilibrium with positive entries ˉx, and then applies the extreme pathways technique in order to design suitable fluxes ensuring that ˉx is an equilibrium of the system. Under some assumptions on the graph, the above propositions then ensure uniqueness of the equilibrium, and its global asymptotic stability.

    In the remainder of this section, we consider Assumption (C). In this case, recall that the dynamics can be re-written as ˙x=J(f)h(x)+ϕ. Different from linear systems, h(x) can contain non-linearities, but the matrix J(f) has the same definition as for linear systems, so that its spectral properties are described by Propositions 9 and 11. Also notice that equilibria, i.e. solutions of J(f)h(x)=ϕ, can be found by solving J(f)h=ϕ, where h is an unknown vector in (R+)n, and then solving h(x)=h. The latter is equivalent to solving Hvi(xvi)=hi, for i=1,,n.

    We start by studying the case with no intakes nor excretions. In this case, by Proposition 9, J(f)h=0 means hN(J(f)), where the nullspace N(J(f)) is generated by π1,,πk, associated with terminal components G1,,Gk. Without loss of generality, we re-label vertices so that the first n1 vertices encompass the first terminal component G1, the following n2 vertices encompass the second terminal component G2, and so on, up to the last nk vertices encompass the last terminal component Gk, and finally the remaining vertices are not in any terminal component (say there are n0 of them). Denote the corresponding subblocks of vector h as h(1),,h(k) and h(0) for the non-terminal ones, and similarly define x(1),,x(k) and x(0) for vector x. By Proposition 9, the nullspace of J(f) is N(J(f))={h(R+)n:h(j)=αj˜πj,αjR and h(0)=0, for j{1,,k}}, where ˜πi is the restriction of π to the component Gi, i.e. is a vector of size ni with strictly positive entries, representing the stationary distribution of the Markov chain whose generator is the transpose of the restriction of J(f) to Gi. Now we need to characterize the set of equilibria with positive entries ˉX:={x(R+)n such that h(x)N(J(f))}. Recall that Hv(xv) are strictly increasing functions, being 0 when xv=0. Denote by Rv the range of Hv (for xv0), notice that either Rv=[0,hmax), or Rv=[0,+). Now denote by Hj the set of vectors αj˜πj such that αj0 and [αj˜πj]vRv for all v in Gj. Then denote by H1jh(j) the vector obtained from h(j)Hj by applying entry-wise the inverse functions H1v. Finally we obtain ˉX={x(R+)n such that x(j)=H1jh(j),h(j)Hj and h(0)=0}.

    Having characterized the set of equilibria with positive entries, now recall that Proposition 13 applies, and trajectories remain bounded, with total mass constant in time, and approach the above-described equilibrium set. We focus on the case with intakes and/or excretions. The definition of the reduced system and the partitioning in two subsystems introduced for linear systems applies also to special LIFE systems under Assumption (C), with the only difference that now the dynamics are non-linear:

    {˙x1=J1h1(x1)+ϕ1˙x2=J2h2(x2)+J21h1(x1)+ϕ2.

    vectors h1(x1) and h2(x2) having replaced x1 and x2 in (16).

    Hence, one can obtain the following analogous of Proposition 12.

    Proposition 16. Consider a weakly connected special LIFE system satisfying Assumption (C). Denote by Rv the range of the function Hv(xv) (for xv0) and define ˉh=J11ϕ1. If ˉhvRv for all vV1, then from any initial condition, the reduced system converges to its unique equilibrium ˉx1 defined by [ˉx1]v=H1v(ˉhv). Otherwise, the system has no equilibrium, and limtx(t)=+.

    In the case where x1(t) converges to ˉx1, if there are some vertices not connected to X, then:

    If there exists a terminal component GT with no excretions such that there is a path from I to GT, then the mass in GT grows unbounded, and hence also limtx2(t)=+ and limtx(t)=+;

    If for all terminal components with no excretions G1,,Gk there is no path from I to Gi, then the mass of the system remains bounded, and moreover entries of x2(t) corresponding to non-terminal components converge to zero, while the restriction of ˉx2 to a terminal component Gi approaches the equilibrium set constructed as follows. The restriction of J2 to vertices in Gi has a nullspace generated by a single positive vector ˜πi; let Hi denote the subset of such nullspace given by vectors h=α˜πi such that hvRv for all vertices v of Gi; the equilibrium set is given by x's such that there exist hHi verifying xv=H1v(hi) for all vertices v of Gi.

    Proof. Recall that x1 is an equilibrium for ˙x1=J1h1(x1)+ϕ1 if and only if J1h1(x1)=ϕ1. By Proposition 11, J1 is invertible and J11 has positive entries; define ˉh to be the unique solution to J1h1=ϕ1. Notice that ˉh has all positive entries. If all entries of ˉh are within the range of the corresponding function Hv(xv), then we have a unique equilibrium with positive entries ˉx1 obtained as in the statement of the proposition, and hence by Proposition 14 x1(t) converges to ˉx. Otherwise, there exists no equilibrium, and hence, by Proposition 14, the mass of system grows indefinitely.

    If there is a terminal component with no excretions but connected to some intakes, then by Proposition 12 there is no equilibrium, and hence by Proposition 14 the mass of system grows indefinitely. If all terminal components with no excretions are not connected to intakes, then the equilibrium set is obtained from the properties of the nullspace of J2, given by Proposition 9. Moreover, Proposition 14 ensures that trajectories remain bounded and approach the equilibrium set.

    In this section, we shortly recall the zero-deficiency theory, and compare it with our results on equilibria of LIFE systems. Zero-deficiency theory arises in the literature on chemical reaction networks, a seminal paper is [8]. Our short overview is based on [7].

    A free closed chemical reaction network with m reactions between p complexes involving n species can be described by ˙x=SΓR(x) where SMn×p, ΓMp×m is the incidence matrix of the network, and R(x) a column vector of size m. The deficiency of the chemical reaction network is the difference of the dimensions of the nullspaces of SΓ and Γ: δ=dimker(SΓ)dimker(Γ), which is equivalent to δ=rank(Γ)rank(SΓ) and to δ=prank(SΓ), where is the number of weakly connected components. The first equivalence is due to rank-nullity Theorem, and the second to the fact that rank(Γ)=p, since Γ is an incidence matrix (Proposition 4.3 of [2]). The reaction rate vector R(x) is governed by the 'mass-action kinetics': R(x)=KΨ(x), where KMm×p and Ψ(x) is a vector of size p, defined as follows: Kej=ke if complex j is the reactant complex of reaction e, Kij=0 otherwise; Ψi(x)=nj=1xsjij. Hence, the dynamics can be equivalently re-written as ˙x=SΓKΨ(x). The following results hold:

    Proposition 17 (Zero-deficiency Theorem). Consider a free closed chemical reaction network with mass-action kinetics. If its deficiency is zero, then: there exists an equilibrium with strictly positive entries if and only if the system is weakly reversible (i.e. each weakly connected component is also strongly connected).

    Moreover, this strictly positive equilibrium is unique in each stoichiometric class (i.e. each weakly connected component has a space of equilibria which has dimension one, so that its equilibrium is unique up to a multiplicative constant representing the total mass in the component), and it is locally asymptotically stable.

    Notice the close resemblance with the results for linear LIFE systems in the case with no intakes nor excretions, described in Propositions 9 and 10. It turns out that a particular class of closed free chemical reaction networks with mass-action kinetics and zero deficiency exactly coincides with linear LIFE systems with no intakes nor excretions. Indeed, let p=n, let S be the identity matrix of size n, let K be defined by Kej=fe if edge e is of the form e=(vj,w) and Kej=0 otherwise; moreover choose exponents sji to be zero and ones so that Ψ(x)=x. This gives exactly a linear LIFE system with no intakes nor excretions. It is immediate to see that the deficiency is zero, since S=I implies rank(SΓ)=rank(Γ).

    In the more general case, free closed chemical reaction networks with mass-action kinetics and zero deficiency are a class of LIFE systems with different assumptions than those considered in this paper, e.g., they might not even satisfy Assumption (A), due to the presence of a rectangular matrix S left-multiplying Γ, and often do not satisfy Assumption (B), due to the dependence of entries of Ψ(x) on various entries of x.

    For general LIFE systems (Assumption (A)), we show the existence of positive solutions. Stability criteria have been shown for the structure of a graph associated to LIFE systems. We show that further conclusions cannot be drawn in the general case for LIFE systems because arbitrary dynamics may be defined.

    For the problem of understanding equilibria for a fix set of metabolites, we show that the rank of the stoichiometric matrix (and the nullspace of this matrix) are determined by structural properties of the graph associated to the system, discuss the effect of the associated graph having intake vertices and excretion vertices on the rank of the stoichiometric matrix, and give necessary conditions on the structure of the graph associated to LIFE systems for the existence of equilibria. More biologically relevant equilibria of LIFE systems are those with all positive metabolites and fluxes and we prove that the extreme pathways method for finding a positive basis describes all such equilibria.

    The field of network flows contributes to LIFE systems by way of the min cut max flow theorem. We show the capacity of edges in network flows relates to saturation of functions corresponding to edges of LIFE systems. The method of extreme pathways to calculate a positive basis is proven to include the entire intersection of the nullspace with the positive orthant. This basis is essential to describing equilibria of LIFE systems.

    Equilibria and asymptotic behavior of metabolites for fixed fluxes are studied. We show that under stricter assumptions (Assumption (C)) we determine the eigenvalues of the jacobian of the LIFE system as well as structure of the graph which admit certain equilibria. We analyze the case with intakes and excretions versus the cases without intakes or excretions.

    Equilibria for LIFE systems with terminal components exist, but these equilibria are dependent upon initial mass of the system. Furthermore, conditions are given for LIFE systems to tend to equilibrium from non-negative initial conditions. LIFE systems under Assumption (B) have bounded solutions if and only if there exits a non-negative equilibrium. Then for Assumption (B), we give conditions on the structure of the associated graph for which nontrivial equilibria exist.

    We show that LIFE systems with no intakes or excretions have "zero deficiency." The rank of the stoichiometric matrix (defined in this work) gives information about the structure of the associated graph, specifically the connectivity of the graph and the existence of strongly connected components. Zero deficiency theory also provides information about the existence of equilibria with respect to this structure.

    The structural conditions of graphs discussed in this work have implications on metabolic networks. We have included results about networks which are considered non-biological for the sake of completeness. These results allow a clear picture of the structure of metabolic systems which are capable of admitting a biologically relevant equilibrium. With this clear picture, we are able to determine metabolic networks for which it is most advantageous to analyze using Linear-In-Flux-Expressions.

    This Appendix contains the calculations for examples 2 and 3.2 Example 3.1 illustrates finding the extreme pathways of RCT using the algorithm described in Section 3.2. Example 3.2 compares a standard basis and positive basis for RCT and shows that when only the positive orthant is considered they describe the same space.

    Example 3.1 continued. As shown previously, the stoichiometric matrix for the RCT is given by:

    S(x)=(100x10000000100x20000000100x30000000x1x2x3x4x400000000x40x500000000x4x5x6)

    First let V0 be the m×m identity matrix where m is the number of columns in S. Then define C0 as follows, V0ST=C0. To better illustrate the row combinations in the algorithm a column is appended to C0 which labels the rows.

    We have

    (1000000000010000000000100000000001000000000010000000000100000000001000000000010000000000100000000001)ST=(100000R1010000R2001000R3x100x100R40x20x200R500x3x300R6000x4x40R7000x40x4R80000x5x5R900000x6R10)

    The next step is to identify any columns of C0 that have no sources or sinks. For this example there are two columns, column 4 and column 5. Then we modify both V0 and C0 by first copying all the rows with a zero in the entry of the first identified column (column 4), then of the remaining rows add all possible positive combinations of two rows which produce a zero in this column. Note that V0 is 10×10 and C0 is 10×6.

    This gives us

    (100000000001000000000010000000000x4x1001000000x4x10001000000x4x2010000000x4x20010000000x4x3100000000x4x3010000000000100000000001)ST=(100000R1010000R2001000R3x4000x40x4x1R4+R7x40000x4x4x1R4+R80x400x40x4x2R5+R70x4000x4x4x2R5+R800x40x40x4x3R6+R700x400x4x4x3R6+R80000x5x5R900000x6R10)

    Note that V1 is 11×10 and C1 is 11×6.

    Now doing the same for column 5 gives

    (100000000001000000000010000000000x4x10010x4x50000x4x10001000000x4x2010x4x500000x4x20010000000x4x310x4x5000000x4x301000000000001)ST=(100000R1010000R2001000R3x40000x4x4x1R4+R7+x4x5R9x40000x4x4x1R4+R80x4000x4x4x2R5+R7+x4x5R90x4000x4x4x2R5+R800x400x4x4x3R6+R7+x4x5R900x400x4x4x3R6+R800000x6R10)

    Note that V2 is 10×10 and C2 is 10×6.

    The same process is then followed for columns containing sources and sinks. Starting at the rightmost column add all possible positive combinations of two rows that create a 0 entry in this column. Do the same operations to both V2 and C2. Doing this for just the rightmost column we have

    (100000000001000000000010000000000x4x10010x4x5x4x6000x4x100010x4x60000x4x2010x4x5x4x60000x4x20010x4x600000x4x310x4x5x4x600000x4x3010x4x6)ST=(100000R1010000R2001000R3x400000x4x1R4+R7+x4x5R9+x4x6R10x400000x4x1R4+R8+x4x6R100x40000x4x2R5+R7+x4x5R9+x4x6R100x40000x4x2R5+R8+x4x6R1000x4000x4x3R6+R7+x4x5R9+x4x6R1000x4000x4x3R6+R8+x4x6R10)

    V3 is 9×10 and C3 is 9×6. After the third column we have the following,

    (10000000000100000000000x4x10010x4x5x4x6000x4x100010x4x60000x4x2010x4x5x4x60000x4x20010x4x600x400x4x310x4x5x4x600x400x4x3010x4x6)ST=(100000R1010000R2x400000x4x1R4+R7+x4x5R9+x4x6R10x400000x4x1R4+R8+x4x6R100x40000x4x2R5+R7+x4x5R9+x4x6R100x40000x4x2R5+R8+x4x6R10000000x4R3+x4x3R6+R7+x4x5R9+x4x6R10000000x4R3+x4x3R6+R8+x4x6R10)

    V4 is 8×10 and C4 is 8×6. After the second column we have the following,

    (1000000000000x4x10010x4x5x4x6000x4x100010x4x60x400x4x2010x4x5x4x60x400x4x20010x4x600x400x4x310x4x5x4x600x400x4x3010x4x6)ST=(100000R1x400000x4x1R4+R7+x4x5R9+x4x6R10x400000x4x1R4+R8+x4x6R10000000x4R2+x4x2R5+R7+x4x5R9+x4x6R10000000x4R2+x4x2R5+R8+x4x6R10000000x4R3+x4x3R6+R7+x4x5R9+x4x6R10000000x4R3+x4x3R6+R8+x4x6R10)

    V5 is 7×10 and C5 is 7×6. And finally,

    (x400x4x10010x4x5x4x6x400x4x100010x4x60x400x4x2010x4x5x4x60x400x4x20010x4x600x400x4x310x4x5x4x600x400x4x3010x4x6)ST=(000000x4R1+x4x1R4+R7+x4x5R9+x4x6R10000000x4R1+x4x1R4+R8+x4x6R10000000x4R2+x4x2R5+R7+x4x5R9+x4x6R10000000x4R2+x4x2R5+R8+x4x6R10000000x4R3+x4x3R6+R7+x4x5R9+x4x6R10000000x4R3+x4x3R6+R8+x4x6R10)

    V6 is 6×10 and C6 is 6×6. The rows of V6 are the basis for N(S(x)). While the dimN(S(x))=4, there are 6 extreme pathway vectors. These vectors are positive, and positively linearly independent,

    V6=(x400x4x10010x4x5x4x6x400x4x100010x4x60x400x4x2010x4x5x4x60x400x4x20010x4x600x400x4x310x4x5x4x600x400x4x3010x4x6).

    Example 3.2 continued. The positive basis of S was found in example 3.1. This basis has six vectors despite the four dimensional nullspace. Also note that the basis vectors are not linearly independent, but positively linearly independent. The span of the positive basis vectors is shown below,

    b1(x600x6x1000x6x401)+b2(x600x6x100x6x40x6x51)+b3(0x600x6x200x6x401)+b4(0x600x6x20x6x40x6x51)+b5(00x600x6x30x6x401)+b6(00x600x6x3x6x40x6x51). (19)

    Next the positive basis vectors for the nullspace must be intersected with the positive orthant N(S(x))(R+)m. When this span is intersected with the positive orthant it is clear there is only one condition,

    For all i, bi0biR. (20)

    We refer to this span intersected with the positive orthant as B.

    For the stoichiometric matrix S(x) the basis for the null space was found. The span of the four basis vectors is shown below,

    a1(x600x6x1000x6x401)+a2(000000x5x4x5x410)+a3(x30x3x3x1010000)+a4(x2x20x2x1100000). (21)

    To find the intersection of this span with the positive orthant, three conditions on the ai must hold.

    {ai0, for i{1,2,3,4}a1x6a3x3+a4x2,a1x6a2x5. (22)

    We refer to the intersection of this span with the positive orthant as C. We show that under the conditions given, B=C.

    First it is shown that BC by showing that for an arbitrary set of bi0, ai's can be chosen to reach the same vector. Using the following substitutions for ai it is clear that if bi0 the inequalities of (22) are satisfied

    {a1=b1+b2+b3+b4+b5+b6,a2=(b2+b4+b6)x6x5,a3=(b5+b6)x6x3,a4=(b3+b4)x6x2. (23)

    This shows that any vector in B can be represented by vectors in C i.e. BC.

    Now it is shown that CB. For arbitrary a's which satisfy (22), the following substitutions for b are used and the conditions of (20) are checked.

    {b1=a1a3x3x6a4x2x6a2x5x6+b4+b6,b2=a2x5x6b4b6,b3=a4x2x6b4,b5=a3x3x6b6. (24)

    To insure that (20) is satisfied the following inequalities must hold,

    a1+b4+b6a3x3x6+a4x2x6+a2x5x6, (25)
    a2x5x6b4+b6, (26)
    a4x2x6b4, (27)
    a3x3x6b6. (28)

    If choices for b4 and b6 can be found which satisfy these inequalities then we will have CB.

    We have the following two conditions:

    Condition 1. a4x2x6+a3x3x6a2x5x6. Setting b4=a4x2x6 and b6=a3x3x6 immediately satisfies (26), (27), and (28). From (22) we have a1x6a2x5 which means that inequality (25) is satisfied.

    Condition 2. a4x2x6+a3x3x6>a2x5x6.

    Under Condition 2 three cases must be considered.

    Case 1. a4x2x6<a2x5x6. In this case we set b4=a4x2x6 and b6=a2x5x6a4x2x6. This immediately satisfies (26) and (27). Since a4x2x6+a3x3x6>a2x5x6 (28) is satisfied as well. Then from (22) we have a1x6a3x3+a4x2, which means that (25) is satisfied.

    Case 2. a4x2x6>a2x5x6 and a3x3x6<a2x5x6. In this case we set b6=a3x3x6 and b4=a2x5x6a3x3x6, this satisfies (26), (27), (28). And from (22) we have a1x6a3x3+a4x2, which means that (25) is also satisfied.

    Case 3. a4x2x6>a2x5x6 and a3x3x6>a2x5x6. In this case we set b4=b6=12a2x5x6 which satisfies (26), (27), (28). Then from (22) we have a1x6a3x3+a4x2, which means that (25) is also satisfied.

    Since the as were arbitrary and bs are found which satisfy (20) this gives us that CB as desired.

    The authors acknowledge the support of the Joseph and Loretta Lopez Chair Professorship endowment, Sanofi via the project "Optimization and Simulation Approaches or Systems Pharmacology in the Pharmaceutical Industry" and the NSF Grant # 1107444 KI-Net "Kinetic description of emerging challenges in multiscale problems of natural sciences".



    [1] C. T. Aage, J. N. Salunke, Fixed points of weak contractions in cone metric spaces, Ann. Funct. Anal., 2 (2011), 71.
    [2] M. Abbas, G. Jungck, Common fixed point results for noncommuting mappings without continuity in cone metric spaces, J. Math. Anal. Appl., 341 (2008), 416-420.
    [3] A. N. Abdou, M. A. Khamsi, Fixed point results of pointwise contractions in modular metric spaces, Fixed Point Theory Appl., (2013), 163.
    [4] S. Aleksic, Z. Kadelburg, Z. D. Mitrovic, et al. A new survey: Cone metric spaces, J. Int. Math. Virtual Institute., 9 (2019), 93-121.
    [5] H. Baghani, M. E. Gordji, M. Ramezani, Orthogonal sets: The axiom of choice and proof of a fixed point theorem, J. Fixed Point Theory Appl., 18 (2016), 465-477. doi: 10.1007/s11784-016-0297-9
    [6] S. Banach, Sur les operations dans les ensembles abstrits et leur applications aux equations integrals, Fund. Math., 3 (1922), 133-181.
    [7] L. E. Brouwer, Uber Abbildung von Mannigfaltigkeiten, Math Ann., 71 (1912), 97-115.
    [8] V. V. Chistyakov, Modular metric spaces, I: basic concepts, Nonlinear Anal., 72 (2010), 1-14.
    [9] V. V. Chistyakov, Modular metric spaces, II: application to superposition operators, Nonlinear Anal., 72 (2010), 15-30.
    [10] S. H. Cho, J. S. Bae, Fixed point theorems for multivalued maps in cone metric spaces, Fixed Point Theory Appl., 2011 (2011), Article number: 87.
    [11] L. Ciric, Some Recent Results in Metrical Fixed Point Theory, University of Belgrade, Beograd, 2003.
    [12] G. Deng, H. Huang, M. Cvetkovic, et al. Cone valued measure of noncompactness and related fixed point theorems, Bull. Int. Math. Virtual Inst., 8 (2018), 233-243.
    [13] M. Eshaghi Gordji, H Habibi, Fixed point theory in generalized orthogonal metric space, J. Linear Topol. Algeb., 6 (2017), 251-260.
    [14] M. Eshaghi Gordji, M. Ramezani, M. De La Sen, et al. On orthogonal sets and Banach fixed point theorem, Fixed Point Theory., 18 (2017), 569-578. doi: 10.24193/fpt-ro.2017.2.45
    [15] L. G. Huang, X. Zhang, Cone metric spaces and fixed point theorems of contractive mappings, J. Math. Anal. Appl., 332 (2007), 1468-1476.
    [16] D. Ilic, V. Rakocevic, Common fixed points for maps on cone metric space, J. Math. Anal. Appl., 341 (2008), 876-882. doi: 10.1016/j.jmaa.2007.10.065
    [17] S. Kakutani, A generalization of Brouwer fixed point theorem, Duke Math. J., 8 (1941), 457-459. doi: 10.1215/S0012-7094-41-00838-4
    [18] M. A. Khamsi, W. K. Kozlowski, S. Reich, Fixed point theory in modular function spaces, Nonlinear Anal., 14 (1990), 935-953.
    [19] M. A. Khamsi, N. Hussain, KKM mappings in metric type spaces, Nonlinear Anal., 73 (2010), 3123-3129.
    [20] W. Kirk, N. Shahzad, Fixed Point Theory in Disatnce Spaces, Springer International Publishing Switzeralan, 2014.
    [21] N. Mehmood, A. Al Rawashdeh, S. Radenovic, New fixed point results for E-metric spaces, Positivity, 23 (2019), 1101-1111.
    [22] S. B. Nadler, Multi-valued contraction mappings, Pacific J. Mafh., 30 (1969), 475-488.
    [23] F. J. Nash, Equilibrium points in n-person game, Proc. Natl. Acad. Sci. USA., 36 (1950), 48-49. doi: 10.1073/pnas.36.1.48
    [24] J. J. Nieto, R. Rodri Guez-Lo'Pez, Contractive Mapping Theorems in Partially Ordered Sets and Applications to Ordinary Differential Equations, Order 22 (2005), 223-239.
    [25] M. A. Noor, An iterative algorithm for variational inequalities, J. Math. Anal. Appl., 158 (1991), 448-455.
    [26] Z. Pales, I-R. Petre, Iterative fixed point theorems in E-metric spaces, Acta Math. Hung., 140 (2013), 134-144.
    [27] E. De Pascale, G. Marino, P. Pietramala, The use of the E-metric spaces in the search for fixed points, Matematiche, 48 (1993), 367-376.
    [28] M. Ramezani, Orthogonal metric space and convex contractions, Int. J. Nonlinear Anal. Appl., 6 (2015), 127-132.
    [29] M. Ramezani, H. Baghani, Contractive gauge functions in strongly orthogonal metric spaces, Int. J. Nonlinear Anal. Appl., 8 (2017), 23-28.
    [30] A. C. M. Ran, M. C. B. Reurings, A fixed point theorem in partially ordered sets and some applications to matrix equations, Proc. Am. Math. Soc., 132 (2004), 1435-1443.
    [31] A. Al-Rawashdeh, W. Shatanawi, M. Khandaqji, Normed ordered and E-metric spaces, Int. J. Math. Sci., 2012 (2012), ID 272137.
    [32] W. Rudin, Principles of mathematical analysis, thired edition, McGraw-Hill, Inc, 1976.
    [33] V. Todorcevic, Harmonic Quasiconformal Mappings and Hyperbolic Type Metrics, Springer Nature Switzerland AG, 2019.
    [34] D. Turkoglu, M. Abuloha, A. bdeljawad T: KKM mappings in cone metric spaces and some fixed point theorems, Nonlinear Anal., 72 (2010), 348-353.
  • This article has been cited by:

    1. Christopher Denaro, Nathaniel J. Merrill, Sean T. McQuade, Logan Reed, Chanchala Kaddi, Karim Azer, Benedetto Piccoli, A pipeline for testing drug mechanism of action and combination therapies: From microarray data to simulations via Linear-In-Flux-Expressions, 2023, 00255564, 108983, 10.1016/j.mbs.2023.108983
    2. Karim Azer, Irina Leaf, Systems biology platform for efficient development and translation of multitargeted therapeutics, 2023, 3, 2674-0702, 10.3389/fsysb.2023.1229532
  • Reader Comments
  • © 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(8216) PDF downloads(1133) Cited by(22)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog