In silico prediction of a disease-associated STIL mutant and its affect on the recruitment of centromere protein J (CENPJ)

Exosomes derived from human umbilical cord mesenchymal stem cells decrease neuroinflammation and facilitate the restoration of nerve function in rats suffering from intracerebral hemorrhage

Exosomes derived from human umbilical cord mesenchymal stem cells (hUCMSC-ex) have become a hopeful substitute for whole-cell therapy due to their minimal immunogenicity and tumorigenicity. The present study aimed to investigate the hypothesis that hUCMSC-ex can alleviate excessive inflammation resulting from intracerebral hemorrhage (ICH) and facilitate the rehabilitation of the nervous system in rats. In vivo, hemorrhagic stroke was induced by injecting collagenase IV into the striatum of rats using stereotactic techniques. hUCMSC-ex were injected via the tail vein at 6 h after ICH model establishment at a dosage of 200 µg. In vitro, astrocytes were pretreated with hUCMSC-ex and then stimulated with hemin (20 μmol/mL) to establish an ICH cell model. The expression of TLR4/NF-κB signaling pathway proteins and inflammatory factors, including TNF-α, IL-1β, and IL-10, was assessed both in vivo and in vitro to investigate the impact of hUCMSC-ex on inflammation. The neurological function of the ICH rats was evaluated using the corner turn test, forelimb placement test, Longa score, and Bederson score on the 1st, 3rd, and 5th day. Additionally, RT-PCR was employed to examine the mRNA expression of TLR4 following hUCMSC-ex treatment. The findings demonstrated that hUCMSC-ex downregulated the protein expression of TLR4, NF-κB/P65, and p-P65, reduced the levels of pro-inflammatory cytokines TNF-α and IL-1β, and increased the expression of the anti-inflammatory cytokine IL-10. Ultimately, the administration of hUCMSC-ex improved the behavioral performance of the ICH rats. However, the results of PT-PCR indicated that hUCMSC-ex did not affect the expression of TLR4 mRNA induced by ICH, suggesting that hUCMSCs-ex may inhibit TLR4 translation rather than transcription, thereby suppressing the TLR4/NF-κB signaling pathway. We can conclude that hUCMSC-ex mitigates hyperinflammation following ICH by inhibiting the TLR4/NF-κB signaling pathway. This study provides preclinical evidence for the potential future application of hUCMSC-ex in the treatment of cerebral injury.

Discovery of Palbociclib as a potent c-Myc G4 stabilizer for lung cancer treatment using molecular docking, molecular dynamics simulation, and in vitro activity evaluation

Despite significant progress in lung cancer treatment, this disease remains a prevalent and serious global malignancy, leading to high rates of illness and death. Urgent research is needed to discover new or alternative therapies that can improve clinical outcomes for lung cancer patients. In our study, we successfully demonstrated the effectiveness of Palbociclib, a CDK4/6 inhibitor, in suppressing the growth of lung cancer cells. The IC50 values obtained were 11.00 μM and 11.74 μM for H1299 and A549 cells, respectively. Furthermore, our findings indicate that Palbociclib may possess strong c-Myc G4 stabilizing properties by significantly reducing both protein and mRNA expression levels of c-Myc. Additionally, Palbociclib induces apoptosis and causes cell cycle arrest at the G2/M phase in two cells. Through circular dichroism (CD), molecular docking, and molecular dynamics (MD) simulation, we have provided evidence that Palbociclib enhances the structural stability of c-Myc G4 while exhibiting a high binding affinity to its ligand's binding site on c-Myc G4. These results suggest that Palbociclib holds promise as a novel c-Myc G4 stabilizer for treating cancers associated with abnormal c-Myc activity; further optimization and development are warranted.

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.

In silico identification of potential protein kinase C alpha inhibitors from phytochemicals from IMPPAT database for anticancer therapeutics: a virtual screening approach

Protein Kinase C alpha (PKCα) is a critical signaling molecule that plays a crucial role in various physiological processes, including cell growth, differentiation, and survival. Over the years, there has been a growing interest in targeting PKCα as a promising drug target for the treatment of various diseases, including cancer. Targeting PKCα can, therefore, serve as a potential strategy to prevent cancer progression and enhance the efficacy of conventional anticancer therapies. We conducted a systematic search for promising compounds for their anticancer potential that target PKCα using natural compounds from the IMPPAT database. The initial compounds were screened through various tests, including analysis of their physical and chemical properties, PAINS filter, ADMET analysis, PASS analysis, and specific interaction analysis. We selected those that showed high binding affinity and specificity to PKCα from the screened compounds, and we further analyzed them using molecular dynamics simulations (MDS) and principal component analysis (PCA). Various systematic parameters from the MDS analyses suggested that the protein-ligand complexes were stabilized throughout the simulation trajectories of 100 nanoseconds (ns). Our findings indicated that compounds Nicandrenone and Withaphysalin D bind to PKCα with high stability and affinity, making them potential candidates for further research in cancer therapeutics innovation in clinical contexts.Communicated by Ramaswamy H. Sarma.

Atractylodes macrocephala Koidz. and Cuscuta chinensis Lam. extract relieves insulin resistance via PI3K/Akt signalling in diabetic Drosophila

Background and aim

Type-2 diabetes mellitus (T2DM) is mainly characterized by insulin resistance (IR) induced by hyperglycaemia and insufficient insulin secretion. We employed a diabetic fly model to examine the effect and molecular mechanism of Atractylodes macrocephala Koidz. and Cuscuta chinensis Lam. (AMK-CCL) extract as traditional Chinese medicine in treating IR and T2DM.

Experimental procedure

The contents of the active ingredients (rhamnose, xylose, mannose, and hyperoside) in AMK-CCL extract were determined by high-performance liquid chromatography. Wild-type (Cg-GAL4/+) or diabetic (Cg > InRK1409A) Drosophila flies were divided into the control group or metformin group and AMK-CCL (0.0125, 0.025, 0.05, 0.1 g/ml) groups. Food intake, haemolymph glucose and trehalose, protein, weight, triglycerides (TAG), and glycogen were measured to assess glycolipid metabolism. Phosphatidylinositol-3-kinase (PI3K)/Akt signalling was detected using fluorescent reporters [tGPH, Drosophila forkhead box O (dFoxO)-green fluorescent protein (GFP), Glut1-GFP, 2-NBDG] in vivo. Glut1/3 mRNA levels and Akt phosphorylation levels were detected by quantitative polymerase chain reaction and western blotting, respectively, in vitro.

Results

