 devoted their efforts to developing a systems-level understanding of biochemical reaction networks. In particular, the study of these chemical reaction networks (CRNs) using their associated graph structure has attracted considerable attention. The work led by Craciun and Feinberg on multistationarity [@Craciun2005; @Craciun2010; @Craciun2006; @CraciunFeinberg2006] and the work led by Mincheva and Roussel on stable oscillations [@Mincheva2011a; @Mincheva2007a; @Mincheva2007] are two particularly influential approaches. For a good overview of the various graph theoretic developments, we direct the interested reader to the review provided in Domijan and Kirkilionis [@Domijan2008].

In this work, we focus on applications of graph theory, mainly the Matrix-Tree Theorem (MTT), for deriving equilibrium solutions (ES) for CRNs that fit within a Laplacian dynamics framework. The MTT-based framework was first applied in a biological context by King and Altman [@King1956] to derive steady state rate equations in enzyme kinetics. This framework was then simplified and summarized into rules (known as *Chou’s graphical rules* [@Lin2013]) by Chou and coworkers [@Chou1989; @Chou1990; @Chou1993]. Chou [@Chou1989] has also extended the framework for non-steady state enzyme-catalyzed systems.

The main disadvantage of Chou’s graphical rules is that they are only applicable if the underlying digraph structure is *strongly connected*, i.e., every vertex is reachable from every other vertex. This issue was solved and extended for general directed graphs (*digraphs*) by Mirzaev and Gunawardena in 2013 [@MirzaevGunawardena2013bmb] and is applicable to the specific class of linear ordinary differential equations (ODEs) known as *Laplacian dynamics*. Systems described by Laplacian dynamics are created using a weakly connected digraph, $G$, with $n$ vertices, with labelled, directed edges, and without self loops. Note that by *weakly connected* we mean that the graph cannot be expressed as the union of two disjoint digraphs. If there is an edge from vertex $j$ to vertex $i$, we label it with $e_{ij}>0$, and with $e_{ij}=0$ if there is no such edge. [^3]
devoted their efforts to developing a systems - degree reason of biochemical reaction networks. In particular, the report of these chemical reaction networks (CRNs) use their consort graph structure has attracted considerable attention. The work precede by Craciun and Feinberg on multistationarity [ @Craciun2005; @Craciun2010; @Craciun2006; @CraciunFeinberg2006 ] and the work led by Mincheva and Roussel on static oscillations [ @Mincheva2011a; @Mincheva2007a; @Mincheva2007 ] are two particularly influential approaches. For a good overview of the various graph theoretic developments, we mastermind the interested reader to the review put up in Domijan and Kirkilionis [ @Domijan2008 ]. 

 In this work, we focus on applications of graph theory, chiefly the Matrix - Tree Theorem (MTT), for deriving balance solutions (ES) for CRNs that fit within a Laplacian dynamics model. The MTT - based framework was first applied in a biological context by King and Altman [ @King1956 ] to derive steady department of state rate equations in enzyme kinetics. This framework was then simplified and summarized into rules (known as * Chou ’s graphical rules * [ @Lin2013 ]) by Chou and coworkers [ @Chou1989; @Chou1990; @Chou1993 ]. Chou [ @Chou1989 ] has besides extend the framework for non - regular state enzyme - catalyzed systems. 

 The main disadvantage of Chou ’s graphical dominion is that they are only applicable if the underlying digraph structure is * strongly connect *, i.e., every vertex is reachable from every other vertex. This issue was solved and extended for general direct graphs (* digraphs *) by Mirzaev and Gunawardena in 2013 [ @MirzaevGunawardena2013bmb ] and is applicable to the specific class of linear ordinary differential equations (ODEs) know as * Laplacian moral force *. Systems described by Laplacian dynamics are created using a weakly connected digraph, $ G$, with $ n$ vertices, with pronounce, directed edges, and without self loops. Note that by * weakly connected * we intend that the graph cannot be expressed as the union of two disjoint digraph. If there is an edge from vertex $ j$ to vertex $ i$, we label it with $ e_{ij}>0 $, and with $ e_{ij}=0 $ if there is no such boundary. [ ^3 ]
 degoted their efforts to dtveloping a systems-level underatanding of biochemical reaction netxorkw. In kcrticular, the study uf these bhemical eeacuion networks (CRNs) using bkeir zdsocnaved graph strucjure has attsacted considesaclz attention. The work led by Craciun wnd Feimbfrg on multistwtiomwritg [@Craciun2005; @Craciun2010; @Craciun2006; @CraciunFsinberg2006] and the work led by Mincheva and Roussfl oj stable oscillatilns [@Minchevq2011a; @Mybcheva2007a; @Mincfeva2007] are two particulatly influential approaches. For a good overview if thf various greph thvoretic developments, fe direvt the interesbed rxadee to the review provived in Domijan and Kyrkilioniv [@Bomijan2008].

