Search

Virus assembly, budding, or surface proteins play important roles such as viral attachment to cells, fusion, and entry into cells. The present study aimed to identify potential antigenic proteins and epitopes that could be used to develop a vaccine or diagnostic assay against the Monkeypox virus (MPXV) which may cause a potential epidemic. To do this, 39 MPXV proteins (including assembly, budding, and surface proteins) were analyzed using an in silico approach. Of these 39 proteins, the F5L virus protein was found to be the best vaccine candidate due to its signal peptide properties, negative GRAVY value, low transmembrane helix content, moderate aliphatic index, large molecular weight, long-estimated half-life, beta wrap motifs, and being stable, soluble, and containing non-allergic features. Moreover, the F5L protein exhibited alpha-helical secondary structures, making it a potential "structural antigen" recognized by antibodies. The other viral protein candidates were A9 and A43, but A9 lacked beta wrap motifs, while A43 had a positive GRAVY value and was insoluble. These two proteins were not as suitable candidates as the F5L protein. The KRVNISLTCL epitope from the F5L protein demonstrated the highest antigen score (2.4684) for MHC-I, while the GRFGYVPYVGYKCI epitope from the A9 protein exhibited the highest antigenicity (1.754) for MHC-II. Both epitopes met the criteria for high antigenicity, non-toxicity, solubility, non-allergenicity, and the presence of cleavage sites. Molecular docking and dynamics (MD) simulations further validated their potential, revealing stable and energetically favorable interactions with MHC molecules. The immunogenicity assessment showed that GRFGYVPYVGYKCI could strongly induce immune responses through both IFN-γ and IL-4 pathways, suggesting its capacity to provoke a balanced Th1 and Th2 response. In contrast, KRVNISLTCL exhibited limited immunostimulatory potential. Overall, these findings lay the groundwork for future vaccine development, indicating that F5L, particularly the GRFGYVPYVGYKCI epitope, may serve as an effective candidate for peptide-based vaccine design against MPXV., Competing Interests: Declarations Conflict of interest The authors declare no competing interests., (© 2024. The Author(s), under exclusive licence to Springer Nature Switzerland AG.)

Background: The human gut microbiota is an extensive population of microorganisms, and it shows significant variations between periods of optimal health and periods of illness. Vancomycin-resistant Enterococcus (VRE) and carbapenem-resistant Klebsiella pneumoniae (CRKP) are both pathogenic agents (BPAs) that can colonize in the gut after dysbiosis of microbiotal composition following antibiotic treatment. Methods: This study aimed to investigate the impact of antibiotics on the microbiotal composition of the gut. For this purpose, the first pass meconiums of 20 patients and the first rectal swabs containing BPAs of the same patients after antibiotic treatment were studied using next-generation sequencing-based 16S rRNA gene analysis. The V1-V9 region of 16S rRNA was sequenced with Oxford Nanopore. Results: Twenty-five phyla were detected in the meconiums, and 12 of them were absent after antibiotic treatment. The four most prevalent phyla in meconiums were Bacillota, Pseudomonadota, Bacteroidota, and Actinomycetota. Only the relative abundance of Pseudomonadota was increased, while a significant decrease was observed in the other three phyla ( p < 0.05). A significant decrease was observed in alpha-diversity in rectal swabs containing BPAs versus meconiums ( p = 0.00408), whereas an increased variance was observed in beta-diversity in all samples ( p < 0.05). As a result of a LEfSe analysis, Pseudomonadota was found to have a higher relative abundance in rectal swabs, and Bacillota was significantly higher in the meconiums of the twins. Conclusions: Our study strongly verified the relationship between the administration of antibiotics, dysbiosis, and colonization of BPAs in the infants' gut microbiota. Further research would be beneficial and needed, comprising the natural development process of the infants' gut microbiota.

The assessment of kidney function within the first year following transplantation is crucial for predicting long-term graft survival. This study aimed to develop a robust and accurate model using metabolite profiles to predict early long-term outcomes in patient groups at the highest risk of early graft loss. A group of 61 kidney transplant recipients underwent thorough monitoring during a one-year follow-up period, which included a one-week hospital stay and follow-up assessments at three and six months. Based on their 12-month follow-up serum creatinine levels: Group 2 had levels exceeding 1.5 mg/dl, while Group 1 had levels below 1.5 mg/dl. Metabolites were detected by mass spectrometer and first pre-processed. Univariate and multivariate statistical analyses were employed to identify significant differences between the two groups. Nineteen metabolites were found to differ significantly in the 1 st week, and seventeen metabolites in the 3 rd month (adjusted p-value < 0.05, quality control (QC) < 30, a fold change (FC) > 1.1 or a FC < 0.91, Variable Influence on Projection (VIP) > 1). However, no significant differences were observed in the 6 th month. These distinctive metabolites mainly belonged to lipid, fatty acid, and amino acid categories. Ten models were constructed using a backward conditional approach, with the best performance seen in model 5 for Group 2 at the 1st-week mark (AUC 0.900) and model 3 at the 3 rd -month mark (AUC 0.924). In conclusion, the models developed in the early stages may offer potential benefits in the management of kidney transplant patients.

MYC dysregulation is pivotal in the onset and progression of IDH-mutant gliomas, mostly driven by copy-number alterations, regulatory element alterations, or epigenetic changes. Our pilot analysis uncovered instances of relative MYC overexpression without alterations in the proximal MYC network (PMN), prompting a deeper investigation into potential novel oncogenic mechanisms. Analysing comprehensive genomics profiles of 236 "IDH-mutant 1p/19q non-co-deleted" lower-grade gliomas from The Cancer Genome Atlas, we identified somatic genomic alterations within the PMN. In tumours without PMN-alterations but with MYC-overexpression, genes correlated with MYC-overexpression were identified. Our analyses yielded that 86/236 of astrocytomas exhibited no PMN-alterations, a subset of 21/86 displaying relative MYC overexpression. Within this subset, we discovered 42 genes inversely correlated with relative MYC expression, all on 19q. Further analysis pinpointed a minimal common region at 19q13.43, encompassing 15 genes. The inverse correlations of these 15 genes with relative MYC overexpression were re-confirmed using independent scRNAseq data. Further, the micro-deleted astrocytoma subset displayed significantly higher genomic instability compared to WT cases, but lower instability compared to PMN-hit cases. This newly identified 19q micro-deletion represents a potential novel mechanism underlying MYC dysregulation in astrocytomas. Given the prominence of 19q loss in IDH-mutant gliomas, our findings bear significant implications for understanding gliomagenesis., (© 2024. The Author(s).)

Solid cancers like pancreatic ductal adenocarcinoma (PDAC), a type of pancreatic cancer, frequently exploit nerves for rapid dissemination. This neural invasion (NI) is an independent prognostic factor in PDAC, but insufficiently modeled in genetically engineered mouse models (GEMM) of PDAC. Here, we systematically screened for human-like NI in Europe's largest repository of GEMM of PDAC, comprising 295 different genotypes. This phenotype screen uncovered 2 GEMMs of PDAC with human-like NI, which are both characterized by pancreas-specific overexpression of transforming growth factor α (TGF-α) and conditional depletion of p53. Mechanistically, cancer-cell-derived TGF-α upregulated CCL2 secretion from sensory neurons, which induced hyperphosphorylation of the cytoskeletal protein paxillin via CCR4 on cancer cells. This activated the cancer migration machinery and filopodia formation toward neurons. Disrupting CCR4 or paxillin activity limited NI and dampened tumor size and tumor innervation. In human PDAC, phospho-paxillin and TGF-α-expression constituted strong prognostic factors. Therefore, we believe that the TGF-α-CCL2-CCR4-p-paxillin axis is a clinically actionable target for constraining NI and tumor progression in PDAC.

