Search

Background: Colorectal cancer (CRC) is the leading cancer among Saudis, and mutations in BRAF , KRAS , and NRAS genes are therapeutically significant due to their association with pathways critical for cell cycle regulation. This study evaluates the prevalence and frequency of somatic mutations in these actionable genes in Saudi CRC patients and assesses their pathogenicity with bioinformatics methods., Methodology: The study employed the TruSight Tumor 15 next-generation sequencing (NGS) panel on 86 colorectal cancer (CRC) samples to detect somatic mutations in BRAF , KRAS , and NRAS genes. Bioinformatic analyses of NGS sequences included variant annotation with ANNOVAR, pathogenicity prediction, variant reclassification with CancerVar, and extensive structural analysis. Additionally, molecular docking assessed the binding of Encorafenib to wild-type and mutant BRAF proteins, providing insights into the therapeutic relevance of pathogenic variants., Results: Out of 86 tumor samples, 40 (46.5%) harbored somatic mutations within actionable genes ( BRAF : 2.3%, KRAS : 43%, NRAS : 2.3%). Fourteen missense variants were identified ( BRAF : n = 1, KRAS : n = 11, NRAS : n = 2). Variants with strong clinical significance included BRAF V600E (2.32%) and KRAS G12D (18.60%). Variants with potential clinical significance included several KRAS and an NRAS mutation, while variants of unknown significance included KRAS E49K and NRAS R102Q. One variant was novel: NRAS R102Q, and two were rare: KRAS E49K and G138E. We further extended the CancerVar prediction capability by adding new pathogenicity prediction tools. Molecular docking demonstrated that Encorafenib inhibits the V600E variant BRAF protein less effectively compared to its wild-type counterpart., Conclusion: Overall, this study highlights the importance of comprehensive molecular screening and bioinformatics in understanding the mutational landscape of CRC in the Saudi population, ultimately improving targeted drug treatments., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024 Aljuhani, Shaik, Alqawas, Bokhary, Al-Mutadares, Al Mahdi, Al-Rayes, El-Harouni, Elango, Banaganapalli and Awan.)

Background: Necrotic enteritis (NE) is an economically important disease of broiler chickens caused by Clostridium perfringens (CP). The pathogenesis, or disease process, of NE is still not clear. This study aimed to identify the alterations of metabolites and metabolic pathways associated with subclinical or clinical NE in CP infected birds and to investigate the possible variations in the metabolic profile of birds infected with different isolates of CP., Methodology: Using a well-established NE model, the protein content of feed was changed abruptly before exposing birds to CP isolates with different toxin genes combinations (cpa, cpb2, netB, tpeL; cpa, cpb2, netB; or cpa, cpb2). Metabolomics analysis of jejunal contents was performed by a targeted, fully quantitative LC-MS/MS based assay., Results: This study detected statistically significant differential expression of 34 metabolites including organic acids, amino acids, fatty acids, and biogenic amines, including elevation of butyric acid at onset of NE in broiler chickens. Subsequent analysis of broilers infected with CP isolates with different toxin gene combinations confirmed an elevation of butyric acid consistently among 21 differentially expressed metabolites including organic acids, amino acids, and biogenic amines, underscoring its potential role during the development of NE. Furthermore, protein-metabolite network analysis revealed significant alterations in butyric acid and arginine-proline metabolisms., Conclusion: This study indicates a significant metabolic difference between CP-infected and non-infected broiler chickens. Among all the metabolites, butyric acid increased significantly in CP-infected birds compared to non-infected healthy broilers. Logistic regression analysis revealed a positive association between butyric acid (coefficient: 1.23, P < 0.01) and CP infection, while showing a negative association with amino acid metabolism. These findings suggest that butyric acid could be a crucial metabolite linked to the occurrence of NE in broiler chickens and may serve as an early indicator of the disease at the farm level. Further metabolomic experiments using different NE animal models and field studies are needed to determine the specificity and to validate metabolites associated with NE, regardless of predisposing factors., (© 2024. The Author(s).)

Endometrial Cancer (EC) is one of the most common gynecological malignancies. Despite its prevalence, molecular pathways, such as the Sonic Hedgehog (SHH) pathway, have not been extensively studied in the context of EC. This study aims to explore the clinical implications of SHH expression in EC, potentially uncovering new insights into the disease's pathogenesis and offering valuable insights for therapeutic strategies in EC. We utilized data from The Cancer Genome Atlas (TCGA) to divide the dataset into ' High SHH ' and ' Low SHH ' groups based on a gene signature score derived from SHH pathway-related genes. We explored the clinical and tumor characteristics of these groups, focusing on key cancer hallmarks, including stemness, proliferation, cytolytic activity, tumor micro-environment, and genomic instability. ' High SHH ' tumors emerged as a distinct category with favorable clinical and molecular features. These tumors exhibited lower proliferation rates, reduced angiogenesis, and diminished genomic instability, indicating a controlled and less aggressive tumor growth pattern. Moreover, ' High SHH ' tumors displayed lower stemness, highlighting a less invasive phenotype. The immune micro-environment in ' High SHH ' tumors was enriched with immune cell types, such as macrophage M0, monocytes, B cells, CD8 T cells, CD4 T cells, follicular helper T cells, and natural killer cells. This immune enrichment, coupled with higher cytolytic activity, suggested an improved anti-tumor immune response. Our study sheds light on the clinical significance of Sonic signaling in EC. ' High SHH ' tumors exhibit a unique molecular and clinical profile associated with favorable cancer hallmarks, lower grades, and better survival. These findings underscore the potential utility of SHH expression as a robust prognostic biomarker, offering valuable insights for tailored therapeutic strategies in EC. Understanding the SHH pathway's role in EC contributes to our growing knowledge of this cancer and may pave the way for more effective treatment strategies in the future.

Background: Familial hypercholesterolemia (FH) is a globally underdiagnosed inherited metabolic disorder. Owing to limited published data from Arab world, this study was conducted with the aim of identifying the genetic and molecular basis of FH in highly consanguineous Saudi population., Methods: We performed clinical screening, biochemical profiling, whole exome sequencing and variant segregation analysis of two Saudi FH families. Additionally, 500 normolipic individuals were screened to ensure the absence of FH variant in general Saudi population. Functional characterization of FH variants on secondary structure characteristics of RNA and protein molecules was performed using different bioinformatics modelling approaches., Results: WES analysis identified two independent rare LDLR gene stop gain variants (p.C231* and p.R744*) consistent to the clinical presentation of FH patients from two different families. RNAfold analysis has shown that both variants were predicted to disturb the free energy dynamics of LDLR mRNA molecule and destabilize its folding pattern and function. PSIPRED based structural modelling analysis has suggested that both variants bring drastic changes disturbing the secondary structural elements of LDLR molecule. The p.C231* and p.R744* variants are responsible for partial or no protein product, thus they are class 1 variants causing loss of function (LoF) LDLR variants., Conclusions: This study highlights the effectiveness of the WES, sanger sequencing, and computational analysis in expanding FH variant spectrum in culturally distinct populations like Saudi Arabia. Genetic testing of FH patients is very essential in better clinical diagnosis, screening, treatment, and management and prevention of cardiovascular disease burden in the society.

Familial hypercholesterolemia (FH) is a globally underdiagnosed genetic condition associated with premature cardiovascular death. The genetic etiology data on Arab FH patients is scarce. Therefore, this study aimed to identify the genetic basis of FH in a Saudi family using whole exome sequencing (WES) and multidimensional bioinformatic analysis. Our WES findings revealed a rare heterozygous gain-of-function variant (R496W) in the exon 9 of the PCSK9 gene as a causal factor for FH in this family. This variant was absent in healthy relatives of the proband and 200 healthy normolipidemic controls from Saudi Arabia. Furthermore, this variant has not been previously reported in various regional and global population genomic variant databases. Interestingly, this variant is classified as "likely pathogenic" (PP5) based on the variant interpretation guidelines of the American College of Medical Genetics (ACMG). Computational functional characterization suggested that this variant could destabilize the native PCSK9 protein and alter its secondary and tertiary structural features. In addition, this variant was predicted to negatively influence its ligand-binding ability with LDLR and Alirocumab antibody molecules. This rare PCSK9 (R496W) variant is likely to expand our understanding of the genetic basis of FH in Saudi Arabia. This study also provides computational structural insights into the genotype-protein phenotype relationship of PCSK9 pathogenic variants and contributes to the development of personalized medicine for FH patients in the future., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Shaik, Al-Shehri, Athar, Awan, Khalili, Al Mahadi, Hejazy, Saadah, Al-Harthi, Elango, Banaganapalli, Alefishat and Awan.)

