Missense mutations in dystrophin that trigger muscular dystrophy decrease protein stability and lead to cross-beta aggregates

Biophysical characterization of the dystrophin C-terminal domain: Dystrophin interacts differentially with dystrobrevin isoforms

Duchenne muscular dystrophy (DMD) gene encodes dystrophin, a large multidomain protein. Its nonfunctionality leads to dystrophinopathies like DMD and Becker muscular dystrophy, for which no cure is yet available. A few therapies targeted towards specific mutations can extend the lifespan of patients, although with limited efficacy and high costs, emphasizing the need for more general treatments. Dystrophin's complex structure with poorly understood domains and the presence of multiple isoforms with varied expression patterns in different tissues pose challenges in therapeutic development. The C-terminal (CT) domain of dystrophin is less understood in terms of its structure-function, although it has been shown to perform important functional roles by interacting with another cytosolic protein, dystrobrevin. Dystrophin and dystrobrevin stabilize the sarcolemma membrane by forming a multiprotein complex called dystrophin-associated glycoprotein complex that is destabilized in DMD. Dystrobrevin has two major isoforms, alpha and beta, with tissue-specific expression patterns. Here, we characterize the CT domain of dystrophin and its interactions with the two dystrobrevin isoforms. We show that the CT domain is nonglobular and shows reversible urea denaturation as well as thermal denaturation. It interacts with dystrobrevin isoforms differentially, with differences in binding affinity and the mode of interaction. We further show that the amino acid differences in the CT region of dystrobrevin isoforms contribute to these differences. These results have implications for the stability of dystrophin-associated glycoprotein complex in different tissues and explain the differing symptoms associated with DMD patients affecting organs beyond the skeletal muscles.

Elucidating the role of MLL1 nsSNPs: Structural and functional alterations and their contribution to leukemia development

(1) BACKGROUND

The Mixed lineage leukemia 1 (MLL1) gene, located on chromosome 11q23, plays a pivotal role in histone lysine-specific methylation and is consistently associated with various types of leukemia. Non-synonymous Single Nucleotide Polymorphisms (nsSNPs) have been tied to numerous diseases, including cancers, and have become valuable cancer biomarkers. There's a notable gap in studies probing the influence of SNPs on MLL1 protein structure, function, and subsequent modifications.

(2) METHODS

We utilized an array of bioinformatics tools, including PredictSNP, InterPro, ConSurf, I-Mutant2.0, MUpro, Musitedeep, Project HOPE, RegulomeDB, Mutpred2, and both CScape and CScape Somatic, to meticulously analyze the consequences of nsSNPs in the MLL1 gene.

(3) RESULTS

Out of 2,097 nsSNPs analyzed, 62 were determined to be significantly pathogenic by the PredictSNP tool, with ten crucial MLL1 functional domains identified using InterPro. Additionally, 50 of these nsSNPs had high conservation scores, hinting at potential effects on protein structure and function, while 32 were found to undermine MLL1 protein stability. Notably, four nsSNPs were deemed oncogenic, with two identified as cancer drivers. The nsSNP, D2724G, between the MLL1 protein's FY-rich domains, could disrupt proteolytic cleavage, altering gene expression patterns and potentially promoting cancer.

(4) CONCLUSIONS

Our research provides a comprehensive assessment of nsSNPs' impact in the MLL1 protein structure and function and consequently on leukemia development, suggesting potential avenues for personalized treatment, early detection, improved prognosis, and a deeper understanding of hematological malignancy genesis.

Core and accessory genomic traits of Vibrio cholerae O1 drive lineage transmission and disease severity

In Bangladesh, Vibrio cholerae lineages are undergoing genomic evolution, with increased virulence and spreading ability. However, our understanding of the genomic determinants influencing lineage transmission and disease severity remains incomplete. Here, we developed a computational framework using machine-learning, genome scale metabolic modelling (GSSM) and 3D structural analysis, to identify V. cholerae genomic traits linked to lineage transmission and disease severity. We analysed in-patients isolates from six Bangladeshi regions (2015-2021), and uncovered accessory genes and core SNPs unique to the most recent dominant lineage, with virulence, motility and bacteriophage resistance functions. We also found a strong correlation between V. cholerae genomic traits and disease severity, with some traits overlapping those driving lineage transmission. GSMM and 3D structure analysis unveiled a complex interplay between transcription regulation, protein interaction and stability, and metabolic networks, associated to lifestyle adaptation, intestinal colonization, acid tolerance and symptom severity. Our findings support advancing therapeutics and targeted interventions to mitigate cholera spread.

Unraveling the potential effects of non-synonymous single nucleotide polymorphisms (nsSNPs) on the Protein structure and function of the human SLC30A8 gene on type 2 diabetes and colorectal cancer: An In silico approach

Background and aims

The single nucleotide polymorphisms (SNPs) in SLC30A8 gene have been recognized as contributing to type 2 diabetes (T2D) susceptibility and colorectal cancer. This study aims to predict the structural stability, and functional impacts on variations in non-synonymous SNPs (nsSNPs) in the human SLC30A8 gene using various computational techniques.

Materials and methods

Several in silico tools, including SIFT, Predict-SNP, SNPs&GO, MAPP, SNAP2, PhD-SNP, PANTHER, PolyPhen-1,PolyPhen-2, I-Mutant 2.0, and MUpro, have been used in our study.

Results

After data analysis, out of 336 missenses, the eight nsSNPs, namely R138Q, I141N, W136G, I349N, L303R, E140A, W306C, and L308Q, were discovered by ConSurf to be in highly conserved regions, which could affect the stability of their proteins. Project HOPE determines any significant molecular effects on the structure and function of eight mutated proteins and the three-dimensional (3D) structures of these proteins. The two pharmacologically significant compounds, Luzonoid B and Roseoside demonstrate strong binding affinity to the mutant proteins, and they are more efficient in inhibiting them than the typical SLC30A8 protein using Autodock Vina and Chimera. Increased binding affinity to mutant SLC30A8 proteins has been determined not to influence drug resistance. Ultimately, the Kaplan-Meier plotter study revealed that alterations in SLC30A8 gene expression notably affect the survival rates of patients with various cancer types.

Conclusion

Finally, the study found eight highly deleterious missense nsSNPs in the SLC30A8 gene that can be helpful for further proteomic and genomic studies for T2D and colorectal cancer diagnosis. These findings also pave the way for personalized treatments using biomarkers and more effective healthcare strategies.

Structural impact, ligand-protein interactions, and molecular phenotypic effects of TGF-β1 gene variants: In silico analysis with implications for idiopathic pulmonary fibrosis

BACKGROUND

Idiopathic Pulmonary Fibrosis (IPF) is a chronic interstitial lung disease resulting in progressively deteriorating lung function. Transforming growth factor-β1 (TGF-β1) belongs to the TGF superfamily and exerts a profibrotic role in promoting lung fibrosis by facilitating fibroblast infiltration and activity, extracellular matrix deposition, and inhibition of collagen breakdown, thus promoting tissue remodelling and IPF.

MATERIALS AND METHODS

We evaluated the link between pathogenic TGF-β1 SNPs and IPF pathogenesis and the structure-activity functional consequences of those SNPs on the TGF-β1 protein. Several computational algorithms were merged to address the functional consequences of TGF-β1 gene mutations to protein stability, putative post-translational modification sites, ligand-protein interactions, and molecular phenotypic effects. These included FATHMM, POLYPHEN2, PROVEAN, and SIFT tools (identifying deleterious nsSNPs in the TGF-β1 gene), along with Pmut, PhD-SNP, SNAP, MutPred and the related TMHMM, MARCOIL, and DisProt algorithms (predicting structural disorders). INPS-MD was also used to evaluate the mutation-induced TGF-β1 protein's stability and MODPRED for recognition of post-translational TGF-β1 modification.

RESULTS