AMK-CCL extract contained 0.038 % rhamnose, 0.017 % xylose, 0.69 % mannose, and 0.039 % hyperoside. AMK-CCL at 0.0125 g/mL significantly suppressed the increase in circulating glucose, and the decrease in body weight, TAG, and glycogen contents of diabetic flies. AMK-CCL improved PI3K activity, Akt phosphorylation, Glut1/3 expression, and glucose uptake in diabetic flies, and also rescued diabetes-induced dFoxO nuclear localisation.

Conclusions

These findings indicate that AMK-CCL extract ameliorates IR-induced diabetes via the PI3K/Akt signalling pathway, providing an experimental basis for clinical treatment.

Acute and sub-chronic toxicity study of novel polyherbal formulation in non-alcoholic fatty liver using Wistar rats

Aim: This study assessed the acute and sub-chronic toxicity of a novel polyherbal formulation tablet in Wistar rats Materials & methods: Acute toxicity and sub-chronic toxicity was assessed following OECD (Organisation for the Economic Co-operation and Development) guidelines based on 423 and 408. Results & conclusion: No mortality and toxicity showed in rats during acute toxicity. The LD50 of the extract was at 2000 mg/kg. In sub-chronic study, both sex rats were orally administered at 250, 500,1000 and 2000 mg/kg for 90 days and revealed no significant difference (p < 0.05) in hematological and other parameters compared with the control. Histopathology study did not reveal morphological alteration. The No observed adverse effect level of the tablet was observed until 2000 mg/kg.

Discovery of potential TAAR1 agonist targeting neurological and psychiatric disorders: An in silico approach

Trace amine-associated receptor 1 (TAAR1) is a G-protein-coupled receptor which is primarily expressed in the brain. It is activated by trace amines which play a role in regulating neurotransmitters like dopamine, serotonin and norepinephrine. TAAR1 agonists have potential applications in the treatment of neurological and psychiatric disorders, especially schizophrenia. In this study, we have used a structure-based virtual screening approach to identify potential TAAR1 agonist(s). We have modelled the structure of TAAR1 and predicted the binding pocket. Further, molecular docking of a few well-known antipsychotic drugs was carried out with TAAR1 model, which showed key interactions with the binding pocket. From screening a library of 5 million compounds from the Enamine REAL Database using structure-based virtual screening method, we shortlisted 12 compounds which showed good docking score, glide energy and interactions with the key residues. One lead compound (Z31378290) was finally selected. The lead compound showed promising binding affinity and stable interactions with TAAR1 during molecular dynamics simulations and demonstrated better van der Waals and binding energy than the known agonist, ulotaront. Our findings suggest that the lead compound may serve as a potential TAAR1 agonist, offering a promising avenue for the development of new therapies for neurological and psychiatric disorders.

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.

Structural insights into the pathogenicity of point mutations in human acyl-CoA dehydrogenase homotetramers

Acyl-CoA dehydrogenase deficiency (ACAD) is an inherited and potentially fatal disorder with variable clinical symptoms. The relationship between pathogenicity and deleterious point mutations is investigated here in ACAD structures of short (SCAD) and medium-chain (MCAD) types. Structures and dynamic features of native and mutant forms of enzymes models were compared. A total of 2.88 µs molecular dynamics simulations were performed at four different temperatures. Total energy, RMSD, protein ligand interactions and affinity, RMSF measures, secondary structure changes, and important interactions were studied. Mutations in the three main domains of ACADs are pathogenic, while those located at linker turns are not. Mutations affect mostly tetramer formations, secondary structures, and many contacts and interactions. In R206H (MCAD mutant) which is experimentally known to cause a huge turnover decrease, the lack of a single H-bond between substrate and FAD was observed. Secondary structures showed temperature-dependent changes, and SCAD activity was found to be highly correlated to the enzyme helix 3-10 content. Finally, RMSF patterns pointed to one important loop that maintains the substrate close to the active site and is a cause of substrate wobbling upon mutation. Despite similar structure, function, and cellular location, SCAD and MCAD may have different optimum temperatures that are related to the structure taken at that specific temperature. In conclusion, new insight has been provided on the effect of various SCAD and MCAD pathogenic mutations on the structure and dynamical features of the enzymes.

Genetic analysis of albinism caused by compound heterozygous mutations of the OCA2 gene in a Chinese family

BACKGROUND

Oculocutaneous albinism (OCA) is a group of rare genetic disorders characterized by a reduced or complete lack of melanin in the skin, hair, and eyes. Patients present with colorless retina, pale pink iris, and pupil, and fear of light. The skin, eyebrows, hair, and other body hair are white or yellowish-white. These conditions are caused by mutations in specific genes necessary for the production of melanin. OCA is divided into eight clinical types (OCA1-8), each with different clinical phenotypes and potential genetic factors. This study aimed to identify the genetic causes of non-syndromic OCA in a Chinese Han family.

METHODS

We performed a comprehensive clinical examination of family members, screened for mutation loci using whole exome sequencing (WES) technology, and predicted mutations using In silico tools.

RESULTS

The patient's clinical manifestations were white skin, yellow hair, a few freckles on the cheeks and bridge of the nose, decreased vision, blue iris, poorly defined optic disk borders, pigmentation of the fundus being insufficient, and significant vascular exposure. The WES test results indicate that the patient has compound heterozygous mutations in the OCA2 gene (c.1258G > A (p.G420R), c.1441G > A (p.A481T), and c.2267-2 A > C), respectively, originating from her parents. Among them, c.1258G > A (p.G420R) is a de novo mutation with pathogenic. Our analysis suggests that compound heterozygous mutations in the OCA2 gene are the primary cause of the disease in this patient.

CONCLUSIONS

The widespread application of next-generation sequencing technologies such as WES in clinical practice can effectively replace conventional detection methods and assist in the diagnosis of clinical diseases more quickly and accurately. The newly discovered c.1258G > A (p.G420R) mutation can update and expand the gene mutation spectrum of OCA2-type albinism.

Genes adaptability and NOL6 protein inhibition studies of fabricated flavan-3-ols lead skeleton intended to treat breast carcinoma

Breast cancer invasive 2.3 million women worldly and second prominent factor of cancer-related mortality. Finding a new site-specific and safe small molecule is a current need in this field. With the aid of deep learning Algorithms, we analyzed the published big database from cancer CBioportal to find the best target protein. Further, Multi-omics analysis such as enrichment analysis, scores of molecular, RNA biological function at a cellular level, and protein domain were obtained and matched to find the better hit molecules. The gene analysis output shows nucleolar protein 6 plays a significant responsibility in breast carcinoma and 354 natural and synthetic lead molecules are docked inside the active site. Docking result gave the output hit molecule falavan-3-ols with a binding score of -5.325 (Kcal/mol) and interaction analysis illustrates, 13 active amino acids favoring the binding interaction with functional groups of the hit molecule compared to the standard molecule Abemacilib (-2.857 (Kcal/mol)). Best docked complex of flavan-3-ols and NOL6 protein subjected to dynamic simulation 100 ns to study the stability. The results proved that π-π stacked, carbon‑hydrogen and electrostatic interactions are stable throughout the 100 ns simulation. The overall results conclude the hit molecule flavan-3-ol will be a safe and potent lead molecule to generate and treat breast carcinoma patients.

