- 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)
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.
Florian Augustin on 15 Jun 2022
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
More Answers (0)
More Answers in the SimBiology Community