Resolving artefacts in voltage-clamp experiments with computational modelling: an application to fast sodium current recordings.

Cellular electrophysiology is the foundation of many fields, from basic science in neurology, cardiology, oncology to safety critical applications for drug safety testing, clinical phenotyping, etc. Patch-clamp voltage clamp is the gold standard technique for studying cellular electrophysiology. Yet, the quality of these experiments is not always transparent, which may lead to erroneous conclusions for studies and applications. Here, we have developed a new computational approach that allows us to explain and predict the experimental artefacts in voltage-clamp experiments. The computational model captures the experimental procedure and its inadequacies, including: voltage offset, series resistance, membrane capacitance and (imperfect) amplifier compensations, such as series resistance compensation and supercharging. The computational model was validated through a series of electrical model cell experiments. Using this computational approach, the artefacts in voltage-clamp experiments of cardiac fast sodium current, one of the most challenging currents to voltage clamp, were able to be resolved and explained through coupling the observed current and the simulated membrane voltage, including some typically observed shifts and delays in the recorded currents. We further demonstrated that the typical way of averaging data for current-voltage relationships would lead to biases in the peak current and shifts in the peak voltage, and such biases can be in the same order of magnitude as those differences reported for disease-causing mutations. Therefore, the presented new computational pipeline will provide a new standard of assessing the voltage-clamp experiments and interpreting the experimental data, which may be able to rectify and provide a better understanding of ion channel mutations and other related applications.

Recording ten-fold larger IKr conductances with automated patch clamping using equimolar Cs+ solutions.

Background: The rapid delayed rectifier potassium current (IKr) is important for cardiac repolarization and is most often involved in drug-induced arrhythmias. However, accurately measuring this current can be challenging in human-induced pluripotent stem cell (hiPSC)-derived cardiomyocytes because of its small current density. Interestingly, the ion channel conducting IKr, hERG channel, is not only permeable to K+ ions but also to Cs+ ions when present in equimolar concentrations inside and outside of the cell. Methods: In this study, IhERG was measured from Chinese hamster ovary (CHO)-hERG cells and hiPSC-CM using either Cs+ or K+ as the charge carrier. Equimolar Cs+ has been used in the literature in manual patch-clamp experiments, and here, we apply this approach using automated patch-clamp systems. Four different (pre)clinical drugs were tested to compare their effects on Cs+- and K+-based currents. Results: Using equimolar Cs+ solutions gave rise to approximately ten-fold larger hERG conductances. Comparison of Cs+- and K+-mediated currents upon application of dofetilide, desipramine, moxifloxacin, or LUF7244 revealed many similarities in inhibition or activation properties of the drugs studied. Using equimolar Cs+ solutions gave rise to approximately ten-fold larger hERG conductances. In hiPSC-CM, the Cs+-based conductance is larger compared to the known K+-based conductance, and the Cs+ hERG conductance can be inhibited similarly to the K+-based conductance. Conclusion: Using equimolar Cs+ instead of K+ for IhERG measurements in an automated patch-clamp system gives rise to a new method by which, for example, quick scans can be performed on effects of drugs on hERG currents. This application is specifically relevant when such experiments are performed using cells which express small IKr current densities in combination with small membrane capacitances.

Inflammation-induced mitochondrial and metabolic disturbances in sensory neurons control the switch from acute to chronic pain.

Pain often persists in patients with an inflammatory disease, even when inflammation has subsided. The molecular mechanisms leading to this failure in pain resolution and the transition to chronic pain are poorly understood. Mitochondrial dysfunction in sensory neurons links to chronic pain, but its role in resolution of inflammatory pain is unclear. Transient inflammation causes neuronal plasticity, called hyperalgesic priming, which impairs resolution of pain induced by a subsequent inflammatory stimulus. We identify that hyperalgesic priming in mice increases the expression of a mitochondrial protein (ATPSc-KMT) and causes mitochondrial and metabolic disturbances in sensory neurons. Inhibition of mitochondrial respiration, knockdown of ATPSCKMT expression, or supplementation of the affected metabolite is sufficient to restore resolution of inflammatory pain and prevents chronic pain development. Thus, inflammation-induced mitochondrial-dependent disturbances in sensory neurons predispose to a failure in resolution of inflammatory pain and development of chronic pain.

