Abstract
The aim of this perspective article is to share with the community of ion channel scientists our thoughts and expectations regarding the increasing role that computational tools will play in the future of our field. The opinions and comments detailed here are the result of a 3-day long international exploratory workshop that took place in October 2013 and that was supported by the Swiss National Science Foundation.
Introduction and current situation
Ion channels are essential proteins that are located at the plasma membrane of virtually all cells. They allow for the passage of ions across the plasma membrane according to their electro-chemical gradients (Hille, 2001). Most of these ion channels show selectivity in their permeability properties (i.e., to sodium, potassium, calcium, or chloride ions found in physiologic fluids). Ion channel gating (opening and closing features) plays a key role in the regulation of ion channel activity. Many ion channel proteins are able to sense the electrical field across the cellular membrane, which is negative at a resting state and inverts upon depolarization of the membrane (Armstrong and Hille, 1998; Bezanilla, 2005). Such channels are said to be voltage-gated, since they open (activate) upon depolarization (Loussouarn and Tarek, 2012). Following depolarization-dependent opening, some of these channels will inactivate and prevent the permeation of further ions.
The essential role of ion channels, particularly in the control of cellular excitability and epithelial transport, was demonstrated at the beginning of the 20th century, even before much was known of their molecular nature (Hille, 2001). In the 1950s, it became clear that ion channels are the targets of very potent drugs that are used to treat a broad range of human pathologies, such as pain, cardiac arrhythmias, hypertension, and epilepsy.
Thanks to the major progress in molecular and human genetics in the 1980s, a myriad of genetic mutations were found in the genes coding for different ion channel subunits. This led to the definition of genetic channelopathies, which groups together all disorders that are primarily caused by dysfunctional mutant ion channels (Ashcroft, 2006). Channelopathies are known to cause epilepsy, pain syndromes, migraines, periodic paralysis, myotonia, cardiac arrhythmias, hypertension, hypotension, cystic fibrosis, and many other rare genetic disorders (Ashcroft, 2006).
Mutations in ion channel subunit genes have been shown to alter ion channel function in many different ways (Cannon, 2007). Mutations can lead to either a loss-of-function (i.e., resulting from a defect in biosynthesis), or to a gain-of-function (i.e., resulting from the alteration of inactivation properties). Peroz et al. (2008), however, demonstrated that many of the molecular dysfunctions do not abide by this simple binary classification. There are, for example, dysfunctions that affect protein targeting, the binding of partner proteins, and the regulation by phospholipids. The demonstration of the important roles of ion channels in numerous diseases, coupled with the fact that ion channels are excellent drug targets because of their location at the cellular membrane, have stimulated the field of ion channel research over the last three decades (Ranjan et al., 2011). It is estimated that about 10,000–15,000 original scientific articles are published every year on the topic of ion channels. Most of the research activity is focused on ion channel-related pathophysiological mechanisms, as well as on finding new diagnostic, preventive, and therapeutic strategies to treat and/or prevent disease.
A broad community of scientists is working to address specific questions about the function and roles of ion channels in disease. Included in this community are biophysicists, biochemists, geneticists, physiologists, pharmacologists, computational biologists, and clinical scientists. They are all tackling these questions with their specific expertise.
Because of the overwhelming amount of data, attempts to integrate this highly fragmented knowledge from vastly different disciplines have been proven to be difficult. One major objective of summarizing the information is to enable its application in clinical practice. Biophysicists, computational biologists, and biochemists are all playing an important role to achieve this objective. The seminal work of Hodgkin and Huxley (1952) on the giant axon of the squid (Nobel Prize 1963) initiated this approach. Applying the strategy of Hodgkin & Huxley, Noble and his colleagues began developing very soon thereafter, a number of models for various types of cardiac cells and tissues (for a review, see Noble et al., 2012). As a result of the developments in information technology over the last few decades, the modeling and simulation of biophysical and biological processes has consequently become much more amenable. Mathematical models of ion channels and excitable tissues have become increasingly refined and complex, while constantly taking into account new experimental findings.
The current trend in ion channel modeling is to represent ion channels as Markovian systems, which generalizes the framework established by Hodgkin and Huxley (1952). Such models recapitulate single channel behavior, whole cell currents, and the action of drugs (Irvine et al., 1999; Clancy et al., 2007; Milescu et al., 2008; Bett et al., 2011). Many models have been developed at the cellular level, notably cardiac cells [e.g., the review by Wilhelms et al. on human atrial cell models (Wilhelms et al., 2013)]. Notably there is an increasing level of structural detail in cell models. For example, cardiac models have started to incorporate the fine structure of the sarcoplasmic reticulum, the T-tubule system and the dyadic spaces (Williams et al., 2011; Nivala et al., 2012). Simulation, using computational methods, is becoming an integral part of the research process in the biomedical sciences, as illustrated by a recent study by Silva et al. (2009). At this stage, these computational approaches have mainly been used to investigate the possible consequences of the mutation-induced alterations of voltage-gated ion channels on the excitability of cells in neurologic disorders such as myotonia (Cannon, 2007), cardiac arrhythmias, and conduction defects (Rudy, 2012). Simulations of more integrative signals have also been performed, such as the recorded cardiac electrical activity on the body surface (ECG) (Rudy, 2012) and the electrical activity of circuits with hundreds of neurons (Hill et al., 2012). An outstanding example of such a current project is the “Human Brain Project,” which has recently been selected as one of the research flagship projects of the European commission (Markram, 2012).
Another important trend is the emerging and very promising field of research using computational frameworks to model molecular function (molecular dynamics), for which the Nobel Prize of Chemistry was awarded to Karplus, Levitt, and Warshel in 2013. These computational approaches are now feasible for very complex molecules, such as proteins and multiprotein complexes, i.e., ion channels (Gosselin-Badaroudine et al., 2012; Tarek and Delemotte, 2013). These tools will most likely play an increasing role for membrane proteins, as the rate of newly identified structures from crystallography experiments is growing exponentially (http://blanco.biomol.uci.edu/mpstruc/), with 433 unique proteins identified as of November 2013.
The future of modeling to study ion channels and channelopathies
A model can be defined as a set of rules, ideally based on the laws of physics and expressed with the formalism of mathematical equations, to describe an observed system. The aim of computer simulation is to implement such a model in a computer program with the goal of replicating the behaviors observed in the real world. To be as close as possible to physical reality, one should strive to formulate models based on established physical principles as long as the experimental data allows for this. However, in practice, this is not always possible, and one often recourses to phenomenological models to describe ion channel behavior. Nevertheless, such models are still highly valuable.
Future modeling developments will most likely be performed using multidisciplinary and multiscale approaches that encompass the modeling of molecules (molecular dynamics), the modeling of ion channel microscopic kinetics, the modeling of cellular electrical activity, and finally, the modeling of the electrical properties of tissues and organs. There is a need to bridge the gap between molecular dynamics simulations and kinetic models of ion channels. The derivation of energy landscapes from molecular dynamics data, and thus, the rate constants in Markovian kinetic models, is an excellent approach that was recently pioneered by Rudy's group, (Silva et al., 2009; Nekouzadeh and Rudy, 2011) where effects of mutations in the KCNQ1 -mediated current were accurately predicted, purely based upon protein sequence. Whether or not such modeling should be performed at all levels remains uncertain. Multiscale modeling is a synergetic process, thus, it seems reasonable that all levels should be modeled and integrated to warrant an optimal yield of information.
Collaboration and multidisciplinarity are essential
There is a tremendously large amount of data available on ion channel biology and function, as well as on channelopathies. There is a critical need to regroup this information into well-constructed databases, which incorporate as much detailed and quality information about ion channels and their generated currents as possible. Collaboration between experimentalists and modeling scientists must be encouraged in order to obtain an optimal yield of information to produce relevant new knowledge. Both scientific individuals, if not united in the same person, and groups should work closely together on projects. Collaborative platforms and interfaces could also encourage a structured exchange of information in both directions. Psychological barriers that may be present on either side could be relieved by working together in the same research group or consortium.
Current cellular electrophysiology protocols
Currently, most of the electrophysiology protocols that are used to record the activity of ion channels from whole-cell recordings do not have the most optimized design to enable the collection of all the intrinsic dynamics of ion channels. There is a need to generate, refine, and use non-conventional cellular electrophysiology protocols which subject the ion channels to continuous perturbations that increase the information obtained from such experiments. An inspiring example is the “dichotomous noise” protocol proposed by Millonas and Hanck (1998). Compared to models based on the Hodgkin-Huxley formalism, Markov models also bring additional challenges in terms of parameter estimation, identifiability, and computational speed. These challenges call for the development of innovative methods to estimate model parameters, to formulate Markov models that are mathematically identifiable, and to allow automatized identification from experimental data that are obtained with optimized voltage clamp protocols. Several research groups have already started to address these different issues (Bruno et al., 2005; Fink and Noble, 2009; Siekmann et al., 2012). All of these approaches could be aptly named “Systems Biophysics,” and should be further developed. Furthermore, the experimentalists should consider “full ion channel (kinetic) models” when communicating their data. Raw data should also be made available, rather than only a summary of activation and inactivation curves, for example. If this had been done in past publications, it may be postulated that there would not have been a need for the “channelome” project of the Swiss Institute of Technology of Lausanne (Ranjan, 2011).
“variability” as a parameter to take into account
Experimental results show “variability” in responses following administration of a given stimulus to biological samples. The sources of this variability are diverse and may reflect either intrinsic biological properties or may result from the recording set-up. With few exceptions (Sato et al., 2010; Lemay et al., 2011), the variability and its consequences have not been adequately addressed by the current models which oftentimes calculate averaged values. This problem must be addressed in the future.
Genetic complexity
It is conceptually simple to model the consequences of mutations in ion channel genes which cause well-defined phenotypes, such as inherited epilepsy or congenital long QT syndrome, that are currently described as monogenic disorders. However, this simple paradigm of “one mutation causing one phenotype” does not hold true in many cases (Probst et al., 2009; Klassen et al., 2011; Bezzina et al., 2013). As a consequence, these observations should be used in future models by incorporating gene-gene and gene-environment interactions.
Conclusions
In conclusion, we are convinced that modeling, through the use of more and more sophisticated computational tools, will prove to be an essential component for the investigation of channelopathies. In some fields of experimental science such as chemistry, computer modeling is already an essential part of the publication process. In the near future, no one will be able to simply say, “I do not believe in computer models!”
Conflict of interest statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Statements
Acknowledgments
This perspective article was drafted during an international exploratory workshop that was partially financed by the Swiss National Science Foundation (IZ32Z0_150895) held on October 9–12, 2013, in Beatenberg (Switzerland). Gildas Loussouarn is supported by a Marie Curie International Outgoing Fellowship within the 7th European Community Framework Programme, and the AFMTéléthon. The authors participated in the workshop and would like to thank the following young scientists for their active participation: Mohamed-Yassine Amarouch (University of Bern), Fayal Abderemane-Ali (University of Nantes), Séverine Petitprez and Rajnish Ranjan (both Swiss Federal Institute of Technology, Lausanne), Marina Kasimova and Magdalena Rucinska (both CNRS-University of Lorraine). We also thank Dr. A. Felley and D. Shy for their comments on the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
ArmstrongC. M.HilleB. (1998). Voltage-gated ion channels and electrical excitability. Neuron20, 371–380. 10.1016/S0896-6273(00)80981-2
2
AshcroftF. M. (2006). From molecule to malady. Nature440, 440–447. 10.1038/nature04707
3
BettG. C.ZhouQ.RasmussonR. L. (2011). Models of HERG gating. Biophys. J. 101, 631–642. 10.1016/j.bpj.2011.06.050
4
BezanillaF. (2005). Voltage-gated ion channels. IEEE Trans. Nanobiosci. 4, 34–48. 10.1109/TNB.2004.842463
5
BezzinaC. R.BarcJ.MizusawaY.RemmeC. A.GourraudJ. B.SimonetF.et al. (2013). Common variants at SCN5A-SCN10A and HEY2 are associated with Brugada syndrome, a rare disease with high risk of sudden cardiac death. Nat. Genet. 45, 1044–1049. 10.1038/ng.2712
6
BrunoW. J.YangJ.PearsonJ. E. (2005). Using independent open-to-closed transitions to simplify aggregated Markov models of ion channel gating kinetics. Proc. Natl. Acad. Sci. U.S.A. 102, 6326–6331. 10.1073/pnas.0409110102
7
CannonS. C. (2007). Physiologic principles underlying ion channelopathies. Neurotherapeutics4, 174–183. 10.1016/j.nurt.2007.01.015
8
ClancyC. E.ZhuZ. I.RudyY. (2007). Pharmacogenetics and anti-arrhythmic drug therapy: a theoretical investigation. Am. J. Physiol. Heart Circ. Physiol. 292, H66–H75. 10.1152/ajpheart.00312.2006
9
FinkM.NobleD. (2009). Markov models for ion channels: versatility versus identifiability and speed. Philos. Trans. A Math. Phys. Eng. Sci. 367, 2161–2179. 10.1098/rsta.2008.0301
10
Gosselin-BadaroudineP.DelemotteL.MoreauA.KleinM. L.ChahineM. (2012). Gating pore currents and the resting state of Nav1.4 voltage sensor domains. Proc. Natl. Acad. Sci. U.S.A. 109, 19250–19255. 10.1073/pnas.1217990109
11
HillS. L.WangY.RiachiI.SchürmannF.MarkramH. (2012). Statistical connectivity provides a sufficient foundation for specific functional connectivity in neocortical neural microcircuits. PNAS109, E2885–E2894. 10.1073/pnas.1202128109
12
HilleB. (2001). Ion Channels of Excitable Membranes. 3rd Edn. Massachusetts: Sinauer, Sunderland.
13
HodgkinA. L.HuxleyA. F. (1952). A quantitative description of membrane current and its application to conduction and excitation in nerve. J. Physiol. 117, 500–544.
14
IrvineL. A.JafriM. S.WinslowR. L. (1999). Cardiac sodium channel Markov model with temperature dependence and recovery from inactivation. Biophys. J. 76, 1868–1885. 10.1016/S0006-3495(99)77346-7
15
KlassenT.DavisC.GoldmanA.BurgessD.ChenT.WheelerD.et al. (2011). Exome sequencing of ion channel genes reveals complex profiles confounding personal risk assessment in epilepsy. Cell145, 1036–1048. 10.1016/j.cell.2011.05.025
16
LemayM.DeL. E.KuceraJ. P. (2011). Effects of stochastic channel gating and distribution on the cardiac action potential. J. Theor. Biol. 281, 84–96. 10.1016/j.jtbi.2011.04.019
17
LoussouarnG.TarekM. (2012). Mechanisms of ion channels voltage-dependency: all about molecular sensors, gates, levers, locks, and grease. Front. Pharmacol. 3:174. 10.3389/fphar.2012.00174
18
MarkramH. (2012). The human brain project. Sci. Am. 306, 50–55. 10.1038/scientificamerican0612-50
19
MilescuL. S.YamanishiT.PtakK.MogriM. Z.SmithJ. C. (2008). Real-time kinetic modeling of voltage-gated ion channels using dynamic clamp. Biophys. J. 95, 66–87. 10.1529/biophysj.107.118190
20
MillonasM. M.HanckD. A. (1998). Non-equilibrium response spectroscopy of voltage-sensitive ion channel gating. Biophys. J. 74, 210–229. 10.1016/S0006-3495(98)77781-1
21
NekouzadehA.RudyY. (2011). Continuum molecular simulation of large conformational changes during ion-channel gating. PLoS ONE6:e20186. 10.1371/journal.pone.0020186
22
NivalaM.DeL. E.RovettiR.QuZ. (2012). Computational modeling and numerical methods for spatiotemporal calcium cycling in ventricular myocytes. Front. Physiol. 3:114. 10.3389/fphys.2012.00114
23
NobleD.GarnyA.NobleP. J. (2012). How the Hodgkin-Huxley equations inspired the cardiac physiome project. J. Physiol. 590, 2613–2628. 10.1113/jphysiol.2011.224238
24
PerozD.RodriguezN.ChoveauF.BaroI.MerotJ.LoussouarnG. (2008). Kv7.1 (KCNQ1) properties and channelopathies. J. Physiol. 586, 1785–1789. 10.1113/jphysiol.2007.148254
25
ProbstV.WildeA. A. M.BarcJ.SacherF.BabutyD.MaboP.et al. (2009). SCN5A mutations and the role of genetic background in the pathophysiology of Brugada syndrome. Circ. Cardiovasc. Genet. 2, 552–557. 10.1161/CIRCGENETICS.109.853374
26
RanjanR. (2011). Engineering Neuron Models: from Ion Channels to Electrical Behavior. Available online at: biblion.epfl.ch/EPFL/theses/2011/5129/5129_abs.pdf
27
RanjanR.KhazenG.GambazziL.RamaswamyS.HillS. L.SchürmannF.et al. (2011). Channelpedia: an integrative and interactive database for ion channels. Front. Neuroinform. 5:36. 10.3389/fninf.2011.00036
28
RudyY. (2012). From genes and molecules to organs and organisms: heart, in Comprehensive Biophysics, ed WeinsteinH. (Oxford: Academic Press), 268–327. 10.1016/B978-0-12-374920-8.00924-3
29
SatoD.XieL. H.NguyenT. P.WeissJ. N.QuZ. (2010). Irregularly appearing early afterdepolarizations in cardiac myocytes: random fluctuations or dynamical chaos?Biophys. J. 99, 765–773. 10.1016/j.bpj.2010.05.019
30
SiekmannI.SneydJ.CrampinE. J. (2012). MCMC can detect non-identifiable models. Biophys. J. 103, 2275–2286. 10.1016/j.bpj.2012.10.024
31
SilvaJ. R.PanH.WuD.NekouzadehA.DeckerK. F.CuiJ.et al. (2009). A multiscale model linking ion-channel molecular dynamics and electrostatics to the cardiac action potential. Proc. Natl. Acad. Sci. U.S.A. 106, 11102–11106. 10.1073/pnas.0904505106
32
TarekM.DelemotteL. (2013). Omega currents in voltage-gated ion channels: what can we learn from uncovering the voltage-sensing mechanism using MD simulations?Acc. Chem. Res. 46, 2755–2762. 10.1021/ar300290u
33
WilhelmsM.HettmannH.MaleckarM. M. C.KoivumäkiJ. T.DösselO.SeemannG. (2013). Benchmarking electrophysiological models of human atrial myocytes. Front. Physiol. 3:487. 10.3389/fphys.2012.00487
34
WilliamsG. S.ChikandoA. C.TuanH. T.SobieE. A.LedererW. J.JafriM. S. (2011). Dynamics of calcium sparks and calcium leak in the heart. Biophys. J. 101, 1287–1296. 10.1016/j.bpj.2011.07.021
Summary
Keywords
ion channels, channelopathies, computer models, simulation, molecular dynamics
Citation
Abriel H, de Lange E, Kucera JP, Loussouarn G and Tarek M (2013) Computational tools to investigate genetic cardiac channelopathies. Front. Physiol. 4:390. doi: 10.3389/fphys.2013.00390
Received
16 November 2013
Accepted
10 December 2013
Published
26 December 2013
Volume
4 - 2013
Edited by
Marcel van der Heyden, University Medical Center, Netherlands
Reviewed by
Anna Weinzinger, University of Vienna, Austria; Gary R. Mirams, University of Oxford, UK
Copyright
© 2013 Abriel, de Lange, Kucera, Loussouarn and Tarek.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hugues Abriel, Department of Clinical Research, University of Bern, Murtenstrasse 35, 3010 Bern, Switzerland e-mail: hugues.abriel@dkf.unibe.ch
This article was submitted to Cardiac Electrophysiology, a section of the journal Frontiers in Physiology.
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.