In total, 14 major pathogenic variants markedly impact the destabilization of the TGF-β1 protein, with most of these high-risk mutations associated with decreased stability of the TGF-β1 protein as per the I-Mutant, MUpro, and INPS-MD tools. R205W, R185W, R180Q, D86Y, and I300T variants were proposed to participate in the post-translational modifications, thus affecting affect protein-ligand interactions. Furthermore, at-risk genetic variants appear to target conserved regions in the alpha helices, random coils, and extracellular loops, resulting in a varied composition of amino acids, charge, hydrophobicity, and spatial architecture.

CONCLUSIONS

This study manuscript comprehensively analyzes gene variants within the TGF-β1 gene, offering novel insights into their structural and functional implications in interacting with target sites. This study is significant for the development of targeted therapeutic strategies and personalized treatment approaches for patients with inflammatory lung diseases such as IPF.

Prediction of the most deleterious non-synonymous SNPs in the human IL1B gene: evidence from bioinformatics analyses

BACKGROUND

Polymorphisms in IL1B play a significant role in depression, multiple inflammatory-associated disorders, and susceptibility to infection. Functional non-synonymous SNPs (nsSNPs) result in changes in the encoded amino acids, potentially leading to structural and functional alterations in the mutant proteins. So far, most genetic studies have concentrated on SNPs located in the IL1B promoter region, without addressing nsSNPs and their association with multifactorial diseases. Therefore, this study aimed to explore the impact of deleterious nsSNPs retrieved from the dbSNP database on the structure and functions of the IL1B protein.

RESULTS

Six web servers (SIFT, PolyPhen-2, PROVEAN, SNPs&GO, PHD-SNP, PANTHER) were used to analyze the impact of 222 missense SNPs on the function and structure of IL1B protein. Five novel nsSNPs (E100K, T240I, S53Y, D128Y, and F228S) were found to be deleterious and had a mutational impact on the structure and function of the IL1B protein. The I-mutant v2.0 and MUPro servers predicted that these mutations decreased the stability of the IL1B protein. Additionally, these five mutations were found to be conserved, underscoring their significance in protein structure and function. Three of them (T240I, D128Y, and F228S) were predicted to be cancer-causing nsSNPs. To analyze the behavior of the mutant structures under physiological conditions, we conducted a 50 ns molecular dynamics simulation using the WebGro online tool. Our findings indicate that the mutant values differ from those of the IL1B wild type in terms of RMSD, RMSF, Rg, SASA, and the number of hydrogen bonds.

CONCLUSIONS

This study provides valuable insights into nsSNPs located in the coding regions of IL1B, which lead to direct deleterious effects on the functional and structural aspects of the IL1B protein. Thus, these nsSNPs could be considered significant candidates in the pathogenesis of disorders caused by IL1B dysfunction, contributing to effective drug discovery and the development of precision medications. Thorough research and wet lab experiments are required to verify our findings. Moreover, bioinformatic tools were found valuable in the prediction of deleterious nsSNPs.

Conformational fingerprinting with Raman spectroscopy reveals protein structure as a translational biomarker of muscle pathology

Neuromuscular disorders are a group of conditions that can result in weakness of skeletal muscles. Examples include fatal diseases such as amyotrophic lateral sclerosis and conditions associated with high morbidity such as myopathies (muscle diseases). Many of these disorders are known to have abnormal protein folding and protein aggregates. Thus, easy to apply methods for the detection of such changes may prove useful diagnostic biomarkers. Raman spectroscopy has shown early promise in the detection of muscle pathology in neuromuscular disorders and is well suited to characterising the conformational profiles relating to protein secondary structure. In this work, we assess if Raman spectroscopy can detect differences in protein structure in muscle in the setting of neuromuscular disease. We utilise in vivo Raman spectroscopy measurements from preclinical models of amyotrophic lateral sclerosis and the myopathy Duchenne muscular dystrophy, together with ex vivo measurements of human muscle samples from individuals with and without myopathy. Using quantitative conformation profiling and matrix factorisation we demonstrate that quantitative 'conformational fingerprinting' can be used to identify changes in protein folding in muscle. Notably, myopathic conditions in both preclinical models and human samples manifested a significant reduction in α-helix structures, with concomitant increases in β-sheet and, to a lesser extent, nonregular configurations. Spectral patterns derived through non-negative matrix factorisation were able to identify myopathy with a high accuracy (79% in mouse, 78% in human tissue). This work demonstrates the potential of conformational fingerprinting as an interpretable biomarker for neuromuscular disorders.

Impact of nsSNPs in human AIM2 and IFI16 gene: a comprehensive in silico analysis

AIM2 and IFI16 are the most studied members of AIM2-like receptors (ALRs) in humans and share a common N-Terminal PYD domain and C-terminal HIN domain. The HIN domain binds to dsDNA in response to the invasion of bacterial and viral DNA, and the PYD domain directs apoptosis-associated speck-like protein via protein-protein interactions. Hence, activation of AIM2 and IFI16 is crucial for protection against pathogenic assaults, and any genetic variation in these inflammasomes can dysregulate the human immune system. In this study, different computational tools were used to identify the most deleterious and disease-causing non-synonymous single nucleotide polymorphisms (nsSNPs) in AIM2 and IFI16 proteins. Molecular dynamic simulation was performed for the top damaging nsSNPs to study single amino acid substitution-induced structural alterations in AIM2 and IFI16. The observed results suggest that the variants G13V, C304R, G266R, and G266D for AIM2, and G13E and C356F are deleterious and affect structural integrity. We hope that the suggested deleterious nsSNPs and structural dynamics of AIM2 and IFI16 variants will guide future research to better understand the function of these variants with large-scale studies and may assist in fresher therapeutics focusing on these polymorphisms.Communicated by Ramaswamy H. Sarma.

A computational analysis to evaluate deleterious SNPs of GSK3β, a multifunctional and regulatory protein, for metabolism, wound healing, and migratory processes

This study focused on GSK-3β, a critical serine/threonine kinase with diverse cellular functions. However, there is limited understanding of the impact of non-synonymous single nucleotide polymorphisms (nsSNPs) on its structure and function. Through an exhaustive in-silico investigation 12 harmful nsSNPs were predicted from a pool of 172 acquired from the NCBI dbSNP database using 12 established tools that detects deleterious SNPs. Consistently, these nsSNPs were discovered in locations with high levels of conservation. Notably, the three harmful nsSNPs F67C, A83T, and T138I were situated in the active/binding site of GSK-3β, which may affect the protein's capacity to bind to substrates and other proteins. Molecular dynamics simulations revealed that the F67C and T138I mutants had stable structures, indicating rigidness, whereas the A83T mutant was unstable. Analysis of secondary structures revealed different modifications in all mutant forms, which may affect the stability, functioning, and interactions of the protein. These mutations appear to alter the structural dynamics of GSK-3β, which may have functional ramifications, such as the formation of novel secondary structures and variations in coil-to-helix transitions. In conclusion, this study illuminates the possible structural and functional ramifications of these GSK-3 nsSNPs, revealing how protein compactness, stiffness, and interactions may affect biological activities.

In silico analysis of missense SNPs in TNFR1a and their possible therapeutic or pathogenic role in immune diseases

Tumor necrosis factor alpha (TNFa) is an inflammatory cytokine that is involved in the pathogenesis of various inflammatory disorders including rheumatoid arthritis. TNF-alpha receptor I (TNFR1a) is one of the receptors TNFa binds with for its activation. Any variation in this receptor might affect the role of TNFa in successive events. Amino acid residue substitutions might happen in TNFR1a through non-synonymous single nucleotide polymorphisms (nsSNPs) which may alter the functioning of TNFa, hence, identifying any such substitutions is of paramount significance. In this study, six nsSNPs at five different evolutionary conserved regions are predicted to be detrimental to the structure and/or function of TNFR1a by using numerous computational tools. Their 3D models are also proposed in this study. Besides, they were found to reduce the stability and affect the molecular mechanisms of this protein. Two contrasting possibilities might happen because of these substitutions. One, they might reduce the production of TNFa which is overexpressed in inflammatory diseases, hence can play therapeutic role in such diseases. Second, they might possibly hinder the apoptosis to occur which can effectuate the uncontrolled division of cells, hence can be pathogenic in diseases like cancer. Further investigations on these nsSNPs using animal models and at cellular level will open doors to understand the underlying mechanisms behind various diseases.

