Allosteric mechanism for KCNE1 modulation of KCNQ1 potassium channel activation

Dynamic protein-protein interactions of KCNQ1 and KCNE1 measured by EPR line shape analysis

KCNQ1, also known as Kv7.1, is a voltage gated potassium channel that associates with the KCNE protein family. Mutations in this protein has been found to cause a variety of diseases including Long QT syndrome, a type of cardiac arrhythmia where the QT interval observed on an electrocardiogram is longer than normal. This condition is often aggravated during strenuous exercise and can cause fainting spells or sudden death. KCNE1 is an ancillary protein that interacts with KCNQ1 in the membrane at varying molar ratios. This interaction allows for the flow of potassium ions to be modulated to facilitate repolarization of the heart. The interaction between these two proteins has been studied previously with cysteine crosslinking and electrophysiology. In this study, electron paramagnetic resonance (EPR) spectroscopy line shape analysis in tandem with site directed spin labeling (SDSL) was used to observe changes in side chain dynamics as KCNE1 interacts with KCNQ1. KCNE1 was labeled at different sites that were found to interact with KCNQ1 based on previous literature, along with sites outside of that range as a control. Once labeled KCNE1 was incorporated into vesicles, KCNQ1 (helices S1-S6) was titrated into the vesicles. The line shape differences observed upon addition of KCNQ1 are indicative of an interaction between the two proteins. This method provides a first look at the interactions between KCNE1 and KCNQ1 from a dynamics perspective using the full transmembrane portion of KCNQ1.

High-throughput functional mapping of variants in an arrhythmia gene, KCNE1, reveals novel biology

BACKGROUND

KCNE1 encodes a 129-residue cardiac potassium channel (IKs) subunit. KCNE1 variants are associated with long QT syndrome and atrial fibrillation. However, most variants have insufficient evidence of clinical consequences and thus limited clinical utility.

METHODS

In this study, we leveraged the power of variant effect mapping, which couples saturation mutagenesis with high-throughput sequencing, to ascertain the function of thousands of protein-coding KCNE1 variants.

RESULTS

We comprehensively assayed KCNE1 variant cell surface expression (2554/2709 possible single-amino-acid variants) and function (2534 variants). Our study identified 470 loss- or partial loss-of-surface expression and 574 loss- or partial loss-of-function variants. Of the 574 loss- or partial loss-of-function variants, 152 (26.5%) had reduced cell surface expression, indicating that most functionally deleterious variants affect channel gating. Nonsense variants at residues 56-104 generally had WT-like trafficking scores but decreased functional scores, indicating that the latter half of the protein is dispensable for protein trafficking but essential for channel function. 22 of the 30 KCNE1 residues (73%) highly intolerant of variation (with > 70% loss-of-function variants) were in predicted close contact with binding partners KCNQ1 or calmodulin. Our functional assay data were consistent with gold standard electrophysiological data (ρ =  - 0.64), population and patient cohorts (32/38 presumed benign or pathogenic variants with consistent scores), and computational predictors (ρ =  - 0.62). Our data provide moderate-strength evidence for the American College of Medical Genetics/Association of Molecular Pathology functional criteria for benign and pathogenic variants.

CONCLUSIONS

Comprehensive variant effect maps of KCNE1 can both provide insight into I Ks channel biology and help reclassify variants of uncertain significance.

Evaluating sequential and allosteric activation models in IKs channels with mutated voltage sensors

The ion-conducting IKs channel complex, important in cardiac repolarization and arrhythmias, comprises tetramers of KCNQ1 α-subunits along with 1-4 KCNE1 accessory subunits and calmodulin regulatory molecules. The E160R mutation in individual KCNQ1 subunits was used to prevent activation of voltage sensors and allow direct determination of transition rate data from complexes opening with a fixed number of 1, 2, or 4 activatable voltage sensors. Markov models were used to test the suitability of sequential versus allosteric models of IKs activation by comparing simulations with experimental steady-state and transient activation kinetics, voltage-sensor fluorescence from channels with two or four activatable domains, and limiting slope currents at negative potentials. Sequential Hodgkin-Huxley-type models approximately describe IKs currents but cannot explain an activation delay in channels with only one activatable subunit or the hyperpolarizing shift in the conductance-voltage relationship with more activatable voltage sensors. Incorporating two voltage sensor activation steps in sequential models and a concerted step in opening via rates derived from fluorescence measurements improves models but does not resolve fundamental differences with experimental data. Limiting slope current data that show the opening of channels at negative potentials and very low open probability are better simulated using allosteric models of activation with one transition per voltage sensor, which implies that movement of all four sensors is not required for IKs conductance. Tiered allosteric models with two activating transitions per voltage sensor can fully account for IKs current and fluorescence activation kinetics in constructs with different numbers of activatable voltage sensors.

