1. #WRITE THE MODEL STRING HERE AND LOAD IT TO TELLURIUM
2. model_stringC4='''
3.   # Identify variables and specify initial values:
4.   AcP = 0; HOAc = 30; ACoA = 0; Pta = 0; Acs = 0; LacI = 0;
5.  
6.   # Reactions with rate laws:
7.   R1: -> ACoA ; k1  
8.   R2: ACoA -> ; k2*ACoA
9.   R3: ACoA -> AcP; k3*Pta*ACoA/(Kpta + ACoA)
10.   R4: AcP -> HOAc; kf*AcP
11.   R5: HOAc -> AcP; kr*HOAc
12.   R6: HOAc -> ACoA; k4*Acs*HOAc/(Kacs+HOAc)
13.   R7: HOAc ->; k5*HOAc
14.   R8: -> Acs; k6*(AcP^n1)/(Kacp^n1 + AcP^n1)
15.   R9: -> LacI; k7*(AcP^n2)/(Kacp^n2 + AcP^n2)
16.   R10: -> Pta; Vmax*(1-(LacI^n3/(LacI^n3 + Krpta^n3)))
17.   R11: Pta -> ; k9*Pta
18.   R12: Acs -> ; k9*Acs
19.   R13: LacI ->; k9*LacI
20.   R14: -> HOAc; k10
21.   # Kinetic parameters:
22.  k1 = 0.5; k2 = 10; k3 = 80000; kf = 1; kr = 1; k4 = 800; k5 = 0.005; k6 = 0.002; k7 = 0.0001; Vmax = 0.002; k9 = 0.06; k10 = 0.005 ;n1 = 2; n2 = 2; n3 = 2;
23.  Kpta = 0.06; Kacs = 0.1; Kacp = 10; Krpta = 0.000001
24.   '''
25.  
26. modelC4 = te.loada(model_stringC4)
27.  
28. #SIMULATE THE MODEL HERE
29. modelC4.resetAll()
30. resultC4 = modelC4.simulate(start=0, end = 3000, points=10000)
31.  
32. #REPLACE NAMES OF THE VARIABLES USED HERE TO ONES YOU USED IN THE SIMULATION
33.  
34. plt.figure(1, (8, 8))
35.  
36. ax = plt.subplot(2, 2, 1)
37. plt.plot(resultC4['time'], resultC4['[AcP]'],label='AcP')
38. plt.plot(resultC4['time'], resultC4['[HOAc]'], label='HOAc')
39.  
40. plt.legend()
41. plt.xlabel('Time, min')
42. plt.ylabel('Concentration, mM')
43.  
44. ax = plt.subplot(2, 2, 2)
45. plt.plot(resultC4['time'], resultC4['[ACoA]'],label='ACoA',color="red")
46. plt.legend()
47. plt.xlabel('Time, min')
48. plt.ylabel('Concentration, mM')
49.  
50. ax = plt.subplot(2, 2, 3)
51. plt.plot(resultC4['time'], resultC4['[Pta]'],label='Pta',color="purple")
52.  
53. plt.legend()
54. plt.xlabel('Time, min')
55. plt.ylabel('Concentration, mM')
56.  
57. ax = plt.subplot(2, 2, 4)
58. plt.plot(resultC4['time'], resultC4['[Acs]'],label='Acs',color="violet")
59. plt.plot(resultC4['time'], resultC4['[LacI]'], label='LacI',color="darkred")
60.  
61. plt.legend()
62. plt.xlabel('Time, min')
63. plt.ylabel('Concentration, mM')
64.  
65. plt.tight_layout()