CRISPR-Cas9 correction in the DMD mouse model is accompanied by upregulation of Dp71f protein

Duchenne muscular dystrophy (DMD) is a severe hereditary disease caused by a deficiency in the dystrophin protein. The most frequent types of disease-causing mutations in the DMD gene are frameshift deletions of one or more exons. Precision genome editing systems such as CRISPR-Cas9 have shown potential to restore open reading frames in numerous animal studies. Here, we applied an AAV-CRISPR double-cut strategy to correct a mutation in the DMD mouse model with exon 8-34 deletion, encompassing the N-terminal actin-binding domain. We report successful excision of the 100-kb genomic sequence, which includes exons 6 and 7, and partial improvement in cardiorespiratory function. While corrected mRNA was abundant in muscle tissues, only a low level of truncated dystrophin was produced, possibly because of protein instability. Furthermore, CRISPR-Cas9-mediated genome editing upregulated the Dp71f dystrophin isoform on the sarcolemma. Given the previously reported Dp71-associated muscle pathology, our results question the applicability of genome editing strategies for some DMD patients with N-terminal mutations. The safety and efficacy of CRISPR-Cas9 constructs require rigorous investigation in patient-specific animal models.

Identification genetic variations in some heat shock protein genes of Tali goat breed and study their structural and functional effects on relevant proteins

BACKGROUND

Animals of different regions have adapted to adverse environmental conditions by modifying their phenotypic and genotypic characteristics in the long run.

OBJECTIVES

In this study, the effect of genetic variations of 10 heat shock protein (HSP) genes (HSP70A4, HSP70A9, HSP40C17, HSP40C27, HSP90AA1, HSP90AB1, HSPB7, HSPB11, HSPD1 and HSPE1) on the three-dimensional protein structure and function of proteins in Tali goat (a tropical breed) were studied and were compared with Saanen goat (as a sensitive breed).

METHODS

A pooled DNA of 15 samples from blood was sequenced and mapped to the goat reference sequence. The bioinformatics analysis was used to identify nsSNPs in the Tali breed and was compared with the Saanen goat. Four online bioinformatics tools (Sorting Intolerant from Tolerant, Protein Variation Effect Analyzer, Polymorphism Phenotyping version2 and Single Nucleotide Polymorphism Database and Gene Ontology) showed three deleterious missense nsSNPs and seven natural missense SNPs in these HSPs genes of Tali goat.

RESULTS

Out of 10 reported nsSNPs, 5 nsSNPs in HSP70A4, 1 nsSNP inHSP70A9, 2 nsSNPs in HSP40C17, 1 nsSNP in HSP40C27 and 1 nsSNP in HSPD1 were detected. ConSurf tools showed that the majority of the predicted nsSNPs occur in conserved sites. Moreover, several post-translational modification (PTM) predictors computed the probability of post-translation change of nsSNPs. The putative phosphorylation and glycosylation sites in HSPs proteins were substitutions rs669769139 and rs666336692 of the Tali goat breed.

CONCLUSION

These results on the effect of type of genetic variants on the function of HSP proteins will assist to predict the resistance to hard conditions in goat breeds. Considering that the identified SNPid rs669769139 (S248) which is located on the N-terminal ATPase domain of HSP70A4 is a PTM site with a highly conserved score and a natural substitution on changing the stability and benign protein that can affect the functional and structural characterization of HSPs protein for adaptation to the local climate.

In-Frame Deletion of Dystrophin Exons 8-50 Results in DMD Phenotype

Mutations that prevent the production of proteins in the DMD gene cause Duchenne muscular dystrophy. Most frequently, these are deletions leading to reading-frame shift. The "reading-frame rule" states that deletions that preserve ORF result in a milder Becker muscular dystrophy. By removing several exons, new genome editing tools enable reading-frame restoration in DMD with the production of BMD-like dystrophins. However, not every truncated dystrophin with a significant internal loss functions properly. To determine the effectiveness of potential genome editing, each variant should be carefully studied in vitro or in vivo. In this study, we focused on the deletion of exons 8-50 as a potential reading-frame restoration option. Using the CRISPR-Cas9 tool, we created the novel mouse model DMDdel8-50, which has an in-frame deletion in the DMD gene. We compared DMDdel8-50 mice to C57Bl6/CBA background control mice and previously generated DMDdel8-34 KO mice. We discovered that the shortened protein was expressed and correctly localized on the sarcolemma. The truncated protein, on the other hand, was unable to function like a full-length dystrophin and prevent disease progression. On the basis of protein expression, histological examination, and physical assessment of the mice, we concluded that the deletion of exons 8-50 is an exception to the reading-frame rule.

Identification of SNPs in hMSH3/MSH6 interaction domain affecting the structure and function of MSH2 protein

MutS homolog 2 (MSH2) is a mismatch repair gene that plays a critical role in DNA repair pathways, and its mutations are associated with different cancers. The present study aimed to find out the single nucleotide polymorphisms (SNPs) of MSH2 protein associated with causing structural and functional changes leading to the development of cancer with the help of computational tools. Four different tools for predicting deleterious SNPs (SIFT, PROVEAN, PANTHER, and PolyPhen), two tools each for identifying disease association (PhD-SNP and SNP&GO) and estimating stability (I-mutant and MUPro) were employed. Homology modeling, energy minimization, and root mean square deviation calculation were used to estimate structural variations. Twenty-seven SNPs and five SNPs (double amino acid change) were identified based on a consensus approach that might be associated with the structural and functional change in MSH2 protein. Molecular docking reveals that six SNPs affect the interaction of MSH2 and MSH6. Twelve identified SNPs were reported to be linked with hereditary nonpolyposis, colorectal cancer, and Lynch syndrome. Further, selected SNPs need to be validated in an in vitro system for their precise association with cancer predisposition.

Fibre optic Raman spectroscopy for the evaluation of disease state in Duchenne muscular dystrophy: an assessment using the mdx model and human muscle

Introduction/Aims

Raman spectroscopy is an emerging technique for the evaluation of muscle disease. In this study we evaluate the ability of in vivo intramuscular Raman spectroscopy to detect the effects of voluntary running in the mdx model of Duchenne muscular dystrophy (DMD). We also compare mdx data with muscle spectra from human DMD patients.

Methods

Thirty 90‐day‐old mdx mice were randomly allocated to an exercised group (48‐hour access to a running wheel) and an unexercised group (n = 15 per group). In vivo Raman spectra were collected from both gastrocnemius muscles and histopathological assessment subsequently performed. Raman data were analyzed using principal component analysis–fed linear discriminant analysis (PCA‐LDA). Exercised and unexercised mdx muscle spectra were compared with human DMD samples using cosine similarity.

Results

Exercised mice ran an average of 6.5 km over 48 hours, which induced a significant increase in muscle necrosis (P = .03). PCA‐LDA scores were significantly different between the exercised and unexercised groups (P < .0001) and correlated significantly with distance run (P = .01). Raman spectra from exercised mice more closely resembled human spectra than those from unexercised mice.

Discussion

Raman spectroscopy provides a readout of the biochemical alterations in muscle in both the mdx mouse and human DMD muscle.

In silico analysis of a novel pathogenic variant c.7G > A in C14orf39 gene identified by WES in a Pakistani family with azoospermia

Infertility is a multifactorial disorder that affects approximately 12% of couples of childbearing ages worldwide. Few studies have been conducted to understand the genetic causes of infertility in depth. The synaptonemal complex (SC), which is essential for the progression of meiosis, is a conserved tripartite structure that binds homologous chromosomes together and is thus required for fertility. This study investigated genetic causes of infertility in a Pakistani consanguineous family containing two patients suffering from non-obstructive azoospermia (NOA). We performed whole-exome sequencing, followed by Sanger sequencing, and identified a novel pathogenic variant (c.7G > A [p.D3N]) in the SC coding gene C14orf39, which was recessively co-segregated with NOA. In silico analysis revealed that charges on wild-type residues were lost, which may result in loss of interactions with other molecules and residues, and a reduction in protein stability occurred, which was caused by the p.D3N mutation. The novel variant generated the mutant protein C14ORF39^D3N, and homozygous mutations in C14orf39 resulted in NOA. The transcriptome profile of C14ORF39 shows that it is specifically expressed in early brain development, which suggests that research in this area is required to study other functions of C14ORF39 in addition to its role in the germline. This research highlights the conserved role of C14orf39/SIX6OS1 in assembly of the SC and its indispensable role in facilitating genetic diagnosis in patients with infertility, which may enable the development of future treatments.