Background: Inflammatory bowel disease (IBD) is a chronic autoimmune disorder characterized by severe inflammation and mucosal destruction of the intestine. The specific, complex molecular processes underlying IBD pathogenesis are not well understood. Therefore, this study is aimed at identifying and uncovering the role of key genetic factors in IBD., Method: The whole exome sequences (WESs) of three consanguineous Saudi families having many siblings with IBD were analyzed to discover the causal genetic defect. Then, we used a combination of artificial intelligence approaches, such as functional enrichment analysis using immune pathways and a set of computational functional validation tools for gene expression, immune cell expression analyses, phenotype aggregation, and the system biology of innate immunity, to highlight potential IBD genes that play an important role in its pathobiology., Results: Our findings have shown a causal group of extremely rare variants in the LILRB1 (Q53L, Y99N, W351G, D365A, and Q376H) and PRSS3 (F4L and V25I) genes in IBD-affected siblings. Findings from amino acids in conserved domains, tertiary-level structural deviations, and stability analysis have confirmed that these variants have a negative impact on structural features in the corresponding proteins. Intensive computational structural analysis shows that both genes have very high expression in the gastrointestinal tract and immune organs and are involved in a variety of innate immune system pathways. Since the innate immune system detects microbial infections, any defect in this system could lead to immune functional impairment contributing to IBD., Conclusion: The present study proposes a novel strategy for unraveling the complex genetic architecture of IBD by integrating WES data of familial cases, with computational analysis., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Jan, Al-Numan, Al-Twaty, Alrayes, Alsufyani, Alaifan, Alhussaini, Shaik, Awan, Qari, Saadah, Banaganapalli, Mosli and Elango.)

The Parkinson disease (PD) is the second most common neurodegenerative disorder affecting the central nervous system and motor functions. The biological complexity of PD is yet to reveal potential targets for intervention or to slow the disease severity. Therefore, this study aimed to compare the fidelity of blood to substantia nigra (SN) tissue gene expression from PD patients to provide a systematic approach to predict role of the key genes of PD pathobiology. Differentially expressed genes (DEGs) from multiple microarray data sets of PD blood and SN tissue from GEO database are identified. Using the theoretical network approach and variety of bioinformatic tools, we prioritized the key genes from DEGs. A total of 540 and 1024 DEGs were identified in blood and SN tissue samples, respectively. Functional pathways closely related to PD such as ERK1 and ERK2 cascades, mitogen-activated protein kinase (MAPK) signaling, Wnt, nuclear factor-κB (NF-κB), and PI3K-Akt signaling were observed by enrichment analysis. Expression patterns of 13 DEGs were similar in both blood and SN tissues. Comprehensive network topological analysis and gene regulatory networks identified additional 10 DEGs functionally connected with molecular mechanisms of PD through the mammalian target of rapamycin (mTOR), autophagy, and AMP-activated protein kinase (AMPK) signaling pathways. Potential drug molecules were identified by chemical-protein network and drug prediction analysis. These potential candidates can be further validated in vitro/in vivo to be used as biomarkers and/or novel drug targets for the PD pathology and/or to arrest or delay the neurodegeneration over the years, respectively., Competing Interests: The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article., (© The Author(s) 2023.)

Paget's disease of bone (PDB) is the second most prevalent metabolic bone disorder worldwide, with a prevalence rate of 1.5%-8.3%. It is characterized by localized areas of accelerated, disorganized, and excessive bone production and turnover. Typically, PDB develops in the later stages of life, particularly in the late 50s, and affects men more frequently than women. PDB is a complex disease influenced by both genetic and environmental factors. PDB has a complex genetic basis involving multiple genes, with SQSTM1 being the gene most frequently associated with its development. Mutations affecting the UBA domain of SQSTM1 have been detected in both familial and sporadic PDB cases, and these mutations are often associated with severe clinical expression. Germline mutations in other genes such as TNFRSF11A, ZNF687 and PFN1, have also been associated with the development of the disease. Genetic association studies have also uncovered several PDB predisposing risk genes contributing to the disease pathology and severity. Epigenetic modifications of genes involved in bone remodelling and regulation, including RANKL, OPG, HDAC2, DNMT1, and SQSTM1, have been implicated in the development and progression of Paget's disease of bone, providing insight into the molecular basis of the disease and potential targets for therapeutic intervention. Although PDB has a tendency to cluster within families, the variable severity of the disease across family members, coupled with decreasing incidence rates, indicates that environmental factors may also play a role in the pathophysiology of PDB. The precise nature of these environmental triggers and how they interact with genetic determinants remain poorly understood. Fortunately, majority of PDB patients can achieve long-term remission with an intravenous infusion of aminobisphosphonates, such as zoledronic acid. In this review, we discuss aspects like clinical characteristics, genetic foundation, and latest updates in PDB research., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Banaganapalli, Fallatah, Alsubhi, Shetty, Awan, Elango and Shaik.)

Introduction: PIM kinases are targets for therapeutic intervention since they are associated with a number of malignancies by boosting cell survival and proliferation. Over the past years, the rate of new PIM inhibitors discovery has increased significantly, however, new generation of potent molecules with the right pharmacologic profiles were in demand that can probably lead to the development of Pim kinase inhibitors that are effective against human cancer. Method: In the current study, a machine learning and structure based approaches were used to generate novel and effective chemical therapeutics for PIM-1 kinase. Four different machine learning methods, namely, support vector machine, random forest, k-nearest neighbour and XGBoost have been used for the development of models. Total, 54 Descriptors have been selected using the Boruta method. Results: SVM, Random Forest and XGBoost shows better performance as compared to k-NN. An ensemble approach was implemented and, finally, four potential molecules (CHEMBL303779, CHEMBL690270, MHC07198, and CHEMBL748285) were found to be effective for the modulation of PIM-1 activity. Molecular docking and molecular dynamic simulation corroborated the potentiality of the selected molecules. The molecular dynamics (MD) simulation study indicated the stability between protein and ligands. Discussion: Our findings suggest that the selected models are robust and can be potentially useful for facilitating the discovery against PIM kinase., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Almukadi, Jadkarim, Mohammed, Almansouri, Sultana, Shaik and Banaganapalli.)

Background: Tetralogy of Fallot (TOF) is a rare, complex congenital heart defect caused by genetic and environmental interactions that results in abnormal heart development during the early stages of pregnancy. Genetic basis of TOF in Saudi populations is not yet studied. Therefore, the objective of this study is to screen for the molecular defects causing TOF in Saudi patients., Methods: A family with non-syndromic TOF was recruited from the Western region of Saudi Arabia. Whole exome sequencing (WES) was performed on the proband and her parents. The identified candidate variant was verified by sanger sequencing. Also, different computational biology tools were used to figure out how candidate variants affect the structure and function of candidate protein involved in TOF., Results: A novel heterozygous de novo mutation in LRP1 (p. G3311D) gene was identified in the index case. Also, this variant was absent in the in-house exome sequencing data of 80 healthy Saudi individuals. This variant was predicted to be likely pathogenic, as it negatively affects the biophysical chemical properties and stability of the protein. Furthermore, functional biology data from knock out mouse models confirms that molecular defects in LRP1 gene leads to cardiac defects and lethality. This variant was not previously reported in both Arab and global population genetic databases., Conclusion: The findings in this study postulate that the LRP1 variant has a role in TOF pathogenesis and facilitate accurate diagnosis as well as the understanding of underlying molecular mechanisms and pathophysiology of the disease., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 Elsevier B.V. All rights reserved.)

Endometrial cancer (EC) is a urogenital cancer affecting millions of post-menopausal women, globally. This study aims to identify key miRNAs, target genes, and drug targets associated with EC metastasis. The global miRNA and mRNA expression datasets of endometrial tissue biopsies (24 tumors +3 healthy tissues for mRNA and 18 tumor +4 healthy tissues for miRNAs), were extensively analyzed by mapping of DEGs, DEMi, biological pathway enrichment, miRNA-mRNA networking, drug target identification, and survival curve output for differentially expressed genes. Our results reveal the dysregulated expression of 26 miRNAs and their 66 target genes involved in focal adhesions, p53 signaling pathway, ECM-receptor interaction, Hedgehog signaling pathway, fat digestion and absorption, glioma as well as retinol metabolism involved in cell growth, migration, and proliferation of endometrial cancer cells. The subsequent miRNA-mRNA network and expression status analysis have narrowed down to 2 hub miRNAs (hsa-mir-200a, hsa-mir-429) and 6 hub genes ( PTCH1, FOSB, PDGFRA, CCND2, ABL1, ALDH1A1 ). Further investigations with different systems biology methods have prioritized ALDH1A1, ABL1 and CCND2 as potential genes involved in endometrial cancer metastasis owing to their high mutation load and expression status. Interestingly, overexpression of PTCH1 , ABL1 and FOSB genes are reported to be associated with a low survival rate among cancer patients. The upregulated hsa-mir-200a-b is associated with the decreased expression of the PTCH1, CCND2, PDGFRA, FOSB and ABL1 genes in endometrial cancer tissue while hsa-mir-429 is correlated with the decreased expression of the ALDH1A1 gene, besides some antibodies, PROTACs and inhibitory molecules. In conclusion, this study identified key miRNAs (hsa-mir-200a, hsa-mir-429) and target genes ALDH1A1 , ABL1 and CCND2 as potential biomarkers for metastatic endometrial cancers from large-scale gene expression data using systems biology approaches., (Copyright © 2023 Ajabnoor, Alsubhi, Shinawi, Habhab, Albaqami, Alqahtani, Nasief, Bondagji, Elango, Shaik and Banaganapalli.)

