A review of novel coronavirus disease (COVID-19): based on genomic structure, phylogeny, current shreds of evidence, candidate vaccines, and drug repurposing

Coronavirus disease (COVID-19) pandemic is instigated by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). As of March 13, 2021, more than 118.9 million cases were infected with COVID-19 worldwide. SARS-CoV-2 is a positive-sense single-stranded RNA beta-CoV. Most COVID-19 infected individuals recover within 1-3 weeks. Nevertheless, approximately 5% of patients develop acute respiratory distress syndrome and other systemic complications, leading to death. Structural genetic analyses of SARS-CoV-2 have shown genomic resemblances but a low evolutionary correlation to SARS-CoV-1 responsible for the 2002-2004 outbreak. The S glycoprotein is critical for cell adhesion and the entrance of the virus into the host. The process of cell entry uses the cellular receptor named angiotensin-converting enzyme 2. Recent evidence proposed that the CD147 as a SARS-CoV-2's potential receptor. The viral genome is mainly held by two non-structural proteins (NSPs), ORF1a and ORF1ab, along with structural proteins. Although NSPs are conserved among the βCoVs, mutations in NSP2 and NSP3 may play critical roles in transmitting the virus and cell tropism. To date, no specific/targeted anti-viral treatments exist. Notably, more than 50 COVID-19 candidate vaccines in clinical trials, and a few being administered. Preventive precautions are the primary strategy to limit the viral load transmission and spread, emphasizing the urgent need for developing significant drug targets and vaccines against COVID-19. This review provides a cumulative overview of the genomic structure, transmission, phylogeny of SARS-CoV-2 from Indian clusters, treatment options, updated discoveries, and future standpoints for COVID-19.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02749-0.

A systemic approach to explore the mechanisms of drug resistance and altered signaling cascades in extensively drug-resistant tuberculosis

BACKGROUND AND AIM

The persistence of extensively drug-resistant (XDR) strains of Mycobacterium tuberculosis (MTB) continue to pose a significant challenge to the treatment and control of tuberculosis infections worldwide. XDR-MTB strains exhibit resistance against first-line anti-TB drugs, fluoroquinolones, and second-line injectable drugs. The mechanisms of drug resistance of MTB remains poorly understood. Our study aims at identifying the differentially expressed genes (DEGs), associated gene networks, and signaling cascades involved in rendering this pathogen resistant to multiple drugs, namely, isoniazid, rifampicin, and capreomycin.

METHODS

We used the microarray dataset GSE53843. The GEO2R tool was used to prioritize the most significant DEGs (top 250) of each drug exposure sample between XDR strains and non-resistant strains. The validation of the 250 DEGs was performed using volcano plots. Protein-protein interaction networks of the DEGs were created using STRING and Cytoscape tools, which helped decipher the relationship between these genes. The significant DEGs were functionally annotated using DAVID and ClueGO. The concomitant biological processes (BP) and molecular functions (MF) were represented as dot plots.

RESULTS AND CONCLUSION

We identified relevant molecular pathways and biological processes, such as cell wall biogenesis, lipid metabolic process, ion transport, phosphopantetheine binding, and triglyceride lipase activity. These processes indicated the involvement of multiple interconnected mechanisms in drug resistance. Our study highlighted the impact of cell wall permeability, with the dysregulation of the mur family of proteins, as essential factors in the inference of resistance. Additionally, upregulation of genes responsible for ion transport such as ctpF, arsC, and nark3, emphasizes the importance of transport channels and efflux pumps in potentially driving out stress-inducing compounds. This study investigated the upregulation of the Lip family of proteins, which play a crucial role in triglyceride lipase activity. Thereby illuminating the potential role of drug-induced dormancy and subsequent resistance in the mycobacterial strains. Multiple mechanisms such as carboxylic acid metabolic process, NAD biosynthetic process, triglyceride lipase activity, phosphopantetheine binding, organic acid biosynthetic process, and growth of symbiont in host cell were observed to partake in resistance of XDR-MTB. This study ultimately provides a platform for important mapping targets for potential therapeutics against XDR-MTB.

Integrated approach on azo dyes degradation using laccase enzyme and Cul nanoparticle