In-silico analysis unravels the structural and functional consequences of non-synonymous SNPs in the human IL-10 gene

Background

Interleukin-10 (IL-10) is an anti-inflammatory cytokine that affects different immune cells. It is also associated with the stimulation of the T and B cells for the production of antibodies. Several genetic polymorphisms in the IL-10 gene have been reported to cause or aggravate certain diseases like inflammatory bowel disease, rheumatoid arthritis, systemic sclerosis, asthma, etc. However, the disease susceptibility and abnormal function of the mutated IL-10 variants remain obscure.

Results

In this study, we used seven bioinformatics tools (SIFT, PROVEAN, PMut, PANTHER, PolyPhen-2, PHD-SNP, and SNPs&GO) to predict the disease susceptible non-synonymous SNPs (nsSNPs) of IL-10. Nine nsSNPs of IL-10 were predicted to be potentially deleterious: R42G, R45Q, F48L, E72G, M95T, A98D, R125S, Y155C, and I168T. Except two, all of the putative deleterious mutations are found in the highly conserved region of IL-10 protein structure, thus affecting the protein's stability. The 3-D structure of mutant proteins was modeled by project HOPE, and the protein–protein interactions were assessed with STRING. The predicted nsSNPs: R42Q, R45Q, F48L, E72G, and I168T are situated in the binding site region of the IL-10R1 receptor. Disruption of binding affinity with its receptor leads to deregulation of the JAK-STAT pathway and results in enhanced inflammation that imbalance in cellular signaling. Finally, Kaplan–Meier Plotter analysis displayed that deregulation of IL-10 expression affects gastric and ovarian cancer patients' survival rate. Thus, IL-10 could be useful as a potential prognostic marker gene for some cancers.

Conclusion

This study has determined the deleterious nsSNPs of IL-10 that might contribute to the malfunction of IL-10 protein and ultimately lead to the IL-10 associated diseases.

In-silico analysis reveals druggable single nucleotide polymorphisms in angiotensin 1 converting enzyme involved in the onset of blood pressure

OBJECTIVE

The Angiotensin 1 converting enzyme (ACE1) gene plays a critical role in regulating blood pressure and thus, it has become a major therapeutic target of antihypertensives. Single nucleotide polymorphisms (SNPs) occurring within a gene most especially at the functional segment of the genes alter the structure-function relationship of that gene.

RESULTS

Our study revealed that five nsSNPs of the ACE1 gene were found to be potentially deleterious and damaging and they include rs2229839, rs14507892, rs12709442, and rs4977 at point mutations P351R, R953Q, I1018T, F1051V, and T1187M. The protein stability predictive tools revealed that all the nsSNPs decreased stability of the protein and the Consurf server which estimates the evolutionary conservation profile of a protein showed that three mutants were in the highly conserved region. In conclusion, this study predicted potential druggable deleterious mutants that can be further explored to understand the pathological basis of cardiovascular disease.

Pathogenic genetic variants from highly connected cancer susceptibility genes confer the loss of structural stability

Genetic polymorphisms in DNA damage repair and tumor suppressor genes have been associated with increasing the risk of several types of cancer. Analyses of putative functional single nucleotide polymorphisms (SNP) in such genes can greatly improve human health by guiding choice of therapeutics. In this study, we selected nine genes responsible for various cancer types for gene enrichment analysis and found that BRCA1, ATM, and TP53 were more enriched in connectivity. Therefore, we used different computational algorithms to classify the nonsynonymous SNPs which are deleterious to the structure and/or function of these three proteins. The present study showed that the major pathogenic variants (V1687G and V1736G of BRCA1, I2865T and V2906A of ATM, V216G and L194H of TP53) might have a greater impact on the destabilization of the proteins. To stabilize the high-risk SNPs, we performed mutation site-specific molecular docking analysis and validated using molecular dynamics (MD) simulation and molecular mechanics/Poisson Boltzmann surface area (MM/PBSA) studies. Additionally, SNPs of untranslated regions of these genes affecting miRNA binding were characterized. Hence, this study will assist in developing precision medicines for cancer types related to these polymorphisms.

Influence of genetic variants in asparaginase pathway on the susceptibility to asparaginase-related toxicity and patients' outcome in childhood acute lymphoblastic leukemia

BACKGROUND

Asparaginase (ASNase) is a key component in the treatment protocols of childhood acute lymphoblastic leukemia (ALL). Asparagine synthetase (ASNS) and the basic region leucine zipper activating transcription factor 5 (ATF5) mediate the anti-leukemic effect of ASNase. Only a few reports studied the association between polymorphisms in these genes and treatment-related toxicity and response. Therefore, the current study aimed to investigate the association of ASNS and ATF5 polymorphisms with the susceptibility to ASNase-related toxicity and disease outcome in a population of childhood ALL Egyptian patients.

METHODS

In this study, 88 children with ALL were enrolled and genotyped for ASNS T629A and ATF5 C362T polymorphisms using allelic discrimination assay.

RESULTS

The studied polymorphisms did not associate with hypersensitivity or thrombosis, while the ATF5 C362T polymorphism was associated significantly with decreased ASNase-associated pancreatitis (AAP) risk under the dominant model. Patients carrying TT/CT genotypes of ATF5 C362T polymorphism had a significantly better overall survival (OS) and longer event-free survival (EFS) compared to patients with CC genotype. Multivariate analysis confirmed the independent prognostic value of the ATF5 C362T dominant model.

CONCLUSION

ATF5 362TT and CT genotypes were associated with decreased risk to develop AAP and better disease outcome demonstrating a low risk for events and superior survival.

The investigation of nonsynonymous SNPs of human SLC6A4 gene associated with depression: An in silico approach

Genetic polymorphism of the SLC6A4 gene is associated with several behavioral disorders, including depression. Since studying the total nonsynonymous single nucleotide polymorphisms (nsSNPs) of the SLC6A4 gene at the population level is a difficult task, we aim to utilize in silico approach to detect the most deleterious nsSNPs of the SLC6A4 gene. In our study, 7 computational tools were used in the initial stage, including SIFT, Polyphen-2, PROVEAN, SNAP2, PhD-SNP, PANTHER, and SNPs&GO to find out the most damaging nsSNPs. In the second phase, we performed structural, functional, and stability analysis of SLC6A4 protein by popular computation tools, including I-Mutant 2.0 and MutPred2. Also, the ConSurf server was utilized to find the conserved region of the SLC6A4 protein to determine the relationship between these conserved regions with high-risk nsSNPs. Based on these analyses, 5 high-risk mutations of the SLC6A4 protein were selected. Then, we carried out comparative modeling by using the Robetta server and aligned the mutant protein model with the native protein structure. Later, we performed the post-translational modification and functional domain analysis of the SLC6A4 protein. This study concludes that Arginine → Tryptophan at position 79 and Arginine → Cysteine at position 104 are the two significant mutations in SLC6A4 protein which might play an essential role in causing diseases. Future studies should take these high-risk nsSNPs (rs1221448303, rs200953188) into consideration while exploring diseases related to the SLC6A4 gene. Besides, our research is the first-ever comprehensive in silico investigation of the SLC6A4 gene. Thus, the findings of this study could be beneficial for developing precision medicines against diseases caused by SLC6A4 malfunction. Furthermore, extensive wet-lab research and experiments on various model organisms might be helpful to investigate the precise role of these damaging nsSNPs of the SLC6A4 gene.

In silico screening and exploration into phenotypic alterations of deleterious oncogenic single nucleotide polymorphisms in HSPB1 gene

