Elucidating the functional impact of G137V and G144R variants in Maroteaux Lamy's Syndrome by Molecular Dynamics Simulation.

Mucopolysaccharidoses VI (Maroteaux Lamy syndrome) is a metabolic disorder due to the loss of enzyme activity of N-acetyl galactosamine-4-sulphatase arising from mutations in the ARSB gene. The mutated ARSB is the origin for the accumulation of GAGs within the lysosome leading to severe growth deformities, causing lysosomal storage disease. The main focus of this study is to identify the deleterious variants by applying bioinformatics tools to predict the conservation, pathogenicity, stability, and effect of the ARSB variants. We examined 170 missense variants, of which G137V and G144R were the resultant variants predicted detrimental to the progression of the disease. The native along with G137V and G144R structures were fixed as the receptors and subjected to Molecular docking with the small molecule Odiparcil to analyze the binding efficiency and the varied interactions of the receptors towards the drug. The interaction resulted in similar docking scores of - 7.3 kcal/mol indicating effective binding and consistent interactions of the drug with residues CYS117, GLN118, THR182, and GLN517 for native, along with G137V and G144R structures. Molecular Dynamics were conducted to validate the stability and flexibility of the native and variant structures on ligand binding. The overall study indicates that the drug has similar therapeutic towards the native and variant based on the higher binding affinity and also the complexes show stability with an average of 0.2 nm RMS value. This can aid in the future development therapeutics for the Maroteaux Lamy syndrome.

Computational screening and structural analysis of Gly201Arg and Gly201Asp missense mutations in human cyclin-dependent kinase 4 protein.

The regulatory proteins, cyclins, and cyclin-dependent kinases (CDKs) control the cell cycle progression. CDK4 gene mutations are associated with certain cancers such as melanoma, breast cancer, and rhabdomyosarcoma. Therefore, understanding the mechanisms of cell cycle control and cell proliferation is essential in developing cancer treatment regimens. In this study, we obtained cancer-causing CDK4 mutations from the COSMIC database and subjected them to a series of in silico analyses to identify the most significant mutations. An overall of 238 mutations (119 missense mutations) retrieved from the COSMIC database were investigated for the pathogenic and destabilizing properties using the PredictSNP and iStable algorithms. Further, the amino acid position of the most pathogenic and destabilizing mutations were analyzed to understand the nature of amino acid conservation across the species during the evolution. We observed that the missense mutations G201R and G201D were more significant and the Glycine at position 201 was found to highly conserved. These significant mutations were subjected to molecular dynamics simulation analysis to understand the protein's structural changes. The results from molecular dynamics simulations revealed that both G201R and G201D of CDK4 are capable of altering the protein's native form. On comparison among the most significant mutations, G201R disrupted the protein structure higher than the protein with G201D.

Controlling cell proliferation by targeting cyclin-dependent kinase 6 using drug repurposing approach.

Cyclin-dependent kinase 6 (CDK6) is an essential kinase in cell cycle progression, which is a viable target for inhibitors in various malignancies, including breast cancer. This study aimed to virtually screen efficient compounds as new leads in treating breast cancer using a drug repurposing approach. Apoptosis regulatory compounds were taken from the seleckchem database. Molecular docking experiments were carried out in the presence of abemaciclib, a routinely used FDA drug. Compared to conventional drugs, the two compounds demonstrated a higher binding affinity for CDK6. Compounds (N-benzyl-6-[(4-hydroxyphenyl)methyl]-8-(naphthalen-1-ylmethyl)-4,7-dioxo-3,6,9,9a-tetrahydro-2H-pyrazino[1,2-a]pyrimidine-1-carboxamide) and (1'-[4-[1-(4-fluorophenyl)indol-3-yl]butyl]spiro[1H-2-benzofuran-3,4'-piperidine]) were discovered to have an inhibitory effect against CDK6 at -8.49 and -6.78kcal/mol, respectively, compared to -8.09kcal/mol of the control molecule, the interacting residues of these two new compounds were found to fall within the binding site of the CDK6 molecule. Both compounds exhibited equal ADME features compared with abemaciclib and would be well distributed and metabolized by the body with an appropriate druglikeness range. Lastly, molecular dynamics was initiated for 200ns for the selected potent inhibitors and abemaciclib as complexed with CDK6. The RMSD, RMSF, Rg, H-Bond interactions, SASA, PCA, FEL, and MM/PBSA analysis were performed for the complexes to assess the stability, fluctuations, radius of gyration, hydrogen bond interaction, solvent accessibility, essential dynamics, free energy landscape, and MM/PBSA. The selected two compounds are small molecules in the appropriate druglikeness range. The results observed in molecular docking and molecular dynamics simulations were most promising for two compounds, suggesting their potent inhibitory effect against CDK6. We propose that these candidate compounds can undergo in vitro validation and in vivo testing for their further use against cancer.

