Clinical interpretation of KCNH2 variants using a robust PS3/BS3 functional patch-clamp assay

Long QT syndrome (LQTS), caused by the dysfunction of cardiac ion channels, increases the risk of sudden death in otherwise healthy young people. For many variants in LQTS genes, there is insufficient evidence to make a definitive genetic diagnosis. We have established a robust functional patch-clamp assay to facilitate classification of missense variants in KCNH2, one of the key LQTS genes. A curated set of 30 benign and 30 pathogenic missense variants were used to establish the range of normal and abnormal function. The extent to which variants reduced protein function was quantified using Z scores, the number of standard deviations from the mean of the normalized current density of the set of benign variant controls. A Z score of -2 defined the threshold for abnormal loss of function, which corresponds to 55% wild-type function. More extreme Z scores were observed for variants with a greater loss-of-function effect. We propose that the Z score for each variant can be used to inform the application and weighting of abnormal and normal functional evidence criteria (PS3 and BS3) within the American College of Medical Genetics and Genomics variant classification framework. The validity of this approach was demonstrated using a series of 18 KCNH2 missense variants detected in a childhood onset LQTS cohort, where the level of function assessed using our assay correlated to the Schwartz score (a scoring system used to quantify the probability of a clinical diagnosis of LQTS) and the length of the corrected QT (QTc) interval.

Prognostic Value of Multiplexed Assays of Variant Effect and Automated Patch-clamping for KCNH2-LQTS Risk Stratification

Background

Long QT syndrome (LQTS) is a lethal arrhythmia condition, frequently caused by rare loss-of-function variants in the cardiac potassium channel encoded by KCNH2. Variant-based risk stratification is complicated by heterogenous clinical data, incomplete penetrance, and low-throughput functional data.

Objective

To test the utility of variant-specific features, including high-throughput functional data, to predict cardiac events among KCNH2 variant heterozygotes.

Methods

We quantified cell-surface trafficking of 18,323 variants in KCNH2 and recorded potassium current densities for 506 KCNH2 variants. Next, we deeply phenotyped 1150 KCNH2 missense variant patients, including ECG features, cardiac event history (528 total cardiac events), and mortality. We then assessed variant functional, in silico, structural, and LQTS penetrance data to stratify event-free survival for cardiac events in the study cohort.

Results

Variant-specific current density (HR 0.28 [0.13-0.60]) and estimates of LQTS penetrance incorporating MAVE data (HR 3.16 [1.59-6.27]) were independently predictive of severe cardiac events when controlling for patient-specific features. Risk prediction models incorporating these data significantly improved prediction of 20 year cardiac events (AUC 0.79 [0.75-0.82]) over patient-only covariates (QTc and sex) (AUC 0.73 [0.70-0.77]).

Conclusion

We show that high-throughput functional data, and other variant-specific features, meaningfully contribute to both diagnosis and prognosis of a clinically actionable monogenic disease.

Multicenter clinical and functional evidence reclassifies a recurrent noncanonical filamin C splice-altering variant

BACKGROUND

Truncating variants in filamin C (FLNC) can cause arrhythmogenic cardiomyopathy (ACM) through haploinsufficiency. Noncanonical splice-altering variants may contribute to this phenotype.

OBJECTIVE

The purpose of this study was to investigate the clinical and functional consequences of a recurrent FLNC intronic variant of uncertain significance (VUS), c.970-4A>G.

METHODS

Clinical data in 9 variant heterozygotes from 4 kindreds were obtained from 5 tertiary health care centers. We used in silico predictors and functional studies with peripheral blood and patient-specific induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs). Isolated RNA was studied by reverse transcription polymerase chain reaction. iPSC-CMs were further characterized at baseline and after nonsense-mediated decay (NMD) inhibition, using quantitative polymerase chain reaction (qPCR), RNA-sequencing, and cellular electrophysiology. American College of Medical Genetics and Genomics (ACMG) criteria were used to adjudicate variant pathogenicity.

RESULTS

Variant heterozygotes displayed a spectrum of disease phenotypes, spanning from mild ventricular dysfunction with palpitations to severe ventricular arrhythmias requiring device shocks or progressive cardiomyopathy requiring heart transplantation. Consistent with in silico predictors, the c.970-4A>G FLNC variant activated a cryptic splice acceptor site, introducing a 3-bp insertion containing a premature termination codon. NMD inhibition upregulated aberrantly spliced transcripts by qPCR and RNA-sequencing. Patch clamp studies revealed irregular spontaneous action potentials, increased action potential duration, and increased sodium late current in proband-derived iPSC-CMs. These findings fulfilled multiple ACMG criteria for pathogenicity.

CONCLUSION

Clinical, in silico, and functional evidence support the prediction that the intronic c.970-4A>G VUS disrupts splicing and drives ACM, enabling reclassification from VUS to pathogenic.

Continuous Bayesian variant interpretation accounts for incomplete penetrance among Mendelian cardiac channelopathies

PURPOSE

The congenital Long QT Syndrome (LQTS) and Brugada Syndrome (BrS) are Mendelian autosomal dominant diseases that frequently precipitate fatal cardiac arrhythmias. Incomplete penetrance is a barrier to clinical management of heterozygotes harboring variants in the major implicated disease genes KCNQ1, KCNH2, and SCN5A. We apply and evaluate a Bayesian penetrance estimation strategy that accounts for this phenomenon.

METHODS

We generated Bayesian penetrance models for KCNQ1-LQT1 and SCN5A-LQT3 using variant-specific features and clinical data from the literature, international arrhythmia genetic centers, and population controls. We analyzed the distribution of posterior penetrance estimates across 4 genotype-phenotype relationships and compared continuous estimates with ClinVar annotations. Posterior estimates were mapped onto protein structure.

