Cargando…
Regression Analysis for Constraining Free Parameters in Electrophysiological Models of Cardiac Cells
A major challenge in computational biology is constraining free parameters in mathematical models. Adjusting a parameter to make a given model output more realistic sometimes has unexpected and undesirable effects on other model behaviors. Here, we extend a regression-based method for parameter sens...
Autores principales: | , |
---|---|
Formato: | Texto |
Lenguaje: | English |
Publicado: |
Public Library of Science
2010
|
Materias: | |
Acceso en línea: | https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2932676/ https://www.ncbi.nlm.nih.gov/pubmed/20824123 http://dx.doi.org/10.1371/journal.pcbi.1000914 |
_version_ | 1782186087260618752 |
---|---|
author | Sarkar, Amrita X. Sobie, Eric A. |
author_facet | Sarkar, Amrita X. Sobie, Eric A. |
author_sort | Sarkar, Amrita X. |
collection | PubMed |
description | A major challenge in computational biology is constraining free parameters in mathematical models. Adjusting a parameter to make a given model output more realistic sometimes has unexpected and undesirable effects on other model behaviors. Here, we extend a regression-based method for parameter sensitivity analysis and show that a straightforward procedure can uniquely define most ionic conductances in a well-known model of the human ventricular myocyte. The model's parameter sensitivity was analyzed by randomizing ionic conductances, running repeated simulations to measure physiological outputs, then collecting the randomized parameters and simulation results as “input” and “output” matrices, respectively. Multivariable regression derived a matrix whose elements indicate how changes in conductances influence model outputs. We show here that if the number of linearly-independent outputs equals the number of inputs, the regression matrix can be inverted. This is significant, because it implies that the inverted matrix can specify the ionic conductances that are required to generate a particular combination of model outputs. Applying this idea to the myocyte model tested, we found that most ionic conductances could be specified with precision (R(2) > 0.77 for 12 out of 16 parameters). We also applied this method to a test case of changes in electrophysiology caused by heart failure and found that changes in most parameters could be well predicted. We complemented our findings using a Bayesian approach to demonstrate that model parameters cannot be specified using limited outputs, but they can be successfully constrained if multiple outputs are considered. Our results place on a solid mathematical footing the intuition-based procedure simultaneously matching a model's output to several data sets. More generally, this method shows promise as a tool to define model parameters, in electrophysiology and in other biological fields. |
format | Text |
id | pubmed-2932676 |
institution | National Center for Biotechnology Information |
language | English |
publishDate | 2010 |
publisher | Public Library of Science |
record_format | MEDLINE/PubMed |
spelling | pubmed-29326762010-09-07 Regression Analysis for Constraining Free Parameters in Electrophysiological Models of Cardiac Cells Sarkar, Amrita X. Sobie, Eric A. PLoS Comput Biol Research Article A major challenge in computational biology is constraining free parameters in mathematical models. Adjusting a parameter to make a given model output more realistic sometimes has unexpected and undesirable effects on other model behaviors. Here, we extend a regression-based method for parameter sensitivity analysis and show that a straightforward procedure can uniquely define most ionic conductances in a well-known model of the human ventricular myocyte. The model's parameter sensitivity was analyzed by randomizing ionic conductances, running repeated simulations to measure physiological outputs, then collecting the randomized parameters and simulation results as “input” and “output” matrices, respectively. Multivariable regression derived a matrix whose elements indicate how changes in conductances influence model outputs. We show here that if the number of linearly-independent outputs equals the number of inputs, the regression matrix can be inverted. This is significant, because it implies that the inverted matrix can specify the ionic conductances that are required to generate a particular combination of model outputs. Applying this idea to the myocyte model tested, we found that most ionic conductances could be specified with precision (R(2) > 0.77 for 12 out of 16 parameters). We also applied this method to a test case of changes in electrophysiology caused by heart failure and found that changes in most parameters could be well predicted. We complemented our findings using a Bayesian approach to demonstrate that model parameters cannot be specified using limited outputs, but they can be successfully constrained if multiple outputs are considered. Our results place on a solid mathematical footing the intuition-based procedure simultaneously matching a model's output to several data sets. More generally, this method shows promise as a tool to define model parameters, in electrophysiology and in other biological fields. Public Library of Science 2010-09-02 /pmc/articles/PMC2932676/ /pubmed/20824123 http://dx.doi.org/10.1371/journal.pcbi.1000914 Text en Sarkar, Sobie. http://creativecommons.org/licenses/by/4.0/ This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are properly credited. |
spellingShingle | Research Article Sarkar, Amrita X. Sobie, Eric A. Regression Analysis for Constraining Free Parameters in Electrophysiological Models of Cardiac Cells |
title | Regression Analysis for Constraining Free Parameters in Electrophysiological Models of Cardiac Cells |
title_full | Regression Analysis for Constraining Free Parameters in Electrophysiological Models of Cardiac Cells |
title_fullStr | Regression Analysis for Constraining Free Parameters in Electrophysiological Models of Cardiac Cells |
title_full_unstemmed | Regression Analysis for Constraining Free Parameters in Electrophysiological Models of Cardiac Cells |
title_short | Regression Analysis for Constraining Free Parameters in Electrophysiological Models of Cardiac Cells |
title_sort | regression analysis for constraining free parameters in electrophysiological models of cardiac cells |
topic | Research Article |
url | https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2932676/ https://www.ncbi.nlm.nih.gov/pubmed/20824123 http://dx.doi.org/10.1371/journal.pcbi.1000914 |
work_keys_str_mv | AT sarkaramritax regressionanalysisforconstrainingfreeparametersinelectrophysiologicalmodelsofcardiaccells AT sobieerica regressionanalysisforconstrainingfreeparametersinelectrophysiologicalmodelsofcardiaccells |