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.

High throughput measurement of hERG drug block kinetics using the CiPA dynamic protocol.

The Comprehensive in vitro Proarrhythmic Assay (CiPA) has promoted use of in silico models of drug effects on cardiac repolarization to improve proarrhythmic risk prediction. These models contain a pharmacodynamic component describing drug binding to hERG channels that required in vitro data for kinetics of block, in addition to potency, to constrain them. To date, development and validation has been undertaken using data from manual patch-clamp. The application of this approach at scale requires the development of a high-throughput, automated patch-clamp (APC) implementation. Here, we present a comprehensive analysis of the implementation of the Milnes, or CiPA dynamic protocol, on an APC platform, including quality control and data analysis. Kinetics and potency of block were assessed for bepridil, cisapride, terfenadine and verapamil with data retention/QC pass rate of 21.8% overall, or as high as 50.4% when only appropriate sweep lengths were considered for drugs with faster kinetics. The variability in IC50 and kinetics between manual and APC was comparable to that seen between sites/platforms in previous APC studies of potency. Whilst the experimental success is less than observed in screens of potency alone, it is still significantly greater than manual patch. With the modifications to protocol design, including sweep length, number of repetitions, and leak correction recommended in this study, this protocol can be applied on APC to acquire data comparable to manual patch clamp.

Metabolic and electrolyte abnormalities as risk factors in drug-induced long QT syndrome.

Drug-induced long QT syndrome (diLQTS) is the phenomenon by which the administration of drugs causes prolongation of cardiac repolarisation and leads to an increased risk of the ventricular tachycardia known as torsades de pointes (TdP). In most cases of diLQTS, the primary molecular target is the human ether-à-go-go-related gene protein (hERG) potassium channel, which carries the rapid delayed rectifier current (IKr) in the heart. However, the proarrhythmic risk associated with drugs that block hERG can be modified in patients by a range of environmental- and disease-related factors, such as febrile temperatures, alterations in pH, dyselectrolytaemias such as hypokalaemia and hypomagnesemia and coadministration with other drugs. In this review, we will discuss the clinical occurrence of drug-induced LQTS in the context of these modifying factors as well as the mechanisms by which they contribute to altered hERG potency and proarrhythmic risk.

Pathophysiological metabolic changes associated with disease modify the proarrhythmic risk profile of drugs with potential to prolong repolarisation.

BACKGROUND AND PURPOSE: Hydroxychloroquine, chloroquine and azithromycin are three drugs that were proposed to treat coronavirus disease 2019 (COVID-19). While concern already existed around their proarrhythmic potential, there are little data regarding how altered physiological states encountered in patients such as febrile state, electrolyte imbalances or acidosis might change their risk profiles. EXPERIMENTAL APPROACH: Potency of human ether-à-go-go related gene (hERG) block was measured using high-throughput electrophysiology in the presence of variable environmental factors. These potencies informed simulations to predict population risk profiles. Effects on cardiac repolarisation were verified in human induced pluripotent stem cell-derived cardiomyocytes from multiple individuals. KEY RESULTS: Chloroquine and hydroxychloroquine blocked hERG with IC50 of 1.47 ± 0.07 and 3.78 ± 0.17 µM, respectively, indicating proarrhythmic risk at concentrations effective against severe acute respiratory syndrome-coronovirus-2 (SARS-CoV-2) in vitro. Hypokalaemia and hypermagnesaemia increased potency of chloroquine and hydroxychloroquine, indicating increased proarrhythmic risk. Acidosis significantly reduced potency of all drugs, whereas increased temperature decreased potency of chloroquine and hydroxychloroquine against hERG but increased potency for azithromycin. In silico simulations demonstrated that proarrhythmic risk was increased by female sex, hypokalaemia and heart failure and identified specific genetic backgrounds associated with emergence of arrhythmia. CONCLUSION AND IMPLICATIONS: Our study demonstrates how proarrhythmic risk can be exacerbated by metabolic changes and pre-existing disease. More broadly, the study acts as a blueprint for how high-throughput in vitro screening, combined with in silico simulations, can help guide both preclinical screening and clinical management of patients in relation to drugs with potential to prolong repolarisation.

Heterozygous KCNH2 variant phenotyping using Flp-In HEK293 and high-throughput automated patch clamp electrophysiology.

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.

A constricted opening in Kir channels does not impede potassium conduction.

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.

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.

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.

Experimentally Validated Pharmacoinformatics Approach to Predict hERG Inhibition Potential of New Chemical Entities.

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.

Evaluation of Chemical Strategies for Improving the Stability and Oral Toxicity of Insecticidal Peptides.

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

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.

In Vitro and In Silico Risk Assessment in Acquired Long QT Syndrome: The Devil Is in the Details.

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?

Lethal effects of an insecticidal spider venom peptide involve positive allosteric modulation of insect nicotinic acetylcholine receptors.

κ-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.'

Measuring kinetics and potency of hERG block for CiPA.

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 ~500ms and ~15s. 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.

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.

In silico assessment of kinetics and state dependent binding properties of drugs causing acquired LQTS.