RESULTS

Bayesian penetrance estimates of KCNQ1-LQT1 and SCN5A-LQT3 are empirically equivalent to 10 and 5 clinically phenotype heterozygotes, respectively. Posterior penetrance estimates were bimodal for KCNQ1-LQT1 and KCNH2-LQT2, with a higher fraction of missense variants with high penetrance among KCNQ1 variants. There was a wide distribution of variant penetrance estimates among identical ClinVar categories. Structural mapping revealed heterogeneity among "hot spot" regions and featured high penetrance estimates for KCNQ1 variants in contact with calmodulin and the S6 domain.

CONCLUSIONS

Bayesian penetrance estimates provide a continuous framework for variant interpretation.

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.

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.

Arrhythmia Variant Associations and Reclassifications in the eMERGE-III Sequencing Study

BACKGROUND

Sequencing Mendelian arrhythmia genes in individuals without an indication for arrhythmia genetic testing can identify carriers of pathogenic or likely pathogenic (P/LP) variants. However, the extent to which these variants are associated with clinically meaningful phenotypes before or after return of variant results is unclear. In addition, the majority of discovered variants are currently classified as variants of uncertain significance, limiting clinical actionability.

METHODS

The eMERGE-III study (Electronic Medical Records and Genomics Phase III) is a multicenter prospective cohort that included 21 846 participants without previous indication for cardiac genetic testing. Participants were sequenced for 109 Mendelian disease genes, including 10 linked to arrhythmia syndromes. Variant carriers were assessed with electronic health record-derived phenotypes and follow-up clinical examination. Selected variants of uncertain significance (n=50) were characterized in vitro with automated electrophysiology experiments in HEK293 cells.

RESULTS

As previously reported, 3.0% of participants had P/LP variants in the 109 genes. Herein, we report 120 participants (0.6%) with P/LP arrhythmia variants. Compared with noncarriers, arrhythmia P/LP carriers had a significantly higher burden of arrhythmia phenotypes in their electronic health records. Fifty-four participants had variant results returned. Nineteen of these 54 participants had inherited arrhythmia syndrome diagnoses (primarily long-QT syndrome), and 12 of these 19 diagnoses were made only after variant results were returned (0.05%). After in vitro functional evaluation of 50 variants of uncertain significance, we reclassified 11 variants: 3 to likely benign and 8 to P/LP.

CONCLUSIONS

Genome sequencing in a large population without indication for arrhythmia genetic testing identified phenotype-positive carriers of variants in congenital arrhythmia syndrome disease genes. As the genomes of large numbers of people are sequenced, the disease risk from rare variants in arrhythmia genes can be assessed by integrating genomic screening, electronic health record phenotypes, and in vitro functional studies.

REGISTRATION

URL: https://www.

CLINICALTRIALS

gov; Unique identifier; NCT03394859.

Veratridine Can Bind to a Site at the Mouth of the Channel Pore at Human Cardiac Sodium Channel NaV1.5

The cardiac sodium ion channel (Na_V1.5) is a protein with four domains (DI-DIV), each with six transmembrane segments. Its opening and subsequent inactivation results in the brief rapid influx of Na⁺ ions resulting in the depolarization of cardiomyocytes. The neurotoxin veratridine (VTD) inhibits Na_V1.5 inactivation resulting in longer channel opening times, and potentially fatal action potential prolongation. VTD is predicted to bind at the channel pore, but alternative binding sites have not been ruled out. To determine the binding site of VTD on Na_V1.5, we perform docking calculations and high-throughput electrophysiology experiments in the present study. The docking calculations identified two distinct binding regions. The first site was in the pore, close to the binding site of Na_V1.4 and Na_V1.5 blocking drugs in experimental structures. The second site was at the “mouth” of the pore at the cytosolic side, partly solvent-exposed. Mutations at this site (L409, E417, and I1466) had large effects on VTD binding, while residues deeper in the pore had no effect, consistent with VTD binding at the mouth site. Overall, our results suggest a VTD binding site close to the cytoplasmic mouth of the channel pore. Binding at this alternative site might indicate an allosteric inactivation mechanism for VTD at Na_V1.5.

Estimating the Posttest Probability of Long QT Syndrome Diagnosis for Rare KCNH2 Variants

BACKGROUND

The proliferation of genetic profiling has revealed many associations between genetic variations and disease. However, large-scale phenotyping efforts in largely healthy populations, coupled with DNA sequencing, suggest variants currently annotated as pathogenic are more common in healthy populations than previously thought. In addition, novel and rare variants are frequently observed in genes associated with disease both in healthy individuals and those under suspicion of disease. This raises the question of whether these variants can be useful predictors of disease. To answer this question, we assessed the degree to which the presence of a variant in the cardiac potassium channel gene KCNH2 was diagnostically predictive for the autosomal dominant long QT syndrome.

METHODS

We estimated the probability of a long QT diagnosis given the presence of each KCNH2 variant using Bayesian methods that incorporated variant features such as changes in variant function, protein structure, and in silico predictions. We call this estimate the posttest probability of disease. Our method was applied to over 4000 individuals heterozygous for 871 missense or in-frame insertion/deletion variants in KCNH2 and validated against a separate international cohort of 933 individuals heterozygous for 266 missense or in-frame insertion/deletion variants.

RESULTS

Our method was well-calibrated for the observed fraction of heterozygotes diagnosed with long QT syndrome. Heuristically, we found that the innate diagnostic information one learns about a variant from 3-dimensional variant location, in vitro functional data, and in silico predictors is equivalent to the diagnostic information one learns about that same variant by clinically phenotyping 10 heterozygotes. Most importantly, these data can be obtained in the absence of any clinical observations.

CONCLUSIONS

We show how variant-specific features can inform a prior probability of disease for rare variants even in the absence of clinically phenotyped heterozygotes.