Whole-exome sequencing of pathogenic genes in a family with congenital heart disease: A case report

RATIONALE

Congenital heart disease (CHD) is the most common birth defect and an important cause of noninfectious deaths in infants and children. It has high prevalence globally, placing an enormous burden on society and families. Studies of individuals with hereditary or sporadic CHD have provided strong evidence for its genetic basis. The aim of this study was to identify causative gene variants in a Chinese family with congenital heart disease.

PATIENT CONCERNS AND DIAGNOSES

Three generations of a CHD family were recruited. Proband III.9 was diagnosed with congenital heart disease at age 11 months, and the echocardiogram showed arterial ductus arteriosus, with a left-to-right shunt at the level of the arteries. Precedent III.10 was a twin of Proband III.9 who was diagnosed with congenital heart disease at age 11 months, in whom the echocardiogram revealed an arterial ductus arteriosus, an unenclosed patent ductus arteriosus, and a left to right shunt at the level of the arteries (second figure). III.8 was diagnosed with congenital heart disease at age 15, but echocardiography in this study showed no abnormalities. No cardiac abnormalities were detected in any of his parents, grandparents, or maternal grandparents. We performed whole-exome sequencing on CHD sufferers and their unexpressing family members to investigate the genetic causes of CHD in this family line. Exome sequencing identified 4 mutation sites in this family line. The variant c.3245A>G (p.His1082Arg) of the AMER1 gene was consistent with concomitant X-chromosome recessive inheritance, the variant c.238G>C (p.Val80Leu) of the KCNE1 gene was consistent with autosomal accessory inheritance, and the other 2 variants did not conform to the law of the mode of inheritance of the disease.

OUTCOMES

The first identified variant, c.3245A>G (p.His1082Arg) of the AMER1 gene, with X-chromosome recessive inheritance, and the variant c.238G>C (p.Val80Leu) of the KCNE1 gene, which has been reported as autosomal dominant, may be the causative agent of CHD in this family line. These findings broaden the genetic scope of congenital heart disease and could help in the development of targeted drugs for the treatment of congenital heart disease.

Dark nanodiscs for evaluating membrane protein thermostability by differential scanning fluorimetry

Measuring protein thermostability provides valuable information on the biophysical rules that govern the structure-energy relationships of proteins. However, such measurements remain a challenge for membrane proteins. Here, we introduce a new experimental system to evaluate membrane protein thermostability. This system leverages a recently developed nonfluorescent membrane scaffold protein to reconstitute proteins into nanodiscs and is coupled with a nano-format of differential scanning fluorimetry (nanoDSF). This approach offers a label-free and direct measurement of the intrinsic tryptophan fluorescence of the membrane protein as it unfolds in solution without signal interference from the "dark" nanodisc. In this work, we demonstrate the application of this method using the disulfide bond formation protein B (DsbB) as a test membrane protein. NanoDSF measurements of DsbB reconstituted in dark nanodiscs loaded with 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and 1,2-dimyristoyl-sn-glycero-3-phosphorylglycerol (DMPG) lipids show a complex biphasic thermal unfolding pattern with a minor unfolding transition followed by a major transition. The inflection points of the thermal denaturation curve reveal two distinct unfolding midpoint melting temperatures (Tm) of 70.5°C and 77.5°C, consistent with a three-state unfolding model. Further, we show that the catalytically conserved disulfide bond between residues C41 and C130 drives the intermediate state of the unfolding pathway for DsbB in a DMPC and DMPG nanodisc. To extend the utility of this method, we evaluate and compare the thermostability of DsbB in different lipid environments. We introduce this method as a new tool that can be used to understand how compositionally and biophysically complex lipid environments drive membrane protein stability.

A generic binding pocket for small molecule IKs activators at the extracellular inter-subunit interface of KCNQ1 and KCNE1 channel complexes

The cardiac IKs ion channel comprises KCNQ1, calmodulin, and KCNE1 in a dodecameric complex which provides a repolarizing current reserve at higher heart rates and protects from arrhythmia syndromes that cause fainting and sudden death. Pharmacological activators of IKs are therefore of interest both scientifically and therapeutically for treatment of IKs loss-of-function disorders. One group of chemical activators are only active in the presence of the accessory KCNE1 subunit and here we investigate this phenomenon using molecular modeling techniques and mutagenesis scanning in mammalian cells. A generalized activator binding pocket is formed extracellularly by KCNE1, the domain-swapped S1 helices of one KCNQ1 subunit and the pore/turret region made up of two other KCNQ1 subunits. A few residues, including K41, A44 and Y46 in KCNE1, W323 in the KCNQ1 pore, and Y148 in the KCNQ1 S1 domain, appear critical for the binding of structurally diverse molecules, but in addition, molecular modeling studies suggest that induced fit by structurally different molecules underlies the generalized nature of the binding pocket. Activation of IKs is enhanced by stabilization of the KCNQ1-S1/KCNE1/pore complex, which ultimately slows deactivation of the current, and promotes outward current summation at higher pulse rates. Our results provide a mechanistic explanation of enhanced IKs currents by these activator compounds and provide a map for future design of more potent therapeutically useful molecules.