Azo dyes released by the textile industries cause severe damage to the environment and living organisms. The degradation of azo dyes is widely studied using enzymatic methods. Laccase is a copper-containing enzyme that degrades the azo dyes into less toxic compounds. In this work, the crude laccase enzyme produced by the alkaliphile Pseudomonas mendocina in the degradation of mixed azo dye showed 0.386 U/mL activity at pH 8.5. A combination of enzymatic and green synthesized nanoparticles was used in the degradation of mixed azo dye. Laccase used in the degradation of mixed azo dyes showed 58.4% in 72 h, while the photocatalytic degradation of mixed azo dyes showed 15.9%. The degradation of azo dyes using copper iodide nanoparticles resulted in 15.8% degradation. However, it was noticed that the combined method of degradation of azo dyes involving both crude laccase and CuI nanoparticles gave a degradation of 62.3% in 60 min. Interaction of laccase enzyme with azo dyes using in silico analysis predicted the binding energy with reactive red (−7.19 kcal/mol), reactive brown (−8.57 kcal/mol), and reactive black dyes (−9.17 kcal/mol) respectively.

Molecular dynamics, residue network analysis, and cross-correlation matrix to characterize the deleterious missense mutations in GALE causing galactosemia III

Epimerase-deficiency galactosemia (EDG) is caused by mutations in the UDP-galactose 4'-epimerase enzyme, encoded by gene GALE. Catalyzing the last reaction in the Leloir pathway, UDP-galactose-4-epimerase catalyzes the interconversion of UDP-galactose and UDP-glucose. This study aimed to use in-depth computational strategies to prioritize the pathogenic missense mutations in GALE protein and investigate the systemic behavior, conformational spaces, atomic motions, and cross-correlation matrix of the GALE protein. We searched four databases (dbSNP, ClinVar, UniProt, and HGMD) and major biological literature databases (PubMed, Science Direct, and Google Scholar), for missense mutations that are associated with EDG patients, our search yielded 190 missense mutations. We applied a systematic computational prediction pipeline, including pathogenicity, stability, biochemical, conservational, protein residue contacts, and structural analysis, to predict the pathogenicity of these mutations. We found three mutations (p.K161N, p.R239W, and p.G302D) with a severe phenotype in patients with EDG that correlated with our computational prediction analysis; thus, they were selected for further structural and simulation analyses to compute the flexibility and stability of the mutant GALE proteins. The three mutants were subjected to molecular dynamics simulation (MDS) with native protein for 200 ns using GROMACS. The MDS demonstrated that these mutations affected the beta-sheets and helical region that are responsible for the catalytic activity; subsequently, affects the stability and flexibility of the mutant proteins along with a decrease and more deviations in compactness when compared to that of a native. Also, three mutations created major variations in the combined atomic motions of the catalytic and C-terminal regions. The network analysis of the residues in the native and three mutant protein structures showed disturbed residue contacts occurred owing to the missense mutations. Our findings help to understand the structural behavior of a protein owing to mutation and are intended to serve as a platform for prioritizing mutations, which could be potential targets for drug discovery and development of targeted therapeutics.

Network analysis of transcriptomics data for the prediction and prioritization of membrane-associated biomarkers for idiopathic pulmonary fibrosis (IPF) by bioinformatics approach

Idiopathic pulmonary fibrosis (IPF) is a rare yet crucial persistent lung disorder that actuates scarring of lung tissues, which makes breathing difficult. Smoking, environmental pollution, and certain viral infections could initiate lung scarring. However, the molecular mechanism involved in IPF remains elusive. To develop an efficient therapeutic arsenal against IPF, it is vital to understand the pathology and deviations in biochemical pathways that lead to disorder. In this study, we availed network analysis and other computational pipelines to delineate the prominent membrane proteins as diagnostic biomarkers and therapeutic targets for IPF. This study yielded a significant role of glycosaminoglycan binding, endothelin, and GABA-B receptor signaling pathway in IPF pathogenesis. Furthermore, ADCY8, CRH, FGB, GPR17, MCHR1, NMUR1, and SAA1 genes were found to be immensely involved with IPF, and the enrichment pathway analysis suggests that most of the pathways were corresponding to membrane transport and signal transduction functionalities. This analysis could help in better understanding the molecular mechanism behind IPF to develop an efficient therapeutic target or biomarkers for IPF.

Investigating mutations at the hotspot position of the ERBB2 and screening for the novel lead compound to treat breast cancer - a computational approach

