Computational SNP Analysis and Molecular Simulation Revealed the Most Deleterious Missense Variants in the NBD1 Domain of Human ABCA1 Transporter

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

BACKGROUND

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

MATERIALS AND METHODS

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

RESULTS

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

CONCLUSIONS

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

Computational exploration of SLC14A1 genetic variants through structure modeling, protein-ligand docking, and molecular dynamics simulation

The urea transporter UT-B1, encoded by the SLC14A1 gene, has been hypothesized to be a significant protein whose deficiency and dysfunction contribute to the pathogenesis of bladder cancer and many other diseases. Several studies reported the association of genetic alterations in the SLC14A1 (UT-B1) gene with bladder carcinogenesis, suggesting a need for thorough characterization of the UT-B1 protein's coding and non-coding variants. This study used various computational techniques to investigate the commonly occurring germ-line missense and non-coding SNPs (ncSNPs) of the SLC14A1 gene (UT-B1) for their structural, functional, and molecular implications for disease susceptibility and dysfunctionality. SLC14A1 missense variants, primarily identified from the ENSEMBL genome browser, were screened through twelve functionality prediction tools leading to two variants D280Y (predicted detrimental by maximum tools) and D280N (high global MAF) for rs1058396. Subsequently, the ConSurf and NetSurf tools revealed the D280 residue to be in a variable site and exposed on the protein surface. According to I-Mutant2.0 and MUpro, both variants are predicted to cause a significant effect on protein stability. Analysis of molecular docking anticipated these two variants to decrease the binding affinity of UT-B1 protein for the examined ligands to a significant extent. Molecular dynamics also disclosed the possible destabilization of the UT-B1 protein due to single nucleotide polymorphism compared to wild-type protein which may result in impaired protein function. Furthermore, several non-coding SNPs were estimated to affect transcription factor binding and regulation of SLC14A1 gene expression. Additionally, two ncSNPs were found to affect miRNA-based post-transcriptional regulation by creating new seed regions for miRNA binding. This comprehensive in-silico study of SLC14A1 gene variants may serve as a springboard for future large-scale investigations examining SLC14A1 polymorphisms.

Comprehensive bioinformatics analysis of structural and functional consequences of deleterious missense mutations in the human QDPR gene

Quinonoid dihydropteridine reductase (QDPR) is an enzyme that regulates tetrahydrobiopterin (BH4), a cofactor for enzymes involved in neurotransmitter synthesis and blood pressure regulation. Reduced QDPR activity can cause dihydrobiopterin (BH2) accumulation and BH4 depletion, leading to impaired neurotransmitter synthesis, oxidative stress, and increased risk of Parkinson's disease. A total of 10,236 SNPs were identified in the QDPR gene, with 217 being missense SNPs. Over 18 different sequence-based and structure-based tools were employed to assess the protein's biological activity, with several computational tools identifying deleterious SNPs. Additionally, the article provides detailed information about the QDPR gene and protein structure and conservation analysis. The results showed that 10 mutations were harmful and linked to brain and central nervous system disorders, and were predicted to be oncogenic by Dr. Cancer and CScape. Following conservation analysis, the HOPE server was used to analyse the effect of six selected mutations (L14P, V15G, G23S, V54G, M107K, G151S) on the protein structure. Overall, the study provides insights into the biological and functional impact of nsSNPs on QDPR activity and the potential induced pathogenicity and oncogenicity. In the future, research can be conducted to systematically evaluate QDPR gene variation through clinical studies, investigate mutation prevalence across different geographical regions, and validate computational results with conclusive experiments.Communicated by Ramaswamy H. Sarma.

Computational Analysis of Non-synonymous SNPs in ATM Kinase: Structural Insights, Functional Implications, and Inhibitor Discovery

