A Modular Workflow for Model Building, Analysis, and Parameter Estimation in Systems Biology and Neuroscience.

João P G Santos, Kadri Pajo, Daniel Trpevski, Andrey Stepaniuk, Olivia Eriksson, Anu G Nair, Daniel Keller, Jeanette Hellgren Kotaleski, Andrei Kramer
Author Information
  1. João P G Santos: Science for Life Laboratory, School of Electrical Engineering and Computer Science, KTH Royal Institute of Technology, 10044, Stockholm, Sweden. ORCID
  2. Kadri Pajo: Department of Neuroscience, Karolinska Institute, 17165, Stockholm, Sweden. ORCID
  3. Daniel Trpevski: Science for Life Laboratory, School of Electrical Engineering and Computer Science, KTH Royal Institute of Technology, 10044, Stockholm, Sweden. ORCID
  4. Andrey Stepaniuk: Blue Brain Project, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland. ORCID
  5. Olivia Eriksson: Science for Life Laboratory, School of Electrical Engineering and Computer Science, KTH Royal Institute of Technology, 10044, Stockholm, Sweden. ORCID
  6. Anu G Nair: Department of Neuroscience, Karolinska Institute, 17165, Stockholm, Sweden. ORCID
  7. Daniel Keller: Blue Brain Project, École Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland. ORCID
  8. Jeanette Hellgren Kotaleski: Science for Life Laboratory, School of Electrical Engineering and Computer Science, KTH Royal Institute of Technology, 10044, Stockholm, Sweden. jeanette@kth.se. ORCID
  9. Andrei Kramer: Science for Life Laboratory, School of Electrical Engineering and Computer Science, KTH Royal Institute of Technology, 10044, Stockholm, Sweden. andreikr@kth.se. ORCID

Abstract

Neuroscience incorporates knowledge from a range of scales, from single molecules to brain wide neural networks. Modeling is a valuable tool in understanding processes at a single scale or the interactions between two adjacent scales and researchers use a variety of different software tools in the model building and analysis process. Here we focus on the scale of biochemical pathways, which is one of the main objects of study in systems biology. While systems biology is among the more standardized fields, conversion between different model formats and interoperability between various tools is still somewhat problematic. To offer our take on tackling these shortcomings and by keeping in mind the FAIR (findability, accessibility, interoperability, reusability) data principles, we have developed a workflow for building and analyzing biochemical pathway models, using pre-existing tools that could be utilized for the storage and refinement of models in all phases of development. We have chosen the SBtab format which allows the storage of biochemical models and associated data in a single file and provides a human readable set of syntax rules. Next, we implemented custom-made MATLAB® scripts to perform parameter estimation and global sensitivity analysis used in model refinement. Additionally, we have developed a web-based application for biochemical models that allows simulations with either a network free solver or stochastic solvers and incorporating geometry. Finally, we illustrate convertibility and use of a biochemical model in a biophysically detailed single neuron model by running multiscale simulations in NEURON. Using this workflow, we can simulate the same model in three different simulators, with a smooth conversion between the different model formats, enhancing the characterization of different aspects of the model.

Keywords

References

  1. Nat Methods. 2011 Feb;8(2):177-83 [PMID: 21186362]
  2. Bioinformatics. 2006 Dec 15;22(24):3067-74 [PMID: 17032683]
  3. Bioinformatics. 2018 Oct 15;34(20):3591-3593 [PMID: 29762723]
  4. Bioinformatics. 2018 Nov 1;34(21):3702-3710 [PMID: 29790940]
  5. BMC Syst Biol. 2012 May 10;6:36 [PMID: 22574658]
  6. Bioinformatics. 2008 Sep 15;24(18):2037-43 [PMID: 18614583]
  7. PLoS Comput Biol. 2010 Jun 17;6(6):e1000815 [PMID: 20585541]
  8. IET Syst Biol. 2011 Nov;5(6):336-6 [PMID: 22129029]
  9. BMC Bioinformatics. 2016 Apr 05;17:154 [PMID: 27044654]
  10. Bioinformatics. 2006 Feb 15;22(4):514-5 [PMID: 16317076]
  11. Nucleic Acids Res. 2016 Jul 8;44(W1):W3-W10 [PMID: 27137889]
  12. J Math Neurosci. 2020 Apr 6;10(1):5 [PMID: 32253526]
  13. Science. 1999 Jan 15;283(5400):381-7 [PMID: 9888852]
  14. Wiley Interdiscip Rev Syst Biol Med. 2012 Mar-Apr;4(2):129-40 [PMID: 22139996]
  15. Neuroinformatics. 2007 Summer;5(2):127-38 [PMID: 17873374]
  16. Bioinformatics. 2019 Jan 15;35(2):284-292 [PMID: 30010712]
  17. Bioinformatics. 2016 Aug 15;32(16):2559-61 [PMID: 27153616]
  18. Front Neuroinform. 2013 Nov 15;7:28 [PMID: 24298253]
  19. Science. 2014 Sep 26;345(6204):1616-20 [PMID: 25258080]
  20. Front Neuroinform. 2018 Jun 26;12:38 [PMID: 29997492]
  21. Nat Rev Neurosci. 2010 Apr;11(4):239-51 [PMID: 20300102]
  22. Bioinformatics. 2009 Aug 1;25(15):1923-9 [PMID: 19505944]
  23. Bioinformatics. 2009 Jun 1;25(11):1453-4 [PMID: 19304877]
  24. J Neurophysiol. 2020 Mar 1;123(3):1042-1051 [PMID: 31851573]
  25. Front Neuroinform. 2008 Dec 19;2:6 [PMID: 19129924]
  26. Bioinformatics. 2003 Mar 1;19(4):524-31 [PMID: 12611808]
  27. Neural Comput. 1997 Aug 15;9(6):1179-209 [PMID: 9248061]
  28. J Comput Neurosci. 2009 Dec;27(3):471-91 [PMID: 19533315]
  29. Curr Opin Genet Dev. 2004 Aug;14(4):375-81 [PMID: 15261653]
  30. Bioinformatics. 2019 Sep 1;35(17):3073-3082 [PMID: 30624608]

MeSH Term

Humans
Models, Biological
Neurons
Neurosciences
Software
Systems Biology
Workflow

Word Cloud

Created with Highcharts 10.0.0modeldifferentbiochemicalsinglemodelstoolsanalysisbiologyNeurosciencescalesscaleusebuildingsystemsconversionformatsinteroperabilitydatadevelopedworkflowstoragerefinementSBtaballowsestimationsensitivitysimulationsParameterSystemsincorporatesknowledgerangemoleculesbrainwideneuralnetworksModelingvaluabletoolunderstandingprocessesinteractionstwoadjacentresearchersvarietysoftwareprocessfocuspathwaysonemainobjectsstudyamongstandardizedfieldsvariousstillsomewhatproblematicoffertaketacklingshortcomingskeepingmindFAIRfindabilityaccessibilityreusabilityprinciplesanalyzingpathwayusingpre-existingutilizedphasesdevelopmentchosenformatassociatedfileprovideshumanreadablesetsyntaxrulesNextimplementedcustom-madeMATLAB®scriptsperformparameterglobalusedAdditionallyweb-basedapplicationeithernetworkfreesolverstochasticsolversincorporatinggeometryFinallyillustrateconvertibilitybiophysicallydetailedneuronrunningmultiscaleNEURONUsingcansimulatethreesimulatorssmoothenhancingcharacterizationaspectsModularWorkflowModelBuildingAnalysisEstimationBiologyGlobalInteroperabilityMultiscalemodeling

Similar Articles

Cited By (1)