Membrane proteins are the most common types of cancer that are active in the prognosis. Membrane proteins are a distinguishing characteristic of a cancer cell. In tumor cell therapy, the overexpressed membrane proteins are becoming ever more relevant. The 3-kinase (PI3K)/AKT phosphatidylinositol pathway is downstream triggered by different extracellular signals, and this signaling pathway activation impacts a variety of proliferation of the cellular processes like cell growth and surviving. Frequent PI3K/AKT dysregulation in human cancer has rendered proteins of this pathway desirable for diagnostic markers. Members of the ERBB family-like ERBB2 and ERBB3 activate intracellular signaling pathways such as PI3K/AKT. The mutations in these proteins dysfunctions the proteins in the downstream. Considering this importance, we have developed a computational pipeline to identify the mutation position with a highest number of mutations and to screen them for pathogenicity, stability, conservation, and structural changes using PredictSNP, iStable, ConSurf, and GROMACS simulation software respectively. Further, a virtual screening approach was initiated to find the most similar non-toxic lead compound, which could be an alternative to the currently used lapatinib. To conclude, protein-ligand dynamics were undertaken to study the actions of native and mutants with the lapatinib and the lead compound. From the overall analysis, we identified position 755 with leucine in the native condition is prone to frequent mutations. The leucine at 755th position is more prone to mutate as serine and tryptophan. Further from the computational analysis, we identified that the mutation L755S is more significant than the L755W mutation. We have witnessed CID140590176 be a potential lead compound with no toxicity. The behavior of the lead compound has shown more compactness with an increased number of intermolecular hydrogen bonds in the ERBB2 with L755S. This lead compound can be further taken for experimental validations, and we believe that this lead compound could be a potent ERBB2 inhibitor.

Mutational landscape of K-Ras substitutions at 12th position-a systematic molecular dynamics approach

K-Ras is a small GTPase and acts as a molecular switch by recruiting GEFs and GAPs, and alternates between the inert GDP-bound and the dynamic GTP-bound forms. The amino acid at position 12 of K-Ras is a hot spot for oncogenic mutations (G12A, G12C, G12D, G12R, G12S, and G12V), disturbing the active fold of the protein, leading to cancer development. This study aimed to investigate the potential conformational changes induced by these oncogenic mutations at the 12^th position, impairing GAP-mediated GTP hydrolysis. Comprehensive computational tools (iStable, FoldX, SNPeffect, DynaMut, and CUPSAT) were used to evaluate the effect of these six mutations on the stability of wild type K-Ras protein. The docking of GTP with K-Ras was carried out using AutoDock4.2, followed by molecular dynamics simulations. Furthermore, on comparison of binding energies between the wild type K-Ras and the six mutants, we have demonstrated that the G12A and G12V mutants exhibited the strongest binding efficiency compared to the other four mutants. Trajectory analyses of these mutations revealed that G12A encountered the least deviation, fluctuation, intermolecular H-bonds, and compactness compared to the wildtype, which was supported by the lower Gibbs free energy value. Our study investigates the molecular dynamics simulations of the mutant K-Ras forms at the 12^th position, which expects to provide insights about the molecular mechanisms involved in cancer development, and may serve as a platform for targeted therapies against cancer. Communicated by Ramaswamy H. Sarma.

Analysis of Differentially Expressed Genes and Molecular Pathways in Familial Hypercholesterolemia Involved in Atherosclerosis: A Systematic and Bioinformatics Approach

