Search

The canonical mechanistic model explaining potassium channel gating is of a conformational change that alternately dilates and constricts a collar-like intracellular entrance to the pore. It is based on the premise that K+ ions maintain a complete hydration shell while passing between the transmembrane cavity and cytosol, which must be accommodated. To put the canonical model to the test, we locked the conformation of a Kir K+ channel to prevent widening of the narrow collar. Unexpectedly, conduction was unimpaired in the locked channels. In parallel, we employed all-atom molecular dynamics to simulate K+ ions moving along the conduction pathway between the lower cavity and cytosol. During simulations, the constriction did not significantly widen. Instead, transient loss of some water molecules facilitated K+ permeation through the collar. The low free energy barrier to partial dehydration in the absence of conformational change indicates Kir channels are not gated by the canonical mechanism., The transition between conducting and non-conducting states of K+ channels has been explained by conformational changes at the intracellular entrance to the conduction pathway. Here authors demonstrate that control over K+ currents in Kir channels is not explained by the canonical pore-gating model, as conduction is not impaired by a constricted inner helix bundle.

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.

KCNH2 is one of the 59 medically actionable genes recommended by the American College of Medical Genetics for reporting of incidental findings from clinical genomic sequencing. However, half of the reported KCNH2 variants in the ClinVar database are classified as variants of uncertain significance. In the absence of strong clinical phenotypes, there is a need for functional phenotyping to help decipher the significance of variants identified incidentally. Here, we report detailed methods for assessing the molecular phenotype of any KCNH2 missense variant. The key components of the assay include quick and cost-effective generation of a bi-cistronic vector to co-express Wild-type (WT) and any KCNH2 variant allele, generation of stable Flp-In HEK293 cell lines and high-throughput automated patch clamp electrophysiology analysis of channel function. Stable cell lines take 3–4 weeks to produce and can be generated in bulk, which will then allow up to 30 variants to be phenotyped per week after 48 h of channel expression. This high-throughput functional genomics assay will enable a much more rapid assessment of the extent of loss of function of any KCNH2 variant.

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-a-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.

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 Ca2+ containing external solution — so that a CaF2 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 CaF2 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.

The hERG (human ether-a-go-go-related gene) encoded potassium ion (K+) channel plays a major role in cardiac repolarization. Drug-induced blockade of hERG has been a major cause of potentially lethal ventricular tachycardia termed Torsades de Pointes (TdPs). Therefore, we presented a pharmacoinformatics strategy using combined ligand and structure based models for the prediction of hERG inhibition potential (IC50) of new chemical entities (NCEs) during early stages of drug design and development. Integrated GRid-INdependent Descriptor (GRIND) models, and lipophilic efficiency (LipE), ligand efficiency (LE) guided template selection for the structure based pharmacophore models have been used for virtual screening and subsequent hERG activity (pIC50) prediction of identified hits. Finally selected two hits were experimentally evaluated for hERG inhibition potential (pIC50) using whole cell patch clamp assay. Overall, our results demonstrate a difference of less than ±1.6 log unit between experimentally determined and predicted hERG inhibition potential (IC50) of the selected hits. This revealed predictive ability and robustness of our models and could help in correctly rank the potency order (lower μM to higher nM range) against hERG.

Spider venoms are a rich source of insecticidal peptide toxins. Their development as bioinsecticides has, however, been hampered due to concerns about potential lack of stability and oral bioactivity. We therefore systematically evaluated several synthetic strategies to increase the stability and oral potency of the potent insecticidal spider-venom peptide &omega, HXTX-Hv1a (Hv1a). Selective chemical replacement of disulfide bridges with diselenide bonds and N- to C-terminal cyclization were anticipated to improve Hv1a resistance to proteolytic digestion, and thereby its activity when delivered orally. We found that native Hv1a is orally active in blowflies, but 91-fold less potent than when administered by injection. Introduction of a single diselenide bond had no effect on the susceptibility to scrambling or the oral activity of Hv1a. N- to C-terminal cyclization of the peptide backbone did not significantly improve the potency of Hv1a when injected into blowflies and it led to a significant decrease in oral activity. We show that this is likely due to a dramatically reduced rate of translocation of cyclic Hv1a across the insect midgut, highlighting the importance of testing bioavailability in addition to toxin stability.

