A massively parallel assay accurately discriminates between functionally normal and abnormal variants in a hotspot domain of KCNH2.

Many genes, including KCNH2, contain "hotspot" domains associated with a high density of variants associated with disease. This has led to the suggestion that variant location can be used as evidence supporting classification of clinical variants. However, it is not known what proportion of all potential variants in hotspot domains cause loss of function. Here, we have used a massively parallel trafficking assay to characterize all single-nucleotide variants in exon 2 of KCNH2, a known hotspot for variants that cause long QT syndrome type 2 and an increased risk of sudden cardiac death. Forty-two percent of KCNH2 exon 2 variants caused at least 50% reduction in protein trafficking, and 65% of these trafficking-defective variants exerted a dominant-negative effect when co-expressed with a WT KCNH2 allele as assessed using a calibrated patch-clamp electrophysiology assay. The massively parallel trafficking assay was more accurate (AUC of 0.94) than bioinformatic prediction tools (REVEL and CardioBoost, AUC of 0.81) in discriminating between functionally normal and abnormal variants. Interestingly, over half of variants in exon 2 were found to be functionally normal, suggesting a nuanced interpretation of variants in this "hotspot" domain is necessary. Our massively parallel trafficking assay can provide this information prospectively.

A calibrated functional patch-clamp assay to enhance clinical variant interpretation in KCNH2-related long QT syndrome.

Modern sequencing technologies have revolutionized our detection of gene variants. However, in most genes, including KCNH2, the majority of missense variants are currently classified as variants of uncertain significance (VUSs). The aim of this study was to investigate the utility of an automated patch-clamp assay for aiding clinical variant classification in KCNH2. The assay was designed according to recommendations proposed by the Clinical Genome Sequence Variant Interpretation Working Group. Thirty-one variants (17 pathogenic/likely pathogenic, 14 benign/likely benign) were classified internally as variant controls. They were heterozygously expressed in Flp-In HEK293 cells for assessing the effects of variants on current density and channel gating in order to determine the sensitivity and specificity of the assay. All 17 pathogenic variant controls had reduced current density, and 13 of 14 benign variant controls had normal current density, which enabled determination of normal and abnormal ranges for applying evidence of moderate or supporting strength for VUS reclassification. Inclusion of functional assay evidence enabled us to reclassify 6 out of 44 KCNH2 VUSs as likely pathogenic. The high-throughput patch-clamp assay can provide moderate-strength evidence for clinical interpretation of clinical KCNH2 variants and demonstrates the value of developing automated patch-clamp assays for functional characterization of ion channel gene variants.

Domain reorientation and rotation of an intracellular assembly regulate conduction in Kir potassium channels.

Potassium channels embedded in cell membranes employ gates to regulate K+ current. While a specific constriction in the permeation pathway has historically been implicated in gating, recent reports suggest that the signature ion selectivity filter located in the outer membrane leaflet may be equally important. Inwardly rectifying K+ channels also control the directionality of flow, using intracellular polyamines to stem ion efflux by a valve-like action. This study presents crystallographic evidence of interdependent gates in the conduction pathway and reveals the mechanism of polyamine block. Reorientation of the intracellular domains, concomitant with activation, instigates polyamine release from intracellular binding sites to block the permeation pathway. Conformational adjustments of the slide helices, achieved by rotation of the cytoplasmic assembly relative to the pore, are directly correlated to the ion configuration in the selectivity filter. Ion redistribution occurs irrespective of the constriction, suggesting a more expansive role of the selectivity filter in gating than previously appreciated.

Multitarget nociceptor sensitization by a promiscuous peptide from the venom of the King Baboon spider.

