Publications
Journal Article
Arrhythmogenic influence of mutations in a myocyte‑based computational model of the pulmonary vein sleeve
Nature Scientific Reports 12 (2022): 7040.Status: Published
Arrhythmogenic influence of mutations in a myocyte‑based computational model of the pulmonary vein sleeve
In the heart, electrophysiological dysregulation arises from defects at many biological levels (from point mutations in ion channel proteins to gross structural abnormalities). These defects disrupt the normal pattern of electrical activation, producing ectopic activity and reentrant arrhythmia. To interrogate mechanisms that link these primary biological defects to macroscopic electrophysiologic dysregulation most prior computational studies have utilized either (i) detailed models of myocyte ion channel dynamics at limited spatial scales, or (ii) homogenized models of action potential conduction that reproduce arrhythmic activity at tissue and organ levels. Here we apply our recent model (EMI), which integrates electrical activation and propagation across these scales, to study human atrial arrhythmias originating in the pulmonary vein (PV) sleeves. These small structures initiate most supraventricular arrhythmias and include pronounced myocyte‑to‑myocyte heterogeneities in ion channel expression and intercellular coupling. To test EMI’s cell‑based architecture in this physiological context we asked whether ion channel mutations known to underlie atrial fibrillation are capable of initiating arrhythmogenic behavior via increased excitability or reentry in a schematic PV sleeve geometry. Our results illustrate that EMI’s improved spatial resolution can directly interrogate how electrophysiological changes at the individual myocyte level manifest in tissue and as arrhythmia in the PV sleeve.
Afilliation | Scientific Computing |
Project(s) | Department of Computational Physiology |
Publication Type | Journal Article |
Year of Publication | 2022 |
Journal | Nature Scientific Reports |
Volume | 12 |
Pagination | 7040 |
Publisher | Springer Nature |
URL | https://rdcu.be/cMntU |
DOI | 10.1038/s41598-022-11110-1 |
Deriving the Bidomain Model of Cardiac Electrophysiology From a Cell-Based Model; Properties and Comparisons
Frontiers in Physiology 12 (2022): 811029.Status: Published
Deriving the Bidomain Model of Cardiac Electrophysiology From a Cell-Based Model; Properties and Comparisons
The bidomain model is considered to be the gold standard for numerical simulation of the electrophysiology of cardiac tissue. The model provides important insights into the conduction properties of the electrochemical wave traversing the cardiac muscle in every heartbeat. However, in normal resolution, the model represents the average over a large number of cardiomyocytes, and more accurate models based on representations of all individual cells have therefore been introduced in order to gain insight into the conduction properties close to the myocytes. The more accurate model considered here is referred to as the EMI model since both the extracellular space (E), the cell membrane (M) and the intracellular space (I) are explicitly represented in the model. Here, we show that the bidomain model can be derived from the cell-based EMI model and we thus reveal the close relation between the two models, and obtain an indication of the error introduced in the approximation. Also, we present numerical simulations comparing the results of the two models and thereby highlight both similarities and differences between the models. We observe that the deviations between the solutions of the models become larger for larger cell sizes. Furthermore, we observe that the bidomain model provides solutions that are very similar to the EMI model when conductive properties of the tissue are in the normal range, but large deviations are present when the resistance between cardiomyocytes is increased.
Afilliation | Scientific Computing |
Project(s) | Department of Computational Physiology |
Publication Type | Journal Article |
Year of Publication | 2022 |
Journal | Frontiers in Physiology |
Volume | 12 |
Pagination | 811029 |
Publisher | Frontiers |
URL | https://www.frontiersin.org/articles/10.3389/fphys.2021.811029/full |
DOI | 10.3389/fphys.2021.811029 |
Journal Article
A computational method for identifying an optimal combination of existing drugs to repair the action potentials of SQT1 ventricular myocytes
PLoS Computational Biology 17 (2021): e1009233.Status: Published
A computational method for identifying an optimal combination of existing drugs to repair the action potentials of SQT1 ventricular myocytes
Mutations are known to cause perturbations in essential functional features of integral membrane proteins, including ion channels. Even restricted or point mutations can result in substantially changed properties of ion currents. The additive effect of these alterations for a specific ion channel can result in significantly changed properties of the action potential (AP). Both AP shortening and AP prolongation can result from known mutations, and the consequences can be life-threatening. Here, we present a computational method for identifying new drugs utilizing combinations of existing drugs. Based on the knowledge of theoretical effects of existing drugs on individual ion currents, our aim is to compute optimal combinations that can ‘repair’ the mutant AP waveforms so that the baseline AP-properties are restored. More specifically, we compute optimal, combined, drug concentrations such that the waveforms of the transmembrane potential and the cytosolic calcium concentration of the mutant cardiomyocytes (CMs) becomes as similar as possible to their wild type coun- terparts after the drug has been applied. In order to demonstrate the utility of this method, we address the question of computing an optimal drug for the short QT syndrome type 1 (SQT1). For the SQT1 mutation N588K, there are available data sets that describe the effect of various drugs on the mutated K+ channel. These published findings are the basis for our computational analysis which can identify optimal compounds in the sense that the AP of the mutant CMs resembles essential biomarkers of the wild type CMs. Using recently developed insights regarding electrophysiological properties among myocytes from different species, we compute optimal drug combinations for hiPSC-CMs, rabbit ventricular CMs and adult human ventricular CMs with the SQT1 mutation. Since the ‘composition’ of ion channels that form the AP is different for the three types of myocytes under consideration, so is the composition of the optimal drug.
Afilliation | Scientific Computing |
Project(s) | Department of Computational Physiology |
Publication Type | Journal Article |
Year of Publication | 2021 |
Journal | PLoS Computational Biology |
Volume | 17 |
Number | 8 |
Pagination | e1009233 |
Publisher | Public Library of Science |
DOI | 10.1371/journal.pcbi.1009233 |
Computational prediction of drug response in short QT syndrome type 1 based on measurements of compound effect in stem cell-derived cardiomyocytes
PLoS Computational Biology 17, no. 2 (2021): e1008089.Status: Published
Computational prediction of drug response in short QT syndrome type 1 based on measurements of compound effect in stem cell-derived cardiomyocytes
Short QT (SQT) syndrome is a genetic cardiac disorder characterized by an abbreviated QT interval of the patient’s electrocardiogram. The syndrome is associated with increased risk of arrhythmia and sudden cardiac death and can arise from a number of ion channel mutations. Cardiomyocytes derived from induced pluripotent stem cells generated from SQT patients (SQT hiPSC-CMs) provide promising platforms for testing pharmacological treatments directly in human cardiac cells exhibiting mutations specific for the syndrome. How- ever, a difficulty is posed by the relative immaturity of hiPSC-CMs, with the possibility that drug effects observed in SQT hiPSC-CMs could be very different from the corresponding drug effect in vivo. In this paper, we apply a multistep computational procedure for translating measured drug effects from these cells to human QT response. This process first detects drug effects on individual ion channels based on measurements of SQT hiPSC-CMs and then uses these results to estimate the drug effects on ventricular action potentials and QT intervals of adult SQT patients. We find that the procedure is able to identify IC50 values in line with measured values for the four drugs quinidine, ivabradine, ajmaline and mexiletine. In addition, the predicted effect of quinidine on the adult QT interval is in good agreement with measured effects of quinidine for adult patients. Consequently, the computational procedure appears to be a useful tool for helping predicting adult drug responses from pure in vitro measurements of patient derived cell lines.
Afilliation | Scientific Computing |
Project(s) | Department of Computational Physiology |
Publication Type | Journal Article |
Year of Publication | 2021 |
Journal | PLoS Computational Biology |
Volume | 17 |
Issue | 2 |
Pagination | e1008089 |
Publisher | Public Library of Science |
URL | https://journals.plos.org/ploscompbiol/article/file?id=10.1371/journal.p... |
DOI | 10.1371/journal.pcbi.1008089 |
Efficient numerical solution of the EMI model representing the extracellular space (E), cell membrane (M) and intracellular space (I) of a collection of cardiac cells
Frontiers in Physics 8 (2021): 579461.Status: Published
Efficient numerical solution of the EMI model representing the extracellular space (E), cell membrane (M) and intracellular space (I) of a collection of cardiac cells
The EMI model represents excitable cells in a more accurate manner than traditional homogenized models at the price of increased computational complexity. The increased complexity of solving the EMI model stems from a significant increase in the number of computational nodes and from the form of the linear systems that need to be solved. Here, we will show that the latter problem can be solved by careful use of operator splitting of the spatially coupled equations. By using this method, the linear systems can be broken into sub-problems that are of the classical type of linear, elliptic boundary value problems. Therefore, the vast collection of methods for solving linear, elliptic partial differential equations can be used. We demonstrate that this enables us to solve the systems using shared-memory parallel computers. The computing time scales perfectly with the number of physical cells. For a collection of 512×256 cells, we manage to solve linear systems with about 2.5×10^8 unknows. Since the computational effort scales linearly with the number of physical cells, we believe that larger computers can be used to simulate millions of excitable cells and thus allow careful analysis of physiological systems of great importance.
Afilliation | Scientific Computing |
Project(s) | Department of Computational Physiology, Department of High Performance Computing |
Publication Type | Journal Article |
Year of Publication | 2021 |
Journal | Frontiers in Physics |
Volume | 8 |
Pagination | 579461 |
Publisher | Frontiers |
URL | https://www.frontiersin.org/articles/10.3389/fphy.2020.579461/full |
DOI | 10.3389/fphy.2020.579461 |
From Millimeters to Micrometers; Re-introducing Myocytes in Models of Cardiac Electrophysiology
Frontiers in Physiology 12 (2021): 763584.Status: Published
From Millimeters to Micrometers; Re-introducing Myocytes in Models of Cardiac Electrophysiology
Computational modeling has contributed significantly to present understanding of cardiac electrophysiology including cardiac conduction, excitation-contraction coupling, and the effects and side-effects of drugs. However, the accuracy of in silico analysis of electrochemical wave dynamics in cardiac tissue is limited by the homogenization procedure (spatial averaging) intrinsic to standard continuum models of conduction. Averaged models cannot resolve the intricate dynamics in the vicinity of individual cardiomyocytes simply because the myocytes are not present in these models. Here we demonstrate how recently developed mathematical models based on representing every myocyte can significantly increase the accuracy, and thus the utility of modeling electrophysiological function and dysfunction in collections of coupled cardiomyocytes. The present gold standard of numerical simulation for cardiac electrophysiology is based on the bidomain model. In the bidomain model, the extracellular (E) space, the cell membrane (M) and the intracellular (I) space are all assumed to be present everywhere in the tissue. Consequently, it is impossible to study biophysical processes taking place close to individual myocytes. The bidomain model represents the tissue by averaging over several hundred myocytes and this inherently limits the accuracy of the model. In our alternative approach both E, M, and I are represented in the model which is therefore referred to as the EMI model. The EMI model approach allows for detailed analysis of the biophysical processes going on in functionally important spaces very close to individual myocytes, although at the cost of significantly increased CPU-requirements.
Afilliation | Scientific Computing |
Project(s) | Department of Computational Physiology |
Publication Type | Journal Article |
Year of Publication | 2021 |
Journal | Frontiers in Physiology |
Volume | 12 |
Pagination | 763584 |
Publisher | Frontiers |
URL | https://www.frontiersin.org/articles/10.3389/fphys.2021.763584/full |
DOI | 10.3389/fphys.2021.763584 |
Identifying drug response by combining measurements of the membrane potential, the cytosolic calcium concentration, and the extracellular potential in microphysiological systems
Frontiers in Pharmacology 11 (2021): 569489.Status: Published
Identifying drug response by combining measurements of the membrane potential, the cytosolic calcium concentration, and the extracellular potential in microphysiological systems
Cardiomyocytes derived from human induced pluripotent stem cells (hiPSC-CMs) offer a new means to study and understand the human cardiac action potential, and can give key insight into how compounds may interact with important molecular pathways to destabilize the electrical function of the heart. Important features of the action potential can be readily measured using standard experimental techniques, such as the use of voltage sensitive dyes and fluorescent genetic reporters to estimate transmembrane potentials and cytosolic calcium concentrations. Using previously introduced computational procedures, such measurements can be used to estimate the current density of major ion channels present in hiPSC-CMs, and how compounds may alter their behavior. However, due to the limitations of optical recordings, resolving the sodium current remains difficult from these data. Here we show that if these optical measurements are complemented with observations of the extracellular potential using multi electrode arrays (MEAs), we can accurately estimate the current density of the sodium channels. This inversion of the sodium current relies on observation of the conduction velocity which turns out to be straightforwardly computed using measurements of extracellular waves across the electrodes. The combined data including the membrane potential, the cytosolic calcium concentration and the extracellular potential further opens up for the possibility of accurately estimating the effect of novel drugs applied to hiPSC-CMs.
Afilliation | Scientific Computing |
Project(s) | Department of Computational Physiology |
Publication Type | Journal Article |
Year of Publication | 2021 |
Journal | Frontiers in Pharmacology |
Volume | 11 |
Pagination | 569489 |
Publisher | Frontiers |
URL | https://www.frontiersin.org/articles/10.3389/fphar.2020.569489/full |
DOI | 10.3389/fphar.2020.569489 |
Book Chapter
Derivation of a Cell-Based Mathematical Model of Excitable Cells
In Modeling Excitable Tissue: The EMI Framework, 1-13. Vol. 7. Cham: Springer International Publishing, 2021.Status: Published
Derivation of a Cell-Based Mathematical Model of Excitable Cells
Excitable cells are of vital importance in biology, and mathematical models have contributed significantly to understand their basic mechanisms. However, classical models of excitable cells are based on severe assumptions that may limit the accuracy of the simulation results. Here, we derive a more detailed approach to modeling that has recently been applied to study the electrical properties of both neurons and cardiomyocytes. The model is derived from first principles and opens up possibilities for studying detailed properties of excitable cells. We refer to the model as the EMI model because both the extracellular space (E), the cell membrane (M) and the intracellular space (I) are explicitly represented in the model, in contrast to classical spatial models of excitable cells. Later chapters of the present text will focus on numerical methods and software for solving the model. Also, in the next chapter, the model will be extended to account for ionic concentrations in the intracellular and extracellular spaces.
Afilliation | Scientific Computing |
Project(s) | Department of Computational Physiology |
Publication Type | Book Chapter |
Year of Publication | 2021 |
Book Title | Modeling Excitable Tissue: The EMI Framework |
Volume | 7 |
Chapter | 1 |
Pagination | 1-13 |
Publisher | Springer International Publishing |
Place Published | Cham |
ISBN Number | 978-3-030-61157-6 |
ISBN | 2512-1677 |
URL | https://link.springer.com/content/pdf/10.1007%2F978-3-030-61157-6_1.pdf |
DOI | 10.1007/978-3-030-61157-6_1 |
Operator Splitting and Finite Difference Schemes for Solving the EMI Model
In Modeling Excitable Tissue: The EMI Framework, 44-55. Vol. 7. Cham: Springer International Publishing, 2021.Status: Published
Operator Splitting and Finite Difference Schemes for Solving the EMI Model
We want to be able to perform accurate simulations of a large number of cardiac cells based on mathematical models where each individual cell is represented in the model. This implies that the computational mesh has to have a typical resolution of a few µm leading to huge computational challenges. In this paper we use a certain operator splitting of the coupled equations and showthat this leads to systems that can be solved in parallel. This opens up for the possibility of simulating large numbers of coupled cardiac cells.
Afilliation | Scientific Computing |
Project(s) | Department of Computational Physiology, Department of High Performance Computing |
Publication Type | Book Chapter |
Year of Publication | 2021 |
Book Title | Modeling Excitable Tissue: The EMI Framework |
Volume | 7 |
Chapter | 4 |
Pagination | 44 - 55 |
Publisher | Springer International Publishing |
Place Published | Cham |
ISBN Number | 978-3-030-61156-9 |
ISBN | 2512-1677 |
URL | http://link.springer.com/content/pdf/10.1007/978-3-030-61157-6_4 |
DOI | 10.1007/978-3-030-61157-6_4 |
Journal Article
Computational translation of drug effects from animal experiments to human ventricular myocytes
Nature Scientific Reports (2020): 10537.Status: Published
Computational translation of drug effects from animal experiments to human ventricular myocytes
Using animal cells and tissues as precise measuring devices for developing new drugs presents a long- standing challenge for the pharmaceutical industry. Despite the very significant resources that continue to be dedicated to animal testing of new compounds, only qualitative results can be obtained. This often results in both false positives and false negatives. Here, we show how the effect of drugs applied to animal ventricular myocytes can be translated, quantitatively, to estimate a number of different effects of the same drug on human cardiomyocytes. We illustrate and validate our methodology by translating, from animal to human, the effect of dofetilide applied to dog cardiomyocytes, the effect of E-4031 applied to zebrafish cardiomyocytes, and, finally, the effect of sotalol applied to rabbit cardiomyocytes. In all cases, the accuracy of our quantitative estimates are demonstrated. Our computations reveal that, in principle, electrophysiological data from testing using animal ventricular myocytes, can give precise, quantitative estimates of the effect of new compounds on human cardiomyocytes.
Afilliation | Scientific Computing |
Project(s) | Department of Computational Physiology |
Publication Type | Journal Article |
Year of Publication | 2020 |
Journal | Nature Scientific Reports |
Number | 10 |
Pagination | 10537 |
Publisher | Nature Publishing Group |
URL | https://rdcu.be/b5iBl |
DOI | 10.1038/s41598-018-35858-7 |
Improved computational identification of drug response using optical measurements of human stem cell derived cardiomyocytes in microphysiological systems
Frontiers in Pharmacology 10 (2020).Status: Published
Improved computational identification of drug response using optical measurements of human stem cell derived cardiomyocytes in microphysiological systems
Cardiomyocytes derived from human induced pluripotent stem cells (hiPSC-CMs) hold great potential for drug screening applications. However, their usefulness is limited by the relative immaturity of the cells’ electrophysiological properties as compared to native cardiomyocytes in the adult human heart. In this work, we extend and improve on methodology to address this limitation, building on previously introduced computational procedures which predict drug effects for adult cells based on changes in optical measurements of action potentials and calcium transients made in stem cell derived cardiac microtissues. This methodology quantifies ion channel changes through the inversion of data into a mathematical model, and maps this response to an adult phenotype through the assumption of functional invariance of fundamental intracellular and membrane channels during maturation. Here we utilize an updated action potential model to represent both hiPSC-CMs and adult cardiomyocytes, apply an IC50-based model of dose-dependent drug effects, and introduce a continuation-based optimization algorithm for analysis of dose escalation measurements using five drugs with known effects. The improved methodology can identify drug induced changes more efficiently, and quantitate important metrics such as IC50 in line with published values. Consequently, the updated methodology is a step towards employing computational procedures to elucidate drug effects in adult cardiomyocytes for new drugs using stem cell-derived experimental tissues.
Afilliation | Scientific Computing |
Project(s) | Department of Computational Physiology, IdentiPhy |
Publication Type | Journal Article |
Year of Publication | 2020 |
Journal | Frontiers in Pharmacology |
Volume | 10 |
Number | 1648 |
Date Published | 02/2020 |
Publisher | Frontiers |
DOI | 10.3389/fphar.2019.01648 |
PhD Thesis
Cell-Based Mathematical Models of Small Collections of Excitable Cells
In The University of Oslo. Vol. PhD. Department of Informatics, University of Oslo, 2019.Status: Published
Cell-Based Mathematical Models of Small Collections of Excitable Cells
Afilliation | Scientific Computing |
Project(s) | Department of Computational Physiology |
Publication Type | PhD Thesis |
Year of Publication | 2019 |
Degree awarding institution | The University of Oslo |
Degree | PhD |
Publisher | Department of Informatics, University of Oslo |
Journal Article
Detecting undetectables: Can conductances of action potential models be changed without appreciable change in the transmembrane potential?
Chaos 29 (2019).Status: Published
Detecting undetectables: Can conductances of action potential models be changed without appreciable change in the transmembrane potential?
Mathematical models describing the dynamics of the cardiac action potential are of great value for understanding how changes to the system can disrupt the normal electrical activity of cells and tissue in the heart. However, to represent specific data, these models must be parameterized, and adjustment of the maximum conductances of the individual contributing ionic currents is a commonly used method. Here, we present a method for investigating the uniqueness of such resulting parameterizations. Our key question is: Can the maximum conductances of a model be changed without giving any appreciable changes in the action potential? If so, the model parameters are not unique and this poses a major problem in using the models to identify changes in parameters from data, for instance, to evaluate potential drug effects. We propose a method for evaluating this uniqueness, founded on the singular value decomposition of a matrix consisting of the individual ionic currents. Small singular values of this matrix signify lack of parameter uniqueness and we show that the conclusion from linear analysis of the matrix carries over to provide insight into the uniqueness of the parameters in the nonlinear case. Using numerical experiments, we quantify the identifiability of the maximum conductances of well-known models of the cardiac action potential. Furthermore, we show how the identifiability depends on the time step used in the observation of the currents, how the application of drugs may change identifiability, and, finally, how the stimulation protocol can be used to improve the identifiability of a model.
Afilliation | Scientific Computing |
Project(s) | Department of Computational Physiology |
Publication Type | Journal Article |
Year of Publication | 2019 |
Journal | Chaos |
Volume | 29 |
Number | 073102 |
Date Published | 07/2019 |
Publisher | AIP |
URL | https://aip.scitation.org/doi/10.1063/1.5087629 |
DOI | 10.1063/1.5087629 |
How does the presence of neural probes affect extracellular potentials?
Journal of Neural Engineering 16 (2019): 026030.Status: Published
How does the presence of neural probes affect extracellular potentials?
{Objective. Mechanistic modeling of neurons is an essential component of computational neuroscience that enables scientists to simulate, explain, and explore neural activity. The conventional approach to simulation of extracellular neural recordings first computes transmembrane currents using the cable equation and then sums their contribution to model the extracellular potential. This two-step approach relies on the assumption that the extracellular space is an infinite and homogeneous conductive medium, while measurements are performed using neural probes. The main purpose of this paper is to assess to what extent the presence of the neural probes of varying shape and size impacts the extracellular field and how to correct for them. Approach. We apply a detailed modeling framework allowing explicit representation of the neuron and the probe to study the effect of the probes and thereby estimate the effect of ignoring it. We use meshes with simplified neurons and different types of probe and compare the extracellular action potentials with and without the probe in the extracellular space. We then compare various solutions to account for the probes’ presence and introduce an efficient probe correction method to include the probe effect in modeling of extracellular potentials. Main results. Our computations show that microwires hardly influence the extracellular electric field and their effect can therefore be ignored. In contrast, multi-electrode arrays (MEAs) significantly affect the extracellular field by magnifying the recorded potential. While MEAs behave similarly to infinite insulated planes, we find that their effect strongly depends on the neuron-probe alignment and probe orientation. Significance. Ignoring the probe effect might be deleterious in some applications, such as neural localization and parameterization of neural models from extracellular recordings. Moreover, the presence of the probe can improve the interpretation of extracellular recordings, by providing a more accurate estimation of the extracellular potential generated by neuronal models.
Afilliation | Scientific Computing |
Project(s) | No Simula project |
Publication Type | Journal Article |
Year of Publication | 2019 |
Journal | Journal of Neural Engineering |
Volume | 16 |
Number | 2 |
Pagination | 026030 |
Date Published | feb |
Publisher | IOP} Publishing |
URL | https://doi.org/10.1088%2F1741-2552%2Fab03a1 |
DOI | 10.1088/1741-2552/ab03a1 |
Properties of cardiac conduction in a cell-based computational model
PLoS Computational Biology 15, no. 5 (2019).Status: Published
Properties of cardiac conduction in a cell-based computational model
The conduction of electrical signals through cardiac tissue is essential for maintaining the function of the heart, and conduction abnormalities are known to potentially lead to life-threatening arrhythmias. The properties of cardiac conduction have therefore been the topic of intense study for decades, but a number of questions related to the mechanisms of conduction still remain unresolved. In this paper, we demonstrate how the so-called EMI model may be used to study some of these open questions. In the EMI model, the extracellular space, the cell membrane, the intracellular space and the cell connections are all represented as separate parts of the computational domain, and the model therefore allows for study of local properties that are hard to represent in the classical homogenized bidomain or monodomain models commonly used to study cardiac conduction. We conclude that a non-uniform sodium channel distribution increases the conduction velocity and decreases the time delays over gap junctions of reduced coupling in the EMI model simulations. We also present a theoretical optimal cell length with respect to conduction velocity and consider the possibility of ephaptic coupling (i.e. cell-to-cell coupling through the extracellular potential) acting as an alternative or supporting mechanism to gap junction coupling. We conclude that for a non-uniform distribution of sodium channels and a sufficiently small intercellular distance, ephaptic coupling can influence the dynamics of the sodium channels and potentially provide cell-to-cell coupling when the gap junction connection is absent.
Afilliation | Scientific Computing |
Project(s) | Department of Computational Physiology |
Publication Type | Journal Article |
Year of Publication | 2019 |
Journal | PLoS Computational Biology |
Volume | 15 |
Issue | 5 |
Number | e1007042 |
Date Published | 05/2019 |
Publisher | Public Library of Science |
URL | https://journals.plos.org/ploscompbiol/article?id=10.1371/journal.pcbi.1... |
DOI | 10.1371/journal.pcbi.1007042 |
Proceedings, refereed
Investigating the Pro- and Anti-Arrhythmic Properties of Human Induced Pluripotent Stem Cell-Derived Cardiomyocytes in Post-Infarction Patient Hearts: A Modeling Study
In 6th International Conference on Computational and Mathematical Biomedical Engineering. Vol. 1. Zeta Computational Resources Ltd, 2019.Status: Published
Investigating the Pro- and Anti-Arrhythmic Properties of Human Induced Pluripotent Stem Cell-Derived Cardiomyocytes in Post-Infarction Patient Hearts: A Modeling Study
Afilliation | Scientific Computing |
Project(s) | Department of Computational Physiology, MI-RISK: Risk factors for sudden cardiac death during acute myocardial infarction |
Publication Type | Proceedings, refereed |
Year of Publication | 2019 |
Conference Name | 6th International Conference on Computational and Mathematical Biomedical Engineering |
Volume | 1 |
Pagination | 357-360 |
Publisher | Zeta Computational Resources Ltd. |
ISSN Number | 2227-9385 |
Poster
In Silico - Augmented Cardiac Microphysiological Systems for Evaluating Cardiac Drug Effects
Keystone Organ on Chip Symposium, Montana, USA, 2018.Status: Published
In Silico - Augmented Cardiac Microphysiological Systems for Evaluating Cardiac Drug Effects
Afilliation | Scientific Computing |
Publication Type | Poster |
Year of Publication | 2018 |
Place Published | Keystone Organ on Chip Symposium, Montana, USA |
In Silico - Augmented Cardiac Microphysiological Systems for Evaluating Cardiac Drug Effects
BMES Conference, Atlanta, USA, 2018.Status: Published
In Silico - Augmented Cardiac Microphysiological Systems for Evaluating Cardiac Drug Effects
Afilliation | Scientific Computing |
Publication Type | Poster |
Year of Publication | 2018 |
Place Published | BMES Conference, Atlanta, USA |
Talks, contributed
In Silico Modeling of Cardiac Microphysiological Systems for Evaluating Drug Side Effects
In Heart By Numbers Conference, Berlin, Germany, 2018.Status: Published
In Silico Modeling of Cardiac Microphysiological Systems for Evaluating Drug Side Effects
Afilliation | Scientific Computing |
Publication Type | Talks, contributed |
Year of Publication | 2018 |
Location of Talk | Heart By Numbers Conference, Berlin, Germany |
Journal Article
Inversion and computational maturation of drug response using human stem cell derived cardiomyocytes in microphysiological systems
Nature Scientific Reports 8 (2018).Status: Published
Inversion and computational maturation of drug response using human stem cell derived cardiomyocytes in microphysiological systems
While cardiomyocytes differentiated from human induced pluripotent stems cells (hiPSCs) hold great promise for drug screening, the electrophysiological properties of these cells can be variable and immature, producing results that are significantly different from their human adult counterparts. Here, we describe a computational framework to address this limitation, and show how in silico methods, applied to measurements on immature cardiomyocytes, can be used to both identify drug action and to predict its effect in mature cells. Our synthetic and experimental results indicate that optically obtained waveforms of voltage and calcium from microphysiological systems can be inverted into information on drug ion channel blockage, and then, through assuming functional invariance of proteins during maturation, this data can be used to predict drug induced changes in mature ventricular cells. Together, this pipeline of measurements and computational analysis could significantly improve the ability of hiPSC derived cardiomycocytes to predict dangerous drug side effects.
Afilliation | Scientific Computing |
Project(s) | Center for Biomedical Computing (SFF) |
Publication Type | Journal Article |
Year of Publication | 2018 |
Journal | Nature Scientific Reports |
Volume | 8 |
Number | 17626 |
Date Published | 12/2018 |
Publisher | Springer Nature |
URL | https://doi.org/10.1038/s41598-018-35858-7 |
DOI | 10.1038/s41598-018-35858-7 |
Metabolically-Driven Maturation of hiPSC-Cell Derived Heart-on-a-Chip
{bioRxiv (2018).Status: Submitted
Metabolically-Driven Maturation of hiPSC-Cell Derived Heart-on-a-Chip
Human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CM) are a promising in vitro tool for drug development and disease modeling, but their immature electrophysiology limits their diagnostic utility. Tissue engineering approaches involving aligned and 3D culture enhance hiPSC-CM maturation but are insufficient to induce electrophysiological maturation. We hypothesized that recapitulating post-natal switching of the heart’s primary adenosine triphosphate source from glycolysis to fatty acid oxidation could enhance maturation of hiPSC-CM. We combined hiPSC-CM with microfabrication to create 3D cardiac microphysiological systems (MPS) that enhanced immediate microtissue alignment and tissue specific extracellular matrix production. Using Robust Experimental design, we identified a maturation media that allowed the cardiac MPS to correctly assess false positive and negative drug response. Finally, we employed mathematical modeling and gene expression data to explain the observed changes in electrophysiology and pharmacology of MPS exposed to maturation media. In contrast, the same media had no effects on 2D hiPSC-CM monolayers. These results suggest that systematic combination of biophysical stimuli and metabolic cues can enhance the electrophysiological maturation of hiPSC-derived cardiomyocytes.
Afilliation | Scientific Computing |
Publication Type | Journal Article |
Year of Publication | 2018 |
Journal | {bioRxiv |
Publisher | Cold Spring Harbor Laboratory |
URL | https://www.biorxiv.org/content/early/2018/12/10/485169 |
DOI | 10.1101/485169 |
Journal Article
A cell-based framework for numerical modelling of electrical conduction in cardiac tissue
Frontiers in Physics, Computational Physics 5 (2017).Status: Published
A cell-based framework for numerical modelling of electrical conduction in cardiac tissue
Afilliation | Scientific Computing |
Project(s) | Center for Biomedical Computing (SFF) |
Publication Type | Journal Article |
Year of Publication | 2017 |
Journal | Frontiers in Physics, Computational Physics |
Volume | 5 |
Date Published | 10/2017 |
Publisher | Frontiers |
URL | https://www.frontiersin.org/articles/10.3389/fphy.2017.00048/full?&utm_s... |
DOI | 10.3389/fphy.2017.00048 |
An Evaluation of the Accuracy of Classical Models for Computing the Membrane Potential and Extracellular Potential for Neurons
Frontiers in Computational Neuroscience 11 (2017): 27.Status: Published
An Evaluation of the Accuracy of Classical Models for Computing the Membrane Potential and Extracellular Potential for Neurons
Two mathematical models are part of the foundation of Computational neurophysiology; a) the Cable equation is used to compute the membrane potential of neurons, and, b) volume-conductor theory describes the extracellular potential around neurons. In the standard procedure for computing extracellular potentials, the transmembrane currents are computed by means of a) and the extracellular potentials are computed using an explicit sum over analytical point-current source solutions as prescribed by volume conductor theory. Both models are extremely useful as they allow huge simplifications of the computational efforts involved in computing extracellular potentials. However, there are more accurate, though
computationally very expensive, models available where the potentials inside and outside the neurons are computed simultaneously in a self-consistent scheme. In the present work we explore the accuracy of the classical models a) and b) by comparing them to these more accurate schemes.
The main assumption of a) is that the ephaptic current can be ignored in the derivation of the Cable equation. We find, however, for our examples with stylized neurons, that the ephaptic current is comparable in magnitude to other currents involved in the computations, suggesting that it may be significant – at least in parts of the simulation. The magnitude of the error introduced in the membrane potential is several millivolts, and this error also translates into errors in the predicted extracellular potentials. While the error becomes negligible if we assume the extracellular conductivity to be very large, this assumption is, unfortunately, not easy to justify a priori for all situations of interest.
Afilliation | Scientific Computing |
Project(s) | Center for Biomedical Computing (SFF) |
Publication Type | Journal Article |
Year of Publication | 2017 |
Journal | Frontiers in Computational Neuroscience |
Volume | 11 |
Pagination | 27 |
Publisher | Frontiers Media SA |
ISSN | 1662-5188 |
URL | http://journal.frontiersin.org/article/10.3389/fncom.2017.00027 |
DOI | 10.3389/fncom.2017.00027 |
Master's thesis
An Investigation of Necessary Grid Resolution for Numerical Simulations of Calcium Dynamics in Cardiac Cells
In The University of Oslo. Department of Mathematics, University of Oslo, 2015.Status: Published
An Investigation of Necessary Grid Resolution for Numerical Simulations of Calcium Dynamics in Cardiac Cells
One of the challenges related to modelling calcium dynamics in cardiac cells is the large difference in the length scales involved. The dyad, where important processes take place, is very small compared to the whole cell. Therefore, resolutions fine enough to capture the details of what happens in the dyad result in huge computational problems for whole-cell simulations, and the exact choice of resolution has a substantial effect on the problem size.
In this thesis, we investigate what grid resolution is necessary to capture the details of what happens in the dyad. We study simple mathematical models of calcium dynamics in the dyad and find analytical solutions to some of these simple models. Numerical simulations of the models are carried out for different resolutions using finite difference methods in 1D and 2D and a finite volume method in 3D. The accuracy of the numerical simulations is then studied by comparing the numerical solutions to analytical solutions and fine-grid numerical solutions, and the results suggest necessary resolutions in the nanometre range.
Afilliation | Scientific Computing |
Project(s) | Center for Biomedical Computing (SFF) |
Publication Type | Master's thesis |
Year of Publication | 2015 |
Degree awarding institution | The University of Oslo |
Publisher | Department of Mathematics, University of Oslo |
URL | https://www.duo.uio.no/handle/10852/45312 |