In silico screening of deleterious single nucleotide polymorphisms (SNPs) and molecular dynamics simulation of disease associated mutations in gene responsible for oculocutaneous albinism type 6 (OCA 6) disorder

Impact of nsSNPs in human AIM2 and IFI16 gene: a comprehensive in silico analysis

AIM2 and IFI16 are the most studied members of AIM2-like receptors (ALRs) in humans and share a common N-Terminal PYD domain and C-terminal HIN domain. The HIN domain binds to dsDNA in response to the invasion of bacterial and viral DNA, and the PYD domain directs apoptosis-associated speck-like protein via protein-protein interactions. Hence, activation of AIM2 and IFI16 is crucial for protection against pathogenic assaults, and any genetic variation in these inflammasomes can dysregulate the human immune system. In this study, different computational tools were used to identify the most deleterious and disease-causing non-synonymous single nucleotide polymorphisms (nsSNPs) in AIM2 and IFI16 proteins. Molecular dynamic simulation was performed for the top damaging nsSNPs to study single amino acid substitution-induced structural alterations in AIM2 and IFI16. The observed results suggest that the variants G13V, C304R, G266R, and G266D for AIM2, and G13E and C356F are deleterious and affect structural integrity. We hope that the suggested deleterious nsSNPs and structural dynamics of AIM2 and IFI16 variants will guide future research to better understand the function of these variants with large-scale studies and may assist in fresher therapeutics focusing on these polymorphisms.Communicated by Ramaswamy H. Sarma.

Comprehensive analysis of non-synonymous missense SNPs of human galactose mutarotase (GALM) gene: an integrated computational approach

Missense Non-synonymous single nucleotide polymorphisms (nsSNPs) of Galactose Mutarotase (GALM) are associated with the Novel type of Galactosemia (Galactosemia type 4) together with symptoms such as high blood galactose levels and eye cataracts. The objective of the present study was to identify deleterious nsSNPs of GALM recorded on the dbSNP database through comprehensive insilico analysis. Among the 319 missense nsSNPs reported, various insilco tools predicted R78S, R82G, A163E, P210S, Y281C, E307G and F339C as the most deleterious mutations. Structural analysis, PTM analysis and molecular dynamics simulations (MDS) were carried out to understand the effect of these mutations on the structural and physicochemical properties of the GALM protein. The residues R82G and E307G were found to be part of the binding site that resulted in decreased surface accessibility. Replacing the charged wild-type residue with a neutral mutant type affected its substrate binding. All 7 mutations were found to increase the rigidity of the protein structure, which is unfavorable during ligand binding. The mutation F339E made the protein structure more rigid than all the other mutations. Y281 is a phosphorylated site, and therefore, less significant structural changes were observed when compared to other mutations; however, it may have significant differences in the usual functioning of the protein. In summary, the structural and functional analysis of missense SNPs of GALM is important to reduce the number of potential mutations to be evaluated in vitro to understand the association with some genetic diseases.Communicated by Ramaswamy H. Sarma.

Modified electrochemical aptasensor for ultrasensitive detection of tetracycline: In silico and in vitro studies

An ultrasensitive electrochemical aptasensor based on a glassy carbon electrode, modified by carbon nanofibers and carboxylated multi-walled carbon nanotubes was fabricated to detect tetracycline in food samples. The affinity of antibiotics, including kanamycin, tetracycline, ampicillin, and sulfadimethoxine toward desired sequences of aptamers and the stability of antibiotic-aptamer complexes were studied using molecular docking and molecular dynamic simulations. Moreover, the highest affinity and most stable complex were observed for tetracycline in complex with kanamycin-specific aptamer (KAP). Finally, KAP was used to develop an aptasensor. The central composite design (CCD) was used to optimize effective parameters. The biosensor achieved a wide dynamic linear range (1.0 × 10^-17-1.0 × 10^-5 M) and a low limit of detection (2.28 × 10^-18 M) under optimized conditions using differential pulse voltammetry. Using the developed aptasensor, tetracycline residues in milk samples were determined.

Regulation of K+-Dependent Na+/Ca2+-Exchangers (NCKX)