The King Baboon spider, Pelinobius muticus, is a burrowing African tarantula. Its impressive size and appealing coloration are tempered by reports describing severe localized pain, swelling, itchiness, and muscle cramping after accidental envenomation. Hyperalgesia is the most prominent symptom after bites from P. muticus, but the molecular basis by which the venom induces pain is unknown. Proteotranscriptomic analysis of P. muticus venom uncovered a cysteine-rich peptide, δ/κ-theraphotoxin-Pm1a (δ/κ-TRTX-Pm1a), that elicited nocifensive behavior when injected into mice. In small dorsal root ganglion neurons, synthetic δ/κ-TRTX-Pm1a (sPm1a) induced hyperexcitability by enhancing tetrodotoxin-resistant sodium currents, impairing repolarization and lowering the threshold of action potential firing, consistent with the severe pain associated with envenomation. The molecular mechanism of nociceptor sensitization by sPm1a involves multimodal actions over several ion channel targets, including NaV1.8, KV2.1, and tetrodotoxin-sensitive NaV channels. The promiscuous targeting of peptides like δ/κ-TRTX-Pm1a may be an evolutionary adaptation in pain-inducing defensive venoms.

Functional coupling between Piezo1 and TRPM4 influences the electrical activity of HL-1 atrial myocytes.

The transient receptor potential melastatin 4 (TRPM4) channel contributes extensively to cardiac electrical activity, especially cardiomyocyte action potential formation. Mechanical stretch can induce changes in heart rate and rhythm, and the mechanosensitive channel Piezo1 is expressed in many cell types within the myocardium. Our previous study showed that TRPM4 and Piezo1 are closely co-localized in the t-tubules of ventricular cardiomyocytes and contribute to the Ca2+ -dependent signalling cascade that underlies hypertrophy in response to mechanical pressure overload. However, there was no direct evidence showing that Piezo1 activation was related to TRPM4 activation in situ. In the present study, we employed the HL-1 mouse atrial myocyte-like cell line as an in vitro model to investigate whether Piezo1-TRPM4 coupling can affect action potential properties. We used the small molecule Piezo1 agonist, Yoda1, as a surrogate for mechanical stretch to activate Piezo1 and detected the action potential changes in HL-1 cells using FluoVolt, a fluorescent voltage sensitive dye. Our results demonstrate that Yoda1-induced activation of Piezo1 changes the action potential frequency in HL-1 cells. This change in action potential frequency is reduced by Piezo1 knockdown using small intefering RNA. Importantly knockdown or pharmacological inhibition of TRPM4 significantly affected the degree to which Yoda1-evoked Piezo1 activation influenced action potential frequency. Thus, the present study provides in vitro evidence of a functional coupling between Piezo1 and TRPM4 in a cardiomyocyte-like cell line. The coupling of a mechanosensitive Ca2+ permeable channel and a Ca2+ -activated TRP channel probably represents a ubiquitous model for the role of TRP channels in mechanosensory transduction. KEY POINTS: The transient receptor potential melastatin 4 (TRPM4) and Piezo1 channels have been confirmed to contribute to the Ca2+ -dependent signalling cascade that underlies cardiac hypertrophy in response to mechanical pressure overload. However, there was no direct evidence showing that Piezo1 activation was related to TRPM4 activation in situ. We employed the HL-1 mouse atrial myocyte-like cell line as an in vitro model to investigate the effect of Piezo1-TRPM4 coupling on cardiac electrical properties. The results show that both pharmacological and genetic inhibition of TRPM4 significantly affected the degree to which Piezo1 activation influenced action potential frequency in HL-1 cells. Our findings provide in vitro evidence of a functional coupling between Piezo1 and TRPM4 in a cardiomyocyte-like cell line. The coupling of a mechanosensitive Ca2+ permeable channel and a Ca2+ -activated TRP channel probably represents a ubiquitous model for the role of TRP channels in mechanosensory transduction in various (patho)physiological processes.

Empirical Quantification of Predictive Uncertainty Due to Model Discrepancy by Training with an Ensemble of Experimental Designs: An Application to Ion Channel Kinetics.

