Bioinformatic Analysis of GJB2 Gene Missense Mutations

Computational analysis and molecular dynamics simulation of high-risk single nucleotide polymorphisms of the ADAM10 gene

Polymorphisms of the disintegrin and metalloproteinase domain-containing protein 10 (ADAM10) are linked to pathophysiological changes in lung inflammation, cancer, Alzheimer's disease (AD), encephalopathy, liver fibrosis, and cardiovascular diseases. In this study, we predicted the pathogenicity of ADAM10 non-synonymous single nucleotide polymorphisms (nsSNPs) in a wide array of mutation analyzing bioinformatics tools. We retrieved 423 nsSNPs from dbSNP-NCBI for the analysis, and 13 were predicted deleterious by each of the ten tools: SIFT, PROVEAN, CONDEL, PANTHER-PSEP, SNAP2, SuSPect, PolyPhen-2, Meta-SNP, Mutation Assessor and Predict-SNP. Further assessment of amino acid sequences, homology models, conservation profiles, and inter-atomic interactions identified C222G, G361E and C639Y as the most pathogenic mutations. We validated this prediction through structural stability analysis using DUET, I-Mutant Suite, SNPeffect and Dynamut. Molecular dynamics simulations and principal component analysis also indicated considerable instability of the C222G, G361E and C639Y variants. Therefore, these ADAM10 nsSNPs could be candidates for diagnostic genetic screening and therapeutic molecular targeting.Communicated by Ramaswamy H. Sarma.

Functional Consequences of Pathogenic Variants of the GJB2 Gene (Cx26) Localized in Different Cx26 Domains

One of the most common forms of genetic deafness has been predominantly associated with pathogenic variants in the GJB2 gene, encoding transmembrane protein connexin 26 (Cx26). The Cx26 molecule consists of an N-terminal domain (NT), four transmembrane domains (TM1-TM4), two extracellular loops (EL1 and EL2), a cytoplasmic loop, and a C-terminus (CT). Pathogenic variants in the GJB2 gene, resulting in amino acid substitutions scattered across the Cx26 domains, lead to a variety of clinical outcomes, including the most common non-syndromic autosomal recessive deafness (DFNB1A), autosomal dominant deafness (DFNA3A), as well as syndromic forms combining hearing loss and skin disorders. However, for rare and poorly documented variants, information on the mode of inheritance is often lacking. Numerous in vitro studies have been conducted to elucidate the functional consequences of pathogenic GJB2 variants leading to amino acid substitutions in different domains of Cx26 protein. In this work, we summarized all available data on a mode of inheritance of pathogenic GJB2 variants leading to amino acid substitutions and reviewed published information on their functional effects, with an emphasis on their localization in certain Cx26 domains.

Comprehensive interpretation of single-nucleotide substitutions in GJB2 reveals the genetic and phenotypic landscape of GJB2-related hearing loss

Genetic variants in GJB2 are the most frequent cause of congenital and childhood hearing loss worldwide. The purpose of this study was to delineate the genetic and phenotypic landscape of GJB2 SNV variants. All possible single-nucleotide substitution variants of the coding region of GJB2 (N = 2043) were manually curated following the ACMG/AMP hearing loss guidelines. As a result, 60 (2.9%), 177 (8.7%), 1499 (73.4%), 301 (14.7%) and 6 (0.3%) of the variants were classified as pathogenic, likely pathogenic, variant of uncertain significance, likely benign, and benign, respectively. 53% (84/158) of the pathogenic/likely pathogenic missense variants were not present in ClinVar. The second transmembrane domain and the 3₁₀ helix were highly enriched for pathogenic missense variants, while the intracellular loops were tolerant to variation. The N-terminal tail and the extracellular loop showed high clustering of variants that are associated with syndromic or dominant non-syndromic hearing loss. In conclusion, our study interpreted all possible single-nucleotide substitution coding variants, characterized novel clinically significant variants in GJB2, and revealed significant genotype-phenotype correlations at this common hearing loss locus. Our work provides a prototype for other genes with similarly high genetic and phenotypic heterogeneity.