Cuproptosis gene-related, neural network-based prognosis prediction and drug-target prediction for KIRC

BACKGROUND

Kidney renal clear cell carcinoma (KIRC), as a common case in renal cell carcinoma (RCC), has the risk of postoperative recurrence, thus its prognosis is poor and its prognostic markers are usually based on imaging methods, which have the problem of low specificity. In addition, cuproptosis, as a novel mode of cell death, has been used as a biomarker to predict disease in many cancers in recent years, which also provides an important basis for prognostic prediction in KIRC. For postoperative patients with KIRC, an important means of preventing disease recurrence is pharmacological treatment, and thus matching the appropriate drug to the specific patient's target is also particularly important. With the development of neural networks, their predictive performance in the field of medical big data has surpassed that of traditional methods, and this also applies to the field of prognosis prediction and drug-target prediction.

OBJECTIVE

The purpose of this study is to screen for cuproptosis genes related to the prognosis of KIRC and to establish a deep neural network (DNN) model for patient risk prediction, while also developing a personalized nomogram model for predicting patient survival. In addition, sensitivity drugs for KIRC were screened, and a graph neural network (GNN) model was established to predict the targets of the drugs, in order to discover potential drug action sites and provide new treatment ideas for KIRC.

METHODS

We used the Cancer Genome Atlas (TCGA) database, International Cancer Genome Consortium (ICGC) database, and DrugBank database for our study. Differentially expressed genes (DEGs) were screened using TCGA data, and then a DNN-based risk prediction model was built and validated using ICGC data. Subsequently, the differences between high- and low-risk groups were analyzed and KIRC-sensitive drugs were screened, and finally a GNN model was trained using DrugBank data to predict the relevant targets of these drugs.

RESULTS

A prognostic model was built by screening 10 significantly different cuproptosis-related genes, the model had an AUC of 0.739 on the training set (TCGA data) and an AUC of 0.707 on the validation set (ICGC data), which demonstrated a good predictive performance. Based on the prognostic model in this paper, patients were also classified into high- and low-risk groups, and functional analyses were performed. In addition, 251 drugs were screened for sensitivity, and four drugs were ultimately found to have high sensitivity, with 5-Fluorouracil having the best inhibitory effect, and subsequently their corresponding targets were also predicted by GraphSAGE, with the most prominent targets including Cytochrome P450 2D6, UDP-glucuronosyltransferase 1A, and Proto-oncogene tyrosine-protein kinase receptor Ret. Notably, the average accuracy of GraphSAGE was 0.817 ± 0.013, which was higher than that of GAT and GTN.

CONCLUSION

Our KIRC risk prediction model, constructed using 10 cuproptosis-related genes, had good independent prognostic ability. In addition, we screened four highly sensitive drugs and predicted relevant targets for these four drugs that might treat KIRC. Finally, literature research revealed that four drug-target interactions have been demonstrated in previous studies and the remaining targets are potential sites of drug action for future research.

Structure-based functional analysis of a novel NADPH-producing glyceraldehyde-3-phosphate dehydrogenase from Corynebacterium glutamicum

Corynebacterium glutamicum is an industrial workhorse applied in the production of valuable biochemicals. In the process of bio-based chemical production, improving cofactor recycling and mitigating cofactor imbalance are considered major solutions for enhancing the production yield and efficiency. Although, glyceraldehyde-3-phosphate dehydrogenase (GapDH), a glycolytic enzyme, can be a promising candidate for a sufficient NADPH cofactor supply, however, most microorganisms have only NAD-dependent GapDHs. In this study, we performed functional characterization and structure determination of novel NADPH-producing GapDH from C. glutamicum (CgGapX). Based on the crystal structure of CgGapX in complex with NADP cofactor, the unique structural features of CgGapX for NADP stabilization were elucidated. Also, N-terminal additional region (Auxiliary domain, AD) appears to have an effect on enzyme stabilization. In addition, through structure-guided enzyme engineering, we developed a CgGapX variant that exhibited 4.3-fold higher kcat, and 1.2-fold higher kcat/KM values when compared with wild-type. Furthermore, a bioinformatic analysis of 100 GapX-like enzymes from 97 microorganisms in the KEGG database revealed that the GapX-like enzymes possess a variety of AD, which seem to determine enzyme stability. Our findings are expected to provide valuable information for supplying NADPH cofactor pools in bio-based value-added chemical production.

In-silico bioprospecting of secondary metabolites from endophytic Streptomyces spp. against Magnaporthe oryzae, a cereal killer fungus

Rice blast disease, caused by Magnaporthe oryzae, is the most devastating cereal killer worldwide. Note that melanin pigment is an essential factor of M. oryzae virulence, thus fungicides interfering with melanin biosynthesizing enzymes would reduce the pathogenicity. Scytalone dehydratase (SDH) is the key target for commercial fungicides, like carpropamid, due to its role in the dehydration reaction of the fungal melanin pathway. However, a single-point mutation (V75M) in SDH elicits resistance to carpropamid. A lack of effective fungicides against this resistant strain expedited the quest for novel bioactive inhibitors. Currently, bacterial endophytes like Streptomyces have been heralded for synthesizing bioactive metabolites to protect plants from phytopathogens. The literature search led to the identification of 21 Streptomyces spp. symbionts of paddy that can suppress M. oryzae growth. An antiSMASH server was used to explore Streptomyces spp. gene clusters and found 4463 putative metabolites. Besides, 745 unique metabolites were subjected to a series of virtual screening techniques. Ideally, this process identified five potential SDH inhibitors. The docking result highlights that the metabolite pseudopyronine A interacted hydrophobically with both Val75 of SDHWT and Met75 of SDHV75M targets. Moreover, pseudopyronine A has a higher binding free energy with SDHWT (- 89.94 kcal/mol) and SDHV75M (- 71.95 kcal/mol). Interestingly, the pyranones scaffold of pseudopyronine A was reported for antifungal activity against phytopathogens. Dynamic behavior confirms that pseudopyronine A has excellent conformational states with both SDHWT and SDHV75M. Altogether, we hope that this study creates a new avenue for the discovery of novel phytopathogen inhibitors from endophytes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03859-7.

In Vitro Antioxidant and In Vivo Hepatoprotective Properties of Wissadula periplocifolia Extract

