function dxdt=MichaelisGSGO3(t,x)

global CGS DGS CGO DGO kfGLU krGLU kfNH3 krNH3 kfGLN krGLN kfALPK krALPK kcatGS kcatGO Con1 Con2;

GS = x(1);
GLU =x(2);
GSGLU=x(3);
GLN=x(4);
GO =x(5);
GOGLN=x(6);
GLNGOALPK=x(7);
GLUGSNH3=x(8);

%CONSERVATION CONDITIONS
%GSGLU + GLU +  NH3 + 2*(GLUGSNH3 + GLN + GOGLN + GLNGOALPK) = Con2;
%GLUGSNH3 + (GLU + GLNGOALPK + ALPK + GOGLN + GLN + GSGLU) = Con1 ;
NH3= Con2-2*(GLUGSNH3 + GLN + GOGLN + GLNGOALPK)-GSGLU - GLU;
ALPK = Con1 - (GLUGSNH3 + GLU + GLNGOALPK + GOGLN + GLN + GSGLU);

%%TRANSCRIPTION TERMS
TranscriptionGS= CGS-DGS*GS;
TranscriptionGO= CGO-DGO*GO;

%%ENZYMATIC REACTIONS
VGSGLU=  kfGLU*GS*GLU  -krGLU*GSGLU;
VGSNH = kfNH3*GSGLU*NH3 -krNH3*GLUGSNH3;
VGOGLN = kfGLN*GO*GLN-krGLN*GOGLN;
VGOAL= kfALPK*GOGLN*ALPK -krALPK*GLNGOALPK;


%1
dGS= TranscriptionGS - VGSGLU + kcatGS*GLUGSNH3;
%2
dGLU= 2*kcatGO*GLNGOALPK   -VGSGLU; %added gogat term
%3
dGSGLU= VGSGLU - VGSNH;
%4
dGLN = kcatGS*GLUGSNH3 -VGOGLN ; %added GOGAT term
%5
dGO = TranscriptionGO -VGOGLN + kcatGO*GLNGOALPK; 
%6
dGOGLN= VGOGLN -VGOAL;
%7
dGLNGOALPK = VGOAL - kcatGO*GLNGOALPK;
%8
dGLUGSNH3= VGSNH - kcatGS*GLUGSNH3;

dxdt=[dGS; dGLU; dGSGLU; dGLN; dGO; dGOGLN; dGLNGOALPK; dGLUGSNH3];