Helicobacter pylori OipA (Outer Inflammatory Protein A) is an outer membrane protein that takes role in the adherence and colonization to the stomach. oipA gene expression is regulated by the slipped-strand mispairing mechanism through a hypermutable CT dinucleotide repeat motif in the 5΄ region. Alterations in the CT number repeats cause frame-shift mutations to result in phase variation of oipA expression. While a functional "On" status has been recognized as a risk factor for peptic ulcer diseases and gastric cancer in many studies, some controversial findings still exist. To this end, this study compiled the sequence data of oipA from 10 different studies between 2000-2019 and 50 oipA DNA sequences from our own research that examined the relationship between the phase On/Off status of oipA and gastric diseases based on CT repeat number. Overall, we have reached 536 oipA DNA sequences from patients. This large collection of oipA sequences first clarified the absolute conservation of the peptide-pentamer of FWLHA for phase ''On'' status, suggesting this pentamer as a superior marker for the determination of oipA status than counting the number of CT repeats. Combining the sequence and patient data, we have re-analyzed the association between the ''On'' status of oipA and gastric diseases. Our results showed a strong association between oipA ''On'' status and gastric cancer supporting previous findings. We also investigated the AlphaFold2 computed structure of OipA that adopts a beta-barrel fold closely resembling to the autotransporter family of H. pylori. Altogether, this study confirms a strong association between oipA ''On'' statuses and severe gastrointestinal diseases like cancer and provides useful insights into the FWLHA pentamer as an indicator of "On" status of oipA putative autotransporter function rather than CT repeats number., (© 2023. Springer-Verlag GmbH Germany, part of Springer Nature.)

Purpose: To investigate the corneal bacterial microbiome in patients with keratoconus using next-generation sequencing and develop a new perspective on the pathogenesis of the disease., Methods: This prospective observational study included 10 patients with keratoconus who underwent corneal crosslinking procedure and 10 healthy controls who underwent photorefractive keratectomy. Patients included in the study were aged 18 years or older. The demographic and clinical characteristics of participants were recorded. Corneal epithelial samples were collected between March 2021 and June 2021. Isolated bacterial DNA from corneal epithelial samples was analyzed using 16 S ribosomal RNA gene analysis. The relative abundance rates at the phylum and genus levels were calculated. Alpha diversity parameters were assessed., Results: Eleven phyla and 521 genera of bacteria were identified in all participants. At the phylum level, Proteobacteria, Firmicutes, Actinobacteria, and Bacteroidetes were most abundant in both groups. There were no statistical differences between the two groups except Bacteriodetes (p < 0.05). At the genus level, the relative abundance rates of twenty bacteria were significantly different between keratoconus and healthy corneas (p < 0.05). Aquabacterium was the most abundant genus in patients with keratoconus, while Shigella was the most abundant genus in healthy controls. Alpha diversity parameters were lower in patients with keratoconus, although the difference did not reach statistical significance (p > 0.05)., Conclusions: Our preliminary study revealed that there are similarities and differences in the corneal microbiome between keratoconus and healthy individuals. Further research is required on the relationship between the abnormal corneal microbiome composition and the pathogenesis of keratoconus., Competing Interests: Declaration of competing interest None., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Cystinosis is a rare, devastating hereditary disease secondary to recessive CTNS gene mutations. The most commonly used diagnostic method is confirmation of an elevated leukocyte cystine level; however, this method is expensive and difficult to perform. This study aimed to identify candidate biomarkers for the diagnosis and follow-up of cystinosis based on multiomics studies. The study included three groups: newly-diagnosed cystinosis patients (patient group, n = 14); cystinosis patients under treatment (treatment group, n = 19); and healthy controls (control group, n = 30). Plasma metabolomics analysis identified 10 metabolites as candidate biomarkers that differed between the patient and control groups [L-serine, taurine, lyxose, 4-trimethylammoniobutanoic acid, orotic acid, glutathione, PE(O-18:1(9Z)/0:0), 2-hydroxyphenyl acetic acid, acetyl-N-formil-5-metoxikinuramine, 3-indoxyl sulphate]. As compared to the healthy control group, in the treatment group, hypotaurine, phosphatidylethanolamine, N-acetyl-d-mannosamine, 3-indolacetic acid, p-cresol, phenylethylamine, 5-aminovaleric acid, glycine, creatinine, and saccharic acid levels were significantly higher, and the metabolites quinic acid, capric acid, lenticin, xanthotoxin, glucose-6-phosphate, taurine, uric acid, glyceric acid, alpha-D-glucosamine phosphate, and serine levels were significantly lower. Urinary metabolomic analysis clearly differentiated the patient group from the control group by means of higher allo-inositol, talose, glucose, 2-hydroxybutiric acid, cystine, pyruvic acid, valine, and phenylalanine levels, and lower metabolite (N-acetyl-L-glutamic acid, 3-aminopropionitrile, ribitol, hydroquinone, glucuronic acid, 3-phosphoglycerate, xanthine, creatinine, and 5-aminovaleric acid) levels in the patient group. Urine metabolites were also found to be significantly different in the treatment group than in the control group. Thus, this study identified candidate biomarkers that could be used for the diagnosis and follow-up of cystinosis.

Non-ionic surfactants were shown to stabilize the active conformation of thermoalkalophilic lipases by mimicking the lipid substrate while the catalytic interactions formed by anionic surfactants have not been well characterized. In this study, we combined μs-scale molecular dynamics (MD) simulations and lipase activity assays to analyze the effect of ionic surfactant, sodium dodecyl sulfate (SDS), on the structure and activity of thermoalkalophilic lipases. Both the open and closed lipase conformations that differ in geometry were recruited to the MD analysis to provide a broader understanding of the molecular effect of SDS on the lipase structure. Simulations at 298 K showed the potential of SDS for maintaining the active lipase through binding to the sn -1 acyl-chain binding pocket in the open conformation or transforming the closed conformation to an open-like state. Consistent with MD findings, experimental analysis showed increased lipase activity upon SDS incubation at ambient temperature. Notably, the lipase cores stayed intact throughout 2 μs regardless of an increase in the simulation temperature or SDS concentration. However, the surface structures were unfolded in the presence of SDS and at elevated temperature for both conformations. Simulations of the dimeric lipase were also carried out and showed reduced flexibility of the surface structures which were unfolded in the monomer, indicating the insulating role of dimer interactions against SDS. Taken together, this study provides insights into the possible substrate mimicry by the ionic surfactant SDS for the thermoalkalophilic lipases without temperature elevation, underscoring SDS's potential for interfacial activation at ambient temperatures.