Ataxia telangiectasia-mutated (ATM) protein kinase, a key player in cellular integrity regulation, is known for its role in DNA damage response. This study investigates the broader impact of ATM on cellular processes and potential clinical manifestations arising from mutations, aiming to expand our understanding of ATM's diverse functions beyond conventional roles. The research employs a comprehensive set of computational techniques for a thorough analysis of ATM mutations. The mutation data are curated from dbSNP and HuVarBase databases. A meticulous assessment is conducted, considering factors such as deleterious effects, protein stability, oncogenic potential, and biophysical characteristics of the identified mutations. Conservation analysis, utilizing diverse computational tools, provides insights into the evolutionary significance of these mutations. Molecular docking and dynamic simulation analyses are carried out for selected mutations, investigating their interactions with Y2080D, AZD0156, and quercetin inhibitors to gauge potential therapeutic implications. Among the 419 mutations scrutinized, five (V1913C, Y2080D, L2656P, C2770G, and C2930G) are identified as both disease causing and protein destabilizing. The study reveals the oncogenic potential of these mutations, supported by findings from the COSMIC database. Notably, Y2080D is associated with haematopoietic and lymphoid cancers, while C2770G shows a correlation with squamous cell carcinomas. Molecular docking and dynamic simulation analyses highlight strong binding affinities of quercetin for Y2080D and AZD0156 for C2770G, suggesting potential therapeutic options. In summary, this computational analysis provides a comprehensive understanding of ATM mutations, revealing their potential implications in cellular integrity and cancer development. The study underscores the significance of Y2080D and C2770G mutations, offering valuable insights for future precision medicine targeting-specific ATM. Despite informative computational analyses, a significant research gap exists, necessitating essential in vitro and in vivo studies to validate the predicted effects of ATM mutations on protein structure and function.

Functional characterization of missense variants affecting the extracellular domains of ABCA1 using a fluorescence-based assay

Excess cholesterol originating from nonhepatic tissues is transported within HDL particles to the liver for metabolism and excretion. Cholesterol efflux is initiated by lipid-free or lipid-poor apolipoprotein A1 interacting with the transmembrane protein ABCA1, a key player in cholesterol homeostasis. Defective ABCA1 results in reduced serum levels of HDL cholesterol, deposition of cholesterol in arteries, and an increased risk of early onset CVD. Over 300 genetic variants in ABCA1 have been reported, many of which are associated with reduced HDL cholesterol levels. Only a few of these have been functionally characterized. In this study, we have analyzed 51 previously unclassified missense variants affecting the extracellular domains of ABCA1 using a sensitive, easy, and low-cost fluorescence-based assay. Among these, only 12 variants showed a distinct loss-of-function phenotype, asserting their direct association with severe HDL disorders. These findings emphasize the crucial role of functional characterization of genetic variants in pathogenicity assessment and precision medicine. The functional rescue of ABCA1 loss-of-function variants through proteasomal inhibition or by the use of the chemical chaperone 4-phenylbutyric acid was genotype specific. Genotype-specific responses were also observed for the ability of apolipoprotein A1 to stabilize the different ABCA1 variants. In view of personalized medicine, this could potentially form the basis for novel therapeutic strategies.

Identification of Important Genes Associated with the Development of Atherosclerosis

Atherosclerosis is one of the most important medical problems due to its prevalence and significant contribution to the structure of temporary and permanent disability and mortality. Atherosclerosis is a complex chain of events occurring in the vascular wall over many years. Disorders of lipid metabolism, inflammation, and impaired hemodynamics are important mechanisms of atherogenesis. A growing body of evidence strengthens the understanding of the role of genetic and epigenetic factors in individual predisposition and development of atherosclerosis and its clinical outcomes. In addition, hemodynamic changes, lipid metabolism abnormalities, and inflammation are closely related and have many overlapping links in regulation. A better study of these mechanisms may improve the quality of diagnosis and management of such patients.

Paradoxical effect of Aβ on protein levels of ABCA1 in astrocytes, microglia, and neurons isolated from C57BL/6 mice: an in vitro and in silico study to elucidate the effect of Aβ on ABCA1 in the brain cells

Impaired cholesterol metabolism has been reported in Alzheimer's disease. Since ABCA1 is one of the main players in the brain's cholesterol homeostasis, here we used the in-vitro and in-silico experiments to investigate the effect of Aβ on ABCA1 protein levels in microglia, astrocytes, and neurons in mice. Microglia, astrocytes, and neurons were cultured and exposed to beta amyloid. ABCA1 in cell lysates was determined by Western blotting, and cholesterol efflux was measured in the conditioned media. Molecular docking, molecular dynamics simulations, and MM-GBSA analysis were conducted to gain a better understanding of the effects of Aβ on ABCA1. In response to Aβ, the protein levels of ABCA1 increase significantly in microglia, astrocytes, and neurons; however, its ability to enhance cholesterol efflux is diminished. Aβ inhibited the function of ABCA1 by obstructing the extracellular tunnel that transports lipids outside the cell, as determined by molecular docking. MD simulation analysis validated these findings. Our results demonstrated that Aβ could increase ABCA1 protein levels in various brain cells, regardless of cell type. Molecular docking, molecular dynamics simulation, and MM-GBSA studies indicate that Aβ has a significant effect on the structural conformation of ABCA1, possibly interfering with its function. We believe that the conformational changes of ABCA1 will inhibit its ability to subsequently release cellular cholesterol. Aβ may obstruct the extracellular tunnel of ABCA1, rendering it less accessible to proteases such as the calpain family, which may explain the increase in ABCA1 levels but decrease in its function.Communicated by Ramaswamy H. Sarma.