Background: Coronavirus disease (COVID-19) infection is known for its severe clinical pathogenesis among individuals with pre-existing comorbidities. However, the molecular basis of this observation remains elusive. Thus, this study aimed to map key genes and pathway alterations in patients with COVID-19 and comorbidities using robust systems biology approaches. Methods: The publicly available genome-wide transcriptomic datasets from 120 COVID-19 patients, 281 patients suffering from different comorbidities (like cardiovascular diseases, atherosclerosis, diabetes, and obesity), and 252 patients with different infectious diseases of the lung (respiratory syncytial virus, influenza, and MERS) were studied using a range of systems biology approaches like differential gene expression, gene ontology (GO), pathway enrichment, functional similarity, mouse phenotypic analysis and drug target identification. Results: By cross-mapping the differentially expressed genes (DEGs) across different datasets, we mapped 274 shared genes to severe symptoms of COVID-19 patients or with comorbidities alone. GO terms and functional pathway analysis highlighted genes in dysregulated pathways of immune response, interleukin signaling, FCGR activation, regulation of cytokines, chemokines secretion, and leukocyte migration. Using network topology parameters, phenotype associations, and functional similarity analysis with ACE2 and TMPRSS2 -two key receptors for this virus-we identified 17 genes with high connectivity ( CXCL10, IDO1, LEPR, MME, PTAFR, PTGS2, MAOB, PDE4B, PLA2G2A, COL5A1, ICAM1, SERPINE1, ABCB1, IL1R1, ITGAL, NCAM1 and PRKD1 ) potentially contributing to the clinical severity of COVID-19 infection in patients with comorbidities. These genes are predicted to be tractable and/or with many existing approved inhibitors, modulators, and enzymes as drugs. Conclusion: By systemic implementation of computational methods, this study identified potential candidate genes and pathways likely to confer disease severity in COVID-19 patients with pre-existing comorbidities. Our findings pave the way to develop targeted repurposed therapies in COVID-19 patients., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Mujalli, Alghamdi, Nasser, Al-Rayes, Banaganapalli, Shaik and Elango.)

Background: Alpha-1 antitrypsin deficiency (A1ATD) is a progressive lung disease caused by inherited pathogenic variants in the SERPINA1 gene. However, their actual role in maintenance of structural and functional characteristics of the corresponding α-1 anti-trypsin (A1AT) protein is not well characterized. Methods: The A1ATD causative SERPINA1 missense variants were initially collected from variant databases, and they were filtered based on their pathogenicity potential. Then, the tertiary protein models were constructed and the impact of individual variants on secondary structure, stability, protein-protein interactions, and molecular dynamic (MD) features of the A1AT protein was studied using diverse computational methods. Results: We identified that A1ATD linked SERPINA1 missense variants like F76S, S77F, L278P, E288V, G216C, and H358R are highly deleterious as per the consensual prediction scores of SIFT, PolyPhen, FATHMM, M-CAP and REVEL computational methods. All these variants were predicted to alter free energy dynamics and destabilize the A1AT protein. These variants were seen to cause minor structural drifts at residue level (RMSD = <2Å) of the protein. Interestingly, S77F and L278P variants subtly alter the size of secondary structural elements like beta pleated sheets and loops. The residue level fluctuations at 100 ns simulation confirm the highly damaging structural consequences of all the six missense variants on the conformation dynamics of the A1AT protein. Moreover, these variants were also predicted to cause functional deformities by negatively impacting the binding energy of A1AT protein with NE ligand molecule. Conclusion: This study adds a new computational biology dimension to interpret the genotype-protein phenotype relationship between S ERPINA1 pathogenic variants with its structural plasticity and functional behavior with NE ligand molecule contributing to the Alpha-1-antitrypsin deficiency. Our results support that A1ATD complications correlates with the conformational flexibility and its propensity of A1AT protein polymerization when misfolded., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Shaik, Saud Al-Saud, Abdulhamid Aljuhani, Jamil, Alnuman, Aljeaid, Sultana, El-Harouni, Awan, Elango and Banaganapalli.)

Background: Prostate cancer (PC) is a fatally aggressive urogenital cancer killing millions of men, globally. Thus, this study aims to identify key miRNAs, target genes, and drug targets associated with prostate cancer metastasis. Methods: The miRNA and mRNA expression datasets of 148 prostate tissue biopsies (39 tumours and 109 normal tissues), were analysed by differential gene expression analysis, protein interactome mapping, biological pathway analysis, miRNA-mRNA networking, drug target analysis, and survival curve analysis. Results: The dysregulated expression of 53 miRNAs and their 250 target genes involved in Hedgehog, ErbB, and cAMP signalling pathways connected to cell growth, migration, and proliferation of prostate cancer cells was detected. The subsequent miRNA-mRNA network and expression status analysis have helped us in narrowing down their number to 3 hub miRNAs (hsa-miR-455-3p, hsa-miR-548c-3p, and hsa-miR-582-5p) and 9 hub genes ( NFIB, DICER1, GSK3B, DCAF7, FGFR1OP, ABHD2, NACC2, NR3C1 , and FGF2 ). Further investigations with different systems biology methods have prioritized NR3C1, ABHD2, and GSK3B as potential genes involved in prostate cancer metastasis owing to their high mutation load and expression status. Interestingly, down regulation of NR3C1 seems to improve the prostate cancer patient survival rate beyond 150 months. The NR3C1, ABHD2, and GSK3B genes are predicted to be targeted by hsa-miR-582-5p, besides some antibodies, PROTACs and inhibitory molecules. Conclusion: This study identified key miRNAs (miR-548c-3p and miR-582-5p) and target genes ( NR3C1, ABHD2, and GSK3B ) as potential biomarkers for metastatic prostate cancers from large-scale gene expression data using systems biology approaches., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Shinawi, Nasser, Moradi, Mujalli, Albaqami, Almukadi, Elango, Shaik and Banaganapalli.)

Asthma is a life-threatening and chronic inflammatory lung disease that is posing a true global health challenge. The genetic basis of the disease is fairly well examined. However, the molecular crosstalk between microRNAs (miRNAs), target genes, and transcription factors (TFs) networks and their contribution to disease pathogenesis and progression is not well explored. Therefore, this study was aimed at dissecting the molecular network between mRNAs, miRNAs, and TFs using robust computational biology approaches. The transcriptomic data of bronchial epithelial cells of severe asthma patients and healthy controls was studied by different systems biology approaches like differentially expressed gene detection, functional enrichment, miRNA-target gene pairing, and mRNA-miRNA-TF molecular networking. We detected the differential expression of 1703 (673 up-and 1030 down-regulated) genes and 71 (41 up-and 30 down-regulated) miRNAs in the bronchial epithelial cells of asthma patients. The DEGs were found to be enriched in key pathways like IL-17 signaling (KEGG: 04657), Th1 and Th2 cell differentiation (KEGG: 04658), and the Th17 cell differentiation (KEGG: 04659) (p-values = 0.001). The results from miRNAs-target gene pairs-transcription factors (TFs) have detected the key roles of 3 miRs (miR-181a-2-3p; miR-203a-3p; miR-335-5p), 6 TFs (TFAM, FOXO1, GFI1, IRF2, SOX9, and HLF) and 32 miRNA target genes in eliciting autoimmune reactions in bronchial epithelial cells of the respiratory tract. Through systemic implementation of comprehensive system biology tools, this study has identified key miRNAs, TFs, and miRNA target gene pairs as potential tissue-based asthma biomarkers., Competing Interests: The authors have declared that no competing interests exist.

Chronic obstructive pulmonary disease (COPD) is a multifactorial progressive airflow obstruction in the lungs, accounting for high morbidity and mortality across the world. This study aims to identify potential COPD blood-based biomarkers by analyzing the dysregulated gene expression patterns in blood and lung tissues with the help of robust computational approaches. The microarray gene expression datasets from blood (136 COPD and 6 controls) and lung tissues (16 COPD and 19 controls) were analyzed to detect shared differentially expressed genes (DEGs). Then these DEGs were used to construct COPD protein network-clusters and functionally enrich them against gene ontology annotation terms. The hub genes in the COPD network clusters were then queried in GWAS catalog and in several cancer expression databases to explore their pathogenic roles in lung cancers. The comparison of blood and lung tissue datasets revealed 63 shared DEGs. Of these DEGs, 12 COPD hub gene-network clusters (SREK1, TMEM67, IRAK2, MECOM, ASB4, C1QTNF2, CDC42BPA, DPF3, DET1, CCDC74B, KHK, and DDX3Y) connected to dysregulations of protein degradation, inflammatory cytokine production, airway remodeling, and immune cell activity were prioritized with the help of protein interactome and functional enrichment analysis. Interestingly, IRAK2 and MECOM hub genes from these COPD network clusters are known for their involvement in different pulmonary diseases. Additional COPD hub genes like SREK1, TMEM67, CDC42BPA, DPF3, and ASB4 were identified as prognostic markers in lung cancer, which is reported in 1% of COPD patients. This study identified 12 gene network- clusters as potential blood based genetic biomarkers for COPD diagnosis and prognosis., Competing Interests: The authors have declared that no competing interests exist.

Background: Autoimmune diseases (AIDs) share a common molecular etiology and often present overlapping clinical presentations. Thus, this study aims to explore the complex molecular basis of AID by whole exome sequencing and computational biology analysis., Methods: Molecular screening of the consanguineous AID family and the computational biology characterization of the potential variants were performed. The potential variants were searched against the exome data of 100 healthy individuals and 30 celiac disease patients., Result: A complex inheritance pattern of PAK2 (V43A), TAP2 (F468Y), and PLCL1 (V473I) genetic variants was observed in the three probands of the AID family. The PAK2 variant (V43A) is a novel one, but TAP2 (F468Y) and PLCL1 (V473I) variants are extremely rare in local Arab (SGHP and GME) and global (gnomAD) databases. All these variants were localized in functional domains, except for the PAK2 variant (V43A) and were predicted to alter the structural (secondary structure elements, folding, active site confirmation, stability, and solvent accessibility) and functional (gene expression) features. Therefore, it is reasonable to postulate that the dysregulation of PAK2 , TAP2 , and PLCL1 genes is likely to elicit autoimmune reactions by altering antigen processing and presentation, T cell receptor signaling, and immunodeficiency pathways., Conclusion: Our findings highlight the importance of exploring the alternate inheritance patterns in families presenting complex autoimmune diseases, where classical genetic models often fail to explain their molecular basis. These findings may have potential implications for developing personalized therapies for complex disease patients., Competing Interests: The authors declare that the research was constructed in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Alharthi, Banaganapalli, Hassan, Rashidi, Al-Shehri, Alaifan, Alhussaini, Alsufyani, Alghamdi, Nasser, Bin-Taleb, Elango, Shaik and Saadah.)

