Computational tools to investigate genetic cardiac channelopathies
- 1Department of Clinical Research, University of Bern, Bern, Switzerland
- 2Department of Knowledge Engineering, Maastricht University, Maastricht, Netherlands
- 3Department of Physiology, University of Bern, Bern, Switzerland
- 4INSERM, UMR 1087, l'Institut du thorax, Nantes, France
- 5Centre National de la Recherche Scientifique, L'Institut du Thorax, UMR 6921, Nantes, France
- 6L'Institut du Thorax, UMR 6921, Université de Nantes, Nantes, France
- 7Theory, Modeling and Simulations, UMR 7565, Université de Lorraine, Vandoeuvre-lés-Nancy, France
- 8Theory, Modeling and Simulations, UMR 7565, Centre National de la Recherche Scientifique, Vandoeuvre-lés-Nancy, France
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.
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.
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.
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.
Bezzina, C. R., Barc, J., Mizusawa, Y., Remme, C. A., Gourraud, J. B., Simonet, F., 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. doi: 10.1038/ng.2712
Bruno, W. J., Yang, J., and Pearson, J. 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. doi: 10.1073/pnas.0409110102
Clancy, C. E., Zhu, Z. I., and Rudy, Y. (2007). Pharmacogenetics and anti-arrhythmic drug therapy: a theoretical investigation. Am. J. Physiol. Heart Circ. Physiol. 292, H66–H75. doi: 10.1152/ajpheart.00312.2006
Gosselin-Badaroudine, P., Delemotte, L., Moreau, A., Klein, M. L., and Chahine, M. (2012). Gating pore currents and the resting state of Nav1.4 voltage sensor domains. Proc. Natl. Acad. Sci. U.S.A. 109, 19250–19255. doi: 10.1073/pnas.1217990109
Hill, S. L., Wang, Y., Riachi, I., Schürmann, F., and Markram, H. (2012). Statistical connectivity provides a sufficient foundation for specific functional connectivity in neocortical neural microcircuits. PNAS 109, E2885–E2894. doi: 10.1073/pnas.1202128109
Irvine, L. A., Jafri, M. S., and Winslow, R. L. (1999). Cardiac sodium channel Markov model with temperature dependence and recovery from inactivation. Biophys. J. 76, 1868–1885. doi: 10.1016/S0006-3495(99)77346-7
Klassen, T., Davis, C., Goldman, A., Burgess, D., Chen, T., Wheeler, D., et al. (2011). Exome sequencing of ion channel genes reveals complex profiles confounding personal risk assessment in epilepsy. Cell 145, 1036–1048. doi: 10.1016/j.cell.2011.05.025
Loussouarn, G., and Tarek, M. (2012). Mechanisms of ion channels voltage-dependency: all about molecular sensors, gates, levers, locks, and grease. Front. Pharmacol. 3:174. doi: 10.3389/fphar.2012.00174
Milescu, L. S., Yamanishi, T., Ptak, K., Mogri, M. Z., and Smith, J. C. (2008). Real-time kinetic modeling of voltage-gated ion channels using dynamic clamp. Biophys. J. 95, 66–87. doi: 10.1529/biophysj.107.118190
Nivala, M., De, L. E., Rovetti, R., and Qu, Z. (2012). Computational modeling and numerical methods for spatiotemporal calcium cycling in ventricular myocytes. Front. Physiol. 3:114. doi: 10.3389/fphys.2012.00114
Probst, V., Wilde, A. A. M., Barc, J., Sacher, F., Babuty, D., Mabo, P., et al. (2009). SCN5A mutations and the role of genetic background in the pathophysiology of Brugada syndrome. Circ. Cardiovasc. Genet. 2, 552–557. doi: 10.1161/CIRCGENETICS.109.853374
Ranjan, R. (2011). Engineering Neuron Models: from Ion Channels to Electrical Behavior. Available online at: biblion.epfl.ch/EPFL/theses/2011/5129/5129_abs.pdf
Ranjan, R., Khazen, G., Gambazzi, L., Ramaswamy, S., Hill, S. L., Schürmann, F., et al. (2011). Channelpedia: an integrative and interactive database for ion channels. Front. Neuroinform. 5:36. doi: 10.3389/fninf.2011.00036
Sato, D., Xie, L. H., Nguyen, T. P., Weiss, J. N., and Qu, Z. (2010). Irregularly appearing early afterdepolarizations in cardiac myocytes: random fluctuations or dynamical chaos? Biophys. J. 99, 765–773. doi: 10.1016/j.bpj.2010.05.019
Silva, J. R., Pan, H., Wu, D., Nekouzadeh, A., Decker, K. F., Cui, J., 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. doi: 10.1073/pnas.0904505106
Tarek, M., and Delemotte, L. (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. doi: 10.1021/ar300290u
Wilhelms, M., Hettmann, H., Maleckar, M. M. C., Koivumäki, J. T., Dössel, O., and Seemann, G. (2013). Benchmarking electrophysiological models of human atrial myocytes. Front. Physiol. 3:487. doi: 10.3389/fphys.2012.00487
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 online: 26 December 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: firstname.lastname@example.org