Wissadula periplocifolia (L.) Thwaites is a traditional medicinal plant belonging to the family Malvaceae, used in folk medicine for inflamed snake bites and bee stings. The current study was designed to investigate the in vitro antioxidant and in vivo anti-inflammatory and hepatoprotective activities of 80% ethanol extract of W. periplocifolia and its different fractions. The crude ethanolic extract (CEE) was then serially fractionated with petroleum ether fraction (PEF), chloroform fraction (CHF), and aqueous fraction (AQF). The antioxidant activity was assessed using 1,1-diphenyl-2-picrylhydrazyl (DPPH) radical assay, anti-inflammatory activity was determined in the xylene-induced ear edema model, and hepatoprotective activity was measured in the paracetamol-induced hepatic injury model. PEF showed a significant scavenging effect with an IC50 value of 33.5 μg/mL, followed by CEE (IC50 = 42.2 μg/mL), CHF (IC50 = 77 μg/mL), and AQF (IC50 = 80 μg/mL), compared to standard butylated hydroxytoluene (IC50 = 14.8 μg/mL). Both doses of CEE (250 and 500 mg/kg) could reduce ear edema by 41.3 and 50%, respectively, compared to standard diclofenac sodium (76.09%). Moreover, CEE significantly reduces the elevated liver enzymes (ALT, AST, and ALP), compared to control. Nevertheless, it elevated blood protein and reduced the blood bilirubin level (p < 0.01), compared to control. Histopathological studies also indicated significant protection of the liver from paracetamol-induced liver damage. In conclusion, W. periplocifolia could be a good source of antioxidant and hepatoprotective phytochemicals; meanwhile, toxicological and pharmacokinetic studies are recommended.

Computational modeling of in-stent restenosis: Pharmacokinetic and pharmacodynamic evaluation

Persistence of the pathology of in-stent restenosis even with the advent of drug-eluting stents warrants the development of highly resolved in silico models. These computational models assist in gaining insights into the transient biochemical and cellular mechanisms involved and thereby optimize the stent implantation parameters. Within this work, an already established fully-coupled Lagrangian finite element framework for modeling the restenotic growth is enhanced with the incorporation of endothelium-mediated effects and pharmacological influences of rapamycin-based drugs embedded in the polymeric layers of the current generation drug-eluting stents. The continuum mechanical description of growth is further justified in the context of thermodynamic consistency. Qualitative inferences are drawn from the model developed herein regarding the efficacy of the level of drug embedment within the struts as well as the release profiles adopted. The framework is then intended to serve as a tool for clinicians to tune the interventional procedures patient-specifically.

Decoding the deactivation mechanism of R192W mutation of ZAP-70 using molecular dynamics simulations and binding free energy calculations

CONTEXT

ZAP-70 (zeta-chain-associated protein of 70 kDa), serving as a critical regulator for T cell antigen receptor signaling, represents an attractive therapeutic target for autoimmunity disease. How the mechanistical mechanism of ZAP-70 to a human autoimmune syndrome-associated R192W mutation remains unclear. The results indicated that the R192W mutation of ZAP-70 clearly affected the conformational flexibility of the N-terminal ITAM-Y2P. Structural analysis unveiled that the R192W mutation of ZAP-70 caused the exposure of the N-terminal ITAM-Y2P to the solvent. MM-GBSA binding free energy calculations exhibited that the R192W mutation decreased the binding affinity of ITAM-Y2P to the ZAP-70 mutant. Residue-based free energy decomposition further revealed that the protein-peptide interaction networks involving electrostatic interactions provide significant contributions for complex formation. The energy unfavorable residues include Arg43, Arg192, Tyr240, and Lys244 from ZAP-70 and Asn301, Leu303, pY304, and pY315 from ITAM-Y2P in the R192W mutant. Our obtained results may help the understanding of the deactivation mechanism of ZAP-70 induced by the R192W mutation.

METHODS

In the work, multiple replica molecular dynamics simulations and molecular mechanics-generalized Born surface area (MM-GBSA) method were performed to reveal the doubly phosphorylated ITAMs (ITAM-Y2P)-mediated deactivation mechanism of ZAP-70 induced by the R192W mutation.

Free energy landscape and thermodynamics properties of novel mutations in PncA of pyrazinamide resistance isolates of Mycobacterium tuberculosis

Pyrazinamide (PZA) is one of the first-line antituberculosis therapy, active against non-replicating Mycobacterium tuberculosis (Mtb). The conversion of PZA into pyrazinoic acid (POA), the active form, required the activity of pncA gene product pyrazinamidase (PZase) activity. Mutations occurred in pncA are the primary cause behind the PZA resistance. However, the resistance mechanism is important to explore using high throughput computational approaches. Here we aimed to explore the mechanism of PZA resistance behind novel P62T, L120R, and V130M mutations in PZase using 200 ns molecular dynamics (MD) simulations. MD simulations were performed to observe the structural changes for these three mutants (MTs) compared to the wild types (WT). Root means square fluctuation, the radius of gyration, free energy landscape, root means square deviation, dynamic cross-correlation motion, and pocket volume were found in variation between WT and MTs, revealing the effects of P62T, L120R, and V130M. The free energy conformational landscape of MTs differs significantly from the WT system, lowering the binding of PZA. The geometric shape complementarity of the drug (PZA) and target protein (PZase) further confirmed that P62T, L120R, and V130M affect the protein structure. These effects on PZase may cause vulnerability to convert PZA into POA.Communicated by Ramaswamy H. Sarma.

Design of siRNA Bioconjugates for Efficient Control of Cancer-Associated Membrane Receptors

Research on siRNA delivery has seen tremendous growth over the past few decades. As one of the major delivery strategies, siRNA bioconjugates offer the potential to enhance and extend the pharmacological properties of siRNAs while minimizing toxicity. In this paper, we suggest the development of a siRNA conjugate platform with peptides and proteins that are ligands of target receptors for cancer treatment. The siRNA bioconjugates target and block the receptor membrane proteins, enter the cells through receptor-mediated endocytosis, and inhibit the expression of that same target membrane receptor, thereby doubly controlling the function of the membrane proteins. The three kinds of bioconjugates targeting CD47, PD-L1, and EGFR were synthesized via two different copper-free click chemistry reactions. Results showed the cellular uptake of each conjugate, reduction of target gene expression, and efficient functional control of receptor proteins. This platform provides an effective approach for regulating membrane proteins in various diseases beyond cancer.

Active site determination of novel plant versatile peroxidase extracted from Citrus sinensis and bioconversion of β-naphthol