Dark nanodiscs as a model membrane for evaluating membrane protein thermostability by differential scanning fluorimetry

Measuring protein thermostability provides valuable information on the biophysical rules that govern structure-energy relationships of proteins. However, such measurements remain a challenge for membrane proteins. Here, we introduce a new experimental system to evaluate membrane protein thermostability. This system leverages a recently-developed non-fluorescent membrane scaffold protein (MSP) to reconstitute proteins into nanodiscs and is coupled with a nano-format of differential scanning fluorimetry (nanoDSF). This approach offers a label-free and direct measurement of the intrinsic tryptophan fluorescence of the membrane protein as it unfolds in solution without signal interference from the "dark" nanodisc. In this work, we demonstrate the application of this method using the disulfide bond formation protein B (DsbB) as a test membrane protein. NanoDSF measurements of DsbB reconstituted in dark nanodiscs show a complex biphasic thermal unfolding pattern in the presence of lipids with a minor unfolding transition followed by a major transition. The inflection points of the thermal denaturation curve reveal two distinct unfolding midpoint melting temperatures (Tm) of 70.5 °C and 77.5 °C, consistent with a three-state unfolding model. Further, we show that the catalytically conserved disulfide bond between residues C41 and C130 drives the intermediate state of the unfolding pathway for DsbB in a nanodisc. We introduce this method as a new tool that can be used to understand how compositionally, and biophysically complex lipid environments drive membrane protein stability.

High-throughput functional mapping of variants in an arrhythmia gene, KCNE1, reveals novel biology

Background

KCNE1 encodes a 129-residue cardiac potassium channel (IKs) subunit. KCNE1 variants are associated with long QT syndrome and atrial fibrillation. However, most variants have insufficient evidence of clinical consequences and thus limited clinical utility.

Results

Here, we demonstrate the power of variant effect mapping, which couples saturation mutagenesis with high-throughput sequencing, to ascertain the function of thousands of protein coding KCNE1 variants. We comprehensively assayed KCNE1 variant cell surface expression (2,554/2,709 possible single amino acid variants) and function (2,539 variants). We identified 470 loss-of-surface expression and 588 loss-of-function variants. Out of the 588 loss-of-function variants, only 155 had low cell surface expression. The latter half of the protein is dispensable for protein trafficking but essential for channel function. 22 of the 30 KCNE1 residues (73%) highly intolerant of variation were in predicted close contact with binding partners KCNQ1 or calmodulin. Our data were highly concordant with gold standard electrophysiological data (ρ = -0.65), population and patient cohorts (32/38 concordant variants), and computational metrics (ρ = -0.55). Our data provide moderate-strength evidence for the ACMG/AMP functional criteria for benign and pathogenic variants.

Conclusions

Comprehensive variant effect maps of KCNE1 can both provide insight into IKs channel biology and help reclassify variants of uncertain significance.

Long-QT mutations in KCNE1 modulate the 17β-estradiol response of Kv7.1/KCNE1

Estradiol (17[Formula: see text]-E2) is implicated in higher arrhythmia risk of women with congenital or acquired long-QT syndrome (LQTS) compared to men. However, the underlying mechanisms remain poorly understood, and little is known about the impact of LQTS-associated mutations. We show that 17[Formula: see text]-E2 inhibits the human cardiac Kv7.1/KCNE1 channel expressed in Xenopus oocytes. We find that the 17[Formula: see text]-E2 effect depends on the Kv7.1 to KCNE1 stoichiometry, and we reveal a critical function of the KCNE1 carboxyl terminus for the effect. LQTS-associated mutations in the KCNE1 carboxyl terminus show a range of responses to 17[Formula: see text]-E2, from a wild-type like response to impaired or abolished response. Together, this study increases our understanding of the mechanistic basis for 17[Formula: see text]-E2 inhibition of Kv7.1/KCNE1 and demonstrates mutation-dependent responses to 17[Formula: see text]-E2. These findings suggest that the 17[Formula: see text]-E2 effect on Kv7.1/KCNE1 might contribute to the higher arrhythmia risk of women, particularly in carriers with specific LQTS-associated mutations.

Simulation and Machine Learning Methods for Ion-Channel Structure Determination, Mechanistic Studies and Drug Design

