630 lines
39 KiB
HTML
630 lines
39 KiB
HTML
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN"
|
|
"http://www.w3.org/TR/REC-html40/loose.dtd">
|
|
<html>
|
|
<meta name="GENERATOR" content="TtH 3.81">
|
|
<style type="text/css"> div.p { margin-top: 7pt;}</style>
|
|
<style type="text/css"><!--
|
|
td div.comp { margin-top: -0.6ex; margin-bottom: -1ex;}
|
|
td div.comb { margin-top: -0.6ex; margin-bottom: -.6ex;}
|
|
td div.hrcomp { line-height: 0.9; margin-top: -0.8ex; margin-bottom: -1ex;}
|
|
td div.norm {line-height:normal;}
|
|
span.roman {font-family: serif; font-style: normal; font-weight: normal;}
|
|
span.overacc2 {position: relative; left: .8em; top: -1.2ex;}
|
|
span.overacc1 {position: relative; left: .6em; top: -1.2ex;} --></style>
|
|
|
|
|
|
|
|
<title> iBioSim Version 2.8 \ Systems Biology Tutorial</title>
|
|
|
|
<h1 align="center">iBioSim Version 2.8 <br />Systems Biology Tutorial </h1>
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<h3 align="center">Chris J. Myers </h3>
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<div class="p"><!----></div>
|
|
<h2><a name="tth_sEc1">
|
|
1</a> Introduction</h2>
|
|
|
|
<div class="p"><!----></div>
|
|
<tt>iBioSim</tt> has been developed for the modeling, analysis, and design of genetic circuits. While <tt>iBioSim</tt> primarily targets models of genetic circuits, models representing metabolic networks, cell-signaling pathways, and other biological and chemical systems can also be analyzed. <tt>iBioSim</tt> includes the following components:
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<ul>
|
|
<li> Model Editor - a tool to create a model of a genetic circuit or other biological system.
|
|
<div class="p"><!----></div>
|
|
</li>
|
|
|
|
<li> Analysis Tool - an abstraction-based ODE, Monte Carlo, and Markov analysis tool.
|
|
<div class="p"><!----></div>
|
|
</li>
|
|
|
|
<li> Learn Tool - a tool to learn a model from <em>time series data</em> (TSD).
|
|
<div class="p"><!----></div>
|
|
</li>
|
|
|
|
<li> TSD Graph Editor- a tool to visualize TSD files.
|
|
<div class="p"><!----></div>
|
|
</li>
|
|
|
|
<li> Histogram Graph Editor - a tool to visualize probability data.
|
|
<div class="p"><!----></div>
|
|
</li>
|
|
</ul>
|
|
|
|
<div class="p"><!----></div>
|
|
This tutorial illustrates each of these features of <tt>iBioSim</tt> using a simple model for the <em>cI</em> and <em>cII</em> genes and the P<sub>R</sub> and P<sub>RE</sub> promoters from the phage λ decision circuit.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<div class="p"><!----></div>
|
|
<h2><a name="tth_sEc2">
|
|
2</a> Project Management</h2>
|
|
|
|
<div class="p"><!----></div>
|
|
Within <tt>iBioSim</tt>, all files are collected within projects. A project is a collection of models, analysis views, learn views, and graphs. As shown below, <tt>iBioSim</tt> displays all project files on the left, the open models, views, and graphs on the right, and a log of all external commands on the bottom. The menu bar is located on the top of the window in the Windows and Linux versions. It is located on the top of the screen in the MacOS version.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/iBioSim.png" alt="screenshots/iBioSim.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
To create a new project, select New → Project from the File menu as shown below. You will then be prompted to browse to a desired location and to give a name to the project directory. Enter the name <tt>SysBioTutorial</tt>. After you do this, click the new button and a new project directory will be created.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/project.png" alt="screenshots/project.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
<h2><a name="tth_sEc3">
|
|
3</a> Model Editor</h2>
|
|
|
|
<div class="p"><!----></div>
|
|
After you have created a project, you can create a new model to add to the project by selecting New → Model from the File menu as shown below. You will then be prompted to enter a model ID. Enter <tt>lambda</tt>. At this point, a Model editor will open in a new tab.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/newModel.png" alt="screenshots/newModel.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/ModelId.png" alt="screenshots/ModelId.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/ModelEditor.png" alt="screenshots/ModelEditor.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Compartments are the membrane-enclosed regions where species can be found and reactions take place. <tt>iBioSim</tt> creates a default compartment initially with the ID of Cell. If you click on the schematic within the Cell compartment, it brings up the compartment editor. Uncheck the "Is Mapped to a Port". This indicates that this compartment should be enclosing this model and not replaced when instantiated in a larger model. Once you press OK, you will notice that the compartment now has rounded corners to indicate that this is membrane enclosed by the compartment Cell.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/compartment.png" alt="screenshots/compartment.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
To add a chemical species, select the Add Species icon and click on the schematic canvas. This will drop a new species with default ID and other values. You may change these defaults by clicking on the selection icon
|
|
, and
|
|
double-clicking on the species to open the Species Editor. In this case, let us change the ID to CI, and click on degrades checkbox. We will leave all other values at their default values. One thing that is important to note is that when this model is analyzed a default degradation reaction will be created which has a rate of 0.0075.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/species.png" alt="screenshots/species.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Add another species for the CI dimer molecule. The next step is to add a complex-formation reaction to convert CI monomers into CI dimers. Select the complex formation icon , highlight the CI species, and, while holding the mouse button, stretch the complex formation arc to the S1 species. Next, edit this species to set its ID to CI2 and select that it degrades. This species is created using a complex-formation reaction with an equilibrium constant of 0.1. Change default to custom for the complex formation equilibrium and set it to this value as shown below.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/species2.png" alt="screenshots/species2.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
If you double-click on the complex formation arc, an influence editor will open which indicates that this is a complex formation arc and the stoichiometry of binding (i.e., the number of molecules of the source species used to construct the sink species) is 2. The default in this case is correct as it does take two molecules of CI to make CI2.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/complex.png" alt="screenshots/complex.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Next, let's add the PR promoter, which initiates transcription of the gene that produces the protein CII. To do this, select the promoter icon and click on the schematic canvas to drop the promoter with a default ID and parameter values. Double-click on the promoter to bring up the promoter editor. Change the ID to PR, customize the RNAP binding equilibrium to be 0.69422, and set the open complex production rate to be 0.014.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/promoter.png" alt="screenshots/promoter.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
The PR promoter is repressed by the CI2 species. To create this relationship, select the repression arc icon
|
|
, highlight the CI2 species, and, while holding the mouse button, stretch the repression arc to the PR promoter. Next, double-click on the repression arc to bring up the influence editor. In this editor, customize the stoichiometry of binding to 1, indicating that just one CI dimer is necessary to repress this promoter, and change the repression binding equilibrium to 0.2165.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/repression.png" alt="screenshots/repression.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
As mentioned earlier, the PR promoter initiates the production of the CII species. Add the CII species following the steps given earlier for adding a species (be sure to mark that it degrades). Then, highlight the PR promoter and, while holding the mouse button, stretch the production arc to the CII species. Note that the icons selected for this are not important because all arcs from promoters to species are always production arcs.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/production.png" alt="screenshots/production.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Finally, CII species activates the production of the CI species from the PRE promoter. Promoters do not need to always be drawn. They can also be implicit on an influence. To add an activation arc with an implicit promoter, select the activation arc icon , highlight the CII species, and, while holding the mouse button, stretch the activation arc to the CI species. This creates not only the influence but also a default promoter. Double-click on the activation arc to bring up the influence editor. In this editor, customize the stoichiometry of binding to 1, indicating that just one CII molecule is necessary to activate this promoter, and change the activation binding equilibrium to 0.00161. Finally, click on the Edit Promoter button and change the ID of this promoter to PRE. Also, customize the RNAP binding equilibrium to be 0.01, the basal production rate to be 0.00004, and the activated production rate to be 0.015.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/activation.png" alt="screenshots/activation.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/activated_promoter.png" alt="screenshots/activated_promoter.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Next, let's go look at some of the other model tabs. Click on the Functions tab. Here one can add/remove/edit definitions for functions. Click on Add Function, and let's create a function to compute the total number of molecules where x is the number of monomers and y is the number of dimers. Let's give this function an ID of total with arguments x and y, and a definition of x + 2 * y.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/function.png" alt="screenshots/function.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Click on the Units tab. Unit definitions allow the user to create custom units for use in models. As an example, click Add Unit and enter the ID nanoMole. A unit is defined using pre-defined base units. Click on Add to List and find mole in the Kind combo box, and change the scale to −9 to define a nanoMole in terms of moles. The exponent and multiplier are real numbers, and the scale is an integer that specifies the relationship between the derived unit and the base unit using the relation below:
|
|
|
|
<br clear="all" /><table border="0" width="100%"><tr><td>
|
|
<table border="0" cellspacing="0" cellpadding="0">
|
|
<tr><td width="50%"></td><td nowrap="nowrap" align="right" colspan="1"><table border="0" cellspacing="0" cellpadding="2"><tr><td nowrap="nowrap" align="left">
|
|
<span class="roman">unit</span> </td></tr></table></td><td nowrap="nowrap" align="left">
|
|
<table border="0" cellspacing="0" cellpadding="2"><tr><td nowrap="nowrap" align="left">
|
|
= </td></tr></table></td><td nowrap="nowrap" align="left">
|
|
<table><tr><td nowrap="nowrap" align="right" colspan="1">(<span class="roman">multiplier</span> * 10<sup><span class="roman">scale</span></sup> * <span class="roman">baseUnit</span>)<sup><span class="roman">exponent</span></sup></td></tr></table></td><td width="50%"></td></tr></table>
|
|
</td></tr></table>
|
|
|
|
|
|
Follow the same steps to create a unit for microLitre.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/units.png" alt="screenshots/units.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Next, click on the Schematic tab. Then, click on the Model button to bring up a window of default units for this model. These are the units to be used for various things when no units are provided. For this model, let's make the Substance Units be in nanoMoles, Time Units be in seconds, and Volume Units be in microLitres. The Extent Units indicate the units of change for reactions. Let's make that also be nanoMoles.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/ModelUnits.png" alt="screenshots/ModelUnits.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Next, click on the Constants tab. This tab includes the default model generation parameters as well as additional global variables. The model generation parameters are used for default values when creating a reaction-based model from our higher level model. These parameters can be edited by selecting them and entering a new initial value. They should not be removed. Note that new global constants can also be added by the user, if desired.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/GCMparam.png" alt="screenshots/GCMparam.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
While you can add a global variable here, the preferred way is to add it on the schematic. Click on the Schematic tab. Select the add variable icon , and click a location to add the variable to the schematic. Then, click on the selection icon
|
|
, and double click on the variable to open the Parameter Editor. Change the ID to CI_total, for the total amount of the species CI in both monomer and dimer form. Its units are in nanoMole and it is not constant.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/parameter.png" alt="screenshots/parameter.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<div class="p"><!----></div>
|
|
Now, let's create a rule to compute the value of CI_total. There are three types of rules: algebraic, assignment, and rate. Algebraic rules are used to specify relationships that must be maintained. Assignment rules are used to define one variable in terms of a mathematical expression. Finally, rate rules are used to indicate a differential equation to govern the evolution of a variable in terms of a mathematical expression on other variables. To add a rule, select the add rule icon
|
|
. In the Rule editor,
|
|
select Assignment, select the variable CI_total, and enter the expression total(CI,CI2). This uses our new function to compute that the amount of CI in total is the number of monomer molecules plus two times the number of dimer molecules.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/rule.png" alt="screenshots/rule.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<div class="p"><!----></div>
|
|
Next, let's add a constraint. A constraint is a condition that must be satisfied or simulation should terminate. To add a constraint, select the add constraint icon
|
|
. Add a constraint that states that CI remains less than or equal to 20 molecules as shown below. Add a second constraint that states that CII remains less than or equal to 50 molecules.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/constraint.png" alt="screenshots/constraint.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Finally, let's add an event. Events are used to specify discrete state changes. For example, let's describe an event for cell division. Click on the Add Event button, enter the ID CellDivision, trigger true, and delay of 2100. The trigger specifies the condition that should be satisfied to enable this event. In this case, we want it to be enabled initially. The delay indicates the time after which it is enabled that it should execute. This combination results in a cell division event at 2100 seconds after simulation begins. The event assignments specify the state change(s) for this event. Click on the Add Assignment button and select the CI variable, and enter the assignment CI / 2. In other words, the number of molecules is being divided between the two daughter cells with a random distribution. Similarly, add event assignments to divide up CI2 and CII.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/event.png" alt="screenshots/event.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
At this point, you should have a model that looks like the one below (though locations of elements may be different). Now, let's make sure the model is saved by either clicking on the Save icon or selecting the Save option from the File menu.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/save.png" alt="screenshots/save.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
<h2><a name="tth_sEc4">
|
|
4</a> ODE Analysis</h2>
|
|
|
|
<div class="p"><!----></div>
|
|
This section describes how to analyze the model just created using <em>ordinary differential equation</em> (ODE) simulation. The first step is to create an analysis view. To do this, right-click on the model file and select Create Analysis View. Enter the analysis ID <tt>lambda</tt> or just press enter. At this point, a new analysis view should open. You should also notice that an icon appears next to your model file. When you click on this, it will show you all of the analysis and learn views associated with this model.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/GCMAnalysis.png" alt="screenshots/GCMAnalysis.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Next, click on the Schematic tab. This tab allows you to select which SBML model elements to include in your analysis, modify initial conditions, and visualize simulation results on your schematic. For now, let us remove both constraints and the event from our analysis by clicking on them to make them invisible.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/SBMLElements.png" alt="screenshots/SBMLElements.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Now, go back to the simulation options tab. Make sure that the simulation type is set to ODE and the simulator selected is Runge-Kutta-Fehlberg. Change the time limit to 2100.0, change the print interval to 10.0, and enter a Simulation ID of <tt>ode</tt>. Then, either press the Save and Run icon or select the Save and Run option from the File menu.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/analysisView.png" alt="screenshots/analysisView.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
After the simulation completes, click on the TSD Graph tab. Double-click on the graph to bring up the graph editor. Open the <tt>ode</tt> simulation, highlight Average, select CI_total and CII, change the Title to "ODE Simulation Results", change the X-Axis Label to "Time (seconds)", and change the Y-Axis Label to "Number of Molecules".
|
|
Press the OK button.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/odeResults.png" alt="screenshots/odeResults.png" /><br />
|
|
<img src="screenshots/odeSimResults.png" alt="screenshots/odeSimResults.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Graphs can be exported in a variety of formats including:
|
|
|
|
<ul>
|
|
<li> Time series data format (tsd).
|
|
<div class="p"><!----></div>
|
|
</li>
|
|
|
|
<li> Comma separated value (csv).
|
|
<div class="p"><!----></div>
|
|
</li>
|
|
|
|
<li> Column separated data (dat).
|
|
<div class="p"><!----></div>
|
|
</li>
|
|
|
|
<li> Encapsulated postscript (eps).
|
|
<div class="p"><!----></div>
|
|
</li>
|
|
|
|
<li> Joint Photographic Experts Group (jpg).
|
|
<div class="p"><!----></div>
|
|
</li>
|
|
|
|
<li> Portable document format (pdf).
|
|
<div class="p"><!----></div>
|
|
</li>
|
|
|
|
<li> Portable network graphics (png).
|
|
<div class="p"><!----></div>
|
|
</li>
|
|
|
|
<li> Scalable vector graphics (svg).
|
|
<div class="p"><!----></div>
|
|
</li>
|
|
</ul>
|
|
In order to export a graph, you can either click on the Export icon or select one of the graph export options from the File menu. When using the Export icon, the type of file exported will depend on the extension provided to the file name. Click on the Export icon, browse to a location on your file system, and enter the file name of <tt>ode.pdf</tt> to create a PDF file for your graph.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/exportTSD.png" alt="screenshots/exportTSD.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Genetic production and its regulation are modeled as a single complex reaction. If the expand reactions abstraction option is selected, the binding reactions for transcription factors are separated into their own reactions providing a more detailed model. To see the effect of this, select Expand Reactions, change the Simulation ID to <tt>ode2</tt> and rerun the simulation. Then, go to the TSD Graph tab, and add CI_total and CII from this simulation to the graph. You should see that there is not a significant difference though the simulation time is slower due to these extra reactions. If you would like to see the actual model that is simulated in this case, you can save the reaction-based model by selecting Model as your simulation type. Then, either press the Save and Run icon or select the Save and Run option from the File menu. In this case, you must provide a new model ID. This new model will appear in your project and it can be opened in the Model Editor. Since this model does not include any layout information, you will need to either lay it out by hand or using one of the default layout routines selectable using the Apply Layout icon ,
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/reactionModel.png" alt="screenshots/reactionModel.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
There are two additional ways to see reaction-based models. If GraphViz is installed on your computer, you can select Network for your Simulation Type to open the reaction-based model in GraphViz (the example below is when abstraction is set to none). If it does not open in GraphViz, make sure that you have files with the <tt>.dot</tt> file extension associated with GraphViz on your computer.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/viewNetwork.png" alt="screenshots/viewNetwork.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
You can also view the model in a web browser by selecting Browser for your simulation type. In this case, you should ensure that you have files with the <tt>.xhtml</tt> extension associated with your favorite browser.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/viewBrowser.png" alt="screenshots/viewBrowser.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Parameters and initial amounts can be modified for a simulation. To try this, select the schematic tab, and click on CII to bring up the parameter editor. Change the degradation rate from default to modified, and change its value to 0.075. Go back to the simulation options tab, select ODE simulation, change abstraction back to none, and change the simulation id to <tt>ode3</tt>. After running simulation, add CI_total and CII from the <tt>ode3</tt> simulation results to the graph. You should observe that the amounts have reduced significantly.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/paramEdit.png" alt="screenshots/paramEdit.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Go back to the schematic and change the degradation rate for CII back to default. Now, click on the Parameters tab. Here, click on kd to open the parameter editor for this constant. Let us sweep this parameter. Change the value type to Modified, and click on the sweep button. Enter a start value of 0.005, stop value of 0.01, and step value of 0.0025. Go back to the simulation tab, and remove the simulation id, and run the simulation. Now, go to the TSD graph tab, and edit the graph. First, press the Deselect All button to clear the selections. You should now have three new simulation folders, one for each value. Go into each and select CI_total from the average for each. This will plot the value of CI_total for each degradation rate.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/parameterEditor.png" alt="screenshots/parameterEditor.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
<h2><a name="tth_sEc5">
|
|
5</a> Stochastic Analysis</h2>
|
|
|
|
<div class="p"><!----></div>
|
|
This section describes how to analyze the model using stochastic, Monte Carlo, methods. First, go back to the parameter tab and change kd back to Original. Then, select the Simulation Options tab again, select <tt>Monte Carlo</tt>, change the number of runs to 100, set the simulation ID to <tt>ssa</tt>, and click on the Save and Run icon. Click on the TSD Graph tab. Double-click on the graph to bring up the graph editor. Open the <tt>ssa</tt> simulation directory, and highlight <tt>run-1</tt>. Select CI_Total and CII, change the title to "SSA Simulation Results", change the X-Axis Label to "Time (seconds)", and change the Y-Axis Label to "Number of Molecules". Press the OK button. Click on the Export icon and enter the file name <tt>ssa-1.pdf</tt>. Repeat these steps to generate graphs for the average (<tt>average.pdf</tt>) and standard deviation (<tt>stddev.pdf</tt>). Note that you can use the "Deselect All" button to remove all items from the graph.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/ssaSim.png" alt="screenshots/ssaSim.png" /><br />
|
|
<img src="screenshots/ssaResults.png" alt="screenshots/ssaResults.png" /><br />
|
|
<img src="screenshots/ssaSimResults.png" alt="screenshots/ssaSimResults.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Another way to view simulation results is on the schematic. To do this, click on the schematic tab. At the bottom of the window, select the Choose Simulation button, which brings up a window with all the simulations in this analysis view. Open the <tt>ssa</tt> directory, select <tt>run-1.tsd</tt>, and press OK.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/chooseSim.png" alt="screenshots/chooseSim.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Now, click on the CI species, which brings up the Edit Species window. Select the Appearance tab. Here you can select how you want the species to appear as you playback the simulation. You can have it change color, size, and/or opacity on a gradient. You can also select the range of molecule counts to specify the ends of the gradient(s). Finally, you can indicate that these selections are either for this species or all species in the model. For our example, let's make CI follow a green color gradient, CI2 follow a red color gradient, and CII follow a blue color gradient.
|
|
|
|
<center>
|
|
<img src="screenshots/editSpeciesAppearance.png" alt="screenshots/editSpeciesAppearance.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Once you have made your selections, you can now play back the simulation. You can either single-step the simulation by pressing the icon or play continuously by pressing the icon. The playback can also be paused by pressing the icon and restarted by pressing the icon.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/movieView.png" alt="screenshots/movieView.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
<h2><a name="tth_sEc6">
|
|
6</a> Reaction-based Abstraction</h2>
|
|
|
|
<div class="p"><!----></div>
|
|
The efficiency of simulation can be improved by employing various automatic abstraction techniques. Go back to the Schematic tab and click on the Total_CI rule to exclude it from the analysis. To activate abstraction, click on the Simulation Options tab, select Reaction-based Abstraction and change the simulation ID to <tt>abs</tt>. Press the Save and Run icon and note that the simulation time is substantially faster. Plot both the SSA results for CI_total and CII with the abstraction results for CI (note this is now equivalent to CI_total after abstraction) and CII.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/absResults.png" alt="screenshots/absResults.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
One way to understand why abstraction is so much faster is by looking at the complexity of the reaction-based model before and after abstraction. To see the difference, select the Network option, then generate the graph for None, Expand Reactions, and Reaction-based Abstraction options. The reaction-based model after abstraction is shown below which is clearly much simpler than the full model shown earlier.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/viewNetworkAbs.png" alt="screenshots/viewNetworkAbs.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
<h2><a name="tth_sEc7">
|
|
7</a> State-based Abstraction</h2>
|
|
|
|
<div class="p"><!----></div>
|
|
Next, let's try checking some properties. To do this, go to the SBML Elements tab and check the boxes next to the constraints. Recall that these constraints terminate simulation whenever CI goes above 20 molecules or CII goes above 50 molecules. Go back to the Simulation Options tab and change abstraction back to none, the Simulation Type to Monte Carlo, and Simulation ID to prob, then press the Save and Run icon. Now, let's plot the results on a histogram by clicking on the Histogram tab and then double-clicking on the graph to bring up the histogram graph editor shown below. Open the prob folder, select the sim-rep file, and check the Use check box to get all fields.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/editProbGraph.png" alt="screenshots/editProbGraph.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
The histogram shown here indicates that CI goes above 20 molecules first about 21 percent of the time, CII goes above 50 molecules first about 74 percent of the time, and the simulation terminates before either happens about 5 percent of the time.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center> <img src="screenshots/probResults.png" alt="screenshots/probResults.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
These results can also be visualized using the TSD graph tool. Click on the TSD graph tab, click on the graph, Deselect All, open the prob folder, select the Percent Termination file, and add both constraints to the graph. The result, shown below is the probability of each constraint terminating the simulation as time evolves.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center> <img src="screenshots/probResultsTSD.png" alt="screenshots/probResultsTSD.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
<h2><a name="tth_sEc8">
|
|
8</a> Learn Tool</h2>
|
|
|
|
<div class="p"><!----></div>
|
|
This section describes how a model can be learned from time series data using <tt>iBioSim</tt>'s Learn Tool. To demonstrate the Learn Tool, first create a simple model, <tt>learnModel</tt>, which just includes the two species CI and CII as shown below. This model represents any background knowledge you have about the model that you wish to learn. Any influences that you add between species will be assumed as part of the final learned model. You can also indicate definite knowledge about there not being any influence between two species by selecting the No Influence icon
|
|
and connecting the species with a no influence arc.
|
|
|
|
<div class="p"><!----></div>
|
|
Now, create a learn view by right-clicking on this model file and selecting Create Learn View. Give this learn view the ID <tt>learnLambda</tt>. At this point, a new learn view should open. You should also notice that an icon appears next to your model file. When you click on this, it will show you all of the analysis and learn views associated with this model.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/createLearn.png" alt="screenshots/createLearn.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
The next step is to add some experimental data from which you wish to learn a model. In this demo, we will just utilize our simulation data as synthetic experimental data. To do this, click Copy From View, and select <tt>lambda/abs</tt>. Highlight <tt>lambda/abs/run-1.tsd</tt> and you should see the simulation data for CI and CII appear on the right in the data editor.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/dataManager.png" alt="screenshots/dataManager.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Now, click on the Learn tab. Here you can edit the various learning options. For example, you can either use auto-generated levels or user-generated levels for your data encoding. Select Use User Generated Levels, which will make the levels below editable. At this point, you can ask the tool to suggest levels by clicking on the Suggest Levels button. Finally, click on the Save and Run icon, which if successful, will ask you for a model ID for the generated model. Enter <tt>learnedModel</tt>, then open this model in the model editor which should show you the model shown below.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/learn.png" alt="screenshots/learn.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/learnResult.png" alt="screenshots/learnResult.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
<h2><a name="tth_sEc9">
|
|
9</a> Advanced Modeling</h2>
|
|
|
|
<div class="p"><!----></div>
|
|
This section is less detailed than the others but it gives some intuition about modeling using reactions, components, and grids.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<div class="p"><!----></div>
|
|
First, let's consider an alternative model of CII degradation which we are going to model using chemical reactions. To do this, create a new model named <tt>CII_degradation</tt>. In this model, create species CII, P1, and CII_P1, making CII have the input type so that we can connect to it later. Set an initial amount of 35 molecules for P1.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/P1.png" alt="screenshots/P1.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Now, select the Add Reaction icon and click on the schematic canvas to drop a reaction. This creates a reaction with a default ID and parameter values that we can change later, if we wish.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/reaction.png" alt="screenshots/reaction.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Now, let's connect up the reactant species. To do this, select the Reaction icon , select the reactant species CII, and, while holding the mouse button, drag the reaction edge to the reaction R1. Similarly, add P1 as a reactant as well. If you double-click on a reactant edge, it brings up a Reactant editor where you can change the stoichiometry, if desired.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/reactant.png" alt="screenshots/reactant.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Connecting product species is accomplished in much the same way, except in this case you select reaction R1 and drag the reaction edge to the product CII_P1. Again, there is a Product editor for changing the stoichiometry. Note that modifiers (i.e., species that are neither produced nor consumed by a reaction but simply catalyze a reaction) can be added in a similar way using the Modifier icon instead.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/product.png" alt="screenshots/product.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Now, let's adjust the parameters for this reaction by clicking on it to open the Reaction Editor. Press the Use Mass Action button to automatically create a kinetic law for this reaction. Then, make this reaction reversible and adjust its forward reaction rate to be 1.0.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/localParam.png" alt="screenshots/localParam.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Follow the same steps to add another reaction that degrades CII in the CII_P1 form and releases the protease molecule P1. This reaction is not reversible and it should have a forward rate of 0.002.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/kineticLaw.png" alt="screenshots/kineticLaw.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Let's now go and add this new degradation mechanism to our lambda model (you might actually want to copy your old model before you do this, which you can do by highlighting the file and selecting Edit → Copy or using the right mouse button menu). To simplify things, remove the rule, constraints, and event. Next, open the Species Editor on CII and deselect the degrades option. Finally, select the Add Component icon and click on the Schematic canvas opening the Add Component(s) window. In this window, browse the combo box to find your CII_degradation model. Pressing OK will then add it to your schematic.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/addComponent.png" alt="screenshots/addComponent.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
Now, let's connect CII to this new component to relate the CII within the component to the outer CII species. To do this, select the CII species and, while holding the mouse button, drag a connection to the component connecting CII to the CII port on the component.
|
|
You may want to now go and try simulating this model, if you like.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/addPort.png" alt="screenshots/addPort.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
In the last example, we will build a model with a grid. First, edit the CII species and make it diffusible, and save the model. Now, create a grid model using the File → New → Grid Model menu, and name the new model <tt>Population</tt>. In the create grid window shown below, select your copy of your lambda model and change the number of rows and columns to 5.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/createGrid.png" alt="screenshots/createGrid.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
The schematic in a grid model is a bit different. It includes a grid in which each location can be empty or contain exactly one component. Only components can be added to grids.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/gridModel.png" alt="screenshots/gridModel.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
When you create a reaction-based model for a grid during analysis, reactions are created to move the diffusible species between the grid locations to provide a coarse form of spatial modeling. If the component within a grid location is enclosed in a compartment membrane (indicated by the rounded corners), the model generated also includes reactions to diffuse the species in and out of the compartment. In the analysis schematic, you can visualize your grid models by clicking on the component in the grid and selecting the species that you would like to see. For each such species, you can set its color, size, and/or opacity gradient. You can also copy these settings to all like models in your grid. Finally, you can click on the area outside of the component within the grid to allow you to also visualize the species that are in the medium.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<center><img src="screenshots/compAppearances.png" alt="screenshots/compAppearances.png" />
|
|
</center>
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<div class="p"><!----></div>
|
|
A more detailed
|
|
<a href="iBioSim_Grid_Tutorial.html">
|
|
grid tutorial
|
|
</a>
|
|
is available in the <tt>docs</tt> directory that comes with the distribution.
|
|
|
|
<div class="p"><!----></div>
|
|
|
|
<br /><br /><hr /><small>File translated from
|
|
T<sub><font size="-1">E</font></sub>X
|
|
by <a href="http://hutchinson.belmont.ma.us/tth/">
|
|
T<sub><font size="-1">T</font></sub>H</a>,
|
|
version 3.81.<br />On 11 Sep 2017, 09:31.</small>
|
|
</html>
|