A ligninolytic peroxidase called versatile peroxidase, VP, (EC 1.11.1.16) is an iron-containing metalloenzyme. The most distinctive feature of this enzyme is its composite molecular framework, which combines lignin peroxidase's capacity to oxidize compounds with high-redox potential with manganese peroxidase's capacity to oxidize Mn2+ to Mn3+. In this study, we have extracted amino acid sequences from the Citrus sinensis source and subjected them to various computation tools to visualize the insight secondary and 3D structure, physicochemical properties, and validation of the structure which have not been studied so far to further investigate the catalytic efficiency and effectiveness of VP. The binding energies of HEME and HEME C (HEC) ligands with produced PDB (6rqf.1. A) have been also assessed, analyzed, and confirmed utilizing AutoDock. Binding energies were calculated using the AutoDock and validated by MD simulation using SCHRODINGER DESMOND. Most stable confirmation was achieved through a protein-ligand interaction study. Bio-technological use of VP in the biotransformation of β-naphthol has also been studied. The findings in the current study will have a substantial impact on proteomics, biochemistry, biotechnology, and possible uses of versatile peroxidase in the bio-remediation of different toxic organic compounds.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03758-x.

3D based generative PROTAC linker design with reinforcement learning

Proteolysis targeting chimera (PROTAC), has emerged as an effective modality to selectively degrade disease-related proteins by harnessing the ubiquitin-proteasome system. Due to PROTACs' hetero-bifunctional characteristics, in which a linker joins a warhead binding to a protein of interest (POI), conferring specificity and a E3-ligand binding to an E3 ubiquitin ligase, this could trigger the ubiquitination and transportation of POI to the proteasome, followed by degradation. The rational PROTAC linker design is challenging due to its relatively large molecular weight and the complexity of maintaining the binding mode of warhead and E3-ligand in the binding pockets of counterpart. Conventional linker generation method can only generate linkers in either 1D SMILES or 2D graph, without taking into account the information of ternary structures. Here we propose a novel 3D linker generative model PROTAC-INVENT which can not only generate SMILES of PROTAC but also its 3D putative binding conformation coupled with the target protein and the E3 ligase. The model is trained jointly with the RL approach to bias the generation of PROTAC structures toward pre-defined 2D and 3D based properties. Examples were provided to demonstrate the utility of the model for generating reasonable 3D conformation of PROTACs. On the other hand, our results show that the associated workflow for 3D PROTAC conformation generation can also be used as an efficient docking protocol for PROTACs.

Elucidation of Pharmacological Mechanism Underlying the Anti-Alzheimer's Disease Effects of Evodia rutaecarpa and Discovery of Novel Lead Molecules: An In Silico Study

Alzheimer's disease (AD) is a brain disease with a peculiarity of multiformity and an insidious onset. Multiple-target drugs, especially Chinese traditional medicine, have achieved a measure of success in AD treatment. Evodia rutaecarpa (Juss.) Benth. (Wuzhuyu, WZY, i.e., E. rutaecarpa), a traditional Chinese herb, has been identified as an effective drug to cure migraines. To our surprise, our in silico study showed that rather than migraines, Alzheimer's disease was the primary disease to which the E. rutaecarpa active compounds were targeted. Correspondingly, a behavioral experiment showed that E. rutaecarpa extract could improve impairments in learning and memory in AD model mice. However, the mechanism underlying the way that E. rutaecarpa compounds target AD is still not clear. For this purpose, we employed methods of pharmacology networking and molecular docking to explore this mechanism. We found that E. rutaecarpa showed significant AD-targeting characteristics, and alkaloids of E. rutaecarpa played the main role in binding to the key nodes of AD. Our research detected that E. rutaecarpa affects the pathologic development of AD through the serotonergic synapse signaling pathway (SLC6A4), hormones (PTGS2, ESR1, AR), anti-neuroinflammation (SRC, TNF, NOS3), transcription regulation (NR3C1), and molecular chaperones (HSP90AA1), especially in the key nodes of PTGS2, AR, SLCA64, and SRC. Graveoline, 5-methoxy-N, N-dimethyltryptamine, dehydroevodiamine, and goshuyuamide II in E. rutaecarpa show stronger binding affinities to these key proteins than currently known preclinical and clinical drugs, showing a great potential to be developed as lead molecules for treating AD.

Evaluation of Antioxidant Activities from a Sustainable Source of Okara Protein Hydrolysate Using Enzymatic Reaction

Okara is a solid byproduct created during the processing of soy milk. The production of protein hydrolysates utilizing enzymatic tests such as papain can result in the production of bioactive peptides (BPs), which are amino acid sequences that can also be produced from the okara protein by hydrolysis. The objective of this study was to investigate the antioxidant activities of okara hydrolysates using papain, based on the in silico and in vitro assays using the papain enzyme. We found that using the in silico assessment, the antioxidant peptides can be found from the precursor (glycinin and conglycinin) in okara. When used as a protease, papain provides the maximum degree of hydrolysis for antioxidative peptides. The highest-peptide-rank peptide sequence was predicted using peptide ranks such as proline-histidine-phenylalanine (PHF), alanine-aspartic acid-phenylalanine (ADF), tyrosine-tyrosine-leucine (YYL), proline-histidine-histidine (PHH), isoleucine-arginine (IR), and serine-valine-leucine (SVL). Molecular docking studies revealed that all peptides generated from the parent protein impeded substrate access to the active site of xanthine oxidase (XO). They have antioxidative properties and are employed in the in silico approach to the XO enzyme. We also use papain to evaluate the antioxidant activity by using in vitro tests for protein hydrolysate following proteolysis. The antioxidant properties of okara protein hydrolysates have been shown in vitro, utilizing DPPH and FRAP experiments. This study suggests that okara hydrolysates generated by papain can be employed as natural antioxidants in food and for further applications, such as active ingredients for antioxidants in packaging.

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.

Evolutionary Pattern of Interferon Alpha Genes in Bovidae and Genetic Diversity of IFNAA in the Bovine Genome