Ion channels are expressed in almost all living cells, controlling the in-and-out communications, making them ideal drug targets, especially for central nervous system diseases. However, owing to their dynamic nature and the presence of a membrane environment, ion channels remain difficult targets for the past decades. Recent advancement in cryo-electron microscopy and computational methods has shed light on this issue. An explosion in high-resolution ion channel structures paved way for structure-based rational drug design and the state-of-the-art simulation and machine learning techniques dramatically improved the efficiency and effectiveness of computer-aided drug design. Here we present an overview of how simulation and machine learning-based methods fundamentally changed the ion channel-related drug design at different levels, as well as the emerging trends in the field.

Cholesterol-induced suppression of Kir2 channels is mediated by decoupling at the inter-subunit interfaces

Cholesterol is a major regulator of multiple types of ion channels. Although there is increasing information about cholesterol binding sites, the molecular mechanisms through which cholesterol binding alters channel function are virtually unknown. In this study, we used a combination of Martini coarse-grained simulations, a network theory-based analysis, and electrophysiology to determine the effect of cholesterol on the dynamic structure of the Kir2.2 channel. We found that increasing membrane cholesterol reduced the likelihood of contact between specific regions of the cytoplasmic and transmembrane domains of the channel, most prominently at the subunit-subunit interfaces of the cytosolic domains. This decrease in contact was mediated by pairwise interactions of specific residues and correlated to the stoichiometry of cholesterol binding events. The predictions of the model were tested by site-directed mutagenesis of two identified residues-V265 and H222-and high throughput electrophysiology.

High-throughput characterization of photocrosslinker-bearing ion channel variants to map residues critical for function and pharmacology

Incorporation of noncanonical amino acids (ncAAs) can endow proteins with novel functionalities, such as crosslinking or fluorescence. In ion channels, the function of these variants can be studied with great precision using standard electrophysiology, but this approach is typically labor intensive and low throughput. Here, we establish a high-throughput protocol to conduct functional and pharmacological investigations of ncAA-containing human acid-sensing ion channel 1a (hASIC1a) variants in transiently transfected mammalian cells. We introduce 3 different photocrosslinking ncAAs into 103 positions and assess the function of the resulting 309 variants with automated patch clamp (APC). We demonstrate that the approach is efficient and versatile, as it is amenable to assessing even complex pharmacological modulation by peptides. The data show that the acidic pocket is a major determinant for current decay, and live-cell crosslinking provides insight into the hASIC1a–psalmotoxin 1 (PcTx1) interaction. Further, we provide evidence that the protocol can be applied to other ion channels, such as P2X2 and GluA2 receptors. We therefore anticipate the approach to enable future APC-based studies of ncAA-containing ion channels in mammalian cells.

This study describes a method to rapidly screen hundreds of ion channel variants containing non-canonical amino acids. A proof-of-principle introducing photocrosslinking non-canonical amino acids into the human ion channel hASIC1a shows how this approach can provide insights into function and pharmacology.

A general mechanism of KCNE1 modulation of KCNQ1 channels involving non-canonical VSD-PD coupling

Voltage-gated KCNQ1 channels contain four separate voltage-sensing domains (VSDs) and a pore domain (PD). KCNQ1 expressed alone opens when the VSDs are in an intermediate state. In cardiomyocytes, KCNQ1 co-expressed with KCNE1 opens mainly when the VSDs are in a fully activated state. KCNE1 also drastically slows the opening of KCNQ1 channels and shifts the voltage dependence of opening by >40 mV. We here show that mutations of conserved residues at the VSD-PD interface alter the VSD-PD coupling so that the mutant KCNQ1/KCNE1 channels open in the intermediate VSD state. Using recent structures of KCNQ1 and KCNE beta subunits in different states, we present a mechanism by which KCNE1 rotates the VSD relative to the PD and affects the VSD-PD coupling of KCNQ1 channels in a non-canonical way, forcing KCNQ1/KCNE1 channels to open in the fully-activated VSD state. This would explain many of the KCNE1-induced effects on KCNQ1 channels.

Allosteric mechanism for KCNE1 modulation of KCNQ1 potassium channel activation

The function of the voltage-gated KCNQ1 potassium channel is regulated by co-assembly with KCNE auxiliary subunits. KCNQ1-KCNE1 channels generate the slow delayed rectifier current, I_Ks, which contributes to the repolarization phase of the cardiac action potential. A three amino acid motif (F57-T58-L59, FTL) in KCNE1 is essential for slow activation of KCNQ1-KCNE1 channels. However, how this motif interacts with KCNQ1 to control its function is unknown. Combining computational modeling with electrophysiological studies, we developed structural models of the KCNQ1-KCNE1 complex that suggest how KCNE1 controls KCNQ1 activation. The FTL motif binds at a cleft between the voltage-sensing and pore domains and appears to affect the channel gate by an allosteric mechanism. Comparison with the KCNQ1-KCNE3 channel structure suggests a common transmembrane-binding mode for different KCNEs and illuminates how specific differences in the interaction of their triplet motifs determine the profound differences in KCNQ1 functional modulation by KCNE1 versus KCNE3.