When using mathematical models to make quantitative predictions for clinical or industrial use, it is important that predictions come with a reliable estimate of their accuracy (uncertainty quantification). Because models of complex biological systems are always large simplifications, model discrepancy arises-models fail to perfectly recapitulate the true data generating process. This presents a particular challenge for making accurate predictions, and especially for accurately quantifying uncertainty in these predictions. Experimentalists and modellers must choose which experimental procedures (protocols) are used to produce data used to train models. We propose to characterise uncertainty owing to model discrepancy with an ensemble of parameter sets, each of which results from training to data from a different protocol. The variability in predictions from this ensemble provides an empirical estimate of predictive uncertainty owing to model discrepancy, even for unseen protocols. We use the example of electrophysiology experiments that investigate the properties of hERG potassium channels. Here, 'information-rich' protocols allow mathematical models to be trained using numerous short experiments performed on the same cell. In this case, we simulate data with one model and fit it with a different (discrepant) one. For any individual experimental protocol, parameter estimates vary little under repeated samples from the assumed additive independent Gaussian noise model. Yet parameter sets arising from the same model applied to different experiments conflict-highlighting model discrepancy. Our methods will help select more suitable ion channel models for future studies, and will be widely applicable to a range of biological modelling problems.

The diagnostic role of T wave morphology biomarkers in congenital and acquired long QT syndrome: A systematic review.

INTRODUCTION: QTc prolongation is key in diagnosing long QT syndrome (LQTS), however 25%-50% with congenital LQTS (cLQTS) demonstrate a normal resting QTc. T wave morphology (TWM) can distinguish cLQTS subtypes but its role in acquired LQTS (aLQTS) is unclear. METHODS: Electronic databases were searched using the terms "LQTS," "long QT syndrome," "QTc prolongation," "prolonged QT," and "T wave," "T wave morphology," "T wave pattern," "T wave biomarkers." Whole text articles assessing TWM, independent of QTc, were included. RESULTS: Seventeen studies met criteria. TWM measurements included T-wave amplitude, duration, magnitude, Tpeak-Tend, QTpeak, left and right slope, center of gravity (COG), sigmoidal and polynomial classifiers, repolarizing integral, morphology combination score (MCS) and principal component analysis (PCA); and vectorcardiographic biomarkers. cLQTS were distinguished from controls by sigmoidal and polynomial classifiers, MCS, QTpeak, Tpeak-Tend, left slope; and COG x axis. MCS detected aLQTS more significantly than QTc. Flatness, asymmetry and notching, J-Tpeak; and Tpeak-Tend correlated with QTc in aLQTS. Multichannel block in aLQTS was identified by early repolarization (ERD30% ) and late repolarization (LRD30% ), with ERD reflecting hERG-specific blockade. Cardiac events were predicted in cLQTS by T wave flatness, notching, and inversion in leads II and V5 , left slope in lead V6 ; and COG last 25% in lead I. T wave right slope in lead I and T-roundness achieved this in aLQTS. CONCLUSION: Numerous TWM biomarkers which supplement QTc assessment were identified. Their diagnostic capabilities include differentiation of genotypes, identification of concealed LQTS, differentiating aLQTS from cLQTS; and determining multichannel versus hERG channel blockade.

Arrhythmogenic effects of ultra-long and bistable cardiac action potentials.

Contemporary accounts of the initiation of cardiac arrhythmias typically rely on after-depolarizations as the trigger for reentrant activity. The after-depolarizations are usually triggered by calcium entry or spontaneous release within the cells of the myocardium or the conduction system. Here we propose an alternative mechanism whereby arrhythmias are triggered autonomously by cardiac cells that fail to repolarize after a normal heartbeat. We investigated the proposal by representing the heart as an excitable medium of FitzHugh-Nagumo cells where a proportion of cells were capable of remaining depolarized indefinitely. As such, those cells exhibit bistable membrane dynamics. We found that heterogeneous media can tolerate a surprisingly large number of bistable cells and still support normal rhythmic activity. Yet there is a critical limit beyond which the medium is persistently arrhythmogenic. Numerical analysis revealed that the critical threshold for arrhythmogenesis depends on both the strength of the coupling between cells and the extent to which the abnormal cells resist repolarization. Moreover, arrhythmogenesis was found to emerge preferentially at tissue boundaries where cells naturally have fewer neighbors to influence their behavior. These findings may explain why atrial fibrillation typically originates from tissue boundaries such as the cuff of the pulmonary vein.