Background: Obesity is associated with the quantitative changes in miRNAs and their target genes. However, the molecular basis of their dysregulation and expression status correlations is incompletely understood. Therefore, this study aims to examine the shared differentially expressed miRNAs and their target genes between blood and adipose tissues of obese individuals to identify potential blood-based biomarkers., Methods: In this study, 3 gene expression datasets (two mRNA and one miRNA), generated from blood and adipose tissues of 68 obese and 39 lean individuals, were analyzed by a series of robust computational concepts, like protein interactome mapping, functional enrichment of biological pathways and construction of miRNA-mRNA and transcription factor gene networks., Results: The comparison of blood versus tissue datasets has revealed the shared differential expression of 210 genes (59.5% upregulated) involved in lipid metabolism and inflammatory reactions. The blood miRNA (GSE25470) analysis has identified 79 differentially expressed miRNAs (71% downregulated). The miRNA-target gene scan identified regulation of 30 shared genes by 22miRNAs. The gene network analysis has identified the inverse expression correlation between 8 target genes (TP53, DYSF, GAB2, GFRA2, NACC2, FAM53C, JNK and GAB2) and 3 key miRNAs (hsa-mir-940, hsa-mir-765, hsa-mir-612), which are further regulated by 24 key transcription factors., Conclusions: This study identifies potential obesity related blood biomarkers from large-scale gene expression data by computational miRNA-target gene interactome and transcription factor network construction methods.

Background: Molecular diagnosis of early onset inflammatory bowel disease (IBD) is very important for adopting suitable treatment strategies. Owing to the sparse data available, this study aims to identify the molecular basis of early onset IBD in Arab patients., Methods: A consanguineous Arab family with monozygotic twins presenting early onset IBD was screened by whole exome sequencing (WES). The variants functional characterization was performed by a series of computational biology methods. The IBD variants were further screened in in-house whole exome data of 100 Saudi cohorts ensure their rare prevalence in the population., Results: Genetic screening has identified the digenic autosomal recessive mode of inheritance of ITGAV (G58V) and FN1 (G313V) variants in IBD twins with early onset IBD. Findings from pathogenicity predictions, stability and molecular dynamics have confirmed the deleterious nature of both variants on structural features of the corresponding proteins. Functional biology data suggested that both genes show abundant expression in gastrointestinal tract and immune organs, involved in immune cell restriction, regulation of different immune related pathways. Data from knockout mouse models for ITGAV gene has revealed that the dysregulated expression of this gene impacts intestinal immune homeostasis. The defective ITGAV and FN1 involved in integrin pathway, are likely to induce intestinal inflammation by disturbing immune homeostasis., Conclusions: Our findings provide novel insights into the molecular etiology of pediatric onset IBD and may likely pave way in developing genomic medicine., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Al-Numan, Jan, Al-Saud, Rashidi, Alrayes, Alsufyani, Mujalli, Shaik, Mosli, Elango, Saadah and Banaganapalli.)

Familial hypercholesterolemia (FH) is a monogenic lipid disorder which promotes atherosclerosis and cardiovascular diseases. Owing to the lack of sufficient published information, this study aims to identify the potential genetic biomarkers for FH by studying the global gene expression profile of blood cells. The microarray expression data of FH patients and controls was analyzed by different computational biology methods like differential expression analysis, protein network mapping, hub gene identification, functional enrichment of biological pathways, and immune cell restriction analysis. Our results showed the dysregulated expression of 115 genes connected to lipid homeostasis, immune responses, cell adhesion molecules, canonical Wnt signaling, mucin type O-glycan biosynthesis pathways in FH patients. The findings from expanded protein interaction network construction with known FH genes and subsequent Gene Ontology (GO) annotations have also supported the above findings, in addition to identifying the involvement of dysregulated thyroid hormone and ErbB signaling pathways in FH patients. The genes like CSNK1A1, JAK3, PLCG2, RALA, and ZEB2 were found to be enriched under all GO annotation categories. The subsequent phenotype ontology results have revealed JAK3I , PLCG2, and ZEB2 as key hub genes contributing to the inflammation underlying cardiovascular and immune response related phenotypes. Immune cell restriction findings show that above three genes are highly expressed by T-follicular helper CD4 + T cells, naïve B cells, and monocytes, respectively. These findings not only provide a theoretical basis to understand the role of immune dysregulations underlying the atherosclerosis among FH patients but may also pave the way to develop genomic medicine for cardiovascular diseases., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2022 The Author(s).)

Celiac disease (CeD) is a multifactorial autoimmune enteropathy characterized by the overactivation of the immune system in response to dietary gluten. The molecular etiology of CeD is still not well-understood. Therefore, this study aims to identify potential candidate genes involved in CeD pathogenesis by applying multilayered system biology approaches. Initially, we identified rare coding variants shared between the affected siblings in two rare Arab CeD families by whole-exome sequencing (WES). Then we used the STRING database to construct a protein network of rare variants and genome-wide association study (GWAS) loci to explore their molecular interactions in CeD. Furthermore, the hub genes identified based on network topology parameters were subjected to a series of computational validation analyses like pathway enrichment, gene expression, knockout mouse model, and variant pathogenicity predictions. Our findings have shown the absence of rare variants showing classical Mendelian inheritance in both families. However, interactome analysis of rare WES variants and GWAS loci has identified a total of 11 hub genes. The multidimensional computational analysis of hub genes has prioritized IL1R1 for family A and CD3E for family B as potential genes. These genes were connected to CeD pathogenesis pathways of T-cell selection, cytokine signaling, and adaptive immune response. Future multi-omics studies may uncover the roles of IL1R1 and CD3E in gluten sensitivity. The present investigation lays forth a novel approach integrating next-generation sequencing (NGS) of familial cases, GWAS, and computational analysis for solving the complex genetic architecture of CeD., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Mansour, Banaganapalli, Nasser, Al-Aama, Shaik, Saadah and Elango.)

Obesity and type 2 and diabetes mellitus (T2D) are two dual epidemics whose shared genetic pathological mechanisms are still far from being fully understood. Therefore, this study is aimed at discovering key genes, molecular mechanisms, and new drug targets for obesity and T2D by analyzing the genome wide gene expression data with different computational biology approaches. In this study, the RNA-sequencing data of isolated primary human adipocytes from individuals who are lean, obese, and T2D was analyzed by an integrated framework consisting of gene expression, protein interaction network (PIN), tissue specificity, and druggability approaches. Our findings show a total of 1932 unique differentially expressed genes (DEGs) across the diabetes versus obese group comparison (p≤0.05). The PIN analysis of these 1932 DEGs identified 190 high centrality network (HCN) genes, which were annotated against 3367 GO terms and functional pathways, like response to insulin signaling, phosphorylation, lipid metabolism, glucose metabolism, etc. (p≤0.05). By applying additional PIN and topological parameters to 190 HCN genes, we further mapped 25 high confidence genes, functionally connected with diabetes and obesity traits. Interestingly, ERBB2, FN1, FYN, HSPA1A, HBA1 , and ITGB1 genes were found to be tractable by small chemicals, antibodies, and/or enzyme molecules. In conclusion, our study highlights the potential of computational biology methods in correlating expression data to topological parameters, functional relationships, and druggability characteristics of the candidate genes involved in complex metabolic disorders with a common etiological basis.

Laterality defects (LDs) or asymmetrically positioned organs are a group of rare developmental disorders caused by environmental and/or genetic factors. However, the exact molecular pathophysiology of LD is not yet fully characterised. In this context, studying Arab population presents an ideal opportunity to discover the novel molecular basis of diseases owing to the high rate of consanguinity and genetic disorders. Therefore, in the present study, we studied the molecular basis of LD in Arab patients, using next-generation sequencing method. We discovered an extremely rare novel missense variant in MYO1D gene (Pro765Ser) presenting with visceral heterotaxy and left isomerism with polysplenia syndrome. The proband in this index family has inherited this homozygous variant from her heterozygous parents following the autosomal recessive pattern. This is the first report to show MYO1D genetic variant causing left-right axis defects in humans, besides previous known evidence from zebrafish, frog and Drosophila models. Moreover, our multilevel bioinformatics-based structural (protein variant structural modelling, divergence, and stability) analysis has suggested that Ser765 causes minor structural drifts and stability changes, potentially affecting the biophysical and functional properties of MYO1D protein like calmodulin binding and microfilament motor activities. Functional bioinformatics analysis has shown that MYO1D is ubiquitously expressed across several human tissues and is reported to induce severe phenotypes in knockout mouse models. In conclusion, our findings show the expanded genetic spectrum of LD, which could potentially pave way for the novel drug target identification and development of personalised medicine for high-risk families., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Alsafwani, Nasser, Shinawi, Banaganapalli, ElSokary, Zaher, Shaik, Abdelmohsen, Al-Aama, Shapiro, O. Al-Radi, Elango and Alahmadi.)

