Structural genomics approach to investigate deleterious impact of nsSNPs in conserved telomere maintenance component 1

Probing the effects of single point mutations in the GKWWRPS motif on the PNAIG motif within Loop 2 of sclerostin (SOST) using in-silico techniques

Sclerostin (SOST), a Wnt signaling pathway inhibitor, is involved in the pathogenesis of skeletal disorders. This study investigated the impact of the GKWWRPS motif on the PNAIG motif in Loop 2 of SOST, which is accountable for the interactions with the LRP6 protein that triggers the down-regulation of the Wnt signaling pathway. Single amino acid mutations on the GKWWRPS motif, hypothesized to have a probable stabilization effect towards the PNAIG motif, led to a significant reduction in the primary interactions between the SOST and LRP6 proteins. Protein-protein docking and molecular dynamic studies were conducted to investigate the role of the motif. The study found that a solitary mutation in the GKWWRPS motif significantly reduced the primary interactions between SOST and LRP6 proteins, except for probable cold-spot residues. The study's findings establish the GKWWRPS motif as a promising target for therapeutic interventions. Based on the obtained results, it can be inferred that alterations implemented within the GKWWRPS motif could lead to the destabilization of the PNAIG motif, which would directly modulate the interactions between the SOST and LRP6 proteins. The present investigation thus presents novel opportunities in the field of anti-sclerostin interventions.

De novo structure prediction of meteorin and meteorin-like protein for identification of domains, functional receptor binding regions, and their high-risk missense variants

Meteorin (Metrn) and Meteorin-like (Metrnl) are homologous secreted proteins involved in neural development and metabolic regulation. In this study, we have performed de novo structure prediction and analysis of both Metrn and Metrnl using Alphafold2 (AF2) and RoseTTAfold (RF). Based on the domain and structural homology analysis of the predicted structures, we have identified that these proteins are composed of two functional domains, a CUB domain and an NTR domain, connected by a hinge/loop region. We have identified the receptor binding regions of Metrn and Metrnl using the machine-learning tools ScanNet and Masif. These were further validated by docking Metrnl with its reported KIT receptor, thus establishing the role of each domain in the receptor interaction. Also, we have studied the effect of non-synonymous SNPs on the structure and function of these proteins using an array of bioinformatics tools and selected 16 missense variants in Metrn and 10 in Metrnl that can affect the protein stability. This is the first study to comprehensively characterize the functional domains of Metrn and Metrnl at their structural level and identify the functional domains, and protein binding regions. This study also highlights the interaction mechanism of the KIT receptor and Metrnl. The predicted deleterious SNPs will allow further understanding of the role of these variants in modulating the plasma levels of these proteins in disease conditions such as diabetes.Communicated by Ramaswamy H. Sarma.

In silico exploration and molecular dynamics of deleterious SNPs on the human TERF1 protein triggering male infertility

By limiting chromosome erosion and end-to-end fusions, telomere integrity is critical for chromosome stability and cell survival. During mitotic cycles or due to environmental stresses, telomeres become progressively shorter and dysfunctional, thus triggering cellular senescence, genomic instability and cell death. To avoid such consequences, the telomerase action, as well as the Shelterin and CST complexes, assure the telomere's protection. Telomeric repeat binding factor 1 (TERF1), which is one of the primary components of the Shelterin complex, binds directly to the telomere and controls its length and function by regulating the telomerase activity. Several reports about TERF1 gene variations have been associated with different diseases, and some of them have linked these variations to male infertility. Hence, this paper can be advantageous to investigate the association between the missense variants of the TERF1 gene and the susceptibility to male infertility. The stepwise prediction of SNPs pathogenicity followed in this study was based on stability and conservation analysis, post-translational modification, secondary structure, functional interaction prediction, binding energy evaluation and finally molecular dynamic simulation. Prediction matching among the tools revealed that out of 18 SNPs, only four (rs1486407144, rs1259659354, rs1257022048 and rs1320180267) were predicted as the most damaging and highly deleterious SNPs affecting the TERF1 protein and its molecular dynamics when interacting with the TERB1 protein by influencing the function, structural stability, flexibility and compaction of the overall complex. Interestingly, these polymorphisms should be considered during genetic screening so they can be used effectively as genetic biomarkers for male infertility diagnosis.Communicated by Ramaswamy H. Sarma.