Background and Aims: Familial hypercholesterolemia (FH) is one of the major risk factor for the progression of atherosclerosis and coronary artery disease. This study focused on identifying the dysregulated molecular pathways and core genes that are differentially regulated in FH and to identify the possible genetic factors and potential underlying mechanisms that increase the risk to atherosclerosis in patients with FH. Methods: The Affymetrix microarray dataset (GSE13985) from the GEO database and the GEO2R statistical tool were used to identify the differentially expressed genes (DEGs) from the white blood cells (WBCs) of five heterozygous FH patients and five healthy controls. The interaction between the DEGs was identified by applying the STRING tool and visualized using Cytoscape software. MCODE was used to determine the gene cluster in the interactive networks. The identified DEGs were subjected to the DAVID v6.8 webserver and ClueGo/CluePedia for functional annotation, such as gene ontology (GO) and enriched molecular pathway analysis of DEGs. Results: We investigated the top 250 significant DEGs (p-value < 0.05; fold two change ≥ 1 or ≤ -1). The GO analysis of DEGs with significant differences revealed that they are involved in critical biological processes and molecular pathways, such as myeloid cell differentiation, peptidyl-lysine modification, signaling pathway of MyD88-dependent Toll-like receptor, and cell-cell adhesion. The analysis of enriched KEGG pathways revealed the association of the DEGs in ubiquitin-mediated proteolysis and cardiac muscle contraction. The genes involved in the molecular pathways were shown to be differentially regulated by either activating or inhibiting the genes that are essential for the canonical signaling pathways. Our study identified seven core genes (UQCR11, UBE2N, ADD1, TLN1, IRAK3, LY96, and MAP3K1) that are strongly linked to FH and lead to a higher risk of atherosclerosis. Conclusion: We identified seven core genes that represent potential molecular biomarkers for the diagnosis of atherosclerosis and might serve as a platform for developing therapeutics against both FH and atherosclerosis. However, functional studies are further needed to validate their role in the pathogenesis of FH and atherosclerosis.

Enzyme therapy: a forerunner in catalyzing a healthy society?

INTRODUCTION

The use of enzymes in various industries has been prevalent for centuries. However, their potency as therapeutics remained latent until the late 1950 s, when scientists finally realized the gold mine they were sitting on. Enzyme therapy has seen rapid development over the past few decades and has been widely used for the therapy of myriad diseases, including lysosomal storage disorders, cancer, Alzheimer's disease, irritable bowel syndrome, exocrine pancreatic insufficiency, and hyperuricemia. Enzymes are also used for wound healing, the treatment of microbial infections, and gene therapy.

AREAS COVERED

This is a comprehensive review of the therapeutic use of enzymes that can act as a guidepost for researchers and academicians and presents a general overview of the developments in enzyme therapy over the years, along with updates on recent advancements in enzyme therapy research.

EXPERT OPINION

Although enzyme therapy is immensely beneficial and induces little auxiliary damage, it has several drawbacks, ranging from high cost, low stability, low production, and hyperimmune responses to the failure to cure a variety of the problems associated with a disease. Further fine-tuning and additional clinical efficacy studies are required to establish enzyme therapy as a forerunner to catalyzing a healthy society.

The Rise and Impact of COVID-19 in India

The coronavirus disease (COVID-19) pandemic, which originated in the city of Wuhan, China, has quickly spread to various countries, with many cases having been reported worldwide. As of May 8th, 2020, in India, 56,342 positive cases have been reported. India, with a population of more than 1.34 billion-the second largest population in the world-will have difficulty in controlling the transmission of severe acute respiratory syndrome coronavirus 2 among its population. Multiple strategies would be highly necessary to handle the current outbreak; these include computational modeling, statistical tools, and quantitative analyses to control the spread as well as the rapid development of a new treatment. The Ministry of Health and Family Welfare of India has raised awareness about the recent outbreak and has taken necessary actions to control the spread of COVID-19. The central and state governments are taking several measures and formulating several wartime protocols to achieve this goal. Moreover, the Indian government implemented a 55-days lockdown throughout the country that started on March 25th, 2020, to reduce the transmission of the virus. This outbreak is inextricably linked to the economy of the nation, as it has dramatically impeded industrial sectors because people worldwide are currently cautious about engaging in business in the affected regions.

Dysregulation of Signaling Pathways Due to Differentially Expressed Genes From the B-Cell Transcriptomes of Systemic Lupus Erythematosus Patients - A Bioinformatics Approach