Motivation: Identifying appropriate pharmacotherapy options from genomics results is a significant challenge in personalized oncology. However, computational methods for prioritizing drugs are underdeveloped. With the hypothesis that network-based approaches can improve the performance by extending the use of potential drug targets beyond direct interactions, we devised two network-based methods for personalized pharmacotherapy prioritization in cancer., Results: We developed novel personalized drug prioritization approaches, PANACEA: PersonAlized Network-based Anti-Cancer therapy EvaluAtion. In PANACEA, initially, the protein interaction network is extended with drugs, and a driverness score is assigned to each altered gene. For scoring drugs, either (i) the 'distance-based' method, incorporating the shortest distance between drugs and altered genes, and driverness scores, or (ii) the 'propagation' method involving the propagation of driverness scores via a random walk with restart framework is performed. We evaluated PANACEA using multiple datasets, and demonstrated that (i) the top-ranking drugs are relevant for cancer pharmacotherapy using TCGA data; (ii) drugs that cancer cell lines are sensitive to are identified using GDSC data; and (iii) PANACEA can perform adequately in the clinical setting using cases with known drug responses. We also illustrate that the proposed methods outperform iCAGES and PanDrugs, two previous personalized drug prioritization approaches., Availability and Implementation: The corresponding R package is available on GitHub. (https://github.com/egeulgen/PANACEA.git)., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2023. Published by Oxford University Press.)

Background Knowledge about cutaneous microbiota in psoriasis vulgaris and seborrheic dermatitis is limited, and a comparison of microbiota in the two diseases was not yet previously undertaken. Aims/Objectives This study aimed to compare the scalp lesional and non-lesional microbiota in psoriasis vulgaris and seborrheic dermatitis with that in a healthy control group. Methods Fifty samples were taken with sterile swabs from patients' and controls' scalps, and 16S rRNA gene sequencing analyses were performed. Results Alpha and beta diversity analyses showed that bacterial load and diversity were significantly increased in psoriasis vulgaris and seborrheic dermatitis lesions compared to the controls. As phyla, Actinobacteria decreased and Firmicutes increased, while as genera, Propionibacterium decreased; Staphylococcus, Streptococcus, Aquabacterium, Neisseria and Azospirillum increased in lesions of both diseases. Specifically, Mycobacterium, Finegoldia, Haemophilus and Ezakiella increased in psoriasis vulgaris and Enhydrobacter, Micromonospora and Leptotrichia increased in seborrheic dermatitis lesions. Mycobacterium, Ezakiella and Peptoniphilus density were higher in psoriasis vulgaris compared to seborrheic dermatitis lesions. The bacterial diversity and load values of non-lesional scalp in psoriasis vulgaris and seborrheic dermatitis lay between those of lesional areas and controls. Limitations The small sample size is the main limitation of this study. Conclusion Higher bacterial diversity was detected in lesions of both psoriasis and seborrheic dermatitis compared to the controls, but similar alterations were observed when the two diseases were compared. Although these differences could be a result rather than a cause of the two diseases, there is a need to analyze all members of the microbiota and microbiota-host interactions.

Modulating the activity of human soluble guanylate cyclase ( h sGC) through allosteric regulation of the βH-NOX domain has been considered as an immediate treatment for cardiovascular disorder (CVDs). Currently available βH-NOX domain-specific agonists including cinaciguat are unable to deal with the conundrum raised due to oxidative stress in the case of CVDs and their associated comorbidities. Therefore, the idea of investigating novel compounds for allosteric regulation of h sGC activation has been rekindled to circumvent CVDs. Current study aims to identify novel βH-NOX domain-specific compounds that can selectively turn on sGC functions by modulating the conformational dynamics of the target protein. Through a comprehensive computational drug-discovery approach, we first executed a target-based performance assessment of multiple docking (PLANTS, QVina, LeDock, Vinardo, Smina) scoring functions based on multiple performance metrices. QVina showed the highest capability of selecting true-positive ligands over false positives thus, used to screen 4.8 million ZINC15 compounds against βH-NOX domain. The docked ligands were further probed in terms of contact footprint and pose reassessment through clustering analysis and PLANTS docking, respectively. Subsequently, energy-based AMBER rescoring of top 100 low-energy complexes, per-residue energy decomposition analysis, and ADME-Tox analysis yielded the top three compounds i.e. ZINC000098973660, ZINC001354120371, and ZINC000096022607. The impact of three selected ligands on the internal structural dynamics of the βH-NOX domain was also investigated through molecular dynamics simulations. The study revealed potential electrostatic interactions for better conformational dialogue between βH-NOX domain and allosteric ligands that are critical for the activation of h sGC as compared to the reference compound.

Recent studies have focused on the early detection of ovarian cancer (OC) using tumor materials by liquid biopsy. The mechanisms of microRNAs (miRNAs) to impact OC and signaling pathways are still unknown. This study aims to reliably perform functional analysis of previously validated circulating miRNAs' target genes by using pathfindR. Also, overall survival and pathological stage analyses were evaluated with miRNAs' target genes which are common in the The Cancer Genome Atlas and GTEx datasets. Our previous studies have validated three downregulated miRNAs (hsa-miR-885-5p, hsa-miR-1909-5p, and hsalet7d-3p) having a diagnostic value in OC patients' sera, with high-throughput techniques. The predicted target genes of these miRNAs were retrieved from the miRDB database (v6.0). Active-subnetwork-oriented Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was conducted by pathfindR using the target genes. Enrichment of KEGG pathways assessed by the analysis of pathfindR indicated that 24 pathways were related to the target genes. Ubiquitin-mediated proteolysis, spliceosome and Notch signaling pathway were the top three pathways with the lowest p-values (p < 0.001). Ninety-three common genes were found to be differentially expressed (p < 0.05) in the datasets. No significant genes were found to be significant in the analysis of overall survival analyses, but 24 genes were found to be significant with pathological stages analysis (p < 0.05). The findings of our study provide in-silico evidence that validated circulating miRNAs' target genes and enriched pathways are related to OC and have potential roles in theranostics applications. Further experimental investigations are required to validate our results which will ultimately provide a new perspective for translational applications in OC management.

Background: Periventricular nodular heterotopia (PNH) is a cell migration disorder associated with mutations in Filamin-A (FLNA) gene on chromosome X. Majority of the individuals with PNH-associated FLNA mutations are female whereas liveborn males with FLNA mutations are very rare. Fetal viability of the males seems to depend on the severity of the variant. Splicing or severe truncations presumed loss of function of the protein product, lead to male lethality and only partial-loss-of-function variants are reported in surviving males. Those variants mostly manifest milder clinical phenotypes in females and thus avoid detection of the disease in females., Methods: We describe a novel p.Arg484Gln variant in the FLNA gene by performing whole exome analysis on the index case, his one affected brother and his healthy non-consanguineous parents. The transmission of PNH from a clinically asymptomatic mother to two sons is reported in a fully penetrant classical X-linked dominant mode. The variant was verified via Sanger sequencing. Additionally, we investigated the impact of missense mutations reported in affected males on the FLNa protein structure, dynamics and interactions by performing molecular dynamics (MD) simulations to examine the disease etiology and possible compensative mechanisms allowing survival of the males., Results: We observed that p.Arg484Gln disrupts the FLNa by altering its structural and dynamical properties including the flexibility of certain regions, interactions within the protein, and conformational landscape of FLNa. However, these impacts existed for only a part the MD trajectories and highly similar patterns observed in the other 12 mutations reported in the liveborn males validated this mechanism., Conclusion: It is concluded that the variants seen in the liveborn males result in transient pathogenic effects, rather than persistent impairments. By this way, the protein could retain its function occasionally and results in the survival of the males besides causing the disease., Competing Interests: The authors have declared that no competing interests exist.

Purpose: Blink reflex (BR) under general anesthesia as an intraoperative neuromonitoring method was used to monitor facial nerves in few studies. This study aimed to test the utility of intraoperative BR during cerebellopontine angle and skull base surgeries, assess its prognostic value for facial nerve functions, and compare it with facial corticobulbar motor evoked potentials (CoMEPs)., Methods: Blink reflex and facial CoMEPs were recorded from 40 patients undergoing skull base surgeries. Subdermal needles were placed in the supraorbital notch for stimulation and in the orbicularis oculi muscle for recording the BR. A double train of 20 to 40 V intensity with an intertrain interval of 40 to 60 milliseconds, an interstimulus interval of 2.5 milliseconds, and a stimulus duration of 0.5 milliseconds were applied. Facial nerve functions were assessed with the House-Brackmann grading system in the postoperative day 1 and third-month period and correlated with intraoperative BR and CoMEPs measurements., Results: Of 40 patients, BR was recordable on the affected side in 32 (80%) and contralateral side in 35 (87.5%) patients. According to our statistical results, BR had a slightly better sensitivity than facial CoMEPs in predicting impairment of facial nerve functions for both postoperative and third-month time points. Blink reflex showed better accuracy for predicting postoperative nerve functions, whereas CoMEPs correlated better in predicting third-month outcome., Conclusions: We suggest that BR is a valuable intraoperative neuromonitoring method that can be used in addition to facial CoMEPs during skull base surgeries to assess real-time facial nerve integrity and predict prognosis., Competing Interests: The authors have no funding or conflicts of interest to disclose., (Copyright © 2020 by the American Clinical Neurophysiology Society.)

Background: Impaired meiosis can result in absence of sperm in the seminal fluid. This condition, namely non-obstructive azoospermia (NOA), is one of the reasons of male infertility. Despite the low number of studies on meiosis 1-associated protein ( M1AP ) in the literature, M1AP is known to be crucial for spermatogenesis. Recently, seven variants (five missense, one frameshift, one splice-site) have been reported in the M1AP gene as associated with NOA, cryptozoospermia and oligozoospermia in two separate studies. However, all missense variants were evaluated as variant of uncertain significance by these studies. Therefore, we aimed to analyze their structural impacts on the M1AP protein that could lead to NOA., Methods: We firstly performed an evolutionary conservation analysis for the variant positions. Afterwards, a comprehensive molecular modelling study was performed for the M1AP structure. By utilizing this model, protein dynamics were sampled for the wild-type and variants by performing molecular dynamics (MD) simulations., Results: All variant positions are highly conserved, indicating that they are potentially important for function. In MD simulations, none of the variants led to a general misfolding or loss of stability in the protein structure, but they did cause severe modifications in the conformational dynamics of M1AP, particularly through changes in local interactions affecting flexibility, hinge and secondary structure., Conclusions: Due to critical perturbations in protein dynamics, we propose that these variants may cause NOA by affecting important interactions regulating meiosis, particularly in wild-type M1AP deficiency since the variants are reported to be homozygous or bi-allelic in the infertile individuals. Our results provided reasonable insights about the M1AP structure and the effects of the variants to the structure and dynamics, which should be further investigated by experimental studies to validate., Competing Interests: The authors declare that they have no competing interests., (© 2022 Gerlevik et al.)

Recessive mutations in the genes encoding the four subunits of the tRNA splicing endonuclease complex (TSEN54, TSEN34, TSEN15, and TSEN2) cause various forms of pontocerebellar hypoplasia, a disorder characterized by hypoplasia of the cerebellum and the pons, microcephaly, dysmorphisms, and other variable clinical features. Here, we report an intronic recessive founder variant in the gene TSEN2 that results in abnormal splicing of the mRNA of this gene, in six individuals from four consanguineous families affected with microcephaly, multiple craniofacial malformations, radiological abnormalities of the central nervous system, and cognitive retardation of variable severity. Remarkably, unlike patients with previously described mutations in the components of the TSEN complex, all the individuals that we report developed atypical hemolytic uremic syndrome (aHUS) with thrombotic microangiopathy, microangiopathic hemolytic anemia, thrombocytopenia, proteinuria, severe hypertension, and end-stage kidney disease (ESKD) early in life. Bulk RNA sequencing of peripheral blood cells of four affected individuals revealed abnormal tRNA transcripts, indicating an alteration of the tRNA biogenesis. Morpholino-mediated skipping of exon 10 of tsen2 in zebrafish produced phenotypes similar to human patients. Thus, we have identified a novel syndrome accompanied by aHUS suggesting the existence of a link between tRNA biology and vascular endothelium homeostasis, which we propose to name with the acronym TRACK syndrome (TSEN2 Related Atypical hemolytic uremic syndrome, Craniofacial malformations, Kidney failure)., (© 2021 John Wiley & Sons A/S . Published by John Wiley & Sons Ltd.)

Phenylketonuria is an autosomal recessive disorder caused by mutations in the phenylalanine hydroxylase gene. In phenylketonuria causes various symptoms including severe mental retardation. PAH gene of a classical Phenylketonuria patient was sequenced, and two novel heterozygous mutations, p.Y198N and p.Y204F, were found. This study aimed to reveal the impacts of these variants on the structural stability of the PAH enzyme. In-silico analyses using prediction tools and molecular dynamics simulations were performed. Mutations were introduced to the wild type catalytic monomer and full length tetramer crystal structures. Variant pathogenicity analyses predicted p.Y198N to be damaging, and p.Y204F to be benign by some prediction tools and damaging by others. Simulations suggested p.Y198N mutation cause significant fluctuations in the spatial organization of two catalytic residues in the temperature accelerated MD simulations with the monomer and increased root-mean-square deviations in the tetramer structure. p.Y204F causes noticeable changes in the spatial positioning of T278 suggesting a possible segregation from the catalytic site in temperature accelerated MD simulations with the monomer. This mutation also leads to increased root-mean-square fluctuations in the regulatory domain which may lead to conformational change resulting in inhibition of dimerization and enzyme activation. Our study reports two novel mutations in the PAH gene and gives insight to their effects on the PAH activity. MD simulations did not yield conclusive results that explains the phenotype but gave plausible insight to possible effects which should be investigated further with in-silico and in-vitro studies to assess the roles of these mutations in etiology of PKU. Communicated by Ramaswamy H. Sarma.

Background: Artificial intelligence (AI) has the potential to transform our healthcare systems significantly. New AI technologies based on machine learning approaches should play a key role in clinical decision-making in the future. However, their implementation in health care settings remains limited, mostly due to a lack of robust validation procedures. There is a need to develop reliable assessment frameworks for the clinical validation of AI. We present here an approach for assessing AI for predicting treatment response in triple-negative breast cancer (TNBC), using real-world data and molecular -omics data from clinical data warehouses and biobanks., Methods: The European "ITFoC (Information Technology for the Future Of Cancer)" consortium designed a framework for the clinical validation of AI technologies for predicting treatment response in oncology., Results: This framework is based on seven key steps specifying: (1) the intended use of AI, (2) the target population, (3) the timing of AI evaluation, (4) the datasets used for evaluation, (5) the procedures used for ensuring data safety (including data quality, privacy and security), (6) the metrics used for measuring performance, and (7) the procedures used to ensure that the AI is explainable. This framework forms the basis of a validation platform that we are building for the "ITFoC Challenge". This community-wide competition will make it possible to assess and compare AI algorithms for predicting the response to TNBC treatments with external real-world datasets., Conclusions: The predictive performance and safety of AI technologies must be assessed in a robust, unbiased and transparent manner before their implementation in healthcare settings. We believe that the consideration of the ITFoC consortium will contribute to the safe transfer and implementation of AI in clinical settings, in the context of precision oncology and personalized care., (© 2021. The Author(s).)

Background: Multiple sclerosis (MS) has a complex pathophysiology, variable clinical presentation, and unpredictable prognosis; understanding the underlying mechanisms requires combinatorial approaches that warrant the integration of diverse molecular omics data., Methods: Here, we combined genomic and proteomic data of the same individuals among a Turkish MS patient group to search for biologically important networks. We previously identified differentially-expressed proteins by cerebrospinal fluid proteome analysis of 179 MS patients and 42 non-MS controls. Among this study group, 11 unrelated MS patients and 60 independent, healthy controls were subjected to whole-genome SNP genotyping, and genome-wide associations were assessed. Pathway enrichment analyses of MS-associated SNPs and differentially-expressed proteins were conducted using the functional enrichment tool, PANOGA., Results: Nine shared pathways were detected between the genomic and proteomic datasets after merging and clustering the enriched pathways. Complement and coagulation cascade was the most significantly associated pathway (hsa04610, P  = 6.96 × 10 -30 ). Other pathways involved in neurological or immunological mechanisms included adherens junctions (hsa04520, P  = 6.64 × 10 -25 ), pathogenic Escherichia coli infection (hsa05130, P  = 9.03 × 10 -14 ), prion diseases (hsa05020, P  = 5.13 × 10 -13 )., Conclusion: We conclude that integrating multiple datasets of the same patients helps reducing false negative and positive results of genome-wide SNP associations and highlights the most prominent cellular players among the complex pathophysiological mechanisms., Competing Interests: The authors declare there are no competing interests., (©2021 Everest et al.)

Metallation status of human Cu/Zn superoxide dismutase 1 (SOD1) plays a pivotal role in the pathogenesis of amyotrophic lateral sclerosis (ALS). All of the amino acids found in the bimetallic center have been associated with ALS except for two positions. H63 which forms the bridging imidazolate ion in the bimetallic center and K136 which is not directly involved in coordination but located in the bimetallic center were not reported to be mutated in any of the identified ALS cases. In this study, we investigated the structure and flexibility of five SOD1 variants by using classical molecular dynamics simulations. These variants include three substitutions on the non-ALS-linked positions; H63A, H63R, K136A and ALS-linked positions; G37R, H46R/H48D. We have generated four systems for each variant differing in metallation and presence of the intramolecular disulfide bond. Overall, a total of 24 different dimers including the wild-type were generated and simulated at two temperatures, 298 and 400 K. We have monitored backbone mobility, fluctuations and compactness of the dimer structures to assess whether the hypothetical mutations would behave similar to the ALS-linked variants. Results showed that particularly two mutants, H63R and K136A, drastically affected the dimer dynamics by increasing the fluctuations of the metal binding loops compared with the control mutations. Further, these variants resulted in demetallation of the dimers, highlighting probable ALS toxicity that could be elicited by the SOD1 variants of H63R and K136A. Overall, this study bridges two putative SOD1 positions in the metallic center and ALS, underlining the potential use of atomistic simulations for studying disease variants., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Timucin, Cinaroglu, Sezerman and Timucin.)

