The ReactionRate property defines the reaction rate equation. You can define a ReactionRate with or without the KineticLaw property. KineticLaw defines the type of reaction rate. The addkineticlaw function configures the ReactionRate based on the KineticLaw and the species and parameters specified in the kinetic law object properties SpeciesVariableNames and ParameterVariableNames.

The reaction takes place in the reverse direction if the Reversible property is true. This is reflected in ReactionRate. The ReactionRate includes the forward and reverse rate if reversible.

You can specify ReactionRate without KineticLaw. Use the set function to specify the reaction rate equation. SimBiology software adds species variables while creating reactionObj using the addreaction method. You must add the parameter variables (to the modelObj in this case). See the example below.

After you specify the ReactionRate without KineticLaw and you later configure the reactionObj to use KineticLaw, the ReactionRate is unset until you specify SpeciesVariableNames and ParameterVariableNames.

See How Reaction Rates Are Evaluated in the SimBiology User's Guide for information on dimensional analysis for reaction rates.

Note If you set the ReactionRate property to an expression that is not continuous and differentiable, see Using Events to Address Discontinuities in Rule and Reaction Rate Expressions before simulating your model.

Characteristics

Applies to	Object: reaction
Data type	`char string`
Data values	Reaction rate string. Default is `''` (empty).
Access	Read/write

Examples

Example 1

Create a model, add a reaction, and assign the expression for the reaction rate equation.

Create a model object, and then add a reaction object.

modelObj = sbiomodel('my_model');
reactionObj = addreaction(modelObj, 'a -> c + d');

Create a kinetic law object for the reaction object of the type 'Henri-Michaelis-Menten'.
```
kineticlawObj = addkineticlaw(reactionObj, 'Henri-Michaelis-Menten');
```
reactionObj KineticLaw property is configured to kineticlawObj.
The 'Henri-Michaelis-Menten' kinetic law has two parameter variables (Vm and Km) and one species variable (S) that you should set. To set these variables, first create the parameter variables as parameter objects (parameterObj1, parameterObj2) with names Vm_d and Km_d and assign them to kineticlawObj.
```
parameterObj1 = addparameter(kineticlawObj, 'Vm_d');
parameterObj2 = addparameter(kineticlawObj, 'Km_d');
```

Set the variable names for the kinetic law object.

set(kineticlawObj,'ParameterVariableNames', {'Vm_d' 'Km_d'});
set(kineticlawObj,'SpeciesVariableNames', {'a'});

Verify that the reaction rate is expressed correctly in the reaction object ReactionRate property.
```
get (reactionObj, 'ReactionRate')
```
MATLAB returns:
```
ans =

Vm_d*a/(Km_d + a)
```

Example 2

Create a model, add a reaction, and specify ReactionRate without a kinetic law.

Create a model object, and then add a reaction object.

modelObj = sbiomodel('my_model');
reactionObj = addreaction(modelObj, 'a + b -> c + d');

Specify ReactionRate and verify the assignment.

set (reactionObj, 'ReactionRate', 'k*a');
get(reactionObj, 'ReactionRate')

MATLAB returns:

ans =

k*a

You cannot simulate the model until you add the parameter k to the modelObj.
```
parameterObj = addparameter(modelObj, 'k');
```
SimBiology adds the parameter to the modelObj with default Value = 1.0 for the parameter.

Free Computational Biology Interactive Kit

See how to analyze, visualize, and model biological data and systems using MathWorks products.

Get free kit

Trials Available

Try the latest computational biology products.

Get trial software

ReactionRate - Reaction rate equation in reaction object

Description