Potassium-dependent sodium-calcium exchangers (NCKX) have emerged as key determinants of calcium (Ca²⁺) signaling and homeostasis, especially in environments where ion concentrations undergo large changes, such as excitatory cells and transport epithelia. The regulation of NCKX transporters enables them to respond to the changing cellular environment thereby helping to shape the extent and kinetics of Ca²⁺ signals. This review examines the current knowledge of the different ways in which NCKX activity can be modulated. These include (i) cellular and dynamic subcellular location (ii); changes in protein expression mediated at the gene, transcript, or protein level (iii); genetic changes resulting in altered protein structure or expression (iv); regulation via changes in substrate concentration (v); and post-translational modification, partner protein interactions, and allosteric regulation. Detailed mechanistic understanding of NCKX regulation is an emerging area of research with the potential to provide important new insights into transporter function, the control of Ca²⁺ signals, and possible interventions for dysregulated Ca²⁺ homeostasis.

Identification of Natural Products as SENP2 Inhibitors for Targeted Therapy in Heart Failure

Aims: Sentrin-specific protease -2 (SENP2) is involved in deSUMOylation. Increased deSUMOylation in murine hearts by SENP2 upregulation resulted in cardiac dysfunction and congenital heart defects. Natural compounds via regulating cell proliferation and survival, induce cell cycle cessation, cell death, apoptosis, and producing reactive oxygen species and various enzyme systems cause disease prevention. Then, natural compounds can be suitable inhibitors and since SENP2 is a protein involved in heart disease, so our aim was inhibition of SENP2 by natural products for heart disease treatment. Material and methods: Molecular docking and molecular dynamics simulation of natural products i.e. Gallic acid (GA), Caffeic acid (CA), Thymoquinone (TQ), Betanin, Betanidin, Fisetin, and Ebselen were done to evaluate the SENP2 inhibitory effect of these natural products. The toxicity of compounds was also predicted. Results: The results showed that Betanin constituted a stable complex with SENP2 active site as it revealed low RMSD, high binding energy, and hydrogen bonds. Further, as compared to Ebselen, Betanin demonstrated low toxicity, formed a stable complex with SENP2 via four to seven hydrogen bonds, and constituted more stable MD plots. Therefore, depending upon the outcomes presented herein, Betanin significantly inhibited SENP2 and hence may be considered as a suitable natural compound for the treatment of heart failure. Further clinical trials must be conducted to validate its use as a potential SENP2 inhibitor.

PSnpBind: a database of mutated binding site protein-ligand complexes constructed using a multithreaded virtual screening workflow

A key concept in drug design is how natural variants, especially the ones occurring in the binding site of drug targets, affect the inter-individual drug response and efficacy by altering binding affinity. These effects have been studied on very limited and small datasets while, ideally, a large dataset of binding affinity changes due to binding site single-nucleotide polymorphisms (SNPs) is needed for evaluation. However, to the best of our knowledge, such a dataset does not exist. Thus, a reference dataset of ligands binding affinities to proteins with all their reported binding sites' variants was constructed using a molecular docking approach. Having a large database of protein-ligand complexes covering a wide range of binding pocket mutations and a large small molecules' landscape is of great importance for several types of studies. For example, developing machine learning algorithms to predict protein-ligand affinity or a SNP effect on it requires an extensive amount of data. In this work, we present PSnpBind: A large database of 0.6 million mutated binding site protein-ligand complexes constructed using a multithreaded virtual screening workflow. It provides a web interface to explore and visualize the protein-ligand complexes and a REST API to programmatically access the different aspects of the database contents. PSnpBind is open source and freely available at https://psnpbind.org .

In silico assessment of DNA damage response gene variants associated with head and neck cancer

Head and neck cancer (HNC), the sixth most common cancer globally, stands first in India, especially Northeast India, where tobacco usage is predominant, which introduces various carcinogens leading to malignancies by accumulating DNA damages. Consequently, the present work aimed to predict the impact of significant germline variants in DNA repair and Tumour Suppressor genes on HNC development. WES in Ion Proton^TM platform on 'discovery set' (n = 15), followed by recurrence assessment of the observed variants on 'confirmation set' (n = 40) using Sanger Sequencing was performed on the HNC-prevalent NE Indian populations. Initially, 53 variants were identified, of which seven HNC-linked DNA damage response gene variants were frequent in the studied populations. Different tools ascertained the biological consequences of these variants, of which the non-coding variants viz. EXO1_rs4150018, RAD52_rs6413436, CHD5_rs2746066, HACE1_rs6918700 showed risk, while FLT3_rs2491227 and BMPR1A_rs7074064 conferred protection against HNC by affecting transcriptional regulation and splicing mechanism. Molecular Dynamics Simulation of the full-length p53 model predicted that the observed coding TP53_rs1042522 variant conferred HNC-risk by altering the structural dynamics of the protein, which displayed difficulty in the transition between active and inactive conformations due to high-energy barrier. Subsequent pathway and gene ontology analysis revealed that EXO1, RAD52 and TP53 variants affected the Double-Strand Break Repair pathway, whereas CHD5 and HACE1 variants inactivated DNA repair cascade, facilitating uncontrolled cell proliferation, impaired apoptosis and malignant transformation. Conversely, FLT3 and BMPR1A variants protected against HNC by controlling tumorigenesis, which requires experimental validation. These findings may serve as prognostic markers for developing preventive measures against HNC.