Delineating the impact of pathogenic mutations on the conformational dynamics of HDL's vital protein ApoA1: a combined computational and molecular dynamic simulation approach

Apolipoprotein A1 (ApoA1), is the important component of high-density lipoproteins (HDL), that has key role in HDL biogenesis, cholesterol trafficking, and reverse cholesterol transport (RCT). Non-synonymous Single Nucleotide Polymorphisms (nsSNPs) in ApoA1 have been linked to cardiovascular diseases and amyloidosis as they alter the protein's native structure and function. Therefore in this study, we attempted to understand the molecular pathogenicity profile of nsSNPs of ApoA1 using various computational approaches. We used state-of-the-art computational methods to thoroughly investigate the 295 ApoA1 nsSNPs at sequence and structural levels. Seven nsSNPs (L13R, L84R, L84P, L99P, R173P, L187P, and L238P) out of 295 were classified as the most deleterious and destabilizing. In order to estimate the effect of such destabilizing mutations on the protein conformation, all-atom molecular dynamics simulations (MDS) of ApoA1 wild-type (WT), L99P and R173P for 100 ns, was carried out using GROMACS 5.0.1 package. The MD simulation investigation revealed significant structural alterations in L99P and R173P. In addition, they had changed principal component analysis and electrostatic surface potential, decreased structural compactness, and intramolecular hydrogen bonds, which supported the rationale underpinning ApoA1 dysfunction with such mutations. This work sheds light on ApoA1 dysfunction due to single amino acid alterations, and offers new insight into the molecular basis of ApoA1-related diseases progression.Communicated by Ramaswamy H. Sarma.

Deleterious single nucleotide polymorphisms (SNPs) of human IFNAR2 gene facilitate COVID-19 severity in patients: a comprehensive in silico approach

In humans, the dimeric receptor complex IFNAR2-IFNAR1 accelerates cellular response triggered by type I interferon (IFN) family proteins in response to viral infection including Coronavirus infection. Studies have revealed the association of the IFNAR2 gene with severe illness in Coronavirus infection and indicated the association of genomic variants, i.e. single nucleotide polymorphisms (SNPs). However, comprehensive analysis of SNPs of the IFNAR2 gene has not been performed in both coding and non-coding region to find the causes of loss of function of IFNAR2 in COVID-19 patients. In this study, we have characterized coding SNPs (nsSNPs) of IFNAR2 gene using different bioinformatics tools and identified deleterious SNPs. We found 9 nsSNPs as pathogenic and disease-causing along with a decrease in protein stability. We employed molecular docking analysis that showed 5 nsSNPs to decrease binding affinity to IFN. Later, MD simulations showed that P136R mutant may destabilize crucial binding with the IFN molecule in response to COVID-19. Thus, P136R is likely to have a high impact on disrupting the structure of the IFNAR2 protein. GTEx portal analysis predicted 14 sQTLs and 5 eQTLs SNPs in lung tissues hampering the post-transcriptional modification (splicing) and altering the expression of the IFNAR2 gene. sQTLs and eQTLs SNPs potentially explain the reduced IFNAR2 production leading to severe diseases. These mutants in the coding and non-coding region of the IFNAR2 gene can help to recognize severe illness due to COVID 19 and consequently assist to develop an effective drug against the infection.Communicated by Ramaswamy H. Sarma.

Whole exome sequencing of high-risk neuroblastoma identifies novel non-synonymous variants

Neuroblastoma (NBL), one of the main death-causing cancers in children, is known for its remarkable genetic heterogeneity and varied patient outcome spanning from spontaneous regression to widespread disease. Specific copy number variations and single gene rearrangements have been proven to be associated with biological behavior and prognosis; however, there is still an unmet need to enlarge the existing armamentarium of prognostic and therapeutic targets. We performed whole exome sequencing (WES) of samples from 18 primary tumors and six relapse samples originating from 18 NBL patients. Our cohort consists of 16 high-risk, one intermediate, and one very low risk patient. The obtained results confirmed known mutational hotspots in ALK and revealed other non-synonymous variants of NBL-related genes (TP53, DMD, ROS, LMO3, PRUNE2, ERBB3, and PHOX2B) and of genes cardinal for other cancers (KRAS, PIK3CA, and FLT3). Beyond, GOSeq analysis determined genes involved in biological adhesion, neurological cell-cell adhesion, JNK cascade, and immune response of cell surface signaling pathways. We were able to identify novel coding variants present in more than one patient in nine biologically relevant genes for NBL, including TMEM14B, TTN, FLG, RHBG, SHROOM3, UTRN, HLA-DRB1, OR6C68, and XIRP2. Our results may provide novel information about genes and signaling pathways relevant for the pathogenesis and clinical course in high-risk NBL.