A small heat shock protein, HSP27, encoded by HSPB1 gene strongly favors survival, proliferation and metastasis of cancer cells and its expression is dependent on post-translational modifications like phosphorylation. This study performed an extensive in silico screening of 20 deleterious non-synonymous SNPs in the coding region of HSPB1 gene, among which four were identified to be cancer associated. The SNP variant I181S introduced a new phosphorylation site in position 181, which might elevate the protein's activation potential. Emergence of other post-translational modifications was also observed in SNP variants: L144P and E130K.Significant conformational changes were observed in I181S, L144P and E130K SNP variants with respect to wild-type HSP27. These SNPs appear in one among 10⁵ individuals, making them more susceptible towards cancer. This study would therefore, instigate development of novel biomarkers for cancer risk detection and would provide a detailed understanding towards varied cancer susceptibility of human population.

Ixazomib, an oral proteasome inhibitor, exhibits potential effect in dystrophin-deficient mdx mice

Dystrophin deficiency makes the sarcolemma fragile and susceptible to degeneration in Duchenne muscular dystrophy. The proteasome is a multimeric protease complex and is central to the regulation of cellular proteins. Previous studies have shown that proteasome inhibition improved pathological changes in mdx mice. Ixazomib is the first oral proteasome inhibitor used as a therapy in multiple myeloma. This study investigated the effects of ixazomib on the dystrophic muscle of mdx mice. MDX mice were treated with ixazomib (7.5 mg/kg/wk by gavage) or 0.2 mL of saline for 12 weeks. The Kondziela test was performed to measure muscle strength. The tibialis anterior (TA) and diaphragm (DIA) muscles were used for morphological analysis, and blood samples were collected for biochemical measurement. We observed maintenance of the muscle strength in the animals treated with ixazomib. Treatment with ixazomib had no toxic effect on the mdx mouse. The morphological analysis showed a reduction in the inflammatory area and fibres with central nuclei in the TA and DIA muscles and an increase in the number of fibres with a diameter of 20 µm2 in the DIA muscle after treatment with ixazomib. There was an increase in the expression of dystrophin and utrophin in the TA and DIA muscles and a reduction in the expression of osteopontin and TGF-β in the DIA muscle of mdx mice treated with ixazomib. Ixazomib was thus shown to increase the expression of dystrophin and utrophin associated with improved pathological and functional changes in the dystrophic muscles of mdx mice.

CARD8 polymorphism rs2043211 protects against noise-induced hearing loss by causing the dysfunction of CARD8 protein

Inflammation is a complicated process and is considered to be responsible for the development of noise-induced hearing loss (NIHL). CARD8 is an important component of inflammasome that has been implicated in inflammation. To decide the relationship between the polymorphisms of CARD8 gene and NIHL risk and deduce the potential mechanism, three SNPs (rs2043211, rs1062808, and rs12459322) were genotyped in a Chinese population consisting of 610 NIHL cases and 612 normal hearing controls. The possible impacts of SNPs on CARD8 structure and function were assessed using a variety of bioinformatics tools. Plasmids expressing wild-type and/or mutated CARD8 were transfected into HEK293 cells to verify the effect of SNPs on CARD8 protein expression level by western blot. The results revealed that rs2043211 AA genotype and A allele were associated with decreased risk of NIHL. Stratified analysis found that the male, drinking and exposed to noise ≥ 92 dB, subjects harboring rs2043211 A allele had a low risk of NIHL. The haplotype AGG (rs2043211-rs1062808-rs12459322) was significantly associated with a decreased risk of NIHL. SNP rs2043211 was predicted to be deleterious and affects CARD8 protein structure and stability. Furthermore, the functional experiment showed the mutant CARD8 could significantly decrease the CARD8 protein expression level. This study confirms that rs2043211 A allele may reduce NIHL risk by causing the loss of PPI combined with the decreased CARD8 expression level leading to CARD8 functional changes, and it may be one valuable genetic biomarker of NIHL susceptibility for Chinese noise-exposed workers.

Comprehensive in-silico analysis of damage associated SNPs in hOCT1 affecting Imatinib response in chronic myeloid leukemia

Non-synonymous single nucleotide polymorphisms (nsSNPs) in hOCT1 (encoded by SLC22A1 gene) are expected to affect Imatinib uptake in chronic myeloid leukemia (CML). In this study, sequence homology-based genetic analysis of a set of 270 coding SNPs identified 18 nsSNPs to be putatively damaging/deleterious using eight different algorithms. Subsequently, based on conservation of amino acid residues, stability analysis, posttranscriptional modifications, and solvent accessibility analysis, the possible structural-functional relationship was established for high-confidence nsSNPs. Furthermore, based on the modeling results, some dissimilarities of mutant type amino acids from wild-type amino acids such as size, charge, interaction and hydrophobicity were revealed. Three highly deleterious mutations consisting of P283L, G401S and R402G in SLC22A1 gene that may alter the protein structure, function and stability were identified. These results provide a filtered data to explore the effect of uncharacterized nsSNP and find their association with Imatinib resistance in CML.

Effects of Mini-Dystrophin on Dystrophin-Deficient, Human Skeletal Muscle-Derived Cells

BACKGROUND

We are developing a novel therapy for Duchenne muscular dystrophy (DMD), involving the transplantation of autologous, skeletal muscle-derived stem cells that have been genetically corrected to express dystrophin. Dystrophin is normally expressed in activated satellite cells and in differentiated muscle fibres. However, in past preclinical validation studies, dystrophin transgenes have generally been driven by constitutive promoters that would be active at every stage of the myogenic differentiation process, including in proliferating muscle stem cells. It is not known whether artificial dystrophin expression would affect the properties of these cells.

AIMS

Our aims are to determine if mini-dystrophin expression affects the proliferation or myogenic differentiation of DMD skeletal muscle-derived cells.

METHODS

Skeletal muscle-derived cells from a DMD patient were transduced with lentivirus coding for mini-dystrophins (R3-R13 spectrin-like repeats (ΔR3R13) or hinge2 to spectrin-like repeats R23 (ΔH2R23)) with EGFP (enhanced green fluorescence protein) fused to the C-terminus, driven by a constitutive promoter, spleen focus-forming virus (SFFV). Transduced cells were purified on the basis of GFP expression. Their proliferation and myogenic differentiation were quantified by ethynyl deoxyuridine (EdU) incorporation and fusion index. Furthermore, dystrophin small interfering ribonucleic acids (siRNAs) were transfected to the cells to reverse the effects of the mini-dystrophin. Finally, a phospho-mitogen-activated protein kinase (MAPK) array assay was performed to investigate signalling pathway changes caused by dystrophin expression.

RESULTS

Cell proliferation was not affected in cells transduced with ΔR3R13, but was significantly increased in cells transduced with ΔH2R23. The fusion index of myotubes derived from both ΔR3R13- and ΔH2R23 -expressing cells was significantly compromised in comparison to myotubes derived from non-transduced cells. Dystrophin siRNA transfection restored the differentiation of ΔH2R23-expressing cells. The Erk1/2- signalling pathway is altered in cells transduced with mini-dystrophin constructs.

CONCLUSIONS

Ectopic expression of dystrophin in cultured human skeletal muscle-derived cells may affect their proliferation and differentiation capacity. Caution should be taken when considering genetic correction of autologous stem cells to express dystrophin driven by a constitutive promoter.

Pathological evaluation of rats carrying in-frame mutations in the dystrophin gene: a new model of Becker muscular dystrophy