Deciphering the effect of mutations in MMAA protein causing methylmalonic acidemia-A computational approach.

Methylmalonic acidemia (MMA) is a rare genetic disorder affecting multiple body systems. We aimed to investigate the pathogenic mutations in MMAA that are associated with isolated methylmalonic acidemia to identify the structural behavior of MMAA upon mutation. The algorithms such as PredictSNP, iStable, ConSurf, and Align GVGD were employed to analyze the consequence of the mutations. Molecular docking was carried out for the native MMAA, L89P, G274D, and R359G to interpret its interactions with the GDP substrate. The docked complexes were simulated for 200ns aiding GROMACS in apprehending the behavior of MMAA upon mutation and GDP binding. After simulation, cα disruptions were observed using the RMSF plot, which indicated that several regions of mutant MMAAs have highly fluctuated. The gyration and H-bond plots were used to understand the compactness and intermolecular interaction with the GDP molecule. The MDS analysis showed that the mutations L89P, G274D, and R359G are highly unstable even after GDP binding, with the least compactness, fewer H-bonds, and larger conformational cα motions. Our study provided structural and dynamic insights into MMAA protein, which further helps to characterize these mutants and provide potential treatment strategies for MMA patients.

Understanding the activating mechanism of the immune system against COVID-19 by Traditional Indian Medicine: Network pharmacology approach.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) transmissions are occurring rapidly; it is raising the alarm around the globe. Though vaccines are currently available, the evolution and mutations in the SARS-CoV-2 threaten available vaccines' significance. The drugs are still undergoing clinical trials, and certain medications are approved for "emergency use" or as an "off-label" drug during the pandemic. These drugs have been effective yet accommodating side effects, which also can be lethal. Complementary and alternative medicine is highly demanded since it embraces a holistic approach. Since ancient times, natural products have been used as drugs to treat various diseases in the medical field and are still widely practiced. Medicinal plants contain many active compounds that serve as the key to an effective drug design. The Kabasura kudineer and Nilavembu kudineer are the two most widely approved formulations to treat COVID-19. However, the mechanism of these formulations is not well known. The proposed study used a network pharmacology approach to understand the immune-boosting mechanism by the Kabasura kudineer, Nilavembu kudineer, and JACOM in treating COVID-19. The plants and phytochemical chemical compounds in the Kabasura kudineer, Nilavembu kudineer, and JACOM were obtained from the literature. The Swiss target prediction algorithm was used to predict the targets for these phytochemical compounds. The common genes for the COVID-19 infection and the drug targets were identified. The gene-gene interaction network was constructed to understand the interactions between these common genes and enrichment analyses to determine the biological process, molecular functions, cellular functions, pathways involved, etc. Finally, virtual screening and molecular docking studies were performed to identify the most potential targets and significant phytochemical compounds to treat the COVID-19. The present study identified potential targets as ACE, Cathepsin L, Cathepsin B, Cathepsin K, DPP4, EGFR, HDAC2, IL6, RIPK1, and VEGFA. Similarly, betulinic acid, 5″-(2⁗-Hydroxybenzyl) uvarinol, antofine, (S)-1'-methyloctyl caffeate, (Z)-3-phenyl-2-propenal, 7-oxo-10α-cucurbitadienol, and PLX-4720 collectively to be potential treatment agents for COVID-19.

An integrative analysis to distinguish between emphysema (EML) and alpha-1 antitrypsin deficiency-related emphysema (ADL)-A systems biology approach.

