UniProt: the universal protein knowledgebase in 2021

Integrated Genomic Analysis Identifies ANKRD36 Gene as a Novel and Common Biomarker of Disease Progression in Chronic Myeloid Leukemia

Simple Summary

Chronic myeloid leukemia is a type of blood cancer that is regarded as a success story in determining the exact biological origin, pathogenesis and development of a molecularly targeted (mutation-specific) therapy that has led to successful treatment of this fatal cancer. It is caused by the BCR-ABL fusion gene, which is formed from the translocation between chromosomes 9 and 22. Anti-BCR-ABL drugs, known as tyrosine kinase inhibitors (TKIs), have led to long-term remissions in more than 80% of CML patients and even cure in about one-third of patients. Nevertheless, many patients face drug resistance, and disease progression occurs in about 30% of CML patients, leading to morbidities and mortality. Unfortunately, no biomarkers of CML progression are available due to a poor understanding of the mechanism of progression. Therefore, finding reliable molecular biomarkers of CML progression is one of the most attractive research areas in 21st-century cancer research. In this study, we report novel genomic variants exclusively found in all our advanced-phase CML patients. This study will help in identifying CML patients at risk of disease progression and timely therapeutic interventions to avoid or at least delay fatal disease progression in this cancer.

Abstract

Background: Chronic myeloid leukemia (CML) is initiated in bone marrow due to chromosomal translocation t(9;22) leading to fusion oncogene BCR-ABL. Targeting BCR-ABL by tyrosine kinase inhibitors (TKIs) has changed fatal CML into an almost curable disease. Despite that, TKIs lose their effectiveness due to disease progression. Unfortunately, the mechanism of CML progression is poorly understood and common biomarkers for CML progression are unavailable. This study was conducted to find novel biomarkers of CML progression by employing whole-exome sequencing (WES). Materials and Methods: WES of accelerated phase (AP) and blast crisis (BC) CML patients was carried out, with chronic-phase CML (CP-CML) patients as control. After DNA library preparation and exome enrichment, clustering and sequencing were carried out using Illumina platforms. Statistical analysis was carried out using SAS/STAT software version 9.4, and R package was employed to find mutations shared exclusively by all AP-/BC-CML patients. Confirmation of mutations was carried out using Sanger sequencing and protein structure modeling using I-TASSER followed by mutant generation and visualization using PyMOL. Results: Three novel genes (ANKRD36, ANKRD36B and PRSS3) were mutated exclusively in all AP-/BC-CML patients. Only ANKRD36 gene mutations (c.1183_1184 delGC and c.1187_1185 dupTT) were confirmed by Sanger sequencing. Protein modeling studies showed that mutations induce structural changes in ANKRD36 protein. Conclusions: Our studies show that ANKRD36 is a potential common biomarker and drug target of early CML progression. ANKRD36 is yet uncharacterized in humans. It has the highest expression in bone marrow, specifically myeloid cells. We recommend carrying out further studies to explore the role of ANKRD36 in the biology and progression of CML.

Verification of the Potential Targets of the Herbal Prescription Sochehwan for Drug Repurposing Processes as Deduced by Network Pharmacology

Network pharmacology (NP) is a useful, emerging means of understanding the complex pharmacological mechanisms of traditional herbal medicines. Sochehwan (SCH) is a candidate herbal prescription for drug repurposing as it has been suggested to have beneficial effects on metabolic syndrome. In this study, NP was adopted to complement the shortcomings of literature-based drug repurposing strategies in traditional herbal medicine. We conducted in vitro studies to confirm the effects of SCH on potential pharmacological targets identified by NP analysis. Herbal compounds and molecular targets of SCH were explored and screened from a traditional Chinese medicine systems pharmacology database and analysis platform (TCMSP) and an oriental medicine advanced searching integrated system (OASIS). Forty-seven key targets selected from a protein-protein interaction (PPI) network were analyzed with gene ontology (GO) term enrichment and KEGG pathway enrichment analysis to identify relevant categories. The tumor necrosis factor (TNF) and mitogen-activated protein kinase (MAPK) signaling pathways were presented as significant signaling pathways with lowest p-values by NP analysis, which were downregulated by SCH treatment. The signal transducer and activator of transcription 3 (STAT3) was identified as a core key target by NP analysis, and its phosphorylation ratio was confirmed to be significantly suppressed by SCH. In conclusion, the NP-based approach used for target prediction and experimental data obtained from Raw 264.7 cells strongly suggested that SCH can attenuate inflammatory status by modulating the phosphorylation status of STAT3.

PDBsum extras: SARS-CoV-2 and AlphaFold models

The PDBsum web server provides structural analyses of the entries in the Protein Data Bank (PDB). Two recent additions are described here. The first is the detailed analysis of the SARS-CoV-2 virus protein structures in the PDB. These include the variants of concern, which are shown both on the sequences and 3D structures of the proteins. The second addition is the inclusion of the available AlphaFold models for human proteins. The pages allow a search of the protein against existing structures in the PDB via the Sequence Annotated by Structure (SAS) server, so one can easily compare the predicted model against experimentally determined structures. The server is freely accessible to all at http://www.ebi.ac.uk/pdbsum.

Ensembl Genomes 2022: an expanding genome resource for non-vertebrates