MutTMPredictor: Robust and accurate cascade XGBoost classifier for prediction of mutations in transmembrane proteins

Transmembrane proteins have critical biological functions and play a role in a multitude of cellular processes including cell signaling, transport of molecules and ions across membranes. Approximately 60% of transmembrane proteins are considered as drug targets. Missense mutations in such proteins can lead to many diverse diseases and disorders, such as neurodegenerative diseases and cystic fibrosis. However, there are limited studies on mutations in transmembrane proteins. In this work, we first design a new feature encoding method, termed weight attenuation position-specific scoring matrix (WAPSSM), which builds upon the protein evolutionary information. Then, we propose a new mutation prediction algorithm (cascade XGBoost) by leveraging the idea learned from consensus predictors and gcForest. Multi-level experiments illustrate the effectiveness of WAPSSM and cascade XGBoost algorithms. Finally, based on WAPSSM and other three types of features, in combination with the cascade XGBoost algorithm, we develop a new transmembrane protein mutation predictor, named MutTMPredictor. We benchmark the performance of MutTMPredictor against several existing predictors on seven datasets. On the 546 mutations dataset, MutTMPredictor achieves the accuracy (ACC) of 0.9661 and the Matthew's Correlation Coefficient (MCC) of 0.8950. While on the 67,584 dataset, MutTMPredictor achieves an MCC of 0.7523 and area under curve (AUC) of 0.8746, which are 0.1625 and 0.0801 respectively higher than those of the existing best predictor (fathmm). Besides, MutTMPredictor also outperforms two specific predictors on the Pred-MutHTP datasets. The results suggest that MutTMPredictor can be used as an effective method for predicting and prioritizing missense mutations in transmembrane proteins. The MutTMPredictor webserver and datasets are freely accessible at http://csbio.njust.edu.cn/bioinf/muttmpredictor/ for academic use.

Molecular alteration in the Gap Junction Beta 2 (GJB2) gene associated with non-syndromic sensorineural hearing impairment

Non-syndromic sensory neural hearing defect is one of the genetic diseases inherited from parents to offerings. The autosomal recessive form affects a large population worldwide and has become a major concern in the social and professional lives of many people. There are many factors and genes which are involved in hearing loss but the Gap Junction Beta 2 (GJB2) gene which encodes the connexin 26 protein, is a major cause of non-syndromic recessive deafness (NSRD). This study aims to record and analyze GJB2 gene mutations in the hearing-impaired population of North Karnataka, India. In this study, we included 368 congenitally hearing-impaired children from North Karnataka, India, under 18 years of age. After thorough clinical examinations, patient's history and proper audiological results, peripheral blood samples were collected and subjected to genetic analysis. We recorded that 54.8% of the NSRD cases have an autosomal recessive mutation in the coding region of the GJB2 gene. The frequency of W24X (25%) mutation was found to be high in the present study population. From this study we can suggest that, identifying this mutation in new-borns definitely helps in the early diagnosis of hearing loss.

In Silico Prediction of the Effects of Nonsynonymous Single Nucleotide Polymorphisms in the Human Catechol-O-Methyltransferase (COMT) Gene

Catechol-O-methyltransferase (COMT) enzyme performs transfer of methyl group to endogenous and exogenous catechol substrates. The COMT enzyme draws interest because of its association with psychiatric, neurological and cardiovascular diseases, and several cancers. Moreover, many prescribed drugs, supplements, and their metabolites are used as substrates of COMT enzyme. The human COMT gene has 226 nonsynonymous single nucleotide polymorphisms (nsSNPs) according to public databases. Uncovering of the molecular impacts of nsSNPs on COMT enzyme function and structure may provide standpoint on how COMT nsSNPs affect enzyme activity and contribute to disease development. Therefore, we aimed in this study to predict possible structural and functional damaging effects of all knowns nsSNPs in COMT gene by applying various bioinformatics tools. Two hundred and twenty-six nsSNPs were obtained from Ensembl, HGMD, ClinVar, and dbSNP databases. Twenty-eight nsSNPs were found to be high-risk changes for protein structure. Some of them were detected in extremely conserved sequences have functional and structural properties. Besides, high-risk nsSNPs were also uncovered to change properties of native COMT protein. Our findings demonstrated the significance of COMT high-risk nsSNPs on protein structure and function. We expect that our results will be helpful in future studies concerning experimental evaluation of the COMT gene polymorphisms and/or the association between COMT polymorphisms and disease development.