Protocol-Dependent Differences in IC50 Values Measured in Human Ether-Á-Go-Go-Related Gene Assays Occur in a Predictable Way and Can Be Used to Quantify State Preference of Drug Binding.

Current guidelines around preclinical screening for drug-induced arrhythmias require the measurement of the potency of block of voltage-gated potassium channel subtype 11.1 (Kv11.1) as a surrogate for risk. A shortcoming of this approach is that the measured IC50 of Kv11.1 block varies widely depending on the voltage protocol used in electrophysiological assays. In this study, we aimed to investigate the factors that contribute to these differences and to identify whether it is possible to make predictions about protocol-dependent block that might facilitate the comparison of potencies measured using different assays. Our data demonstrate that state preferential binding, together with drug-binding kinetics and trapping, is an important determinant of the protocol dependence of Kv11.1 block. We show for the first time that differences in IC50 measured between protocols occurs in a predictable way, such that machine-learning algorithms trained using a selection of simple voltage protocols can indeed predict protocol-dependent potency. Furthermore, we also show that the preference of a drug for binding to the open versus the inactivated state of Kv11.1 can also be inferred from differences in IC50 values measured between protocols. Our work therefore identifies how state preferential drug binding is a major determinant of the protocol dependence of IC50 values measured in preclinical Kv11.1 assays. It also provides a novel method for quantifying the state dependence of Kv11.1 drug binding that will facilitate the development of more complete models of drug binding to Kv11.1 and improve our understanding of proarrhythmic risk associated with compounds that block Kv11.1.

hERG K(+) channels: structure, function, and clinical significance.

The human ether-a-go-go related gene (hERG) encodes the pore-forming subunit of the rapid component of the delayed rectifier K(+) channel, Kv11.1, which are expressed in the heart, various brain regions, smooth muscle cells, endocrine cells, and a wide range of tumor cell lines. However, it is the role that Kv11.1 channels play in the heart that has been best characterized, for two main reasons. First, it is the gene product involved in chromosome 7-associated long QT syndrome (LQTS), an inherited disorder associated with a markedly increased risk of ventricular arrhythmias and sudden cardiac death. Second, blockade of Kv11.1, by a wide range of prescription medications, causes drug-induced QT prolongation with an increase in risk of sudden cardiac arrest. In the first part of this review, the properties of Kv11.1 channels, including biogenesis, trafficking, gating, and pharmacology are discussed, while the second part focuses on the pathophysiology of Kv11.1 channels.

The Temperature Dependence of Kinetics Associated with Drug Block of hERG Channels Is Compound-Specific and an Important Factor for Proarrhythmic Risk Prediction.

Current mandated preclinical tests for drug-induced proarrhythmia are very sensitive, but not sufficiently specific. This has led to concern that there is a high attrition rate of potentially safe drugs that could have been beneficial to patients. The comprehensive in vitro proarrhythmia initiative has proposed new metrics based around in silico risk predictions, which are informed, among other things, by measures of human ether-à-go-go-related gene channel (hERG) block kinetics. However, high-throughput patch-clamp systems set to collect these data largely operate at ambient temperature, whereas the simulations for risk prediction are carried out at physiologic temperature. The aims of this study were to: 1) determine to what degree kinetics of drug block of hERG are temperature-dependent, 2) assess the impact of any temperature dependence of drug binding kinetics on repolarization in silico, and 3) identify whether a common set of Q10 scalars can be used to extrapolate kinetic data gathered at ambient to physiologic temperatures for use in in silico proarrhythmic risk prediction. We show that, for a range of drugs, kinetics of block are temperature-dependent and, furthermore, that the degree of temperature dependence is different for each drug. As a result, no common set of Q10 scalars could describe the observed range of temperature dependencies. These results suggest that if accurate physiologic temperature models of the kinetics of drug binding are important for in silico risk prediction, the in vitro data should be acquired at physiologic temperature.

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.