We present a global atlas of 4,728 metagenomic samples from mass-transit systems in 60 cities over 3 years, representing the first systematic, worldwide catalog of the urban microbial ecosystem. This atlas provides an annotated, geospatial profile of microbial strains, functional characteristics, antimicrobial resistance (AMR) markers, and genetic elements, including 10,928 viruses, 1,302 bacteria, 2 archaea, and 838,532 CRISPR arrays not found in reference databases. We identified 4,246 known species of urban microorganisms and a consistent set of 31 species found in 97% of samples that were distinct from human commensal organisms. Profiles of AMR genes varied widely in type and density across cities. Cities showed distinct microbial taxonomic signatures that were driven by climate and geographic differences. These results constitute a high-resolution global metagenomic atlas that enables discovery of organisms and genes, highlights potential public health and forensic applications, and provides a culture-independent view of AMR burden in cities., Competing Interests: Declaration of interests C.E.M. is co-founder of Biotia and Onegevity Health. D.B. is co-founder and CSO of Poppy Health Inc. The other authors declare they have no competing interests that impacted this study., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

Background: Cancer develops due to "driver" alterations. Numerous approaches exist for predicting cancer drivers from cohort-scale genomics data. However, methods for personalized analysis of driver genes are underdeveloped. In this study, we developed a novel personalized/batch analysis approach for driver gene prioritization utilizing somatic genomics data, called driveR., Results: Combining genomics information and prior biological knowledge, driveR accurately prioritizes cancer driver genes via a multi-task learning model. Testing on 28 different datasets, this study demonstrates that driveR performs adequately, achieving a median AUC of 0.684 (range 0.651-0.861) on the 28 batch analysis test datasets, and a median AUC of 0.773 (range 0-1) on the 5157 personalized analysis test samples. Moreover, it outperforms existing approaches, achieving a significantly higher median AUC than all of MutSigCV (Wilcoxon rank-sum test p < 0.001), DriverNet (p < 0.001), OncodriveFML (p < 0.001) and MutPanning (p < 0.001) on batch analysis test datasets, and a significantly higher median AUC than DawnRank (p < 0.001) and PRODIGY (p < 0.001) on personalized analysis datasets., Conclusions: This study demonstrates that the proposed method is an accurate and easy-to-utilize approach for prioritizing driver genes in cancer genomes in personalized or batch analyses. driveR is available on CRAN: https://cran.r-project.org/package=driveR .

