forked from BIMIB-DISCo/scFBA
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathsetScFBABound.asv
124 lines (109 loc) · 6.22 KB
/
setScFBABound.asv
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
function [model] = setScFBABound(model,nPops)
% Compute some statistic over the single cell expression profile and
% generate the field geneNameOff. May export the resoult in excel and draw
% an histogram of the distribution for each significative genes.
%
% USAGE:
%
% [model] = setScFBABound(model,nPops)
%
%
% INPUT:
% model: metabolic model in COBRA format
%
% OPTIONAL INPUTS:
% CutOff: A gene is defined significative if the mean of its expression
% profile have a ratio less then CutOff between the bulk
% expression level. (Default = 0.1)
% Hist: If TRUE draw a histogram of each significative genes
% expression profile. (Default = FALSE)
% ExpXLS: If TRUE export the result in a excel file. (Default = FALSE)
% Epsilon: Value added to single cell transcript levels. (Default = minimum of all expression level)
% nBins: If Hist = TRUE states the number of bins to draw.
% FitFunction: If Hist = TRUE states the function with which fits the distribution of
% each significative genes and draw this curve on the histogram.
% (Default = 'kernel')
%
% OUTPUT:
% scStruct: Single Cells dataset with fields added.
%
%
% .. Author:
% - Chiara Damiani 30/01/2018
%LOOPS
model = EditBoundaries(model, 'SLC25A5_SLC25A4_SLC25A6', 0,10000, false, false);
%*********************************
%%% vincoli exchange e cooperation reactions
% set boundaries - EX reaction
model = EditBoundaries(model, 'Ex_biomass', 0,10000, false, false);
model = EditBoundaries(model, 'Ex_glucose', -100*nPops,0, false, false);
model = EditBoundaries(model, 'Ex_O2', -100*nPops,0, false, false);
model = EditBoundaries(model, 'Ex_glutamine', -100*nPops,0, false, false);
model = EditBoundaries(model, 'Ex_folate', -100*nPops,0, false, false);
model = EditBoundaries(model, 'Ex_lactateL', 0,0, false, false);
model = EditBoundaries(model, 'Ex_urea', 0,0, false, false);
model = EditBoundaries(model, 'Ex_glutamate', 0,0, false, false);
model = EditBoundaries(model, 'Ex_NH3', 0,0, false, false);
model = EditBoundaries(model, 'Ex_arginines', -100*nPops,0 , false, false);
model = EditBoundaries(model, 'Ex_citratec', 0,0, false, false);
model = EditBoundaries(model, 'Ex_Carnitine', -100*nPops,0, false, false);
model = EditBoundaries(model, 'Ex_Palmitate', 0,0, false, false);
model = EditBoundaries(model, 'Ex_AKG', 0,0, false, false);
model = EditBoundaries(model, 'Ex_proline', 0,0, false, false);
model = EditBoundaries(model, 'Ex_alanine', 0,0, false, false);
model = EditBoundaries(model, 'Ex_aspartate', 0,0, false, false);
model = EditBoundaries(model, 'Ex_asparagine', 0,0, false, false);
model = EditBoundaries(model, 'Ex_serine', 0,0, false, false);
model = EditBoundaries(model, 'Ex_glycine', 0,0, false, false);
model = EditBoundaries(model, 'Ex_cystine', 0,0, false, false);
%set boundaries EX_SC reactions
model = EditBoundaries(model, 'Ex_SC_mercaptopyruvatec', -10000,0, false, false);
model = EditBoundaries(model, 'Ex_SC_H2O', -10000,10000, false, false);
model = EditBoundaries(model, 'Ex_SC_Pi', -10000,10000, false, false);
model = EditBoundaries(model, 'Ex_SC_H', -10000,10000, false, false);
model = EditBoundaries(model, 'DM_glucose', 0,0, false, false);
model = EditBoundaries(model, 'DM_glutamine', 0,0, false, false);
model = EditBoundaries(model, 'DM_folate', 0,0, false, false);
model = EditBoundaries(model, 'DM_lactateL', 0,10000, false, false);
model = EditBoundaries(model, 'DM_urea', 0,10000, false, false);
model = EditBoundaries(model, 'DM_glutamate', 0,10000, false, false);
model = EditBoundaries(model, 'DM_NH3', 0,10000, false, false);
model = EditBoundaries(model, 'DM_arginines', 0,0, false, false);
model = EditBoundaries(model, 'DM_Carnitine', 0,0, false, false);
model = EditBoundaries(model, 'DM_Palmitate', 0,10000, false, false);
model = EditBoundaries(model, 'DM_CO2', 0,10000, false, false);
model = EditBoundaries(model, 'DM_AKG', 0,0, false, false);
model = EditBoundaries(model, 'DM_proline', 0,0, false, false);
model = EditBoundaries(model, 'DM_alanine', 0,0, false, false);
model = EditBoundaries(model, 'DM_aspartate', 0,0, false, false);
model = EditBoundaries(model, 'DM_asparagine', 0,0, false, false);
model = EditBoundaries(model, 'DM_serine', 0,0, false, false);
model = EditBoundaries(model, 'DM_glycine', 0,0, false, false);
% set boundaries - coop reaction (siccome ultimo argomento ? false la
% cambia automatiamente in tutte le cellule)
model = EditBoundaries(model, 'Glucose_DM_COOP', -10000, 0, false, false);
model = EditBoundaries(model, 'Glutamine_DM_COOP', -10000, 0, false, false);
model = EditBoundaries(model, 'Glutamate_DM_COOP', -10000,10000, false, false); %la reazione con consumo di ATP per il demand non viene utilizzata
model = EditBoundaries(model, 'NH3_DM_COOP', -10000,10000, false, false);
model = EditBoundaries(model, 'Palmitate_UP_COOP', 0, 10000, false, false);
model = EditBoundaries(model, 'Palmitate_DM_COOP', 0, 10000, false, false);
model = EditBoundaries(model, 'LactateL_DM_COOP', -10000, 10000, false, false);
model = EditBoundaries(model, 'Oxygen_DM_COOP', -10000,0, false, false);
model = EditBoundaries(model, 'Arginine_DM_COOP', -10000, 0, false, false);
model = EditBoundaries(model, 'Folate_DM_COOP', -10000, 0, false, false);
model= EditBoundaries(model, 'Urea_DM_COOP', 0,0, false, false); %urea
model = EditBoundaries(model, 'Carnitine_DM_COOP', -10000,0, false, false);
model = EditBoundaries(model, 'AKG_DM_COOP', 0,0, false, false);
model = EditBoundaries(model, 'proline_DM_COOP', 0,0, false, false);
model = EditBoundaries(model, 'alanine_DM_COOP', 0,0, false, false);
model = EditBoundaries(model, 'aspartate_DM_COOP', 0,0, false, false);
model = EditBoundaries(model, 'asparagine_DM_COOP', 0,0, false, false);
model = EditBoundaries(model, 'serine_DM_COOP', 0,0, false, false);
model = EditBoundaries(model, 'glycine_DM_COOP', 0,0, false, false);
model = EditBoundaries(model, 'TranspCystineGlu', 0,0, false, false);
model = EditBoundaries(model, 'TranspCystineSer', 0,0, false, false);
%*********************************
%%%reazioni da silenziare
model = EditBoundaries(model, 'Complex1', 0,0, false, false);
model = EditBoundaries(model, 'Complex1ROS', 0,10000, false, false);
end