Convergence of models of human ventricular myocyte electrophysiology after global optimization to recapitulate clinical long QT phenotypes.

In-silico models of human cardiac electrophysiology are now being considered for prediction of cardiotoxicity as part of the preclinical assessment phase of all new drugs. We ask the question whether any of the available models are actually fit for this purpose. We tested three models of the human ventricular action potential, the O'hara-Rudy (ORD11), the Grandi-Bers (GB10) and the Ten Tusscher (TT06) models. We extracted clinical QT data for LQTS1 and LQTS2 patients with nonsense mutations that would be predicted to cause 50% loss of function in IKs and IKr respectively. We also obtained clinical QT data for LQTS3 patients. We then used a global optimization approach to improve the existing in silico models so that they reproduced all three clinical data sets more closely. We also examined the effects of adrenergic stimulation in the different LQTS subsets. All models, in their original form, produce markedly different and unrealistic predictions of QT prolongation for LQTS1, 2 and 3. After global optimization of the maximum conductances for membrane channels, all models have similar current densities during the action potential, despite differences in kinetic properties of the channels in the different models, and more closely reproduce the prolongation of repolarization seen in all LQTS subtypes. In-silico models of cardiac electrophysiology have the potential to be tremendously useful in complementing traditional preclinical drug testing studies. However, our results demonstrate they should be carefully validated and optimized to clinical data before they can be used for this purpose.

Temperature Effects on Kinetics of KV11.1 Drug Block Have Important Consequences for In Silico Proarrhythmic Risk Prediction.

Drug block of voltage-gated potassium channel subtype 11.1 human ether-a-go-go related gene (Kv11.1) (hERG) channels, encoded by the KCNH2 gene, is associated with reduced repolarization of the cardiac action potential and is the predominant cause of acquired long QT syndrome that can lead to fatal cardiac arrhythmias. Current safety guidelines require that potency of KV11.1 block is assessed in the preclinical phase of drug development. However, not all drugs that block KV11.1 are proarrhythmic, meaning that screening on the basis of equilibrium measures of block can result in high attrition of potentially low-risk drugs. The basis of the next generation of drug-screening approaches is set to be in silico risk prediction, informed by in vitro mechanistic descriptions of drug binding, including measures of the kinetics of block. A critical issue in this regard is characterizing the temperature dependence of drug binding. Specifically, it is important to address whether kinetics relevant to physiologic temperatures can be inferred or extrapolated from in vitro data gathered at room temperature in high-throughout systems. Here we present the first complete study of the temperature-dependent kinetics of block and unblock of a proarrhythmic drug, cisapride, to KV11.1. Our data highlight a complexity to binding that manifests at higher temperatures and can be explained by accumulation of an intermediate, non-blocking encounter-complex. These results suggest that for cisapride, physiologically relevant kinetic parameters cannot be simply extrapolated from those measured at lower temperatures; rather, data gathered at physiologic temperatures should be used to constrain in silico models that may be used for proarrhythmic risk prediction.

Computational cardiology and risk stratification for sudden cardiac death: one of the grand challenges for cardiology in the 21st century.

Risk stratification in the context of sudden cardiac death has been acknowledged as one of the major challenges facing cardiology for the past four decades. In recent years, the advent of high performance computing has facilitated organ-level simulation of the heart, meaning we can now examine the causes, mechanisms and impact of cardiac dysfunction in silico. As a result, computational cardiology, largely driven by the Physiome project, now stands at the threshold of clinical utility in regards to risk stratification and treatment of patients at risk of sudden cardiac death. In this white paper, we outline a roadmap of what needs to be done to make this translational step, using the relatively well-developed case of acquired or drug-induced long QT syndrome as an exemplar case.

Rescue of protein expression defects may not be enough to abolish the pro-arrhythmic phenotype of long QT type 2 mutations.