Lung Emphysema is an abnormal enlargement of the air sacs followed by the destruction of alveolar walls without any prominent fibrosis. This study primarily identifies the differentially expressed genes (DEGs), interactions between them, and their significant involvement in the activated signaling cascades. The dataset with ID GSE1122 (five normal lung tissue samples, five of usual emphysema, and five of alpha-1 antitrypsin deficiency-related emphysema) from the gene expression omnibus (GEO) was analyzed using the GEO2R tool. The physical association between the DEGs were mapped using the STRING tool and was visualized in the Cytoscape software. The enriched functional processes were identified with the ClueGO plugin's help from Cytoscape. Further integrative functional annotation was performed by implying the GeneGo Metacore™ to distinguish the enriched pathway maps, process networks, and GO processes. The results from this analysis revealed the critical signaling cascades that have been either activated or inhibited due to identified DEGs. We found the activated pathways such as immune response IL-1 signaling pathway, positive regulation of smooth muscle migration, BMP signaling pathway, positive regulation of leukocyte migration, NIK/NF-kappB signaling, and cytochrome-c oxidase activity. Finally, we mapped four crucial genes (CCL5, ALK, TAC1, CD74, and HLA-DOA) by comparing the functional annotations that could be significantly influential in emphysema molecular pathogenesis. Our study provides insights into the pathogenesis of emphysema and helps in developing potential drug targets against emphysema.

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.

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.

A computational approach to analyse the amino acid variants of GLB1 protein causing GM1 Gangliosidosis.

Lysosomal storage diseases comprise different forms of autosomal recessive disorders from which GM1 gangliosidosis has categorized by the accumulation of complex glycolipids associated with a range of progressive neurologic phenotypes. GM1 gangliosidosis is an inherited disorder that progressively destroys nerve cells (neurons) in the brain and spinal cord. GM1 has three main types of onsets, namely infantile (type I), juvenile (type II), and adult (type III) forms. This study provides a series of computational methods that examine the mutations that occurred in GLB1 protein. Initially, the mutational analysis started with 689 amino acid variants for a sequence-based screening and it was done with quite a few In-silico tools to narrow down the most significant variants by utilizing the standard tools; namely, Evolutionary analysis (77 variants), Pathogenicity prediction (44 variants), Stability predictions (30 variants), Biophysical functions (19 variants) and according to the binding site of protein structure with PDB ID 3THC, seven variants (Y83D, Y83H, Y270S, Y270D, W273R, W273D, and Y333H) were narrowed down. Structure based analysis was performed to understand the interacting profile of the native protein and variants with Miglustat; which is the currently used FDA drug as an alternative to enzyme replacement therapy. Molecular Docking study was done to analyze the protein interaction with Miglustat (ligand), as a result native (3THC) structure had a binding affinity of -8.18 kcal/mol and two variant structures had an average binding affinities of -2.61 kcal/mol (Y83D) and - 7.63 kcal/mol (Y270D). Finally, Molecular Dynamics Simulation was performed to know the mutational activity of the protein structures on Miglustat for 50,000 ps. The Y83D variant showed higher deviation than native protein and Y270D in all trajectory analysis. The analysis was done to the protein structures to check the structural variations happened through simulations. This study aids to understand the most deleterious mutants, the activity of the drug to the protein structure and also gives an insight on the stability of the drug with the native and selected variants.

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.

Computational model to analyze and characterize the functional mutations of NOD2 protein causing inflammatory disorder - Blau syndrome.

