embedded_kinetic_velocities

0001 function v = embedded_kinetic_velocities(c,network_combined,kinetics)
0002 
0003 eval(default('kinetics','network_combined.kinetics'));
0004 
0005 % evaluate rates for the network model
0006 
0007 v = network_velocities(c,network_combined,kinetics.kinetics_network);
0008 
0009 % evaluate rates for the kinetic models and insert them
0010 % if there are overlaps between kinetic models, the first one of them wins
0011 
0012 for it = length(network_combined.kinetics.kinetic_models):-1:1,
0013   kk = network_combined.kinetics.kinetic_models{it};
0014   my_c = c(network_combined.kinetics.mapping_metabolites{it});
0015   vv = network_velocities(my_c,kk);
0016   if isfield(kinetics,'enzyme_adjustment'),
0017     vv =  vv .* kinetics.enzyme_adjustment{it};
0018   end
0019   v(network_combined.kinetics.mapping_reactions{it}) = vv;
0020 end
0021