Computational insights into NIMA-related kinase 6: unraveling mutational effects on structure and function

The NEK6 (NIMA-related kinase 6) serine/threonine kinase is a pivotal player in a multitude of cellular processes, including the regulation of the cell cycle and the response to DNA damage. Its significance extends to disease pathogenesis, as changes in NEK6 activity have been linked to the development of cancer. Non-synonymous single nucleotide polymorphisms (nsSNPs) in NEK6 have been linked to cancer as they alter the protein's native structure and function. The association between NEK6 activity and cancer development has prompted researchers to explore the effects of genetic variations within the NEK6 gene. Therefore, we utilized advanced computational tools to analyze 155 high-confidence nsSNPs in the NEK6 gene. From this analysis, 21 nsSNPs were identified as potentially harmful, raising concerns about their impact on NEK6 activity and cancer risk. These 21 mutations were then examined for structural alterations, and eight of nsSNPs (I51M, V76A, I134N, Y152D, R171Q, V186G, L237R, and C285S) were found to destabilize the protein. Among the destabilizing mutations screened, a specific mutation, R171Q, stood out due to its conserved nature. To understand its impact on the protein and conformation, all-atom molecular dynamics simulations (MDS) for 100 ns were performed for both Wildtype NEK6 (WT-NEK6) and R171Q. The simulations revealed that the R171Q variant was unstable and led to significant conformational changes in NEK6. This study provides valuable insights into NEK6 dysfunction caused by single amino acid alterations, offering a novel understanding of the molecular mechanisms underlying NEK6-related cancer progression.

Deciphering structural aspects of reverse cholesterol transport: mapping the knowns and unknowns

Atherosclerosis is a major contributor to the onset and progression of cardiovascular disease (CVD). Cholesterol-loaded foam cells play a pivotal role in forming atherosclerotic plaques. Induction of cholesterol efflux from these cells may be a promising approach in treating CVD. The reverse cholesterol transport (RCT) pathway delivers cholesteryl ester (CE) packaged in high-density lipoproteins (HDL) from non-hepatic cells to the liver, thereby minimising cholesterol load of peripheral cells. RCT takes place via a well-organised interplay amongst apolipoprotein A1 (ApoA1), lecithin cholesterol acyltransferase (LCAT), ATP binding cassette transporter A1 (ABCA1), scavenger receptor-B1 (SR-B1), and the amount of free cholesterol. Unfortunately, modulation of RCT for treating atherosclerosis has failed in clinical trials owing to our lack of understanding of the relationship between HDL function and RCT. The fate of non-hepatic CEs in HDL is dependent on their access to proteins involved in remodelling and can be regulated at the structural level. An inadequate understanding of this inhibits the design of rational strategies for therapeutic interventions. Herein we extensively review the structure-function relationships that are essential for RCT. We also focus on genetic mutations that disturb the structural stability of proteins involved in RCT, rendering them partially or completely non-functional. Further studies are necessary for understanding the structural aspects of RCT pathway completely, and this review highlights alternative theories and unanswered questions.

Non-synonymous SNPs variants of PRKCG and its association with oncogenes predispose to hepatocellular carcinoma

BACKGROUND

PRKCG encodes PKC γ, which is categorized under the classical protein kinase C family. No studies have specifically established the relationship between PRKCG nsSNPs with structural and functional variations in PKC γ in the context of hepatocellular carcinoma (HCC). The present study aims to uncover this link through in-silico and experimental studies.

METHODS