In this work, ww docus on dpplkxatkona pf graph thxory, mainly the Matrix-Rree Theorem (MTT), fot qvtiving equiljbrium sjlutions (ES) for CRNs that fit within a Napmacian dynamics framewoek. The MTT-based frameaork was sirst applied in a biological context by King and Altmen [@Kiun1956] gi ferive steady state rate equations in enzyme hjnttibs. This framework was then somolosied and summxrized info rules (known as *Fhou’s gtaphicql rules* [@Jin2013]) ny Chou and coworkers [@Chou1989; @Xhou1990; @Chou1993]. Chjy [@Chou1989] has also excended the fxamewotk for non-steady state enzyme-eatalyaed systems.

Hhe main sksadvantage of Cfou’x craphical rules is that thqy are onoy a'plicabld if the ugderlying figraiv structure is *strlngly cmnnected*, i.f., every vertex is reachable from every other festeq. This isfue wss solved and extended for yeneral birectdd graphs (*sigraphv*) by Mirzaed and Gunawarggna in 2013 [@MirzaxvGunawarqena2013vmb] qnd is xoplicable to tne specifpc claws of linear ordinery ajfferential equcuiobs (ODEs) known ax *Lxplwcpan drtamics*. Systeks ddscfobed cy Laplaciak dhnamocs are created usinc a seakly connected dogvaph, $G$, wijh $n$ vertyces, with labrlled, directed edgts, and wivhout xels loops. Note that by *weakly cohnected* wf mcan that the drapm cagnot be ex'ressed as the union of two disjoint digcaphs. If there is an edte from vertex $j$ to yertex $i$, we kabel it with $a_{ij}>0$, and with $e_{ij}=0$ if rhere is no such tdge. [^3]
devoted their efforts to developing a systems-level biochemical networks. In the study of using associated graph structure attracted considerable attention. work led by Craciun and Feinberg multistationarity [@Craciun2005; @Craciun2010; @Craciun2006; @CraciunFeinberg2006] and the work led by Mincheva and Roussel stable oscillations [@Mincheva2011a; @Mincheva2007a; @Mincheva2007] are two particularly influential approaches. For a good of various theoretic we direct the interested reader to the review provided in Domijan and Kirkilionis [@Domijan2008]. In this we focus on applications of graph theory, mainly Matrix-Tree Theorem (MTT), for equilibrium solutions (ES) for CRNs fit a Laplacian framework. MTT-based was first applied a biological context by King and Altman [@King1956] to derive steady state rate equations in enzyme kinetics. framework was and summarized rules as graphical rules* [@Lin2013]) and coworkers [@Chou1989; @Chou1990; @Chou1993]. Chou extended the framework for non-steady state enzyme-catalyzed systems. main disadvantage Chou’s graphical rules is that they only applicable if the underlying digraph structure is connected*, i.e., every vertex is reachable from every other vertex. This issue was solved and general directed graphs (*digraphs*) Mirzaev and Gunawardena 2013 and applicable the specific of linear ordinary differential equations (ODEs) known as *Laplacian dynamics*. Systems by Laplacian dynamics are created using a weakly connected digraph, $n$ with labelled, directed and without self loops. that *weakly connected* we mean graph be union two digraphs. If there is edge from vertex $j$ to $i$, we label it if there is no such edge. [^3]
 devoted their efforts to deveLoping a sysTems-lEveL unDeRstaNdinG of biochemical REactIon networks. In particulaR, the sTuDY of tHEsE chemIcal reaCTiON NetWoRkS (CRns) USiNg theIr aSsociatEd graph strUctUrE has attracteD CoNsiderable AttEntion. The worK leD by CraCiUn aND FeinBerG on muLtistaTIonariTy [@Craciun2005; @crACiun2010; @CrACiun2006; @CraCIUnfeinBerg2006] and the work led BY MINcheva and RoussEl on stAbLE oSCIllAtiOns [@MinchevA2011a; @minchEVa2007a; @MincHEvA2007] ARE twO Particularly iNfluential aPProAches. FOr A goOD overvIew of ThE VarIous graph thEoreTic develoPments, WE direct THe interEsted rEadEr tO the REvIeW prOvIDed IN DOmiJAn aNd KirkilIoNiS [@DomiJan2008].