Interferons are secretory proteins induced in response to specific extracellular stimuli which stimulate intra- and intercellular networks for regulating innate and acquired immunity, resistance to viral infections, and normal and tumor cell survival and death. Type 1 interferons plays a major role in the CD8 T-cell response to viral infection. The genomic analysis carried out here for type I interferons within Bovidae family shows that cattle, bison, water buffalo, goat, and sheep (all Bovidae), have different number of genes of the different subtypes, with a large increase in the numbers, compared to human and mouse genomes. A phylogenetic analysis of the interferon alpha (IFNA) proteins in this group shows that the genes do not follow the evolutionary pattern of the species, but rather a cycle of duplications and deletions in the different species. In this study we also studied the genetic diversity of the bovine interferon alpha A (IFNAA), as an example of the IFNA genes in cattle, sequencing a fragment of the coding sequence in 18 breeds of cattle from Pakistan, Nigeria and USA. Similarity analysis allowed the allocation of sequences into 22 haplotypes. Bhagnari, Brangus, Sokoto Gudali, and White Fulani, had the highest number of haplotypes, while Angus, Hereford and Nari Master had the least. However, when analyzed by the average haplotype count, Angus, Bhagnari, Hereford, Holstein, Muturu showed the highest values, while Cholistani, Lohani, and Nari Master showed the lowest values. Haplotype 4 was found in the highest number of individuals (74), and in 15 breeds. Sequences for yak, bison, and water buffalo, were included within the bovine haplotypes. Medium Joining network showed that the sequences could be divided into 4 groups: one with highly similar haplotypes containing mostly Asian and African breeds, one with almost all of the Bos taurus American breeds, one mid-diverse group with mostly Asian and African sequences, and one group with highly divergent haplotypes with five N'Dama sequences and one from each of White Fulani, Dhanni, Tharparkar, and Bhagnari. The large genetic diversity found in IFNAA could be a very good indication of the genetic variation among the different genes of IFNA and could be an adaptation for these species in response to viral challenges they face.

Investigating the pathogenic SNPs in BLM helicase and their biological consequences by computational approach

The BLM helicase protein plays a vital role in DNA replication and the maintenance of genomic integrity. Variation in the BLM helicase gene resulted in defects in the DNA repair mechanism and was reported to be associated with Bloom syndrome (BS) and cancer. Despite extensive investigation of helicase proteins in humans, no attempt has previously been made to comprehensively analyse the single nucleotide polymorphism (SNPs) of the BLM gene. In this study, a comprehensive analysis of SNPs on the BLM gene was performed to identify, characterize and validate the pathogenic SNPs using computational approaches. We obtained SNP data from the dbSNP database version 150 and mapped these data to the genomic coordinates of the "NM_000057.3" transcript expressing BLM helicase (P54132). There were 607 SNPs mapped to missense, 29 SNPs mapped to nonsense, and 19 SNPs mapped to 3'-UTR regions. Initially, we used many consensus tools of SIFT, PROVEAN, Condel, and PolyPhen-2, which together increased the accuracy of prediction and identified 18 highly pathogenic non-synonymous SNPs (nsSNPs) out of 607 SNPs. Subsequently, these 18 high-confidence pathogenic nsSNPs were analysed for BLM protein stability, structure-function relationships and disease associations using various bioinformatics tools. These 18 mutants of the BLM protein along with the native protein were further investigated using molecular dynamics simulations to examine the structural consequences of the mutations, which might reveal their malfunction and contribution to disease. In addition, 28 SNPs were predicted as "stop gained" nonsense SNPs and one SNP was predicted as "start lost". Two SNPs in the 3'UTR were found to abolish miRNA binding and thus may enhance the expression of BLM. Interestingly, we found that BLM mRNA overexpression is associated with different types of cancers. Further investigation showed that the dysregulation of BLM is associated with poor overall survival (OS) for lung and gastric cancer patients and hence led to the conclusion that BLM has the potential to be used as an important prognostic marker for the detection of lung and gastric cancer.

Structural prediction, whole exome sequencing and molecular dynamics simulation confirms p.G118D somatic mutation of PIK3CA as functionally important in breast cancer patients

To understand the structural and functional importance of PIK3CA somatic mutations, whole exome sequencing, molecular dynamics simulation techniques in combination with in silico prediction algorithms such as SIFT, PolyPhen, Provean and CADD were employed. Twenty out of eighty missense somatic mutations in PIK3CA gene were found to be pathogenic by all the four algorithms. Most recurrent mutations found were known hotspot PIK3CA mutations with known clinical significance like p.E545 K, p.E545A, p.E545 G and p.C420R. A missense mutation p.G118D was found to be recurrently mutated in 5 cases. Interestingly, this mutation was observed in one of the patients who underwent whole exome sequencing and was completely absent from the controls. To see the effect of this mutation on the structure of PIK3CA protein, molecular dynamics simulation was performed. By molecular dynamics approach, we have shown that p.G118D mutation deviated from the native structure which was supported by the decrease in the number of hydrogen bonds, difference in hydrogen bond distance and angle, difference in root mean square deviation between the native and the mutant structures.

Novel STIL Compound Heterozygous Mutations Cause Severe Fetal Microcephaly and Centriolar Lengthening

STIL (SCL/TAL1 interrupting locus) is a core component of the centriole duplication process. STIL mutations have been associated with both autosomal recessive primary microcephaly (MCPH) and holoprosencephaly. In this report, we describe a family with multiple miscarriages and 2 terminations of pregnancy due to marked fetal microcephaly, delayed cortical gyrification, and dysgenesis of the corpus callosum. Whole exome sequencing allowed us to identify novel compound heterozygous mutations in STIL. The mutations lie, respectively, in the CPAP/CENPJ and the hsSAS6 interacting domains of STIL. M-phase synchronized amniocytes from both affected fetuses did not display an aberrant number of centrioles, as shown previously for either STIL-depleted or overexpressing cells. However, we observed an elongation of at least 1 centriole for each duplicated centrosome. These preliminary results may point to a novel mechanism causing MCPH and embryonic lethality in humans.

Prediction of a highly deleterious mutation E17K in AKT-1 gene: An in silico approach

The AKT1 (v-akt murine thymoma viral oncogene homologue 1) kinase is a member of most frequently activated proliferation and survival signaling pathway in cancer. Recently, hyperactivation of AKT1, due to functional point mutation in the pleckstrin homology (PH) domain of AKT1 gene, has been found to be associated with human colorectal, breast and ovarian cancer. Thus, considering its crucial role in cellular signaling pathway, a functional analysis of missense mutations of AKT1 gene was undertaken in this study. Twenty nine nsSNPs (non-synonymous single nucleotide polymorphism) within coding region of AKT1 gene were selected for our investigation and six SNPs were found to be deleterious by combinatorial predictions of various computational tools. RMSD values were calculated for the mutant models which predicted four substitutions (E17K, E319G, D32E and A255T) to be highly deleterious. The insight of the structural attribute was gained through analysis of, secondary structures, solvent accessibility and intermolecular hydrogen bond analysis which confirmed one missense mutation (E17K) to be highly deleterious nsSNPs. In conclusion, the investigated gene AKT1 has twenty nine SNPs in the coding region and through progressive analysis using different bioinformatics tools one highly deleterious SNP with rs121434592 was profiled. Thus, results of this study can pave a new platform to sort nsSNPs for several important regulatory genes that can be undertaken for the confirmation of their phenotype and their correlation with diseased status in case control studies.

•

We have added a small portion of text in introduction part as per reviewers comment.

•

We have added a schematic representation of methodology used (Fig. 1).

•