Blau syndrome (BS), which affects the eyes, skin, and joints, is an autosomal dominant genetic inflammatory disorder. BS is caused by mutations in the NOD2 gene. However, there are no direct treatments, and treatment with conventional anti-inflammatory drugs such as adrenal glucocorticoids, anti-metabolites, and biological agents such as anti-TNF and infliximab have all been attempted with varying degrees of success. In this study, we tried to identify all the reported mutations in the NOD2 protein that cause BS. Collectively, 114 missense mutations were extracted from the UniProt, ClinVar, and HGMD databases. The mutations were further subjected to pathogenic, stability, and conservation analyses. According to these computational analyses, six missense mutations (R334Q, R334W, E383G, E383K, R426H, and T605P) were found to be highly deleterious, destabilizing, and positioned in the conserved position. ADP to ATP conversion plays a crucial role in switching the closed-form of NOD2 protein to the open-form, thus activating the protein. Accordingly, the mutations in the ADP binding sites have received more attention in comparison to the mutations in the non-ADP binding positions. Interestingly, the W490L mutation is positioned in the ADP binding site and exhibits highly deleterious and destabilizing properties. Additionally, W490L was also found to be conserved, with a ConSurf score of 7. Therefore, we further performed homology modeling to determine the 3D structure of native NOD2 and the W490L mutant. Molecular docking analysis was carried out to understand the change in the interaction of ADP with the NOD2 protein. We observed that ADP had a stronger interaction with the native NOD2 protein compared to the W490L mutant. Finally, ADP complexed with native NOD2 and W490L mutant were subjected to molecular dynamics simulations, and the trajectories were analyzed. In the simulations, we observed decreased deviation and fluctuations in native NOD2, whereas decreased compactness and inter- and intramolecular hydrogen bonds were observed in the W490L mutant. This study is expected to serve as a platform for developing targeted drug therapy for BS.

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 hadIC50 values of 22.9 ±â€¯5.8 µM for cancerous A549 cell lines, whilst not targeting non-cancerous vero cells [IC50: 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.

Molecular insights of the G2019S substitution in LRRK2 kinase domain associated with Parkinson's disease: A molecular dynamics simulation approach.

The G2019S substitution in the Leucine-rich repeat kinase 2 (LRRK2) is significantly associated with Parkinson's disease (PD). This substitution was identified in both familial and sporadic forms of PD with a higher frequency. Few computational studies have reported the impact of G2019S substitution on inhibitors of the kinase domain of LRRK2. However, no computational study deeply investigated the possible impact of the G2019S substitution on the kinase domain in its Apo conformation. Therefore, in this study, we used 200 ns molecular dynamic simulation using the GROMACS 5.1.4 package software to investigate the impact of the G2019S substitution on the structure of the kinase domain of LRRK2. Our results indicate that the G2019S substitution affects the dynamics and stability of LRRK2 by decreasing the flexibility and increasing the compactness of the kinase domain and showing its tendency to be in an active conformation for long time interval because of the high energy barrier between active and inactive conformation. This study predicts the molecular pathogenicity mechanism of the G2019S on patients with PD and provides a potential platform for developing therapeutics for patients with PD that harbor this amino acid substitution.

Analyzing the Effect of V66M Mutation in BDNF in Causing Mood Disorders: A Computational Approach.

Mental disorders or mood disorders are prevalent globally irrespective of region, race, and ethnic groups. Of the types of mood disorders, major depressive disorder (MDD) and bipolar disorder (BPD) are the most prevalent forms of psychiatric condition. A number of preclinical studies emphasize the essential role of brain-derived neurotrophic factor (BDNF) in the pathophysiology of mood disorders. Additionally, BDNF is the most common growth factor in the central nervous system along with their essential role during the neural development and the synaptic elasticity. A malfunctioning of this protein is associated with many types of mood disorders. The variant methionine replaces valine at 66th position is strongly related to BPD, and an individual with a homozygous condition of this allele is at a greater risk of developing MDD. There are very sparse reports suggesting the structural changes of the protein occurring upon the mutation. Consequently, in this study, we applied a computational pipeline to understand the effects caused by the mutation on the protein's structure and function. With the use of in silico tools and computational macroscopic methods, we identified a decrease in the alpha-helix nature, and an overall increase in the random coils that could have probably resulted in deformation of the protein.

Influence of V54M mutation in giant muscle protein titin: a computational screening and molecular dynamics approach.

Recent genetic studies have revealed the impact of mutations in associated genes for cardiac sarcomere components leading to dilated cardiomyopathy (DCM). The cardiac sarcomere is composed of thick and thin filaments and a giant muscle protein known as titin or connectin. Titin interacts with T-cap/telethonin in the Z-line region and plays a vital role in regulating sarcomere assembly. Initially, we screened all the variants associated with giant protein titin and analyzed their impact with the aid of pathogenicity and stability prediction methods. V54M mutation found in the hydrophobic core region of the protein associated with abnormal clinical phenotype leads to DCM was selected for further analysis. To address this issue, we mapped the deleterious mutant V54M, modeled the mutant protein complex, and deciphered the impact of mutation on binding with its partner telethonin in the titin crystal structure of PDB ID: 1YA5 with the aid of docking analysis. Furthermore, two run molecular dynamics simulation was initiated to understand the mechanistic action of V54M mutation in altering the protein structure, dynamics, and stability. According to the results obtained from the repeated 50 ns trajectory files, the overall effect of V54M mutation was destabilizing and transition of bend to coil in the secondary structure was observed. Furthermore, MMPBSA elucidated that V54M found in the Z-line region of titin decreases the binding affinity of titin to Z-line proteins T-cap/telethonin thereby hindering the protein-protein interaction.

Computational identification of pathogenic associated nsSNPs and its structural impact in UROD gene: a molecular dynamics approach.

Uroporphyrinogen decarboxylase is a cytosolic enzyme involved in the biosynthetic pathway of heme production. Decreased activity of this enzyme results in porphyria cutanea tarda and hepato erythropoietic porphyria. Nonsynonymous single nucleotide polymorphisms (nsSNPs) alter protein sequence and can cause disease. Identifying the deleterious nsSNPs that contribute to disease is an important task. We used five different in silico tools namely SIFT, PANTHER, PolyPhen2, SNPs&GO, and I-mutant3 to identify deleterious nsSNPs in UROD gene. Further, we used molecular dynamic (MD) approach to evaluate the impact of deleterious mutations on UROD protein structure. By comparing the results of all the five prediction results, we screened 35 (51.47 %) nsSNPs as highly deleterious. MD analysis results show that all the three L161Q, L282R, and I334T deleterious variants were affecting the UROD protein structural stability and flexibility. Our findings provide strong evidence on the effect of deleterious nsSNPs in UROD gene. A detailed MD study provides a new insight in the conformational changes occurred in the mutant structures of UROD protein.

In silico profiling and structural insights of missense mutations in RET protein kinase domain by molecular dynamics and docking approach.

A major challenge remaining in drug design efforts towards protein kinase is due to the development of drug resistance initiated by the missense mutations in the kinase catalytic domain. Gain or loss of function mutations in the REarranged during Transfection (RET) tyrosine kinase gene have been associated with the development of a wide range of human associated cancers and Hirschsprung's disease. However, to what extent these mutations might affect bio-molecular functions remains unclear. In this article, the functionally significant mutations in RET were screened with the aid of various sequence and structure based in silico prediction methods. We mapped the deleterious mutants, modelled mutant proteins and deciphered the impact of mutations on drug binding mechanisms in the RET crystal structure of PDB ID: with the potential inhibitor vandetanib by docking analysis. Furthermore, molecular dynamics simulations were undertaken to understand the mechanistic action of cancer associated mutations in altering the protein kinase structure, dynamics, and stability. According to our results, the overall effect of V804M, M918T and S922Y were destabilizing and mostly alter the electrostatic component of the binding energy. Specifically, the mutation of gatekeeper residue valine 804 present in the ATP binding pocket affects the protein stability and confers resistance to the drug vandetanib, which was consistent with previously published experimental results. Overall, our findings may provide useful structural insights for in-depth understanding of the molecular mechanism underlying RET mutation and developing effective drugs.

Computational pipeline to identify and characterize functional mutations in ornithine transcarbamylase deficiency.

Ornithine transcarbamylase (OTC) (E.C. 2.1.3.3) is one of the enzymes in the urea cycle, which involves in a sequence of reactions in the liver cells. During protein assimilation in our body surplus nitrogen is made, this open nitrogen is altered into urea and expelled out of the body by kidneys, in this cycle OTC helps in the conversion of free toxic nitrogen into urea. Ornithine transcarbamylase deficiency (OTCD: OMIM#311250) is triggered by mutation in this OTC gene. To date more than 200 mutations have been noted. Mutation in OTC gene indicates alteration in enzyme production, which upsets the ability to carry out the chemical reaction. The computational analysis was initiated to identify the deleterious nsSNPs in OTC gene in causing OTCD using five different computational tools such as SIFT, PolyPhen 2, I-Mutant 3, SNPs&Go, and PhD-SNP. Studies on the molecular basis of OTC gene and OTCD have been done partially till date. Hence, in silico categorization of functional SNPs in OTC gene can provide valuable insight in near future in the diagnosis and treatment of OTCD.

Computational methods to work as first-pass filter in deleterious SNP analysis of alkaptonuria.

A major challenge in the analysis of human genetic variation is to distinguish functional from nonfunctional SNPs. Discovering these functional SNPs is one of the main goals of modern genetics and genomics studies. There is a need to effectively and efficiently identify functionally important nsSNPs which may be deleterious or disease causing and to identify their molecular effects. The prediction of phenotype of nsSNPs by computational analysis may provide a good way to explore the function of nsSNPs and its relationship with susceptibility to disease. In this context, we surveyed and compared variation databases along with in silico prediction programs to assess the effects of deleterious functional variants on protein functions. In other respects, we attempted these methods to work as first-pass filter to identify the deleterious substitutions worth pursuing for further experimental research. In this analysis, we used the existing computational methods to explore the mutation-structure-function relationship in HGD gene causing alkaptonuria.