Consensus interpretation of the p.Met34Thr and p.Val37Ile variants in GJB2 by the ClinGen Hearing Loss Expert Panel

PURPOSE

Pathogenic variants in GJB2 are the most common cause of autosomal recessive sensorineural hearing loss. The classification of c.101T>C/p.Met34Thr and c.109G>A/p.Val37Ile in GJB2 are controversial. Therefore, an expert consensus is required for the interpretation of these two variants.

METHODS

The ClinGen Hearing Loss Expert Panel collected published data and shared unpublished information from contributing laboratories and clinics regarding the two variants. Functional, computational, allelic, and segregation data were also obtained. Case-control statistical analyses were performed.

RESULTS

The panel reviewed the synthesized information, and classified the p.Met34Thr and p.Val37Ile variants utilizing professional variant interpretation guidelines and professional judgment. We found that p.Met34Thr and p.Val37Ile are significantly overrepresented in hearing loss patients, compared with population controls. Individuals homozygous or compound heterozygous for p.Met34Thr or p.Val37Ile typically manifest mild to moderate hearing loss. Several other types of evidence also support pathogenic roles for these two variants.

CONCLUSION

Resolving controversies in variant classification requires coordinated effort among a panel of international multi-institutional experts to share data, standardize classification guidelines, review evidence, and reach a consensus. We concluded that p.Met34Thr and p.Val37Ile variants in GJB2 are pathogenic for autosomal recessive nonsyndromic hearing loss with variable expressivity and incomplete penetrance.

Comparison of Predictive In Silico Tools on Missense Variants in GJB2, GJB6, and GJB3 Genes Associated with Autosomal Recessive Deafness 1A (DFNB1A)

In silico predictive software allows assessing the effect of amino acid substitutions on the structure or function of a protein without conducting functional studies. The accuracy of in silico pathogenicity prediction tools has not been previously assessed for variants associated with autosomal recessive deafness 1A (DFNB1A). Here, we identify in silico tools with the most accurate clinical significance predictions for missense variants of the GJB2 (Cx26), GJB6 (Cx30), and GJB3 (Cx31) connexin genes associated with DFNB1A. To evaluate accuracy of selected in silico tools (SIFT, FATHMM, MutationAssessor, PolyPhen-2, CONDEL, MutationTaster, MutPred, Align GVGD, and PROVEAN), we tested nine missense variants with previously confirmed clinical significance in a large cohort of deaf patients and control groups from the Sakha Republic (Eastern Siberia, Russia): Сх26: p.Val27Ile, p.Met34Thr, p.Val37Ile, p.Leu90Pro, p.Glu114Gly, p.Thr123Asn, and p.Val153Ile; Cx30: p.Glu101Lys; Cx31: p.Ala194Thr. We compared the performance of the in silico tools (accuracy, sensitivity, and specificity) by using the missense variants in GJB2 (Cx26), GJB6 (Cx30), and GJB3 (Cx31) genes associated with DFNB1A. The correlation coefficient (r) and coefficient of the area under the Receiver Operating Characteristic (ROC) curve as alternative quality indicators of the tested programs were used. The resulting ROC curves demonstrated that the largest coefficient of the area under the curve was provided by three programs: SIFT (AUC = 0.833, p = 0.046), PROVEAN (AUC = 0.833, p = 0.046), and MutationAssessor (AUC = 0.833, p = 0.002). The most accurate predictions were given by two tested programs: SIFT and PROVEAN (Ac = 89%, Se = 67%, Sp = 100%, r = 0.75, AUC = 0.833). The results of this study may be applicable for analysis of novel missense variants of the GJB2 (Cx26), GJB6 (Cx30), and GJB3 (Cx31) connexin genes.

