A single Markov-type kinetic model accounting for the macroscopic currents of all human voltage-gated sodium channel isoforms

Pietro Balbi, Paolo Massobrio, Jeanette Hellgren Kotaleski

Research output: Contribution to journalArticle

3 Citations (Scopus)

Abstract

Modelling ionic channels represents a fundamental step towards developing biologically detailed neuron models. Until recently, the voltage-gated ion channels have been mainly modelled according to the formalism introduced by the seminal works of Hodgkin and Huxley (HH). However, following the continuing achievements in the biophysical and molecular comprehension of these pore-forming transmembrane proteins, the HH formalism turned out to carry limitations and inconsistencies in reproducing the ion-channels electrophysiological behaviour. At the same time, Markov-type kinetic models have been increasingly proven to successfully replicate both the electrophysiological and biophysical features of different ion channels. However, in order to model even the finest non-conducting molecular conformational change, they are often equipped with a considerable number of states and related transitions, which make them computationally heavy and less suitable for implementation in conductance-based neurons and large networks of those. In this purely modelling study we develop a Markov-type kinetic model for all human voltage-gated sodium channels (VGSCs). The model framework is detailed, unifying (i.e., it accounts for all ion-channel isoforms) and computationally efficient (i.e. with a minimal set of states and transitions). The electrophysiological data to be modelled are gathered from previously published studies on whole-cell patch-clamp experiments in mammalian cell lines heterologously expressing the human VGSC subtypes (from NaV1.1 to NaV1.9). By adopting a minimum sequence of states, and using the same state diagram for all the distinct isoforms, the model ensures the lightest computational load when used in neuron models and neural networks of increasing complexity. The transitions between the states are described by original ordinary differential equations, which represent the rate of the state transitions as a function of voltage (i.e., membrane potential). The kinetic model, developed in the NEURON simulation environment, appears to be the simplest and most parsimonious way for a detailed phenomenological description of the human VGSCs electrophysiological behaviour.

Original languageEnglish
Pages (from-to)e1005737
JournalPLoS Computational Biology
Volume13
Issue number9
DOIs
Publication statusPublished - Sep 2017

Fingerprint

Voltage-Gated Sodium Channels
sodium channels
Kinetic Model
Ion Channels
Sodium
Protein Isoforms
Voltage
sodium
kinetics
Kinetics
Electric potential
Neuron Model
ion channels
Neurons
NAV1.1 Voltage-Gated Sodium Channel
Patch Clamp
ion
Ions
Membrane Potential
Neural Networks (Computer)

Keywords

  • Cell Line
  • Computational Biology
  • Computer Simulation
  • Humans
  • Kinetics
  • Markov Chains
  • Models, Neurological
  • Protein Isoforms
  • Voltage-Gated Sodium Channels
  • Journal Article

Cite this

A single Markov-type kinetic model accounting for the macroscopic currents of all human voltage-gated sodium channel isoforms. / Balbi, Pietro; Massobrio, Paolo; Hellgren Kotaleski, Jeanette.

In: PLoS Computational Biology, Vol. 13, No. 9, 09.2017, p. e1005737.

Research output: Contribution to journalArticle

@article{2b7e08aaa6a04ffcaa45943e7af7e71c,
title = "A single Markov-type kinetic model accounting for the macroscopic currents of all human voltage-gated sodium channel isoforms",
abstract = "Modelling ionic channels represents a fundamental step towards developing biologically detailed neuron models. Until recently, the voltage-gated ion channels have been mainly modelled according to the formalism introduced by the seminal works of Hodgkin and Huxley (HH). However, following the continuing achievements in the biophysical and molecular comprehension of these pore-forming transmembrane proteins, the HH formalism turned out to carry limitations and inconsistencies in reproducing the ion-channels electrophysiological behaviour. At the same time, Markov-type kinetic models have been increasingly proven to successfully replicate both the electrophysiological and biophysical features of different ion channels. However, in order to model even the finest non-conducting molecular conformational change, they are often equipped with a considerable number of states and related transitions, which make them computationally heavy and less suitable for implementation in conductance-based neurons and large networks of those. In this purely modelling study we develop a Markov-type kinetic model for all human voltage-gated sodium channels (VGSCs). The model framework is detailed, unifying (i.e., it accounts for all ion-channel isoforms) and computationally efficient (i.e. with a minimal set of states and transitions). The electrophysiological data to be modelled are gathered from previously published studies on whole-cell patch-clamp experiments in mammalian cell lines heterologously expressing the human VGSC subtypes (from NaV1.1 to NaV1.9). By adopting a minimum sequence of states, and using the same state diagram for all the distinct isoforms, the model ensures the lightest computational load when used in neuron models and neural networks of increasing complexity. The transitions between the states are described by original ordinary differential equations, which represent the rate of the state transitions as a function of voltage (i.e., membrane potential). The kinetic model, developed in the NEURON simulation environment, appears to be the simplest and most parsimonious way for a detailed phenomenological description of the human VGSCs electrophysiological behaviour.",
keywords = "Cell Line, Computational Biology, Computer Simulation, Humans, Kinetics, Markov Chains, Models, Neurological, Protein Isoforms, Voltage-Gated Sodium Channels, Journal Article",
author = "Pietro Balbi and Paolo Massobrio and {Hellgren Kotaleski}, Jeanette",
year = "2017",
month = "9",
doi = "10.1371/journal.pcbi.1005737",
language = "English",
volume = "13",
pages = "e1005737",
journal = "PLoS Computational Biology",
issn = "1553-734X",
publisher = "Public Library of Science",
number = "9",

}