Leak current, even with gigaohm seals, can cause misinterpretation of stem cell-derived cardiomyocyte action potential recordings.

AIMS: Human-induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) have become an essential tool to study arrhythmia mechanisms. Much of the foundational work on these cells, as well as the computational models built from the resultant data, has overlooked the contribution of seal-leak current on the immature and heterogeneous phenotype that has come to define these cells. The aim of this study is to understand the effect of seal-leak current on recordings of action potential (AP) morphology. METHODS AND RESULTS: Action potentials were recorded in human iPSC-CMs using patch clamp and simulated using previously published mathematical models. Our in silico and in vitro studies demonstrate how seal-leak current depolarizes APs, substantially affecting their morphology, even with seal resistances (Rseal) above 1 GΩ. We show that compensation of this leak current is difficult due to challenges with obtaining accurate measures of Rseal during an experiment. Using simulation, we show that Rseal measures (i) change during an experiment, invalidating the use of pre-rupture values, and (ii) are polluted by the presence of transmembrane currents at every voltage. Finally, we posit that the background sodium current in baseline iPSC-CM models imitates the effects of seal-leak current and is increased to a level that masks the effects of seal-leak current on iPSC-CMs. CONCLUSION: Based on these findings, we make recommendations to improve iPSC-CM AP data acquisition, interpretation, and model-building. Taking these recommendations into account will improve our understanding of iPSC-CM physiology and the descriptive ability of models built from such data.

Interplay between calcium and sarcomeres directs cardiomyocyte maturation during regeneration.

Zebrafish hearts can regenerate by replacing damaged tissue with new cardiomyocytes. Although the steps leading up to the proliferation of surviving cardiomyocytes have been extensively studied, little is known about the mechanisms that control proliferation and redifferentiation to a mature state. We found that the cardiac dyad, a structure that regulates calcium handling and excitation-contraction coupling, played a key role in the redifferentiation process. A component of the cardiac dyad called leucine-rich repeat-containing 10 (Lrrc10) acted as a negative regulator of proliferation, prevented cardiomegaly, and induced redifferentiation. We found that its function was conserved in mammalian cardiomyocytes. This study highlights the importance of the underlying mechanisms required for heart regeneration and their application to the generation of fully functional cardiomyocytes.

PITX2 induction leads to impaired cardiomyocyte function in arrhythmogenic cardiomyopathy.

Arrhythmogenic cardiomyopathy (ACM) is an inherited progressive disease characterized by electrophysiological and structural remodeling of the ventricles. However, the disease-causing molecular pathways, as a consequence of desmosomal mutations, are poorly understood. Here, we identified a novel missense mutation within desmoplakin in a patient clinically diagnosed with ACM. Using CRISPR-Cas9, we corrected this mutation in patient-derived human induced pluripotent stem cells (hiPSCs) and generated an independent knockin hiPSC line carrying the same mutation. Mutant cardiomyocytes displayed a decline in connexin 43, NaV1.5, and desmosomal proteins, which was accompanied by a prolonged action potential duration. Interestingly, paired-like homeodomain 2 (PITX2), a transcription factor that acts a repressor of connexin 43, NaV1.5, and desmoplakin, was induced in mutant cardiomyocytes. We validated these results in control cardiomyocytes in which PITX2 was either depleted or overexpressed. Importantly, knockdown of PITX2 in patient-derived cardiomyocytes is sufficient to restore the levels of desmoplakin, connexin 43, and NaV1.5.

The zebrafish grime mutant uncovers an evolutionarily conserved role for Tmem161b in the control of cardiac rhythm.