Glycoprotein M6A, a stress related gene, plays an important role in synapse and filopodia formation. Filopodia formation is vital for development, immunity, angiogenesis, wound healing and metastasis. In this study, structural and functional analysis of high-risk SNPs associated with Glycoprotein M6-A were evaluated using six different bioinformatics tools. Results classified T210I, T134I, Y153H, I215T, F156L, T160I, I226T, R247W, R178C, W159R, N157S and P151L as deleterious mutants that are crucial for the structure and function of the protein causing malfunction of M6-a and ultimately leads to disease development. The three-dimensional structure of wild-type M6-a and mutant M6-a were also predicted. Furthermore, the effects of high risk substitutions were also analyzed with interaction with valproic acid. Based on structural models obtained, the binding pocket of ligand bound glycoprotein M6-A structure showed few core interacting residues which are different in the mutant models. Among all substitutions, F156L showed complete loss of binding pocket when interacting with valproic acid as compared to the wild type model. Up to the best of our knowledge this is the first comprehensive study where GPM6A mutations were analyzed. The mechanism of action of GPM6A is still not fully defined which limits the understanding of functional details encoding M6-A. Our results may help enlighten some molecular aspects underlying glycoprotein M6-A. Communicated by Ramaswamy H. Sarma.

Background: In the clinical setting, workflows for analyzing individual genomics data should be both comprehensive and convenient for clinical interpretation. In an effort for comprehensiveness and practicality, we attempted to create a clinical individual whole exome sequencing (WES) analysis workflow, allowing identification of genomic alterations and presentation of neurooncologically-relevant findings., Methods: The analysis workflow detects germline and somatic variants and presents: (1) germline variants, (2) somatic short variants, (3) tumor mutational burden (TMB), (4) microsatellite instability (MSI), (5) somatic copy number alterations (SCNA), (6) SCNA burden, (7) loss of heterozygosity, (8) genes with double-hit, (9) mutational signatures, and (10) pathway enrichment analyses. Using the workflow, 58 WES analyses from matched blood and tumor samples of 52 patients were analyzed: 47 primary and 11 recurrent diffuse gliomas., Results: The median mean read depths were 199.88 for tumor and 110.955 for normal samples. For germline variants, a median of 22 (14-33) variants per patient was reported. There was a median of 6 (0-590) reported somatic short variants per tumor. A median of 19 (0-94) broad SCNAs and a median of 6 (0-12) gene-level SCNAs were reported per tumor. The gene with the most frequent somatic short variants was TP53 (41.38%). The most frequent chromosome-/arm-level SCNA events were chr7 amplification, chr22q loss, and chr10 loss. TMB in primary gliomas were significantly lower than in recurrent tumors (p = 0.002). MSI incidence was low (6.9%)., Conclusions: We demonstrate that WES can be practically and efficiently utilized for clinical analysis of individual brain tumors. The results display that NOTATES produces clinically relevant results in a concise but exhaustive manner.

Little is known about the mutational processes that shape the genetic landscape of gliomas. Numerous mutational processes leave marks on the genome in the form of mutations, copy number alterations, rearrangements or their combinations. To explore gliomagenesis, we hypothesized that gliomas with different underlying oncogenic mechanisms would have differences in the burden of various forms of these genomic alterations. This was an analysis on adult diffuse gliomas, but IDH -mutant gliomas as well as diffuse midline gliomas H3 -K27M were excluded to search for the possible presence of new entities among the very heterogenous group of IDH -WT glioblastomas. The cohort was divided into two molecular subsets: (1) Molecularly-defined GBM (mGBM) as those that carried molecular features of glioblastomas (including TERT promoter mutations, 7/10 pattern, or EGFR -amplification), and (2) those who did not (others). Whole exome sequencing was performed for 37 primary tumors and matched blood samples as well as 8 recurrences. Single nucleotide variations (SNV), short insertion or deletions (indels) and copy number alterations (CNA) were quantified using 5 quantitative metrics (SNV burden, indel burden, copy number alteration frequency-wGII, chromosomal arm event ratio-CAER, copy number amplitude) as well as 4 parameters that explored underlying oncogenic mechanisms (chromothripsis, double minutes, microsatellite instability and mutational signatures). Findings were validated in the TCGA pan-glioma cohort. mGBM and "Others" differed significantly in their SNV (only in the TCGA cohort) and CNA metrics but not indel burden. SNV burden increased with increasing age at diagnosis and at recurrences and was driven by mismatch repair deficiency. On the contrary, indel and CNA metrics remained stable over increasing age at diagnosis and with recurrences. Copy number alteration frequency (wGII) correlated significantly with chromothripsis while CAER and CN amplitude correlated significantly with the presence of double minutes, suggesting separate underlying mechanisms for different forms of CNA.

Formation of atherosclerotic plaques, called atherogenesis, is a complex process affected by genetic and environmental factors. It was proposed that endoplasmic reticulum (ER) stress is an important factor in the pathogenesis of atherosclerosis and that vitamin E affects atherosclerotic plaque formation via its antioxidant properties. Here, we investigated ER stress-related molecular mechanisms in high-cholesterol diet (HCD, 2%)-induced atherosclerosis model and the role of vitamin E supplementation in it, beyond its antioxidant properties. The consequences of HCD and vitamin E supplementation were examined by determining protein levels of ER stress markers in aortic tissues. As vitamin E supplementation acts on several unfolded protein response (UPR) factors, it decreased ER stress induced by HCD. To elucidate the associated pathways, gene expression profiling was performed, revealing differentially expressed genes enriched in ER stress-related pathways such as the proteasome and the apoptosis pathways. We further assessed the proteasomal activity impaired by HCD in the aorta and showed that vitamin E reversed it to that of control animals. Overall, the study characterized the effects of HCD and vitamin E on ER stress-related gene expression, revealing the role of proteolytic systems during atherogenesis., (© 2020 International Union of Biochemistry and Molecular Biology.)