Systemic lupus erythematosus (SLE) is an autoimmune inflammatory disorder that is clinically complex and has increased production of autoantibodies. Via emerging technologies, researchers have identified genetic variants, expression profiling of genes, animal models, and epigenetic findings that have paved the way for a better understanding of the molecular and genetic mechanisms of SLE. Our current study aimed to illustrate the essential genes and molecular pathways that are potentially involved in the pathogenesis of SLE. This study incorporates the gene expression profiling data of the microarray dataset GSE30153 from the Gene Expression Omnibus (GEO) database, and differentially expressed genes (DEGs) between the B-cell transcriptomes of SLE patients and healthy controls were screened using the GEO2R web tool. The identified DEGs were subjected to STRING analysis and Cytoscape to explore the protein-protein interaction (PPI) networks between them. The MCODE (Molecular Complex Detection) plugin of Cytoscape was used to screen the cluster subnetworks that are highly interlinked between the DEGs. Subsequently, the clustered DEGs were subjected to functional annotation with ClueGO/CluePedia to identify the significant pathways that were enriched. For integrative analysis, we used GeneGo MetacoreTM, a Cortellis Solution software, to exhibit the Gene Ontology (GO) and enriched pathways between the datasets. Our study identified 4 upregulated and 13 downregulated genes. Analysis of GO and functional enrichment using ClueGO revealed the pathways that were statistically significant, including pathways involving T-cell costimulation, lymphocyte costimulation, negative regulation of vascular permeability, and B-cell receptor signaling. The DEGs were mainly enriched in metabolic networks such as the phosphatidylinositol-3,4,5-triphosphate pathway and the carnitine pathway. Additionally, potentially enriched pathways, such as the signaling pathways induced by oxidative stress and reactive oxygen species (ROS), chemotaxis and lysophosphatidic acid signaling induced via G protein-coupled receptors (GPCRs), and the androgen receptor activation pathway, were identified from the DEGs that were mainly associated with the immune system. Four genes (EGR1, CD38, CAV1, and AKT1) were identified to be strongly associated with SLE. Our integrative analysis using a multitude of bioinformatics tools might promote an understanding of the dysregulated pathways that are associated with SLE development and progression. The four DEGs in SLE patients might shed light on the pathogenesis of SLE and might serve as potential biomarkers in early diagnosis and as therapeutic targets for SLE.

Comprehensive in silico screening and molecular dynamics studies of missense mutations in Sjogren-Larsson syndrome associated with the ALDH3A2 gene

Sjögren-Larsson syndrome (SLS) is an autoimmune disorder inherited in an autosomal recessive pattern. To date, 80 missense mutations have been identified in association with the Aldehyde Dehydrogenase 3 Family Member A2 (ALDH3A2) gene causing SLS. Disruption of the function of ALDH3A2 leads to excessive accumulation of fat in the cells, which interferes with the normal function of protective membranes or materials that are necessary for the body to function normally. We retrieved 54 missense mutations in the ALDH3A2 from the OMIM, UniProt, dbSNP, and HGMD databases that are known to cause SLS. These mutations were examined with various in silico stability tools, which predicted that the mutations p.S308N and p.R423H that are located at the protein-protein interaction domains are the most destabilizing. Furthermore, to determine the atomistic-level differences within the protein-protein interactions owing to mutations, we performed macromolecular simulation (MMS) using GROMACS to validate the motion patterns and dynamic behavior of the biological system. We found that both mutations (p.S380N and p.R423H) had significant effects on the protein-protein interaction and disrupted the dimeric interactions. The computational pipeline provided in this study helps to elucidate the potential structural and functional differences between the ALDH3A2 native and mutant homodimeric proteins, and will pave the way for drug discovery against specific targets in the SLS patients.

Identification of potential inhibitors against pathogenic missense mutations of PMM2 using a structure-based virtual screening approach

The autosomal recessive phosphomannomutase 2-congenital disorder of glycosylation (PMM2-CDG) is characterized by defective functioning of the PMM2 enzyme, which is necessary for the conversion of mannose-6-phosphate into mannose-1-phosphate. Here, a computational pipeline was drawn to identify the most significant mutations, and further, we used a virtual screening approach to identify a new lead compound to treat the identified significant mutations. We searched for missense mutation data related to PMM2-CDG in HGMD®, UniProt, and ClinVar. Our search yielded a total of 103 mutations, of which 91 are missense mutations. The D65Y, I132N, I132T, and F183S mutations were classified as deleterious, destabilizing, and altering the biophysical properties using the PredictSNP, iStable, and Align GVGD in silico prediction tools. Additionally, we applied a multistep protocol to screen for an alternative lead compound to the existing CID2876053 (1-(3-chlorophenyl)-3,3-bis(pyridine-2-yl)urea) with affinity to these identified significant mutants. Two compounds, CHEMBL1491007 (6-chloro-4-phenyl-3-(4-pyridin-2-ylpiperazin-1-yl)-1H-quinolin-2-one) and CHEMBL3653029 (5-chloro-4-[6-[(3-fluorophenyl)methylamino]pyridin-2-yl]-N-(piperidin-4-ylmethyl)pyridin-2-amine), exhibited the highest binding affinity with the selected mutants and were chosen for further analysis. Through molecular docking, molecular dynamics simulation, and MMPBSA analysis, we found that the known compound, i.e. CID2876053, has stronger interaction with the D65Y mutant. The newly identified lead compound CHEMBL1491007 showed stronger interaction with the I132N and I132T mutants, whereas the most deleterious mutant, F183S, showed stronger interaction with CHEMBL3653029. This study is expected to aid in the field of precision medicine, and further to in vivo and in vitro analysis of these lead compounds might shed light on the treatment of PMM2-CDG. Communicated by Ramaswamy H. Sarma.

