Using an integrative machine learning approach utilising homology modelling to clinically interpret genetic variants: CACNA1F as an exemplar

Truncated variants of thyroid hormone receptor beta display disease-inflicting malfunctioning at cellular level

Thyroid hormone receptor β (THRβ) is a member of the nuclear receptor superfamily of ligand-modulated transcription factors. Upon ligand binding, THRβ sequentially recruits the components of transcriptional machinery to modulate target gene expression. In addition to regulating diverse physiological processes, THRβ plays a crucial role in hypothalamus-pituitary-thyroid axis feedback regulation. Anomalies in THRβ gene/protein structure are associated with onset of diverse disease states. In this study, we investigated disease-inflicting truncated variants of THRβ using in-silico analysis and cell-based assays. We examined the THRβ truncated variants on multiple test parameters, including subcellular localization, ligand-receptor interactions, transcriptional functions, interaction with heterodimeric partner RXR, and receptor-chromatin interactions. Moreover, molecular dynamic simulation approaches predicted that shortened THRβ-LBD due to point mutations contributes proportionally to the loss of structural integrity and receptor stability. Deviant subcellular localization and compromised transcriptional function were apparent with these truncated variants. Present study shows that 'mitotic bookmarking' property of some THRβ variants is also affected. The study highlights that structural and conformational attributes of THRβ are necessary for normal receptor functioning, and any deviations may contribute to the underlying cause of the inflicted diseases. We anticipate that insights derived herein may contribute to improved mechanistic understanding to assess disease predisposition.

Protein modeling and in silico analysis to assess pathogenicity of ABCA4 variants in patients with inherited retinal disease

Purpose

The retina-specific ABCA transporter, ABCA4, plays an essential role in translocating retinoids required by the visual cycle. ABCA4 genetic variants are known to cause a wide range of inherited retinal disorders, including Stargardt disease and cone-rod dystrophy. More than 1,400 ABCA4 missense variants have been identified; however, more than half of these remain variants of uncertain significance (VUS). The purpose of this study was to employ a predictive strategy to assess the pathogenicity of ABCA4 variants in inherited retinal diseases using protein modeling and computational approaches.

Methods

We studied 13 clinically well-defined patients with ABCA4 retinopathies and identified the presence of 10 missense variants, including one novel variant in the ABCA4 gene, by next-generation sequencing (NGS). All variants were structurally analyzed using AlphaFold2 models and existing experimental structures of human ABCA4 protein. The results of these analyses were compared with patient clinical presentations to test the effectiveness of the methods employed in predicting variant pathogenicity.

Results

We conducted a phenotype-genotype comparison of 13 genetically and phenotypically well-defined retinal disease patients. The in silico protein structure analyses we employed successfully detected the deleterious effect of missense variants found in this affected patient cohort. Our study provides American College of Medical Genetics and Genomics (ACMG)-defined supporting evidence of the pathogenicity of nine missense ABCA4 variants, aligning with the observed clinical phenotypes in this cohort.

Conclusions

In this report, we describe a systematic approach to predicting the pathogenicity of ABCA4 variants by means of three-dimensional (3D) protein modeling and in silico structure analysis. Our results demonstrate concordance between disease severity and structural changes in protein models induced by genetic variations. Furthermore, the present study suggests that in silico protein structure analysis can be used as a predictor of pathogenicity and may facilitate the assessment of genetic VUS.

In silico versus functional characterization of genetic variants: lessons from muscle channelopathies

Accurate determination of the pathogenicity of missense genetic variants of uncertain significance is a huge challenge for implementing genetic data in clinical practice. In silico predictive tools are used to score variants' pathogenicity. However, their value in clinical settings is often unclear, as they have not usually been validated against robust functional assays. We compared nine widely used in silico predictive tools, including more recently developed tools (EVE and REVEL) with detailed cell-based electrophysiology, for 126 CLCN1 variants discovered in patients with the skeletal muscle channelopathy myotonia congenita. We found poor accuracy for most tools. The highest accuracy was obtained with MutationTaster (84.58%) and REVEL (82.54%). Both of these scores showed poor specificity, although specificity was better using EVE. Combining methods based on concordance improved performance overall but still lacked specificity. Our calculated statistics for the predictive tools were different to reported values for other genes in the literature, suggesting that the utility of the tools varies between genes. Overall, current predictive tools for this chloride channel are not reliable for clinical use, and tools with better specificity are urgently required. Improving the accuracy of predictive tools is a wider issue and a huge challenge for effective clinical implementation of genetic data.

Cav1.4 congenital stationary night blindness is associated with an increased rate of proteasomal degradation

Pathogenic, generally loss-of-function, variants in CACNA1F, encoding the Ca_v1.4α₁ calcium channel, underlie congenital stationary night blindness type 2 (CSNB2), a rare inherited retinal disorder associated with visual disability. To establish the underlying pathomechanism, we investigated 10 clinically derived CACNA1F missense variants located across pore-forming domains, connecting loops, and the carboxy-tail domain of the Ca_v1.4α subunit. Homology modeling showed that all variants cause steric clashes; informatics analysis correctly predicted pathogenicity for 7/10 variants. In vitro analyses demonstrated that all variants cause a decrease in current, global expression, and protein stability and act through a loss-of-function mechanism and suggested that the mutant Ca_v1.4α proteins were degraded by the proteasome. We showed that the reduced current for these variants could be significantly increased through treatment with clinical proteasome inhibitors. In addition to facilitating clinical interpretation, these studies suggest that proteasomal inhibition represents an avenue of potential therapeutic intervention for CSNB2.