Dystrophin, encoded by the DMD gene on the X chromosome, stabilizes the sarcolemma by linking the actin cytoskeleton with the dystrophin-glycoprotein complex (DGC). In-frame mutations in DMD cause a milder form of X-linked muscular dystrophy, called Becker muscular dystrophy (BMD), characterized by the reduced expression of truncated dystrophin. So far, no animal model with in-frame mutations in Dmd has been established. As a result, the effect of in-frame mutations on the dystrophin expression profile and disease progression of BMD remains unclear. In this study, we established a novel rat model carrying in-frame Dmd gene mutations (IF rats) and evaluated the pathology. We found that IF rats exhibited reduced expression of truncated dystrophin in a proteasome-independent manner. This abnormal dystrophin expression caused dystrophic changes in muscle tissues but did not lead to functional deficiency. We also found that the expression of additional dystrophin named dpX, which forms the DGC in the sarcolemma, was associated with the appearance of truncated dystrophin. In conclusion, the outcomes of this study contribute to the further understanding of BMD pathology and help elucidate the efficiency of dystrophin recovery treatments in Duchenne muscular dystrophy, a more severe form of X-linked muscular dystrophy.

Morpholino Oligomer-Induced Dystrophin Isoforms to Map the Functional Domains in the Dystrophin Protein

Dystrophin plays a crucial role in maintaining sarcolemma stability during muscle contractions, and mutations that prevent the expression of a functional protein cause Duchenne muscular dystrophy (DMD). Antisense oligonucleotide-mediated manipulation of pre-messenger RNA splicing to bypass Duchenne-causing mutations and restore functional dystrophin expression has entered the clinic for the most common DMD mutations. The rationale of "exon skipping" is based upon genotype-phenotype correlations observed in Becker muscular dystrophy, a milder allelic disorder generally characterized by in-frame deletions and internally truncated but semi-functional dystrophin isoforms. However, there is a lack of genotype-phenotype correlations downstream of DMD exon 55, as deletions in this region are rare and most single exon deletions would disrupt the reading frame. Consequently, the amenability of mutations in this region of the DMD gene to exon skipping strategies remains unknown. Here, we induced "Becker muscular dystrophy-like" in-frame dystrophin isoforms in vivo by intraperitoneal injection of peptide-conjugated phosphorodiamidate morpholino oligomers targeting selected exons. The dystrophin isoform encoded by the transcript lacking exons 56+57 appears to be more functional than that encoded by the 58+59-deleted transcript, as determined by higher dystrophin expression, stabilized β-dystroglycan, and less severe dystrophic pathology, indicating some potential for the strategy to address Duchenne-causing mutations affecting these exons.

The chaperonin TRiC is blocked by native and glycated prion protein

Eukaryotic double-ring chaperonin TRiC is an ATP-dependent protein-folding machine. Most of its substrates are known to form large ordered structures from multiple polypeptide chains. Since these structures are similar to fibrillar and oligomeric forms of amyloidogenic proteins, we hypothesized that TRiC may play a role in the development of neurodegenerative diseases of amyloid nature including prion diseases. Enzyme-linked immunosorbent assay showed that monomeric, oligomeric and fibrillar forms of prion protein (PrP) bind strongly to chaperonin TRiC, whereas glycation reduces the prion protein affinity for chaperonin. Nevertheless, dynamic light scattering, electron microscopy and thioflavin T fluorescence confirmed that all studied forms of PrP undergo an amyloid transformation after interaction with chaperonin, but different forms of prion protein are capable of having different effects on the functional state of TRiC. For example, prion protein monomers completely block its ability to reactivate the chaperonin's natural substrate - sperm-specific glyceraldehyde-3-phosphate dehydrogenase (GAPDS). At the same time, PrP oligomers and fibrils only partially prevent the reactivation of GAPDS upon the action of TRiC. The monomeric forms of prion protein glycated by methylglyoxal do not inhibit, but only slow down the chaperone-dependent reactivation of GAPDS. Thus, the interaction of amyloidogenic proteins with chaperonins could cause cell malfunction.

Tissue-Specificity of Dystrophin-Actin Interactions: Isoform-Specific Thermodynamic Stability and Actin-Binding Function of Tandem Calponin-Homology Domains

Genetic mutations in Duchenne muscular dystrophy (DMD) gene affecting the expression of dystrophin protein lead to a number of muscle disorders collectively called dystrophinopathies. In addition to muscle dystrophin, mutations in brain-specific dystrophin isoforms, in particular those that are expressed in the brain cortex and Purkinje neurons, result in cognitive impairment associated with DMD. These isoforms carry minor variations in the flanking region of the N-terminal actin-binding domain (ABD1) of dystrophin, which is composed of two calponin-homology (CH) domains in tandem. Determining the effect of these sequence variations is critical for understanding the mechanisms that govern varied symptoms of the disease. We studied the impact of differences in the N-terminal flanking region on the structure and function of dystrophin tandem CH domain isoforms. The amino acid changes did not affect the global structure of the protein but drastically affected the thermodynamic stability, with the muscle isoform more stable than the brain and Purkinje isoforms. Actin binding investigated with actin from different sources (skeletal muscle, smooth muscle, cardiac muscle, and platelets) revealed that the muscle isoform binds to filamentous actin (F-actin) with a lower affinity compared to the brain and Purkinje isoforms, and a similar trend was observed with actin from different sources. In addition, all isoforms showed a higher affinity to smooth muscle actin in comparison to actin from other sources. In conclusion, tandem CH domain isoforms might be using minor sequence variations in the N-terminal flanking regions to modulate their thermodynamic stability and actin-binding function, thus leading to specificity in dystrophin-actin interactions in various tissues.

Molecular Modelling and Dynamics Study of nsSNP in STXBP1 Gene in Early Infantile Epileptic Encephalopathy Disease

Early Infantile Epileptic Encephalopathy (known as Ohtahara Syndrome) is one of the most severe and earliest forms of epilepsy, characterized by early seizures onset. It affects newborns and children between two and six years old. Among the genes that have been associated with early infantile epileptic encephalopathy, the STXBP1 gene, which encodes the Syntaxin binding protein1a that is involved in SNARE complex formation, contributes to synaptic vesicles exocytosis. The aim of this study was to identify the most pathogenic polymorphisms of STXBP1 gene and determine their impact on the structure and stability of Stxbp1 protein. The high-risk nonsynonymous single nucleotide polymorphisms (nsSNPs) in the STXBP1 gene were predicted using 13 bioinformatics tools. The conservation analysis was realized by CONSURF web server. The analysis of the impact of the pathogenic SNPs on the structure of Stxbp1 protein was realized using YASARA software, and the molecular dynamics simulation was performed using GROMACS software. Out of 245 nsSNPs, we identified 11 (S42P, H103D R190W, R235G, D238E, L256P, P335S, C354Y, L365V, R406C, and G544D) as deleterious using in silico prediction tools. Conservation analysis results revealed that all these nsSNPs were located in conserved regions. The comparison of the hydrogen and hydrophobic interactions in the wild type Stxbp1 structure and its mutant forms showed that all these nsSNPs affect the protein structure on different levels. The molecular dynamics simulations revealed that the total of nsSNPs affect the protein stability, residual fluctuation, and the compaction at different levels. This study provides helpful information on high risk nsSNPs that may affect the Stxbp1 protein structure and function. Thus, these variants should be taken into consideration during the genetic screening of patients suffering from early infantile epileptic encephalopathy.

Large in-frame 5' deletions in DMD associated with mild Duchenne muscular dystrophy: Two case reports and a review of the literature

Duchenne muscular dystrophy is caused by mutations in the dystrophin-encoding DMD gene. While Duchenne is most commonly caused by large intragenic deletions that cause frameshift and complete loss of dystrophin expression, in-frame deletions in DMD can result in the expression of internally truncated dystrophin proteins and may be associated with a milder phenotype. In this study, we describe two individuals with large in-frame 5' deletions (exon 3-23 and exon 3-28) that remove the majority of the N-terminal region, including part of the actin binding and central rod domains. Both patients had progressive muscle weakness during childhood but are observed to have a relatively mild disease course compared to typical Duchenne. We show that in muscle biopsies from both patients, truncated dystrophin is expressed at the sarcolemma. We have additionally developed a patient-specific fibroblast-derived cell model, which can be inducibly reprogrammed to form myotubes that largely recapitulate biopsy findings for the patient with the exon 3-23 deletion, providing a culture model for future investigation of this unusual case. We discuss these mutations in the context of previously reported 5' in-frame DMD deletions and relevant animal models, and review the spectrum of phenotypes associated with these deletions.

