Difference between revisions of "Optimization examples"
(→Testing the methods convergence rate for models with the different number of parameters) |
|||
| (26 intermediate revisions by one user not shown) | |||
| Line 1: | Line 1: | ||
| + | <font size=3> | ||
| + | |||
Here we give some examples of the [[BioUML]] usage for solving the problem of parameter estimation applied to the models of biochemical pathways. | Here we give some examples of the [[BioUML]] usage for solving the problem of parameter estimation applied to the models of biochemical pathways. | ||
For details about creation your oun optimization document in BioUML, see the chapter [[Optimization document]]. All information about the optimization methods implemented in BioUML is done in the chapter [[Optimization problem]]. | For details about creation your oun optimization document in BioUML, see the chapter [[Optimization document]]. All information about the optimization methods implemented in BioUML is done in the chapter [[Optimization problem]]. | ||
| Line 9: | Line 11: | ||
<li>'''Experimental data''': ''data'' > ''Examples'' > ''Optimization'' > ''Data'' > ''Experiments'' > ''exp_data_1''</li> | <li>'''Experimental data''': ''data'' > ''Examples'' > ''Optimization'' > ''Data'' > ''Experiments'' > ''exp_data_1''</li> | ||
</ul> | </ul> | ||
| + | |||
To analyze a convergence rate of the optimization methods implemented in [[BioUML]] [1], we considered a reaction chain extracted from the model by Neumann et al. [2] and representing activation of caspase-8 triggered by the receptor | To analyze a convergence rate of the optimization methods implemented in [[BioUML]] [1], we considered a reaction chain extracted from the model by Neumann et al. [2] and representing activation of caspase-8 triggered by the receptor | ||
| Line 238: | Line 241: | ||
</table> | </table> | ||
| − | ==Testing the | + | ==Testing the computational speed of the optimization methods== |
<ul> | <ul> | ||
| − | <li>'''Optimization documents''': ''data'' > ''Examples'' > ''Optimization'' > ''Data'' > ''Documents'' ( | + | <li>'''Optimization documents''': ''data'' > ''Examples'' > ''Optimization'' > ''Data'' > ''Documents'' <div></div> (test cases 1A, 1B, 1C, 2, 3) </li> |
| − | <li>'''Models''': ''data'' > ''Examples'' > ''Optimization'' > ''Data'' > ''Diagrams'' ( | + | <li>'''Models''': ''data'' > ''Examples'' > ''Optimization'' > ''Data'' > ''Diagrams'' <div></div> (diagrams 1A, 1B, 1C, 2, 3 for the corresponding test cases) </li> |
| − | <li>'''Experimental data''': ''data'' > ''Examples'' > ''Optimization'' > ''Data'' > ''Experiments'' (''exp_data_1'' for the test cases 1A, 1B, 1C; ''exp_data_2'' for the test case 2; ''exp_data_3'' for the test case 3)</li> | + | <li>'''Experimental data''': ''data'' > ''Examples'' > ''Optimization'' > ''Data'' > ''Experiments'' <div></div> (''exp_data_1'' for the test cases 1A, 1B, 1C; ''exp_data_2'' for the test case 2; ''exp_data_3'' for the test case 3) </li> |
</ul> | </ul> | ||
| + | |||
| + | |||
| + | We tested a computational speed of such optimization methods in BioUML as particle swarm optimization, adaptive simulated annealing, stochastic ranking evolution strategy (SRES), and cellular genetic algorithm. | ||
| + | For this purpose, we used biochemical models with the different number of parameters and species introduced in the following test cases. | ||
| + | Firstly, we derived three models of CD95-induced activation of caspase-8 from the model by Neumann et al. [2] with varying degrees of detail (the test cases 1A, 1B and 1C). | ||
| + | Secondly, we took the test case proposed by Mendes et al. [3] for the MAP kinase cascade model developed by Kholodenko et al. [4] (the test case 2). | ||
| + | Finally, we tested the model by Bagci et al. [5] representing the mitochondria-depended apoptosis (the test case 3). | ||
| + | |||
| + | As expected, the computational speed for all test cases directly depended on the number of parameters and species in the model. | ||
| + | The greatest computational speed was shown by the method of particle swarm optimization. | ||
| + | The other methods showed about the same speed for running in one core, wherease, for running in several cores, the computational speed of PSO, SRES and MOCell was evidently higher compared to the simulated annealing. | ||
| + | |||
| + | <font size=2> | ||
| + | <gallery> | ||
| + | File:optimization_examples_model_1.png|'''Test case 1A''': 7 species and 5 reaction rate parameters. | ||
| + | File:optimization_examples_model_2.png|'''Test case 1B''': 13 species and 10 reaction rate parameters. | ||
| + | File:optimization_examples_model_3.png|'''Test case 1C''': 18 species and 12 reaction rate parameters. | ||
| + | File:optimization_examples_model_4.png|'''Test case 2''': 8 species and 22 reaction rate parameters. | ||
| + | File:optimization_examples_model_5.png|'''Test case 3''': 32 species and 57 reaction rate parameters. | ||
| + | File:optimization_examples_figure_3.png|The number of objective function evaluations per second for the different test cases in BioUML. | ||
| + | </gallery> | ||
| + | </font> | ||
==References== | ==References== | ||
| − | # Kutumova E., Ryabova A., Valeev T., Kolpakov F. BioUML plug-in for nonlinear parameter estimation using multiple experimental data. ''Virtual | + | # Kutumova E., Ryabova A., Valeev T., Kolpakov F. BioUML plug-in for nonlinear parameter estimation using multiple experimental data. ''Virtual Biology''. 2013. 1:47-58. |
# Neumann L., Pforr C., Beaudouin J., Pappa A., Fricker N., Krammer P.H., Lavrik I.N., Eils R. Dynamics within the CD95 death-inducing signaling complex decide life and death of cells. ''Molecular Systems Biology''. 2010. 6:352. | # Neumann L., Pforr C., Beaudouin J., Pappa A., Fricker N., Krammer P.H., Lavrik I.N., Eils R. Dynamics within the CD95 death-inducing signaling complex decide life and death of cells. ''Molecular Systems Biology''. 2010. 6:352. | ||
| + | # Mendes P., Hoops S., Sahle S., Gauges R., Dada J., Kummer U. Computational modeling of biochemical networks using COPASI. ''Methods in Molecular Biology''. 2009. 500:17–59. | ||
| + | # Kholodenko B.N. Negative feedback and ultrasensitivity can bring about oscillations in the mitogenactivated protein kinase cascades. ''European Journal of Biochemistry''. 2000. 267(6):1583–1588. | ||
| + | # Bagci E.Z., Vodovotz Y., Billiar T.R., Ermentrout G.B., Bahar I. Bistability in apoptosis: roles of bax, bcl-2, and mitochondrial permeability transition pores. ''Biophysical Journal''. 2006. 90(5):1546–1559. | ||
| + | |||
| + | </font> | ||
Latest revision as of 12:07, 16 March 2022
Here we give some examples of the BioUML usage for solving the problem of parameter estimation applied to the models of biochemical pathways. For details about creation your oun optimization document in BioUML, see the chapter Optimization document. All information about the optimization methods implemented in BioUML is done in the chapter Optimization problem.
[edit] Testing the convergence rate of the optimization methods
- Optimization document: data > Examples > Optimization > Data > Documents > test_case_1A
- Model: data > Examples > Optimization > Data > Diagrams > diagram_1A
- Experimental data: data > Examples > Optimization > Data > Experiments > exp_data_1
To analyze a convergence rate of the optimization methods implemented in BioUML [1], we considered a reaction chain extracted from the model by Neumann et al. [2] and representing activation of caspase-8 triggered by the receptor
CD95 (APO-1/Fas).
 |
|
We performed estimation of parameters using the search space defined as:
where upper bounds were chosen based on the order of magnitude of parameter values proposed in [2].
Estimation was based on the experimental data obtained by Neumann et al. [2] for procaspase-8 and its cleaved products p43/p41 and caspase-8.
| Time (min-1) | p43/p41 (nM) | pro-8 (nM) | casp-8 (nM) |
| 0.0 | 0.058 | 59.963 | 0.000 |
| 10.0 | 0.268 | 57.565 | 0.041 |
| 20.0 | 4.760 | 58.590 | 0.316 |
| 30.0 | 8.252 | 59.422 | 1.397 |
| 45.0 | 16.144 | 48.190 | 3.520 |
| 60.0 | 17.021 | 38.950 | 3.947 |
| 90.0 | 15.269 | 23.502 | 4.871 |
| 120.0 | 12.530 | 13.127 | 4.878 |
| 150.0 | 10.335 | 10.703 | 4.228 |
We reviewed solutions obtained by all optimization methods for 100 runs. Each run was based on the generation of 107 different guesses. The best result was obtained by the particle swarm optimization (PSO) and the cellular genetic algorithm (MOCell). Methods SRES, MOCell and PSO found similar solutions. Methods ASA and glbSolve found other values for parameters k1 and k2 showing lower efficiency.
 |
The best guesses obtained by optimization methods for 100 runs
|
Values of the objective function for 100 runs
|
[edit] Testing the computational speed of the optimization methods
- Optimization documents: data > Examples > Optimization > Data > Documents (test cases 1A, 1B, 1C, 2, 3)
- Models: data > Examples > Optimization > Data > Diagrams (diagrams 1A, 1B, 1C, 2, 3 for the corresponding test cases)
- Experimental data: data > Examples > Optimization > Data > Experiments (exp_data_1 for the test cases 1A, 1B, 1C; exp_data_2 for the test case 2; exp_data_3 for the test case 3)
We tested a computational speed of such optimization methods in BioUML as particle swarm optimization, adaptive simulated annealing, stochastic ranking evolution strategy (SRES), and cellular genetic algorithm.
For this purpose, we used biochemical models with the different number of parameters and species introduced in the following test cases.
Firstly, we derived three models of CD95-induced activation of caspase-8 from the model by Neumann et al. [2] with varying degrees of detail (the test cases 1A, 1B and 1C).
Secondly, we took the test case proposed by Mendes et al. [3] for the MAP kinase cascade model developed by Kholodenko et al. [4] (the test case 2).
Finally, we tested the model by Bagci et al. [5] representing the mitochondria-depended apoptosis (the test case 3).
As expected, the computational speed for all test cases directly depended on the number of parameters and species in the model. The greatest computational speed was shown by the method of particle swarm optimization. The other methods showed about the same speed for running in one core, wherease, for running in several cores, the computational speed of PSO, SRES and MOCell was evidently higher compared to the simulated annealing.
Test case 1A: 7 species and 5 reaction rate parameters.
Test case 1B: 13 species and 10 reaction rate parameters.
Test case 1C: 18 species and 12 reaction rate parameters.
Test case 2: 8 species and 22 reaction rate parameters.
Test case 3: 32 species and 57 reaction rate parameters.
The number of objective function evaluations per second for the different test cases in BioUML.
[edit] References
- Kutumova E., Ryabova A., Valeev T., Kolpakov F. BioUML plug-in for nonlinear parameter estimation using multiple experimental data. Virtual Biology. 2013. 1:47-58.
- Neumann L., Pforr C., Beaudouin J., Pappa A., Fricker N., Krammer P.H., Lavrik I.N., Eils R. Dynamics within the CD95 death-inducing signaling complex decide life and death of cells. Molecular Systems Biology. 2010. 6:352.
- Mendes P., Hoops S., Sahle S., Gauges R., Dada J., Kummer U. Computational modeling of biochemical networks using COPASI. Methods in Molecular Biology. 2009. 500:17–59.
- Kholodenko B.N. Negative feedback and ultrasensitivity can bring about oscillations in the mitogenactivated protein kinase cascades. European Journal of Biochemistry. 2000. 267(6):1583–1588.
- Bagci E.Z., Vodovotz Y., Billiar T.R., Ermentrout G.B., Bahar I. Bistability in apoptosis: roles of bax, bcl-2, and mitochondrial permeability transition pores. Biophysical Journal. 2006. 90(5):1546–1559.
