Control of Metabolic Systems Modeled with Timed Continuous Petri Nets Roberto Ross-León1 , Antonio Ramirez-Treviño1 , José Alejandro Morales2 , and Javier Ruiz-León1 1 Centro de Investigaciones y Estudios Avanzados del I.P.N. Unidad Guadalajara {rross,art,jruiz}@gdl.cinvestav.mx 2 Centro Universitario de Ciencias Exactas e Ingenierías, Universidad de Guadalajara alejandro.morales@cucei.udg.mx Abstract. This paper is concerned with the control problem of biolog- ical systems modeled with Timed Continuous Petri Nets under in…nite server semantics. This work introduces two main contributions. The …rst one is a bottom-up modeling methodology that uses T CP N to represent cell metabolism. The second contribution is the control wich solves the Regulation Control Problem (RCP ) (to reach a required state and maintain it). The control is based on a Lyapunov criterion that ensures reaching the required state. Key words: Cell metabolome, Petri nets, Controllability, Stability. This work was supported by project No. 23-2008-335, COETCYJAL-CEMUE, Mexico. R. Ross-León was supported by CONACYT, grant No. 13527. 1 Introduction Petri nets P N [1], [2], [3] are a formal paradigm for modelling and analysis of sys- tems that can be seen as discrete dynamical systems. Unfortunately, due to state explosion problem, most of the analysis techniques cannot be applied in heavy marked Petri nets. In order to overcome this problem, the Petri net community developed the Timed Continuous Petri Nets (T CP N ) [4], [5], a relaxation of the Petri Nets where the marking becomes continuous and the state equation is represented by a positive, bounded set of linear di¤erential equations. The main T CP N characteristics such as the nice pictorially representation, the mathematical background, the synchronization of several products to start an activity and the representation of causal relationship make T CP N amenable to represent biochemical reactions and cell metabolism. In fact T CP N mark- ing captures the concentration of molecular species while di¤erential equations together with the …ring vectors represent the reaction velocity and the graph cap- tures the metabolic pathways. The entire T CP N captures the cell metabolome. Several works model [6], [7], analyse [8], [9] and control [10], [11] metabolic pathways. Most of them deal with pseudo-steady states of the biochemical reac- tion dynamic. Nowadays, the scienti…c community is exploring the use of P N and Recent Advances in Petri Nets and Concurrency, S. Donatelli, J. Kleijn, R.J. Machado, J.M. Fernandes (eds.), CEUR Workshop Proceedings, ISSN 1613-0073, Jan/2012, pp. 87–102. 88 Petri Nets & Concurrency Ross-León et al. 2 their extensions [12], [13] to model biological systems since the former are able to capture the compounds ‡ow, the reaction velocity, the enabling/inhibiting reactions and both the transitory and steady states of reaction dynamic into a single formalism. This work is concerned on how to model the entire metabolome with T CP N . It proposes a bottom-up modeling methodology where biochemical reactions are modeled through elementary modules, and shows how these modules are merged to form metabolic pathways, and at the end the cell metabolism. The resulting model captures both, the transitory and steady state metabolome dynamics. It is worth noticing that the derived T CP N model condenses several particular be- haviors represented by the set of di¤erential equations generated by the T CP N itself. For instance, a single transition with four input places (a reaction needing four substrates) generates a set of four possible di¤erential equations while two transitions with four input places each will generate a set of sixteen possible di¤erential equations. Therefore highly complex behaviors emerging from few compounds interacting can be captured by T CP N . This work also presents the control problem of reaching a required state (marking) representing a certain metabolite concentration. In order to solve this problem, an error equation is stated and stabilized using a Lyapunov approach. The solution is the reaction rate vector which is greater or equal to zero and lower or equal to the maximum settled by the kinetics of Michaelis-Menten for the current enzyme concentration. Thus, if a solution exists, it could be imple- mented in vivo by directed genetic mutation, knock-in (or knock-out) strategies or pharmacological e¤ects. Present paper is organized as follows. Section 2 gives T CP N basic de…ni- tions, controllability and cell metabolic concepts. Next section introduces the proposed metabolome modeling methodology. Section 4 presents the problem of reaching a required state and synthesizes Lyapunov like transition ‡ow for solving this problem. Following section presents an illustrative example to show the performance of the computed control law. In the last section the conclusions and future work are presented. 2 Basic De…nitions This section presents brie‡y the basic concepts related with P N , Continuous P N and T CP N . An interested reader can review [3], [14], [15] and [16] for further information. At the end of this section a useful form of the state equation for T CP N under in…nite server semantics is presented. 2.1 Petri Net concepts De…nition 1. A Continuous Petri Net (ContP N ) system is a pair (N; m0 ), where N = (P; T; P re; P ost) is a Petri net structure (P N ) and m0 2 fR+ [0gjP j is the initial marking. P = fp1 ; :::; pn g and T = ft1 ; :::; tk g are …nite sets of jP j jT j elements named places and transitions, respectively. P re; P ost 2 fN [ 0g Control of metabolic systems Petri Nets & Concurrency – 89 3 are the Pre and Post incidence matrices, respectively, where P re[i; j], P ost[i; j] represent the weights of the arcs from pi to tj and from tj to pi , respectively. The Incidence matrix denoted by C is de…ned by C = P ost P re: Each place pi has a marking denoted by mi 2 fR+ [ 0g. The set ti = fpj j P re[j; i] > 0g ; (ti = fpj j P ost[j; i] > 0g) is the preset (postset) of ti : Sim- ilarly the set pi = ftj j P ost[i; j] > 0g ; (pi = ftj j P ost[i; j] > 0g) is the preset (postset) of pi . A transition tj 2 T is enabled at marking m i¤ 8pi 2 tj , mi > 0. Its enabling degree is: mi enab(tj ; m) = min (1) pi 2 tj P re [i; j] and it is said that mi constraints the …ring of tj . Equation (1) denotes the maximum amount that tj can be …red at marking m; indeed tj can …re in any real amount ; where 0 < < enab(tj ; m) leading to a new marking m0 = m + C[ ; j]. If m is reachable from m0 through a …nite sequence of enabled transitions, then m can be computed with the equation: m = m0 + C (2) jT j named the ContP N state equation, where 2 fR+ [ 0g is the …ring count vector, i.e., j is the cumulative amount of …ring of tj in the sequence . The set of all reachable markings from m0 is called the reachability space and it is denoted by RS (N; m0 ). In the case of a ContP N system, RS (N; m0 ) is a convex set [17]. De…nition 2. A contP N is bounded when every place is bounded (8p 2 P; 9bp 2 R with m [p] bp at every reachable marking m). It is live when every transi- tion is live (it can ultimately occur from every reachable marking). Liveness is extended to lim-live when in…nitely long sequence can be …red. A transition t is non lim-live i¤ a sequence of successively reachable markings exists which con- verge to a marking such that none of its successors enables a transition t. 2.2 Timed continuous Petri nets De…nition 3. A timed ContP N is the 3-tuple T CP N = (N; ; m0 ) ; where N is a ContP N , : T ! fR+ gjT j is a function that associates a maximum …ring rate to each transition, and m0 is the initial marking of the net N . The state equation of a T CP N is m( ) = Cf ( ) (3) where f ( ) = ( ) And under the in…nite server semantics, the ‡ow of transition tj is given by fj ( ) = j enab(tj ; m( )) (4) 90 Petri Nets & Concurrency Ross-León et al. 4 where j represents the maximum …ring rate of transition tj . Notice that T CP N under in…nite server semantics is a piecewise linear system (a class of hybrid systems) due to the minimum operator that appears in the enabling function of the ‡ow de…nition. De…nition 4. A con…guration of a T CP N at m is a set of (p; t) arcs describing the e¤ ective ‡ow of all transitions. 1 (m) [i; j] = P re[i;j] if pi is constraining tj (5) 0 otherwise De…nition 5. The maximum …ring rate matrix is denoted by = diag 1; : : : ; jT j : (6) According to previous notation, the state equation and the ‡ow vector are described by: m=C (m) m (7) f= (m) m The only action that can be applied to a T CP N system is to slow down the …ring ‡ow. The forced ‡ow of a controlled transition ti becomes fi ui where fi is the ‡ow of the unforced system (i.e. without control) and u is the control action, with 0 ui fi . The controlled state equation is: m=C[ (m) m u] (8) 0 ui [ (m) m]i (9) In order to obtain a simpli…ed version of the state equation, the input vector u is rewritten as u = Iu (m) m, where Iu = diag Iu1 ; : : : ; IujT j and 0 Iui 1. Then the matrix Ic = I Iu is constructed and the controlled state equation can be rewritten as: m = CIc (m) m (10) Notice that 0 Ici 1. 2.3 Controllability The classical linear systems de…nition of controllability cannot be applied to T CP N systems because the required hypothesis are not ful…lled, that is, the input should be unbounded and the state space should be RjP j . The next de…n- itions are taken from [18]. De…nition 6. Let N be net of a T CP N . The structural admissible states set is jP j de…ned as SASS (N ) = fR+ [ f0gg (all inital markings that can be imposed to a net). Let B be the base of the left annuller of the incidence matrix C. The equivalence relation : SASS (N ) ! SASS (N ) is de…ned as m1 m2 i¤ B T m1 = B T m2 , 8m1 ; m2 2 SASS(N ). The system admissible states set is the equivalent class of the initial marking Class (m0 ) under . Control of metabolic systems Petri Nets & Concurrency – 91 5 In the sequel, let us denote by int (Class (m0 )) the set of relative interior of Class (m0 ). De…nition 7. Let (N; ; m0 ) be a T CP N system. It is fully controllable with bounded input (BIF C) if there is an input such that for any two markings m1 ; m2 2 Class (m0 ), it is possible to transfer the marking from m1 to m2 in …nite or in…nite time, and the input ful…lls (9) along the trajectory, and is controllable with bounded input (BIC) over S Class (m0 ) if there is an input such that for any two markings m1 ; m2 2 S, it is possible to transfer the mark- ing from m1 to m2 in …nite or in…nite time, and the input ful…lls (9) along the trajectory. De…nition 8. Let (N; ; m0 ) be a T CP N system. Let mr 2 RS (N; m0 ) and 0 Icr [i; i] 1. Then (mr ; Icr ) is an equilibrium point if mr = CIcr (mr ) m = 0. Then, the steady state ‡ow for (mr ; Icr ) is fss (mr ; Icr ) = Icr (mr ) mr . An equilibrium point represents a state in which the system can be main- tained using the de…ned control action. Given an initial marking m0 and a re- quired marking mr , one control problem is to reach mr and then keep it. For a further information about equilibrium points an interested reader can review [19]. 2.4 Cell Metabolism For the wellbeing of an given organism, each cell of that organism must transform the substances available in its surroundings to useful molecules. Such transforma- tions take place as chemical reactions catalyzed by enzymes. In these reactions, a substrate tightly binds non-covalently to its enzyme active site to build an enzyme-substrate complex. At that moment, the enzyme chemically changes the substrate into one or more products and then releases it. The enzyme did not su¤er any irreversible alterations in the process, and now is free to accept a new substrate [20]. There is no limit to the number of possible reactions occurring in nature. Nonetheless, after exhaustive analysis certain general patterns had emerged that became useful to describe several characteristics of biochemical reactions. In the case where a sole substrate becomes a single product, the reaction process is represented by the scheme: E+S ES ! EQ E+Q (11) where E is the enzyme, S is the substrate, ES and EQ are the bound complexes and Q is the product. Typically, the rate of these reactions is settled by the kinetics of Michaelis- Menten [21]. Under this kinetic model, the enzyme and substrate react rapidly to form an enzyme-substrate complex while [S] and [ES] are considered to be at concentration equilibrium (the same applies to [EQ] and [Q]), that is, the rate 92 Petri Nets & Concurrency Ross-León et al. 6 at which ES dissociates into E + S is much faster than the rate at which ES brakes down to EQ. Throughout the present work, we will consider a physiological cellular state where [S] >> [E], which means that [S] [ES] equilibrium will always tend to complex formation. Therefore, ES dissociation rate is irrelevant and Scheme 11 can be abbreviated as follows: E+S !E+Q (12) where the association-dissociation is implicit. In reactions with more than one substrate, binding can occur in di¤erent sequences; for instance, the following scheme represents an enzyme system with two substrates and all the possible sequences: 8 > > ES1 + S2 > > < E + S1 + S2 ES1 S2 ! E + Q1 + ::: + Qn (13) > > > > : ES2 + S1 Frequently the product of an enzyme is the substrate of another reaction and so on, to build a chain of reactions called metabolic pathways represented by M P j = 1j 2j j j j n where i is a reaction (12) or (13) of a pathway j and k uses one or more products of im . Notice that j and m may represent di¤erent pathways. Then a (Cell) Metabolome is CM = M P i M P i is a metabolic pathway , and the purpose of CM is to produce a particular set of metabolites in certain concentrations, essential to that cell. 3 Modelling the Metabolome In order to model the metabolome using T CP N it is necessary to identify how the elements involved in it will be represented. The next table relates the meaning of each element of the T CP N with respect to metabolic reactions. TCPN term Molecular interpretation Place Molecular Species Marking Concentration Transition Reaction Firing Rate Rate of Reaction Arc Weights Stoichiometric Coe¢ cients The bottom-up approach herein proposed to model the metabolome consists of: a) representing reactions, the results of this stage are the elementary modules; b) merging elementary modules, where places of elementary modules represent- ing the same molecular species on the same physical space in the cell will merge Control of metabolic systems Petri Nets & Concurrency – 93 7 Fig. 1. Four elementary modules representing four diferent reactions. into a single place. The results of this stage are pathway modules; and c) merg- ing pathways modules, where places of pathways modules representing the same molecular species on the same cellular space will merge into a single place; the result of this stage is the metabolome model. For stages b and c, any specie being protein-mediated transported into a di¤erent organelle shall be modeled through the same elementary module, representing instead of substrate and product the same molecule in di¤erent spaces. Next section describes these stages. 3.1 Representing Reactions In order to represent each reaction i with T CP N elementary modules repre- senting the Scheme (12) or (13) are constructed. There exists one place pj for each molecular species at the same physical space msj and one transition ti to represent the reaction i . There exists one arc (ps ; ti ) if ps represents a substrate. There exists one arc (ti ; pq ) if pq represents a product. Finally, there exists a self- loop around pe and ti if pe represents an enzyme. The initial marking m0 [pj ] is the concentration of the molecular species msj at time = 0. Associated to transition ti is i representing the rate of reaction. Example 1. Let P 1 + E1 ! P 2 + E1 be the 1 reaction. There is one place for each molecular species (P 1, P 2 and E1), and one transition t1 representing 1 . Finally, arcs are …xed in the way depicted in Figure 1a. Assuming that the substrate concentration will remain higher than the en- zyme concentration (this is an expected behavior of the system), the con‡ict 94 Petri Nets & Concurrency Ross-León et al. 8 Fig. 2. Example of a Pathway Module. between substrate and enzyme can be ignored. Hence, if a system has 2n con…g- urations originated by the n number of enzymes in con‡ict with substrates, all those con…gurations are eliminated because min ([E]; [S]) = [E] for all 0. 3.2 Merging Elementary Modules Let N 1 and N 2 be two elementary modules, then the merging N is such that N = (P; T; P re; P ost) where P = P 1 [ P 2 , T = T 1 [ T 2 , P re = P re1 [ P re2 and P ost = P ost1 [ P ost2 . Notice that places representing the same molecular species in the same physical space are merged into a single place. After a merging of elementary modules is made, pathway modules are ob- tained. Example 2. Let N 1 = P 1 ; T 1 ; P re1 ; P ost1 and N 2 = P 2 ; T 2 ; P re2 ; P ost2 be two elementary modules showed in Figure 1a and Figure 1b respectively. Then, the merging is N = (P; T; P re; P ost) where P = P 1 [P 2 = fP 1; :::; P 5; E1; E2g, T = T 1 [ T 2 = ft1 ; t2 g and arcs are …xed in the way depicted in Figure 2, where the merging is showed. 3.3 Merging Pathway Modules Let N 1 be a pathway module and N 2 be a pathway or an elementary module, then the merging N is such that N = (P; T; P re; P ost) where P = P 1 [ P 2 , T = T 1 [ T 2 , P re = P re1 [ P re2 and P ost = P ost1 [ P ost2 . Notice that places representing the same molecular species are merged into a single place. After a merging of pathway modules is made, a metabolic model is obtained. Example 3. Let N 1 = P 1 ; T 1 ; P re1 ; P ost1 be the pathway module showed in Figure 2. Let N 2 = P 2 ; T 2 ; P re2 ; P ost2 and N 3 = P 3 ; T 3 ; P re3 ; P ost3 be two elementary modules showed in Figure 1c and Figure 1d respectively. Then Control of metabolic systems Petri Nets & Concurrency – 95 9 Fig. 3. Metabolic Model. the merging is N = (P; T; P re; P ost), where P = [P i = fP 1; :::; P 5; E1 ; :::; E4 g, T = [T i = ft1 ; t2 ; t3 ; t4 g for i = 1; :::; 4. Arcs are …xed in Figure 3, where the merging is showed. Although obtained metabolic models could be not live, the addition of a virtual transition and arcs going from the last place representing …nal products to the virtual transition and from virtual transition to the places representing initial products with an appropriate virtual reaction velocity will make the metabolic model live. For instance, consider the net of Figure 1a, it is a non-live net, but if we add a virtual input transition tv to the place S and a virtual output transition tv to the place Q the system will gain liveness, see Figure 4. Notice that tv must to be the same transition added to the inital and …nal metabolites, this is because it is necessary to maintain the conservativeness of the matter of the system. This notion is based assuming that each module belongs to a bigger system, therefore, although the real input and output transitions could be not the same, they must have the same …ring ratio. 4 Control Law An important control problem in the metabolic engineering area is to reach a certain metabolome state such that the production of selected metabolites is regulated or particular processes are limited or favored. This problem is captured in T CP N as the reachability problem, i.e. to reach a required state mr from an initial state m0 by means of an appropriate control action. This is formalized as follows. De…nition 9. Let T CP N be a metabolic model. Then the Regulation Control Problem in (mr ; Icr ) (RCP (mr ; Icr ))deals with the computation of a control law 96 Petri Nets & Concurrency Ross-León et al. 10 Fig. 4. Module forced to be live with the addition of a virtual transition tv (in gray). Ic ( ); 0 < f feasible in the T CP N such that m( ss ) = mr and Ic ( ss ) = Icr ; 8 ss f: In order to solve this problem, some extra places are added to the T CP N metabolic model to detect the material passing through transitions. The follow- ing de…nition shows how these places are added to the T CP N: De…nition 10. Let (N; m0 ; ) be a metabolic model T CP N; where N = (P; T; F ). Its extension is de…ned by xT CP N = (Nx ; m0x ; ), where Nx = (P [ Pa ; T; F [ T Fa ); jPa j = jT j, m0x = m0 0jT j ; Fa = f(ti ; pai ) j8ti 2 T and 8pai 2 Pa g : T Then the incidence matrix of xT CP N is Cx = C IjT j . Since x (mx ) = (m) 0jT j jT j ; then the state equation of xT CP N is: " # m CIc (m) m mx = = (14) ma Ic (m) m m(0) = m0 , ma (0) = 0 (15) Remark 1. Notice that the extension has the same dynamic over the metabolic model places and the extra places can only increase its marking. In fact, due to the T CP N is live, then by construction the xT CP N is also live. Then there exists at least one enabled transition. Hence (m) m > 0 (or equivalently ma 0; the zero could be forced by an appropriate control law Ic ). Example 4. An example of an extended net is presented in Figure 5. 4.1 Solution to the RCP (mr ; Icr ) Theorem 1. Let (N; m0 ; ) be a metabolic model T CP N and let xT CP N = (Nx ; m0x ; ) be its extension. If (N; m0 ; ) is BIC over int (Class (m0 )) (notice Control of metabolic systems Petri Nets & Concurrency – 97 11 Fig. 5. Example of an extended net. The gray places are the set of added places Pa and the intermittent arrows are the set of added arcs Fa . that neither the initial marking nor the required marking could be zero compo- nents) and (Icr ; mr ) is an arbitrary equilibrium point, then there exists Ic ( ); 0 f feasible in the T CP N such that m( ss ) = mr , Ic ( ss ) = Icr ; 8 ss f: Proof. If the system is BIC over int (Class (m0 )), then there exists a positive solution r ( ) feasible such that mr = m0 + C r (16) This result was taken from [18]. Thus there exists f ( ) such that: Z f Z f f ( )d = Ic (m) md = r (17) 0 0 From (14): ma ( f ) = r (18) Now, let T ex ( ) = e( ) ea ( ) , 0 f T T (19) e( ) ea ( ) = mr m( ) r ma ( ) and V (ex ) = eTx PL ex (20) where 0 0 PL = (21) 0 IjT j and IjT j is an identity matrix of order jT j jT j. We claim that V (ex ) is a Lyapunov function, i.e it is positive de…nite and its derivative is negative de…nite. 98 Petri Nets & Concurrency Ross-León et al. 12 Since Equation (20) is clearly non-negative de…nite, then we assume that (22) is positive semide…nite, then there exists ex ( 0 ) 6= 0 such that: 0 0 e V (ex ( 0 )) = eT eTa =0 (22) 0 IjT j ea From (22), it is clear that ea ( 0 ) = 0; then from (19) ma ( 0 ) = r : Thus, from (14) and (17) and letting f = 0 : Z f m( )d = C r (23) 0 m( f ) m(0) = C r Thus, from Equation 16 m( f ) = mr ; then ex = 0; a contradiction. Hence V (ex ) is positive de…nite. Now , we prove that V (ex ) is negative de…nite. The di¤erentiate of V (ex ) is: T V (ex ) = 2eTa ea = 2[ r ma ] ma (24) Then, choosing Ic such that: 1 if ma [i] < r [i] Ici = (25) 0 otherwhise we obtain: T [ r ma ] Ic > 0 (26) and T T [ r ma ] Ic = 0 i¤ [ r ma ] = 0 thus V (ex ) < 0 and V (0) = 0: Since ma (0) = 0 and it only increase its value, then Ici is feasible leading from ma (0) = 0 to ma ( f ) = r ; i.e. from m0 to mr . Moreover, assuming mr 2 int (Class (m0 )) it is reached in …nite time because ma [i] = m (min ( ti )) e and m (min ( t)) 6= 0 8 . At f the control law must be switched from Ic ( f ) to Ic ( ss ) = Icr and the regulation control problem is solved. The solution to the RCP (mr ; Icr ) include both, the transitory and steady state control of metabolic systems. It is an improvement to current control solutions, where the biologist and metabolic engineers use stoichiometric non- dynamical approaches such as F BA (Flux Balance Analysis) [22], [23], [24] for the control of metabolic systems. Those are based on a pseudo-stationary state model, represented by the equation: Sv = 0 (27) where S is the matrix of stoichiometry coe¢ cients and the solution v gives the balance of mass for a single equilibrium point at that state (v is the reaction rates vector in a steady state). Control of metabolic systems Petri Nets & Concurrency – 99 13 Fig. 6. Marking evolution of the net of Figure 3 applying RCP (mr ; Icr ) . 5 Illustrative Controlling Metabolic System Example In order to illustrate the RCP (mr ; Icr ) applied to a metabolic system, suppose the pathway module of Figure 2 together with modules c and d of Figure 1 com- prise a cell metabolome. The initial marking used for this example is an arbitrary but physiologically possible initial state for the alleged metabolic model. Example 5. Let the metabolome model of the Figure 3 be the system T CP N = T (N; ; m0 ) with = diag (2; 3; 4; 1) and m0 = 100 80 100 50 70 5 3 2 4 . Let T mr = 95 70 60 65 110 5 3 2 4 be a required marking. We make the extended system like the procedure showed in the Figure 5. We need the solution of r from mr = m0 + C r . Notice that there are a lot of solutions for r but we only focus on the smallest solution of r . For this example the solution is: T r = 30 40 25 0 Solving the RCP (mr ; Icr ) and applying the control (25) to the T CP N = (N; ; m0 ) ; the metabolite concentrations are depicted in Figure 6. The reac- tion velocities (transition ‡ux) is depicted in Figure 7. Notice that from = 0 to = f 4:5 occurs the transitory dynamics, and for > f the steady state is reached. Example 6. In Figure 8 the evolution of marking ma is depicted. When occurs ma [i] = r [i] the control Ici = 0 makes fi = 0 and ma [i] is maintained until = f (ma = r ). Then Ic switches to Icr for the steady state control. 100 Petri Nets & Concurrency Ross-León et al. 14 Fig. 7. Reaction velocities (transition ‡ows) of the controlled metabolic model of the Example 5. Notice that Ic ( ) is applied for 0 < f and Icr ( ) for > f . Fig. 8. Marking ma of the Example 5. Control of metabolic systems Petri Nets & Concurrency – 101 15 6 Conclusions This work presented a model methodology to capture the metabolome behav- ior. It uses a bottom-up approach where each individual biochemical reaction is modeled by elementary T CP N modules and, afterwards, all the modules are merged into a single one to capture the whole metabolome behavior. Such char- acteristic of the methodology makes it simple and easy to use while the complex cell metabolic behavior is captured. This work also presented the problem of reaching a required metabolome state. The solution to this problem are the instantaneous reaction velocities that are realizable in biological system. Present results are being applied to optimize metabolome fermentation in the production of tequila and to biofuels generation. Future perspective involves introduction of stochastic modelling and merging the metabolome with the signaling and genetic networks. References 1. Murata, T.: Petri nets: Properties, analysis and applications. In: Proceedings of IEEE. Volume 77. (1989) 541–580 2. Javier, Esparza; Manuel, S.: Compositional synthesis of live and bounded free choice petri nets. LNCS 571 (1991) 172–187 3. Jorg Dessel, J.E.: Free Choice Petri Nets. Cambridge University Press (1995) 4. Manuel Silva, L.R.: Continuization of timed petri nets: From performance evalu- ation to observation and control. 26th International Conference On Application and Theory of Petri Nets and Other Models of Concurrency (2006) 5. R. David, H.A.: Continuous petri nets. In: In Proc. Of the 8th European Workshop on Application and Theory of Petri Nets. (1987) 275–294 6. Peter J. E. Goss, J.P.: Quantitative modeling of stochastic systems in molecular biology by using stochastic petr nets. In: Natl. Acad. Sci. USA. Volume 95. (1998) 6750–6755 7. Masao Nagasaki, Atsushi Doi, H.M.S.M.: Petri net based description and modeling of biological pathways. Algebraic Biology (2005) 19–31 8. Monika Heiner, David Gilbert, R.D.: Petri nets for systems and syntethic biology. (2008) 215–264 9. David Angeli, Patrick De Leenheer, E.S.: A petri net approach to persistance analy- sis in chemical reaction networks. Biology and Control Theory: Current Challenges (2007) 181–216 10. Riel, N.A.W.V.: Dynamic modelling and analysis of biochemical networks: Mechanism-based models and model-based experiments. Brie…ngs In Bioinfor- matics 7 (2006) 364–374 11. Nevoigt, E.: Progress in metabolic engineering of saccharomyces cerevisiae. Mi- crobiology And Molecular Biology Reviews 72 (2008) 379–412 12. Mor Peleg, Daniel Rubin, R.B.A.: Using petri net tools to study properties and dynamics of biological systems. Journal of the American Medical Informatics As- sociation 12 (2005) 181–199 13. Hiroshi Matsuno, Atsushi Doi, M.N.S.M.: Hybrid petri net representation of gene regulatory network. Paci…c Symposium on Biocomputing 5 (2000) 338–349 102 Petri Nets & Concurrency Ross-León et al. 16 14. Cassandras, C.G.: Discrete Event Systems. Modelling and Performance Analysis. Asken Associates (1993) 15. René David, H.A.: Continuous and hybrid petri nets. Journal of Circuites, Systems and Computers 8 (1998) 159–188 16. R. David, H.A.: Continuous petri nets. Proceedings of the 8th European Workshop on Application and Theory of Petri Nets (1987) 275–294 17. Manuel Silva, L.R.: On ‡uidi…cation of petri net models: From discrete to hybrid and continuous model. Annual Reviews in Control 28 (2004) 253–266 18. C. R. Vázquez, A. Ramírez-Treviño, L.R.M.S.: On controllability of timed contin- uous petri nets. 11th Int. Workshop Hybrid Systems: Computational and Control 4981 (2008) 528–541 19. Cristian Mahulea, Antonio Ramirez, L.R.M.S.: Steady state control, zero valued poles and token conservation law in continuous net systems. Proceedings of the International Workshop on Control of Hybrid and Discrete Event Systems (2005) 20. Laszlo Kurti, B.C.: Strategic Applications of Named Reactions in Organic Synthe- sis. Elsevier Academic Press. USA (2005) 21. Segel, I.H.: Enzyme Kinetics. New York: Wiley-Interscience (1975) 22. J. M. Lee, E. P. Gianchandani, J.A.P.: Flux balance analysis in the era of metabolomics. Briefs in Bioinformatics 7 (2006) 140–150 23. J. S. Edwards, M. Covert, B.P.: Metabolic modelling of microbes: The ‡ux-balance approach. Environmental Microbiology 4 (2002) 133–140 24. Francisco Llaneras, J.P.: Stoichiometric modelling of cell metabolism. Journal of Bioscience and Bioengineering 105 (2008) 1–11