Importance of muscle biopsy to establish pathogenicity of DMD missense and splice variants

A precise genetic diagnosis of a dystrophinopathy has far-reaching implications for affected boys and their families. We present three boys with DMD single nucleotide variants associated with Becker muscular dystrophy presenting with myalgia, reduced exercise capacity, neurodevelopmental symptoms and elevated creatine kinase. The DMD variants were difficult to classify: AIII:1 a synonymous variant in exon 13 c.1602G>A, p.Lys534Lys; BIII:1 an essential splice-site variant in intron 33 c.4674+1G>A, and CII:1 a missense mutation within the cysteine-rich domain, exon 66 c.9619T>C, p.Cys3207Arg. Complementary DNA (cDNA) analysis using muscle-derived mRNA established splice-altering effects of variants for AIII:1 and BIII:1, and normal splicing in CII:1. Western blot analysis demonstrated mildly to moderately reduced dystrophin levels (17.6 - 36.1% the levels of controls), supporting dystrophinopathy as a probable diagnosis. These three cases highlight the diagnostic utility of muscle biopsy for mRNA studies and western blot to investigate DMD variants of uncertain pathogenicity, by exploring effects on splicing and dystrophin protein levels.

Identification of deleterious SNPs and their effects on BCL11A, the master regulator of fetal hemoglobin expression

The B-cell lymphoma/leukemia 11A protein (encoded by BCL11A gene) is a key regulator of fetal-to-adult hemoglobin switching as seen in post-natal life. Although genetic polymorphisms like SNPs in BCL11A gene are expected to affect fetal hemoglobin (HbF) expression levels, yet their implications are poorly studied. This study utilizes a computational approach to identify the deleterious SNPs which may affect the structure and function of BCL11A protein. The study also generated a 3D structure of native and mutants. The analysis identified two SNPs in BCL11A as highly deleterious: N391K and C414S which are expected to affect structure and stability of the protein. According to conservation analysis, both residues N391 and C414 were identified as highly conserved. Additionally, post-translational modification sites were predicted at both sites. Ligand binding sites were also predicted in N391 and C414. Therefore, N391K and C414S in BCL11A can considered as important candidates to mediate HbF variation.

Identification and structural characterization of deleterious non-synonymous single nucleotide polymorphisms in the human SKP2 gene

In SCF (Skp, Cullin, F-box) ubiquitin-protein ligase complexes, S-phase kinase 2 (SKP2) is one of the major players of F-box family, that is responsible for the degradation of several important cell regulators and tumor suppressor proteins. Despite of having significant evidence for the role of SKP2 on tumorgenesis, there is a lack of available data regarding the effect of non-synonymous polymorphisms. In this communication, the structural and functional consequences of non-synonymous single nucleotide polymorphisms (nsSNPs) of SKP2 have been reported by employing various computational approaches and molecular dynamics simulation. Initially, several computational tools like SIFT, PolyPhen-2, PredictSNP, I-Mutant 2.0 and ConSurf have been implicated in this study to explore the damaging SNPs. In total of 172 nsSNPs, 5 nsSNPs were identified as deleterious and 3 of them were predicted to be decreased the stability of protein. Guided from ConSurf analysis, P101L (rs761253702) and Y346C (rs755010517) were categorized as the highly conserved and functional disrupting mutations. Therefore, these mutations were subjected to three dimensional model building and molecular dynamics simulation study for the detailed structural consequences upon the mutations. The study revealed that P101L and Y346C mutations increased the flexibility and changed the structural dynamics. As both these mutations are located in the most functional regions of SKP2 protein, these computational insights might be helpful to consider these nsSNPs for wet-lab confirmatory analysis as well as in rationalizing future population based studies and structure based drug design against SKP2.

Mouse models of two missense mutations in actin-binding domain 1 of dystrophin associated with Duchenne or Becker muscular dystrophy

Missense mutations in the dystrophin protein can cause Duchenne muscular dystrophy (DMD) or Becker muscular dystrophy (BMD) through an undefined pathomechanism. In vitro studies suggest that missense mutations in the N-terminal actin-binding domain (ABD1) cause protein instability, and cultured myoblast studies reveal decreased expression levels that can be restored to wild-type with proteasome inhibitors. To further elucidate the pathophysiology of missense dystrophin in vivo, we generated two transgenic mdx mouse lines expressing L54R or L172H mutant dystrophin, which correspond to missense mutations identified in human patients with DMD or BMD, respectively. Our biochemical, histologic and physiologic analysis of the L54R and L172H mice show decreased levels of dystrophin which are proportional to the phenotypic severity. Proteasome inhibitors were ineffective in both the L54R and L172H mice, yet mice homozygous for the L172H transgene were able to express even higher levels of dystrophin which caused further improvements in muscle histology and physiology. Given that missense dystrophin is likely being degraded by the proteasome but whole body proteasome inhibition was not possible, we screened for ubiquitin-conjugating enzymes involved in targeting dystrophin to the proteasome. A myoblast cell line expressing L54R mutant dystrophin was screened with an siRNA library targeting E1, E2 and E3 ligases which identified Amn1, FBXO33, Zfand5 and Trim75. Our study establishes new mouse models of dystrophinopathy and identifies candidate E3 ligases that may specifically regulate dystrophin protein turnover in vivo.

Dynamics of Dystrophin's Actin-Binding Domain

We have used pulsed electron paramagnetic resonance, calorimetry, and molecular dynamics simulations to examine the structural mechanism of binding for dystrophin's N-terminal actin-binding domain (ABD1) and compare it to utrophin's ABD1. Like other members of the spectrin superfamily, dystrophin's ABD1 consists of two calponin-homology (CH) domains, CH1 and CH2. Several mutations within dystrophin's ABD1 are associated with the development of severe degenerative muscle disorders Duchenne and Becker muscular dystrophies, highlighting the importance of understanding its structural biology. To investigate structural changes within dystrophin ABD1 upon binding to actin, we labeled the protein with spin probes and measured changes in inter-CH domain distance using double-electron electron resonance. Previous studies on the homologous protein utrophin showed that actin binding induces a complete structural opening of the CH domains, resulting in a highly ordered ABD1-actin complex. In this study, double-electron electron resonance shows that dystrophin ABD1 also undergoes a conformational opening upon binding F-actin, but this change is less complete and significantly more structurally disordered than observed for utrophin. Using molecular dynamics simulations, we identified a hinge in the linker region between the two CH domains that grants conformational flexibility to ABD1. The conformational dynamics of both dystrophin's and utrophin's ABD1 showed that compact conformations driven by hydrophobic interactions are preferred and that extended conformations are energetically accessible through a flat free-energy surface. Considering that the binding free energy of ABD1 to actin is on the order of 6-7 kcal/mole, our data are compatible with a mechanism in which binding to actin is largely dictated by specific interactions with CH1, but fine tuning of the binding affinity is achieved by the overlap between conformational ensembles of ABD1 free and bound to actin.

Identification and in silico analysis of functional SNPs of human TAGAP protein: A comprehensive study

Genetic polymorphisms in TAGAP gene have been associated with many diseases including rheumatoid arthritis, multiple sclerosis and other autoimmune disorders. Identifying functional SNPs in such disease associated genes is an uphill task hence before planning larger population study, it is better to scrutinize putative functional SNPs. In this study we used various computational approaches to identify nsSNPs which are deleterious to the structure and/or function of TAGAP protein that might be causing these diseases. Computational analysis was performed by five different in silico tools including SIFT, PROVEAN, PolyPhen-2, PhD-SNP and SNPs&GO. The study concludes that mutations of Glycine → Glutamic Acid at position 120, Glycine → Tryptophan at position 141 and Valine → Methionine at position 151 are major mutations in native TAGAP protein which might contribute to its malfunction and ultimately causing disease. The study also proposed 3D structures of native TAGAP protein and its three mutants. Future studies should consider these nsSNPs as main target mutations in various diseases involving TAGAP malfunction. This is the first comprehensive study, where TAGAP gene variants were analyzed using in silico tools hence will be of great help while considering large scale studies and also in developing precision medicines for cure of diseases related to these polymorphisms. Furthermore, animal models of various autoimmune diseases and having these mutations might be of help in exploring their precise roles.