Comparative functional characterization of novel non-syndromic GJB2 gene variant p.Gly45Arg and lethal syndromic variant p.Gly45Glu

We characterized a novel GJB2 missense variant, c.133G>A, p.Gly45Arg, and compared it with the only other variant at the same amino acid position of the connexin 26 protein (Cx26) reported to date: c.134G>A, p.Gly45Glu. Whereas both variants are associated with hearing loss and are dominantly inherited, p.Gly45Glu has been implicated in the rare fatal keratitis-ichthyosis-deafness (KID) syndrome, which results in cutaneous infections and septicemia with premature demise in the first year of life. In contrast, p.Gly45Arg appears to be non-syndromic. Subcellular localization experiments in transiently co-transfected HeLa cells demonstrated that Cx26-WT (wild-type) and p.Gly45Arg form gap junctions, whereas Cx26-WT with p.Gly45Glu protein does not. The substitution of a nonpolar amino acid glycine in wildtype Cx26 at position 45 with a negatively charged glutamic acid (acidic) has previously been shown to interfere with Ca²⁺ regulation of hemichannel gating and to inhibit the formation of gap junctions, resulting in cell death. The novel variant p.Gly45Arg, however, changes this glycine to a positively charged arginine (basic), resulting in the formation of dysfunctional gap junctions that selectively affect the permeation of negatively charged inositol 1,4,5-trisphosphate (IP₃) and contribute to hearing loss. Cx26 p.Gly45Arg transfected cells, unlike cells transfected with p.Gly45Glu, thrived at physiologic Ca²⁺ concentrations, suggesting that Ca²⁺ regulation of hemichannel gating is unaffected in Cx26 p.Gly45Arg transfected cells. Thus, the two oppositely charged amino acids that replace the highly conserved uncharged glycine in p.Gly45Glu and p.Gly45Arg, respectively, produce strikingly different effects on the structure and function of the Cx26 protein.

In Silico Analysis of FMR1 Gene Missense SNPs

The FMR1 gene, a member of the fragile X-related gene family, is responsible for fragile X syndrome (FXS). Missense single-nucleotide polymorphisms (SNPs) are responsible for many complex diseases. The effect of FMR1 gene missense SNPs is unknown. The aim of this study, using in silico techniques, was to analyze all known missense mutations that can affect the functionality of the FMR1 gene, leading to mental retardation (MR) and FXS. Data on the human FMR1 gene were collected from the Ensembl database (release 81), National Centre for Biological Information dbSNP Short Genetic Variations database, 1000 Genomes Browser, and NHLBI Exome Sequencing Project Exome Variant Server. In silico analysis was then performed. One hundred-twenty different missense SNPs of the FMR1 gene were determined. Of these, 11.66 % of the FMR1 gene missense SNPs were in highly conserved domains, and 83.33 % were in domains with high variety. The results of the in silico prediction analysis showed that 31.66 % of the FMR1 gene SNPs were disease related and that 50 % of SNPs had a pathogenic effect. The results of the structural and functional analysis revealed that although the R138Q mutation did not seem to have a damaging effect on the protein, the G266E and I304N SNPs appeared to disturb the interaction between the domains and affect the function of the protein. This is the first study to analyze all missense SNPs of the FMR1 gene. The results indicate the applicability of a bioinformatics approach to FXS and other FMR1-related diseases. I think that the analysis of FMR1 gene missense SNPs using bioinformatics methods would help diagnosis of FXS and other FMR1-related diseases.

The pathological effects of connexin 26 variants related to hearing loss by in silico and in vitro analysis