The establishment of cardiac function in the developing embryo is essential to ensure blood flow and, therefore, growth and survival of the animal. The molecular mechanisms controlling normal cardiac rhythm remain to be fully elucidated. From a forward genetic screen, we identified a unique mutant, grime, that displayed a specific cardiac arrhythmia phenotype. We show that loss-of-function mutations in tmem161b are responsible for the phenotype, identifying Tmem161b as a regulator of cardiac rhythm in zebrafish. To examine the evolutionary conservation of this function, we generated knockout mice for Tmem161b. Tmem161b knockout mice are neonatal lethal and cardiomyocytes exhibit arrhythmic calcium oscillations. Mechanistically, we find that Tmem161b is expressed at the cell membrane of excitable cells and live imaging shows it is required for action potential repolarization in the developing heart. Electrophysiology on isolated cardiomyocytes demonstrates that Tmem161b is essential to inhibit Ca2+ and K+ currents in cardiomyocytes. Importantly, Tmem161b haploinsufficiency leads to cardiac rhythm phenotypes, implicating it as a candidate gene in heritable cardiac arrhythmia. Overall, these data describe Tmem161b as a highly conserved regulator of cardiac rhythm that functions to modulate ion channel activity in zebrafish and mice.

Is zebrafish heart regeneration "complete"? Lineage-restricted cardiomyocytes proliferate to pre-injury numbers but some fail to differentiate in fibrotic hearts.

Adult zebrafish are frequently described to be able to "completely" regenerate the heart. Yet, the extent to which cardiomyocytes lost to injury are replaced is unknown, since existing evidence for cardiomyocyte proliferation is indirect or non-quantitative. We established stereological methods to quantify the number of cardiomyocytes at several time-points post cryoinjury. Intriguingly, after cryoinjuries that killed about 1/3 of the ventricular cardiomyocytes, pre-injury cardiomyocyte numbers were restored already within 30 days. Yet, many hearts retained small residual scars, and a subset of cardiomyocytes bordering these fibrotic areas remained smaller, lacked differentiated sarcomeric structures, and displayed defective calcium signaling. Thus, a subset of regenerated cardiomyocytes failed to fully mature. While lineage-tracing experiments have shown that regenerating cardiomyocytes are derived from differentiated cardiomyocytes, technical limitations have previously made it impossible to test whether cardiomyocyte trans-differentiation contributes to regeneration of non-myocyte cell lineages. Using Cre responder lines that are expressed in all major cell types of the heart, we found no evidence for cardiomyocyte transdifferentiation into endothelial, epicardial, fibroblast or immune cell lineages. Overall, our results imply a refined answer to the question whether zebrafish can completely regenerate the heart: in response to cryoinjury, preinjury cardiomyocyte numbers are indeed completely regenerated by proliferation of lineage-restricted cardiomyocytes, while restoration of cardiomyocyte differentiation and function, as well as resorption of scar tissue, is less robustly achieved.

Accounting for variability in ion current recordings using a mathematical model of artefacts in voltage-clamp experiments.

Mathematical models of ion channels, which constitute indispensable components of action potential models, are commonly constructed by fitting to whole-cell patch-clamp data. In a previous study, we fitted cell-specific models to hERG1a (Kv11.1) recordings simultaneously measured using an automated high-throughput system, and studied cell-cell variability by inspecting the resulting model parameters. However, the origin of the observed variability was not identified. Here, we study the source of variability by constructing a model that describes not just ion current dynamics, but the entire voltage-clamp experiment. The experimental artefact components of the model include: series resistance, membrane and pipette capacitance, voltage offsets, imperfect compensations made by the amplifier for these phenomena, and leak current. In this model, variability in the observations can be explained by either cell properties, measurement artefacts, or both. Remarkably, by assuming that variability arises exclusively from measurement artefacts, it is possible to explain a larger amount of the observed variability than when assuming cell-specific ion current kinetics. This assumption also leads to a smaller number of model parameters. This result suggests that most of the observed variability in patch-clamp data measured under the same conditions is caused by experimental artefacts, and hence can be compensated for in post-processing by using our model for the patch-clamp experiment. This study has implications for the question of the extent to which cell-cell variability in ion channel kinetics exists, and opens up routes for better correction of artefacts in patch-clamp data. This article is part of the theme issue 'Uncertainty quantification in cardiac and cardiovascular modelling and simulation'.