Integrative Bioinformatics Approaches to Map Potential Novel Genes and Pathways Involved in Ovarian Cancer

Background and aims: Ovarian cancer (OC) is the seventh most commonly detected cancer among women. This study aimed to map the hub and core genes and potential pathways that might be involved in the molecular pathogenesis of OC.

Methods: In the present work, we analyzed a microarray dataset (GSE126519) from the Gene Expression Omnibus (GEO) database and used the GEO2R tool to screen OC cells and ovarian SINE-resistant cancer cells for differentially expressed genes (DEGs). For the functional annotation of the DEGs, we conducted Gene Ontology (GO) and pathway enrichment analyses (KEGG) using the DAVID v6.8 online server and GenoGo Metacore™, Cortellis Solution software. Protein–protein interaction (PPI) networks were constructed using the STRING database, and Cytoscape software was used for visualization. The survival analysis was performed using the online platform GEPIA2 to determine the prognostic value of the expression of hub genes in cell lines from OC patients.

Results: We identified a total of 809 upregulated and 700 downregulated DEGs. GO analysis revealed that the genes with statistically significant differences in expression were mainly associated with biological processes involved in the cell cycle, the mitotic cell cycle, mitotic nuclear division, organ morphogenesis, cell development, and cell morphogenesis. By using the Analyze Networks (AN) algorithm in GeneGo, we identified the most relevant biological networks involving DEGs that were mainly enriched in the cell cycle (in metaphase checkpoints) and revealed the role of APC in cell cycle regulation pathways. We found 10 hub genes and four core genes (FZD6, FZD8, CDK2, and RBBP8) that are strongly linked to OC.

Conclusion: This study sheds light on the molecular pathogenesis of OC and is expected to provide potential molecular biomarkers that are beneficial for the treatment and clinical molecular diagnosis of OC.

Bioinformatics classification of mutations in patients with Mucopolysaccharidosis IIIA

Mucopolysaccharidosis (MPS) IIIA, also known as Sanfilippo syndrome type A, is a severe, progressive disease that affects the central nervous system (CNS). MPS IIIA is inherited in an autosomal recessive manner and is caused by a deficiency in the lysosomal enzyme sulfamidase, which is required for the degradation of heparan sulfate. The sulfamidase is produced by the N-sulphoglucosamine sulphohydrolase (SGSH) gene. In MPS IIIA patients, the excess of lysosomal storage of heparan sulfate often leads to mental retardation, hyperactive behavior, and connective tissue impairments, which occur due to various known missense mutations in the SGSH, leading to protein dysfunction. In this study, we focused on three mutations (R74C, S66W, and R245H) based on in silico pathogenic, conservation, and stability prediction tool studies. The three mutations were further subjected to molecular dynamic simulation (MDS) analysis using GROMACS simulation software to observe the structural changes they induced, and all the mutants exhibited maximum deviation patterns compared with the native protein. Conformational changes were observed in the mutants based on various geometrical parameters, such as conformational stability, fluctuation, and compactness, followed by hydrogen bonding, physicochemical properties, principal component analysis (PCA), and salt bridge analyses, which further validated the underlying cause of the protein instability. Additionally, secondary structure and surrounding amino acid analyses further confirmed the above results indicating the loss of protein function in the mutants compared with the native protein. The present results reveal the effects of three mutations on the enzymatic activity of sulfamidase, providing a molecular explanation for the cause of the disease. Thus, this study allows for a better understanding of the effect of SGSH mutations through the use of various computational approaches in terms of both structure and functions and provides a platform for the development of therapeutic drugs and potential disease treatments.

