I've used Matlab extensively to model biochemical networks with systems of ODEs, though I'm new to using SimBiology. I've imported a calcium channel model into SimBiology and can execute it, plot the results, etc.
What I'd like to do now is take one of the species, Ca — currently with its InitialAmount set to 0.5 millimole/litre, and its BoundaryCondition and ConstantAmount attributes set to on — and control its concentration as a function of time, treating it as an input function to the rest of the system. (I'm investigating how different ways of pulsing Ca concentration will affect downstream production of autonomous CaMKII.)
What is the best way to accomplish building such an input function in SimBiology?
I found this article, which describes how to create a custom user-defined function and use it in reaction, rule, and event expressions. But I don't think this is what I need to do.
I'd appreciate any guidance. Thanks.
Best Answer