We have added text in the discussion portion as per the comment of reviewer.

•

We have also corrected the conclusion as per reviewer's comments.

Bioinformatic Analysis of GJB2 Gene Missense Mutations

Gap junction beta 2 (GJB2) gene is the most commonly mutated connexin gene in patients with autosomal recessive and dominant hearing loss. According to Ensembl (release 74) database, 1347 sequence variations are reported in the GJB2 gene and about 13.5% of them are categorized as missense SNPs or nonsynonymous variant. Because of the high incidence of GJB2 mutations in hearing loss patients, revealing the molecular effect of GJB2 mutations on protein structure may also provide clear point of view regarding the molecular etiology of deafness. Hence, the aim of this study is to analyze structural and functional consequences of all known GJB2 missense variations to the Cx26 protein by applying multiple bioinformatics methods. Two-hundred and eleven nonsynonymous variants were collected from Ensembl release 74, Leiden Open Variation Database (LOVD) and The Human Gene Mutation Database (HGMD). A number of bioinformatic tools were utilized for predicting the effect of GJB2 missense mutations at the sequence, structural, and functional levels. Some of the mutations were found to locate highly conserved regions and have structural and functional properties. Moreover, GJB2 mutations were also found to affect Cx26 protein at the molecular level via loss or gain of disorder, catalytic site, and post-translational modifications, including methylation, glycosylation, and ubiquitination. Findings, presented here, demonstrated the application of bioinformatic algorithms to predict the effects of mutations causing hearing impairment. I expect, this type of analysis will serve as a start point for future experimental evaluation of the GJB2 gene mutations and it will also be helpful in evaluating other deafness-related gene mutations.

Influence of MDM2 polymorphisms on squamous cell carcinoma susceptibility: a meta-analysis

PURPOSE

Controversial associations between single-nucleotide polymorphisms (rs2279744, rs937283, rs3730485) of the MDM2 gene and the etiology of squamous cell carcinomas (SCCs) have been reported. This merits further comprehensive assessment.

MATERIALS AND METHODS

We systematically reviewed the available data and conducted an updated meta-analysis to evaluate the genetic effect of MDM2 polymorphisms in SCC susceptibility, using Stata/SE 12.0 software.

RESULTS

After screening, 7,987 SCC cases and 12,954 controls from 26 eligible case-control studies were enrolled. Overall, compared with the control group, a significantly increased SCC risk was observed for the MDM2 rs2279744 polymorphism in the Asian population (test of association: odds ratio [OR] 1.12, P=0.027 for G vs T; OR 1.26, P=0.016 for GG vs TT; OR 1.25, P<0.001 for GG vs TT + TG; and OR 1.08, P=0.023 for carrier G vs T). In subgroup analysis by SCC type, a similarly increased esophageal SCC risk was detected (OR 1.19, P<0.001 for G vs T; OR 1.46, P<0.001 for GG vs TT; and OR 1.48, P=0.005 for GG vs TT + TG). Furthermore, MDM2-TP53 double mutation was statistically associated with increased SCC susceptibility overall (OR 1.52, P=0.001), especially in the Asian population (OR 1.49, P=0.022). However, no significant difference between the control and case groups was obtained for MDM2 rs937283 or rs3730485 under any genetic model (all P>0.05).

CONCLUSION

Our results highlight a positive association between the GG genotype of MDM2 rs2279744 polymorphism and an increased risk of esophageal SCC in the Asian population, which needs to be clarified by more large-scale studies.

Prediction of the damage-associated non-synonymous single nucleotide polymorphisms in the human MC1R gene

The melanocortin 1 receptor (MC1R) is involved in the control of melanogenesis. Polymorphisms in this gene have been associated with variation in skin and hair color and with elevated risk for the development of melanoma. Here we used 11 computational tools based on different approaches to predict the damage-associated non-synonymous single nucleotide polymorphisms (nsSNPs) in the coding region of the human MC1R gene. Among the 92 nsSNPs arranged according to the predictions 62% were classified as damaging in more than five tools. The classification was significantly correlated with the scores of two consensus programs. Alleles associated with the red hair color (RHC) phenotype and with the risk of melanoma were examined. The R variants D84E, R142H, R151C, I155T, R160W and D294H were classified as damaging by the majority of the tools while the r variants V60L, V92M and R163Q have been predicted as neutral in most of the programs The combination of the prediction tools results in 14 nsSNPs indicated as the most damaging mutations in MC1R (L48P, R67W, H70Y, P72L, S83P, R151H, S172I, L206P, T242I, G255R, P256S, C273Y, C289R and R306H); C273Y showed to be highly damaging in SIFT, Polyphen-2, MutPred, PANTHER and PROVEAN scores. The computational analysis proved capable of identifying the potentially damaging nsSNPs in MC1R, which are candidates for further laboratory studies of the functional and pharmacological significance of the alterations in the receptor and the phenotypic outcomes.

Making an impact: FEBS Open Bio celebrates its 3rd anniversary

As FEBS celebrates its 50th anniversary, it also marks the launch 3 years ago of its newest journal, FEBS Open Bio. This online-only open access journal is now indexed in Web of Science. Papers originally submitted to The FEBS Journal can be transferred to FEBS Open Bio together with any reviewer reports, offering authors a fast-track alternative for publication of their work.

Computational SNP analysis: current approaches and future prospects

The computational approaches in determining disease-associated Non-synonymous single nucleotide polymorphisms (nsSNPs) have evolved very rapidly. Large number of deleterious and disease-associated nsSNP detection tools have been developed in last decade showing high prediction reliability. Despite of all these highly efficient tools, we still lack the accuracy level in determining the genotype-phenotype association of predicted nsSNPs. Furthermore, there are enormous questions that are yet to be computationally compiled before we might talk about the prediction accuracy. Earlier we have incorporated molecular dynamics simulation approaches to foster the accuracy level of computational nsSNP analysis roadmap, which further helped us to determine the changes in the protein phenotype associated with the computationally predicted disease-associated mutation. Here we have discussed on the present scenario of computational nsSNP characterization technique and some of the questions that are crucial for the proper understanding of pathogenicity level for any disease associated mutations.

Evidence of colorectal cancer-associated mutation in MCAK: a computational report

Computational prediction of disease-associated non-synonymous polymorphism (nsSNP) has provided a significant platform to filter out the pathological mutations from large pool of SNP datasets at a very low cost input. Several methodologies and complementary protocols have been previously implemented and has provided significant prediction results. Although the previously implicated prediction methods were capable of investigating the most likely deleterious nsSNPs, but due to the lack of genotype-phenotype association analysis, the prediction results lacked in accuracy level. In this work we implemented the computational compilation of protein conformational changes as well as the probable disease-associated phenotypic outcomes. Our result suggested E403K mutation in mitotic centromere-associated kinesin protein as highly damaging and showed strong concordance to the previously observed colorectal cancer mutations aggregation tendency and energy value changes. Moreover, the molecular dynamics simulation results showed major loss in conformation and stability of mutant N-terminal kinesin-like domain structure. The result obtained in this study will provide future prospect of computational approaches in determining the SNPs that may affect the native conformation of protein structure and lead to cancer-associated disorders.