In Silico Analysis Revealed Five Novel High-Risk Single-Nucleotide Polymorphisms (rs200384291, rs201163886, rs193141883, rs201139487, and rs201723157) in ELANE Gene Causing Autosomal Dominant Severe Congenital Neutropenia 1 and Cyclic Hematopoiesis

Single-nucleotide polymorphisms in the ELANE (Elastase, Neutrophil Expressed) gene are associated with severe congenital neutropenia, while the ELANE gene provides instructions for making a protein called neutrophil elastase. We identified disease susceptibility single-nucleotide polymorphisms (SNPs) in the ELANE gene using several computational tools. We used cutting-edge computational techniques to investigate the effects of ELANE mutations on the sequence and structure of the protein. Our study suggested that eight nsSNPs (rs28931611, rs57246956, rs137854448, rs193141883, rs201723157, rs201139487, rs137854451, and rs200384291) are the most deleterious in ELANE gene and disturb protein structure and function. The mutants F218L, R34W, G203S, R193W, and T175M have not yet been identified in patients suffering from SCN and cyclic hematopoiesis, while C71Y, P139R, C151Y, G214R, and G203C reported in our study are already associated with both of the disorders. These mutations are shown to destabilize structure and disrupt ELANE protein activation, splicing, and folding and might diminish trypsin-like serine protease efficiency. Prediction of posttranslation modifications highlighted the significance of deleterious nsSNPs because some of nsSNPs affect potential phosphorylation sites. Gene-gene interactions showed the relation of ELANE with other genes depicting its importance in numerous pathways and coexpressions. We identified the deleterious nsSNPs, constructed mutant protein structures, and evaluated the impact of mutation by employing molecular docking. This research sheds light on how ELANE failure upon mutation results in disease progression, including congenital neutropenia, and validation of these novel predicted nsSNPs is required through the wet lab.

Comparative analysis of web-based programs for single amino acid substitutions in proteins

Single amino-acid substitution in a protein affects its structure and function. These changes are the primary reasons for the advent of many complex diseases. Analyzing single point mutations in a protein is crucial to see their impact and to understand the disease mechanism. This has given many biophysical resources, including databases and web-based tools to explore the effects of mutations on the structure and function of human proteins. For a given mutation, each tool provides a score-based outcomes which indicate deleterious probability. In recent years, developments in existing programs and the introduction of new prediction algorithms have transformed the state-of-the-art protein mutation analysis. In this study, we have performed a systematic study of the most commonly used mutational analysis programs (10 sequence-based and 5 structure-based) to compare their prediction efficiency. We have carried out extensive mutational analyses using these tools for previously known pathogenic single point mutations of five different proteins. These analyses suggested that sequence-based tools, PolyPhen2, PROVEAN, and PMut, and structure-based web tool, mCSM have a better prediction accuracy. This study indicates that the employment of more than one program based on different approaches should significantly improve the prediction power of the available methods.

Coats plus in prematurity

BACKGROUND

Coats plus syndrome or cerebroretinal microangiopathy with calcifications and cysts (CMCC) is an exceedingly rare autosomal recessive disorder that predominantly affects the microvasculature in the retina, brain, bones, and gastrointestinal system. Unlike Coats disease, CMCC is bilateral and affects multiple organ systems.

MATERIALS AND METHODS

Case report.

RESULTS

We report the case of two brothers with Coats Plus syndrome who presented with variable phenotypic expression. One sibling (Patient 1) was thought to have atypical retinopathy of prematurity and was only diagnosed with Coats plus after his older brother (Patient 2) presented with a seizure and a left upper extremity tremor at 4 years of age. The CTC1 mutation was confirmed in both patients. Aggressive treatment with laser photocoagulation and intravitreal bevacizumab dramatically improved the retinal vascular and exudative changes.