The novel coronavirus pneumonia, which was named later as coronavirus disease 2019 (COVID-19), is caused by the severe acute respiratory syndrome coronavirus 2, namely SARS-CoV-2. It is a positive-strand RNA virus that is the seventh coronavirus known to infect humans. The COVID-19 outbreak presents enormous challenges for global health behind the pandemic outbreak. The first diagnosed patient in Turkey has been reported by the Republic of Turkey Ministry of Health on March 11, 2020. In May, over 150,000 cases in Turkey, and 5.5 million cases around the world have been declared. Due to the urgent need for a vaccine and antiviral drug, isolation of the virus is crucial. Here, we report 1 of the first isolation and characterization studies of SARS-CoV-2 from nasopharyngeal and oropharyngeal specimens of diagnosed patients in Turkey. This study provides an isolation and replication methodology,and cell culture tropism of the virus that will be available to the research communities., Competing Interests: CONFLICT OF INTEREST: none declared, (Copyright © 2020 The Author(s).)

In the nitric oxide (NO) signaling pathway, human soluble guanylate cyclase (hsGC) synthesizes cyclic guanosine monophosphate (cGMP); responsible for the regulation of cGMP-specific protein kinases (PKGs) and phosphodiesterases (PDEs). The crystal structure of the inactive hsGC cyclase dimer is known, but there is still a lack of information regarding the substrate-specific internal motions that are essential for the catalytic mechanism of the hsGC. In the current study, the hsGC cyclase heterodimer complexed with guanosine triphosphate (GTP) and cGMP was subjected to molecular dynamics simulations, to investigate the conformational dynamics that have functional implications on the catalytic activity of hsGC. Results revealed that in the GTP-bound complex of the hsGC heterodimer, helix 1 of subunit α (α:h1) moves slightly inwards and comes close to helix 4 of subunit β (β:h4). This conformational change brings loop 2 of subunit β (β:L2) closer to helix 2 of subunit α (α:h2). Likewise, loop 2 of subunit α (α:L2) comes closer to helix 2 of subunit β (β:h2). These structural events stabilize and lock GTP within the closed pocket for cyclization. In the cGMP-bound complex, α:L2 detaches from β:h2 and establishes interactions with β:L2, which results in the loss of global structure compactness. Furthermore, with the release of pyrophosphate, the interaction between α:h1 and β:L2 weakens, abolishing the tight packing of the binding pocket. This study discusses the conformational changes induced by the binding of GTP and cGMP to the hsGC catalytic domain, valuable in designing new therapeutic strategies for the treatment of cardiovascular diseases.

Bacillus thermocatenulatus lipase (BTL2), a member of the isolated lipase family known as thermoalkalophilic lipases, carries potential for industrial applications owing to its ability to catalyze versatile reactions under extreme conditions. This study investigates the molecular effects of distinct solvents on the stability of BTL2 at different temperatures, aiming to contribute to lipase use in industrial applications. Initially, molecular dynamic (MD) simulations were carried out to address for the molecular impacts of distinct solvents on the structural stability of BTL2 at different temperatures. Two lipase conformations representing the active and inactive forms were simulated in 5 solvents including water, ethanol, methanol, cyclohexane, and toluene. Low temperature simulations showed that polar solvents led to enhanced lid fluctuations compared with non-polar solvents reflecting a more dynamic equilibrium between active and inactive lipase conformations in polar solvents including water, while the overall structure of the lipase in both forms became more rigid in non-polar solvents than they were in polar solvent. Notably, the native lipase fold was maintained in non-polar solvents even at high temperatures, indicating an enhancement of lipase's thermostability in non-polar organic solvents. Next, we conducted experiments for which BTL2 was expressed in a heterologous host and purified to homogeneity, and its thermostability in different solvents was assessed. Parallel to the computational findings, experimental results suggested that non-polar organic solvents contributed to BTL2's thermostability at concentrations as high as 70% (v/v). Altogether, this study provides beneficial insights to the lipase use under extreme conditions. Graphical Abstract.

Precision oncology relies on accurate discovery and interpretation of genomic variants, enabling individualized diagnosis, prognosis and therapy selection. We found that six prominent somatic cancer variant knowledgebases were highly disparate in content, structure and supporting primary literature, impeding consensus when evaluating variants and their relevance in a clinical setting. We developed a framework for harmonizing variant interpretations to produce a meta-knowledgebase of 12,856 aggregate interpretations. We demonstrated large gains in overlap between resources across variants, diseases and drugs as a result of this harmonization. We subsequently demonstrated improved matching between a patient cohort and harmonized interpretations of potential clinical significance, observing an increase from an average of 33% per individual knowledgebase to 57% in aggregate. Our analyses illuminate the need for open, interoperable sharing of variant interpretation data. We also provide a freely available web interface (search.cancervariants.org) for exploring the harmonized interpretations from these six knowledgebases.

Tyrosine kinase inhibitor (TKI) resistance is a major problem in chronic myeloid leukemia (CML). We generated a TKI-resistant K562 sub-population, K562-IR, under selective imatinib-mesylate pressure. K562-IR cells are CD34-/CD38-, BCR-Abl-independent, proliferate slowly, highly adherent and form intact tumor spheroids. Loss of CD45 and other hematopoietic markers reveal these cells have diverged from their hematopoietic origin. CD34 negativity, high expression of E-cadherin and CD44; decreased levels of CD45 and β-catenin do not fully confer with the leukemic stem cell (LSC) phenotype. Expression analyses reveal that K562-IR cells differentially express tissue/organ development and differentiation genes. Our data suggest that the observed phenotypic shift is an adaptive process rendering cells under TKI stress to become oncogene independent. Cells develop transcriptional instability in search for a gene expression framework suitable for new environmental stresses, resulting in an adaptive phenotypic shift in which some cells partially display LSC-like properties. With leukemic/cancer stem cell targeted therapies underway, the difference between treating an entity and a spectrum of dynamic cellular states will have conclusive effects on the outcome., Competing Interests: The authors have declared that no competing interests exist.

Pathway analysis is often the first choice for studying the mechanisms underlying a phenotype. However, conventional methods for pathway analysis do not take into account complex protein-protein interaction information, resulting in incomplete conclusions. Previously, numerous approaches that utilize protein-protein interaction information to enhance pathway analysis yielded superior results compared to conventional methods. Hereby, we present pathfindR, another approach exploiting protein-protein interaction information and the first R package for active-subnetwork-oriented pathway enrichment analyses for class comparison omics experiments. Using the list of genes obtained from an omics experiment comparing two groups of samples, pathfindR identifies active subnetworks in a protein-protein interaction network. It then performs pathway enrichment analyses on these identified subnetworks. To further reduce the complexity, it provides functionality for clustering the resulting pathways. Moreover, through a scoring function, the overall activity of each pathway in each sample can be estimated. We illustrate the capabilities of our pathway analysis method on three gene expression datasets and compare our results with those obtained from three popular pathway analysis tools. The results demonstrate that literature-supported disease-related pathways ranked higher in our approach compared to the others. Moreover, pathfindR identified additional pathways relevant to the conditions that were not identified by other tools, including pathways named after the conditions., (Copyright © 2019 Ulgen, Ozisik and Sezerman.)