The 3D structure of PKC γ was predicted. Molecular Dynamic (MD) Simulations were run and estimates were made for interactions, stability, conservation and post-translational alterations between wild and mutant structures. The association of PRKCG levels with HCC survival rate was determined. Genotyping analyses were conducted to investigate the deleterious PRKCG nsSNP association with HCC. mRNA expression of PKC γ, HIF-1 alpha, AKT, SOCS3 and VEGF in the blood of controls and HCC patients was analyzed and a genetic cascade was constructed depicting these interactions.

RESULTS

The expression level of studied oncogenes was compared to tumour suppressor genes. Through Alphafold, the 3D structure of PKC γ was explored. Fifteen SNPs were narrowed down for in-silico analyses that were identified in exons 5, 10 and 18 and the regulatory and kinase domain of PKC γ. Root mean square deviation and fluctuation along with the radius of gyration unveiled potential changes between the wild and mutated variant structures. Mutant genotype AA (homozygous) corresponding to nsSNP, rs386134171 had more frequency in patients with OR (2.446), RR (1.564) and P-values (< 0.0029) that highlights its significant association with HCC compared to controls in which the wild genotype GG was found more prevalent.

CONCLUSION

nsSNP rs386134171 can be a genetic marker for HCC diagnosis and therapeutic studies. This study has laid down a road map for future studies to be conducted on HCC.

Unveiling the Potentiality of Shikonin Derivatives Inhibiting SARS-CoV-2 Main Protease by Molecular Dynamic Simulation Studies

Shikonin, a phytochemical present in the roots of Lithospermum erythrorhizon, is well-known for its broad-spectrum activity against cancer, oxidative stress, inflammation, viruses, and anti-COVID-19 agents. A recent report based on a crystallographic study revealed a distinct conformation of shikonin binding to the SARS-CoV-2 main protease (M^pro), suggesting the possibility of designing potential inhibitors based on shikonin derivatives. The present study aimed to identify potential shikonin derivatives targeting the M^pro of COVID-19 by using molecular docking and molecular dynamics simulations. A total of 20 shikonin derivatives were screened, of which few derivatives showed higher binding affinity than shikonin. Following the MM-GBSA binding energy calculations using the docked structures, four derivatives were retained with the highest binding energy and subjected to molecular dynamics simulation. Molecular dynamics simulation studies suggested that alpha-methyl-n-butyl shikonin, beta-hydroxyisovaleryl shikonin, and lithospermidin-B interacted with two conserved residues, His41 and Cys145, through multiple bonding in the catalytic sites. This suggests that these residues may effectively suppress SARS-CoV-2 progression by inhibiting M^pro. Taken together, the present in silico study concluded that shikonin derivatives may play an influential role in M^pro inhibition.

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

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

The effects of missense OPN3 mutations in melanocytic lesions on protein structure and light-sensitive function

Opsin 3 (OPN3), a member of the light-sensitive, retinal-dependent opsin family, is widely expressed in a variety of human tissues and plays a multitude of light-dependent and light-independent roles. We recently identified five missense variants of OPN3, including p. I51T, p. V134A, p. V183I, p. M256I and p. C331Y, in human melanocytic tumours. However, it remains unclear how these OPN3 variants affect OPN3 protein structure and function. Herein, we conducted structural and functional studies of these variant proteins in OPN3 by molecular docking and molecular dynamics simulations. Moreover, we performed in vitro fluorescence calcium imaging to assess the functional properties of five single-nucleotide variant (SNV) proteins using a site-directed mutagenesis method. Notably, the p. I51T variant was not able to effectively dock with 11-cis-retinal. Additionally, in vitro, the p. I51T SNVs failed to induce any detectable changes in intracellular Ca²⁺ concentration at room temperature. Taken together, these results reveal that five SNVs in the OPN3 gene have deleterious effects on protein structure and function, suggesting that these mutations, especially the p. I51T variant, significantly disrupt the canonical function of the OPN3 protein. Our findings provide new insight into the role of OPN3 variants in the loss of protein function.

Genetic and Epigenetic Regulation of Lipoxygenase Pathways and Reverse Cholesterol Transport in Atherogenesis