A natural anticancer pigment,Pheophytin a,from a seagrass acts as a high affinity human mitochondrial translocator protein (TSPO) ligand, in silico, to reduce mitochondrial membrane Potential (∆ψmit) in adenocarcinomic A549 cells

BACKGROUND

The present investigation looks at the most likely possibilities of usage of a naturally occurring photosynthetic pigment, Pheophytin a, from the seagrass, Syringodium isoetifolium, for plausible use as human TSPO ligand.

METHODS

Pheophytin a isolated in our laboratory previously was administered to A549 cell lines in vitro to examine its effects on cell migrations, DNA, cell cycle, Mitochondrial Membrane Potential and gene expressions. In silico tools were used to predict the nature of the compound and target binding.

RESULTS

Pheophytin a hadIC₅₀ values of 22.9 ± 5.8 µM for cancerous A549 cell lines, whilst not targeting non-cancerous vero cells [IC₅₀: 183.6 ± 1.92 µM]. Pheophytin a hindered cellular migration, fragmented DNA, arrested cell cycle precisely at S phase, reduced ∆ψ_mit and directed mRNA expressions toward apoptosis. In silico tools indicate that the compound binds to TSPO with high effectiveness to collapse ∆ψ_mit(which is proved using wet lab experiments) to promote mitophagy.

CONCLUSION

Hence Pheophytin a could be seen as a possible TSPO ligand for targeting metastatic alveolar cancers like A549 via intrinsic apoptotic pathway.

GENERAL SIGNIFICANCE

Given the inherent non-toxic nature of the compound and easy extractability from almost all autotrophic eukaryotes, one could be confident to testing in animal models.

Understanding the structure-function relationship of HPRT1 missense mutations in association with Lesch-Nyhan disease and HPRT1-related gout by in silico mutational analysis

The nucleotide salvage pathway is used to recycle degraded nucleotides (purines and pyrimidines); one of the enzymes that helps to recycle purines is hypoxanthine guanine phosphoribosyl transferase 1 (HGPRT1). Therefore, defects in this enzyme lead to the accumulation of DNA and nucleotide lesions and hence replication errors and genetic disorders. Missense mutations in hypoxanthine phosphoribosyl transferase 1 (HPRT1) are associated with deficiencies such as Lesch-Nyhan disease and chronic gout, which have manifestations such as arthritis, neurodegeneration, and cognitive disorders. In the present study, we collected 88 non-synonymous single nucleotide polymorphisms (nsSNPs) from the UniProt, dbSNP, ExAC, and ClinVar databases. We used a series of sequence-based and structure-based in silico tools to prioritize and characterize the most pathogenic and stabilizing or destabilizing nsSNPs. Moreover, to obtain the structural impact of the pathogenic mutations, we mapped the mutations to the crystal structure of the HPRT protein. We further subjected these mutant proteins to a 50 ns molecular dynamics simulation (MDS). The MDS trajectory showed that all mutant proteins altered the structural conformation and dynamic behavior of the HPRT protein and corroborated its association with LND and gout. This study provides essential information regarding the use of HPRT protein mutants as potential targets for therapeutic development.

Computational modelling approaches as a potential platform to understand the molecular genetics association between Parkinson's and Gaucher diseases

Gaucher's disease (GD) is a genetic disorder in which glucocerebroside accumulates in cells and specific organs. It is broadly classified into type I, type II and type III. Patients with GD are at high risk of Parkinson's disease (PD), and the clinical and pathological presentation of GD patients with PD is almost identical to idiopathic PD. Several experimental models like cell culture, animal models, and transgenic mice models were used to understand the molecular mechanism behind GD and PD association; however, such mechanism remains unclear. In this context, based on literature reports, we identified the most common mutations K198T, E326K, T369M, N370S, V394L, D409H, L444P, and R496H, in the Glucosylceramidase (GBA) protein that are known to cause GD1, and represent a risk of developing PD. However, to date, no computational analyses have designed to elucidate the potential functional role of GD mutations with increased risk of PD. The present computational pipeline allows us to understand the structural and functional significance of these GBA mutations with PD. Based on the published data, the most common and severe mutations were E326K, N370S, and L444P, which further selected for our computational analysis. PredictSNP and iStable servers predicted L444P mutant to be the most deleterious and responsible for the protein destabilization, followed by the N370S mutation. Further, we used the structural analysis and molecular dynamics approach to compare the most frequent deleterious mutations (N370S and L444P) with the mild mutation E326K. The structural analysis demonstrated that the location of E326K and N370S in the alpha helix region of the protein whereas the mutant L444P was in the starting region of the beta sheet, which might explain the predicted pathogenicity level and destabilization effect of the L444P mutant. Finally, Molecular Dynamics (MD) at 50 ns showed the highest deviation and fluctuation pattern in the L444P mutant compared to the two mutants E326K and N370S and the native protein. This was consistent with more loss of intramolecular hydrogen bonds and less compaction of the radius of gyration in the L444P mutant. The proposed study is anticipated to serve as a potential platform to understand the mechanism of the association between GD and PD, and might facilitate the process of drug discovery against both GD and PD.