TY - JOUR

T1 - A single Markov-type kinetic model accounting for the macroscopic currents of all human voltage-gated sodium channel isoforms

AU - Balbi, Pietro

AU - Massobrio, Paolo

AU - Hellgren Kotaleski, Jeanette

PY - 2017/9

Y1 - 2017/9

N2 - Modelling ionic channels represents a fundamental step towards developing biologically detailed neuron models. Until recently, the voltage-gated ion channels have been mainly modelled according to the formalism introduced by the seminal works of Hodgkin and Huxley (HH). However, following the continuing achievements in the biophysical and molecular comprehension of these pore-forming transmembrane proteins, the HH formalism turned out to carry limitations and inconsistencies in reproducing the ion-channels electrophysiological behaviour. At the same time, Markov-type kinetic models have been increasingly proven to successfully replicate both the electrophysiological and biophysical features of different ion channels. However, in order to model even the finest non-conducting molecular conformational change, they are often equipped with a considerable number of states and related transitions, which make them computationally heavy and less suitable for implementation in conductance-based neurons and large networks of those. In this purely modelling study we develop a Markov-type kinetic model for all human voltage-gated sodium channels (VGSCs). The model framework is detailed, unifying (i.e., it accounts for all ion-channel isoforms) and computationally efficient (i.e. with a minimal set of states and transitions). The electrophysiological data to be modelled are gathered from previously published studies on whole-cell patch-clamp experiments in mammalian cell lines heterologously expressing the human VGSC subtypes (from NaV1.1 to NaV1.9). By adopting a minimum sequence of states, and using the same state diagram for all the distinct isoforms, the model ensures the lightest computational load when used in neuron models and neural networks of increasing complexity. The transitions between the states are described by original ordinary differential equations, which represent the rate of the state transitions as a function of voltage (i.e., membrane potential). The kinetic model, developed in the NEURON simulation environment, appears to be the simplest and most parsimonious way for a detailed phenomenological description of the human VGSCs electrophysiological behaviour.

AB - Modelling ionic channels represents a fundamental step towards developing biologically detailed neuron models. Until recently, the voltage-gated ion channels have been mainly modelled according to the formalism introduced by the seminal works of Hodgkin and Huxley (HH). However, following the continuing achievements in the biophysical and molecular comprehension of these pore-forming transmembrane proteins, the HH formalism turned out to carry limitations and inconsistencies in reproducing the ion-channels electrophysiological behaviour. At the same time, Markov-type kinetic models have been increasingly proven to successfully replicate both the electrophysiological and biophysical features of different ion channels. However, in order to model even the finest non-conducting molecular conformational change, they are often equipped with a considerable number of states and related transitions, which make them computationally heavy and less suitable for implementation in conductance-based neurons and large networks of those. In this purely modelling study we develop a Markov-type kinetic model for all human voltage-gated sodium channels (VGSCs). The model framework is detailed, unifying (i.e., it accounts for all ion-channel isoforms) and computationally efficient (i.e. with a minimal set of states and transitions). The electrophysiological data to be modelled are gathered from previously published studies on whole-cell patch-clamp experiments in mammalian cell lines heterologously expressing the human VGSC subtypes (from NaV1.1 to NaV1.9). By adopting a minimum sequence of states, and using the same state diagram for all the distinct isoforms, the model ensures the lightest computational load when used in neuron models and neural networks of increasing complexity. The transitions between the states are described by original ordinary differential equations, which represent the rate of the state transitions as a function of voltage (i.e., membrane potential). The kinetic model, developed in the NEURON simulation environment, appears to be the simplest and most parsimonious way for a detailed phenomenological description of the human VGSCs electrophysiological behaviour.

KW - Cell Line

KW - Computational Biology

KW - Computer Simulation

KW - Humans

KW - Kinetics

KW - Markov Chains

KW - Models, Neurological

KW - Protein Isoforms

KW - Voltage-Gated Sodium Channels

KW - Journal Article

U2 - 10.1371/journal.pcbi.1005737

DO - 10.1371/journal.pcbi.1005737

M3 - Article

C2 - 28863150

VL - 13

SP - e1005737

JO - PLoS Computational Biology

JF - PLoS Computational Biology

SN - 1553-734X

IS - 9

ER -