Development, Implementation and Testing of a Multicellular Dynamic Action Potential Clamp Simulator for Drug Cardiac Safety Assessment.

As drugs can be multichannel blockers it is important to assess their cardiac safety taking into account multiple currents. In silico action potential (AP) models have been proposed for being able to integrate drugs effect on ionic currents and generate the resulting AP. However, a mathematical description of drug effects is required, which could be inaccurate. Dynamic Clamp has been proposed for drug cardiac safety assessment. In the dynamic action potential clamp (dAPC) configuration it creates an hybrid model connecting a real cell with a computer simulation. This way, drugs could be administrated directly to real cells, and effects on currents can be taken into account when generating the AP. Here we design and simulate a parallel multichannel dAPC system. The system includes the real cells overexpressing the currents of interest, the voltage clamp acquisition system, and the AP in silico model.

A nonlinear and time-dependent leak current in the presence of calcium fluoride patch-clamp seal enhancer.

Automated patch-clamp platforms are widely used and vital tools in both academia and industry to enable high-throughput studies such as drug screening. A leak current to ground occurs whenever the seal between a pipette and cell (or internal solution and cell in high-throughput machines) is not perfectly insulated from the bath (extracellular) solution. Over 1 GΩ seal resistance between pipette and bath solutions is commonly used as a quality standard for manual patch work. With automated platforms it can be difficult to obtain such a high seal resistance between the intra- and extra-cellular solutions. One suggested method to alleviate this problem is using an F - containing internal solution together with a Ca 2+ containing external solution - so that a CaF 2 crystal forms when the two solutions meet which 'plugs the holes' to enhance the seal resistance. However, we observed an unexpected nonlinear-in-voltage and time-dependent current using these solutions on an automated patch-clamp platform. We performed manual patch-clamp experiments with the automated patch-clamp solutions, but no biological cell, and observed the same nonlinear time-dependent leak current. The current could be completely removed by washing out F - ions to leave a conventional leak current that was linear and not time-dependent. We therefore conclude fluoride ions interacting with the CaF 2 crystal are the origin of the nonlinear time-dependent leak current. The consequences of such a nonlinear and time-dependent leak current polluting measurements should be considered carefully if it cannot be isolated and subtracted.

Required GK1 to Suppress Automaticity of iPSC-CMs Depends Strongly on IK1 Model Structure.