CONCLUSION

Coats Plus syndrome can have a variable phenotypic presentation, including retinal vascular findings. This rare genetic disease should be in the differential diagnosis in patients who present with atypical retinal pathology, including Retinopathy of Prematurity, Familial Exudative Vitreoretinopathy, or Coats disease associated with non-specific multiorgan abnormalities.

Impact of Single Amino Acid Substitutions in Parkinsonism-Associated Deglycase-PARK7 and Their Association with Parkinson's Disease

Parkinsonism-associated deglycase-PARK7/DJ-1 (PARK7) is a multifunctional protein having significant roles in inflammatory and immune disorders and cell protection against oxidative stress. Mutations in PARK7 may result in the onset and progression of a few neurodegenerative disorders such as Parkinson's disease. This study has analyzed the non-synonymous single nucleotide polymorphisms (nsSNPs) resulting in single amino acid substitutions in PARK7 to explore its disease-causing variants and their structural dysfunctions. Initially, we retrieved the mutational dataset of PARK7 from the Ensembl database and performed detailed analyses using sequence-based and structure-based approaches. The pathogenicity of the PARK7 was then performed to distinguish the destabilizing/deleterious variants. Aggregation propensity, noncovalent interactions, packing density, and solvent accessible surface area analyses were carried out on the selected pathogenic mutations. The SODA study suggested that mutations in PARK7 result in aggregation, inducing disordered helix and altering the strand propensity. The effect of mutations alters the number of hydrogen bonds and hydrophobic interactions in PARK7, as calculated from the Arpeggio server. The study indicated that the alteration in the hydrophobic contacts and frustration of the protein could alter the stability of the missense variants of the PARK7, which might result in disease progression. This study provides a detailed understanding of the destabilizing effects of single amino acid substitutions in PARK7.

Comprehensive Characterization of the Coding and Non-Coding Single Nucleotide Polymorphisms in the Tumor Protein p63 (TP63) Gene Using In Silico Tools

Single nucleotide polymorphisms (SNPs) help to understand the phenotypic variations in humans. Genome-wide association studies (GWAS) have identified SNPs located in the tumor protein 63 (TP63) locus to be associated with the genetic susceptibility of cancers. However, there is a lack of in-depth characterization of the structural and functional impacts of the SNPs located at the TP63 gene. The current study was designed for the comprehensive characterization of the coding and non-coding SNPs in the human TP63 gene for their functional and structural significance. The functional and structural effects of the SNPs were investigated using a wide variety of computational tools and approaches, including molecular dynamics (MD) simulation. The deleterious impact of eight nonsynonymous SNPs (nsSNPs) affecting protein stability, structure, and functions was measured by using 13 bioinformatics tools. These eight nsSNPs are in highly conserved positions in protein and were predicted to decrease protein stability and have a deleterious impact on the TP63 protein function. Molecular docking analysis showed five nsSNPs to reduce the binding affinity of TP63 protein to DNA with significant results for three SNPs (R319H, G349E, and C347F). Further, MD simulations revealed the possible disruption of TP63 and DNA binding, hampering the essential protein function. PolymiRTS study found five non-coding SNPs in miRNA binding sites, and the GTEx portal recognized five eQTLs SNPs in single tissue of the lung, heart (LV), and cerebral hemisphere (brain). Characterized nsSNPs and non-coding SNPs will help researchers to focus on TP63 gene loci and ascertain their association with certain diseases.

Impact of Deleterious Mutations on Structure, Function and Stability of Serum/Glucocorticoid Regulated Kinase 1: A Gene to Diseases Correlation

Serum and glucocorticoid-regulated kinase 1 (SGK1) is a Ser/Thr protein kinase involved in regulating cell survival, growth, proliferation, and migration. Its elevated expression and dysfunction are reported in breast, prostate, hepatocellular, lung adenoma, and renal carcinomas. We have analyzed the SGK1 mutations to explore their impact at the sequence and structure level by utilizing state-of-the-art computational approaches. Several pathogenic and destabilizing mutations were identified based on their impact on SGK1 and analyzed in detail. Three amino acid substitutions, K127M, T256A, and Y298A, in the kinase domain of SGK1 were identified and incorporated structurally into original coordinates of SGK1 to explore their time evolution impact using all-atom molecular dynamic (MD) simulations for 200 ns. MD results indicate substantial conformational alterations in SGK1, thus its functional loss, particularly upon T256A mutation. This study provides meaningful insights into SGK1 dysfunction upon mutation, leading to disease progression, including cancer, and neurodegeneration.

