From b58eebbb4a9cf82e54275b2642dbc1b597db242f Mon Sep 17 00:00:00 2001
From: Viktoria Petrova <vipet103@hhu.de>
Date: Sun, 10 Nov 2024 12:34:14 +0100
Subject: [PATCH] add FBA assay and protocol
---
.../README.md | 0
.../dataset/.gitkeep | 0
.../isa.assay.xlsx | Bin 0 -> 6912 bytes
.../protocols/.gitkeep | 0
...luxBalanceAnalysisToDefinePathwaysProtocol.md | 9 +++++++++
5 files changed, 9 insertions(+)
create mode 100644 assays/FluxBalanceAnalysisToDefinePathways/README.md
create mode 100644 assays/FluxBalanceAnalysisToDefinePathways/dataset/.gitkeep
create mode 100644 assays/FluxBalanceAnalysisToDefinePathways/isa.assay.xlsx
create mode 100644 assays/FluxBalanceAnalysisToDefinePathways/protocols/.gitkeep
create mode 100644 assays/FluxBalanceAnalysisToDefinePathways/protocols/FluxBalanceAnalysisToDefinePathwaysProtocol.md
diff --git a/assays/FluxBalanceAnalysisToDefinePathways/README.md b/assays/FluxBalanceAnalysisToDefinePathways/README.md
new file mode 100644
index 0000000..e69de29
diff --git a/assays/FluxBalanceAnalysisToDefinePathways/dataset/.gitkeep b/assays/FluxBalanceAnalysisToDefinePathways/dataset/.gitkeep
new file mode 100644
index 0000000..e69de29
diff --git a/assays/FluxBalanceAnalysisToDefinePathways/isa.assay.xlsx b/assays/FluxBalanceAnalysisToDefinePathways/isa.assay.xlsx
new file mode 100644
index 0000000000000000000000000000000000000000..057ce443ebc9eb1bb3b4a6b28679b48e911b216b
GIT binary patch
literal 6912
zcmai3by(DE)235;=}xJo5tI~Am!-QUmlT!`0Y#({mJ&f43F!u5DJkg=0Tt<%ZqRR$
z7d^*w-tT?>*lVxdxo4j1p4plC-K(mAiiU%PgoK5ZY+$PwJ>870hm3@@jfR9og80-x
z7V6*%ad0)(^m2r_81Z=8KPrk-a%kZr3_Oy!w*Um+Fw)ez*OFP;LzcIXG~y}4O*6fl
z&@}Bq%%0E5kZ(f5l_;|1o!UmJQXL^bP5}~23Rk-Yz@8fBuqUrQPHxJ9F=b$BKXW2;
zDU&6RopTVM@yZHF0xo7o4oS(&Ln820m&EmpHF<;xPi}E%1v_Ko#vQfJKIcw}ur{YM
zAW(aetNXsYZ2{5)JGzTPMclRb@MF0%@)v3&nT0w-gM;80IxI>LcdGoVTx8q;ad~U5
zW#r&Sk5t7$oLI_$V+_24h4Ib^+#cTL_(@1C35Su^D)CD$9lRQg1yuh`D*jt+>*IyU
z_Lvvv%8KO%_eJoX_qV!;PJHLTKlTmi&%p#aSw%TW9QIej3}vMch*0TYBw-$J`p(80
z8OtwzxADW%=(mr{uJ|gxLxgwy!{~eR`)aX+M*)m$J#Qq4Cr{VP%*40?*U${yr3&Y$
zRsWL-q@S5UkpJsoAl?}}Lms*CUfWz{Yd1}PY8COF0tylm?k^jjOUt5|_sXq&H>9_G
z?k$|Xb^u{zWWyP<KC=l)HnTWTS4F&Od;V=JG)nqLgbik-)5fw@&w8JD#9D;<A<Yf$
z0xAG{Pk(TKbI*4hCjw4fj9_L=0mrNxNlTU&ip=1=t_GMMDX0_UarAMzTeXBz0#_uU
zTcyWFhGSrKHN?F+azDR#%k{L&Zz&+Jx?g&#j`i&Y|FG~Q2LtG8^rmShgO&byH@k%*
zZRj(v(XHelouv#;)xl*W?Fsm~L~bf*FeAM>*_ucrPkF9KzsZee+6a4j;Dv%r&5v4|
z*Ma`b9(vg3=K{kwMFxG4t{M3!Yp4iDJRkj{?4JhMR|fyaYyeA{_9e5m2xh7NYUTlT
zwljl5?Jj=_*F@X4@)62yUEq`*>AUNLIig)So*5b(p$1e|D?XC97O}rL(xL{@U}a#Y
zwH1C$`~X_-I}#01Q5?lAN#_I;t|*l@Yyo_Prq8Sh+olX8&|o-hEFVWb)<tN9mT8@R
zS)!H1IUIX;?tumHA1K`mS=m#f9)5efXSC~W-$M^}9momG$S2e_T$Qp4oZ)_LhSwZ(
zPwQo;4>63N-+r-j0uEqF4eI5kO1Wc!{=mqBK5AlKN1?;zOS{WynmB%e!3~0N?XRvf
zXZ;*f+x4L*$7H?pF?`01n_Y~%lnX+q_xe8MQF`snuhIW^qVcZ?x}?U%8Uk_svrKEk
z)an8dM7AIZr}&5MPYC}dpD*Lpph<j$eltg0Q8Uz{A}ow3`R-~}@<~GWDTf;rB!wvV
zt7EIr4&b5+3&rVU1ks|@lK5wh$$Lr5#S!>2)On}=&D@O9E4BF+O<bJ=!Dj`ep)Vl5
z?fPTk{`FY04&4MNvBR(&qL%wuD@u2%?1J5BsK10nCMHtJ!33*@X)7q!j2>YCbr<Q@
zrV`k?n7$=!bcL3VZ<)TWi#L1(nDNN0xfwzRtVn3ISup9ZjBLW3;J<iBU;J(jL+%lT
zYcu?L%-i|(0r*3kYWtmvG}Hd$#7%#ZV`b-}7x|kvoSu^Sp7rS^w}UwMz;om@IugT~
zJe|P8E%vB)y`TGzc)Fl>^5yyti!O>&C9!MmxWB1^nW}1vHrmOZ2!YC^EM5|aYE=%s
zkh<cBw^|$xOZ>A|7fGot^%Ct40JGhw-*%fL&#>$$ym%6a0I^FSyrDZ7Gy)X#-7GKy
z2%XD0<66rbM)9>BlVlIt7Wkp3YqAKzY@fuJ15sYQ*=#zfWWPbcHW?i+?J`i-Qvp^(
zv0NnM$px#SS%%QTzQbO5)sF=(FMT_&eROUVg+}T<r?V_qE1MwID9xv~LzqQa?`)m9
z!q~~QE$km!7t9ddEP~AjruR2QF=k76`vhwk+b7LL1_?_5$0ysQ1sN*_$EUzV_*t)i
zRHHbbDFI%2@af@To+Xh5!NF$6k0}!G2742+otLXA8%m>LQ_mck_w!(Nb_Ez7fm)SX
z<N4ocgg;gXO5T5aw<t~mb0_pn_TvR12ax}+$0>VfdgYj+z7NHRr(DU_BG20_Wh=a@
zb%$@%dXARIz!WKH03W1o8kXr#VIC-GS39Ll2EICc+RM;L)6^|>UFBB|b|<be(E%|I
z?<2+`;Hq+6tWBLE78<V3HV#%UKdUyb57B7}1GdBi7scFGks0H%a0Wr3N;w-*N5dF<
z^De3Ob?LoBQBoP{H*X~#A0>Hi(B6E#Rb@uM0@5ehE9i<76na^dJ=t-d-OB?NWhThM
zfZ%o1%g*gBFfr|NgYeqxcb<4~p=sgSriXscz)p#9qw>+MD@sq1eVcc8Uq5ux^)XNg
z)T#$Dm^YV9avG;P`7(E*!War-i&hyvxZ@iUpu#jn-WV=7SjLJ{g)dAs(CrYMt5ztN
ztEaC8y2zN!AA_#+4;9<5b6T`BQER`{74G1KZ{B}&@x0$jm|db(v$*+ao)SYaFRDnI
z{v@bTiAiJgJ3b%*M*}sVSw7phIYKDmIaJd7KCvSHu|<jX(@GU7^?2T}japrvu!B;V
zcz}Rm8%aPG-sFiaWH`ki2<fAX*e3==_=iWS2!j`tJPGDgzaB3%o!J<{`8J&_E26>v
z1kn_WI0cU93Yc@rakT3}Ol&PMh7*+Q{QTE$?rIIOhx}^#S80st7ryxdF+|%KuC^M)
zRJbzx`TDc>Z0e~%r*8p%aQc5MwlLUDOLs8saUZ#5<KMA4$jm%KW)2ntafYe^D5geu
zP%&Ivc~IS_O)0cZ(6u+p2}si0$9+R2<i9gdX(m<M44FD1C1pL?+WBxqM{nF#FrgxA
zOJmN16|Qf25@t^&x;q&qHy`f>9pn{Olm!9sh#*7`iGtj{>l4M9K+3O#)C{ztz)=BI
zj$DIyRv(J{b2v$JSg8>d1knLJ^aiq#%F=eAqA>f3Lc-`yAtr$TaKli@@ka_4&c5tj
zy1)`pm_>8$d9#XDuR_+x*{tOAki>0vkw_x@M81KjlpEW5LU>^T0hquMWvLvu&HO4r
zX1VmOn-eNhF15AIj_UG@$H$HwE2703cDx5V0vAc|9yaK#kHw#=BklU?+UUKhc({{%
z`xa*RYpTzmAJ^z>Z{NNc!}FO|`VR1M3-q|PebHD`Zi*i#{28x)+Qz*u|A|LniT&4g
zcli1!sSgegdng;V^>mR|L%Z#(D55vh`ko_hrw83MjWT+*E8PpsdH_4dJ7)C+extRZ
zn9isQ4mw79N4m^~C=&DF7pPeE9AOE5rf{>PO-`SiPnD;@5@FSR-SSqljUOq%#}qBg
z7Muo4z(Lz3243TOQSbBfjh!tkmDFU3?{1PuIiGHYb$HcINcrv!5ExTE{c$+Io;$O)
z!vuSJ(4$lG{diAlEn8+CZdH;h<$ZE7@$}+c6BX*>nh;C?3JG+=Or!A7!@?a^^dmTd
z%ilne17HJ?)KNFOq>;ZQGUm6smD;wP(Fg?MdooL^zoSDFkC#KOrwzbye8oRyY6wrB
z?lL@CY#Y**x((ZV2gjmaUpgG8`4pZ@xGR@Cj9V(>lICACMDakEs!Q<JdtK&3Kf?B6
zc-yjKX|xa>;hX0n<1DDo><2Q_END%=u4wPa!tC9*sA=!Jt-+6V@3mA*YvvMbw7(Eh
z<yLzOc_H)q!8ZoV6%bwTyAldR(DpushUI>cE>=n}T^ffgR5oY`?)$1=hO~<{aV1JD
zq)fcLhwY(XOJ2NSsK5AOPJ(m@lX;pv5Ig5vR2IP#X|)Xa776!^y}``$8(9QV3-Y~g
z1}K7VQvEAfrYXC(2wrYXEy+>Srg#<x9O8EXQqmbaIvSGhh$jh1tQsI&ytFT@=mTbT
z*bc)={RBAj?1Um&Ww0kv#Vyf&+^y_Ft%5QMcbNJ-*lauZzs8N^3`?x3#A$*(%WV!=
z4&_4nv)!bm!Md!z8L$aT*cUh0v)&=^&aWQv&)S|iOy3}EKlP(zVCQpPC6hC20F|ex
zm+5kesEcfN<BjV#ki@JMwGZe|&OdowzWjz%oZNT~kR5k`l8bdj&uv%79(l;0%<@XB
z52CgAw9q!}a2En}(d#iB7f36N*L%sU{T-;1&&utyD998?*Gbtf#su_1+qf{L#~4xO
ztn_IW!b`uMru;tQHh3DtY8RL3`Ew*C)7Q>-T6KH1V1fJ)6!6UoRS*5g?Xx`8mTwk!
z=d#(b7;344hwb?B17G*vw>ub{hG(BvijaAsIOOl2vpC1*vE)+`AFQ2#;x@z0hDJ-R
z<7V1gy1mfxZ=$k7Ta9j-qs?M4N+_Pkb@A4F&!+MS<kW41S$Zigw}=--bp|GZnaG2X
zY3SBl$v4(5HMZ8m;_o-?@cY=?FJVySNz9`i8jhh0OTMn|B|0BL>vP08m5a$Eq<sJ!
zdF3M8q+*%q%m@deQF8mn6-luWyF8B}G?UoEpelOzlFRLEk9DF_N~%Y1UAf~8(IMo0
zXFl$L0P9pPJ(ltvZj52f_uFlvS^Q1m&8iL!$rO&D07HEjkSEVO9c@DeoxP&BhSi1D
z?3p7-eJXK20s~)dyyb7!M8>V2=5wcD6MVP`7j>KAsf-uG#gBSyN_OTM>=UV_HPAh_
zRA<@xWRGKK7E>ltUsk0kKUdq;-O%K2R8q0FNewtdL~2_|7pBjpqm$DRVPNm>i55K*
zXe*qfZo&IOIph40>?m4~4fKF~a#53t|0ZxRW3!sJWbQF5eD)BmJXWTFS7QXMGOE%4
zRtteo=ZT}d>15!BPk!<}wo2<?dUd#Z+|lm)j=ienn7dL}AJZVu@v)O}F~Il}H5&I`
zsrIT!L7#HbXkqmkYwom7-0O!$8f)YAk85TQKHabQ(J2(EQ94rmSh;9|eXKiG$Jy3h
zrfG<KECrmsIXWGYI`TxQM!3*t2hIoWUY}~Zo^e;h3ti0bp*v!GMnw!Ty#MSVu3nEI
zKbMc}I3vW$3<%mf;^LmsK*i$@^|Q4J)67fa6pnRY_meL?MWq~xRyx}f!yqPho60%t
z{!ju-+dCDkUCol_?G-XC0AT7!-d7~8Td;jKSCaH8$$V5dAIw>4#vxgXF=Ig=(?N&k
zJtb&KM5q}%*HCp36P5E(P(F#vND!U*<*gX~mDjX3wyoBa@8s2^29L~{ICMN2xFI`I
zgz|kZ2Ht`3!c{&!=rzHqfzjOYM`GD{pK%<g&;|9=c|w#t3kMZS!axrxKCg5?M2=g@
zpF{WFbh}FyKt_UH3S~wo+G#bN)l;{DP;0Je>wWZ17YpvO%5Ez7gm>GivpYrx*FGvx
zy|)X^m9thfot#sh=1zu43&vgNk6{nb`LXUwJs2`?0h0RIy(8&T&45yUb!Jp9mJ!F&
zL!ov^RJ-juNCCwcPEd0c#!yT<f^DQ^1T(zoXWPdc{cNkzqLq9&&khMSc^QZuv~^fq
zRmbKqsQ+<=o8$Ipn|pz>bXG5rO>o)TV&Nq-`MB;njAPw~Q)dG0jys>%V2l)0t;m9|
zG(VDWmCi)mGji*a$~coD-;WLZB%gB}?oDkblWyM=xVT0T^i0n^cH+rs-3MnqbPCo2
zYp#tYp@cb(@4etJ-P>^ooyL{46*<T6Q#1+TRD~4;5t|01RH<`m9tE)Q10%Y}E>DGz
zEGM{k$%bq<?^R5XJx<|0A8fjCJG7nIvHTE>5d~pp=y?%N+b5-@A2x93`NfkT&4;x|
z_bATNGMkTf?WDFX`BS|2j_>U5vSUr`JaQvn97;S|iD(K?@i!5m2LfLyI=1s{Li@Eu
zSysM7h!#AFElwLNX08kL-<L28sJW^$)bY=!z^=tsNC0A&W<xxe6C&#Pmzk-f<K-Tm
z5_1<i&37Yk>x?w$TST%0dyup?e-`y)-7+Nk>G|d^S8mOCv}oa&i_<RrxjXC1aETa#
zfi#Z%Y1!p4I+Dq+TxEuFnbJ@h`-?_m7#ymTcJJ<uMhtRs?YR&ZhXc$ic}6p{aLolY
zm6z6u^I)(5&|NxQ{LpJSU=(>-U=w0ZWw*_fFo7@8vtC{w?P;<vJg(z}1ITUD8u&r2
zk14o4H1dtYB~{zl%yk@TIc=0+JWD!bc5eNBx`wl)=56cnC1E|uHT5p9ylM=ZVwX~f
zW+Kj$an&QN*8W>2BsODv`43tdNiCi(7a4`mkrU#JzFc9{+vZlT*5J1#4~>0-zZ$47
z$Htnw#J>EcwmNO!S;)f4CXNI&{whSQ2Y>ZJ#j~~_3=jK5A7JpCSE5Nh!;0J@9H2kr
zdO}-=CAp=lyoM~k;k2#120!54r|_7?rx;LZe3WmGi@ko%Gs3Kh=ayy@OMSjhc22im
zv3y}78k##YxxUz!d9E8j&zkHFZZvxL9rZe~SDnWLjjF~Nk+FG1_aXXg#>}D4ke?}2
z@40;Rsox`w-pZSge?Xd@hxYa+C*u%KTi?g&aiKDxx5DeD4_-pq{v<UiDU&;r$$c;g
zk?*s;?iEVYHvw|#fze>7l|_Son+CDhg8pucbwV_ic;d27Vrd>9erihDSJh`Q|LDWK
zQo(ZBN{;M@8|e#A0DwDD)brp-<T_xL;l>Hsim1GU4R$B%17+&gkNA;UHF0L8m6$q}
z80xC*>1YwD<wS3+fc{(;Y3RT`Sv}4|;5mmx^z+k|F^0l?TTLZ>o0YVon#B<XG224d
z3QlB}fh<rUYGlWvUBg7qFy%$Gou{U*^p~77KDqDPJJ=}jeIyy~_qe&-L1JzT_n4B+
z3EUvb`_CGr`?5t_#^2Sw>1X%z_cjc2fdV?uPCu<fIoCa=`;Y{1#N^?4sR<KZt@E~b
zcbx<pBPql6%FVcGRYsEFM>)&krxm7?QoGO(v3I3r+M9ZtkpG?yOQ7r%Lqs-Q5!oQP
zngf3oqatFIR1kt+(AIev`%w|_Ir7+)JXVzZT^c{GU3HJ}c;5C^-5(MT526Xld7L)d
zV8&v4MJnSV%=d@CXq6`{tEX!;4ms#_)lB4076Ia#g-z;P!G|DcBjZ9zU-3r`Qds>S
znq9u&j9vkC?@DPG`oSIf^*4cS4xb<ZCK*}pHVY;9#O0n|fA?KQ_W7q2jrLV7m?d|@
zc~TgO@L04}ylJg7^$zHh13gz-6uqb^%Es>JJ_!JuYSm8OrZ=^Fx31wN5$$z`X<rqV
z>XWqe;*{(pcZO*oZ0?hDgKgFOqrFk|fbARZZ_Tp!q*u#^u#>>fO7l`XIt=c!jUz!R
z@OzC>dKY`dylO%pWvqs|3SO_A=kyfRvuAyc!JbzbeV{549Q;X**4Kc)ILT6!#6dF)
zThb5fbcci`l1A8u4jY$4c;-dYk}VR1ME;ptp#>lFYECCFw1Yaf{}fkM0T~5=^xsu-
zL^l3%{zWtY?ckp+^HSWsI+tBB43Xh$LGtfG{}g+#QAkL|zc4=)$^S3dr6ziHE`v2A
zO6jLUdL8MXx4YCxug)c4_6qRt7yMH*y#^xwTR?<G#JPg~wP^q6eJ`cOt8*Ei10lm8
z{QOO5yuSKX%lcDlyhbA-Y5fEHTX_7<*p-(*WuM<(2rz%%;2$E;Zx8>A!V;{X9<F7I
zzhD1Pq2hPo@xQ`;%ND=m{`<1{+d(X%XCSz9`L#O!hW@WB3laZiTK|S!Mf`KD|Bh(z
zSJ;1s{GaGQN6c?D9zt*XwIW;(oIl}z_Qc=tB+_4{axK@ZDxm-Ttq~&Dj}AmdvXEbz
F{~u^-Y-Rue
literal 0
HcmV?d00001
diff --git a/assays/FluxBalanceAnalysisToDefinePathways/protocols/.gitkeep b/assays/FluxBalanceAnalysisToDefinePathways/protocols/.gitkeep
new file mode 100644
index 0000000..e69de29
diff --git a/assays/FluxBalanceAnalysisToDefinePathways/protocols/FluxBalanceAnalysisToDefinePathwaysProtocol.md b/assays/FluxBalanceAnalysisToDefinePathways/protocols/FluxBalanceAnalysisToDefinePathwaysProtocol.md
new file mode 100644
index 0000000..d9d3d1c
--- /dev/null
+++ b/assays/FluxBalanceAnalysisToDefinePathways/protocols/FluxBalanceAnalysisToDefinePathwaysProtocol.md
@@ -0,0 +1,9 @@
+## Flux balance analysis (FBA) to define pathways
+
+As we suspect that the root may convert sucrose into other sugars before secretion, we would like to find out the genes involved in these conversion pathways for further analysis. To assign genes in the metabolic network to these pathways, we applied FBA (Orth, J.D. et al., 2010), which searches for the combination of genes and reactions that results in the highest yield for such conversions, assuming sucrose to be the sole input carbon source of the metabolic network.
+
+We used the AraGEM (de Oliveira Dal’Molin, C.G. et al., 2009) a curated FBA model of *Arabidopsis*, to perform our simulations and define the genes associated with each conversion pathway. The AraGEM model is published with 3 different sets of parameters that correspond to their own case of metabolism: (a) photosynthesis, (b) photorespiration, and (c) non-photosynthetic cell. Case (c) is relevant to the root condition, as it considers the conversion of sucrose into other biomass metabolites. Hence, our simulations are developed from case (c), which uses sucrose as the input carbon source.
+
+We switched off the default objective function of AraGEM, which includes a set of metabolites with weights found in real plants. To investigate the conversion of sucrose into different sugars, we set different monosaccharides, disaccharides and combinations of monosaccharides as the objective; these include all possible monosaccharides and disaccharides in AraGEM. In total, we analyzed 11 different objective functions: (1) maltose, (2) glucose, (3) beta-D-fructose, (4) D-galactose, (5) D-xylose, (6) glucose + beta-D-fructose, (7) glucose + D-galactose, (8) glucose + D-xylose, (9) beta-D-fructose + D-galactose, (10) beta-D-fructose + D-xylose, (11) D-galactose + D-xylose. Note that the molecules within objective functions have equal weights, as some of them have multiple sugar molecules.
+
+We used the python package COBRApy (Ebrahim, A. et al., 2013) along with the GUROBI solver (Gurobi Optimizer Version 3.0. Houston, Texas: Gurobi Optimization, Inc., April 2010. (software program)) to perform FBA simulation. For each of the 11 objective functions, we performed parsimony-FBA simulation (Holzhütter, H.-G., 2004) to find out the most efficient reactions that perform the conversion of sucrose into the objective sugars. We also calculated the shadow price (Reznik, E., 2013) of each metabolite involved in these reactions; a metabolite with negative shadow price is deemed critical, as an extra supply of this metabolite leads to a higher flux of the objective function. Reactions that have critical metabolites both at their input and output are also considered critical. We collected all the critical reactions of an objective function, and the genes associated with any critical reactions are assigned to the pathway (refer to Table S6 for the genes in these pathways).
\ No newline at end of file
--
GitLab