High-Throughput Reclassification of SCN5A Variants

Partial or complete loss-of-function variants in SCN5A are the most common genetic cause of the arrhythmia disorder Brugada syndrome (BrS1). However, the pathogenicity of SCN5A variants is often unknown or disputed; 80% of the 1,390 SCN5A missense variants observed in at least one individual to date are variants of uncertain significance (VUSs). The designation of VUS is a barrier to the use of sequence data in clinical care. We selected 83 variants: 10 previously studied control variants, 10 suspected benign variants, and 63 suspected Brugada syndrome-associated variants, selected on the basis of their frequency in the general population and in individuals with Brugada syndrome. We used high-throughput automated patch clamping to study the function of the 83 variants, with the goal of reclassifying variants with functional data. The ten previously studied controls had functional properties concordant with published manual patch clamp data. All 10 suspected benign variants had wild-type-like function. 22 suspected BrS variants had loss of channel function (<10% normalized peak current) and 22 variants had partial loss of function (10%-50% normalized peak current). The previously unstudied variants were initially classified as likely benign (n = 2), likely pathogenic (n = 10), or VUSs (n = 61). After the patch clamp studies, 16 variants were benign/likely benign, 45 were pathogenic/likely pathogenic, and only 12 were still VUSs. Structural modeling identified likely mechanisms for loss of function including altered thermostability and disruptions to alpha helices, disulfide bonds, or the permeation pore. High-throughput patch clamping enabled reclassification of the majority of tested VUSs in SCN5A.

High-throughput discovery of trafficking-deficient variants in the cardiac potassium channel KV11.1

BACKGROUND

KCHN2 encodes the K_V11.1 potassium channel responsible for I_Kr, a major repolarization current during the cardiomyocyte action potential. Variants in KCNH2 that lead to decreased I_Kr have been associated with long QT syndrome type 2 (LQT2). The mechanism of LQT2 is most often induced loss of K_V11.1 trafficking to the cell surface. Accurately discriminating between variants with normal and abnormal trafficking would aid in understanding the deleterious nature of these variants; however, the volume of reported nonsynonymous KCNH2 variants precludes the use of conventional methods for functional study.

OBJECTIVE

The purpose of this study was to report a high-throughput, multiplexed screening method for KCNH2 genetic variants capable of measuring the cell surface abundance of hundreds of missense variants in the resulting K_V11.1 channel.

METHODS

We developed a method to quantitate K_V11.1 variant trafficking on a pilot region of 11 residues in the S5 helix.

RESULTS

We generated trafficking scores for 220 of 231 missense variants in the pilot region. For 5 of 5 variants, high-throughput trafficking scores validated when tested in single variant flow cytometry and confocal microscopy experiments. We further explored these results with planar patch electrophysiology and found that loss-of-trafficking variants do not produce I_Kr. Conversely, but expectedly, some variants that traffic normally were still functionally compromised.

CONCLUSION

We describe a new method for detecting K_V11.1 trafficking-deficient variants in a multiplexed assay. This new method accurately generated trafficking data for variants in K_V11.1 and is extendable both to all residues in K_V11.1 and to other cell surface proteins.

A Bayesian method to estimate variant-induced disease penetrance

A major challenge emerging in genomic medicine is how to assess best disease risk from rare or novel variants found in disease-related genes. The expanding volume of data generated by very large phenotyping efforts coupled to DNA sequence data presents an opportunity to reinterpret genetic liability of disease risk. Here we propose a framework to estimate the probability of disease given the presence of a genetic variant conditioned on features of that variant. We refer to this as the penetrance, the fraction of all variant heterozygotes that will present with disease. We demonstrate this methodology using a well-established disease-gene pair, the cardiac sodium channel gene SCN5A and the heart arrhythmia Brugada syndrome. From a review of 756 publications, we developed a pattern mixture algorithm, based on a Bayesian Beta-Binomial model, to generate SCN5A penetrance probabilities for the Brugada syndrome conditioned on variant-specific attributes. These probabilities are determined from variant-specific features (e.g. function, structural context, and sequence conservation) and from observations of affected and unaffected heterozygotes. Variant functional perturbation and structural context prove most predictive of Brugada syndrome penetrance.

Deep Mutational Scan of an SCN5A Voltage Sensor

BACKGROUND

Variants in ion channel genes have classically been studied in low throughput by patch clamping. Deep mutational scanning is a complementary approach that can simultaneously assess function of thousands of variants.

METHODS

We have developed and validated a method to perform a deep mutational scan of variants in SCN5A, which encodes the major voltage-gated sodium channel in the heart. We created a library of nearly all possible variants in a 36 base region of SCN5A in the S4 voltage sensor of domain IV and stably integrated the library into HEK293T cells.

RESULTS

In preliminary experiments, challenge with 3 drugs (veratridine, brevetoxin, and ouabain) could discriminate wild-type channels from gain- and loss-of-function pathogenic variants. High-throughput sequencing of the pre- and postdrug challenge pools was used to count the prevalence of each variant and identify variants with abnormal function. The deep mutational scan scores identified 40 putative gain-of-function and 33 putative loss-of-function variants. For 8 of 9 variants, patch clamping data were consistent with the scores.

CONCLUSIONS

These experiments demonstrate the accuracy of a high-throughput in vitro scan of SCN5A variant function, which can be used to identify deleterious variants in SCN5A and other ion channel genes.

SCN5A (NaV1.5) Variant Functional Perturbation and Clinical Presentation: Variants of a Certain Significance

BACKGROUND