Investigating single amino acid substitutions in PIM1 kinase: A structural genomics approach

PIM1, is a serine/threonine proto-oncogene kinase, involved in many biological functions, including cell survival, proliferation, and differentiation, thus play a key role in oncogenesis. It plays a crucial role in the onset and progression of various hematopoietic and non-hematopoietic malignancies, including acute myeloid leukemia and prostate cancer. Mutations in PIM1, especially in its kinase domain, can induce abnormal structural changes and thus alter functionalities that can lead to disease progression and other complexities. Herein, we have performed an extensive analysis of the PIM1 mutations at sequence and structure level while utilizing state-of-the-art computational approaches. Based on the impact on PIM1, numerous pathogenic and destabilizing mutations were identified and subsequently analyzed in detail. Finally, two amino acid substitutions (W109C and F147C) in the kinase domain of PIM1 were selected to explore their impact on the PIM1 structure in a time evolution manner using all-atom molecular dynamics (MD) simulations for 200 ns. MD results indicate significant conformational altercations in the structure of PIM1, especially upon F147C mutation. This study provides a significant insight into the PIM1 dysfunction upon single amino acid substitutions, which can be utilized to get insights into the molecular basis of PIM1-associated disease progression.

Genomic Variations in the Structural Proteins of SARS-CoV-2 and Their Deleterious Impact on Pathogenesis: A Comparative Genomics Approach

A continual rise in severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection causing coronavirus disease (COVID-19) has become a global threat. The main problem comes when SARS-CoV-2 gets mutated with the rising infection and becomes more lethal for humankind than ever. Mutations in the structural proteins of SARS-CoV-2, i.e., the spike surface glycoprotein (S), envelope (E), membrane (M) and nucleocapsid (N), and replication machinery enzymes, i.e., main protease (Mpro) and RNA-dependent RNA polymerase (RdRp) creating more complexities towards pathogenesis and the available COVID-19 therapeutic strategies. This study analyzes how a minimal variation in these enzymes, especially in S protein at the genomic/proteomic level, affects pathogenesis. The structural variations are discussed in light of the failure of small molecule development in COVID-19 therapeutic strategies. We have performed in-depth sequence- and structure-based analyses of these proteins to get deeper insights into the mechanism of pathogenesis, structure-function relationships, and development of modern therapeutic approaches. Structural and functional consequences of the selected mutations on these proteins and their association with SARS-CoV-2 virulency and human health are discussed in detail in the light of our comparative genomics analysis.

Impact of single amino acid substitution on the structure and function of TANK-binding kinase-1

Tank-binding kinase 1 (TBK1) is a serine/threonine protein kinase involved in various signaling pathways and subsequently regulates cell proliferation, apoptosis, autophagy, antiviral and antitumor immunity. Dysfunction of TBK1 can cause many complex diseases, including autoimmunity, neurodegeneration, and cancer. This dysfunction of TBK1 may result from single amino acid substitutions and subsequent structural alterations. This study analyzed the effect of substituting amino acids on TBK1 structure, function, and subsequent disease using advanced computational methods and various tools. In the initial assessment, a total of 467 mutations were found to be deleterious. After that, in detailed structural and sequential analyses, 13 mutations were found to be pathogenic. Finally, based on the functional importance, two variants (K38D and S172A) of the TBK1 kinase domain were selected and studied in detail by utilizing all-atom molecular dynamics (MD) simulation for 200 ns. MD simulation, including correlation matrix and principal component analysis, helps to get deeper insights into the TBK1 structure at the atomic level. We observed a substantial change in variants' conformation, which may be possible for structural alteration and subsequent TBK1 dysfunction. However, substitution S172A shows a significant conformational change in TBK1 structure as compared to K38D. Thus, this study provides a structural basis to understand the effect of mutations on the kinase domain of TBK1 and its function associated with disease progression.