Human-induced pluripotent stem cells derived cardiomyocytes (hiPSC-CMs) are a virtually endless source of human cardiomyocytes that may become a great tool for safety pharmacology; however, their electrical phenotype is immature: they show spontaneous action potentials (APs) and an unstable and depolarized resting membrane potential (RMP) because of lack of IK1. Such immaturity hampers their application in assessing drug safety. The electronic overexpression of IK1 (e.g., through the dynamic clamp (DC) technique) is an option to overcome this deficit. In this computational study, we aim to estimate how much IK1 is needed to bring hiPSC-CMs to a stable and hyperpolarized RMP and which mathematical description of IK1 is most suitable for DC experiments. We compared five mature IK1 formulations (Bett, Dhamoon, Ishihara, O'Hara-Rudy, and ten Tusscher) with the native one (Paci), evaluating the main properties (outward peak, degree of rectification), and we quantified their effects on AP features (RMP, V˙max, APD50, APD90 (AP duration at 50 and 90% of repolarization), and APD50/APD90) after including them in the hiPSC-CM mathematical model by Paci. Then, we automatically identified the critical conductance for IK1 ( GK1, critical), the minimally required amount of IK1 suppressing spontaneous activity. Preconditioning the cell model with depolarizing/hyperpolarizing prepulses allowed us to highlight time dependency of the IK1 formulations. Simulations showed that inclusion of mature IK1 formulations resulted in hyperpolarized RMP and higher V˙max, and observed GK1, critical and the effect on AP duration strongly depended on IK1 formulation. Finally, the Ishihara IK1 led to shorter (-16.3%) and prolonged (+6.5%) APD90 in response to hyperpolarizing and depolarizing prepulses, respectively, whereas other models showed negligible effects. Fine-tuning of GK1 is an important step in DC experiments. Our computational work proposes a procedure to automatically identify how much IK1 current is required to inject to stop the spontaneous activity and suggests the use of the Ishihara IK1 model to perform DC experiments in hiPSC-CMs.

Genetic variation in GNB5 causes bradycardia by augmenting the cholinergic response via increased acetylcholine-activated potassium current (IK,ACh).

Mutations in GNB5, encoding the G-protein ß5 subunit (Gß5), have recently been linked to a multisystem disorder that includes severe bradycardia. Here, we investigated the mechanism underlying bradycardia caused by the recessive p.S81L Gß5 variant. Using CRISPR/Cas9-based targeting, we generated an isogenic series of human induced pluripotent stem cell (hiPSC) lines that were either wild type, heterozygous or homozygous for the GNB5 p.S81L variant. These were differentiated into cardiomyocytes (hiPSC-CMs) that robustly expressed the acetylcholine-activated potassium channel [I(KACh); also known as IK,ACh]. Baseline electrophysiological properties of the lines did not differ. Upon application of carbachol (CCh), homozygous p.S81L hiPSC-CMs displayed an increased acetylcholine-activated potassium current (IK,ACh) density and a more pronounced decrease of spontaneous activity as compared to wild-type and heterozygous p.S81L hiPSC-CMs, explaining the bradycardia in homozygous carriers. Application of the specific I(KACh) blocker XEN-R0703 resulted in near-complete reversal of the phenotype. Our results provide mechanistic insights and proof of principle for potential therapy in patients carrying GNB5 mutations.This article has an associated First Person interview with the first author of the paper.

The influence of hERG1a and hERG1b isoforms on drug safety screening in iPSC-CMs.

The human Ether-à-go-go Related Gene (hERG) encodes the pore forming subunit of the channel that conducts the rapid delayed rectifier potassium current IKr. IKr drives repolarization in the heart and when IKr is dysfunctional, cardiac repolarization delays, the QT interval on the electrocardiogram (ECG) prolongs and the risk of developing lethal arrhythmias such as Torsade de Pointes (TdP) increases. TdP risk is incorporated in drug safety screening for cardiotoxicity where hERG is the main target since the IKr channels appear highly sensitive to blockage. hERG block is also included as an important read-out in the Comprehensive in Vitro Proarrhythmia Assay (CiPA) initiative which aims to combine in vitro and in silico experiments on induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) to screen for cardiotoxicity. However, the hERG channel has some unique features to consider for drug safety screening, which we will discuss in this study. The hERG channel consists of different isoforms, hERG1a and hERG1b, which individually influence the kinetics of the channel and the drug response in the human heart and in iPSC-CMs. hERG1b is often underappreciated in iPSC-CM studies, drug screening assays and in silico models, and the fact that its contribution might substantially differ between iPSC-CM and healthy but also diseased human heart, adds to this problem. In this study we show that the activation kinetics in iPSC-CMs resemble hERG1b kinetics using Cs+ as a charge carrier. Not including hERG1b in drug safety testing might underestimate the actual role of hERG1b in repolarization and drug response, and might lead to inappropriate conclusions. We stress to focus more on including hERG1b in drug safety testing concerning IKr.

Flotillins in the intercalated disc are potential modulators of cardiac excitability.