Lamellar ichthyosis (LI) is a rare inherited disease where affected infants present a extensive skin scaling characterized by hyperkeratosis. Inherited mutations in the Transglutaminase 1 (TGM1) protein is one of the known causative genetic factor for the LI. The main objective of this study is to explore the impact of LI causative missense mutations on the structural and stability aspects of TGM1 protein using structural modeling, molecular docking and molecular dynamics approaches. By testing all LI causative TMG1 mutations against multiple stability prediction methods, we found that L362R and L388P mutations positioned in the Transglut&#95;core domain were most destabilizing to the stability of TGM1 protein. These 2 mutations were 3D protein modeled and further analyzed by molecular docking and dynamic simulation methods. Molecular docking of these TGM1 mutant structures with chitosan, a natural polyphenolic compound and known inducer for transglutaminase enzyme, has shown stable molecular interactions between the native TGM1-chitosan and TGM1(L388P)-chitosan complex, when compared to the TGM1(L362R)-chitosan complex. Interestingly, molecular dynamics analysis have also yielded similar findings, where L388P-chitosan complex is shown to develop B-sheets and attain better stability, whereas TGM1-L362R complex possessed coils over the simulation period, pointing its highly destabilizing behavior on the protein structure. This study concludes that missense mutations in Transglut&#95;core domain of the TGM1 are deleterious to the stability and structural changes of TGM1 protein and also suggest that chitosan molecule could act as a natural activator against few pathogenic TGM1 mutations. Communicated by Ramaswamy H. Sarma.

Background: The spread of a novel severe acute respiratory syndrome corona virus 2 (SARS-CoV-2) has affected both the public health and the global economy. The current study was aimed at analysing the genetic sequence of this highly contagious corona virus from an evolutionary perspective, comparing the genetic variation features of different geographic strains, and identifying the key miRNAs as well as their gene targets from the transcriptome data of infected lung tissues., Methods: A multilevel robust computational analysis was undertaken for viral genetic sequence alignment, phylogram construction, genome-wide transcriptome data interpretation of virus-infected lung tissues, miRNA mapping, and functional biology networking., Results: Our findings show both genetic similarities as well as notable differences in the S protein length among SARS-CoV-1, SARS-CoV-2 and MERS viruses. All SARS-CoV-2 strains showed a high genetic similarity with the parent Wuhan strain, but Saudi Arabian, South African, USA, Russia and New Zealand strains carry 3 additional genetic variations like P333L (RNA -dependant RNA polymerase), D614G (spike), and P4715L (ORF1ab). The infected lung tissues demonstrated the upregulation of 282 (56.51%) antiviral defensive response pathway genes and downregulation of 217 (43.48%) genes involved in autophagy and lung repair pathways. By miRNA mapping, 4 key miRNAs (hsa-miR-342-5p, hsa-miR-432-5p, hsa-miR-98-5p and hsa-miR-17-5p), targeting multiple host genes (MYC, IL6, ICAM1 and VEGFA) as well as SARS-CoV2 gene (ORF1ab) were identified., Conclusion: Systems biology methods offer a new perspective in understanding the molecular basis for the faster spread of SARS-CoV-2 infection. The antiviral miRNAs identified in this study may aid in the ongoing search for novel personalized therapeutic avenues for COVID patients., (Copyright © 2021 Elsevier Ltd. All rights reserved.)

Familial hypercholesterolemia (FH), a well-known lipid disease caused by inherited genetic defects in cholesterol uptake and metabolism is underdiagnosed in many countries including Saudi Arabia. The present study aims to identify the molecular basis of severe clinical manifestations of FH patients from unrelated Saudi consanguineous families. Two Saudi families with multiple FH patients fulfilling the combined FH diagnostic criteria of Simon Broome Register, and the Dutch Lipid Clinic Network (DLCN) were recruited. LipidSeq, a targeted resequencing panel for monogenic dyslipidemias, was used to identify causative pathogenic mutation in these two families and in 92 unrelated FH cases. Twelve FH patients from two unrelated families were sharing a very rare, pathogenic and founder LDLR stop gain mutation i.e., c.2027delG (p.Gly676Alafs * 33) in both the homozygous or heterozygous states, but not in unrelated patients. Based on the variant zygosity, a marked phenotypic heterogeneity in terms of LDL-C levels, clinical presentations and resistance to anti-lipid treatment regimen (ACE inhibitors, β-blockers, ezetimibe, statins) of the FH patients was observed. This loss-of-function mutation is predicted to alter the free energy dynamics of the transcribed RNA, leading to its instability. Protein structural mapping has predicted that this non-sense mutation eliminates key functional domains in LDLR, which are essential for the receptor recycling and LDL particle binding. In conclusion, by combining genetics and structural bioinformatics approaches, this study identified and characterized a very rare FH causative LDLR pathogenic variant determining both clinical presentation and resistance to anti-lipid drug treatment., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Awan, Rashidi, Al-Shehri, Jamil, Elango, Al-Aama, Hegele, Banaganapalli and Shaik.)

Background: Celiac disease (CD) is a genetically complex autoimmune disease which is triggered by dietary gluten. Human leukocyte antigen (HLA) class II genes are known to act as high-risk markers for CD, where >95% of CD patients carry (HLA), DQ2 and/or DQ8 alleles. Therefore, the present study was conducted to investigate the distribution of HLA haplotypes among Saudi CD patients and healthy controls by using the tag single nucleotide polymorphisms (SNP)., Methods: HLA-tag SNPs showing strong linkage value (r2>0.99) were used to predict the HLA DQ2 and DQ8 genotypes in 101 Saudi CD patients and in 103 healthy controls by using real-time polymerase chain reaction technique. Genotype calls were further validated by Sanger sequencing method., Results: A total of 63.7% of CD cases and of 60.2% of controls were predicted to carry HLA-DQ2 and DQ8 heterodimers, either in the homozygous or heterozygous states. The prevalence of DQ8 in our CD patients was predicted to be higher than the patients from other ethnic populations (35.6%). More than 32% of the CD patients were found to be non-carriers of HLA risk haplotypes as predicted by the tag SNPs., Conclusion: The present study highlights that the Caucasian specific HLA-tag SNPs would be of limited value to accurately predict CD specific HLA haplotypes in Saudi population, when compared with the Caucasian groups. Prediction of risk haplotypes by tag SNPs in ethnic groups is a good alternate approach as long as the tag SNPs were identified from the local population genetic variant databases., (© 2021 The Author(s).)

Uterine smooth muscular neoplastic growths like benign leiomyomas (UL) and metastatic leiomyosarcomas (ULMS) share similar clinical symptoms, radiological and histological appearances making their clinical distinction a difficult task. Therefore, the objective of this study is to identify key genes and pathways involved in transformation of UL to ULMS through molecular differential analysis. Global gene expression profiles of 25 ULMS, 25 UL, and 29 myometrium (Myo) tissues generated on Affymetrix U133A 2.0 human genome microarrays were analyzed by deploying robust statistical, molecular interaction network, and pathway enrichment methods. The comparison of expression signals across Myo vs UL, Myo vs ULMS, and UL vs ULMS groups identified 249, 1037, and 716 significantly expressed genes, respectively (p ≤ 0.05). The analysis of 249 DEGs from Myo vs UL confirms multistage dysregulation of various key pathways in extracellular matrix, collagen, cell contact inhibition, and cytokine receptors transform normal myometrial cells to benign leiomyomas (p value ≤ 0.01). The 716 DEGs between UL vs ULMS were found to affect cell cycle, cell division related Rho GTPases and PI3K signaling pathways triggering uncontrolled growth and metastasis of tumor cells (p value ≤ 0.01). Integration of gene networking data, with additional parameters like estimation of mutation burden of tumors and cancer driver gene identification, has led to the finding of 4 hubs (JUN, VCAN, TOP2A, and COL1A1) and 8 bottleneck genes (PIK3R1, MYH11, KDR, ESR1, WT1, CCND1, EZH2, and CDKN2A), which showed a clear distinction in their distribution pattern among leiomyomas and leiomyosarcomas. This study provides vital clues for molecular distinction of UL and ULMS which could further assist in identification of specific diagnostic markers and therapeutic targets. Abbreviations UL: Uterine Leiomyomas; ULMS: Uterine Leiomyosarcoma; Myo: Myometrium; DEGs: Differential Expressed Genes; RMA: Robust Multiarray Average; DC: Degree of Centrality; BC: Betweenness of Centrality; CGC: Cancer Gene Census; FDR: False Discovery Rate; TCGA: Cancer Genome Atlas; BP: Biological Process; CC: Cellular Components; MF: Molecular Function; PPI: Protein-Protein Interaction.