Assessing the Pathogenicity of In-Frame CACNA1F Indel Variants Using Structural Modeling

Small in-frame insertion-deletion (indel) variants are a common form of genomic variation whose impact on rare disease phenotypes has been understudied. The prediction of the pathogenicity of such variants remains challenging. X-linked incomplete congenital stationary night blindness type 2 (CSNB2) is a nonprogressive, inherited retinal disorder caused by variants in CACNA1F, encoding the Ca_v1.4α1 channel protein. Here, structural analysis was used through homology modeling to interpret 10 disease-correlated and 10 putatively benign CACNA1F in-frame indel variants. CSNB2-correlated changes were found to be more highly conserved compared with putative benign variants. Notably, all 10 disease-correlated variants but none of the benign changes were within modeled regions of the protein. Structural analysis revealed that disease-correlated variants are predicted to destabilize the structure and function of the Ca_v1.4α1 channel protein. Overall, the use of structural information to interpret the consequences of in-frame indel variants provides an important adjunct that can improve the diagnosis for individuals with CSNB2.

Computational approaches for predicting variant impact: An overview from resources, principles to applications

One objective of human genetics is to unveil the variants that contribute to human diseases. With the rapid development and wide use of next-generation sequencing (NGS), massive genomic sequence data have been created, making personal genetic information available. Conventional experimental evidence is critical in establishing the relationship between sequence variants and phenotype but with low efficiency. Due to the lack of comprehensive databases and resources which present clinical and experimental evidence on genotype-phenotype relationship, as well as accumulating variants found from NGS, different computational tools that can predict the impact of the variants on phenotype have been greatly developed to bridge the gap. In this review, we present a brief introduction and discussion about the computational approaches for variant impact prediction. Following an innovative manner, we mainly focus on approaches for non-synonymous variants (nsSNVs) impact prediction and categorize them into six classes. Their underlying rationale and constraints, together with the concerns and remedies raised from comparative studies are discussed. We also present how the predictive approaches employed in different research. Although diverse constraints exist, the computational predictive approaches are indispensable in exploring genotype-phenotype relationship.

Improving the clinical interpretation of missense variants in X linked genes using structural analysis

BACKGROUND

Improving the clinical interpretation of missense variants can increase the diagnostic yield of genomic testing and lead to personalised management strategies. Currently, due to the imprecision of bioinformatic tools that aim to predict variant pathogenicity, their role in clinical guidelines remains limited. There is a clear need for more accurate prediction algorithms and this study aims to improve performance by harnessing structural biology insights. The focus of this work is missense variants in a subset of genes associated with X linked disorders.

METHODS

We have developed a protein-specific variant interpreter (ProSper) that combines genetic and protein structural data. This algorithm predicts missense variant pathogenicity by applying machine learning approaches to the sequence and structural characteristics of variants.

RESULTS

ProSper outperformed seven previously described tools, including meta-predictors, in correctly evaluating whether or not variants are pathogenic; this was the case for 11 of the 21 genes associated with X linked disorders that met the inclusion criteria for this study. We also determined gene-specific pathogenicity thresholds that improved the performance of VEST4, REVEL and ClinPred, the three best-performing tools out of the seven that were evaluated; this was the case in 11, 11 and 12 different genes, respectively.

CONCLUSION

ProSper can form the basis of a molecule-specific prediction tool that can be implemented into diagnostic strategies. It can allow the accurate prioritisation of missense variants associated with X linked disorders, aiding precise and timely diagnosis. In addition, we demonstrate that gene-specific pathogenicity thresholds for a range of missense prioritisation tools can lead to an increase in prediction accuracy.

In Silico Tools and Approaches for the Prediction of Functional and Structural Effects of Single-Nucleotide Polymorphisms on Proteins: An Expert Review

Single-nucleotide polymorphisms (SNPs) are single-base variants that contribute to human biological variation and pathogenesis of many human diseases. Among all SNP types, nonsynonymous single-nucleotide polymorphisms (nsSNPs) can alter many structural, biochemical, and functional features of a protein such as folding characteristics, charge distribution, stability, dynamics, and interactions with other proteins/nucleotides. These modifications in the protein structure can lead nsSNPs to be closely associated with many multifactorial diseases such as cancer, diabetes, and neurodegenerative diseases. Predicting structural and functional effects of nsSNPs with experimental approaches can be time-consuming and costly; hence, computational prediction tools and algorithms are being widely and increasingly utilized in biology and medical research. This expert review examines the in silico tools and algorithms for the prediction of functional or structural effects of SNP variants, in addition to the description of the phenotypic effects of nsSNPs on protein structure, association between pathogenicity of variants, and functional or structural features of disease-associated variants. Finally, case studies investigating the functional and structural effects of nsSNPs on selected protein structures are highlighted. We conclude that creating a consistent workflow with a combination of in silico approaches or tools should be considered to increase the performance, accuracy, and precision of the biological and clinical predictions made in silico.