IN THIS worK, we FocuS on apPlications of gRapH theORy, mAinly The MaTrix-trEe TheOrem (MTt), for dErIving equilibriuM solUtions (ES) fOr CrNS thAt Fit wiTHin a LaPlaCiaN dynamiCs frameWOrk. thE mtt-bAsed framework was fiRsT APpLied in a bIologiCAl CoNText by KiNg And altmAN [@king1956] tO derIVe Steady stAte ratE EqUaTions in EnZyme kiNeTicS. ThIs fraMEworK was thEn simpliFied aND summarized intO Rules (known as *CHOu’S GRaPHicaL ruLes* [@Lin2013]) by ChoU and COworKers [@cHoU1989; @ChOU1990; @Chou1993]. chou [@CHoU1989] HaS Also extended the framEwOrk for Non-stEady state enzyMe-catalyzeD SYStems.

The Main DIsADvantage of Chou’S grapHical rules IS that theY are oNly appliCable if thE UNderlyinG diGraPh sTruCTUrE is *strongly coNNEcteD*, i.E., every vErtEx is reaChaBle FroM evErY other verTex. This iSsUe WaS sOlvEd and EXtended fOr GenErAl dIrectED graphS (*digrAphs*) By miRZaeV and GunAWaRDEna iN 2013 [@MIrZaevgunAwArdenA2013bmb] ANd iS applicAble to the SpeCIfic ClAsS of lineAr ordinary difFeRential equAtIonS (ODEs) kNOWn as *LaplAcian dynamics*. Systems desCRibed by lapLaciaN dynAmics are cReaTed usiNg a WEakly cOnnectEd digRaPh, $G$, WITh $n$ veRTIcEs, wItH labelled, dIREctEd edgEs, And wIthout sElf loops. Note that by *WEakLy connected* we MeaN thaT THe GraPH cANnoT bE ExpRESsed as the union oF two disjoiNt DIgRaphs. If theRE is An Edge froM vertex $J$ to veRTex $i$, we lAbel it witH $e_{ij}>0$, and wiTh $E_{ij}=0$ iF THerE is no such eDge. [^3]
 devoted their efforts todeveloping  a sy ste ms- le velunde rstanding of b i oche mical reaction network s. In  p a rtic u la r, th e study of t hes ech emi ca l  r eacti onnetwork s (CRNs) u sin gtheir associ a te d graph st ruc ture has att rac ted co ns ide r ableatt entio n. The work l ed by Cra ci u n andF einberg o nmult istationarity [@C r ac i un2005; @Craci un2010 ;@ Cr a c iun 200 6; @Craciu nF einbe r g2006]a nd t h e w o rk led by Min cheva and R o uss el onst abl e  oscil latio ns [@M incheva2011 a; @ Mincheva2 007a;@ Minchev a 2007] a re two  pa rti cula r ly  i nfl ue n tia l  a ppr o ach es. Forago od ov ervi e w o f th e v ario us gr aph theoretic  de velo p men ts, w e dir ectth e int ereste d rea de r to the review  pro vided inDom ij anan d Kir k ilioni s [ @Do mijan20 08].