Objective: Processes that cause or contribute to cancer, such as aging, exposure to carcinogens, or DNA damage repair deficiency (DDRd), create predictable and traceable nucleotide alterations in one's genetic code (termed "mutational signatures"). Large studies have previously identified various such mutational signatures across cancers that can be attributed to the specific causative processes. To gain further insight into the processes in glioma development, the authors analyzed mutational signatures in adult diffuse gliomas (DGs)., Methods: Twenty-five DGs and paired blood samples were whole exome sequenced. Somatic mutational signatures were identified using 2 different methods. Associations of the signatures with age at diagnosis, molecular subset, and mutational load were investigated. As DDRd-related signatures were frequently observed, germline and somatic DDR gene mutations as well as microsatellite instability (MSI) status were determined for all samples. For validation of signature prevalence, publicly available data from The Cancer Genome Atlas (TCGA) were used., Results: Each tumor had a unique combination of signatures. The most common signatures were signature 1 (88%, aging related), signature 3 (52%, homologous recombination related), and signature 15 (56%, mismatch repair related). Eighty-four percent of the tumors contained at least 1 DDRd signature. The findings were validated using public TCGA data. The weight of signature 1 positively correlated with age (r = 0.43) while cumulative weight of DDRd signatures negatively correlated with age (r = -0.16). Each subject had at least 1 germline/somatic alteration in a DDR gene, the most common being the risk single nucleotide polymorphism rs1800734 in MLH1. The rs1800734-AA genotype had a higher cumulative DDRd weight as well as higher mutational load; TP53 was the most common somatically altered DDR gene. MSI was observed in 24% of the tumors. No significant associations of MSI status with mutational load, rs1800734, or the cumulative weight of DDRd signatures were identified., Conclusions: Current findings suggest that DDRd may act as a fundamental mechanism in gliomagenesis rather than being a random, secondary event.

Antigen delivery platforms based on engineered viruses or virus-like particles are currently developed as vaccines against infectious diseases. As the interaction of vaccines with dendritic cells (DCs) shapes the immunological response, we compared the interaction of a range of virus-based vectors and virus-like particles with DCs in a murine model of systemic administration and transcriptome analyses of splenic DCs. The transcriptome profiles of DCs separated the vaccine vectors into two distinct groups characterized by high- and low-magnitude differential gene expression, which strongly correlated with (1) the surface expression of costimulatory molecules CD40, CD83, and CD86 on DCs, and (2) antigen-specific T-cell responses. Pathway analysis using PANOGA (Pathway and Network-Oriented GWAS Analysis) revealed that the JAK/STAT pathway was significantly activated by both groups of vaccines. In contrast, the oxidative phosphorylation pathway was significantly downregulated only by the high-magnitude DC-stimulating vectors. A gene signature including exclusively chemokine-, cytokine-, and receptor-related genes revealed a vector-specific pattern. Overall, this in vivo DC stimulation model demonstrated a strong relationship between the levels of induced DC maturation and the intensity of T-cell-specific immune responses with a distinct cytokine/chemokine profile, metabolic shifting, and cell surface expression of maturation markers. It could represent an important tool for vaccine design.

Genetic generalized epilepsies (GGE) are genetically determined, as their name implies and they are clinically characterized by generalized seizures involving both sides of the brain in the absence of detectable brain lesions or other known causes. GGEs are yet complex and are influenced by many different genetic and environmental factors. Methylation specific epigenetic marks are one of the players of the complex epileptogenesis scenario leading to GGE. In this study, we have set out to perform genome-wide methylation profiling to analyze GGE trios each consisting of an affected parent-offspring couple along with an unaffected parent. We have developed a novel scoring scheme within trios to categorize each locus analyzed as hypo or hypermethylated. This stringent approach classified differentially methylated genes in each trio and helped us to produce trio specific and pooled gene lists with inherited and aberrant methylation levels. In order to analyze the methylation differences from a boarder perspective, we performed enrichment analysis with these lists using the PANOGA software. This collective effort has led us to detect pathways associated with the GGE phenotype, including the neurotrophin signaling pathway. We have demonstrated a trio based approach to genome-wide DNA methylation analysis that identified individual and possibly minor changes in methylation marks that could be involved in epileptogenesis leading to GGE., Competing Interests: The authors have declared that no competing interests exist.

Leishmania virus (LRV) has previously been identified in different Leishmania species. Host-LRV interaction is associated with exacerbated clinical manifestations of cutaneous leishmaniasis (CL) and may cause poor therapeutic response. CL cases due to L. major with large skin lesions resistant to routine therapy were recently identified in Turkey. Here, we report the first autochthonous case of cutaneous leishmaniasis caused by LRV-positive Leishmania major, using conventional PCR targeting the viral capsid protein of LRV. The lesion of the case was 6 months old, relatively large (4 cm), and did not recover despite three consecutive intralesional applications of glucantime. Assessment of LRV's influence on prognosis and clinical outcomes of leishmaniasis, based on additional studies, is required.

The Drug often binds to more than one targets defined as polypharmacology, one application of which is drug repurposing also referred as drug repositioning or therapeutic switching. The traditional drug discovery and development is a high-priced and tedious process, thus making drug repurposing a popular alternate strategy. We proposed an integrative method based on similarity scheme that predicts approved and novel Drug targets with new disease associations. We combined PPI, biological pathways, binding site structural similarities and disease-disease similarity measures. The results showed 94% Accuracy with 0.93 Recall and 0.94 Precision measure in predicting the approved and novel targets surpassing the existing methods. All these parameters help in elucidating the unknown associations between drug and diseases for finding the new uses for old drugs., (Copyright © 2018 Elsevier Inc. All rights reserved.)

The discovery or development of thermoalkalophilic lipases that show high levels of catalytic activity in organic solvents would have important industrial ramifications. However, this goal is yet to be achieved because organic solvents induce structural changes in lipases that suppress their catalytic abilities. A deep understanding of these structural changes to lipases in the presence of organic solvents is required before strategies can be devised to stop them from occurring. In this work, we investigated the effects of an organic reaction medium, toluene, on the structure of the Bacillus thermocatenulatus lipase BTL2 using MD simulation. The main aims were to identify the regions of the protein that are particularly sensitive to the presence of an organic solvent, and how the presence of a hydrophobic medium affects the overall stability of the enzyme. Upon analyzing how the behavior of the enzyme differed in aqueous and hydrophobic media, it was found that many significant zones of the protein suffer in the presence of an organic solvent, which increases the rigidity of the system. This was readily apparent when we investigated important noncovalent interactions (salt bridges) and probed how distances between the atoms of the catalytic triad Ser114, Asp318, and His359 change in the presence of toluene. Moreover, the high tendency for the system to destabilize in toluene was explained by the results of FoldX calculations. Calculations showed that the addition of a small amount of water to the hydrophobic reaction environment should restore the required flexibility of BTL2. The insights gained from the analysis of our simulations allowed us to propose a modification of BTL2, the G116P mutation, that should result in the structural behavior of BTL2 in organic solvent being closer to that of BTL2 in water.

Drug resistance is a major obstacle faced by therapist in treating HIV infected patients. The reason behind these phenomena is either protein mutation or the changes in gene expression level that induces resistance to drug treatments. These mutations affect the drug binding activity, hence resulting in failure of treatment. Therefore, it is necessary to conduct resistance testing in order to carry out HIV effective therapy. This study combines both sequence and structural features for predicting HIV resistance by applying SVM and Random Forests classifiers. The model was tested on the mutants of HIV-1 protease and reverse transcriptase. Taken together the features we have used in our method, total contact energies among multiple mutations have a strong impact in predicting resistance as they are crucial in understanding the interactions of HIV mutants. The combination of sequence-structure features offers high accuracy with support vector machines as compared to Random Forests classifier. Both single and acquisition of multiple mutations are important in predicting HIV resistance to certain drug treatments. We have discovered the practicality of these features; hence, these can be used in the future to predict resistance for other complex diseases.

The known challenge of underutilization of data and biological material from biorepositories as potential resources for medical research has been the focus of discussion for over a decade. Recently developed guidelines for improved data availability and reusability-entitled FAIR Principles (Findability, Accessibility, Interoperability, and Reusability)-are likely to address only parts of the problem. In this article, we argue that biological material and data should be viewed as a unified resource. This approach would facilitate access to complete provenance information, which is a prerequisite for reproducibility and meaningful integration of the data. A unified view also allows for optimization of long-term storage strategies, as demonstrated in the case of biobanks. We propose an extension of the FAIR Principles to include the following additional components: (1) quality aspects related to research reproducibility and meaningful reuse of the data, (2) incentives to stimulate effective enrichment of data sets and biological material collections and its reuse on all levels, and (3) privacy-respecting approaches for working with the human material and data. These FAIR-Health principles should then be applied to both the biological material and data. We also propose the development of common guidelines for cloud architectures, due to the unprecedented growth of volume and breadth of medical data generation, as well as the associated need to process the data efficiently.