The Paget disease (PDB; OMIM is 167250) is a chronic bone disease caused by pathogenic mutations in Sequestome1/p62 (SQSTM1) gene. This study has aimed to interpret the relationship of PDB linked SQSTM1 mutations with protein structure and its molecular dynamic features. The disease causative missense mutations were initially collected, and then analyzed for their, exonic and domain distribution, impact on secondary and tertiary structures, and their ability on protein-ligand interactions, using a combination of systems biology approaches. Our results show that most PDB linked SQSTM1 missense mutations affect amino acid residues clustered within or near the UBA domain (aa 389-434), which participates in the ubiquitination of substrates. We also report that the majority mutations occurred in α-helices over β-strands but their effects on the secondary structure were mostly neutral. Global tertiary structure deviations were minimal; however, at amino acid residue level minor structural changes were evident. The molecular dynamics simulation analysis showed that both PB1 and UBA domains were under constant structural fluctuations resulting in closed form conformation of SQSMT1 protein structure, when it is bound to PRKCI ligand. We also found salt bridge conformation changes in the UBA domain of SQSTM1 mutants when they bound to the PRKCI interactor protein. This finding suggests the possibility that mutations in SQSTM1 could impair its ability to ubiquitinate the substrates, eventually affecting autophagy and apoptosis, especially in mature osteoclasts. This study presents the additional insight into structure and function relationship between SQSTM1 mutations and PDB pathogenesis. Communicated by Ramaswamy H. Sarma.

Background: Alström syndrome (AS) is a very rare childhood disorder characterized by cardiomyopathy, progressive hearing loss and blindness. Inherited genetic variants of ALMS1 gene are the known molecular cause of this disease. The objective of this study was to characterize the genetic basis and understand the genotype-phenotype relationship in Saudi AS patients. Methods: Clinical phenotyping and whole-exome sequencing (WES) analysis were performed on six AS patients belonging to two unrelated consanguineous Saudi families. Sanger sequencing was performed to determine the mode of inheritance of ALMS1 variant in first-degree family relatives and also to ensure its rare prevalence in 100 healthy population controls. Results: We identified that Alström patients from both the families were sharing a very rare ALMS1, 3'-splice site acceptor (c.11873-2 A>T) variant, which skips entire exon-19 and shortens the protein by 80 amino acids. This disease variant was inherited by AS patients in autosomal recessive mode and is not yet reported in any population-specific genetic databases. AS patients carrying this mutation showed heterogeneity in clinical presentations. Computational analysis of the mutant centroid structure of ALMS1 mRNA revealed that exon-19 skipping enlarges the hairpin loop and decreases the free energy, eventually affecting its folding pattern, stability, and function. Hence, we propose c.11873-2A as an AS causative potential founder mutation in Saudi Arabia because it is found in two families lacking a common lineage. Conclusions: We conclude that WES analysis potentially helps in clinical phenotyping, early diagnosis, and better clinical management of Alström patients showing variable clinical expressivity., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Saadah, Banaganapalli, Kamal, Sahly, Alsufyani, Mohammed, Ahmad, Nasser, Al-Aama, Shaik and Elango.)

Myocardial infarction (MI) is the most prevalent coronary heart disease caused by the complex molecular interactions between multiple genes and environment. Here, we aim to identify potential biomarkers for the disease development and for prognosis of MI. We have used gene expression dataset (GSE66360) generated from 51 healthy controls and 49 patients experiencing acute MI and analyzed the differentially expressed genes (DEGs), protein-protein interactions (PPI), gene network-clusters to annotate the candidate pathways relevant to MI pathogenesis. Bioinformatic analysis revealed 810 DEGs. Their functional annotations have captured several MI targeting biological processes and pathways like immune response, inflammation and platelets degranulation. PPI network identify seventeen hub and bottleneck genes, whose involvement in MI was further confirmed by DisGeNET database. OpenTarget Platform reveal unique bottleneck genes as potential target for MI. Our findings identify several potential biomarkers associated with early stage MI providing a new insight into molecular mechanism underlying the disease., (Copyright © 2020 Elsevier Inc. All rights reserved.)

Celiac disease (CeD) is a gastrointestinal autoimmune disorder, whose specific molecular basis is not yet fully interpreted. Therefore, in this study, we compared the global gene expression profile of duodenum tissues from CeD patients, both at the time of disease diagnosis and after two years of the gluten-free diet. A series of advanced systems biology approaches like differential gene expression, protein-protein interactions, gene network-cluster analysis were deployed to annotate the candidate pathways relevant to CeD pathogenesis. The duodenum tissues from CeD patients revealed the differential expression of 106 up- and 193 down-regulated genes. The pathway enrichment of differentially expressed genes (DEGs) highlights the involvement of biological pathways related to loss of cell division regulation (cell cycle, p53 signalling pathway), immune system processes (NOD-like receptor signalling pathway, Th1, and Th2 cell differentiation, IL-17 signalling pathway) and impaired metabolism and absorption (mineral and vitamin absorptions and drug metabolism) in celiac disease. The molecular dysfunctions of these 3 biological events tend to increase the number of intraepithelial lymphocytes (IELs) and villous atrophy of the duodenal mucosa promoting the development of CeD. For the first time, this study highlights the involvement of aberrant cell division, immune system, absorption, and metabolism pathways in CeD pathophysiology and presents potential novel therapeutic opportunities.

Carbonic anhydrase 2 (CA2) enzyme deficiency caused by CA2 gene mutations is an inherited disorder characterized by symptoms like osteopetrosis, renal tubular acidosis, and cerebral calcification. This study has collected the CA2 deficiency causal missense mutations and assessed their pathogenicity using diverse computational programs. The 3D protein models for all missense mutations were built, and analyzed for structural divergence, protein stability, and molecular dynamics properties. We found M-CAP as the most sensitive prediction method to measure the deleterious potential of CA2 missense mutations. Free energy dynamics of tertiary structure models of CA2 mutants with DUET, mCSM, and SDM based consensus methods predicted only 50% of the variants as destabilizing. Superimposition of native and mutant CA2 models revealed the minor structural fluctuations at the amino acid residue level but not at the whole protein structure level. Near native molecular dynamic simulation analysis indicated that CA2 causative missense variants result in residue level fluctuation pattern in the protein structure. This study expands the understanding of genotype-protein phenotype correlations underlying CA2 variant pathogenicity and presents a potential avenue for modifying the CA2 deficiency by targeting biophysical structural features of CA2 protein. Communicated by Ramaswamy H. Sarma.

Objective: The objective of this study was to investigate the association of rs1051740, rs2234922 (in microsomal epoxide hydrolase 1; EPHX1 ), rs268 (in lipoprotein lipase; LPL ) and rs6025 (in Factor V Leiden; F5 ) genetic variants with the risk of preeclampsia development in Saudi women., Materials and Methods: This case-control study recruited 233 Saudi women (94 preeclampsia cases and 139 healthy controls) who visited the Gynecology and Obstetrics Departments of two hospitals in Jeddah, Saudi Arabia, for routine postpregnancy clinical follow-ups. All the women underwent thorough clinical and biochemical investigations conducted according to the standard clinical guidelines. Genotyping of the study participants was done using real-time polymerase chain reaction-based TaqMan allelic discrimination assay. The strength of the association between genetic variants and disease development was assessed using chi-square, odds ratio, 95% confidence interval and multifactor dimensionality reduction tests., Result: The minor alleles "G" in rs268 ( LPL ) and "A" in rs6025 ( F5 ) were absent in Saudi women. The frequencies of rs1051740 and rs2234922 of EPHX1 , both in the homozygous and allelic forms, were not significantly different between preeclampsia patients and healthy controls (for all tests, P > 0.05). The multifactor dimensionality reduction analysis also indicated that the interaction between the four studied single-nucleotide polymorphisms (SNPs) had no significant association with preeclampsia risk., Conclusion: This study found that none of the studied genetic variants (neither the single SNP nor the SNP-SNP interactions) explain the development of preeclampsia in the Saudi population. These findings not only underscore the disease heterogeneity but also highlight the need to develop population-specific diagnostic genetic biomarkers for preeclampsia., Competing Interests: There are no conflicts of interest., (Copyright: © 2020 Saudi Journal of Medicine & Medical Sciences.)

Background: Familial hypercholesterolemia (FH) is a lipid disorder caused by pathogenic mutations in LDLRAP1 gene. The present study has aimed to deepen our understanding about the pathogenicity predictions of FH causative genetic mutations, as well as their relationship to phenotype changes in LDLRAP1 protein, by utilizing multidirectional computational analysis., Methods: FH linked LDLRAP1 mutations were mined from databases, and the prediction ability of several pathogenicity classifiers against these clinical variants, was assessed through different statistical measures. Furthermore, these mutations were 3D modelled in protein structures to assess their impact on protein phenotype changes., Results: Our findings suggest that Polyphen-2, when compared with SIFT, M-CAP and CADD tools, can make better pathogenicity predictions for FH causative LDLRAP1 mutations. Through, 3D simulation and superimposition analysis of LDLRAP1 protein structures, it was found that missense mutations do not create any gross changes in the protein structure, although they could induce subtle structural changes at the level of amino acid residues. Near native molecular dynamic analysis revealed that missense mutations could induce variable degrees of fluctuation differences guiding the protein flexibility. Stability analysis showed that most missense mutations shifts the free energy equilibrium and hence they destabilize the protein. Molecular docking analysis demonstrates the molecular shifts in hydrogen and ionic bonds and Van der waals bonding properties, which further cause differences in the binding energy of LDLR-LDLRAP1 proteins., Conclusions: The diverse computational approaches used in the present study may provide a new dimension for exploring the structure-function relationship of the novel and deleterious LDLRAP1 mutations linked to FH., (© 2020 John Wiley & Sons, Ltd.)