Accurately predicting the impact of rare nonsynonymous variants on disease risk is an important goal in precision medicine. Variants in the cardiac sodium channel SCN5A (protein Na_V1.5; voltage-dependent cardiac Na+ channel) are associated with multiple arrhythmia disorders, including Brugada syndrome and long QT syndrome. Rare SCN5A variants also occur in ≈1% of unaffected individuals. We hypothesized that in vitro electrophysiological functional parameters explain a statistically significant portion of the variability in disease penetrance.

METHODS

From a comprehensive literature review, we quantified the number of carriers presenting with and without disease for 1712 reported SCN5A variants. For 356 variants, data were also available for 5 Na_V1.5 electrophysiological parameters: peak current, late/persistent current, steady-state V1/2 of activation and inactivation, and recovery from inactivation.

RESULTS

We found that peak and late current significantly associate with Brugada syndrome (P<0.001; ρ=-0.44; Spearman rank test) and long QT syndrome disease penetrance (P<0.001; ρ=0.37). Steady-state V1/2 activation and recovery from inactivation associate significantly with Brugada syndrome and long QT syndrome penetrance, respectively. Continuous estimates of disease penetrance align with the current American College of Medical Genetics classification paradigm.

CONCLUSIONS

Na_V1.5 in vitro electrophysiological parameters are correlated with Brugada syndrome and long QT syndrome disease risk. Our data emphasize the value of in vitro electrophysiological characterization and incorporating counts of affected and unaffected carriers to aid variant classification. This quantitative analysis of the electrophysiological literature should aid the interpretation of Na_V1.5 variant electrophysiological abnormalities and help improve Na_V1.5 variant classification.

Patient-independent human induced pluripotent stem cell model: A new tool for rapid determination of genetic variant pathogenicity in long QT syndrome

BACKGROUND

Commercial genetic testing for long QT syndrome (LQTS) has rapidly expanded, but the inability to accurately predict whether a rare variant is pathogenic has limited its clinical benefit. Novel missense variants are routinely reported as variant of unknown significance (VUS) and cannot be used to screen family members at risk for sudden cardiac death. Better approaches to determine the pathogenicity of VUS are needed.

OBJECTIVE

The purpose of this study was to rapidly determine the pathogenicity of a CACNA1C variant reported by commercial genetic testing as a VUS using a patient-independent human induced pluripotent stem cell (hiPSC) model.

METHODS

Using CRISPR/Cas9 genome editing, CACNA1C-p.N639T was introduced into a previously established hiPSC from an unrelated healthy volunteer, thereby generating a patient-independent hiPSC model. Three independent heterozygous N639T hiPSC lines were generated and differentiated into cardiomyocytes (CM). Electrophysiological properties of N639T hiPSC-CM were compared to those of isogenic and population control hiPSC-CM by measuring the extracellular field potential (EFP) of 96-well hiPSC-CM monolayers and by patch clamp.

RESULTS

Significant EFP prolongation was observed only in optically stimulated but not in spontaneously beating N639T hiPSC-CM. Patch-clamp studies revealed that N639T prolonged the ventricular action potential by slowing voltage-dependent inactivation of Ca_V1.2 currents. Heterologous expression studies confirmed the effect of N639T on Ca_V1.2 inactivation.

CONCLUSION

The patient-independent hiPSC model enabled rapid generation of functional data to support reclassification of a CACNA1C VUS to likely pathogenic, thereby establishing a novel LQTS type 8 mutation. Furthermore, our results indicate the importance of controlling beating rates to evaluate the functional significance of LQTS VUS in high-throughput hiPSC-CM assays.

Protein structure aids predicting functional perturbation of missense variants in SCN5A and KCNQ1

Rare variants in the cardiac potassium channel K_V7.1 (KCNQ1) and sodium channel Na_V1.5 (SCN5A) are implicated in genetic disorders of heart rhythm, including congenital long QT and Brugada syndromes (LQTS, BrS), but also occur in reference populations. We previously reported two sets of Na_V1.5 (n = 356) and K_V7.1 (n = 144) variants with in vitro characterized channel currents gathered from the literature. Here we investigated the ability to predict commonly reported Na_V1.5 and K_V7.1 variant functional perturbations by leveraging diverse features including variant classifiers PROVEAN, PolyPhen-2, and SIFT; evolutionary rate and BLAST position specific scoring matrices (PSSM); and structure-based features including “functional densities” which is a measure of the density of pathogenic variants near the residue of interest. Structure-based functional densities were the most significant features for predicting Na_V1.5 peak current (adj. R² = 0.27) and K_V7.1 + KCNE1 half-maximal voltage of activation (adj. R² = 0.29). Additionally, use of structure-based functional density values improves loss-of-function classification of SCN5A variants with an ROC-AUC of 0.78 compared with other predictive classifiers (AUC = 0.69; two-sided DeLong test p = .01). These results suggest structural data can inform predictions of the effect of uncharacterized SCN5A and KCNQ1 variants to provide a deeper understanding of their burden on carriers.

Multiple mechanisms underlie increased cardiac late sodium current

BACKGROUND

We recently reported a quantitative relationship between the degree of functional perturbation reported in the literature for 356 variants in the cardiac sodium channel gene SCN5A and the penetrance of resulting arrhythmia phenotypes. In the course of that work, we identified multiple SCN5A variants, including R1193Q, that are common in populations but are reported in human embryonic kidney (HEK) cells to generate large late sodium current (I_Na-L).

OBJECTIVE

The purpose of this study was to compare the functional properties of R1193Q with those of the well-studied type 3 long QT syndrome mutation ΔKPQ.

METHODS

We compared functional properties of SCN5A R1193Q with those of ΔKPQ in Chinese hamster ovary (CHO) cells at baseline and after exposure to intracellular phosphatidylinositol (3,4,5)-trisphosphate (PIP₃), which inhibits I_Na-L generated by decreased Phosphoinositide 3-kinase (PI3K) activity. We also used CRISPR/Cas9 editing to generate R1193Q in human-induced pluripotent stem cells differentiated to cardiomyocytes (hiPSC-CMs).