BACKGROUND: The intercalated disc (ID) is important for cardiac remodeling and has become a subject of intensive research efforts. However, as yet the composition of the ID has still not been conclusively resolved and the role of many proteins identified in the ID, like Flotillin-2, is often unknown. The Flotillin proteins are known to be involved in the stabilization of cadherins and desmosomes in the epidermis and upon cancer development. However, their role in the heart has so far not been investigated. Therefore, in this study, we aimed at identifying the role of Flotillin-1 and Flotillin-2 in the cardiac ID. METHODS: Location of Flotillins in human and murine cardiac tissue was evaluated by fluorescent immunolabeling and co-immunoprecipitation. In addition, the effect of Flotillin knockout (KO) on proteins of the ID and in electrical excitation and conduction was investigated in cardiac samples of wildtype (WT), Flotillin-1 KO, Flotilin-2 KO and Flotilin-1/2 double KO mice. Consequences of Flotillin knockdown (KD) on cardiac function were studied (patch clamp and Multi Electrode Array (MEA)) in neonatal rat cardiomyocytes (NRCMs) transfected with siRNAs against Flotillin-1 and/or Flotillin-2. RESULTS: First, we confirmed presence in the ID and mutual binding of Flotillin-1 and Flotillin-2 in murine and human cardiac tissue. Flotillin KO mice did not show cardiac fibrosis, nor hypertrophy or changes in expression of the desmosomal ID proteins. However, protein expression of the cardiac sodium channel NaV1.5 was significantly decreased in Flotillin-1 and Flotillin-1/2 KO mice compared to WT mice. In addition, sodium current density showed a significant decrease upon Flotillin-1/2 KD in NRCMs as compared to scrambled siRNA-transfected NRCMs. MEA recordings of Flotillin-2 KD NRCM cultures showed a significantly decreased spike amplitude and a tendency of a reduced spike slope when compared to control and scrambled siRNA-transfected cultures. CONCLUSIONS: In this study, we demonstrate the presence of Flotillin-1, in addition to Flotillin-2 in the cardiac ID. Our findings indicate a modulatory role of Flotillins on NaV1.5 expression at the ID, with potential consequences for cardiac excitation.

Optogenetic sensors in the zebrafish heart: a novel in vivo electrophysiological tool to study cardiac arrhythmogenesis.

Cardiac arrhythmias are among the most challenging human disorders to diagnose and treat due to their complex underlying pathophysiology. Suitable experimental animal models are needed to study the mechanisms causative for cardiac arrhythmogenesis. To enable in vivo analysis of cardiac cellular electrophysiology with a high spatial and temporal resolution, we generated and carefully validated two zebrafish models, one expressing an optogenetic voltage indicator (chimeric VSFP-butterfly CY) and the other a genetically encoded calcium indicator (GCaMP6f) in the heart. Methods: High-speed epifluorescence microscopy was used to image chimeric VSFP-butterfly CY and GCaMP6f in the embryonic zebrafish heart, providing information about the spatiotemporal patterning of electrical activation, action potential configuration and intracellular Ca2+ dynamics. Plotting VSFP or GCaMP6f signals on a line along the myocardial wall over time facilitated the visualization and analysis of electrical impulse propagation throughout the heart. Administration of drugs targeting the sympathetic nervous system or cardiac ion channels was used to validate sensitivity and kinetics of both zebrafish sensor lines. Using the same microscope setup, we imaged transparent juvenile casper fish expressing GCaMP6f, demonstrating the feasibility of imaging cardiac optogenetic sensors at later stages of development. Results: Isoproterenol slightly increased heart rate, diastolic Ca2+ levels and Ca2+ transient amplitudes, whereas propranolol caused a profound decrease in heart rate and Ca2+ transient parameters in VSFP-Butterfly and GCaMP6f embryonic fish. Ikr blocker E-4031 decreased heart rate and increased action potential duration in VSFP-Butterfly fish. ICa,L blocker nifedipine caused total blockade of Ca2+ transients in GCaMP6f fish and a reduced heart rate, altered ventricular action potential duration and disrupted atrial-ventricular electrical conduction in VSFP-Butterfly fish. Imaging of juvenile animals demonstrated the possibility of employing an older zebrafish model for in vivo cardiac electrophysiology studies. We observed differences in atrial and ventricular Ca2+ recovery dynamics between 3 dpf and 14 dpf casper fish, but not in Ca2+ upstroke dynamics. Conclusion: By introducing the optogenetic sensors chimeric VSFP-butterfly CY and GCaMP6f into the zebrafish we successfully generated an in vivo cellular electrophysiological readout tool for the zebrafish heart. Complementary use of both sensor lines demonstrated the ability to study heart rate, cardiac action potential configuration, spatiotemporal patterning of electrical activation and intracellular Ca2+ homeostasis in embryonic zebrafish. In addition, we demonstrated the first successful use of an optogenetic sensor to study cardiac function in older zebrafish. These models present a promising new research tool to study the underlying mechanisms of cardiac arrhythmogenesis.

