# Time-dependent Growth Rate of Species in Model

4 visualizaciones (últimos 30 días)
mitpi_182 el 15 de Jun. de 2022
Editada: Florian Augustin el 29 de Jul. de 2022
Hello,
I'm trying to implement in Simbiology the model of a system which describes the kinetics of many species. The kinetics of one these species is non-standard in that it's change rate is time dependent with different kinetics in three different periods. In the first period, which occurs in the first five days, the species is able to distribute from its initial compartment to another at a fixed rate. Following the distribution of the species into the second compartment, in the second period which occurs for the next seven days, the species grows at a fixed constant rate. Finally, for the next seven days during the third period, the species decays at a fixed constant rate. Thus, if I'm not mistaken, the species has a time-dependent growth rate that can be described as a step function in time with changes in time = 5, 12, and 19 days. Is there a way implement these kinetics using Simbiology?
##### 0 comentariosMostrar -2 comentarios más antiguosOcultar -2 comentarios más antiguos

Iniciar sesión para comentar.

Florian Augustin el 15 de Jun. de 2022
Hi Marco,
you may be able to use SimBiology Events to switch between reaction rates. The general idea would be as follows:
• Define a reation rate rate1, rate2, rate3, ... for each of the phases you describe.
• Add dummy (phase-indicator) parameters phase1, phase2, phase3, ... to your model. The initial value of phase1 should be 1, all other parameters should have an initial value of 0.
• Add SimBiology Events to switch between different phases.
Below is an example using SimBiology on the MATLAB Command Window. Dependent on your actual rates you may be able to condense/rewrite the events; you may not even need the phase-indicator parameters. But I hope this demonstrates how you could achieve different reation rates in different time intervals.
% Define model
model = sbiomodel("time-dependent reaction rates");
% Define reaction rate: here I am using rate1 = k*A, rate2 = -k, rate3 = -k*A as an example.
reaction = addreaction(model, "A -> null");
reaction.ReactionRate = "phase1*k*A - phase2*k - phase3*k*A";
% Add events to switch between phases:
addevent(model, "time >= 5" , ["phase1 = 0", "phase2 = 1", "k = 0.1"]);
addevent(model, "time >= 12", ["phase2 = 0", "phase3 = 1", "k = 0.4"]);
% Simulate model
configset = getconfigset(model);
configset.StopTime = 19;
configset.RuntimeOptions.StatesToLog = "A";
simData = sbiosimulate(model);
% Plot results
sbioplot(simData)
Best,
Florian
##### 9 comentariosMostrar 7 comentarios más antiguosOcultar 7 comentarios más antiguos
mitpi_182 el 28 de Jul. de 2022
Hello Florian,
Is there any chance you could show me how to do this using the Simbiology GUI? Specifically, how to integrate the sigmoid function in my model in order to use it to define the reaction rates. Many thanks for your time.
Best,
Marco.
Florian Augustin el 29 de Jul. de 2022
Hi Marco,
attached is an sbproj file with the example above. The project is compatible with MATLAB version R2020b and later. I hope this helps you to get started.
Best,
Florian

Iniciar sesión para comentar.

### Más respuestas (0)

Más respuestas en  SimBiology Community

### Categorías

Más información sobre Import Data en Help Center y File Exchange.

R2020b

### Community Treasure Hunt

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

Start Hunting!

Translated by