Celiac disease (CD) is a gastrointestinal disorder whose genetic basis is not fully understood. Therefore, we studied a Saudi family with two CD affected siblings to discover the causal genetic defect. Through whole exome sequencing (WES), we identified that both siblings have inherited an extremely rare and deleterious CPED1 genetic variant (c.241 A > G; p.Thr81Ala) segregating as autosomal recessive mutation, suggesting its putative causal role in the CD. Saudi population specific minor allele frequency (MAF) analysis has confirmed its extremely rare prevalence in homozygous condition (MAF is 0.0004). The Sanger sequencing analysis confirmed the absence of this homozygous variant in 100 sporadic Saudi CD cases. Genotype-Tissue Expression (GTEx) data has revealed that CPED1 is abundantly expressed in gastrointestinal mucosa. By using a combination of systems biology approaches like protein 3D modeling, stability analysis and nucleotide sequence conservation analysis, we have further established that this variant is deleterious to the structural and functional aspects of CPED1 protein. To the best of our knowledge, this variant has not been previously reported in CD or any other gastrointestinal disease. The cell culture and animal model studies could provide further insight into the exact role of CPED1 p.Thr81Ala variant in the pathophysiology of CD. In conclusion, by using WES and systems biology analysis, present study for the first-time reports CPED1 as a potential causative gene for CD in a Saudi family with potential implications to both disease diagnosis and genetic counseling., (© 2020 The Author(s).)

Obesity is a multifactorial disease caused by complex interactions between genes and dietary factors. Salt-rich diet is related to the development and progression of several chronic diseases including obesity. However, the molecular basis of how salt sensitivity genes (SSG) contribute to adiposity in obesity patients remains unexplored. In this study, we used the microarray expression data of visceral adipose tissue samples and constructed a complex protein-interaction network of salt sensitivity genes and their co-expressed genes to trace the molecular pathways connected to obesity. The Salt Sensitivity Protein Interaction Network (SSPIN) of 2691 differentially expressed genes and their 15474 interactions has shown that adipose tissues are enriched with the expression of 23 SSGs, 16 hubs and 84 bottlenecks (p = 2.52 x 10-16) involved in diverse molecular pathways connected to adiposity. Fifteen of these 23 SSGs along with 8 other SSGs showed a co-expression with enriched obesity-related genes (r ≥ 0.8). These SSGs and their co-expression partners are involved in diverse metabolic pathways including adipogenesis, adipocytokine signaling pathway, renin-angiotensin system, etc. This study concludes that SSGs could act as molecular signatures for tracing the basis of adipogenesis among obese patients. Integrated network centered methods may accelerate the identification of new molecular targets from the complex obesity genomics data., Competing Interests: The authors have declared that no competing interests exist.

Alström syndrome (AS, OMIM ID 203800) is a rare childhood multiorgan disorder, which is widely studied in non-Arab ethnic patients. The clinical and molecular basis of AS and the mode of disease inheritance in consanguineous Arab populations is not well investigated. Therefore, to identify the molecular basis of AS in familial forms, the present study performed whole exome sequencing of 5 AS patients belonging to 2 different Bedouin families from Saudi Arabia. The present study identified the AS causative rare biallelic mutations in ALMS gene:T376S in exon 5 and S909* in exon 8 for family A and an R2721* in exon 10 (R2721*) for family B. ALMS1 targeted genetic sequencing of healthy population controls and family members has confirmed its extremely rare frequency and autosomal recessive mode of inheritance. The truncating mutations S909* and R2721* could cause the loss of CC domains and ALMS motif on C-terminal end of the protein and creates unstable protein, which eventually undergoes intracellular degradation. The premature protein truncating mutations described in our study may eventually provide further insight into the functional domains of the ALMS1 protein and contribute to the understanding of the phenotypic spectrum of AS. Whole exome sequencing based molecular diagnosis is expected to rule out ambiguity surrounding clinical diagnosis of suspected AS cases., Competing Interests: Authors declares that they have no conflict of interest., (© 2019 King Saud University.)

Mitochondrial disorders (MIDs) shows overlapping clinical presentations owing to the genetic and metabolic defects of mitochondria. However, specific relationship between inherited mutations in nuclear encoded mitochondrial proteins and their functional impacts in terms of metabolic defects in patients is not yet well explored. Therefore, using high throughput whole exome sequencing (WES), we screened a chronic kidney disease (CKD) and sensorineural hearing loss (SNHL) patient, and her family members to ascertain the mode of inheritance of the mutation, and healthy population controls to establish its rare frequency. The impact of mutation on biophysical characteristics of the protein was further studied by mapping it in 3D structure. Furthermore, LC-MS tandem mass spectrophotometry based untargeted metabolomic profiling was done to study the fluctuations in plasma metabolites relevant to disease causative mutations and kidney damage. We identified a very rare homozygous c.631G > A (p.Val211Met) pathogenic mutation in RMND1 gene in the proband, which is inherited in an autosomal recessive fashion. This gene is involved in the mitochondrial translational pathways and contribute in mitochondrial energy metabolism. The p.Val211Met mutation is found to disturb the structural orientation (RMSD is -2.95 Å) and stability (ΔΔG is -0.552 Kcal/mol) of the RMND1 protein. Plasma metabolomics analysis revealed the aberrant accumulation of metabolites connected to lipid and amino acid metabolism pathways. Of these metabolites, pathway networking has discovered ceramide, a metabolite of sphingolipids, which plays a role in different signaling cascades including mitochondrial membrane biosynthesis, is highly elevated in this patient. This study suggests that genetic defects in RMND1 gene alters the mitochondrial energy metabolism leading to the accumulation of ceramide, and subsequently promote dysregulated apoptosis and tissue necrosis in kidneys., (© 2019 The Authors.)

Rheumatoid arthritis (RA) is a chronic synovial autoinflammatory disease that destructs the cartilage and bone, leading to disability. The functional regulation of major immunity-related pathways like nuclear factor kappa B (NF-κB), which is involved in the chronic inflammatory reactions underlying the development of RA, remains to be explored. Therefore, this study has adopted statistical and knowledge-based systemic investigations (like gene correlation, semantic similarity, and topological parameters based on graph theory) to study the gene expression status of NF-κB protein family (NK PF ) and its regulators in synovial tissues to trace the molecular pathways through which these regulators contribute to RA. A complex protein-protein interaction map (PPIM) of 2,742 genes and 37,032 interactions was constructed from differentially expressed genes ( p ≤ 0.05). PPIM was further decomposed into a Regulator Allied Protein Interaction Network (RA PIN ) based on the interaction between genes (5 NK PF , 31 seeds, 131 hubs, and 652 bottlenecks). Pathway network analysis has shown the RA-specific disturbances in the functional connectivity between seed genes ( RIPK1 , ATG7 , TLR4 , TNFRSF1A , KPNA1 , CFLAR , SNW1 , FOSB , PARVA , CX3CL1 , and TRPC6 ) and NK PF members ( RELA , RELB , NFKB2 , and REL ). Interestingly, these genes are known for their involvement in inflammation and immune system (signaling by interleukins, cytokine signaling in immune system, NOD-like receptor signaling, MAPK signaling, Toll-like receptor signaling, and TNF signaling) pathways connected to RA. This study, for the first time, reports that SNW1, along with other NK regulatory genes, plays an important role in RA pathogenesis and might act as potential biomarker for RA. Additionally, these genes might play important roles in RA pathogenesis, as well as facilitate the development of effective targeted therapies. Our integrative data analysis and network-based methods could accelerate the identification of novel drug targets for RA from high-throughput genomic data., (Copyright © 2019 Sabir, El Omri, Banaganapalli, Al-Shaeri, Alkenani, Sabir, Hajrah, Zrelli, Ciesla, Nasser, Elango, Shaik and Khan.)

Obesity is a multifactorial metabolic disorder characterized by low grade chronic inflammation. Rare and novel mutations in genes which are vital in several key pathways have been reported to alter the energy expenditure which regulates body weight. The TP53 or p53 gene plays a prominent role in regulating various metabolic activities such as glycolysis, lipolysis, and glycogen synthesis. Recent genome-wide association studies reported that tumor suppressor gene p53 variants play a critical role in the predisposition of type 2 diabetes and obesity. Till date, no reports are available from the Arabian population; hence the present study was intended to assess the association between p53 variants with risk of obesity development in the Saudi population. We have selected three p53 polymorphisms, rs1642785 (C > G), and rs9894946 (A > G), and rs1042522 (Pro72Arg; C > G) and assessed their association with obesity risk in the Saudi population. Phenotypic and biochemical parameters were also evaluated to check their association with p53 genotypes and obesity. Genotyping was carried out on 136 obese and 122 normal samples. We observed that there is significantly increased prevalence p52 Pro72Arg (rs1042522) polymorphism in obese persons when compared to controls at GG genotype in overall comparison (OR: 2.169, 95% CI: 1.086-4.334, p = 0.02716). Male obese subjects showed three-fold higher risk at GG genotype (OR: 3.275, 95% CI: 1.230-8.716, p = 0.01560) and two-fold risk at G allele (OR: 1.827, 95% CI: 1.128-2.958, p = 0.01388) of p53 variant Pro72Arg respectively. This variant has also shown significant influence on cholesterol, LDL level, and random insulin levels in obese subjects (p ≤ 0.05). In conclusion, p53 Pro72Arg variant is highly prevalent among obese individuals and may act as a genetic modifier for obesity development among Saudis., (© 2019 Production and hosting by Elsevier B.V. on behalf of King Saud University.)