Cardiac Ca2+ signalling in zebrafish: Translation of findings to man.

Sudden cardiac death is a leading cause of death worldwide, mainly caused by highly disturbed electrical activation patterns in the heart. Currently, murine models are the most popular model to study underlying molecular mechanisms of inherited or acquired cardiac electrical abnormalities, although the numerous electrophysiological discrepancies between mouse and human raise the question whether mice are the optimal model to study cardiac rhythm disorders. Recently it has been uncovered that the zebrafish cardiac electrophysiology seems surprisingly similar to the human heart, mainly because the zebrafish AP contains a clear plateau phase and ECG characteristics show alignment with the human ECG. Although, before using zebrafish as a model to study cardiac arrhythmogenesis, however, it is very important to gain a better insight into the electrophysiological characteristics of the zebrafish heart. In this review we outline the electrophysiological machinery of the zebrafish cardiomyocytes, with a special focus on the intracellular Ca2+ dynamics and excitation-contraction coupling. We debate the potential of zebrafish as a model to study human cardiovascular diseases and postulate steps to employ zebrafish into a more 'humanized' model.

Sinusoidal voltage protocols for rapid characterisation of ion channel kinetics.

KEY POINTS: Ion current kinetics are commonly represented by current-voltage relationships, time constant-voltage relationships and subsequently mathematical models fitted to these. These experiments take substantial time, which means they are rarely performed in the same cell. Rather than traditional square-wave voltage clamps, we fitted a model to the current evoked by a novel sum-of-sinusoids voltage clamp that was only 8 s long. Short protocols that can be performed multiple times within a single cell will offer many new opportunities to measure how ion current kinetics are affected by changing conditions. The new model predicts the current under traditional square-wave protocols well, with better predictions of underlying currents than literature models. The current under a novel physiologically relevant series of action potential clamps is predicted extremely well. The short sinusoidal protocols allow a model to be fully fitted to individual cells, allowing us to examine cell-cell variability in current kinetics for the first time. ABSTRACT: Understanding the roles of ion currents is crucial to predict the action of pharmaceuticals and mutations in different scenarios, and thereby to guide clinical interventions in the heart, brain and other electrophysiological systems. Our ability to predict how ion currents contribute to cellular electrophysiology is in turn critically dependent on our characterisation of ion channel kinetics - the voltage-dependent rates of transition between open, closed and inactivated channel states. We present a new method for rapidly exploring and characterising ion channel kinetics, applying it to the hERG potassium channel as an example, with the aim of generating a quantitatively predictive representation of the ion current. We fitted a mathematical model to currents evoked by a novel 8 second sinusoidal voltage clamp in CHO cells overexpressing hERG1a. The model was then used to predict over 5 minutes of recordings in the same cell in response to further protocols: a series of traditional square step voltage clamps, and also a novel voltage clamp comprising a collection of physiologically relevant action potentials. We demonstrate that we can make predictive cell-specific models that outperform the use of averaged data from a number of different cells, and thereby examine which changes in gating are responsible for cell-cell variability in current kinetics. Our technique allows rapid collection of consistent and high quality data, from single cells, and produces more predictive mathematical ion channel models than traditional approaches.