Gap junctions (GJs) are intercellular channels associated with cell-cell communication. Connexin 26 (Cx26) encoded by the GJB2 gene forms GJs of the inner ear, and mutations of GJB2 cause congenital hearing loss that can be syndromic or non-syndromic. It is difficult to predict pathogenic effects using only genetic analysis. Using ionic and biochemical coupling tests, we evaluated the pathogenic effects of Cx26 variants using computational analyses to predict structural abnormalities. For seven out of ten variants, we predicted the variation would result in a loss of GJ function, whereas the others would completely fail to form GJs. Functional studies demonstrated that, although all variants were able to function normally as hetero-oligomeric GJ channels, six variants (p.E47K, p.E47Q, p.H100L, p.H100Y, p.R127L, and p.M195L) did not function normally as homo-oligomeric GJ channels. Interestingly, GJs composed of the Cx26 variant p.R127H were able to function normally, even as homo-oligomeric GJ channels. This study demonstrates the particular location and property of an amino acid are more important mainly than the domain where they belong in the formation and function of GJ, and will provide information that is useful for the accurate diagnosis of hearing loss.

Regulation of connexin signaling by the epigenetic machinery

Connexins and their channels are involved in the control of all aspects of the cellular life cycle, ranging from cell growth to cell death, by mediating extracellular, intercellular and intracellular communication. These multifaceted aspects of connexin-related cellular signaling obviously require strict regulation. While connexin channel activity is mainly directed by posttranslational modifications, connexin expression as such is managed by classical cis/trans mechanisms. Over the past few years, it has become clear that connexin production is equally dictated by epigenetic actions. This paper provides an overview of the role of major determinants of the epigenome, including DNA methylation, histone acetylation and microRNA species, in connexin expression.

Structural and functional impact of missense mutations in TPMT: An integrated computational approach

BACKGROUND

Thiopurine S-methyltransferase (TPMT) detoxifies thiopurine drugs which are used for treatment of various diseases including inflammatory bowel disease (IBD), and hematological malignancies. Individual variation in TPMT activity results from mutations in TPMT gene. In this study, the effects of all the known missense mutations in TPMT enzyme were studied at the sequence and structural level

METHODS

A broad set of bioinformatic tools was used to assess all the known missense mutations affecting enzyme activity. The effects of these mutations on protein stability, aggregation propensity, and residue interaction network were analyzed.

RESULTS

Our results indicate that the missense mutations have diverse effects on TPMT structure and function. Stability and aggregation propensities are affected by various mutations. Several mutations also affect residues in ligand binding site.

CONCLUSIONS

In vitro study of missense mutation is laborious and time-consuming. However, computational methods can be used to obtain information about effects of missense mutations on protein structure. In this study, the effects of most of the mutations on enzyme activity could be explained by computational methods. Thus, the present approach can be used for understanding the protein structure-function relationships.

The p.Cys169Tyr variant of connexin 26 is not a polymorphism

Mutations in the GJB2 gene, which encodes the gap junction protein connexin 26 (Cx26), are the primary cause of hereditary prelingual hearing impairment. Here, the p.Cys169Tyr missense mutation of Cx26 (Cx26C169Y), previously classified as a polymorphism, has been identified as causative of severe hearing loss in two Qatari families. We have analyzed the effect of this mutation using a combination of confocal immunofluorescence microscopy and molecular dynamics simulations. At the cellular level, our results show that the mutant protein fails to form junctional channels in HeLa transfectants despite being correctly targeted to the plasma membrane. At the molecular level, this effect can be accounted for by disruption of the disulfide bridge that Cys169 forms with Cys64 in the wild-type structure (Cx26WT). The lack of the disulfide bridge in the Cx26C169Y protein causes a spatial rearrangement of two important residues, Asn176 and Thr177. In the Cx26WT protein, these residues play a crucial role in the intra-molecular interactions that permit the formation of an intercellular channel by the head-to-head docking of two opposing hemichannels resident in the plasma membrane of adjacent cells. Our results elucidate the molecular pathogenesis of hereditary hearing loss due to the connexin mutation and facilitate the understanding of its role in both healthy and affected individuals.