Investigation of binding phenomenon of NSP3 and p130Cas mutants and their effect on cell signalling

Members of the novel SH2-containing protein (NSP3) and Crk-associated substrate (p130Cas) protein families form a multi-domain signalling platforms that mediate cell signalling process. We analysed the damaging consequences of three mutations, each from NSP3 (NSP3(L469R), NSP3(L623E), NSP3(R627E)) and p130Cas (p130Cas(F794R), p130Cas(L787E), p130Cas(D797R)) protein with respect to their native biological partners. Mutations depicted notable loss in interaction affinity towards their corresponding biological partners. NSP3(L469R) and p130Cas(D797R) mutations were predicted as most prominent in docking analysis. Molecular dynamics (MD) studies were conducted to evaluate structural consequences of most prominent mutation in NSP3 and p130Cas obtained from the docking analysis. MD analysis confirmed that mutation in NSP3(L469R) and p130Cas(D797R) showed significant structural deviation, changes in conformations and increased flexibility, which in turn affected the binding affinity with their biological partners. Moreover, the root mean square fluctuation has indicated a rise in fluctuation of residues involved in moderate interaction acquired between the NSP3 and p130Cas. It has significantly affected the binding interaction in mutant complexes. The results obtained in this work present a detailed overview of molecular mechanisms involved in the loss of cell signalling associated with NSP3 and p130Cas protein.

Computational investigation of cancer-associated molecular mechanism in Aurora A (S155R) mutation

Centrosomes are the key-regulating element of cell cycle progression. Aberrations in their functional mechanism lead to several cancer-related disorders. Aurora A protein is a centrosome-associated protein that regulates the centriole duplication and its abberations are associated with multiple cases of aneuploidy and cancer-related disorders. S155R mutation in Aurora A is reported to induce cancer like phenotype and disrupt its binding with TPX2 protein. In this study, we have demonstrated the structural consequences of Aurora A S155R mutation and the atomic changes that influenced the loss of TPX2-binding affinity. Docking and molecular dynamics simulation results suggested significant loss in atomic contacts between mutant Aurora A and TPX2 protein. Further, we observed a notable changes in conformation of mutant Aurora A-TPX2 docked complex as compared to the native. Loss of binding affinity rendered the TPX2 domain free which then induced unfolding in its coiled region and enabled the overall expansion of mutant complex as compared to the native. The significant outcomes obtained from this study will facilitate in future cancer researches and in developing the potent drug therapies.

AKT kinase pathway: a leading target in cancer research

AKT1, a serine/threonine-protein kinase also known as AKT kinase, is involved in the regulation of various signalling downstream pathways including metabolism, cell proliferation, survival, growth, and angiogenesis. The AKT kinases pathway stands among the most important components of cell proliferation mechanism. Several approaches have been implemented to design an efficient drug molecule to target AKT kinases, although the promising results have not been confirmed. In this paper we have documented the detailed molecular insight of AKT kinase protein and proposed a probable doxorubicin based approach in inhibiting miR-21 based cancer cell proliferation. Moreover, the inhibition of miR-21 activation by raising the FOXO3A concentration seems promising in reducing miR-21 mediated cancer activation in cell. Furthermore, the use of next generation sequencing and computational drug design approaches will greatly assist in designing a potent drug molecule against the associated cancer cases.

In silico screening and molecular dynamics simulation of disease-associated nsSNP in TYRP1 gene and its structural consequences in OCA3

Oculocutaneous albinism type III (OCA3), caused by mutations of TYRP1 gene, is an autosomal recessive disorder characterized by reduced biosynthesis of melanin pigment in the hair, skin, and eyes. The TYRP1 gene encodes a protein called tyrosinase-related protein-1 (Tyrp1). Tyrp1 is involved in maintaining the stability of tyrosinase protein and modulating its catalytic activity in eumelanin synthesis. Tyrp1 is also involved in maintenance of melanosome structure and affects melanocyte proliferation and cell death. In this work we implemented computational analysis to filter the most probable mutation that might be associated with OCA3. We found R326H and R356Q as most deleterious and disease associated by using PolyPhen 2.0, SIFT, PANTHER, I-mutant 3.0, PhD-SNP, SNP&GO, Pmut, and Mutpred tools. To understand the atomic arrangement in 3D space, the native and mutant (R326H and R356Q) structures were modelled. Finally the structural analyses of native and mutant Tyrp1 proteins were investigated using molecular dynamics simulation (MDS) approach. MDS results showed more flexibility in native Tyrp1 structure. Due to mutation in Tyrp1 protein, it became more rigid and might disturb the structural conformation and catalytic function of the structure and might also play a significant role in inducing OCA3. The results obtained from this study would facilitate wet-lab researches to develop a potent drug therapies against OCA3.

Cancer associated E17K mutation causes rapid conformational drift in AKT1 pleckstrin homology (PH) domain

BACKGROUND

AKT1 (v-akt murine thymoma viral oncogene homologue 1) kinase is one of the most frequently activated proliferated and survival pathway of cancer. Recently it has been shown that E17K mutation in the Pleckstrin Homology (PH) domain of AKT1 protein leads to cancer by amplifying the phosphorylation and membrane localization of protein. The mutant has shown resistance to AKT1/2 inhibitor VIII drug molecule. In this study we have demonstrated the detailed structural and molecular consequences associated with the activity regulation of mutant protein.

METHODS

The docking score exhibited significant loss in the interaction affinity to AKT1/2 inhibitor VIII drug molecule. Furthermore, the molecular dynamics simulation studies presented an evidence of rapid conformational drift observed in mutant structure.

RESULTS

There was no stability loss in mutant as compared to native structure and the major cation-π interactions were also shown to be retained. Moreover, the active residues involved in membrane localization of protein exhibited significant rise in NHbonds formation in mutant. The rise in NHbond formation in active residues accounts for the 4-fold increase in the membrane localization potential of protein.

CONCLUSION

The overall result suggested that, although the mutation did not induce any stability loss in structure, the associated pathological consequences might have occurred due to the rapid conformational drifts observed in the mutant AKT1 PH domain.

GENERAL SIGNIFICANCE

The methodology implemented and the results obtained in this work will facilitate in determining the core molecular mechanisms of cancer-associated mutations and in designing their potential drug inhibitors.