Impact of missense mutations in survival motor neuron protein (SMN1) leading to Spinal Muscular Atrophy (SMA): A computational approach

Spinal muscular atrophy (SMA) is a neuromuscular disorder caused by the mutations in survival motor neuron 1 gene (SMN1). The molecular pathology of missense mutations in SMN1 is not thoroughly investigated so far. Therefore, we collected all missense mutations in the SMN1 protein, using all possible search terms, from three databases (PubMed, PMC and Google Scholar). All missense mutations were subjected to in silico pathogenicity, conservation, and stability analysis tools. We used statistical analysis as a QC measure for validating the specificity and accuracy of these tools. PolyPhen-2 demonstrated the highest specificity and accuracy. While PolyPhen-1 showed the highest sensitivity; overall, PolyPhen2 showed better measures in comparison to other in silico tools. Three mutations (D44V, Y272C, and Y277C) were identified as the most pathogenic and destabilizing. Further, we compared the physiochemical properties of the native and the mutant amino acids and observed loss of H-bonds and aromatic stacking upon the cysteine to tyrosine substitution, which led to the loss of aromatic rings and may reduce protein stability. The three mutations were further subjected to Molecular Dynamics Simulation (MDS) analysis using GROMACS to understand the structural changes. The Y272C and Y277C mutants exhibited maximum deviation pattern from the native protein as compared to D44V mutant. Further MDS analysis predicted changes in the stability that may have been contributed due to the loss of hydrogen bonds as observed in intramolecular hydrogen bond analysis and physiochemical analysis. A loss of function/structural impact was found to be severe in the case of Y272C and Y277C mutants in comparison to D44V mutation. Correlating the results from in silico predictions, physiochemical analysis, and MDS, we were able to observe a loss of stability in all the three mutants. This combinatorial approach could serve as a platform for variant interpretation and drug design for spinal muscular dystrophy resulting from missense mutations.

Computational insights of K1444N substitution in GAP-related domain of NF1 gene associated with neurofibromatosis type 1 disease: a molecular modeling and dynamics approach

The NF1 gene encodes for neurofibromin protein, which is ubiquitously expressed, but most highly in the central nervous system. Non-synonymous SNPs (nsSNPs) in the NF1 gene were found to be associated with Neurofibromatosis Type 1 disease, which is characterized by the growth of tumors along nerves in the skin, brain, and other parts of the body. In this study, we used several in silico predictions tools to analyze 16 nsSNPs in the RAS-GAP domain of neurofibromin, the K1444N (K1423N) mutation was predicted as the most pathogenic. The comparative molecular dynamic simulation (MDS; 50 ns) between the wild type and the K1444N (K1423N) mutant suggested a significant change in the electrostatic potential. In addition, the RMSD, RMSF, Rg, hydrogen bonds, and PCA analysis confirmed the loss of flexibility and increase in compactness of the mutant protein. Further, SASA analysis revealed exchange between hydrophobic and hydrophilic residues from the core of the RAS-GAP domain to the surface of the mutant domain, consistent with the secondary structure analysis that showed significant alteration in the mutant protein conformation. Our data concludes that the K1444N (K1423N) mutant lead to increasing the rigidity and compactness of the protein. This study provides evidence of the benefits of the computational tools in predicting the pathogenicity of genetic mutations and suggests the application of MDS and different in silico prediction tools for variant assessment and classification in genetic clinics.