SimBiology®

KineticLaw object - Kinetic law information for reaction

Description

The kinetic law object holds information about the abstract kinetic law applied to a reaction and provides a template for the reaction rate. In the model, the SimBiology software uses the information you provide in a fully defined kinetic law object to determine the ReactionRate property in the reaction object.

When you first create a kinetic law object, you must specify the name of the abstract kinetic law to use. The SimBiology software fills in the KineticLawName property and the Expression property in the kinetic law object with the name of the abstract kinetic law you specified and the mathematical expression respectively. The software also fills in the ParameterVariables property and the SpeciesVariables property of the kinetic law object with the values found in the corresponding properties of the abstract kinetic law object.

To obtain the reaction rate, you must fully define the kinetic law object:

1. In the ParameterVariableNames property, specify the parameters from the model that you want to substitute in the expression (Expression property).

2. In the SpeciesVariableNames property, specify the species from the model that you want to substitute in the expression.

   The SimBiology software substitutes in the expression, the names of parameter variables and species variables in the order specified in the ParameterVariables and SpeciesVariables properties respectively.

   The software then shows the substituted expression as the reaction rate in the ReactionRate property of the reaction object. If the kinetic law object is not fully defined, the ReactionRate property remains ' ' (empty).

For links to kinetic law object property reference pages, see Property Summary.

Properties define the characteristics of an object. Use the get and set commands to list object properties and change their values at the command line. You can interactively change object properties in the SimBiology desktop.

For an explanation of how relevant properties relate to one another, see Command Line.

The following sections use a kinetic law example to show how you can fully define your kinetic law object to obtain the reaction rate in the SimBiology desktop and at the command line.

The Henri-Michaelis-Menten kinetic law is expressed as follows:

In the SimBiology software Henri-Michaelis-Menten is a built-in abstract kinetic law, where V_m and K_m are defined in the ParameterVariables property of the abstract kinetic law object, and S is defined in the SpeciesVariables property of the abstract kinetic law object.

SimBiology Desktop

To fully define kinetic law, in the SimBiology desktop, define the names of the species variables and parameter variables that participate in the reaction rate in the Project Settings-Reactions pane on the Kinetic Law tab. To add a reaction and set the reaction rate in the SimBiology desktop, see "Adding Reactions to a Model" in the SimBiology Getting Started Guide documentation.

Command Line

To fully define the kinetic law object at the command line, define the names of the parameters in the ParameterVariableNames property of the kinetic law object, and define the species names in the SpeciesVariableNames property of the kinetic law object. For example, to apply the Henri-Michaelis-Menten abstract kinetic law to a reaction

  A -> B
  where Vm = Va, Km = Ka
  and S = A

Define Va and Ka in the ParameterVariableNames property to substitute the variables that are in the ParameterVariables property (Vm and Km). Define A in the SpeciesVariableName property to be used to substitute the species variable in the SpeciesVariables property (S). Specify the order of the model parameters to be used for substitution in the same order that the parameter variables are listed in the ParameterVariables property. Similarly, specify species order if more than one species variable is represented.

% Find the order of the parameter variables 
% in the kinetic law expression. 

get(kineticlawObj, 'ParameterVariables')

ans = 

    'Vm'    'Km'

% Find the species variable in the
% kinetic law expression

get(kineticlawObj, 'SpeciesVariables')
ans = 

    'S'

% Specify the parameters and species variables 
% to be used in the substitution. 
% Remember to specify order, for example Vm = Va 
% Vm is listed first in 'ParameterVariables', 
% therefore list Va first in 'ParameterVariableNames'.

set(kineticlawObj,'ParameterVariableNames', {'Va' 'Ka'});
set(kineticlawObj,'SpeciesVariableNames', {'A'});

The rate equation is assigned in the reaction object as follows:

For a detailed procedure, see Examples.

The following table summarizes the relationships between the properties in the abstract kinetic law object and the kinetic law object in the context of the above example.

Henri-Michaelis
-Menten

Vm and Km

S

Va and Ka

A

Constructor Summary

Method Summary

Property Summary

Examples

This example shows how to define the reaction rate for a reaction.

1. Create a model object, and add a reaction object to the model.

   modelObj = sbiomodel ('my_model');
   reactionObj = addreaction (modelObj, 'A -> B');

2. Define a kinetic law for the reaction object.

   kineticlawObj = addkineticlaw(reactionObj, 'Henri-Michaelis-Menten');
   

3. Query the parameters and species variables defined in the kinetic law.

   get(kineticlawObj, 'ParameterVariables')
   
   ans = 
   
       'Vm'    'Km'
   
   get(kineticlawObj, 'SpeciesVariables')
   ans = 
   
       'S'
   

4. Define Va and Ka as ParameterVariableNames, which correspond to the ParameterVariables Vm and Km. To set these variables, first create the parameter variables as parameter objects (parameterObj1, parameterObj2) with the names Va and Ka, and then add them to kineticlawObj. The species object with Name A is created when reactionObj is created and need not be redefined.

   parameterObj1 = addparameter(kineticlawObj, 'Va');
   parameterObj2 = addparameter(kineticlawObj, 'Ka');

5. Set the variable names for the kinetic law object.

   set(kineticlawObj,'ParameterVariableNames', {'Va' 'Ka'});
   set(kineticlawObj,'SpeciesVariableNames', {'A'});

6. Verify that the reaction rate is expressed correctly in the reaction object ReactionRate property.

   get (reactionObj, 'ReactionRate')
   

   MATLAB returns:

   ans =
   
   Va*A/(Ka+A)

See Also

AbstractKineticLaw object, Configset object, Model object, Parameter object, Reaction object, Root object, Rule object, Species object

SimBiology property Expression

getvariant (model)

Model object