I n  th is w o rk , we focus on appl ic a t io ns of gr aph th e or y, mainly t he  Ma trix - T ree T heor e m(MTT), f or der i vi ng  equili br ium so lu tio ns(ES)f or C RNs th at fit w ithin a Laplacian dy n amics framewo r k. T he MTT- bas ed framewor k wa s  fir st a p pl ied in abiolo gi c al context by King and  A ltman[@Kin g1956] to der ive steady s t ate rate  equ a ti o ns in enzyme k ineti cs. This f r ameworkwas t hen simp lified an d summariz edint o r ule s (k nown as *Chou ’ s  gra ph ical ru les * [@Lin 201 3])  by  Ch ou  and cowo rkers [@ Ch ou 19 89 ; @ Chou1 9 90; @Cho u1 993 ].  Ch ou [@ C hou198 9] ha s al so  e x ten ded the fr a m ewor kfo r no n-s te ady s tate enz yme-cat alyzed sy ste m s.
Th emain di sadvantage of  C hou’s grap hi cal  rules i s that t hey are only applicable if theund erlyi ng d igraph st ruc ture i s * s trongl y conn ected *,  i. e . , eve r y  v ert ex  is reacha b l e f rom e ve ry o ther ve rtex. This issue w a s s olved and ext end ed f o r  g ene r al dir ec t edg r aphs (*digraphs *) by Mirz ae v  a nd Gunawar d ena  i n 2013[@Mirza evGun a wardena 2013bmb]and is ap pl icab l e  to  the speci fic clas s of line a r ord i na ry di ffe rentia lequ ation s (ODE s ) k nownas *La pl aciandynam ic s*. Syst ems described by Laplac ian dy namic s a re create d u s ing  a weakly  con nected dig rap h,$G$,wit h  $n$vert i ce s,w ith l abel l ed, direc t ed  ed g e s,  and withou t s elf  loop s.N ote th at b y *weakly connect e d* we mean tha t th e gra phc anno tbe expressed a s t he u nion oftw o disjointdigraphs .I f the re isan edg e fromv e rt e x $j$to v ert ex $i$, w e l ab e l it wi th  $ e _{ij}> 0$,an d with  $e_{i j }=0$ i f there is no su ch ed g e . [^3 ] 
 devoted_their efforts_to developing a systems-level_understanding of_biochemical_reaction networks._In_particular, the study_of these chemical_reaction networks (CRNs) using_their associated graph_structure_has attracted considerable attention. The work led by Craciun and Feinberg on multistationarity [@Craciun2005;_@Craciun2010;_@Craciun2006; @CraciunFeinberg2006]_and_the_work led by Mincheva and_Roussel on stable oscillations [@Mincheva2011a;_@Mincheva2007a; @Mincheva2007]_are two particularly influential approaches. For a good_overview_of the various_graph theoretic developments, we direct the interested reader to_the review provided in Domijan and_Kirkilionis [@Domijan2008].

In this_work,_we_focus on applications of_graph theory, mainly the Matrix-Tree Theorem_(MTT), for deriving equilibrium solutions (ES)_for CRNs that fit within a Laplacian_dynamics framework. The MTT-based framework was_first applied in a biological_context by_King and Altman [@King1956] to_derive steady state_rate equations_in enzyme kinetics._This framework was then simplified and_summarized into rules_(known as *Chou’s graphical rules* [@Lin2013])_by_Chou and coworkers_[@Chou1989;_@Chou1990;_@Chou1993]. Chou_[@Chou1989] has also_extended_the framework_for_non-steady state enzyme-catalyzed systems.

The main disadvantage_of_Chou’s graphical rules is that they are_only applicable if the_underlying_digraph structure is *strongly_connected*, i.e., every vertex is_reachable from every other vertex. This_issue was_solved and_extended for general directed graphs (*digraphs*) by Mirzaev and Gunawardena in_2013 [@MirzaevGunawardena2013bmb] and is applicable to_the specific class of_linear ordinary_differential_equations (ODEs) known_as_*Laplacian dynamics*._Systems described by Laplacian dynamics are created_using a_weakly connected digraph, $G$, with $n$_vertices, with labelled, directed_edges,_and without self loops. Note that_by *weakly connected* we mean that_the graph cannot be expressed_as_the_union of two disjoint digraphs._If there is an edge from_vertex $j$ to_vertex $i$, we label it with $e_{ij}>0$,_and_with $e_{ij}=0$ if there is no_such_edge. [^3]