Ensembl Genomes (https://www.ensemblgenomes.org) provides access to non-vertebrate genomes and analysis complementing vertebrate resources developed by the Ensembl project (https://www.ensembl.org). The two resources collectively present genome annotation through a consistent set of interfaces spanning the tree of life presenting genome sequence, annotation, variation, transcriptomic data and comparative analysis. Here, we present our largest increase in plant, metazoan and fungal genomes since the project's inception creating one of the world's most comprehensive genomic resources and describe our efforts to reduce genome redundancy in our Bacteria portal. We detail our new efforts in gene annotation, our emerging support for pangenome analysis, our efforts to accelerate data dissemination through the Ensembl Rapid Release resource and our new AlphaFold visualization. Finally, we present details of our future plans including updates on our integration with Ensembl, and how we plan to improve our support for the microbial research community. Software and data are made available without restriction via our website, online tools platform and programmatic interfaces (available under an Apache 2.0 license). Data updates are synchronised with Ensembl's release cycle.

Immunohistochemical Evaluation of Adaptor Protein FAM159B Expression in Normal and Neoplastic Human Tissues

FAM159B is a so-called adaptor protein. These proteins are essential components in numerous cell signalling pathways. However, little is known regarding FAM159B expression in normal and neoplastic human tissues. The commercially available rabbit polyclonal anti-human FAM159B antibody HPA011778 was initially characterised for its specificity using Western blot analyses and immunocytochemistry and then applied to a large series of formalin-fixed, paraffin-embedded normal and neoplastic human tissue samples. Confirmation of FAM159B's predicted size and antibody specificity was achieved in BON-1 cells, a neuroendocrine tumour cell line endogenously expressing FAM159B, using targeted siRNA. Immunocytochemical experiments additionally revealed cytoplasmic expression of the adaptor protein. Immunohistochemical staining detected FAM159B expression in neuronal and neuroendocrine tissues such as the cortex, the trigeminal ganglia, dorsal root and intestinal ganglia, the pancreatic islets and the neuroendocrine cells of the bronchopulmonary and gastrointestinal tract, but also in the syncytiotrophoblasts of the placenta. FAM159B was also expressed in many of the 28 tumour entities investigated, with high levels in medullary and anaplastic thyroid carcinomas, parathyroid adenomas, lung and ovarian carcinomas, lymphomas and neuroendocrine tumours of different origins. The antibody HPA011778 can act as a useful tool for basic research and identifying FAM159B expression in tissue samples.

The high-throughput production of membrane proteins

Membrane proteins, found at the junctions between the outside world and the inner workings of the cell, play important roles in human disease and are used as biosensors. More than half of all therapeutics directly affect membrane protein function while nanopores enable DNA sequencing. The structural and functional characterisation of membrane proteins is therefore crucial. However, low levels of naturally abundant protein and the hydrophobic nature of membrane proteins makes production difficult. To maximise success, high-throughput strategies were developed that rely upon simple screens to identify successful constructs and rapidly exclude those unlikely to work. Parameters that affect production such as expression host, membrane protein origin, expression vector, fusion-tags, encapsulation reagent and solvent composition are screened in parallel. In this way, constructs with divergent requirements can be produced for a variety of structural applications. As structural techniques advance, sample requirements will change. Single-particle cryo-electron microscopy requires less protein than crystallography and as cryo-electron tomography and time-resolved serial crystallography are developed new sample production requirements will evolve. Here we discuss different methods used for the high-throughput production of membrane proteins for structural biology.

Novel Autoantibodies in Idiopathic Small Fiber Neuropathy

Objective

Small fiber neuropathy (SFN) is clinically and etiologically heterogeneous. Although autoimmunity has been postulated to be pathophysiologically important in SFN, few autoantibodies have been described. We aimed to identify autoantibodies associated with idiopathic SFN (iSFN) by a novel high‐throughput protein microarray platform that captures autoantibodies expressed in the native conformational state.

Methods

Sera from 58 SFN patients and 20 age‐ and gender‐matched healthy controls (HCs) were screened against >1,600 immune‐related antigens. Fluorescent unit readout and postassay imaging were performed, followed by composite data normalization and protein fold change (pFC) analysis. Analysis of an independent validation cohort of 33 SFN patients against the same 20 HCs was conducted to identify reproducible proteins in both cohorts.

Results

Nine autoantibodies were screened with statistical significance and pFC criteria in both cohorts, with at least 50% change in serum levels. Three proteins showed consistently high fold changes in main and validation cohorts: MX1 (FC = 2.99 and 3.07, respectively, p = 0.003, q = 0.076), DBNL (FC = 2.11 and 2.16, respectively, p = 0.009, q < 0.003), and KRT8 (FC = 1.65 and 1.70, respectively, p = 0.043, q < 0.003). Further subgroup analysis into iSFN and SFN by secondary causes (secondary SFN) in the main cohort showed that MX1 is higher in iSFN compared to secondary SFN (FC = 1.61 vs 0.106, p = 0.009).

Interpretation

Novel autoantibodies MX1, DBNL, and KRT8 are found in iSFN. MX1 may allow diagnostic subtyping of iSFN patients. ANN NEUROL 2022;91:66–77

Telomerase-positive circulating tumor cells are associated with poor prognosis via a neutrophil-mediated inflammatory immune environment in glioma

BACKGROUND

Gliomas are the most common aggressive cancer in the central nervous system. Considering the difficulty in monitoring glioma response and progression, an approach is needed to evaluate the progression or survival of patients with glioma. We propose an application to facilitate clinical detection and treatment monitoring in glioma patients by using telomerase-positive circulating tumor cells (CTCs) and to further evaluate the relationship between the immune microenvironment and CTCs in glioma patients.

METHODS

From October 2014 to June 2017, 106 patients newly diagnosed with glioma were enrolled. We used the telomerase reverse transcriptase CTC detection method to detect and analyze the CTC statuses of glioma patients before and after surgery. FlowSight and FISH confirmed the CTCs detected by the telomerase-based method. To verify the correlation between CTCs and the immune response, peripheral white blood cell RNA sequencing was performed.

RESULTS

CTCs were common in the peripheral blood of glioma patients and were not correlated with the pathological classification or grade of patients. The results showed that the presence of postoperative CTCs but not preoperative CTCs in glioma patients was a poor prognostic factor. The level of postoperative CTCs, which predicts a poor prognosis after surgery, may be associated with neutrophils. RNA sequencing suggested that postoperative CTCs were positively correlated with innate immune responses, especially the activation of neutrophils and the generation of neutrophil extracellular traps, but negatively correlated with the cytotoxic response.

CONCLUSIONS

Our results showed that telomerase-positive CTCs can predict a poor prognosis of patients with glioma. Our results also showed a correlation between CTCs and the immune macroenvironment, which provides a new perspective for the treatment of glioma.

dbPTM in 2022: an updated database for exploring regulatory networks and functional associations of protein post-translational modifications

Protein post-translational modifications (PTMs) play an important role in different cellular processes. In view of the importance of PTMs in cellular functions and the massive data accumulated by the rapid development of mass spectrometry (MS)-based proteomics, this paper presents an update of dbPTM with over 2 777 000 PTM substrate sites obtained from existing databases and manual curation of literature, of which more than 2 235 000 entries are experimentally verified. This update has manually curated over 42 new modification types that were not included in the previous version. Due to the increasing number of studies on the mechanism of PTMs in the past few years, a great deal of upstream regulatory proteins of PTM substrate sites have been revealed. The updated dbPTM thus collates regulatory information from databases and literature, and merges them into a protein-protein interaction network. To enhance the understanding of the association between PTMs and molecular functions/cellular processes, the functional annotations of PTMs are curated and integrated into the database. In addition, the existing PTM-related resources, including annotation databases and prediction tools are also renewed. Overall, in this update, we would like to provide users with the most abundant data and comprehensive annotations on PTMs of proteins. The updated dbPTM is now freely accessible at https://awi.cuhk.edu.cn/dbPTM/.

The IntAct database: efficient access to fine-grained molecular interaction data

The IntAct molecular interaction database (https://www.ebi.ac.uk/intact) is a curated resource of molecular interactions, derived from the scientific literature and from direct data depositions. As of August 2021, IntAct provides more than one million binary interactions, curated by twelve global partners of the International Molecular Exchange consortium, for which the IntAct database provides a shared curation and dissemination platform. The IMEx curation policy has always emphasised a fine-grained data and curation model, aiming to capture the relevant experimental detail essential for the interpretation of the provided molecular interaction data. Here, we present recent curation focus and progress, as well as a completely redeveloped website which presents IntAct data in a much more user-friendly and detailed way.

Quadruple and Truncated MEK3 Mutants Identified from Acute Lymphoblastic Leukemia Promote Degradation and Enhance Proliferation

Compared to other ethnicities, Hispanic children incur the highest rates of leukemia, and most cases are diagnosed as Acute Lymphoblastic Leukemia (ALL). Despite improved treatment and survival for ALL, disproportionate health outcomes in Hispanics persist. Thus, it is essential to identify oncogenic mutations within this demographic to aid in the development of new strategies to diagnose and treat ALL. Using whole-exome sequencing, five single nucleotide polymorphisms within mitogen-activated protein kinase 3 (MAP2K3) were identified in an ALL cancer patient library from the U.S./Mexico border. MAP2K3 R26T and P11T are located near the substrate-binding site, while R65L and R67W localized to the kinase domain. Truncated-MAP2K3 mutant Q73* was also identified. Transfection in HEK293 cells showed that the quadruple-MEK3 mutant (4M-MEK3) impacted protein stability, inducing degradation and reducing expression. The expression of 4M-MEK3 could be rescued by cysteine/serine protease inhibition, and proteasomal degradation of truncated-MEK3 occurred in a ubiquitin-independent manner. MEK3 mutants displayed reduced auto-phosphorylation and enzymatic activity, as seen by decreases in p38 phosphorylation. Furthermore, uncoupling of the MEK3/p38 signaling pathway resulted in less suppressive activity on HEK293 cell viability. Thus, disruption of MEK3 activation may promote proliferative signals in ALL. These findings suggest that MEK3 represents a potential therapeutic target for treating ALL.

A resource for improved predictions of Trypanosoma and Leishmania protein three-dimensional structure

AlphaFold2 and RoseTTAfold represent a transformative advance for predicting protein structure. They are able to make very high-quality predictions given a high-quality alignment of the protein sequence with related proteins. These predictions are now readily available via the AlphaFold database of predicted structures and AlphaFold or RoseTTAfold Colaboratory notebooks for custom predictions. However, predictions for some species tend to be lower confidence than model organisms. Problematic species include Trypanosoma cruzi and Leishmania infantum: important unicellular eukaryotic human parasites in an early-branching eukaryotic lineage. The cause appears to be due to poor sampling of this branch of life (Discoba) in the protein sequences databases used for the AlphaFold database and ColabFold. Here, by comprehensively gathering openly available protein sequence data for Discoba species, significant improvements to AlphaFold2 protein structure prediction over the AlphaFold database and ColabFold are demonstrated. This is made available as an easy-to-use tool for the parasitology community in the form of Colaboratory notebooks for generating multiple sequence alignments and AlphaFold2 predictions of protein structure for Trypanosoma, Leishmania and related species.

De novo identification of mammalian ciliary motility proteins using cryo-EM

Dynein-decorated doublet microtubules (DMTs) are critical components of the oscillatory molecular machine of cilia, the axoneme, and have luminal surfaces patterned periodically by microtubule inner proteins (MIPs). Here we present an atomic model of the 48-nm repeat of a mammalian DMT, derived from a cryoelectron microscopy (cryo-EM) map of the complex isolated from bovine respiratory cilia. The structure uncovers principles of doublet microtubule organization and features specific to vertebrate cilia, including previously unknown MIPs, a luminal bundle of tektin filaments, and a pentameric dynein-docking complex. We identify a mechanism for bridging 48- to 24-nm periodicity across the microtubule wall and show that loss of the proteins involved causes defective ciliary motility and laterality abnormalities in zebrafish and mice. Our structure identifies candidate genes for diagnosis of ciliopathies and provides a framework to understand their functions in driving ciliary motility.

IID 2021: towards context-specific protein interaction analyses by increased coverage, enhanced annotation and enrichment analysis

Improved bioassays have significantly increased the rate of identifying new protein-protein interactions (PPIs), and the number of detected human PPIs has greatly exceeded early estimates of human interactome size. These new PPIs provide a more complete view of disease mechanisms but precise understanding of how PPIs affect phenotype remains a challenge. It requires knowledge of PPI context (e.g. tissues, subcellular localizations), and functional roles, especially within pathways and protein complexes. The previous IID release focused on PPI context, providing networks with comprehensive tissue, disease, cellular localization, and druggability annotations. The current update adds developmental stages to the available contexts, and provides a way of assigning context to PPIs that could not be previously annotated due to insufficient data or incompatibility with available context categories (e.g. interactions between membrane and cytoplasmic proteins). This update also annotates PPIs with conservation across species, directionality in pathways, membership in large complexes, interaction stability (i.e. stable or transient), and mutation effects. Enrichment analysis is now available for all annotations, and includes multiple options; for example, context annotations can be analyzed with respect to PPIs or network proteins. In addition to tabular view or download, IID provides online network visualization. This update is available at http://ophid.utoronto.ca/iid.

Wnt signaling mediates acquisition of blood-brain barrier properties in naïve endothelium derived from human pluripotent stem cells

Endothelial cells (ECs) in the central nervous system (CNS) acquire their specialized blood-brain barrier (BBB) properties in response to extrinsic signals, with Wnt/β-catenin signaling coordinating multiple aspects of this process. Our knowledge of CNS EC development has been advanced largely by animal models, and human pluripotent stem cells (hPSCs) offer the opportunity to examine BBB development in an in vitro human system. Here we show that activation of Wnt signaling in hPSC-derived naïve endothelial progenitors, but not in matured ECs, leads to robust acquisition of canonical BBB phenotypes including expression of GLUT-1, increased claudin-5, decreased PLVAP and decreased permeability. RNA-seq revealed a transcriptome profile resembling ECs with CNS-like characteristics, including Wnt-upregulated expression of LEF1, APCDD1, and ZIC3. Together, our work defines effects of Wnt activation in naïve ECs and establishes an improved hPSC-based model for interrogation of CNS barriergenesis.

Modeling sequence-space exploration and emergence of epistatic signals in protein evolution

During their evolution, proteins explore sequence space via an interplay between random mutations and phenotypic selection. Here, we build upon recent progress in reconstructing data-driven fitness landscapes for families of homologous proteins, to propose stochastic models of experimental protein evolution. These models predict quantitatively important features of experimentally evolved sequence libraries, like fitness distributions and position-specific mutational spectra. They also allow us to efficiently simulate sequence libraries for a vast array of combinations of experimental parameters like sequence divergence, selection strength, and library size. We showcase the potential of the approach in reanalyzing two recent experiments to determine protein structure from signals of epistasis emerging in experimental sequence libraries. To be detectable, these signals require sufficiently large and sufficiently diverged libraries. Our modeling framework offers a quantitative explanation for different outcomes of recently published experiments. Furthermore, we can forecast the outcome of time- and resource-intensive evolution experiments, opening thereby a way to computationally optimize experimental protocols.

R. Buzalaf M, Crespo-Lopez ME, Lima RR. Maternal methylmercury exposure changes the proteomic profile of the offspring's salivary glands: Prospects on translational toxicology

Background

Methylmercury (MeHg) remains a public health issue since developing organisms are particularly vulnerable to this environmental contaminant. This study investigated the effect of maternal MeHg exposure on the modulation of proteomic profile of parotid (PA), submandibular (SM), and sublingual (SL) glands of offspring rats.

Materials and methods

Pregnant Wistar rats were daily exposed to 40 μg/kg MeHg during both gestational and lactation periods. The proteomic profiles of the major salivary glands of the offspring rats were analyzed through mass spectrometry.

Results

The offspring rats exposed to MeHg showed significant alterations in the proteomic profiles of the PA, SM, and SL glands. Altered proteins were associated with cytoskeleton components, tissue morphogenesis, and response to stimulus and stress.

Conclusion

This original study showed that maternal MeHg exposure significantly modulates the expression of proteins and induces alterations in the proteomic profiles of developing salivary glands.

Hymenoptera Genome Database: new genomes and annotation datasets for improved go enrichment and orthologue analyses

We report an update of the Hymenoptera Genome Database (HGD; http://HymenopteraGenome.org), a genomic database of hymenopteran insect species. The number of species represented in HGD has nearly tripled, with fifty-eight hymenopteran species, including twenty bees, twenty-three ants, eleven wasps and four sawflies. With a reorganized website, HGD continues to provide the HymenopteraMine genomic data mining warehouse and JBrowse/Apollo genome browsers integrated with BLAST. We have computed Gene Ontology (GO) annotations for all species, greatly enhancing the GO annotation data gathered from UniProt with more than a ten-fold increase in the number of GO-annotated genes. We have also generated orthology datasets that encompass all HGD species and provide orthologue clusters for fourteen taxonomic groups. The new GO annotation and orthology data are available for searching in HymenopteraMine, and as bulk file downloads.

VARIDT 2.0: structural variability of drug transporter

The structural variability data of drug transporter (DT) are key for research on precision medicine and rational drug use. However, these valuable data are not sufficiently covered by the available databases. In this study, a major update of VARIDT (a database previously constructed to provide DTs' variability data) was thus described. First, the experimentally resolved structures of all DTs reported in the original VARIDT were discovered from PubMed and Protein Data Bank. Second, the structural variability data of each DT were collected by literature review, which included: (a) mutation-induced spatial variations in folded state, (b) difference among DT structures of human and model organisms, (c) outward/inward-facing DT conformations and (d) xenobiotics-driven alterations in the 3D complexes. Third, for those DTs without experimentally resolved structural variabilities, homology modeling was further applied as well-established protocol to enrich such valuable data. As a result, 145 mutation-induced spatial variations of 42 DTs, 1622 inter-species structures originating from 292 DTs, 118 outward/inward-facing conformations belonging to 59 DTs, and 822 xenobiotics-regulated structures in complex with 57 DTs were updated to VARIDT (https://idrblab.org/varidt/ and http://varidt.idrblab.net/). All in all, the newly collected structural variabilities will be indispensable for explaining drug sensitivity/selectivity, bridging preclinical research with clinical trial, revealing the mechanism underlying drug-drug interaction, and so on.

A Missense Variant in the Bardet-Biedl Syndrome 2 Gene (BBS2) Leads to a Novel Syndromic Retinal Degeneration in the Shetland Sheepdog

Canine progressive retinal atrophy (PRA) describes a group of hereditary diseases characterized by photoreceptor cell death in the retina, leading to visual impairment. Despite the identification of multiple PRA-causing variants, extensive heterogeneity of PRA is observed across and within dog breeds, with many still genetically unsolved. This study sought to elucidate the causal variant for a distinct form of PRA in the Shetland sheepdog, using a whole-genome sequencing approach. Filtering variants from a single PRA-affected Shetland sheepdog genome compared to 176 genomes of other breeds identified a single nucleotide variant in exon 11 of the Bardet-Biedl syndrome-2 gene (BBS2) (c.1222G>C; p.Ala408Pro). Genotyping 1386 canids of 155 dog breeds, 15 cross breeds and 8 wolves indicated the c.1222G>C variant was only segregated within Shetland sheepdogs. Out of 505 Shetland sheepdogs, seven were homozygous for the variant. Clinical history and photographs for three homozygotes indicated the presence of a novel phenotype. In addition to PRA, additional clinical features in homozygous dogs support the discovery of a novel syndromic PRA in the breed. The development and utilization of a diagnostic DNA test aim to prevent the mutation from becoming more prevalent in the breed.

Reversible processes in collagen dehydration: A molecular dynamics study

Collagen dehydration is an unavoidable damaging process that causes the lack of fibers' physical properties and it is usually irreversible. However, the identification of low hydration conditions that permit a recovering of initial collagen features after a rehydration treatment is particularly of interest. Monitoring structural changes by means of MD simulations, we investigated the hydration-dehydration-rehydration cycle of two microfibril models built on different fragments of the sequence of rat tail collagen type I. The microfibrils have different hydropathic features, to investigate the influence of amino acid composition on the whole process. We showed that with low hydration at a level corresponding to the first shell, microfibril gains in compactness and tubularity. Crucially, some water molecules remain trapped inside the fibrils, allowing, by rehydrating, a recovery of the initial collagen structural features. Water rearranges in cluster around the protein, and its first layer is more anchored to the surface. However, these changes in distribution and mobility in low hydration conditions get back with rehydration.

Chromosome Level Genome Assembly and Annotation of Highly Invasive Japanese Stiltgrass (Microstegium vimineum)

The invasive Japanese stiltgrass (Microstegium vimineum) affects a wide range of ecosystems and threatens biodiversity across the eastern USA. However, the mechanisms underlying rapid adaptation, plasticity, and epigenetics in the invasive range are largely unknown. We present a chromosome-level assembly for M. vimineum to investigate genome dynamics, evolution, adaptation, and the genomics of phenotypic plasticity. We generated a 1.12-Gb genome with scaffold N50 length of 53.44 Mb respectively, taking a de novo assembly approach that combined PacBio and Dovetail Genomics Omni-C sequencing. The assembly contains 23 pseudochromosomes, representing 99.96% of the genome. BUSCO assessment indicated that 80.3% of Poales gene groups are present in the assembly. The genome is predicted to contain 39,604 protein-coding genes, of which 26,288 are functionally annotated. Furthermore, 66.68% of the genome is repetitive, of which unclassified (35.63%) and long-terminal repeat (LTR) retrotransposons (26.90%) are predominant. Similar to other grasses, Gypsy (41.07%) and Copia (32%) are the most abundant LTR-retrotransposon families. The majority of LTR-retrotransposons are derived from a significant expansion in the past 1-2 Myr, suggesting the presence of relatively young LTR-retrotransposon lineages. We find corroborating evidence from Ks plots for a stiltgrass-specific duplication event, distinct from the more ancient grass-specific duplication event. The assembly and annotation of M. vimineum will serve as an essential genomic resource facilitating studies of the invasion process, the history and consequences of polyploidy in grasses, and provides a crucial tool for natural resource managers.

Investigating the interaction between organic anion transporter 1 and ochratoxin A: An in silico structural study to depict early molecular events of substrate recruitment and the impact of single point mutations

Organic anion transporters (OATs) belong to a subgroup of the solute carrier 22 transporter family. OATs have a central role in xenobiotic disposition affecting the toxicokinetics of its substrates and inter-individual differences in their expression, activity and function impact both toxicokinetics and toxicodynamics. Amongst OATs, OAT1 (solute carrier family 22 member 6) is involved in the urinary excretion of many xenobiotics bringing substrates into renal proximal tubular cells which can then be secreted across the apical membrane into the tubule lumen. The mycotoxin ochratoxin A has been shown to have a high affinity for OAT1, which is an important renal transporter involved in its urinary excretion. Nowadays, molecular modeling techniques are widely applied to assess protein-ligand interactions and may provide a tool to depict the mechanic of xenobiotic action be it toxicokinetics or toxicodynamics. This work provides a structured pipeline consisting of docking and molecular dynamic simulations to study OAT1-ligand interactions and the impact of OAT1 polymorphisms on such interactions. Such a computational structure-based analytical framework allowed to: i) model OAT1-substrate complex formation and depict the features correlating its sequence, structure and its capability to recruit substrates; and ii) investigate the impact of OAT1 missense mutations on substrate recruitment. Perspectives on applying such a structured pipeline to xenobiotic-metabolising enzymes are discussed.

Transcriptomic Profile of the Cockle Cerastoderma edule Exposed to Seasonal Diarrhetic Shellfish Toxin Contamination

Bivalves constitute an important source of proteins for human consumption, but some accumulate biotoxins such as diarrhetic shellfish toxins (DSTs), constituting a risk to human health. The cockle Cerastoderma edule is one of the most important species harvested in the Portuguese coast but also one of the most affected species due to recurrent DSTs exposure. However, little is known regarding the effects of the toxins produced by blooming dinoflagellates on C. edule. Herein, we explore the Differentially Expressed Genes (DEGs) of two tissues (gills and digestive gland) from wild cockles sampled in Portugal, through their whole transcriptomic response in two different seasons (exposed and not exposed to DSTs). The de novo transcriptome assembly returned 684,723 contigs, N50 of 1049, and 98.53% completeness. Altogether, 1098 DEGs were identified, of which 353 DEGs were exclusive for the digestive gland, 536 unique for the gills and 209 DEGs were common. Among DEGs were identified known DSTs-biomarkers including glutathione peroxidase, glutathione S-transferase, superoxide dismutase, cytochrome P450, ABC transporters, actin and tubulin-related proteins, Heat shock proteins and complement C1Q-like proteins. This study provides the first transcriptomic profile of C. edule, giving new insights about its molecular responses under different environmental conditions of DSTs exposure.

System-Wide Quantitative N-Glycoproteomic Analysis from K562 Cells and Mouse Liver Tissues

As a key regulator of many biological processes, glycosylation is an essential post-translational modification (PTM) in the living system. Over 50% of human proteins are known to be glycosylated. Alterations in glycoproteins are directly linked to many diseases, making it crucial to understand system-wide glycosylation changes. The majority of known glycoproteins are from plasma membrane; however, glycosylation is a dynamic process that occurs throughout multiple subcellular organelles and involves sets of enzymes, chaperones, transporters, and sugar donor molecules. Many glycoproteins are expressed not only in plasma membranes but also in subcellular organelles. Here, we developed a mass-spectrometry-based quantitative workflow for the system-wide N-glycoproteomic analysis of membrane and cytosolic proteins extracted using a MEM-PER kit. The kit facilitates the extraction and solubilization of both membrane and cytosolic proteins in a simple, efficient, and reproducible manner. We analyzed the K562 cell line and mouse liver tissue to evaluate this approach. A total of 934 glycosites, 5154 glycopeptides, and 536 glycoproteins from the K562 cell line and a total of 1449 glycosites, 7549 glycopeptides, and 660 glycoproteins from mouse liver tissue were identified. This simple and reproducible approach provides a unique way to understand system-wide glycosylation in biological processes and enables the identification and quantitation of glycan profiles at glycosylation sites in proteins.

Establishment of Host-Algal Endosymbioses: Genetic Response to Symbiont Versus Prey in a Sponge Host

The freshwater sponge Ephydatia muelleri and its Chlorella-like algal partner is an emerging model for studying animal: algal endosymbiosis. The sponge host is a tractable laboratory organism, and the symbiotic algae are easily cultured. We took advantage of these traits to interrogate questions about mechanisms that govern the establishment of durable intracellular partnerships between hosts and symbionts in facultative symbioses. We modified a classical experimental approach to discern the phagocytotic mechanisms that might be co-opted to permit persistent infections, and identified genes differentially expressed in sponges early in the establishment of endosymbiosis. We exposed algal-free E. muelleri to live native algal symbionts and potential food items (bacteria and native heat-killed algae), and performed RNA-Seq to compare patterns of gene expression among treatments. We found a relatively small but interesting suite of genes that are differentially expressed in the host exposed to live algal symbionts, and a larger number of genes triggered by host exposure to heat-killed algae. The upregulated genes in sponges exposed to live algal symbionts were mostly involved in endocytosis, ion transport, metabolic processes, vesicle-mediated transport, and oxidation–reduction. One of the host genes, an ATP-Binding Cassette transporter that is downregulated in response to live algal symbionts, was further evaluated for its possible role in the establishment of the symbiosis. We discuss the gene expression profiles associated with host responses to living algal cells in the context of conditions necessary for long-term residency within host cells by phototrophic symbionts as well as the genetic responses to sponge phagocytosis and immune-driven pathways.

Nephronophthisis-Pathobiology and Molecular Pathogenesis of a Rare Kidney Genetic Disease

The exponential rise in our understanding of the aetiology and pathophysiology of genetic cystic kidney diseases can be attributed to the identification of cystogenic genes over the last three decades. The foundation of this was laid by positional cloning strategies which gradually shifted towards next-generation sequencing (NGS) based screenings. This shift has enabled the discovery of novel cystogenic genes at an accelerated pace unlike ever before and, most notably, the past decade has seen the largest increase in identification of the genes which cause nephronophthisis (NPHP). NPHP is a monogenic autosomal recessive cystic kidney disease caused by mutations in a diverse clade of over 26 identified genes and is the most common genetic cause of renal failure in children. NPHP gene types present with some common pathophysiological features alongside a diverse range of extra-renal phenotypes associated with specific syndromic presentations. This review provides a timely update on our knowledge of this disease, including epidemiology, pathophysiology, anatomical and molecular features. We delve into the diversity of the NPHP causing genes and discuss known molecular mechanisms and biochemical pathways that may have possible points of intersection with polycystic kidney disease (the most studied renal cystic pathology). We delineate the pathologies arising from extra-renal complications and co-morbidities and their impact on quality of life. Finally, we discuss the current diagnostic and therapeutic modalities available for disease management, outlining possible avenues of research to improve the prognosis for NPHP patients.

The Arabidopsis PeptideAtlas: Harnessing worldwide proteomics data to create a comprehensive community proteomics resource

We developed a resource, the Arabidopsis PeptideAtlas (www.peptideatlas.org/builds/arabidopsis/), to solve central questions about the Arabidopsis thaliana proteome, such as the significance of protein splice forms and post-translational modifications (PTMs), or simply to obtain reliable information about specific proteins. PeptideAtlas is based on published mass spectrometry (MS) data collected through ProteomeXchange and reanalyzed through a uniform processing and metadata annotation pipeline. All matched MS-derived peptide data are linked to spectral, technical, and biological metadata. Nearly 40 million out of ∼143 million MS/MS (tandem MS) spectra were matched to the reference genome Araport11, identifying ∼0.5 million unique peptides and 17,858 uniquely identified proteins (only isoform per gene) at the highest confidence level (false discovery rate 0.0004; 2 non-nested peptides ≥9 amino acid each), assigned canonical proteins, and 3,543 lower-confidence proteins. Physicochemical protein properties were evaluated for targeted identification of unobserved proteins. Additional proteins and isoforms currently not in Araport11 were identified that were generated from pseudogenes, alternative start, stops, and/or splice variants, and small Open Reading Frames; these features should be considered when updating the Arabidopsis genome. Phosphorylation can be inspected through a sophisticated PTM viewer. PeptideAtlas is integrated with community resources including TAIR, tracks in JBrowse, PPDB, and UniProtKB. Subsequent PeptideAtlas builds will incorporate millions more MS/MS data.

Structural Dynamics and Catalytic Mechanism of ATP13A2 (PARK9) from Simulations

ATP13A2 is a gene encoding a protein of the P5B subfamily of ATPases and is a PARK gene. Molecular defects of the gene are mainly associated with variations of Parkinson's disease (PD). Despite the established importance of the protein in regulating neuronal integrity, the three-dimensional structure of the protein currently remains unresolved crystallographically. We have modeled the structure and reactivity of the full-length protein in its E1-ATP state. Using molecular dynamics (MD), quantum cluster, and quantum mechanical/molecular mechanical (QM/MM) methods, we aimed at describing the main catalytic reaction, leading to the phosphorylation of Asp513. Our MD simulations suggest that two positively charged Mg²⁺ cations are present at the active site during the catalytic reaction, stabilizing a specific triphosphate binding mode. Using QM/MM calculations, we subsequently calculated the reaction profiles for the phosphoryl transfer step in the presence of one and two Mg²⁺ cations. The calculated barrier heights in both cases are found to be ∼12.5 and 7.5 kcal mol^-1, respectively. We elucidated details of the catalytically competent ATP conformation and the binding mode of the second Mg²⁺ cofactor. We also examined the role of the conserved Arg686 and Lys859 catalytic residues. We observed that by significantly lowering the barrier height of the ATP cleavage reaction, Arg686 had major effect on the reaction. The removal of Arg686 increased the barrier height for the ATP cleavage by more than 5.0 kcal mol^-1 while the removal of key electrostatic interactions created by Lys859 to the γ-phosphate and Asp513 destabilizes the reactant state. When missense mutations occur in close proximity to an active site residue, they can interfere with the barrier height of the reaction, which can halt the normal enzymatic rate of the protein. We also found large binding pockets in the full-length structure, including a transmembrane domain pocket, which is likely where the ATP13A2 cargo binds.

Epinephrine affects gene expression levels and has a complex effect on biofilm formation in M icrococcus luteus strain C01 isolated from human skin

In this study, the effect of epinephrine on the biofilm formation of Micrococcus luteus C01 isolated from human skin was investigated in depth for the first time. This hormone has a complex effect on biofilms in various systems. In a system with polytetrafluoroethylene (PTFE) cubes, treatment with epinephrine at a physiological concentration of 4.9 × 10^-9 M increased the total amount of 72-h biofilm biomass stained with crystal violet and increased the metabolic activity of biofilms, but at higher and lower concentrations, the treatment had no significant effect. On glass fiber filters, treatment with the hormone decreased the number of colony forming units (CFUs) and changed the aggregation but did not affect the metabolic activity of biofilm cells. In glass bottom plates examined by confocal microscopy, epinephrine notably inhibited the growth of biofilms. RNA-seq analysis and RT-PCR demonstrated reproducible upregulation of genes encoding Fe-S cluster assembly factors and cyanide detoxification sulfurtransferase, whereas genes encoding the co-chaperone GroES, the LysE superfamily of lysine exporters, short-chain alcohol dehydrogenase and the potential c-di-GMP phosphotransferase were downregulated. Our results suggest that epinephrine may stimulate matrix synthesis in M. luteus biofilms, thereby increasing the activity of NAD(H) oxidoreductases. Potential c-di-GMP pathway proteins are essential in these processes.

The Function of the Kynurenine Pathway in the Placenta: A Novel Pharmacotherapeutic Target?

(L-)tryptophan is metabolized via the kynurenine pathway into several kynurenine metabolites with distinct functions. Dysfunction of the kynurenine pathway can lead to impairments in vascular regulation, immune regulation, and tolerance. The first and rate limiting enzyme of this pathway, indoleamine 2,3-dioxygenase (IDO), is highly expressed in the placenta and reduced in placentas from complicated pregnancies. IDO is essential during pregnancy, as IDO inhibition in pregnant mice resulted in fetal loss. However, the exact function of placental IDO, as well as its exact placental localization, remain controversial. This review identified that two isoforms of IDO; IDO1 and IDO2, are differently expressed between placental cells, suggesting spatial segregation. Furthermore, this review summarizes how the placental kynurenine pathway is altered in pregnancy complications, including recurrent miscarriage, preterm birth, preeclampsia, and fetal growth restriction. Importantly, we describe that these alterations do not affect maternally circulating metabolite concentrations, suggesting that the kynurenine pathway functions as a local signaling pathway. In the placenta, it is an important source of de novo placental NAD⁺ synthesis and regulates fetal tryptophan and kynurenine metabolite supply. Therefore, kynurenine pathway interventions might provide opportunities to treat pregnancy complications, and this review discusses how such treatment could affect placental function and pregnancy development.

Changes in the distribution of membrane lipids during growth of Thermotoga maritima at different temperatures: Indications for the potential mechanism of biosynthesis of ether-bound diabolic acid (membrane-spanning) lipids

Membrane-spanning lipids are present in a wide variety of archaea but they are rarely in bacteria. Nevertheless, the (hyper)thermophilic members of the order Thermotogales harbor tetraester, tetraether, and mixed ether/ester membrane-spanning lipids mostly composed of core lipids derived from diabolic acids, C_30, C₃₂ and C₃₄ dicarboxylic acids with two adjacent mid-chain methyl substituents. Lipid analysis of Thermotoga maritima across growth phases revealed a decrease of the relative abundance of fatty acids together with an increase of diabolic acids with independence of growth temperature. We also identified isomers of C₃₀ and C₃₂ diabolic acids, i.e. dicarboxylic acids with only one methyl group at C-15. Their distribution suggests they are products of the condensation reaction but preferably produced when the length of the acyl chains is not optimal. In comparison with growth at the optimal temperature of 80°C, an increase of glycerol ether-derived lipids was observed at 55°C. Besides, our analysis only detected diabolic acid-containing intact polar lipids with phosphoglycerol (PG) headgroups. Considering these findings, we hypothesize a biosynthetic pathway for the synthesis of membrane-spanning lipids based on PG polar lipid formation, suggesting that the protein catalyzing this process could be a membrane protein. We also identified, by genomic and protein domain analyses, a gene coding for a putative plasmalogen synthase homologue in T. maritima, which is also present in other bacteria producing sn1-alkyl ether lipids but not plasmalogens, suggesting it could be involved in the conversion of the ester to ether bond in the diabolic acids bound in membrane-spanning lipids. Importance Membrane-spanning lipids are unique compounds found in most archaeal membranes, but they are also present in specific bacterial groups like the Thermotogales. The synthesis and physiological role of membrane-spanning lipids in bacteria represent an evolutionary and biochemical open question that points to the differentiation of the membrane lipids composition. Understanding the formation of membrane-spanning lipids is crucial to solving this question and identifying the enzymatic and biochemical mechanism performing this procedure. In the present work, we found changes at the core lipid level, and we propose that the growth phase drives the biosynthesis of these lipids rather than temperature. Our results identified physiological conditions influencing the membrane-spanning lipids biosynthetic process which can further clarify the pathway leading to the biosynthesis of these compounds.

Temporal dynamics of base excision/single-strand break repair protein complex assembly/disassembly are modulated by the PARP/NAD+/SIRT6 axis

Assembly and disassembly of DNA repair protein complexes at DNA damage sites are essential for maintaining genomic integrity. Investigating factors coordinating assembly of the base excision repair (BER) proteins DNA polymerase β (Polβ) and XRCC1 to DNA lesion sites identifies a role for Polβ in regulating XRCC1 disassembly from DNA repair complexes and, conversely, demonstrates Polβ’s dependence on XRCC1 for complex assembly. LivePAR, a genetically encoded probe for live-cell imaging of poly(ADP-ribose) (PAR), reveals that Polβ and XRCC1 require PAR for repair-complex assembly, with PARP1 and PARP2 playing unique roles in complex dynamics. Further, BER complex assembly is modulated by attenuation/augmentation of NAD⁺ biosynthesis. Finally, SIRT6 does not modulate PARP1 or PARP2 activation but does regulate XRCC1 recruitment, leading to diminished Polβ abundance at sites of DNA damage. These findings highlight coordinated yet independent roles for PARP1, PARP2, and SIRT6 and their regulation by NAD⁺ bioavailability to facilitate BER.

Koczor et al. use quantitative confocal microscopy to characterize DNA-damage-induced poly(ADP-ribose) (PAR) formation and assembly/disassembly kinetics in human cells. These studies highlight the coordinated yet independent roles for XRCC1, POLΒ, PARP1, PARP2, and SIRT6 (and regulation by NAD⁺) to facilitate BER/SSBR protein complex dynamics.

Switching an active site helix in dihydrofolate reductase reveals limits to subdomain modularity

To what degree are individual structural elements within proteins modular such that similar structures from unrelated proteins can be interchanged? We study subdomain modularity by creating 20 chimeras of an enzyme, Escherichia coli dihydrofolate reductase (DHFR), in which a catalytically important, 10-residue α-helical sequence is replaced by α-helical sequences from a diverse set of proteins. The chimeras stably fold but have a range of diminished thermal stabilities and catalytic activities. Evolutionary coupling analysis indicates that the residues of this α-helix are under selection pressure to maintain catalytic activity in DHFR. Reversion to phenylalanine at key position 31 was found to partially restore catalytic activity, which could be explained by evolutionary coupling values. We performed molecular dynamics simulations using replica exchange with solute tempering. Chimeras with low catalytic activity exhibit nonhelical conformations that block the binding site and disrupt the positioning of the catalytically essential residue D27. Simulation observables and in vitro measurements of thermal stability and substrate-binding affinity are strongly correlated. Several E. coli strains with chromosomally integrated chimeric DHFRs can grow, with growth rates that follow predictions from a kinetic flux model that depends on the intracellular abundance and catalytic activity of DHFR. Our findings show that although α-helices are not universally substitutable, the molecular and fitness effects of modular segments can be predicted by the biophysical compatibility of the replacement segment.

Wilavidin* - a novel member of the avidin family that forms unique biotin-binding hexamers

Nature's optimization of protein functions is a highly intricate evolutionary process. In addition to optimal tertiary folding, the intramolecular recognition among the monomers that generate higher-order quaternary arrangements is driven by stabilizing interactions that have a pivotal role for ideal activity. Homotetrameric avidin and streptavidin are regularly utilized in many applications, whereby their ultra-high affinity toward biotin is dependent on their quaternary arrangements. In recent years, a new subfamily of avidins was discovered that comprises homodimers rather than tetramers, in which the high affinity toward biotin is maintained. Intriguingly, several of the respective dimers have been shown to assemble into higher-order cylindrical hexamers or octamers that dissociate into dimers upon biotin binding. Here, we present wilavidin, a newly discovered member of the dimeric subfamily, forming hexamers in the apo form, which are uniquely maintained upon biotin binding with six high-affinity binding sites. Removal of the short C-terminal segment of wilavidin resulted in the presence of the dimer only, thus emphasizing the role of this segment in stabilizing the hexamer. Utilization of a hexavalent biotin-binding form of avidin would be beneficial for expanding the biotechnological toolbox. Additionally, this unique family of dimeric avidins and their propensity to oligomerize to hexamers or octamers can serve as a basis for protein oligomerization and intermonomeric recognition as well as cumulative interactions that determine molecular assemblies.

Plant Metabolic Network 15: A resource of genome-wide metabolism databases for 126 plants and algae

To understand and engineer plant metabolism, we need a comprehensive and accurate annotation of all metabolic information across plant species. As a step towards this goal, we generated genome-scale metabolic pathway databases of 126 algal and plant genomes, ranging from model organisms to crops to medicinal plants (https://plantcyc.org). Of these, 104 have not been reported before. We systematically evaluated the quality of the databases, which revealed that our semi-automated validation pipeline dramatically improves the quality. We then compared the metabolic content across the 126 organisms using multiple correspondence analysis and found that Brassicaceae, Poaceae, and Chlorophyta appeared as metabolically distinct groups. To demonstrate the utility of this resource, we used recently published sorghum transcriptomics data to discover previously unreported trends of metabolism underlying drought tolerance. We also used single-cell transcriptomics data from the Arabidopsis root to infer cell type-specific metabolic pathways. This work shows the quality and quantity of our resource and demonstrates its wide-ranging utility in integrating metabolism with other areas of plant biology.

O-GlcNAcylation drives calcium signaling toward osteoblast differentiation: A bioinformatics-oriented study

This study aimed to reveal the possible mechanisms by which O-linked-N-acetylglucosaminylation (O-GlcNAcylation) regulates osteoblast differentiation using a series of bioinformatics-oriented experiments. To examine the influence of O-GlcNAcylation levels on osteoblast differentiation, osteoblastic MC3T3-E1 cells were treated with O-GlcNAc transferase (OGT) and O-GlcNAcase (OGA) inhibitors. Correlations between the levels of O-GlcNAcylation and the expression of osteogenic markers as well as OGT were evaluated by qPCR and western blotting. The O-GlcNAcylated proteins assumed to correlate with Runx2 expression were retrieved from several public databases and used for further bioinformatics analysis. Following the findings of the bioinformatics analysis, intracellular calcium ([Ca²⁺ ]_i ) was monitored in the cells treated with OGT and OGA inhibitors using a confocal laser-scanning microscope (CLS). The interaction effect between O-GlcNAcylation and [Ca²⁺ ]_i on osteogenic marker expression was determined using stable OGT knockdown MC3T3-E1 cells. O-GlcNAcylation was positively associated with osteoblast differentiation. The time-course profile of global O-GlcNAcylated proteins showed a distinctive pattern with different molecular weights during osteoblast differentiation. The expression pattern of several O-GlcNAcylated proteins was significantly similar to that of Runx2 expression. Bioinformatic analysis of the retrieved Runx2-related-O-GlcNAcylated-proteins revealed the importance of [Ca²⁺ ]_i . CLS showed that alteration of O-GlcNAcylation rapidly changed [Ca²⁺ ]_i in MC3T3-E1 cells. O-GlcNAcylation and [Ca²⁺ ]_i showed an interaction effect on the expression of osteogenic markers. OGT knockdown disrupted the [Ca²⁺ ]_i -induced expression changes of osteogenic markers. O-GlcNAcylation interacts with [Ca²⁺ ]_i and elicits osteoblast differentiation by regulating the expression of osteogenic markers.

Bakta: rapid and standardized annotation of bacterial genomes via alignment-free sequence identification

Command-line annotation software tools have continuously gained popularity compared to centralized online services due to the worldwide increase of sequenced bacterial genomes. However, results of existing command-line software pipelines heavily depend on taxon-specific databases or sufficiently well annotated reference genomes. Here, we introduce Bakta, a new command-line software tool for the robust, taxon-independent, thorough and, nonetheless, fast annotation of bacterial genomes. Bakta conducts a comprehensive annotation workflow including the detection of small proteins taking into account replicon metadata. The annotation of coding sequences is accelerated via an alignment-free sequence identification approach that in addition facilitates the precise assignment of public database cross-references. Annotation results are exported in GFF3 and International Nucleotide Sequence Database Collaboration (INSDC)-compliant flat files, as well as comprehensive JSON files, facilitating automated downstream analysis. We compared Bakta to other rapid contemporary command-line annotation software tools in both targeted and taxonomically broad benchmarks including isolates and metagenomic-assembled genomes. We demonstrated that Bakta outperforms other tools in terms of functional annotations, the assignment of functional categories and database cross-references, whilst providing comparable wall-clock runtimes. Bakta is implemented in Python 3 and runs on MacOS and Linux systems. It is freely available under a GPLv3 license at https://github.com/oschwengers/bakta. An accompanying web version is available at https://bakta.computational.bio.

Beta turn propensity and a model polymer scaling exponent identify intrinsically disordered phase-separating proteins

The complex cellular milieu can spontaneously demix, or phase separate, in a process controlled in part by intrinsically disordered (ID) proteins. A protein's propensity to phase separate is thought to be driven by a preference for protein-protein over protein-solvent interactions. The hydrodynamic size of monomeric proteins, as quantified by the polymer scaling exponent (v), is driven by a similar balance. We hypothesized that mean v, as predicted by protein sequence, would be smaller for proteins with a strong propensity to phase separate. To test this hypothesis, we analyzed protein databases containing subsets of proteins that are folded, disordered, or disordered and known to spontaneously phase separate. We find that the phase-separating disordered proteins, on average, had lower calculated values of v compared with their non-phase-separating counterparts. Moreover, these proteins had a higher sequence-predicted propensity for β-turns. Using a simple, surface area-based model, we propose a physical mechanism for this difference: transient β-turn structures reduce the desolvation penalty of forming a protein-rich phase and increase exposure of atoms involved in π/sp2 valence electron interactions. By this mechanism, β-turns could act as energetically favored nucleation points, which may explain the increased propensity for turns in ID regions (IDRs) utilized biologically for phase separation. Phase-separating IDRs, non-phase-separating IDRs, and folded regions could be distinguished by combining v and β-turn propensity. Finally, we propose a new algorithm, ParSe (partition sequence), for predicting phase-separating protein regions, and which is able to accurately identify folded, disordered, and phase-separating protein regions based on the primary sequence.

Integration of transcriptomic and proteomic approaches for snake venom profiling

INTRODUCTION

Snake venoms contain many protein and peptide isoforms with high levels of sequence variation, even within a single species.

AREAS COVERED

In this review, we highlight several examples, from both published and unpublished work in our lab, demonstrating how a combined venom gland transcriptome and proteome methodology allows for comprehensive characterization of venoms, including those from understudied rear-fanged snake species, and we provide recommendations for using these approaches.

EXPERT OPINION

When characterizing venoms, peptide mass fingerprinting using databases built predominately from protein sequences originating from model organisms can be disadvantageous, especially when the intention is to document protein diversity. Therefore, the use of species-specific venom gland transcriptomes corrects for the absence of these unique peptide sequences in databases. The integration of transcriptomics and proteomics improves the accuracy of either approach alone for venom profiling.

A comparative study of pan-genome methods for microbial organisms: Acinetobacter baumannii pan-genome reveals structural variation in antimicrobial resistance-carrying plasmids

Microbial organisms have diverse populations, where using a single linear reference sequence in comparative studies introduces reference-bias in downstream analyses, and leads to a failure to account for variability in the population. Recently, pan-genome graphs have emerged as an alternative to the traditional linear reference with many successful applications and a rapid increase in the number of methods available in the literature. Despite this enthusiasm, there has been no attempt at exploring these graph construction methods in depth, demonstrating their practical use. In this study, we aim to develop a general guide to help researchers who may want to incorporate pan-genomes in their analyses of microbial organisms. We evaluated the state-of-the art pan-genome construction tools to model a collection of 70 Acinetobacter baumannii strains. Our results suggest that all tools produced pan-genome graphs conforming to our expectations based on previous literature, and that their approach to homologue detection is likely to be the most influential in determining the final size and complexity of the pan-genome. The graphs overlapped most in the core pan-genome content while the cloud genes varied significantly among tools. We propose an alternative approach for pan-genome construction by combining two of the tools, Panaroo and Ptolemy, to further exploit them in downstream analyses, and demonstrate the effectiveness of our pipeline for structural variant calling in beta-lactam resistance genes in the same set of A. baumannii isolates, identifying various transposon structures for carbapenem resistance in chromosome, as well as plasmids. We identify a novel plasmid structure in two multidrug-resistant clinical isolates that had previously been studied, and which could be important for their resistance phenotypes.

Protein Structure, Models of Sequence Evolution, and Data Type Effects in Phylogenetic Analyses of Mitochondrial Data: A Case Study in Birds

Phylogenomic analyses have revolutionized the study of biodiversity, but they have revealed that estimated tree topologies can depend, at least in part, on the subset of the genome that is analyzed. For example, estimates of trees for avian orders differ if protein-coding or non-coding data are analyzed. The bird tree is a good study system because the historical signal for relationships among orders is very weak, which should permit subtle non-historical signals to be identified, while monophyly of orders is strongly corroborated, allowing identification of strong non-historical signals. Hydrophobic amino acids in mitochondrially-encoded proteins, which are expected to be found in transmembrane helices, have been hypothesized to be associated with non-historical signals. We tested this hypothesis by comparing the evolution of transmembrane helices and extramembrane segments of mitochondrial proteins from 420 bird species, sampled from most avian orders. We estimated amino acid exchangeabilities for both structural environments and assessed the performance of phylogenetic analysis using each data type. We compared those relative exchangeabilities with values calculated using a substitution matrix for transmembrane helices estimated using a variety of nuclear- and mitochondrially-encoded proteins, allowing us to compare the bird-specific mitochondrial models with a general model of transmembrane protein evolution. To complement our amino acid analyses, we examined the impact of protein structure on patterns of nucleotide evolution. Models of transmembrane and extramembrane sequence evolution for amino acids and nucleotides exhibited striking differences, but there was no evidence for strong topological data type effects. However, incorporating protein structure into analyses of mitochondrially-encoded proteins improved model fit. Thus, we believe that considering protein structure will improve analyses of mitogenomic data, both in birds and in other taxa.

Stem-loop binding protein and metal carcinogenesis

Pre-mRNA processing of the replication-dependent canonical histone mRNAs requires an endonucleolytic cleavage immediately after a conserved stem loop structure which occurs before RNA Pol II encounters any poly(A) signal. Thus, in contrast to all other eukaryotic mRNAs, the canonical histone mRNAs are not polyadenylated in their 3' ends. The binding of stem-loop binding protein (SLBP) to the stem loop structure of the histone mRNAs is required for this process. SLBP is also involved in regulation of histone mRNA nuclear export, degradation, and translation. Depletion of SLBP has been shown to induce polyadenylation of histone mRNAs and alteration of histone protein levels, which are considered to contribute to the observed aberrant cell cycle progress and genomic instability resulting from the loss of SLBP function. Recent studies have demonstrated that some heavy metal carcinogens, including arsenic and nickel, can induce the loss of SLBP and the gain of polyadenylation of canonical histone mRNAs. Polyadenylated canonical histone H3 can result in abnormal transcription, cell cycle arrest, genomic instability, and cell transformation, which links SLBP depletion and subsequent histone mRNA misprocessing to cancer. This review seeks to briefly summarize what is known about regulation of SLBP expression, consequences of SLBP depletion, its roles in cancer-related end points, with particular focus on metal-induced SLBP depletion and the potential of SLBP depletion as a new mechanism for metal-induced carcinogenesis.

Toward improvement of screening through mass spectrometry-based proteomics: ovarian cancer as a case study

Ovarian cancer is one of the leading causes of cancer related deaths affecting United States women. Early-stage detection of ovarian cancer has been linked to increased survival, however, current screening methods, such as biomarker testing, have proven to be ineffective in doing so. Therefore, further developments are necessary to be able to achieve positive patient prognosis. Ongoing efforts are being made in biomarker discovery towards clinical applications in screening for early-stage ovarian cancer. In this perspective, we discuss and provide examples for several workflows employing mass spectrometry-based proteomics towards protein biomarker discovery and characterization in the context of ovarian cancer; workflows include protein identification and characterization as well as intact protein profiling. We also discuss the opportunities to merge these workflows for a multiplexed approach for biomarkers. Lastly, we provide our insight as to future developments that may serve to enhance biomarker discovery workflows while also considering translational potential.

RNAInter v4.0: RNA interactome repository with redefined confidence scoring system and improved accessibility

Establishing an RNA-associated interaction repository facilitates the system-level understanding of RNA functions. However, as these interactions are distributed throughout various resources, an essential prerequisite for effectively applying these data requires that they are deposited together and annotated with confidence scores. Hence, we have updated the RNA-associated interaction database RNAInter (RNA Interactome Database) to version 4.0, which is freely accessible at http://www.rnainter.org or http://www.rna-society.org/rnainter/. Compared with previous versions, the current RNAInter not only contains an enlarged data set, but also an updated confidence scoring system. The merits of this 4.0 version can be summarized in the following points: (i) a redefined confidence scoring system as achieved by integrating the trust of experimental evidence, the trust of the scientific community and the types of tissues/cells, (ii) a redesigned fully functional database that enables for a more rapid retrieval and browsing of interactions via an upgraded user-friendly interface and (iii) an update of entries to >47 million by manually mining the literature and integrating six database resources with evidence from experimental and computational sources. Overall, RNAInter will provide a more comprehensive and readily accessible RNA interactome platform to investigate the regulatory landscape of cellular RNAs.

qPTMplants: an integrative database of quantitative post-translational modifications in plants

As a crucial molecular mechanism, post-translational modifications (PTMs) play critical roles in a wide range of biological processes in plants. Recent advances in mass spectrometry-based proteomic technologies have greatly accelerated the profiling and quantification of plant PTM events. Although several databases have been constructed to store plant PTM data, a resource including more plant species and more PTM types with quantitative dynamics still remains to be developed. In this paper, we present an integrative database of quantitative PTMs in plants named qPTMplants (http://qptmplants.omicsbio.info), which hosts 1 242 365 experimentally identified PTM events for 429 821 nonredundant sites on 123 551 proteins under 583 conditions for 23 PTM types in 43 plant species from 293 published studies, with 620 509 quantification events for 136 700 PTM sites on 55 361 proteins under 354 conditions. Moreover, the experimental details, such as conditions, samples, instruments and methods, were manually curated, while a variety of annotations, including the sequence and structural characteristics, were integrated into qPTMplants. Then, various search and browse functions were implemented to access the qPTMplants data in a user-friendly manner. Overall, we anticipate that the qPTMplants database will be a valuable resource for further research on PTMs in plants.

Determining oncogenic patterns and cancer predisposition through the transcriptomic profile in Mitchell-Riley syndrome with heterotopic gastric mucosa and duodenal atresia: a case report

Background

Homozygous mutations in the transcription factor RFX6 are the cause of the Mitchell–Riley syndrome (MRS) associating neonatal diabetes, congenital digestive system, such as biliary atresia, pancreatic hypoplasia, duodenal and/or jejunal atresia, intestinal malrotation, gallbladder aplasia, cholestasis. A constitutive inactivation of RFX6 leads also to gastric heterotopia. Application of RNA-seq in human diseases may help to better understand pathogenic mechanism of diseases and to predict the risk of developing chronic disorders and personalizing their prevention and treatment. We evaluated oncogenic patterns and cancer predisposition using the transcriptomic profile in a case of MRS with neonatal diabetes, duodenal atresia, and extensive intestinal tract gastric heterotopia.

Results

We signalled the interactors of RFX6 with other up and downregulated genes, that may be interested in severity of diabetic condition, in multi-organs impairment and cancer predisposition. Furthermore, several dysregulated genes are involved in biological processes that can lead to promote cancer including “Evading apoptosis” (BAD, BBC3, EGF, FGFR2, FLT3LG, HMOX1, HRAS, IFNAR2, IGF1R, IL12RB1, IL13RA1, IL15, IL2RB, IL2RG, IL6R, KEAP1, MGST1, PDGFA, PDGFRB, PIK3R3, RALB, RALGDS, RASSF1, SOS1, TGFA, TXNRD3), “Proliferation” (APC, BRAF, CCND2, CCND3, CCNE2, FGFR2, FLT3LG, FZD1, FZD6, HMOX1, HRAS, IGF1R, KEAP1, LRP6, MAPK3, MGST1, PDGFA, PDGFB, PDGFRB, RB1, SOS1, TGFA, TXNRD3, WNT10B), “Sustained angiogenesis” (BRAF, FGFR2, FLT3LG, HRAS, IGF1R, JAG1, MAPK3, NOTCH2, PDGFA, PDGFB, PDGFRB, SOS1, TGFA, TGFB1), “Genomic instability” (BAD, BBC3) and “Insensitivity to anti-growth signals” (SMAD2, TGFB1). We also inspected the signalings and their related genes in cancer, such as “PI3K signaling”, “ERK signaling”, “JAK-STAT signaling”, “Calcium signaling”, “Other RAS signaling”, “WNT signaling”.

Conclusions

In our MRS patient, we signaled the interactors of RFX6 with other up- and downregulated genes that may be related to severe diabetic condition, multi-organ impairment, and cancer predisposition. Notably, many dysregulated genes may lead to triggering carcinogenesis. The possibility of the patient developing cancer degeneration in heterotopic gastric mucosa and/or additional long-term tumoral sequelae is not excluded. Personalized prevention and treatment strategies should be proposed.

Therapeutic target database update 2022: facilitating drug discovery with enriched comparative data of targeted agents

Drug discovery relies on the knowledge of not only drugs and targets, but also the comparative agents and targets. These include poor binders and non-binders for developing discovery tools, prodrugs for improved therapeutics, co-targets of therapeutic targets for multi-target strategies and off-target investigations, and the collective structure-activity and drug-likeness landscapes of enhanced drug feature. However, such valuable data are inadequately covered by the available databases. In this study, a major update of the Therapeutic Target Database, previously featured in NAR, was therefore introduced. This update includes (a) 34 861 poor binders and 12 683 non-binders of 1308 targets; (b) 534 prodrug-drug pairs for 121 targets; (c) 1127 co-targets of 672 targets regulated by 642 approved and 624 clinical trial drugs; (d) the collective structure-activity landscapes of 427 262 active agents of 1565 targets; (e) the profiles of drug-like properties of 33 598 agents of 1102 targets. Moreover, a variety of additional data and function are provided, which include the cross-links to the target structure in PDB and AlphaFold, 159 and 1658 newly emerged targets and drugs, and the advanced search function for multi-entry target sequences or drug structures. The database is accessible without login requirement at: https://idrblab.org/ttd/.