Magnesium deficiency increases the susceptibility of plants toward heat stress. The correlation between magnesium levels and stress response has been studied at the physiological level; albeit, the molecular explanation to this relationship remains elusive. Among diverse pathways implicated in the heat stress, the abscisic acid (ABA) signal modulates the heat stress response by magnesium dependent phosphatases (PP2Cs). Exclusively, sequestration of PP2Cs by ABA receptors (PYLs) in the heterodimer form activates the stress response through ABA responsive transcription factors. In this study, the molecular interplay between magnesium levels and ABA related heat stress response was investigated. Molecular dynamics simulations have been applied to two different PYL-PP2C heterodimer systems representing normal and magnesium deficient conditions. The heterodimer conformation and stability were delineated at high temperatures mimicking heat stress. Results showed that the thermostability of the heat stress response heterodimer was significantly dependent on the magnesium. Furthermore, a conserved aromatic cluster at the dimer interface acted synergistically with the metal to confer thermostability to the heterodimer structure. These structural insights into one of the possible links between magnesium levels and stress highlight the importance of metal micronutrients for tuning the stability of the stress-related proteins and optimizing tolerance.

Extracting useful knowledge from an unstructured textual data is a challenging task for biologists, since biomedical literature is growing exponentially on a daily basis. Building an automated method for such tasks is gaining much attention of researchers. ZK DrugResist is an online tool that automatically extracts mutations and expression changes associated with drug resistance from PubMed. In this study we have extended our tool to include semantic relations extracted from biomedical text covering drug resistance and established a server including both of these features. Our system was tested for three relations, Resistance (R), Intermediate (I) and Susceptible (S) by applying hybrid feature set. From the last few decades the focus has changed to hybrid approaches as it provides better results. In our case this approach combines rule-based methods with machine learning techniques. The results showed 97.67% accuracy with 96% precision, recall and F-measure. The results have outperformed the previously existing relation extraction systems thus can facilitate computational analysis of drug resistance against complex diseases and further can be implemented on other areas of biomedicine., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Metal ions play pivotal roles in protein structure, function and stability. The functional and structural diversity of proteins in nature expanded with the incorporation of metal ions or clusters in proteins. Approximately one-third of these proteins in the databases contain metal ions. Many biological and chemical processes in nature involve metal ion-binding proteins, aka metalloproteins. Many cellular reactions that underpin life require metalloproteins. Most of the remarkable, complex chemical transformations are catalysed by metalloenzymes. Realization of the importance of metal-binding sites in a variety of cellular events led to the advancement of various computational methods for their prediction and characterization. Furthermore, as structural and functional knowledgebase about metalloproteins is expanding with advances in computational and experimental fields, the focus of the research is now shifting towards de novo design and redesign of metalloproteins to extend nature's own diversity beyond its limits. In this review, we will focus on the computational toolbox for prediction of metal ion-binding sites, de novo metalloprotein design and redesign. We will also give examples of tailor-made artificial metalloproteins designed with the computational toolbox., (© 2017 The Author(s).)

Allatostatins (AST) are neuropeptides with variable function ranging from regulation of developmental processes to the feeding behavior in insects. They exert their effects by binding to cognate GPCRs, called Allatostatin receptors (AlstR), which emerge as promising targets for pesticide design. However, AlstRs are rarely studied. This study is the first reported structural study on AlstR-AST interaction. In this work, the first C type AlstR from the stick insect Carausius morosus (CamAlstR-C) was identified and its interaction with type C AST peptide was shown to be physically consistent with the experimental results. The proposed structure of CamAlstR-C revealed a conserved motif within the third extracellular loop, which, together with the N-terminus is essential for ligand binding. In this work, computational studies were combined with molecular and nano-scale approaches in order to introduce an unknown GPCR-ligand system. Consequently, the data obtained provided a reliable target region for future agonist/inverse agonist studies on AlstRs.

Autophagy is biological mechanism allowing recycling of long-lived proteins, abnormal protein aggregates, and damaged organelles under cellular stress conditions. Following sequestration in double- or multimembrane autophagic vesicles, the cargo is delivered to lysosomes for degradation. ATG5 is a key component of an E3-like ATG12-ATG5-ATG16 protein complex that catalyzes conjugation of the MAP1LC3 protein to lipids, thus controlling autophagic vesicle formation and expansion. Accumulating data indicate that ATG5 is a convergence point for autophagy regulation. Here, we describe the scaffold protein RACK1 (receptor activated C-kinase 1, GNB2L1) as a novel ATG5 interactor and an autophagy protein. Using several independent techniques, we showed that RACK1 interacted with ATG5. Importantly, classical autophagy inducers (starvation or mammalian target of rapamycin blockage) stimulated RACK1-ATG5 interaction. Knockdown of RACK1 or prevention of its binding to ATG5 using mutagenesis blocked autophagy activation. Therefore, the scaffold protein RACK1 is a new ATG5-interacting protein and an important and novel component of the autophagy pathways., (© 2016 by The American Society for Biochemistry and Molecular Biology, Inc.)

The single nucleotide polymorphism rs55705857, located in a non-coding but evolutionarily conserved region at 8q24.21, is strongly associated with IDH-mutant glioma development and was suggested to be a causal variant. However, the molecular mechanism underlying this association has remained unknown. With a case control study in 285 gliomas, 316 healthy controls, 380 systemic cancers, 31 other CNS-tumors, and 120 IDH-mutant cartilaginous tumors, we identified that the association was specific to IDH-mutant gliomas. Odds-ratios were 9.25 (5.17-16.52; 95% CI) for IDH-mutated gliomas and 12.85 (5.94-27.83; 95% CI) for IDH-mutated, 1p/19q co-deleted gliomas. Decreasing strength with increasing anaplasia implied a modulatory effect. No somatic mutations were noted at this locus in 114 blood-tumor pairs, nor was there a copy number difference between risk-allele and only-ancestral allele carriers. CCDC26 RNA-expression was rare and not different between the two groups. There were only minor subtype-specific differences in common glioma driver genes. RNA sequencing and LC-MS/MS comparisons pointed to significantly altered MYC-signaling. Baseline enhancer activity of the conserved region specifically on the MYC promoter and its further positive modulation by the SNP risk-allele was shown in vitro. Our findings implicate MYC deregulation as the underlying cause of the observed association.

Cellulases have great potential to be widely used for industrial applications. In general, naturally occurring cellulases are not optimized and limited to meet the industrial needs. These limitations lead to demand for novel cellulases with enhanced enzymatic properties. Here, we describe the enzymatic and structural properties of two novel enzymes, EG3&#95;S1 and EG3&#95;S2, obtained through the single-gene shuffling approach of Cel12A(EG3) gene from Trichoderma reseei EG3&#95;S1 and EG3&#95;S2 shuffled enzymes display 59 and 75% identity in protein sequence with respect to native, respectively. Toward 4-MUC, the minimum activity of EG3&#95;S1 was reported as 5.9-fold decrease in native at 35°C, whereas the maximum activity of EG3&#95;S2 was reported as 15.4-fold increase in native activity at 40°C. Also, the diminished enzyme activity of EG3&#95;S1 was reported within range of 0.6- to 0.8-fold of native and within range of 0.5- to 0.7-fold of native toward CMC and Na-CMC, respectively. For EG3&#95;S2 enzyme, the improved enzymatic activities within range of 1.1- to 1.4-fold of native and within range of 1.1- to 1.6-fold of native were reported toward CMC and Na-CMC, respectively. Moreover, we have reported 6.5-fold increase in the kcat/Km ratio of EG3&#95;S2 with respect to native and suggested EG3&#95;S2 enzyme as more efficient catalysis for hydrolysis reactions than its native counterpart., (© The Author 2016. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)