Atherosclerosis is one of the most important medical and social problems of modern society. Atherosclerosis causes a large number of hospitalizations, disability, and mortality. A considerable amount of evidence suggests that inflammation is one of the key links in the pathogenesis of atherosclerosis. Inflammation in the vascular wall has extensive cross-linkages with lipid metabolism, and lipid mediators act as a central link in the regulation of inflammation in the vascular wall. Data on the role of genetics and epigenetic factors in the development of atherosclerosis are of great interest. A growing body of evidence is strengthening the understanding of the significance of gene polymorphism, as well as gene expression dysregulation involved in cross-links between lipid metabolism and the innate immune system. A better understanding of the genetic basis and molecular mechanisms of disease pathogenesis is an important step towards solving the problems of its early diagnosis and treatment.

Identification and functional assays of single-nucleotide variants of opsins genes in melanocytic tumors

Epidermal melanocytes sense solar light via the opsin-coupled signaling pathway which is involved in a range of biological functions, including regulating pigmentation, proliferation, apoptosis, and tumorigenesis. However, it remains unclear whether there are genetic variants within these opsins that affect opsin protein structure and function, and further melanocyte biological behaviors. Here, we examined single-nucleotide variants (SNVs) of five opsin (RGR, OPN1SW, OPN2, OPN4, and OPN5) genes in MM (malignant melanoma; n = 76) and MN (melanocytic nevi; n = 157), using next-generation sequencing. The effects of these pathogenic single-nucleotide variants (SNVs) on opsin structure and function were further investigated using molecular dynamics (MD) simulations, dynamic cross-correlation (DCC), and site-directed mutagenesis. In total, 107 SNV variants were identified. Of these variants, 14 nonsynonymous SNVs (nsSNVs) of opsin genes were detected, including three mutations in the RGR gene, three mutations in the OPN1SW gene, two mutations in the OPN2 gene, and six mutations in the OPN4 gene. The effect of these missense mutations on opsin function was then assessed using eight prediction tools to estimate the potential impact of an amino acid substitution. The impact of each nsSNV was investigated using MD simulations and DCC analysis. Furthermore, we performed in vitro fluorescence calcium imaging to assess the functional properties of nsSNV proteins using a site-directed mutagenesis method. Taken together, these results revealed that p.A103V (RGR), p.T167I (RGR), p.G141S (OPN1SW), p.R144C (OPN1SW), and p.S231F (OPN4) had more deleterious effects on protein structure and function among the 14 nsSNVs. Opsin gene alterations showed the low frequency of missense mutations in melanocytic tumors, and although rare, some mutations in these opsin genes disrupt the canonical function of opsin. Our findings provide new insight into the role of opsin variants in the loss of function.

Dynamic insights into the effects of nonsynonymous polymorphisms (nsSNPs) on loss of TREM2 function

Single nucleotide variations in Triggering Receptor Expressed on Myeloid Cells 2 (TREM2) are associated with many neurodegenerative diseases, including Nasu-Hakola disease (NHD), frontotemporal dementia (FTD), and late-onset Alzheimer's disease because they disrupt ligand binding to the extracellular domain of TREM2. However, the effects of nonsynonymous single nucleotide polymorphisms (nsSNPs) in TREM2 on disease progression remain unknown. In this study, we identified several high-risk nsSNPs in the TREM2 gene using various deleterious SNP predicting algorithms and analyzed their destabilizing effects on the ligand recognizing region of the TREM2 immunoglobulin (Ig) domain by molecular dynamics (MD) simulation. Cumulative prediction by all tools employed suggested the three most deleterious nsSNPs involved in loss of TREM2 function are rs549402254 (W50S), rs749358844 (R52C), and rs1409131974 (D104G). MD simulation showed that these three variants cause substantial structural alterations and conformational remodeling of the apical loops of the TREM2 Ig domain, which is responsible for ligand recognition. Detailed analysis revealed that these variants substantially increased distances between apical loops and induced conformation remodeling by changing inter-loop nonbonded contacts. Moreover, all nsSNPs changed the electrostatic potentials near the putative ligand-interacting region (PLIR), which suggested they might reduce specificity or loss of binding affinity for TREM2 ligands. Overall, this study identifies three potential high-risk nsSNPs in the TREM2 gene. We propose further studies on the molecular mechanisms responsible for loss of TREM2 function and the associations between TREM2 nsSNPs and neurodegenerative diseases.

Structural Consequence of Non-Synonymous Single-Nucleotide Variants in the N-Terminal Domain of LIS1

