Fitting and comparison of calcium-calmodulin kinetic schemes to a common data set using non-linear mixed effects modelling.
Calmodulin is a calcium binding protein that is essential in calcium signalling in the brain. There are many computational models of calcium-calmodulin binding that capture various calmodulin features. However, existing models have generally been fit to different data sets, with some publications no...
Saved in:
Main Authors: | , , , , |
---|---|
Format: | Article |
Language: | English |
Published: |
Public Library of Science (PLoS)
2025-01-01
|
Series: | PLoS ONE |
Online Access: | https://doi.org/10.1371/journal.pone.0318646 |
Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
_version_ | 1823856870871269376 |
---|---|
author | Domas Linkevicius Angus Chadwick Guido C Faas Melanie I Stefan David C Sterratt |
author_facet | Domas Linkevicius Angus Chadwick Guido C Faas Melanie I Stefan David C Sterratt |
author_sort | Domas Linkevicius |
collection | DOAJ |
description | Calmodulin is a calcium binding protein that is essential in calcium signalling in the brain. There are many computational models of calcium-calmodulin binding that capture various calmodulin features. However, existing models have generally been fit to different data sets, with some publications not reporting their training and validation performance. Moreover, there is no model comparison using a common benchmark data set as is common practice in other modeling domains. Finally, some calmodulin models have been fit as a part of a larger kinetic scheme, which may have resulted in parameters being underdetermined. We address these three limitations of previous models by fitting the published calcium-calmodulin schemes to a common calcium-calmodulin data set comprising equilibrium data from Shifman et al. and dynamical data from Faas et al. Due to technical limitations, the amount of uncaged calcium in Faas et al. data could not be predicted with certainty. To find good parameter fits, despite this uncertainty, we used non-linear mixed effects modelling as implemented in the Pumas.jl package. The Akaike information criterion values for our reaction rate constants were significantly lower than for the published parameters, indicating that the published parameters are suboptimal. Moreover, there were significant differences in calmodulin activation, both between the schemes and between our reaction rate and those previously published. A kinetic scheme with independent lobes and unique, rather than identical, binding sites fit the data best. Our results support two hypotheses: (1) partially bound calmodulin is important in cellular signalling; and (2) calcium binding sites within a calmodulin lobe are kinetically distinct rather than identical. We conclude that more attention should be given to validation and comparison of models of individual molecules. |
format | Article |
id | doaj-art-f4d99a577477458e812a8b705532cc11 |
institution | Kabale University |
issn | 1932-6203 |
language | English |
publishDate | 2025-01-01 |
publisher | Public Library of Science (PLoS) |
record_format | Article |
series | PLoS ONE |
spelling | doaj-art-f4d99a577477458e812a8b705532cc112025-02-12T05:31:08ZengPublic Library of Science (PLoS)PLoS ONE1932-62032025-01-01202e031864610.1371/journal.pone.0318646Fitting and comparison of calcium-calmodulin kinetic schemes to a common data set using non-linear mixed effects modelling.Domas LinkeviciusAngus ChadwickGuido C FaasMelanie I StefanDavid C SterrattCalmodulin is a calcium binding protein that is essential in calcium signalling in the brain. There are many computational models of calcium-calmodulin binding that capture various calmodulin features. However, existing models have generally been fit to different data sets, with some publications not reporting their training and validation performance. Moreover, there is no model comparison using a common benchmark data set as is common practice in other modeling domains. Finally, some calmodulin models have been fit as a part of a larger kinetic scheme, which may have resulted in parameters being underdetermined. We address these three limitations of previous models by fitting the published calcium-calmodulin schemes to a common calcium-calmodulin data set comprising equilibrium data from Shifman et al. and dynamical data from Faas et al. Due to technical limitations, the amount of uncaged calcium in Faas et al. data could not be predicted with certainty. To find good parameter fits, despite this uncertainty, we used non-linear mixed effects modelling as implemented in the Pumas.jl package. The Akaike information criterion values for our reaction rate constants were significantly lower than for the published parameters, indicating that the published parameters are suboptimal. Moreover, there were significant differences in calmodulin activation, both between the schemes and between our reaction rate and those previously published. A kinetic scheme with independent lobes and unique, rather than identical, binding sites fit the data best. Our results support two hypotheses: (1) partially bound calmodulin is important in cellular signalling; and (2) calcium binding sites within a calmodulin lobe are kinetically distinct rather than identical. We conclude that more attention should be given to validation and comparison of models of individual molecules.https://doi.org/10.1371/journal.pone.0318646 |
spellingShingle | Domas Linkevicius Angus Chadwick Guido C Faas Melanie I Stefan David C Sterratt Fitting and comparison of calcium-calmodulin kinetic schemes to a common data set using non-linear mixed effects modelling. PLoS ONE |
title | Fitting and comparison of calcium-calmodulin kinetic schemes to a common data set using non-linear mixed effects modelling. |
title_full | Fitting and comparison of calcium-calmodulin kinetic schemes to a common data set using non-linear mixed effects modelling. |
title_fullStr | Fitting and comparison of calcium-calmodulin kinetic schemes to a common data set using non-linear mixed effects modelling. |
title_full_unstemmed | Fitting and comparison of calcium-calmodulin kinetic schemes to a common data set using non-linear mixed effects modelling. |
title_short | Fitting and comparison of calcium-calmodulin kinetic schemes to a common data set using non-linear mixed effects modelling. |
title_sort | fitting and comparison of calcium calmodulin kinetic schemes to a common data set using non linear mixed effects modelling |
url | https://doi.org/10.1371/journal.pone.0318646 |
work_keys_str_mv | AT domaslinkevicius fittingandcomparisonofcalciumcalmodulinkineticschemestoacommondatasetusingnonlinearmixedeffectsmodelling AT anguschadwick fittingandcomparisonofcalciumcalmodulinkineticschemestoacommondatasetusingnonlinearmixedeffectsmodelling AT guidocfaas fittingandcomparisonofcalciumcalmodulinkineticschemestoacommondatasetusingnonlinearmixedeffectsmodelling AT melanieistefan fittingandcomparisonofcalciumcalmodulinkineticschemestoacommondatasetusingnonlinearmixedeffectsmodelling AT davidcsterratt fittingandcomparisonofcalciumcalmodulinkineticschemestoacommondatasetusingnonlinearmixedeffectsmodelling |