Obesity is connected to the activation of chronic inflammatory pathways in both adipocytes and macrophages located in adipose tissues. The nuclear factor (NF)-κB is a central molecule involved in inflammatory pathways linked to the pathology of different complex metabolic disorders. Investigating the gene expression data in the adipose tissue would potentially unravel disease relevant gene interactions. The present study is aimed at creating a signature molecular network and at prioritizing the potential biomarkers interacting with NF-κB family of proteins in obesity using system biology approaches. The dataset GSE88837 associated with obesity was downloaded from Gene Expression Omnibus (GEO) database. Statistical analysis represented the differential expression of a total of 2650 genes in adipose tissues (p = <0.05). Using concepts like correlation, semantic similarity, and theoretical graph parameters we narrowed down genes to a network of 23 genes strongly connected with NF-κB family with higher significance. Functional enrichment analysis revealed 21 of 23 target genes of NF-κB were found to have a critical role in the pathophysiology of obesity. Interestingly, GEM and PPP1R13L were predicted as novel genes which may act as potential target or biomarkers of obesity as they occur with other 21 target genes with known obesity relationship. Our study concludes that NF-κB and prioritized target genes regulate the inflammation in adipose tissues through several molecular signaling pathways like NF-κB, PI3K-Akt, glucocorticoid receptor regulatory network, angiogenesis and cytokine pathways. This integrated system biology approaches can be applied for elucidating functional protein interaction networks of NF-κB protein family in different complex diseases. Our integrative and network-based approach for finding therapeutic targets in genomic data could accelerate the identification of novel drug targets for obesity., Competing Interests: The authors have declared that no competing interests exist.

Rheumatoid arthritis (RA) is a chronic autoimmune disorder of bone joints caused by the complex interplay between several factors like body physiology, the environment with genetic background. The recent meta-analysis of GWAS has expanded the total number of RA-associated loci to more than 100, out of which approximately ∼97% (98 variants) loci are located in non-coding regions, and the other ∼3% (3 variants) are in three different non-HLA genes, i.e., TYK2 (Prp1104Ala), IL6R (Asp358Ala), and PTPN22 (Trp620Arg). However, whether these variants prompt changes in the protein phenotype with regards to its stability, structure, and interaction with other molecules, remains unknown. Thus, we selected the three clinically pathogenic variants described above, as positive controls and applied diverse computational methods to scrutinize if those mutations cause changes in the protein phenotype. Both wild type and mutant protein structures of PTPN22 (W620R), IL6R (D358A), and TYK2 (P1104A) were modeled and studied for structural deviations. Furthermore, we have also studied the secondary structure characteristics, solvent accessibility and stability, and the molecular interaction deformities caused by the amino acid substitutions. We observed that simple nucleotide predictions of SIFT, PolyPhen, CADD and FATHMM yields mixed findings in screening the RA-missense variants which showed a ≥ P -value threshold of 5 × 10 -8 in genome wide association studies. However, structure-based analysis confirms that mutant structures shows subtle but significant changes at their core regions, but their functional domains seems to lose wild type like functional interaction. Our findings suggest that the multidirectional computational analysis of clinically potential RA-mutations could act as a primary screening step before undertaking functional biology assays.

The somatic mutations in ATP binding cleft of the tyrosine kinase binding domain of EGFR are known to occur in 15-40% of non-small cell lung cancer (NSCLC) patients. Although first and second generation anti-EGFR inhibitors are widely used to treat these patients, their therapeutic efficacy is modest and often results in adverse effects or drug resistance. Therefore, there is a need to develop novel as well as safe anti-EGFR drugs. The rapid emergence of computational drug designing provided a great opportunity to both discover and predict the efficacy of novel EGFR inhibitors from plant sources. In the present study, we designed several chemical analogues of edible curcumin (CUCM) compound and assessed their drug likeliness, ADME and toxicity properties using a diverse range of advanced computational methods. We also have examined the structural plasticity and binding characteristics of EGFR wild-type and mutant forms (S769L and K846R) against ligand molecules like Gefitinib, native CUCM, and different CUCM analogues. Through multidimensional experimental approaches, we conclude that CUCM-36 ((1E,4Z,6E)-1-(3,4-Diphenoxyphenyl)-5-hydroxy-7-(4-hydroxy-3-phenoxyphenyl)-1,4,6-heptatrien-3-one) is the best anti-EGFR compound with high drug-likeness, ADME properties, and low toxicity properties. CUCM-36 compound has demonstrated better affinity towards both wild-type (ΔG is -8.5 kcal/Mol) and mutant forms (V769L & K846R; ΔG for both is >-9.20 kcal/Mol) compared to natural CUCM and Gefitinib inhibitor. This study advises the future laboratory assays to develop CUCM-36 as a novel drug compound for treating EGFR positive non-small cell lung cancer patients.

MED12, a subunit of mediator complex genes is known to harbor genetic mutations, (mostly in exon 2), causal to the genesis of uterine leiomyomas among Caucasian, African American, and Asian women. However, the precise relationship between genetic mutations vs. protein or disease phenotype is not well-explained. Therefore, we sought to replicate the MED12 mutation frequency in leiomyomas of Saudi Arabian women, who represents ethnically and culturally distinct population. We performed molecular screening of MED12 gene (in 308 chromosomes belonging to 154 uterine biopsies), analyzed the genotype-disease phenotype correlations and determined the biophysical characteristics of mutated protein through diverse computational approaches. We discovered that >44% (34/77) leiomyomas of Arab women carry a spectrum of MED12 mutations (30 missense, 1 splice site, and 3 indels). In addition to known codon 44, we observed novel somatic mutations in codons 36, 38, and 55. Most genetically mutated tumors (27/30; 90%) demonstrated only one type of genetic change, highlighting that even single allele change in MED12 can have profound impact in transforming the normal uterine myometrium to leiomyomas. An interesting inverse correlation between tumor size and LH is observed when tumor is positive to MED12 mutation ( p < 0.05). Our computational investigations suggest that amino acid substitution mutations in exon-2 region of MED12 might contribute to potential alterations in phenotype as well as the stability of MED12 protein. Our study, being the first one from Arab world, confirms the previous findings that somatic MED12 mutations are critical to development and progression of uterine leiomyomas irrespective of the ethnic background. We recommend that mutation screening, particularly codon 44 of MED12 can assist in molecular diagnostics of uterine leiomyomas in majority of the patients.

The life-threatening group of irregular cardiac rhythmic disorders also known as Cardiac Arrhythmias (CA) are caused by mutations in highly conserved Calmodulin (CALM/CaM) genes. Herein, we present a multidimensional approach to diagnose changes in phenotypic, stability, and Ca 2+ ion binding properties of CA-causing mutations. Mutation pathogenicity was determined by diverse computational machine learning approaches. We further modeled the mutations in 3D protein structure and analyzed residue level phenotype plasticity. We have also examined the influence of torsion angles, number of H-bonds, and free energy dynamics on the stability, near-native simulation dynamic potential of residue fluctuations in protein structures, Ca 2+ ion binding potentials, of CaM mutants. Our study recomends to use M-CAP method for measuring the pathogenicity of CA causing CaM variants. Interestingly, most CA-causing variants we analyzed, exists in either third (V/H-96, S/I-98, V-103) or fourth (G/V-130, V/E/H-132, H-134, P-136, G-141, and L-142) EF-hands located in carboxyl domains of the CaM molecule. We observed that the minor structural fluctuations caused by these variants are likely tolerable owing to the highly flexible nature of calmodulin's globular domains. However, our molecular docking results supports that these variants disturb the affinity of CaM toward Ca 2+ ions and corroborate previous findings from functional studies. Taken together, these computational findings can explain the molecular reasons for subtle changes in structure, flexibility, and stability aspects of mutant CaM molecule. Our comprehensive molecular scanning approach demonstrates the utility of computational methods in quick preliminary screening of CA- CaM mutations before undertaking time consuming and complicated functional laboratory assays., (© 2018 Wiley Periodicals, Inc.)

The deleterious amino acid substitution mutations in IL-10 receptor alpha gene are most frequently reported in several autoimmune diseases including early onset-inflammatory bowel disease (IBD). Despite the important role of IL-10 RA in maintaining immune homeostasis, the specific structural and functional implications of these mutations on protein phenotype, stability, ligand binding and post translational characteristics is not well explored. Therefore, this study performed the multidimensional computational analysis of IL10RA missense variations causative to pediatric or early onset inflammatory bowel disease (<5 years of age). Our computational algorithmic screening identified the deleterious nature of p. W45G, p. Y57C, p. W69G, p.T84I, p.Y91C, p.R101W, p.R117C, and p.R117H, IBD causative IL10-RA mutations. The sensitivity and specificity analysis of different computational methods showed that CADD outperform SIFT, PolyPhen 2.0, FATHMM, LRT, MetaLR, MetaSVM, PROVEAN and Condel in predicting the pathogenicity of IL10RA mutations. Our three-dimensional protein modeling assays showed that the point mutations cause major drifts in the structural plasticity of IL10 RA molecule and negatively influence its stability. Findings from molecular docking analysis have shown that these point mutations decrease the binding affinity of IL10RA toward IL10 and may likely to disturb the IL10 signaling pathway. This study provides an easy frame work for phenotypic characterization of mutant IL10RA molecule in terms of structure, flexibility and stability aspects. Our approach may also add a new dimension to conventional functional biology assays in quickly studying IL10 RA mutations and also for designing and developing inhibitors for mutant IL10RA molecule.