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

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.

Predicting the functional impact of KCNQ1 variants with artificial neural networks

Recent advances in experimental and computational protein structure determination have provided access to high-quality structures for most human proteins and mutants thereof. However, linking changes in structure in protein mutants to functional impact remains an active area of method development. If successful, such methods can ultimately assist physicians in taking appropriate treatment decisions. This work presents three artificial neural network (ANN)-based predictive models that classify four key functional parameters of KCNQ1 variants as normal or dysfunctional using PSSM-based evolutionary and/or biophysical descriptors. Recent advances in predicting protein structure and variant properties with artificial intelligence (AI) rely heavily on the availability of evolutionary features and thus fail to directly assess the biophysical underpinnings of a change in structure and/or function. The central goal of this work was to develop an ANN model based on structure and physiochemical properties of KCNQ1 potassium channels that performs comparably or better than algorithms using only on PSSM-based evolutionary features. These biophysical features highlight the structure-function relationships that govern protein stability, function, and regulation. The input sensitivity algorithm incorporates the roles of hydrophobicity, polarizability, and functional densities on key functional parameters of the KCNQ1 channel. Inclusion of the biophysical features outperforms exclusive use of PSSM-based evolutionary features in predicting activation voltage dependence and deactivation time. As AI is increasingly applied to problems in biology, biophysical understanding will be critical with respect to 'explainable AI', i.e., understanding the relation of sequence, structure, and function of proteins. Our model is available at www.kcnq1predict.org.

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.

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.

Structures Illuminate Cardiac Ion Channel Functions in Health and in Long QT Syndrome

The cardiac action potential is critical to the production of a synchronized heartbeat. This electrical impulse is governed by the intricate activity of cardiac ion channels, among them the cardiac voltage-gated potassium (Kv) channels KCNQ1 and hERG as well as the voltage-gated sodium (Nav) channel encoded by SCN5A. Each channel performs a highly distinct function, despite sharing a common topology and structural components. These three channels are also the primary proteins mutated in congenital long QT syndrome (LQTS), a genetic condition that predisposes to cardiac arrhythmia and sudden cardiac death due to impaired repolarization of the action potential and has a particular proclivity for reentrant ventricular arrhythmias. Recent cryo-electron microscopy structures of human KCNQ1 and hERG, along with the rat homolog of SCN5A and other mammalian sodium channels, provide atomic-level insight into the structure and function of these proteins that advance our understanding of their distinct functions in the cardiac action potential, as well as the molecular basis of LQTS. In this review, the gating, regulation, LQTS mechanisms, and pharmacological properties of KCNQ1, hERG, and SCN5A are discussed in light of these recent structural findings.

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.

Genotype-Phenotype Correlation in a Family with Brugada Syndrome Harboring the Novel p.Gln371* Nonsense Variant in the SCN5A Gene

Brugada syndrome (BrS) is marked by coved ST-segment elevation and increased risk of sudden cardiac death. The genetics of this syndrome are elusive in over half of the cases. Variants in the SCN5A gene are the single most common known genetic unifier, accounting for about a third of cases. Research models, such as animal models and cell lines, are limited. In the present study, we report the novel NM_198056.2:c.1111C>T (p.Gln371*) heterozygous variant in the SCN5A gene, as well as its segregation with BrS in a large family. The results herein suggest a pathogenic effect of this variant. Functional studies are certainly warranted to characterize the molecular effects of this variant.

Computational identification of vesicular transport proteins from sequences using deep gated recurrent units architecture

Protein function prediction is one of the most well-studied topics, attracting attention from countless researchers in the field of computational biology. Implementing deep neural networks that help improve the prediction of protein function, however, is still a major challenge. In this research, we suggested a new strategy that includes gated recurrent units and position-specific scoring matrix profiles to predict vesicular transportation proteins, a biological function of great importance. Although it is difficult to discover its function, our model is able to achieve accuracies of 82.3% and 85.8% in the cross-validation and independent dataset, respectively. We also solve the problem of imbalance in the dataset via tuning class weight in the deep learning model. The results generated showed sensitivity, specificity, MCC, and AUC to have values of 79.2%, 82.9%, 0.52, and 0.861, respectively. Our strategy shows superiority in results on the same dataset against all other state-of-the-art algorithms. In our suggested research, we have suggested a technique for the discovery of more proteins, particularly proteins connected with vesicular transport. In addition, our accomplishment could encourage the use of gated recurrent units architecture in protein function prediction.