Simultant: Simultaneous Curve Fitting of Functions and Differential Equations Using Analytical Gradient Calculations

Overview

Journal BMC Bioinformatics

Publisher Biomed Central

Specialty Biology

Date 2022 May 21

PMID 35597909

Authors

Julius B Kirkegaard

Affiliations

Soon will be listed here.

Abstract

Background: The initial step in comparing mathematical models to experimental data is to do a fit. This process can be complicated when either the mathematical models are not analytically solvable (e.g. because of nonlinear differential equations) or when the relation between data and models is complex (e.g. when some fitting parameters must be shared between many data sets).

Results: We introduce Simultant, a software package that allows complex fitting setups to be easily defined using a simple graphical user interface. Fitting functions can be defined directly as mathematical expressions or indirectly as the solution to specified ordinary differential equations. Analytical gradients of these functions, including the solution of differential equations, are automatically calculated to provide fast fitting even for functions with many parameters. The software enables easy definition of complex fitting setups in which parameters can be shared across both data sets and models to allow simultaneous fits to be performed.

Conclusions: Simultant exploits differentiable programming and simplifies modern fitting approaches in a unified graphical interface.

References

Meisl G, Yang X, Hellstrand E, Frohm B, Kirkegaard J, Cohen S . Differences in nucleation behavior underlie the contrasting aggregation kinetics of the Aβ40 and Aβ42 peptides. Proc Natl Acad Sci U S A. 2014; 111(26):9384-9. PMC: 4084462. DOI: 10.1073/pnas.1401564111. View

Rodrigues J, Bershtein S, Li A, Lozovsky E, Hartl D, Shakhnovich E . Biophysical principles predict fitness landscapes of drug resistance. Proc Natl Acad Sci U S A. 2016; 113(11):E1470-8. PMC: 4801265. DOI: 10.1073/pnas.1601441113. View

Cohen S, Linse S, Luheshi L, Hellstrand E, White D, Rajah L . Proliferation of amyloid-β42 aggregates occurs through a secondary nucleation mechanism. Proc Natl Acad Sci U S A. 2013; 110(24):9758-63. PMC: 3683769. DOI: 10.1073/pnas.1218402110. View

Meisl G, Kirkegaard J, Arosio P, Michaels T, Vendruscolo M, Dobson C . Molecular mechanisms of protein aggregation from global fitting of kinetic models. Nat Protoc. 2016; 11(2):252-72. DOI: 10.1038/nprot.2016.010. View

Elowitz M, Leibler S . A synthetic oscillatory network of transcriptional regulators. Nature. 2000; 403(6767):335-8. DOI: 10.1038/35002125. View

Gimenez B, Dalgaard P . Modelling and predicting the simultaneous growth of Listeria monocytogenes and spoilage micro-organisms in cold-smoked salmon. J Appl Microbiol. 2003; 96(1):96-109. DOI: 10.1046/j.1365-2672.2003.02137.x. View

Kohram M, Vashistha H, Leibler S, Xue B, Salman H . Bacterial Growth Control Mechanisms Inferred from Multivariate Statistical Analysis of Single-Cell Measurements. Curr Biol. 2020; 31(5):955-964.e4. DOI: 10.1016/j.cub.2020.11.063. View

Beechem J . Global analysis of biochemical and biophysical data. Methods Enzymol. 1992; 210:37-54. DOI: 10.1016/0076-6879(92)10004-w. View

Finke M, DEFOLIART G, Benevenga N . Use of simultaneous curve fitting and a four-parameter logistic model to evaluate the nutritional quality of protein sources at growth rates of rats from maintenance to maximum gain. J Nutr. 1987; 117(10):1681-8. DOI: 10.1093/jn/117.10.1681. View

10.

Shen-Orr S, Milo R, Mangan S, Alon U . Network motifs in the transcriptional regulation network of Escherichia coli. Nat Genet. 2002; 31(1):64-8. DOI: 10.1038/ng881. View

11.

Cleland W . Enzyme kinetics. Annu Rev Biochem. 1967; 36:77-112. DOI: 10.1146/annurev.bi.36.070167.000453. View

12.

Guo T, Abed A, Lovell N, Dokos S . Parameter fitting using multiple datasets in cardiac action potential modeling. Annu Int Conf IEEE Eng Med Biol Soc. 2012; 2011:158-61. DOI: 10.1109/IEMBS.2011.6089918. View

13.

Le Marc Y, Valik L, Medvedova A . Modelling the effect of the starter culture on the growth of Staphylococcus aureus in milk. Int J Food Microbiol. 2009; 129(3):306-11. DOI: 10.1016/j.ijfoodmicro.2008.12.015. View

14.

Rial D, Vazquez J, Murado M . Effects of three heavy metals on the bacteria growth kinetics: a bivariate model for toxicological assessment. Appl Microbiol Biotechnol. 2011; 90(3):1095-109. DOI: 10.1007/s00253-011-3138-1. View

15.

Ashyraliyev M, Siggens K, Janssens H, Blom J, Akam M, Jaeger J . Gene circuit analysis of the terminal gap gene huckebein. PLoS Comput Biol. 2009; 5(10):e1000548. PMC: 2760955. DOI: 10.1371/journal.pcbi.1000548. View

16.

Johnson K, Simpson Z, Blom T . Global kinetic explorer: a new computer program for dynamic simulation and fitting of kinetic data. Anal Biochem. 2009; 387(1):20-9. DOI: 10.1016/j.ab.2008.12.024. View

17.

HODGKIN A, HUXLEY A . A quantitative description of membrane current and its application to conduction and excitation in nerve. J Physiol. 1952; 117(4):500-44. PMC: 1392413. DOI: 10.1113/jphysiol.1952.sp004764. View

18.

Beard D . A biophysical model of the mitochondrial respiratory system and oxidative phosphorylation. PLoS Comput Biol. 2005; 1(4):e36. PMC: 1201326. DOI: 10.1371/journal.pcbi.0010036. View