RESULTS

Both R1193Q and ΔKPQ generated robust I_Na-L in CHO cells. PIP₃ abrogated the late current phenotype in R1193Q cells but had no effect on ΔKPQ. Homozygous R1193Q hiPSC-CMs displayed increased I_Na-L and long action potentials with frequent triggered beats, which were reversed with the addition of PIP₃.

CONCLUSION

The consistency between the late current produced in HEK cells, CHO cells, and hiPSC-CMs suggests that the late current is a feature of the SCN5A R1193Q variant in human cardiomyocytes but that the mechanism by which the late current is produced is distinct and indirect, as compared with the more highly penetrant ΔKPQ. These data suggest that observing a late current in an in vitro setting does not necessarily translate to highly pathogenic type 3 long QT syndrome phenotype but depends on the underlying mechanism.

Predicting the Functional Impact of KCNQ1 Variants of Unknown Significance

BACKGROUND

An emerging standard-of-care for long-QT syndrome uses clinical genetic testing to identify genetic variants of the KCNQ1 potassium channel. However, interpreting results from genetic testing is confounded by the presence of variants of unknown significance for which there is inadequate evidence of pathogenicity.

METHODS AND RESULTS

In this study, we curated from the literature a high-quality set of 107 functionally characterized KCNQ1 variants. Based on this data set, we completed a detailed quantitative analysis on the sequence conservation patterns of subdomains of KCNQ1 and the distribution of pathogenic variants therein. We found that conserved subdomains generally are critical for channel function and are enriched with dysfunctional variants. Using this experimentally validated data set, we trained a neural network, designated Q1VarPred, specifically for predicting the functional impact of KCNQ1 variants of unknown significance. The estimated predictive performance of Q1VarPred in terms of Matthew's correlation coefficient and area under the receiver operating characteristic curve were 0.581 and 0.884, respectively, superior to the performance of 8 previous methods tested in parallel. Q1VarPred is publicly available as a web server at http://meilerlab.org/q1varpred.

CONCLUSIONS

Although a plethora of tools are available for making pathogenicity predictions over a genome-wide scale, previous tools fail to perform in a robust manner when applied to KCNQ1. The contrasting and favorable results for Q1VarPred suggest a promising approach, where a machine-learning algorithm is tailored to a specific protein target and trained with a functionally validated data set to calibrate informatics tools.

A Mechanism of Calmodulin Modulation of the Human Cardiac Sodium Channel

The function of the human cardiac sodium channel (Na_V1.5) is modulated by the Ca²⁺ sensor calmodulin (CaM), but the underlying mechanism(s) are controversial and poorly defined. CaM has been reported to bind in a Ca²⁺-dependent manner to two sites in the intracellular loop that is critical for inactivation of Na_V1.5 (inactivation gate [IG]). The affinity of CaM for the complete IG was significantly stronger than that of fragments that lacked both complete binding sites. Structural analysis by nuclear magnetic resonance, crystallographic, and scattering approaches revealed that CaM simultaneously engages both IG sites using an extended configuration. Patch-clamp recordings for wild-type and mutant channels with an impaired CaM-IG interaction revealed CaM binding to the IG promotes recovery from inactivation while impeding the kinetics of inactivation. Models of full-length Na_V1.5 suggest that CaM binding to the IG directly modulates channel function by destabilizing the inactivated state, which would promote resetting of the IG after channels close.

Structural and biochemical differences between the Notch and the amyloid precursor protein transmembrane domains

γ-Secretase cleavage of the Notch receptor transmembrane domain is a critical signaling event for various cellular processes. Efforts to develop inhibitors of γ-secretase cleavage of the amyloid-β precursor C99 protein as potential Alzheimer's disease therapeutics have been confounded by toxicity resulting from the inhibition of normal cleavage of Notch. We present biochemical and structural data for the combined transmembrane and juxtamembrane Notch domains (Notch-TMD) that illuminate Notch signaling and that can be compared and contrasted with the corresponding traits of C99. The Notch-TMD and C99 have very different conformations, adapt differently to changes in model membrane hydrophobic span, and exhibit different cholesterol-binding properties. These differences may be exploited in the design of agents that inhibit cleavage of C99 while allowing Notch cleavage.

Structural basis for KCNE3 modulation of potassium recycling in epithelia

The single-span membrane protein KCNE3 modulates a variety of voltage-gated ion channels in diverse biological contexts. In epithelial cells, KCNE3 regulates the function of the KCNQ1 potassium ion (K(+)) channel to enable K(+) recycling coupled to transepithelial chloride ion (Cl(-)) secretion, a physiologically critical cellular transport process in various organs and whose malfunction causes diseases, such as cystic fibrosis (CF), cholera, and pulmonary edema. Structural, computational, biochemical, and electrophysiological studies lead to an atomically explicit integrative structural model of the KCNE3-KCNQ1 complex that explains how KCNE3 induces the constitutive activation of KCNQ1 channel activity, a crucial component in K(+) recycling. Central to this mechanism are direct interactions of KCNE3 residues at both ends of its transmembrane domain with residues on the intra- and extracellular ends of the KCNQ1 voltage-sensing domain S4 helix. These interactions appear to stabilize the activated "up" state configuration of S4, a prerequisite for full opening of the KCNQ1 channel gate. In addition, the integrative structural model was used to guide electrophysiological studies that illuminate the molecular basis for how estrogen exacerbates CF lung disease in female patients, a phenomenon known as the "CF gender gap."

Documentation of an Imperative To Improve Methods for Predicting Membrane Protein Stability