KEY POINTS: Most missense long QT syndrome type 2 (LQTS2) mutations result in Kv11.1 channels that show reduced levels of membrane expression. Pharmacological chaperones that rescue mutant channel expression could have therapeutic potential to reduce the risk of LQTS2-associated arrhythmias and sudden cardiac death, but only if the mutant Kv11.1 channels function normally (i.e. like WT channels) after membrane expression is restored. Fewer than half of mutant channels exhibit relatively normal function after rescue by low temperature. The remaining rescued missense mutant Kv11.1 channels have perturbed gating and/or ion selectivity characteristics. Co-expression of WT subunits with gating defective missense mutations ameliorates but does not eliminate the functional abnormalities observed for most mutant channels. For patients with mutations that affect gating in addition to expression, it may be necessary to use a combination therapy to restore both normal function and normal expression of the channel protein. ABSTRACT: In the heart, Kv11.1 channels pass the rapid delayed rectifier current (IKr ) which plays critical roles in repolarization of the cardiac action potential and in the suppression of arrhythmias caused by premature stimuli. Over 500 inherited mutations in Kv11.1 are known to cause long QT syndrome type 2 (LQTS2), a cardiac electrical disorder associated with an increased risk of life threatening arrhythmias. Most missense mutations in Kv11.1 reduce the amount of channel protein expressed at the membrane and, as a consequence, there has been considerable interest in developing pharmacological agents to rescue the expression of these channels. However, pharmacological chaperones will only have clinical utility if the mutant Kv11.1 channels function normally after membrane expression is restored. The aim of this study was to characterize the gating phenotype for a subset of LQTS2 mutations to assess what proportion of mutations may be suitable for rescue. As an initial screen we used reduced temperature to rescue expression defects of mutant channels expressed in Xenopus laevis oocytes. Over half (∼56%) of Kv11.1 mutants exhibited functional gating defects that either dramatically reduced the amount of current contributing to cardiac action potential repolarization and/or reduced the amount of protective current elicited in response to premature depolarizations. Our data demonstrate that if pharmacological rescue of protein expression defects is going to have clinical utility in the treatment of LQTS2 then it will be important to assess the gating phenotype of LQTS2 mutations before attempting rescue.

Getting to the heart of hERG K(+) channel gating.

Potassium ion channels encoded by the human ether-a-go-go related gene (hERG) form the ion-conducting subunit of the rapid delayed rectifier potassium current (IKr ). Although hERG channels exhibit a widespread tissue distribution they play a particularly important role in the heart. There has been considerable interest in hERG K(+) channels for three main reasons. First, they have very unusual gating kinetics, most notably rapid and voltage-dependent inactivation coupled to slow deactivation, which has led to the suggestion that they may play a specific role in the suppression of arrhythmias. Second, mutations in hERG are the cause of 30-40% of cases of congenital long QT syndrome (LQTS), the commonest inherited primary arrhythmia syndrome. Third, hERG is the molecular target for the vast majority of drugs that cause drug-induced LQTS, the commonest cause of drug-induced arrhythmias and cardiac death. Drug-induced LQTS has now been reported for a large range of both cardiac and non-cardiac drugs, in which this side effect is entirely undesired. In recent years there have been comprehensive reviews published on hERG K(+) channels (Vandenberg et al. 2012) and we will not re-cover this ground. Rather, we focus on more recent work on the structural basis and dynamics of hERG gating with an emphasis on how the latest developments may facilitate translational research in the area of stratifying risk of arrhythmias.

Multiple interactions between cytoplasmic domains regulate slow deactivation of Kv11.1 channels.