Disruptive neuronal migration during early brain development causes severe brain malformation. Characterized by mislocalization of cortical neurons, this condition is a result of the loss of function of migration regulating genes. One known neuronal migration disorder is lissencephaly (LIS), which is caused by deletions or mutations of the LIS1 (PAFAH1B1) gene that has been implicated in regulating the microtubule motor protein cytoplasmic dynein. Although this class of diseases has recently received considerable attention, the roles of non-synonymous polymorphisms (nsSNPs) in LIS1 on lissencephaly progression remain elusive. Therefore, the present study employed combined bioinformatics and molecular modeling approach to identify potential damaging nsSNPs in the LIS1 gene and provide atomic insight into their roles in LIS1 loss of function. Using this approach, we identified three high-risk nsSNPs, including rs121434486 (F31S), rs587784254 (W55R), and rs757993270 (W55L) in the LIS1 gene, which are located on the N-terminal domain of LIS1. Molecular dynamics simulation highlighted that all variants decreased helical conformation, increased the intermonomeric distance, and thus disrupted intermonomeric contacts in the LIS1 dimer. Furthermore, the presence of variants also caused a loss of positive electrostatic potential and reduced dimer binding potential. Since self-dimerization is an essential aspect of LIS1 to recruit interacting partners, thus these variants are associated with the loss of LIS1 functions. As a corollary, these findings may further provide critical insights on the roles of LIS1 variants in brain malformation.

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

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

Computational Insights into the Deleterious Impacts of Missense Variants on N-Acetyl-d-glucosamine Kinase Structure and Function

An enzyme of the mammalian amino-sugar metabolism pathway, N-acetylglucosamine kinase (NAGK), that synthesizes N-acetylglucosamine (GlcNAc)-6-phosphate, is reported to promote dynein functions during mitosis, axonal and dendritic growth, cell migration, and selective autophagy, which all are unrelated to its enzyme activity. As non-enzymatic structural functions can be altered by genetic variation, we made an effort in this study aimed at deciphering the pathological effect of nonsynonymous single-nucleotide polymorphisms (nsSNPs) in NAGK gene. An integrated computational approach, including molecular dynamics (MD) simulation and protein-protein docking simulation, was used to identify the damaging nsSNPs and their detailed structural and functional consequences. The analysis revealed the four most damaging variants (G11R, G32R, G120E, and A156D), which are highly conserved and functional, positioned in both small (G11R and G32R) and large (G120E and A156D) domains of NAGK. G11R is located in the ATP binding region, while variants present in the large domain (G120E and A156D) were found to induce substantial alterations in the structural organizations of both domains, including the ATP and substrate binding sites. Furthermore, all variants were found to reduce binding energy between NAGK and dynein subunit DYNLRB1, as revealed by protein-protein docking and MM-GBSA binding energy calculation supporting their deleteriousness on non-canonical function. We hope these findings will direct future studies to gain more insight into the role of these variants in the loss of NAGK function and their role in neurodevelopmental disorders.

Molecular simulation studies to reveal the binding mechanisms of shikonin derivatives inhibiting VEGFR-2 kinase

Traditional vascular endothelial growth factor receptor 2 (VEGFR-2) inhibitors can manage angiogenesis; however, severe toxicity and resistance limit their long-term applications in clinical therapy. Shikonin (SHK) and its derivatives could be promising to inhibit the VEGFR-2 mediated angiogenesis, as they are reported to bind in the catalytic kinase domain with low affinity. However, the detailed molecular insights and binding dynamics of these natural inhibitors are unknown, which is crucial for potential SHK based lead design. Therefore, the present study employed molecular modeling and simulations techniques to get insight into the binding behaviors of SHK and its two derivates, β-hydroxyisovalerylshikonin (β-HIVS) and acetylshikonin (ACS). Here the intermolecular interactions between protein and ligands were studied by induced fit docking approach, which were further evaluated by treating QM/MM (quantum mechanics/molecular mechanics) and molecular dynamics (MD) simulation. The result showed that the naphthazarin ring of the SHK derivates is vital for strong binding to the catalytic domain; however, the binding stability can be modulated by the side chain modification. Because of having electrostatic potential, this ring makes essential interactions with the DFG (Asp1046 and Phe1047) motif and also allows interacting with the allosteric binding site. Taken together, the studies will advance our knowledge and scope for the development of new selective VEGFR-2 inhibitors based on SHK and its analogs.