Acquired long QT syndrome, mostly as a result of drug block of the Kv11. 1 potassium channel in the heart, is characterized by delayed cardiac myocyte repolarization, prolongation of the T interval on the ECG, syncope and sudden cardiac death due to the polymorphic ventricular arrhythmia Torsade de Pointes (TdP). In recent years, efforts are underway through the Comprehensive in vitro proarrhythmic assay (CiPA) initiative, to develop better tests for this drug induced arrhythmia based in part on in silico simulations of pharmacological disruption of repolarization. However, drug binding to Kv11.1 is more complex than a simple binary molecular reaction, meaning simple steady state measures of potency are poor surrogates for risk. As a result, there is a plethora of mechanistic detail describing the drug/Kv11.1 interaction-such as drug binding kinetics, state preference, temperature dependence and trapping-that needs to be considered when developing in silico models for risk prediction. In addition to this, other factors, such as multichannel pharmacological profile and the nature of the ventricular cell models used in simulations also need to be considered in the search for the optimum in silico approach. Here we consider how much of mechanistic detail needs to be included for in silico models to accurately predict risk and further, how much of this detail can be retrieved from protocols that are practical to implement in high throughout screens as part of next generation of preclinical in silico drug screening approaches?

Introduction The Comprehensive in vitro Proarrhythmic Assay (CiPA) aims to update current cardiac safety testing to better evaluate arrhythmic risk. A central theme of CiPA is the use of in silico approaches to risk prediction incorporating models of drug binding to hERG. To parameterize these models, accurate in vitro measurement of potency and kinetics of block is required. The Ion Channel Working Group was tasked with: i) selecting a protocol that could measure kinetics of block and was easily implementable on automated platforms for future rollout in industry and ii) acquiring a reference dataset using the standardized protocol. Methods Data were acquired using a ‘step depolarisation’ protocol using manual patch-clamp at ambient temperature. Results Potency, kinetics and trapping characteristics of hERG block for the CiPA training panel of twelve drugs were measured. Timecourse of block and trapping characteristics could be reliably measured if the time constant for onset of block was between ~ 500 ms and ~ 15 s. Seven drugs, however had time courses of block faster than this cut-off. Discussion Here we describe the implementation of the standardized protocol for measurement of kinetics and potency of hERG block for CiPA. The results highlight the challenges in identifying a single protocol to measure hERG block over a range of kinetics. The dataset from this study is being used by the In Silico Working Group to develop models of drug binding for risk prediction and is freely available as a ‘gold standard’ ambient temperature dataset to evaluate variability across high throughput platforms.

© 2017 Elsevier Ltd κ-Hexatoxins (κ-HXTXs) are a family of excitotoxic insect-selective neurotoxins from Australian funnel-web spiders that are lethal to a wide range of insects, but display no toxicity towards vertebrates. The prototypic κ-HXTX-Hv1c selectively blocks native and expressed cockroach large-conductance calcium-activated potassium (BKCa or KCa1.1) channels, but not their mammalian orthologs. Despite this potent and selective action on insect KCa1.1 channels, we found that the classical KCa1.1 blockers paxilline, charybdotoxin and iberiotoxin, which all block insect KCa1.1 channels, are not lethal in crickets. We therefore used whole-cell patch-clamp analysis of cockroach dorsal unpaired median (DUM) neurons to study the effects of κ-HXTX-Hv1c on sodium-activated (KNa), delayed-rectifier (KDR) and ‘A-type’ transient (KA) K+ channels. 1 μM κ-HXTX-Hv1c failed to significantly inhibit cockroach KNa and KDR channels, but did cause a 30 ± 7% saturating inhibition of KA channel currents, possibly via a Kv4 (Shal-like) action. However, this modest action at such a high concentration of κ-HXTX-Hv1c would indicate a different lethal target. Accordingly, we assessed the actions of κ-HXTX-Hv1c on neurotransmitter-gated ion channels in cockroach DUM neurons. We found that κ-HXTX-Hv1c failed to produce any major effects on GABAA or glutamate-Cl receptors but dramatically slowed nicotine-evoked ACh receptor (nAChR) current decay and reversed nAChR desensitization. These actions occurred without any alterations to nAChR current amplitude or the nicotine concentration-response curve, and are consistent with a positive allosteric modulation of nAChRs. κ-HXTX-Hv1c therefore represents the first venom peptide that selectively modulates insect nAChRs with a mode of action similar to the excitotoxic insecticide spinosyn A. This article is part of the Special Issue entitled ‘Venom-derived Peptides as Pharmacological Tools.’

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.