A dystroglycan mutation (p.Cys667Phe) associated to muscle-eye-brain disease with multicystic leucodystrophy results in ER-retention of the mutant protein

Dystroglycan (DG) is a cell adhesion complex composed by two subunits, the highly glycosylated α-DG and the transmembrane β-DG. In skeletal muscle, DG is involved in dystroglycanopathies, a group of heterogeneous muscular dystrophies characterized by a reduced glycosylation of α-DG. The genes mutated in secondary dystroglycanopathies are involved in the synthesis of O-mannosyl glycans and in the O-mannosylation pathway of α-DG. Mutations in the DG gene (DAG1), causing primary dystroglycanopathies, destabilize the α-DG core protein influencing its binding to modifying enzymes. Recently, a homozygous mutation (p.Cys699Phe) hitting the β-DG ectodomain has been identified in a patient affected by muscle-eye-brain disease with multicystic leucodystrophy, suggesting that other mechanisms than hypoglycosylation of α-DG could be implicated in dystroglycanopathies. Herein, we have characterized the DG murine mutant counterpart by transfection in cellular systems and high-resolution microscopy. We observed that the mutation alters the DG processing leading to retention of its uncleaved precursor in the endoplasmic reticulum. Accordingly, small-angle X-ray scattering data, corroborated by biochemical and biophysical experiments, revealed that the mutation provokes an alteration in the β-DG ectodomain overall folding, resulting in disulfide-associated oligomerization. Our data provide the first evidence of a novel intracellular mechanism, featuring an anomalous endoplasmic reticulum-retention, underlying dystroglycanopathy.

Functional correction of dystrophin actin binding domain mutations by genome editing

Dystrophin maintains the integrity of striated muscles by linking the actin cytoskeleton with the cell membrane. Duchenne muscular dystrophy (DMD) is caused by mutations in the dystrophin gene (DMD) that result in progressive, debilitating muscle weakness, cardiomyopathy, and a shortened lifespan. Mutations of dystrophin that disrupt the amino-terminal actin-binding domain 1 (ABD-1), encoded by exons 2-8, represent the second-most common cause of DMD. In the present study, we compared three different strategies for CRISPR/Cas9 genome editing to correct mutations in the ABD-1 region of the DMD gene by deleting exons 3-9, 6-9, or 7-11 in human induced pluripotent stem cells (iPSCs) and by assessing the function of iPSC-derived cardiomyocytes. All three exon deletion strategies enabled the expression of truncated dystrophin protein and restoration of cardiomyocyte contractility and calcium transients to varying degrees. We show that deletion of exons 3-9 by genomic editing provides an especially effective means of correcting disease-causing ABD-1 mutations. These findings represent an important step toward eventual correction of common DMD mutations and provide a means of rapidly assessing the expression and function of internally truncated forms of dystrophin-lacking portions of ABD-1.

The N-Terminal Flanking Region Modulates the Actin Binding Affinity of the Utrophin Tandem Calponin-Homology Domain

Despite sharing a high degree of sequence similarity, the tandem calponin-homology (CH) domain of utrophin binds to actin 30 times stronger than that of dystrophin. We have previously shown that this difference in actin binding affinity could not be ascribed to the differences in inter-CH-domain linkers [Bandi, S., et al. (2015) Biochemistry 54, 5480-5488]. Here, we examined the role of the N-terminal flanking region. The utrophin tandem CH domain contains a 27-residue flanking region before its CH1 domain. We examined its effect by comparing the structure and function of full-length utrophin tandem CH domain Utr(1-261) and its truncated Utr(28-261) construct. Both full-length and truncated constructs are monomers in solution, with no significant differences in their secondary or tertiary structures. Truncated construct Utr(28-261) binds to actin 30 times weaker than that of the full-length Utr(1-261), similar to that of the dystrophin tandem CH domain with a much shorter flanking region. Deletion of the N-terminal flanking region stabilizes the CH1 domain. The magnitude of the change in binding free energy upon truncation is similar to that of the change in thermodynamic stability. The isolated N-terminal peptide by itself is significantly random coil and does not bind to F-actin in the affinity range of Utr(1-261) and Utr(28-261). These results indicate that the N-terminal flanking region significantly affects the actin binding affinity of tandem CH domains. This observation further stresses that protein regions other than the three actin-binding surfaces identified earlier, irrespective of whether they directly bind to actin, also contribute to the actin binding affinity of tandem CH domains.

The Structurally Plastic CH2 Domain Is Linked to Distinct Functions of Fimbrins/Plastins

Fimbrins/plastins have been implicated in the generation of distinct actin structures, which are linked to different cellular processes. Historically, fimbrins/plastins were mainly considered as generating tight actin bundles. Here, we demonstrate that different members of the fimbrin/plastin family have diverged biochemically during evolution to generate either tight actin bundles or loose networks with distinct biochemical and biophysical properties. Using the phylogenetically and functionally distinct Arabidopsis fimbrins FIM4 and FIM5 we found that FIM4 generates both actin bundles and cross-linked actin filaments, whereas FIM5 only generates actin bundles. The distinct functions of FIM4 and FIM5 are clearly observed at single-filament resolution. Domain swapping experiments showed that cooperation between the conformationally plastic calponin-homology domain 2 (CH2) and the N-terminal headpiece determines the function of the full-length protein. Our study suggests that the structural plasticity of fimbrins/plastins has biologically meaningful consequences, and provides novel insights into the structure-function relationship of fimbrins/plastins as well as shedding light on how cells generate distinct actin structures.

The Dystrophin Complex: Structure, Function, and Implications for Therapy

The dystrophin complex stabilizes the plasma membrane of striated muscle cells. Loss of function mutations in the genes encoding dystrophin, or the associated proteins, trigger instability of the plasma membrane, and myofiber loss. Mutations in dystrophin have been extensively cataloged, providing remarkable structure-function correlation between predicted protein structure and clinical outcomes. These data have highlighted dystrophin regions necessary for in vivo function and fueled the design of viral vectors and now, exon skipping approaches for use in dystrophin restoration therapies. However, dystrophin restoration is likely more complex, owing to the role of the dystrophin complex as a broad cytoskeletal integrator. This review will focus on dystrophin restoration, with emphasis on the regions of dystrophin essential for interacting with its associated proteins and discuss the structural implications of these approaches.

Pathogenic mutation R959W alters recognition dynamics of dysferlin inner DysF domain

Dysferlin, a 220 kD protein, plays a major role in regulating plasma membrane repair in muscle cells. Mutations in the dysferlin inner DysF domain are known to cause different types of muscular dystrophy, including limb-girdle muscular dystrophy type 2B (LGMD2B) and Miyoshi myopathy (MM). Replacement of arginine in position 959 by tryptophan has been frequently associated with both LGMD2B and MM, but the molecular mechanisms by which this mutation alters dysferlin function remain unknown. In this study, we have used protein binding site predictions and microsecond molecular dynamics (MD) simulations to determine the effects pathogenic mutation R959W on the structural dynamics of dysferlin inner DysF domain. Analysis of 2 μs long MD trajectories revealed that mutation R959W does not induce local destabilization, unfolding or misfolding of the domain. We used a binding site predictor to discover a protein-binding site (residues T958-I966 and E1031-H1037) that resembles pincers in shape. Cartesian principal component analysis and interresidue distance distributions of the wild-type domain showed that the predicted protein-binding site undergoes a pincer motion, and populates two structural states, open and closed. We found that mutation R959W inhibits the pincer motion of the protein-binding site and completely shifts the equilibrium toward the open state. These differences in the structural dynamics of the predicted binding site suggest that mutation R959W alters recognition dynamics of the inner DysF domain. Based on these findings and on previous experimental studies, we propose a novel role for the inner DysF domain in muscle membrane repair through recruitment of dysferlin to plasma membrane. In conclusion, these findings have important implications for our understanding of the structural aspects of muscular dystrophies in atomic-level resolution.