Longitudinal non-targeted metabolomic profiling of urine samples for monitoring of kidney transplantation patients.

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 1st week, and seventeen metabolites in the 3rd 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 6th 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 3rd-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.

Effect of the switch status of Helicobacter pylori outer inflammatory protein A on gastric diseases.

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.

Phenotype screens of murine pancreatic cancer identify a Tgf-α-Ccl2-paxillin axis driving human-like neural invasion.

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.

Corneal bacterial microbiome in patients with keratoconus using next-generation sequencing-based 16S rRNA gene analysis.

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.

Metabolomic Analyses to Identify Candidate Biomarkers of Cystinosis.

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.

PANACEA: network-based methods for pharmacotherapy prioritization in personalized oncology.

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.

Integrated bioinformatics analysis of validated and circulating miRNAs in ovarian cancer.

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.

Computational analysis of missense filamin-A variants, including the novel p.Arg484Gln variant of two brothers with periventricular nodular heterotopia.

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.

Is cutaneous microbiota a player in disease pathogenesis? Comparison of cutaneous microbiota in psoriasis and seborrheic dermatitis with scalp involvement.

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.

Structural analysis of M1AP variants associated with severely impaired spermatogenesis causing male infertility.

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.

A splice site mutation in the TSEN2 causes a new syndrome with craniofacial and central nervous system malformations, and atypical hemolytic uremic syndrome.

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).

Intraoperative Neuromonitoring of Blink Reflex During Posterior Fossa Surgeries and Its Correlation With Clinical Outcome.

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.

A recursive molecular docking coupled with energy-based pose-rescoring and MD simulations to identify hsGC ßH-NOX allosteric modulators for cardiovascular dysfunctions.

Modulating the activity of human soluble guanylate cyclase (hsGC) 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 hsGC 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 hsGC as compared to the reference compound.

Identifying and elucidating the roles of Y198N and Y204F mutations in the PAH enzyme through molecular dynamic simulations.

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.

Bridging the Bridging Imidazolate in the Bimetallic Center of the Cu/Zn SOD1 and ALS.

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.

A comprehensive study on identifying the structural and functional SNPs of human neuronal membrane glycoprotein M6A (GPM6A).

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.

Probing the Structural Dynamics of the Catalytic Domain of Human Soluble Guanylate Cyclase.

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.

SARS-CoV-2 isolation and propagation from Turkish COVID-19 patients.

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.

Understanding thermal and organic solvent stability of thermoalkalophilic lipases: insights from computational predictions and experiments.

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.

High-throughput profiling reveals perturbation of endoplasmic reticulum stress-related genes in atherosclerosis induced by high-cholesterol diet and the protective role of vitamin E.

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.