How to introduce a delay in a simbiology model?

8 次查看（过去 30 天）
Archishman Ghosh
Archishman Ghosh 2022-4-13
评论： Jeremy Huard 2022-4-14
I have a model here for transcription and translation which I am trying to fit to data containing a column for mRNA and for protein.
TsR and TlR are transcription and translation resources.
The data has a delay and is pretty sigmoidal. The michaelis-menten based ODEs I used don't fit. How do a introduce that delay time factor so that it fits by recognizing the t-t0 time delay? (mS or GDmS is the protein)
Data attached.

请先登录，再回答此问题。

回答（1 个）

Arthur Goldsipe
Arthur Goldsipe 2022-4-13
SimBiology does not directly support delays. One approach would be a more mechanistic model of the delay. For exampe, you might try switching to Hill kinetics to see if that can reproduce the sigmoidal response. Alternatively, you could implement the model with basic MATLAB code and use dde23 to solve the equations.
  2 个评论
Archishman Ghosh
Archishman Ghosh 2022-4-13
I was afraid of this. I really wanted to use the simbiology module. I tried Hill, didn't work. I suppose I will have to use the delay differential equations solver.
Jeremy Huard
Jeremy Huard 2022-4-14
Ran in:
@Archishman Ghosh
One way to do this in SimBiology would be to approximate the delay with a chain of reactions of type with a reaction rate where
sbioloadproject('timeDelay_reactionchain.sbproj')
cs = getconfigset(m1);
cs.SolverOptions.OutputTimes = 0:0.05:25;
simdata = sbiosimulate(m1);
[time,data] = selectbyname(simdata,{'y','ydelay'});
plot(time, data(:,1), time, data(:,2));
xlabel('Time [hrs]');
legend('y','y_{delay}');
title('Delay = 2.5 hrs');
grid on;
set(gca, XLimitMethod="padded", YLimitMethod="padded");
Alternatively, one could write a function that keep tracks of the previous values of y and interpolate between them to find .
I hope this helps.
Jérémy

请先登录，再进行评论。

社区

更多回答在  SimBiology Community

类别

在 Help Center 和 File Exchange 中查找有关 Extend Modeling Environment 的更多信息

标签

产品


版本

R2021a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by