- 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.

# Time-dependent Growth Rate of Species in Model

11 views (last 30 days)

Show older comments

mitpi_182
on 15 Jun 2022

Edited: Florian Augustin
on 29 Jul 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?

Many thanks in advance for any suggestions or comments.

##### 0 Comments

### Accepted Answer

Florian Augustin
on 15 Jun 2022

Hi Marco,

you may be able to use SimBiology Events to switch between reaction rates. The general idea would be as follows:

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");

addspecies(model, "A", 10);

addparameter(model, "k", 1, "Constant", false);

% Add phase-indicator (dummy) parameters:

addparameter(model, "phase1", 1, "Constant", false);

addparameter(model, "phase2", 0, "Constant", false);

addparameter(model, "phase3", 0, "Constant", false);

% 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 Comments

Florian Augustin
on 29 Jul 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

### More Answers (0)

### Communities

More Answers in the SimBiology Community

### See Also

### Categories

### Community Treasure Hunt

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

Start Hunting!