The intracellular domains of many ion channels are important for fine-tuning their gating kinetics. In Kv11.1 channels, the slow kinetics of channel deactivation, which are critical for their function in the heart, are largely regulated by the N-terminal N-Cap and Per-Arnt-Sim (PAS) domains, as well as the C-terminal cyclic nucleotide-binding homology (cNBH) domain. Here, we use mutant cycle analysis to probe for functional interactions between the N-Cap/PAS domains and the cNBH domain. We identified a specific and stable charge-charge interaction between Arg(56) of the PAS domain and Asp(803) of the cNBH domain, as well an additional interaction between the cNBH domain and the N-Cap, both of which are critical for maintaining slow deactivation kinetics. Furthermore, we found that positively charged arginine residues within the disordered region of the N-Cap interact with negatively charged residues of the C-linker domain. Although this interaction is likely more transient than the PAS-cNBD interaction, it is strong enough to stabilize the open conformation of the channel and thus slow deactivation. These findings provide novel insights into the slow deactivation mechanism of Kv11.1 channels.

PhosphOrtholog: a web-based tool for cross-species mapping of orthologous protein post-translational modifications.

BACKGROUND: Most biological processes are influenced by protein post-translational modifications (PTMs). Identifying novel PTM sites in different organisms, including humans and model organisms, has expedited our understanding of key signal transduction mechanisms. However, with increasing availability of deep, quantitative datasets in diverse species, there is a growing need for tools to facilitate cross-species comparison of PTM data. This is particularly important because functionally important modification sites are more likely to be evolutionarily conserved; yet cross-species comparison of PTMs is difficult since they often lie in structurally disordered protein domains. Current tools that address this can only map known PTMs between species based on known orthologous phosphosites, and do not enable the cross-species mapping of newly identified modification sites. Here, we addressed this by developing a web-based software tool, PhosphOrtholog ( www.phosphortholog.com ) that accurately maps protein modification sites between different species. This facilitates the comparison of datasets derived from multiple species, and should be a valuable tool for the proteomics community. RESULTS: Here we describe PhosphOrtholog, a web-based application for mapping known and novel orthologous PTM sites from experimental data obtained from different species. PhosphOrtholog is the only generic and automated tool that enables cross-species comparison of large-scale PTM datasets without relying on existing PTM databases. This is achieved through pairwise sequence alignment of orthologous protein residues. To demonstrate its utility we apply it to two sets of human and rat muscle phosphoproteomes generated following insulin and exercise stimulation, respectively, and one publicly available mouse phosphoproteome following cellular stress revealing high mapping and coverage efficiency. Although coverage statistics are dataset dependent, PhosphOrtholog increased the number of cross-species mapped sites in all our example data sets by more than double when compared to those recovered using existing resources such as PhosphoSitePlus. CONCLUSIONS: PhosphOrtholog is the first tool that enables mapping of thousands of novel and known protein phosphorylation sites across species, accessible through an easy-to-use web interface. Identification of conserved PTMs across species from large-scale experimental data increases our knowledgebase of functional PTM sites. Moreover, PhosphOrtholog is generic being applicable to other PTM datasets such as acetylation, ubiquitination and methylation.

Kinetics of drug interaction with the Kv11.1 potassium channel.

The Kv11.1 potassium channel is the molecular target for the majority of drugs implicated in acquired long QT syndrome, the most common cause of drug-induced sudden cardiac death, and a common reason for drug restriction or withdrawal from the market. While the IC50 for block of Kv11.1 is commonly used to estimate the risk of acquired long QT syndrome, this approach is crude, and it is widely accepted that the kinetics of drug interactions with the channel are a critical component in understanding their mechanism of action and risk profiles. In this study we report the first directly measured kinetics of block and unblock of Kv11.1 by a QT prolonging drug: the antipsychotic clozapine. Our data show that clozapine binding to Kv11.1 is complex. There are at least two kinetically distinct components to both block and unblock, while the kinetics of unblock are dependent on the dose or duration of drug application. Based on these observations, we have proposed a model incorporating kinetically distinct binding to the open and inactivated states of Kv11.1 that can describe the observed kinetic features of clozapine block and correctly predict the overall affinity and apparent nonstate-dependent interaction of clozapine with Kv11.1. Mechanistic insights into drug block of Kv11.1 gained though detailed kinetic analyses such as this have a potential role in development of drugs targeted to specific channel states to reduce unwanted side effects, as well as in the design of better high-throughput preclinical tests for assessing the proarrhythmic effects of QT prolonging drugs.