The Kv11.1 or hERG potassium channel is responsible for one of the major repolarising currents (IKr) in cardiac myocytes. Drug binding to hERG can result in reduction in IKr, action potential prolongation, acquired long QT syndrome and fatal cardiac arrhythmias. The current guidelines for pre-clinical assessment of drugs in development is based on the measurement of the drug concentration that causes 50% current block, i.e., IC50. However, drugs with the same apparent IC50 may have very different kinetics of binding and unbinding, as well as different affinities for the open and inactivated states of Kv11.1. Therefore, IC50 measurements may not reflect the true risk of drug induced arrhythmias. Here we have used an in silico approach to test the hypothesis that drug binding kinetics and differences in state-dependent affinity will influence the extent of cardiac action potential prolongation independent of apparent IC50 values. We found, in general that drugs with faster overall kinetics and drugs with higher affinity for the open state relative to the inactivated state cause more action potential prolongation. These characteristics of drug-hERG interaction are likely to be more arrhythmogenic but cannot be predicted by IC50 measurement alone. Our results suggest that the pre-clinical assessment of Kv11.1-drug interactions should include descriptions of the kinetics and state dependence of drug binding. Further, incorporation of this information into sophisticated in silico models should be able to better predict arrhythmia risk and therefore more accurately assess safety of new drugs in development.

Do vicinal disulfide bridges mediate functionally important redox transformations in proteins?

Vicinal disulfide bridges, in which a disulfide bond is formed between adjacent cysteine residues, constitute an unusual but expanding class of potential allosteric disulfides. Although vicinal disulfide rings (VDRs) are relatively uncommon, they have proven to be functionally critical in almost all proteins in which they have been discovered. However, it has proved difficult to test whether these sterically constrained disulfides participate in functionally important redox transformations. We demonstrate that chemical replacement of VDRs with dicarba or diselenide bridges can be used to assess whether VDRs function as allosteric disulfides. Our approach leads to the hypothesis that not all VDRs participate in functionally important redox reactions.

Spider-venom peptides as bioinsecticides.

Over 10,000 arthropod species are currently considered to be pest organisms. They are estimated to contribute to the destruction of ~14% of the world's annual crop production and transmit many pathogens. Presently, arthropod pests of agricultural and health significance are controlled predominantly through the use of chemical insecticides. Unfortunately, the widespread use of these agrochemicals has resulted in genetic selection pressure that has led to the development of insecticide-resistant arthropods, as well as concerns over human health and the environment. Bioinsecticides represent a new generation of insecticides that utilise organisms or their derivatives (e.g., transgenic plants, recombinant baculoviruses, toxin-fusion proteins and peptidomimetics) and show promise as environmentally-friendly alternatives to conventional agrochemicals. Spider-venom peptides are now being investigated as potential sources of bioinsecticides. With an estimated 100,000 species, spiders are one of the most successful arthropod predators. Their venom has proven to be a rich source of hyperstable insecticidal mini-proteins that cause insect paralysis or lethality through the modulation of ion channels, receptors and enzymes. Many newly characterized insecticidal spider toxins target novel sites in insects. Here we review the structure and pharmacology of these toxins and discuss the potential of this vast peptide library for the discovery of novel bioinsecticides.

A novel family of insect-selective peptide neurotoxins targeting insect large-conductance calcium-activated K+ channels isolated from the venom of the theraphosid spider Eucratoscelus constrictus.

Spider venoms are actively being investigated as sources of novel insecticidal agents for biopesticide engineering. After screening 37 theraphosid spider venoms, a family of three new "short-loop" inhibitory cystine knot insecticidal toxins (κ-TRTX-Ec2a, κ-TRTX-Ec2b, and κ-TRTX-Ec2c) were isolated and characterized from the venom of the African tarantula Eucratoscelus constrictus. Whole-cell patch-clamp recordings from cockroach dorsal unpaired median neurons revealed that, despite significant sequence homology with other theraphosid toxins, these 29-residue peptides lacked activity on insect voltage-activated sodium and calcium channels. It is noteworthy that κ-TRTX-Ec2 toxins were all found to be high-affinity blockers of insect large-conductance calcium-activated K(+) (BK(Ca)) channel currents with IC(50) values of 3 to 25 nM. In addition, κ-TRTX-Ec2a caused the inhibition of insect delayed-rectifier K(+) currents, but only at significantly higher concentrations. κ-TRTX-Ec2a and κ-TRTX-Ec2b demonstrated insect-selective effects, whereas the homologous κ-TRTX-Ec2c also resulted in neurotoxic signs in mice when injected intracerebroventricularly. Unlike other theraphosid toxins, κ-TRTX-Ec2 toxins induce a voltage-independent channel block, and therefore, we propose that these toxins interact with the turret and/or loop region of the external entrance to the channel and do not project deeply into the pore of the channel. Furthermore, κ-TRTX-Ec2a and κ-TRTX-Ec2b differ from other theraphotoxins at the C terminus and positions 5 to 6, suggesting that these regions of the peptide contribute to the phyla selectivity and are involved in targeting BK(Ca) channels. This study therefore establishes these toxins as tools for studying the role of BK(Ca) channels in insects and lead compounds for the development of novel insecticides.