Key criteria for engineering mycotoxin binding aptamers via computational simulations: Aflatoxin B1 as a case study

Due to the difficulties in monoclonal antibody production specific to mycotoxins, aptameric probes have been considered as suitable alternatives. The low efficiency of the SELEX procedure in screening high affinity aptamers for binding mycotoxins as small molecules can be significantly improved through computational techniques. Previously, we designed five new aptamers to aflatoxin B₁ (AFB1) based on a known aptamer sequence (Patent: PCT/CA2010/001 292, Apt1) through a genetic algorithm-based in silico maturation strategy and experimentally measured their affinity to the target toxin. Here, integrated molecular dynamic simulation (MDs) studies with molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) analysis to clarify the binding modes, critical interacting nucleic bases and energy component contributions in the six AFB1-binding aptamers. The aptamer F20, which was selected in the first work, showed the best free binding energy and complex stability compared to other aptamers. The trajectory analysis revealed that AFB1 recognized F20 through the groove binding mode along with precise shape complementarity. The MD simulation results revealed that dynamic water intermediate interactions also play a key role in promoting complex stability. According to the MM-PBSA calculations, van der Waals contacts were identified as dominant energy components in all complexes. Interestingly, a high consistency is observed between the experimentally obtained binding affinities of the six aptamers with their free energy solvation. The computational findings, confirmed via previous experiments, highlighted the binding modes, the dynamic hydration of complex components and the total free interacting energy as the crucial criteria in discovering high functional aptameric probes.

Computational Analysis of Gly482Ser Single-Nucleotide Polymorphism in PPARGC1A Gene Associated with CAD, NAFLD, T2DM, Obesity, Hypertension, and Metabolic Diseases

Peroxisome proliferator-activated receptor-gamma coactivator 1-alpha (PPARGC1A) regulates the expression of energy metabolism's genes and mitochondrial biogenesis. The essential roles of PPARGC1A encouraged the researchers to assess the relation between metabolism-related diseases and its variants. To study Gly482Ser (+1564G/A) single-nucleotide polymorphism (SNP) after PPARGC1A modeling, we substitute Gly482 for Ser482. Stability prediction tools showed that this substitution decreases the stability of PPARGC1A or has a destabilizing effect on this protein. We then utilized molecular dynamics simulation of both the Gly482Ser variant and wild type of the PPARGC1A protein to analyze the structural changes and to reveal the conformational flexibility of the PPARGC1A protein. We observed loss flexibility in the RMSD plot of the Gly482Ser variant, which was further supported by a decrease in the SASA value in the Gly482Ser variant structure of PPARGC1A and an increase of H-bond with the increase of β-sheet and coil and decrease of turn in the DSSP plot of the Gly482Ser variant. Such alterations may significantly impact the structural conformation of the PPARGC1A protein, and it might also affect its function. It showed that the Gly482Ser variant affects the PPARGC1A structure and makes the backbone less flexible to move. In general, molecular dynamics simulation (MDS) showed more flexibility in the native PPARGC1A structure. Essential dynamics (ED) also revealed that the range of eigenvectors in the conformational space has lower extension of motion in the Gly482Ser variant compared with WT. The Gly482Ser variant also disrupts PPARGC1A interaction. Due to this single-nucleotide polymorphism in PPARGC1A, it became more rigid and might disarray the structural conformation and catalytic function of the protein and might also induce type 2 diabetes mellitus (T2DM), coronary artery disease (CAD), and nonalcoholic fatty liver disease (NAFLD). The results obtained from this study will assist wet lab research in expanding potent treatment on T2DM.

Melanins as Sustainable Resources for Advanced Biotechnological Applications