There is a compelling and growing need to accurately predict the impact of amino acid mutations on protein stability for problems in personalized medicine and other applications. Here the ability of 10 computational tools to accurately predict mutation-induced perturbation of folding stability (ΔΔG) for membrane proteins of known structure was assessed. All methods for predicting ΔΔG values performed significantly worse when applied to membrane proteins than when applied to soluble proteins, yielding estimated concordance, Pearson, and Spearman correlation coefficients of <0.4 for membrane proteins. Rosetta and PROVEAN showed a modest ability to classify mutations as destabilizing (ΔΔG < −0.5 kcal/mol), with a 7 in 10 chance of correctly discriminating a randomly chosen destabilizing variant from a randomly chosen stabilizing variant. However, even this performance is significantly worse than for soluble proteins. This study highlights the need for further development of reliable and reproducible methods for predicting thermodynamic folding stability in membrane proteins.

Probing Structural Dynamics and Topology of the KCNE1 Membrane Protein in Lipid Bilayers via Site-Directed Spin Labeling and Electron Paramagnetic Resonance Spectroscopy

KCNE1 is a single transmembrane protein that modulates the function of voltage-gated potassium channels, including KCNQ1. Hereditary mutations in the genes encoding either protein can result in diseases such as congenital deafness, long QT syndrome, ventricular tachyarrhythmia, syncope, and sudden cardiac death. Despite the biological significance of KCNE1, the structure and dynamic properties of its physiologically relevant native membrane-bound state are not fully understood. In this study, the structural dynamics and topology of KCNE1 in bilayered lipid vesicles was investigated using site-directed spin labeling (SDSL) and electron paramagnetic resonance (EPR) spectroscopy. A 53-residue nitroxide EPR scan of the KCNE1 protein sequence including all 27 residues of the transmembrane domain (45-71) and 26 residues of the N- and C-termini of KCNE1 in lipid bilayered vesicles was analyzed in terms of nitroxide side-chain motion. Continuous wave-EPR spectral line shape analysis indicated the nitroxide spin label side-chains located in the KCNE1 TMD are less mobile when compared to the extracellular region of KCNE1. The EPR data also revealed that the C-terminus of KCNE1 is more mobile when compared to the N-terminus. EPR power saturation experiments were performed on 41 sites including 18 residues previously proposed to reside in the transmembrane domain (TMD) and 23 residues of the N- and C-termini to determine the topology of KCNE1 with respect to the 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC)/1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-(1'-rac-glycerol) (POPG) lipid bilayers. The results indicated that the transmembrane domain is indeed buried within the membrane, spanning the width of the lipid bilayer. Power saturation data also revealed that the extracellular region of KCNE1 is solvent-exposed with some of the portions partially or weakly interacting with the membrane surface. These results are consistent with the previously published solution NMR structure of KCNE1 in micelles.

Mapping membrane protein backbone dynamics: a comparison of site-directed spin labeling with NMR 15N-relaxation measurements

The ability to detect nanosecond backbone dynamics with site-directed spin labeling (SDSL) in soluble proteins has been well established. However, for membrane proteins, the nitroxide appears to have more interactions with the protein surface, potentially hindering the sensitivity to backbone motions. To determine whether membrane protein backbone dynamics could be mapped with SDSL, a nitroxide was introduced at 55 independent sites in a model polytopic membrane protein, TM0026. Electron paramagnetic resonance spectral parameters were compared with NMR (15)N-relaxation data. Sequential scans revealed backbone dynamics with the same trends observed for the R1 relaxation rate, suggesting that nitroxide dynamics remain coupled to the backbone on membrane proteins.

Personalized biochemistry and biophysics

Whole human genome sequencing of individuals is becoming rapid and inexpensive, enabling new strategies for using personal genome information to help diagnose, treat, and even prevent human disorders for which genetic variations are causative or are known to be risk factors. Many of the exploding number of newly discovered genetic variations alter the structure, function, dynamics, stability, and/or interactions of specific proteins and RNA molecules. Accordingly, there are a host of opportunities for biochemists and biophysicists to participate in (1) developing tools to allow accurate and sometimes medically actionable assessment of the potential pathogenicity of individual variations and (2) establishing the mechanistic linkage between pathogenic variations and their physiological consequences, providing a rational basis for treatment or preventive care. In this review, we provide an overview of these opportunities and their associated challenges in light of the current status of genomic science and personalized medicine, the latter often termed precision medicine.

Structural investigation of the transmembrane domain of KCNE1 in proteoliposomes

KCNE1 is a single-transmembrane protein of the KCNE family that modulates the function of voltage-gated potassium channels, including KCNQ1. Hereditary mutations in KCNE1 have been linked to diseases such as long QT syndrome (LQTS), atrial fibrillation, sudden infant death syndrome, and deafness. The transmembrane domain (TMD) of KCNE1 plays a key role in mediating the physical association with KCNQ1 and in subsequent modulation of channel gating kinetics and conductance. However, the mechanisms associated with these roles for the TMD remain poorly understood, highlighting a need for experimental structural studies. A previous solution NMR study of KCNE1 in LMPG micelles revealed a curved transmembrane domain, a structural feature proposed to be critical to KCNE1 function. However, this curvature potentially reflects an artifact of working in detergent micelles. Double electron electron resonance (DEER) measurements were conducted on KCNE1 in LMPG micelles, POPC/POPG proteoliposomes, and POPC/POPG lipodisq nanoparticles to directly compare the structure of the TMD in a variety of different membrane environments. Experimentally derived DEER distances coupled with simulated annealing molecular dynamic simulations were used to probe the bilayer structure of the TMD of KCNE1. The results indicate that the structure is helical in proteoliposomes and is slightly curved, which is consistent with the previously determined solution NMR structure in micelles. The evident resilience of the curvature in the KCNE1 TMD leads us to hypothesize that the curvature is likely to be maintained upon binding of the protein to the KCNQ1 channel.

The homology model of PMP22 suggests mutations resulting in peripheral neuropathy disrupt transmembrane helix packing

Peripheral myelin protein 22 (PMP22) is a tetraspan membrane protein strongly expressed in myelinating Schwann cells of the peripheral nervous system. Myriad missense mutations in PMP22 result in varying degrees of peripheral neuropathy. We used Rosetta 3.5 to generate a homology model of PMP22 based on the recently published crystal structure of claudin-15. The model suggests that several mutations known to result in neuropathy act by disrupting transmembrane helix packing interactions. Our model also supports suggestions from previous studies that the first transmembrane helix is not tightly associated with the rest of the helical bundle.

Structure of the Neisserial outer membrane protein Opa₆₀: loop flexibility essential to receptor recognition and bacterial engulfment

The structure and dynamics of Opa proteins, which we report herein, are responsible for the receptor-mediated engulfment of Neisseria gonorrheae or Neisseria meningitidis by human cells and can offer deep understanding into the molecular recognition of pathogen-host receptor interactions. Such interactions are vital to understanding bacterial pathogenesis as well as the mechanism of foreign body entry to a human cell, which may provide insights for the development of targeted pharmaceutical delivery systems. The size and dynamics of the extracellular loops of Opa60 required a hybrid refinement approach wherein membrane and distance restraints were used to generate an initial NMR structural ensemble, which was then further refined using molecular dynamics in a DMPC bilayer. The resulting ensemble revealed that the extracellular loops, which bind host receptors, occupy compact conformations, interact with each other weakly, and are dynamic on the nanosecond time scale. We predict that this conformational sampling is critical for enabling diverse Opa loop sequences to engage a common set of receptors.

Purification and structural study of the voltage-sensor domain of the human KCNQ1 potassium ion channel

KCNQ1 (also known as KV7.1 or KVLQT1) is a voltage-gated potassium channel modulated by members of the KCNE protein family. Among multiple functions, KCNQ1 plays a critical role in the cardiac action potential. This channel is also subject to inherited mutations that cause certain cardiac arrhythmias and deafness. In this study, we report the overexpression, purification, and preliminary structural characterization of the voltage-sensor domain (VSD) of human KCNQ1 (Q1-VSD). Q1-VSD was expressed in Escherichia coli and purified into lyso-palmitoylphosphatidylglycerol micelles, conditions under which this tetraspan membrane protein yields excellent nuclear magnetic resonance (NMR) spectra. NMR studies reveal that Q1-VSD shares a common overall topology with other channel VSDs, with an S0 helix followed by transmembrane helices S1-S4. The exact sequential locations of the helical spans do, however, show significant variations from those of the homologous segments of previously characterized VSDs. The S4 segment of Q1-VSD was seen to be α-helical (with no 310 component) and underwent rapid backbone amide H-D exchange over most of its length. These results lay the foundation for more advanced structural studies and can be used to generate testable hypotheses for future structure-function experiments.

Backbone ¹H, ¹³C and ¹⁵N resonance assignments of the α-helical membrane protein TM0026 from Thermotoga maritima

Critical to the use of solution NMR to describe the structure and flexibility of membrane proteins is the thorough understanding of the degree of perturbation induced by the detergent or other membrane mimetic. To develop a deeper understanding of the interaction between membrane proteins and micelles or bicelles, we will investigate the differences in structure and flexibility of a model membrane protein TM0026 from Thermotoga maritima using solution NMR. A comparison of the structural differences between TM0026 solubilized in different detergent combinations will provide important insight into the degree of modulation of membrane proteins by detergent physical properties. Here we report the nearly complete backbone and Cβ resonance assignments of the two transmembrane helical model protein TM0026. These assignments are the first step to using TM0026 to elucidate the interaction between membrane proteins and membrane mimetics.

DEER EPR measurements for membrane protein structures via bifunctional spin labels and lipodisq nanoparticles

Pulsed EPR DEER structural studies of membrane proteins in a lipid bilayer have often been hindered by difficulties in extracting accurate distances when compared to those of globular proteins. In this study, we employed a combination of three recently developed methodologies, (1) bifunctional spin labels (BSL), (2) SMA-Lipodisq nanoparticles, and (3) Q band pulsed EPR measurements, to obtain improved signal sensitivity, increased transverse relaxation time, and more accurate and precise distances in DEER measurements on the integral membrane protein KCNE1. The KCNE1 EPR data indicated an ∼2-fold increase in the transverse relaxation time for the SMA-Lipodisq nanoparticles when compared to those of proteoliposomes and narrower distance distributions for the BSL when compared to those of the standard MTSL. The certainty of information content in DEER data obtained for KCNE1 in SMA-Lipodisq nanoparticles is comparable to that in micelles. The combination of techniques will enable researchers to potentially obtain more precise distances in cases where the traditional spin labels and membrane systems yield imprecise distance distributions.

Reversible folding of human peripheral myelin protein 22, a tetraspan membrane protein

Misfolding of the α-helical membrane protein peripheral myelin protein 22 (PMP22) has been implicated in the pathogenesis of the common neurodegenerative disease known as Charcot-Marie-Tooth disease (CMTD) and also several other related peripheral neuropathies. Emerging evidence suggests that the propensity of PMP22 to misfold in the cell may be due to an intrinsic lack of conformational stability. Therefore, quantitative studies of the conformational equilibrium of PMP22 are needed to gain insight into the molecular basis of CMTD. In this work, we have investigated the folding and unfolding of wild type (WT) human PMP22 in mixed micelles. Both kinetic and thermodynamic measurements demonstrate that the denaturation of PMP22 by n-lauroyl sarcosine (LS) in dodecylphosphocholine (DPC) micelles is reversible. Assessment of the conformational equilibrium indicates that a significant fraction of unfolded PMP22 persists even in the absence of the denaturing detergent. However, we find the stability of PMP22 is increased by glycerol, which facilitates quantitation of thermodynamic parameters. To our knowledge, this work represents the first report of reversible unfolding of a eukaryotic multispan membrane protein. The results indicate that WT PMP22 possesses minimal conformational stability in micelles, which parallels its poor folding efficiency in the endoplasmic reticulum. Folding equilibrium measurements for PMP22 in micelles may provide an approach to assess the effects of cellular metabolites or potential therapeutic agents on its stability. Furthermore, these results pave the way for future investigation of the effects of pathogenic mutations on the conformational equilibrium of PMP22.

MAPK phosphorylation of connexin 43 promotes binding of cyclin E and smooth muscle cell proliferation

RATIONALE

Dedifferentiation of vascular smooth muscle cells (VSMC) leading to a proliferative cell phenotype significantly contributes to the development of atherosclerosis. Mitogen-activated protein kinase (MAPK) phosphorylation of proteins including connexin 43 (Cx43) has been associated with VSMC proliferation in atherosclerosis.

OBJECTIVE

To investigate whether MAPK phosphorylation of Cx43 is directly involved in VSMC proliferation.

METHODS AND RESULTS

We show in vivo that MAPK-phosphorylated Cx43 forms complexes with the cell cycle control proteins cyclin E and cyclin-dependent kinase 2 (CDK2) in carotids of apolipoprotein-E receptor null (ApoE(-/-)) mice and in C57Bl/6 mice treated with platelet-derived growth factor-BB (PDGF). We tested the involvement of Cx43 MAPK phosphorylation in vitro using constructs for full-length Cx43 (Cx43) or the Cx43 C-terminus (Cx43(CT)) and produced null phosphorylation Ser>Ala (Cx43(MK4A)/Cx43(CTMK4A)) and phospho-mimetic Ser>Asp (Cx43(MK4D)/Cx43(CTMK4D)) mutations. Coimmunoprecipitation studies in primary VSMC isolated from Cx43 wild-type (Cx43(+/+)) and Cx43 null (Cx43(-/-)) mice and analytic size exclusion studies of purified proteins identify that interactions between cyclin E and Cx43 requires Cx43 MAPK phosphorylation. We further demonstrate that Cx43 MAPK phosphorylation is required for PDGF-mediated VSMC proliferation. Finally, using a novel knock-in mouse containing Cx43-MK4A mutation, we show in vivo that interactions between Cx43 and cyclin E are lost and VSMC proliferation does not occur after treatment of carotids with PDGF and that neointima formation is significantly reduced in carotids after injury.

CONCLUSIONS

We identify MAPK-phosphorylated Cx43 as a novel interacting partner of cyclin E in VSMC and show that this interaction is critical for VSMC proliferation. This novel interaction may be important in the development of atherosclerotic lesions.

Identification and removal of nitroxide spin label contaminant: impact on PRE studies of α-helical membrane proteins in detergent

NMR paramagnetic relaxation enhancement (PRE) provides long-range distance constraints (~15-25 Å) that can be critical to determining overall protein topology, especially where long-range NOE information is unavailable such as in the case of larger proteins that require deuteration. However, several challenges currently limit the use of NMR PRE for α-helical membrane proteins. One challenge is the nonspecific association of the nitroxide spin label to the protein-detergent complex that can result in spurious PRE derived distance restraints. The effect of the nitroxide spin label contaminant is evaluated and quantified and a robust method for the removal of the contaminant is provided to advance the application of PRE restraints to membrane protein NMR structure determination.

Structural origins of nitroxide side chain dynamics on membrane protein α-helical sites

Understanding the structure and dynamics of membrane proteins in their native, hydrophobic environment is important to understanding how these proteins function. EPR spectroscopy in combination with site-directed spin labeling (SDSL) can measure dynamics and structure of membrane proteins in their native lipid environment; however, until now the dynamics measured have been qualitative due to limited knowledge of the nitroxide spin label's intramolecular motion in the hydrophobic environment. Although several studies have elucidated the structural origins of EPR line shapes of water-soluble proteins, EPR spectra of nitroxide spin-labeled proteins in detergents or lipids have characteristic differences from their water-soluble counterparts, suggesting significant differences in the underlying molecular motion of the spin label between the two environments. To elucidate these differences, membrane-exposed α-helical sites of the leucine transporter, LeuT, from Aquifex aeolicus, were investigated using X-ray crystallography, mutational analysis, nitroxide side chain derivatives, and spectral simulations in order to obtain a motional model of the nitroxide. For each crystal structure, the nitroxide ring of a disulfide-linked spin label side chain (R1) is resolved and makes contacts with hydrophobic residues on the protein surface. The spin label at site I204 on LeuT makes a nontraditional hydrogen bond with the ortho-hydrogen on its nearest neighbor F208, whereas the spin label at site F177 makes multiple van der Waals contacts with a hydrophobic pocket formed with an adjacent helix. These results coupled with the spectral effect of mutating the i ± 3, 4 residues suggest that the spin label has a greater affinity for its local protein environment in the low dielectric than on a water-soluble protein surface. The simulations of the EPR spectra presented here suggest the spin label oscillates about the terminal bond nearest the ring while maintaining weak contact with the protein surface. Combined, the results provide a starting point for determining a motional model for R1 on membrane proteins, allowing quantification of nitroxide dynamics in the aliphatic environment of detergent and lipids. In addition, initial contributions to a rotamer library of R1 on membrane proteins are provided, which will assist in reliably modeling the R1 conformational space for pulsed dipolar EPR and NMR paramagnetic relaxation enhancement distance determination.