Melanins are a class of biopolymers that are widespread in nature and have diverse origins, chemical compositions, and functions. Their chemical, electrical, optical, and paramagnetic properties offer opportunities for applications in materials science, particularly for medical and technical uses. This review focuses on the application of analytical techniques to study melanins in multidisciplinary contexts with a view to their use as sustainable resources for advanced biotechnological applications, and how these may facilitate the achievement of the United Nations Sustainable Development Goals.

Bioinformatic and computational analysis for predominant mutations of the Nrf2/Keap1 complex in pediatric leukemia

The levels of reactive oxygen species (ROS) are tightly controlled and regulated by Nuclear Factor Erythroid-2-Like 2 (Nrf2) transcription factor, which is the main regulator of antioxidant responses and its suppressor protein Kelch-like ECH-associated protein 1 (Keap1). Our previous study has identified six novel changes in Nrf2/Keap1 pathway in pediatric ALL, which were described for the first time. These changes in the pathway are likely to alter the evolutionary process of amino acids and cause structural changes in the final products of genes. In this study, we aimed to compare the pathogenicity of eight determined mutations reported in our previous study by utilizing different programs with different algorithms and molecular dynamics simulation. Since it is too difficult to handle each existing mutation in a wet laboratory, in silico methods may give suggestion to choose the important mutations for further analysis and to establish the appropriate patient population and conduct wet laboratory studies. For this purpose, four different algorithms were used to evaluate the effects of single amino acid mutation. In addition, root-mean-square deviation, root-mean-square fluctuation and free-energy landscape analyses were performed to observe stability, flexibility and energetically favorable conformations, respectively, for each amino acid mutation. As a result, our study emphasizes the importance of Keap1 mutations in pediatric ALL Nrf2/Keap1 pathway, a total of eight mutations, two of which were shown for the first time in our study. Especially the mutations in the Keap1 Broad-Complex, Tramtrack and Bric-à-brac domain are worthy of attention.Communicated by Ramaswamy H. Sarma.

Computational Screening and Analysis of Lung Cancer Related Non-Synonymous Single Nucleotide Polymorphisms on the Human Kirsten Rat Sarcoma Gene

The human KRAS (Kirsten rat sarcoma) is an oncogene, involved in the regulation of cell growth and division. The mutations in the KRAS gene have the potential to cause normal cells to become cancerous in human lungs. In the present study, we focus on non-synonymous single nucleotide polymorphisms (nsSNPs), which are point mutations in the DNA sequence leading to the amino acid variants in the encoded protein. To begin with, we developed a pipeline to utilize a set of computational tools in order to obtain the most deleterious nsSNPs (Q22K, Q61P, and Q61R) associated with lung cancer in the human KRAS gene. Furthermore, molecular dynamics simulation and structural analyses of the 3D structures of native and mutant proteins confirmed the impact of these nsSNPs on the stability of the protein. Finally, the experimental results demonstrated that the structural stability of the mutant proteins was worse than that of the native protein. This study provides significant guidance for narrowing down the number of KRAS mutations to be screened as potential diagnostic biomarkers and to better understand the structural and functional mechanisms of the KRAS protein.

Discover potential inhibitors for PFKFB3 using 3D-QSAR, virtual screening, molecular docking and molecular dynamics simulation

The 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase-3 (PFKFB3) is a master regulator of glycolysis in cancer cells by synthesizing fructose-2,6-bisphosphate (F-2,6-BP), a potent allosteric activator of phosphofructokinase-1 (PFK-1), which is a rate-limiting enzyme of glycolysis. PFKFB3 is an attractive target for cancer treatment. It is valuable to discover promising inhibitors by using 3D-QSAR pharmacophore modeling, virtual screening, molecular docking and molecular dynamics simulation. Twenty molecules with known activity were used to build 3D-QSAR pharmacophore models. The best pharmacophore model was ADHR called Hypo1, which had the highest correlation value of 0.98 and the lowest RMSD of 0.82. Then, the Hypo1 was validated by cost value method, test set method and decoy set validation method. Next, the Hypo1 combined with Lipinski's rule of five and ADMET properties were employed to screen databases including Asinex and Specs, total of 1,048,159 molecules. The hits retrieved from screening were docked into protein by different procedures including HTVS, SP and XP. Finally, nine molecules were picked out as potential PFKFB3 inhibitors. The stability of PFKFB3-lead complexes was verified by 40 ns molecular dynamics simulation. The binding free energy and the energy contribution of per residue to the binding energy were calculated by MM-PBSA based on molecular dynamics simulation.