LocTree3 prediction of localization

Phenanthrene metabolism in Panicum miliaceum: anatomical adaptations, degradation pathway, and computational analysis of a dioxygenase enzyme

Polycyclic aromatic compounds (PAHs) are persistent organic pollutants of environmental concern due to their potential impacts on food chain, with plants being particularly vulnerable. While plants can uptake, transport, and transform PAHs, the precise mechanisms underlying their localization and degradation are not fully understood. Here, a cultivation experiment conducted with Panicum miliaceum exposed different concentrations of phenanthrene (PHE). Intermediate PHE degradation compounds were identified via GC-MS analysis, leading to the proposal of a phytodegradation pathway featuring three significant benzene ring cleavage steps. Our results showed that P. miliaceum exhibited the ability to effectively degrade high levels of PHE, resulting in the production of various intermediate products through several chemical changes. Examination of the localization and anatomical characteristics revealed structural alterations linked to PHE stress, with an observed enhancement in PHE accumulation density in both roots and shoots as treatment levels increased. Following a 2-week aging period, a decrease in the amount of PHE accumulation was observed, along with a change in its localization. Bioinformatics analysis of the P. miliaceum 2-oxoglutarate-dependent dioxygenase (2-ODD) DAO-like protein revealed a 299 amino acid structure with two highly conserved domains, namely 2OG-FeII_Oxy and DIOX_N. Molecular docking analysis aligned with experimental results, strongly affirming the potential link and direct action of 2-ODD DAO-like protein with PHE. Our study highlights P. miliaceum capacity for PAHs degradation and elucidates the mechanisms behind enhanced degradation efficiency. By integrating experimental evidence with bioinformatics analysis, we offer valuable insights into the potential applications of plant-based remediation strategies for PAHs-contaminated environments.

SCLpred-ECL: Subcellular Localization Prediction by Deep N-to-1 Convolutional Neural Networks

The subcellular location of a protein provides valuable insights to bioinformaticians in terms of drug designs and discovery, genomics, and various other aspects of medical research. Experimental methods for protein subcellular localization determination are time-consuming and expensive, whereas computational methods, if accurate, would represent a much more efficient alternative. This article introduces an ab initio protein subcellular localization predictor based on an ensemble of Deep N-to-1 Convolutional Neural Networks. Our predictor is trained and tested on strict redundancy-reduced datasets and achieves 63% accuracy for the diverse number of classes. This predictor is a step towards bridging the gap between a protein sequence and the protein's function. It can potentially provide information about protein-protein interaction to facilitate drug design and processes like vaccine production that are essential to disease prevention.

Protein Engineering with Lightweight Graph Denoising Neural Networks

Protein engineering faces challenges in finding optimal mutants from a massive pool of candidate mutants. In this study, we introduce a deep-learning-based data-efficient fitness prediction tool to steer protein engineering. Our methodology establishes a lightweight graph neural network scheme for protein structures, which efficiently analyzes the microenvironment of amino acids in wild-type proteins and reconstructs the distribution of the amino acid sequences that are more likely to pass natural selection. This distribution serves as a general guidance for scoring proteins toward arbitrary properties on any order of mutations. Our proposed solution undergoes extensive wet-lab experimental validation spanning diverse physicochemical properties of various proteins, including fluorescence intensity, antigen-antibody affinity, thermostability, and DNA cleavage activity. More than 40% of ProtLGN-designed single-site mutants outperform their wild-type counterparts across all studied proteins and targeted properties. More importantly, our model can bypass the negative epistatic effect to combine single mutation sites and form deep mutants with up to seven mutation sites in a single round, whose physicochemical properties are significantly improved. This observation provides compelling evidence of the structure-based model's potential to guide deep mutations in protein engineering. Overall, our approach emerges as a versatile tool for protein engineering, benefiting both the computational and bioengineering communities.

Surface-Exposed Protein Moieties of Burkholderia cenocepacia J2315 in Microaerophilic and Aerobic Conditions

Burkholderia cepacia complex infections remain life-threatening to cystic fibrosis patients, and due to the limited eradication efficiency of current treatments, novel antimicrobial therapies are urgently needed. Surface proteins are among the best targets to develop new therapeutic strategies since they are exposed to the host's immune system. A surface-shaving approach was performed using Burkholderia cenocepacia J2315 to quantitatively compare the relative abundance of surface-exposed proteins (SEPs) expressed by the bacterium when grown under aerobic and microaerophilic conditions. After trypsin incubation of live bacteria and identification of resulting peptides by liquid chromatography coupled with mass spectrometry, a total of 461 proteins with ≥2 unique peptides were identified. Bioinformatics analyses revealed a total of 53 proteins predicted as localized at the outer membrane (OM) or extracellularly (E). Additionally, 37 proteins were predicted as moonlight proteins with OM or E secondary localization. B-cell linear epitope bioinformatics analysis of the proteins predicted to be OM and E-localized revealed 71 SEP moieties with predicted immunogenic epitopes. The protegenicity higher scores of proteins BCAM2761, BCAS0104, BCAL0151, and BCAL0849 point out these proteins as the best antigens for vaccine development. Additionally, 10 of the OM proteins also presented a high probability of playing important roles in adhesion to host cells, making them potential targets for passive immunotherapeutic approaches. The immunoreactivity of three of the OM proteins identified was experimentally demonstrated using serum samples from cystic fibrosis patients, validating our strategy for identifying immunoreactive moieties from surface-exposed proteins of potential interest for future immunotherapies development.

Chloroplast phosphate transporter CrPHT4-7 regulates phosphate homeostasis and photosynthesis in Chlamydomonas

In eukaryotic cells, phosphorus is assimilated and utilized primarily as phosphate (Pi). Pi homeostasis is mediated by transporters that have not yet been adequately characterized in green algae. This study reports on PHOSPHATE TRANSPORTER 4-7 (CrPHT4-7) from Chlamydomonas reinhardtii, a member of the PHT4 transporter family, which exhibits remarkable similarity to AtPHT4;4 from Arabidopsis (Arabidopsis thaliana), a chloroplastic ascorbate transporter. Using fluorescent protein tagging, we show that CrPHT4-7 resides in the chloroplast envelope membrane. Crpht4-7 mutants, generated by the CRISPR/Cas12a-mediated single-strand templated repair, show retarded growth, especially in high light, reduced ATP level, strong ascorbate accumulation, and diminished non-photochemical quenching in high light. On the other hand, total cellular phosphorous content was unaffected, and the phenotype of the Crpht4-7 mutants could not be alleviated by ample Pi supply. CrPHT4-7-overexpressing lines exhibit enhanced biomass accumulation under high light conditions in comparison with the wild-type strain. Expressing CrPHT4-7 in a yeast (Saccharomyces cerevisiae) strain lacking Pi transporters substantially recovered its slow growth phenotype, demonstrating that CrPHT4-7 transports Pi. Even though CrPHT4-7 shows a high degree of similarity to AtPHT4;4, it does not display any substantial ascorbate transport activity in yeast or intact algal cells. Thus, the results demonstrate that CrPHT4-7 functions as a chloroplastic Pi transporter essential for maintaining Pi homeostasis and photosynthesis in C. reinhardtii.

In-silico identification of potential inhibitors against FabI protein in Klebsiella pneumoniae

The development of new antimicrobial drugs is needed to combat multi-drug resistant and novel hypervirulent strains of Klebsiella pneumoniae (KPN) that are associated with increased morbidity and mortality globally. The FabI protein plays a crucial role in fatty acid biosynthesis and has been identified as an important target for in-silico, in-vitro, and in-vivo drug discovery. In this study we have used computer integrated-drug discovery approaches and binding-free energy calculations to identify three novel inhibitors (21272541, 67724550, and 67724551) of the FabI protein. All inhibitors showed strong affinity including van der Waals energy, electrostatic energy, polar and non-polar energies; however, the 21272541 compound was the most effective inhibitor and bound with the strongest affinity (ΔGbind -59.02 kcal/mol) to the FabI protein. Nevertheless, all three inhibitors are promising targets for new novel antimicrobial drugs that could contribute to the management of antimicrobial resistant KPN infections based on various computational analysis. Additional in-vitro and in-vivo clinical studies will be needed to confirm drug effectiveness for the treatment of KPN infections.

Unveiling the defensive role of Snakin-3, a member of the subfamily III of Snakin/GASA peptides in potatoes

KEY MESSAGE

The first in-depth characterization of a subfamily III Snakin/GASA member was performed providing experimental evidence on promoter activity and subcellular localization and unveiling a role of potato Snakin-3 in defense Snakin/GASA proteins share 12 cysteines in conserved positions in the C-terminal region. Most of them were involved in different aspects of plant growth and development, while a small number of these peptides were reported to have antimicrobial activity or participate in abiotic stress tolerance. In potato, 18 Snakin/GASA genes were identified and classified into three groups based on phylogenetic analysis. Snakin-1 and Snakin-2 are members of subfamilies I and II, respectively, and were reported to be implicated not only in defense against pathogens but also in plant development. In this work, we present the first in-depth characterization of Snakin-3, a member of the subfamily III within the Snakin/GASA gene family of potato. Transient co-expression of Snakin-3 fused to the green fluorescent protein and organelle markers revealed that it is located in the endoplasmic reticulum. Furthermore, expression analyses via pSnakin-3::GUS transgenic plants showed GUS staining mainly in roots and vascular tissues of the stem. Moreover, GUS expression levels were increased after inoculation with Pseudomonas syringae pv. tabaci or Pectobacterium carotovorum subsp. carotovorum and also after auxin treatment mainly in roots and stems. To gain further insights into the function of Snakin-3 in planta, potato overexpressing lines were challenged against P. carotovorum subsp. carotovorum showing enhanced tolerance to this bacterial pathogen. In sum, here we report the first functional characterization of a Snakin/GASA gene from subfamily III in Solanaceae. Our findings provide experimental evidence on promoter activity and subcellular localization and reveal a role of potato Snakin-3 in plant defense.

Protein Sorting Prediction

Many computational methods are available for predicting protein sorting in bacteria. When comparing them, it is important to know that they can be grouped into three fundamentally different approaches: signal-based, global property-based, and homology-based prediction. In this chapter, the strengths and drawbacks of each of these approaches are described through many examples of methods that predict secretion, integration into membranes, or subcellular locations in general. The aim of this chapter is to provide a user-level introduction to the field with a minimum of computational theory.

Unveiling promising immunogenic targets in Coxiella burnetii through in silico analysis: paving the way for novel vaccine strategies

BACKGROUND

Coxiella burnetii, an intracellular pathogen, serves as the causative agent of zoonotic Q fever. This pathogen presents a significant threat due to its potential for airborne transmission, environmental persistence, and pathogenicity. The current whole-cell vaccine (WCV) utilized in Australia to combat Q fever exhibits notable limitations, including severe adverse reactions and limited regulatory approval for human use. This research employed the reverse vaccinology (RV) approach to uncover antigenic proteins and epitopes of C. burnetii, facilitating the development of more potent vaccine candidates.

METHODS

The potential immunogenic proteins derived from C. burnetii RSA493/Nine Mile phase I (NMI) were extracted through manual, automated RV, and virulence factor database (VFDB) methods. Web tools and bioinformatics were used to evaluate physiochemical attributes, subcellular localization, antigenicity, allergenicity, human homology, B-cell epitopes, MHC I and II binding ratios, functional class scores, adhesion probabilities, protein-protein interactions, and molecular docking.

RESULTS

Out of the 1850 proteins encoded by RSA493/NMI, a subset of 178 demonstrated the potential for surface or membrane localization. Following a series of analytical iterations, 14 putative immunogenic proteins emerged. This collection included nine proteins (57.1%) intricately involved in cell wall/membrane/envelope biogenesis processes (CBU_0197 (Q83EW1), CBU_0311 (Q83EK8), CBU_0489 (Q83E43), CBU_0939 (Q83D08), CBU_1190 (P39917), CBU_1829 (Q83AQ2), CBU_1412 (Q83BU0), CBU_1414 (Q83BT8), and CBU_1600 (Q83BB2)). The CBU_1627 (Q83B86 ) (7.1%) implicated in intracellular trafficking, secretion, and vesicular transport, and CBU_0092 (Q83F57) (7.1%) contributing to cell division. Additionally, three proteins (21.4%) displayed uncharacterized functions (CBU_0736 (Q83DJ4), CBU_1095 (Q83CL9), and CBU_2079 (Q83A32)). The congruent results obtained from molecular docking and immune response stimulation lend support to the inclusion of all 14 putative proteins as potential vaccine candidates. Notably, seven proteins with well-defined functions stand out among these candidates.

CONCLUSIONS

The outcomes of this study introduce promising proteins and epitopes for the forthcoming formulation of subunit vaccines against Q fever, with a primary emphasis on cellular processes and the virulence factors of C. burnetii.

Gallionellaceae pangenomic analysis reveals insight into phylogeny, metabolic flexibility, and iron oxidation mechanisms

IMPORTANCE

Neutrophilic iron-oxidizing bacteria (FeOB) produce copious iron (oxyhydr)oxides that can profoundly influence biogeochemical cycles, notably the fate of carbon and many metals. To fully understand environmental microbial iron oxidation, we need a thorough accounting of iron oxidation mechanisms. In this study, we show the Gallionellaceae FeOB genomes encode both characterized iron oxidases as well as uncharacterized multiheme cytochromes (MHCs). MHCs are predicted to transfer electrons from extracellular substrates and likely confer metabolic capabilities that help Gallionellaceae occupy a range of different iron- and mineral-rich niches. Gallionellaceae appear to specialize in iron oxidation, so it would be advantageous for them to have multiple mechanisms to oxidize various forms of iron, given the many iron minerals on Earth, as well as the physiological and kinetic challenges faced by FeOB. The multiple iron/mineral oxidation mechanisms may help drive the widespread ecological success of Gallionellaceae.

Evolutionary analysis and quality assessment of ζ-carbonic anhydrase sequences from environmental microbiome

Carbonic anhydrase (CA) is one of the most vital enzymes in living cells. This study has been performed due to the significance of this metalloenzyme for life and the novelty of some CA families like ζ-CA to evaluate evolutionary processes and quality check their sequences. In this study, bioinformatics methods revealed the presence of ζ-CA in some eukaryotic and prokaryotic microorganisms. Notably, it has not been previously reported in prokaryotes. The coexistence of β- and ζ-CAs in some microorganisms is also a novel finding as well. Also, our analysis identified several CA proteins with 6-14 amino acid intervals between histidine and cysteine in the second highly conserved motif, which can be classified as the novel ζ-CA subfamily members that emerged under the Zn deficiency of aquatic ecosystems and selection pressure in these environments. There is also a possibility that the achieved results are rooted in the contamination of samples from the environmental microbiome genome with genomes of diatom species and the occurrence of errors was observed in the DNA sequencing outcomes. Combining of all results from evolutionary analysis to quality control of ζ-CA DNA sequences is the incentive motivation to explore more the hidden aspects of ζ-CAs.

Hypothetical Proteins of Mycoplasma synoviae Reannotation and Expression Changes Identified via RNA-Sequencing

Mycoplasma synoviae infection rates in chickens are increasing worldwide. Genomic studies have considerably improved our understanding of M. synoviae biology and virulence. However, approximately 20% of the predicted proteins have unknown functions. In particular, the M. synoviae ATCC 25204 genome has 663 encoding DNA sequences, among which 155 are considered encoding hypothetical proteins (HPs). Several of these genes may encode unknown virulence factors. This study aims to reannotate all 155 proteins in M. synoviae ATCC 25204 to predict new potential virulence factors using currently available databases and bioinformatics tools. Finally, 125 proteins were reannotated, including enzymes (39%), lipoproteins (10%), DNA-binding proteins (6%), phase-variable hemagglutinin (19%), and other protein types (26%). Among 155 proteins, 28 proteins associated with virulence were detected, five of which were reannotated. Furthermore, HP expression was compared before and after the M. synoviae infection of cells to identify potential virulence-related proteins. The expression of 14 HP genes was upregulated, including that of five virulence-related genes. Our study improved the functional annotation of M. synoviae ATCC 25204 from 76% to 95% and enabled the discovery of potential virulence factors in the genome. Moreover, 14 proteins that may be involved in M. synoviae infection were identified, providing candidate proteins and facilitating the exploration of the infection mechanism of M. synoviae.

Autophagy-mediated surveillance of Rim4-mRNA interaction safeguards programmed meiotic translation

In yeast meiosis, autophagy is active and essential. Here, we investigate the fate of Rim4, a meiosis-specific RNA-binding protein (RBP), and its associated transcripts during meiotic autophagy. We demonstrate that Rim4 employs a nuclear localization signal (NLS) to enter the nucleus, where it loads its mRNA substrates before nuclear export. Upon reaching the cytoplasm, active autophagy selectively spares the Rim4-mRNA complex. During meiotic divisions, autophagy preferentially degrades Rim4 in an Atg11-dependent manner, coinciding with the release of Rim4-bound mRNAs for translation. Intriguingly, these released mRNAs also become vulnerable to autophagy. In vitro, purified Rim4 and its RRM-motif-containing variants activate Atg1 kinase in meiotic cell lysates and in immunoprecipitated (IP) Atg1 complexes. This suggests that the conserved RNA recognition motifs (RRMs) of Rim4 are involved in stimulating Atg1 and thereby facilitating selective autophagy. Taken together, our findings indicate that autophagy surveils Rim4-mRNA interaction to ensure stage-specific translation during meiosis.

The first multi-tissue genome-scale metabolic model of a woody plant highlights suberin biosynthesis pathways in Quercus suber

Over the last decade, genome-scale metabolic models have been increasingly used to study plant metabolic behaviour at the tissue and multi-tissue level under different environmental conditions. Quercus suber, also known as the cork oak tree, is one of the most important forest communities of the Mediterranean/Iberian region. In this work, we present the genome-scale metabolic model of the Q. suber (iEC7871). The metabolic model comprises 7871 genes, 6231 reactions, and 6481 metabolites across eight compartments. Transcriptomics data was integrated into the model to obtain tissue-specific models for the leaf, inner bark, and phellogen, with specific biomass compositions. The tissue-specific models were merged into a diel multi-tissue metabolic model to predict interactions among the three tissues at the light and dark phases. The metabolic models were also used to analyse the pathways associated with the synthesis of suberin monomers, namely the acyl-lipids, phenylpropanoids, isoprenoids, and flavonoids production. The models developed in this work provide a systematic overview of the metabolism of Q. suber, including its secondary metabolism pathways and cork formation.

A ligand-receptor interactome atlas of the zebrafish

Studies in zebrafish can unravel the functions of cellular communication and thus identify novel bench-to-bedside drugs targeting cellular communication signaling molecules. Due to the incomplete annotation of zebrafish proteome, the knowledge of zebrafish receptors, ligands, and tools to explore their interactome is limited. To address this gap, we de novo predicted the cellular localization of zebrafish reference proteome using deep learning algorithm. We combined the predicted and existing annotations on cellular localization of zebrafish proteins and created repositories of zebrafish ligands, membrane receptome, and interactome as well as associated diseases and targeting drugs. Unlike other tools, our interactome atlas is based on both the physical interaction data of zebrafish proteome and existing human ligand-receptor pair databases. The resources are available as R and Python scripts. DanioTalk provides a novel resource for researchers interested in targeting cellular communication in zebrafish, as we demonstrate in applications studying synapse and axo-glial interactome. DanioTalk methodology can be applied to build and explore the ligand-receptor atlas of other non-mammalian model organisms.

A systematic in silico report on iron and zinc proteome of Zea mays

Zea mays is an essential staple food crop across the globe. Maize contains macro and micronutrients but is limited in essential mineral micronutrients such as Fe and Zn. Worldwide, serious health concerns have risen due to the deficiencies of essential nutrients in human diets, which rigorously jeopardizes economic development. In the present study, the systematic in silico approach has been used to predict Fe and Zn binding proteins from the whole proteome of maize. A total of 356 and 546 putative proteins have been predicted, which contain sequence and structural motifs for Fe and Zn ions, respectively. Furthermore, the functional annotation of these predicted proteins, based on their domains, subcellular localization, gene ontology, and literature support, showed their roles in distinct cellular and biological processes, such as metabolism, gene expression and regulation, transport, stress response, protein folding, and proteolysis. The versatile roles of these shortlisted putative Fe and Zn binding proteins of maize could be used to manipulate many facets of maize physiology. Moreover, in the future, the predicted Fe and Zn binding proteins may act as relevant, novel, and economical markers for various crop improvement programs.

In Planta Study Localizes an Effector Candidate from Austropuccinia psidii Strain MF-1 to the Nucleus and Demonstrates In Vitro Cuticular Wax-Dependent Differential Expression

Austropuccinia psidii is a biotrophic fungus that causes myrtle rust. First described in Brazil, it has since spread to become a globally important pathogen that infects more than 480 myrtaceous species. One of the most important commercial crops affected by A. psidii is eucalypt, a widely grown forestry tree. The A. psidii-Eucalyptus spp. interaction is poorly understood, but pathogenesis is likely driven by pathogen-secreted effector molecules. Here, we identified and characterized a total of 255 virulence effector candidates using a genome assembly of A. psidii strain MF-1, which was recovered from Eucalyptus grandis in Brazil. We show that the expression of seven effector candidate genes is modulated by cell wax from leaves sourced from resistant and susceptible hosts. Two effector candidates with different subcellular localization predictions, and with specific gene expression profiles, were transiently expressed with GFP-fusions in Nicotiana benthamiana leaves. Interestingly, we observed the accumulation of an effector candidate, Ap28303, which was upregulated under cell wax from rust susceptible E. grandis and described as a peptidase inhibitor I9 domain-containing protein in the nucleus. This was in accordance with in silico analyses. Few studies have characterized nuclear effectors. Our findings open new perspectives on the study of A. psidii-Eucalyptus interactions by providing a potential entry point to understand how the pathogen manipulates its hosts in modulating physiology, structure, or function with effector proteins.

Molecular characterization suggests kinetic modulation of expression of accessory viral protein, W, in Newcastle disease virus infected DF1 cells

Viruses adopt strategies to efficiently utilize their compact genome. Members of the family Paramyxoviridae, exhibit a cotranscriptional RNA editing mechanism wherein polymerase stuttering generates accessory proteins from Phosphoprotein (P) gene. Newcastle disease virus (NDV), an avian paramyxovirus, expresses two accessory proteins, V and W, by RNA editing. While P and V proteins are well studied, very little is known about W protein. Recent studies confirmed W protein expression in NDV and the unique subcellular localization of W proteins of virulent and avirulent NDV. We characterized the W protein of NDV strain Komarov, a moderately virulent vaccine strain. W mRNA expression ranged between 7 and 9% of total P gene transcripts similar to virulent NDV. However, W protein expression, detectable by 6 h, peaked at 24 h and dropped by 48 h post infection in DF1 cells indicating a kinetically regulated expression by the virus. The W protein localized in the nucleus and by mutations, a strong nuclear localization signal was identified in the C-terminal region of W protein. The viral growth kinetics study suggested neither supplementation of W protein nor subcellular localization pattern of the supplemented W protein influenced viral replication in vitro similar to that noticed in avirulent NDV. A cytoplasmic mutant of W protein localized in cytoplasm unlike specific mitochondrial colocalization as recorded in velogenic NDV strain SG10 indicating a possible role of W protein in determining the viral pathogenicity. This study describes for the first time, the distinct features of W protein of moderately virulent NDV.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13337-023-00813-2.

Gallionellaceae pangenomic analysis reveals insight into phylogeny, metabolic flexibility, and iron oxidation mechanisms

The iron-oxidizing Gallionellaceae drive a wide variety of biogeochemical cycles through their metabolisms and biominerals. To better understand the environmental impacts of Gallionellaceae, we need to improve our knowledge of their diversity and metabolisms, especially any novel iron oxidation mechanisms. Here, we used a pangenomic analysis of 103 genomes to resolve Gallionellaceae phylogeny and explore the range of genomic potential. Using a concatenated ribosomal protein tree and key gene patterns, we determined Gallionellaceae has four genera, divided into two groups-iron-oxidizing bacteria (FeOB) Gallionella, Sideroxydans, and Ferriphaselus with known iron oxidases (Cyc2, MtoA) and nitrite-oxidizing bacteria (NOB) Candidatus Nitrotoga with nitrite oxidase (Nxr). The FeOB and NOB have similar electron transport chains, including genes for reverse electron transport and carbon fixation. Auxiliary energy metabolisms including S oxidation, denitrification, and organotrophy were scattered throughout the Gallionellaceae FeOB. Within FeOB, we found genes that may represent adaptations for iron oxidation, including a variety of extracellular electron uptake (EEU) mechanisms. FeOB genomes encoded more predicted c-type cytochromes overall, notably more multiheme c-type cytochromes (MHCs) with >10 CXXCH motifs. These include homologs of several predicted outer membrane porin-MHC complexes, including MtoAB and Uet. MHCs are known to efficiently conduct electrons across longer distances and function across a wide range of redox potentials that overlap with mineral redox potentials, which can help expand the range of usable iron substrates. Overall, the results of pangenome analyses suggest that the Gallionellaceae genera Gallionella, Sideroxydans, and Ferriphaselus are primarily iron oxidizers, capable of oxidizing dissolved Fe2+ as well as a range of solid iron or other mineral substrates.

Comparison of PB1-F2 Proximity Interactomes Reveals Functional Differences between a Human and an Avian Influenza Virus

Most influenza viruses express the PB1-F2 protein which is regarded as a virulence factor. However, PB1-F2 behaves differently in avian and mammalian hosts, suggesting that this protein may be involved in the species barrier crossings regularly observed in influenza viruses. To better understand the functions associated with this viral protein, we decided to compare the BioID2-derived proximity interactome of a human PB1-F2 from an H3N2 virus with that of an avian PB1-F2 from an H7N1 strain. The results obtained reveal that the two proteins share only a few interactors and thus common functions. The human virus protein is mainly involved in signaling by Rho GTPases while the avian virus protein is mainly involved in ribonucleoprotein complex biogenesis. PB1-F2 H3N2 interactors include several members of the 14-3-3 protein family, a family of regulatory proteins involved in many signaling pathways. We then validated the interaction with 14-3-3 proteins and were able to show that the association of H3N2-PB1-F2 with YWHAH increased the activity of the antiviral sensor MDA5, while H7N1-PB1-F2 had no effect. Collectively, these results show that PB1-F2 can associate with a large range of protein complexes and exert a wide variety of functions. Furthermore, PB1-F2 interactome differs according to the avian or human origin of the protein.

How Many Sirtuin Genes Are Out There? Evolution of Sirtuin Genes in Vertebrates With a Description of a New Family Member

Studying the evolutionary history of gene families is a challenging and exciting task with a wide range of implications. In addition to exploring fundamental questions about the origin and evolution of genes, disentangling their evolution is also critical to those who do functional/structural studies to allow a deeper and more precise interpretation of their results in an evolutionary context. The sirtuin gene family is a group of genes that are involved in a variety of biological functions mostly related to aging. Their duplicative history is an open question, as well as the definition of the repertoire of sirtuin genes among vertebrates. Our results show a well-resolved phylogeny that represents an improvement in our understanding of the duplicative history of the sirtuin gene family. We identified a new sirtuin gene family member (SIRT3.2) that was apparently lost in the last common ancestor of amniotes but retained in all other groups of jawed vertebrates. According to our experimental analyses, elephant shark SIRT3.2 protein is located in mitochondria, the overexpression of which leads to an increase in cellular levels of ATP. Moreover, in vitro analysis demonstrated that it has deacetylase activity being modulated in a similar way to mammalian SIRT3. Our results indicate that there are at least eight sirtuin paralogs among vertebrates and that all of them can be traced back to the last common ancestor of the group that existed between 676 and 615 millions of years ago.

Regulation of Cell Membrane Permeability Enhanced the Non-Classical Secretion of γ-Cyclodextrin Glycosyltransferase in Bacillus subtilis

γ-Cyclodextrin glycosyltransferase (γ-CGTase, EC 2.4.1.19) is an essential enzyme required in the production of γ-cyclodextrin, which shows huge prospects in the food, medicine, materials, and chemical industries. In this study, γ-CGTase from Bacillus sp. G-825-6 STB17 was successfully cloned and expressed in Bacillus subtilis WB600. The final extracellular activity of γ-CGTase can reach 45.34 U/mL with the deletion of the signal peptide, which was about 11.3 times of the initial level of γ-CGTase secreted by the general pathway. By monitoring the whole cultivation process, secretion was divided into two stages, which were dominated by cell membrane changes and apoptosis. The measurement of lactate dehydrogenase and the results of fluorescence microscopy demonstrated that the cell membrane permeability changed significantly in the middle stage of fermentation, proving that it played a crucial role in the non-classical secretion of γ-CGTase. Furthermore, the addition of Triton X-100, a non-ionic surfactant, remarkably enhanced the secretion level of γ-CGTase by 21.5%, which was caused by the increase in cell membrane permeability. This work is the first to obtain the extracellular expression of CGTase via the non-classical secretion pathway in B. subtilis and provides a new strategy to enhance extracellular expression by regulating the cell membranes.

Genome-wide characterization of FK506-binding proteins, parvulins and phospho-tyrosyl phosphatase activators in wheat and their regulation by heat stress

Peptidyl-prolyl cis-trans isomerases (PPIases) are ubiquitous proteins which are essential for cis-trans isomerisation of peptide bonds preceding the proline residue. PPIases are categorized into four sub-families viz., cyclophilins, FK506-binding proteins (FKBPs), parvulins and protein phosphatase 2A phosphatase activators (PTPAs). Apart from catalysing the cis-trans isomerization, these proteins have also been implicated in diverse cellular functions. Though PPIases have been identified in several important crop plants, information on these proteins, except cyclophilins, is scanty in wheat. In order to understand the role of these genes in wheat, we carried out genome-wide identification using computational approaches. The present study resulted in identification of 71 FKBP (TaFKBP) 12 parvulin (TaPar) and 3 PTPA (TaPTPA) genes in hexaploid wheat genome, which are distributed on different chromosomes with uneven gene densities. The TaFKBP and TaPar proteins, besides PPIase domain, also contain additional domains, indicating functional diversification. In silico prediction also revealed that TaFKBPs are localized to ER, nucleus, chloroplast and cytoplasm, while the TaPars are confined to cytoplasm and nucleus. The TaPTPAs, on the contrary, appear to be present only in the cytoplasm. Evolutionary studies predicted that most of the TaFKBP, TaPar and TaPTPA genes in hexaploid wheat have been derived from their progenitor species, with some events of loss or gain. Syntenic analysis revealed the presence of many collinear blocks of TaFKBP genes in wheat and its sub-genome donors. qRT-PCR analysis demonstrated that expression of TaFKBP and TaPar genes is regulated differentially by heat stress, suggesting their likely involvement in thermotolerance. The findings of this study will provide basis for further functional characterization of these genes and their likely applications in crop improvement.

Lipidome analysis and characterization of Buglossoides arvensis acyltransferases that incorporate polyunsaturated fatty acids into triacylglycerols

Buglossoides arvensis is a burgeoning oilseed crop that contains an unique combination of ω-3 and ω-6 polyunsaturated fatty acids (PUFA), constituting ~80-85% of seed triacylglycerols (TAGs). To uncover the critical TAG biosynthetic pathways contributing for high PUFA accumulation, we performed lipidome of developing seeds and characterized acyltransferases involved in the final step of TAG biosynthesis. During seed development, distribution of lipid molecular species in individual lipid classes showed distinct patterns from an early-stage (6 days after flowering (DAF)) to the middle-stage (12 and 18 DAF) of oil biosynthesis. PUFA-containing TAG species drastically increased from 6 to 12 DAF. The expression profiles of key triacylglycerol biosynthesis genes and patterns of phosphatidylcholine, diacylglycerol and triacylglycerol molecular species during seed development were used to predict the contribution of diacylglycerol acyltransferases (DGAT1 and DGAT2) and phospholipid: diacylglycerol acyltransferases (PDAT1 and PDAT2) to PUFA-rich TAG biosynthesis. Our analysis suggests that DGATs play a crucial role in enriching TAGs with PUFA compared to PDATs. This was further confirmed by fatty acid feeding studies in yeast expressing acyltransferases. BaDGAT2 preferentially incorporated high amounts of PUFAs into TAG, compared to BaDGAT1. Our results provide insight into the molecular mechanisms of TAG accumulation in this plant and identify target genes for transgenic production of SDA in traditional oilseed crops.

Challenges in the Heterologous Production of Furanocoumarins in Escherichia coli

Coumarins and furanocoumarins are plant secondary metabolites with known biological activities. As they are present in low amounts in plants, their heterologous production emerged as a more sustainable and efficient approach to plant extraction. Although coumarins biosynthesis has been positively established, furanocoumarin biosynthesis has been far more challenging. This study aims to evaluate if Escherichia coli could be a suitable host for furanocoumarin biosynthesis. The biosynthetic pathway for coumarins biosynthesis in E. coli was effectively constructed, leading to the production of umbelliferone, esculetin and scopoletin (128.7, 17.6, and 15.7 µM, respectively, from tyrosine). However, it was not possible to complete the pathway with the enzymes that ultimately lead to furanocoumarins production. Prenyltransferase, psoralen synthase, and marmesin synthase did not show any activity when expressed in E. coli. Several strategies were tested to improve the enzymes solubility and activity with no success, including removing potential N-terminal transit peptides and expression of cytochrome P450 reductases, chaperones and/or enzymes to increase dimethylallylpyrophosphate availability. Considering the results herein obtained, E. coli does not seem to be an appropriate host to express these enzymes. However, new alternative microbial enzymes may be a suitable option for reconstituting the furanocoumarins pathway in E. coli. Nevertheless, until further microbial enzymes are identified, Saccharomyces cerevisiae may be considered a preferred host as it has already been proven to successfully express some of these plant enzymes.

Genome-Wide Identification of the A20/AN1 Zinc Finger Protein Family Genes in Ipomoea batatas and Its Two Relatives and Function Analysis of IbSAP16 in Salinity Tolerance

Stress-associated protein (SAP) genes—encoding A20/AN1 zinc-finger domain-containing proteins—play pivotal roles in regulating stress responses, growth, and development in plants. They are considered suitable candidates to improve abiotic stress tolerance in plants. However, the SAP gene family in sweetpotato (Ipomoea batatas) and its relatives is yet to be investigated. In this study, 20 SAPs in sweetpotato, and 23 and 26 SAPs in its wild diploid relatives Ipomoea triloba and Ipomoea trifida were identified. The chromosome locations, gene structures, protein physiological properties, conserved domains, and phylogenetic relationships of these SAPs were analyzed systematically. Binding motif analysis of IbSAPs indicated that hormone and stress responsive cis-acting elements were distributed in their promoters. RT-qPCR or RNA-seq data revealed that the expression patterns of IbSAP, ItbSAP, and ItfSAP genes varied in different organs and responded to salinity, drought, or ABA (abscisic acid) treatments differently. Moreover, we found that IbSAP16 driven by the 35 S promoter conferred salinity tolerance in transgenic Arabidopsis. These results provided a genome-wide characterization of SAP genes in sweetpotato and its two relatives and suggested that IbSAP16 is involved in salinity stress responses. Our research laid the groundwork for studying SAP-mediated stress response mechanisms in sweetpotato.

Proteome and morphological analysis show unexpected differences between promastigotes of Leishmania amazonensis PH8 and LV79 strains

BACKGROUND

Leishmaniases are diseases caused by Leishmania protozoans that affect around 12 million people. Leishmania promastigotes are transmitted to vertebrates by female phlebotomine flies during their blood meal. Parasites attach to phagocytic cells, are phagocytosed and differentiate into amastigotes. We previously showed that PH8 and LV79 strains of Leishmania amazonensis have different virulence in mice and that their amastigotes differ in their proteomes. In this work, we compare promastigotes' infectivity in macrophages, their proteomes and morphologies.

METHODS/PRINCIPAL FINDINGS

Phagocytosis assays showed that promastigotes adhesion to and phagocytosis by macrophages is higher in PH8 than LV79. To identify proteins that differ between the two strains and that may eventually contribute for these differences we used a label-free proteomic approach to compare promastigote´s membrane-enriched fractions. Proteomic analysis enabled precise discrimination of PH8 and LV79 protein profiles and the identification of several differentially abundant proteins. The proteins more abundant in LV79 promastigotes participate mainly in translation and amino acid and nucleotide metabolism, while the more abundant in PH8 are involved in carbohydrate metabolism, cytoskeleton composition and vesicle/membrane trafficking. Interestingly, although the virulence factor GP63 was more abundant in the less virulent LV79 strain, zymography suggests a higher protease activity in PH8. Enolase, which may be related to virulence, was more abundant in PH8 promastigotes. Unexpectedly, flow cytometry and morphometric analysis indicate higher abundance of metacyclics in LV79.

CONCLUSIONS/SIGNIFICANCE

Proteome comparison of PH8 and LV79 promastigotes generated a list of differential proteins, some of which may be further prospected to affect the infectivity of promastigotes. Although proteomic profile of PH8 includes more proteins characteristic of metacyclics, flow cytometry and morphometric analysis indicate a higher abundance of metacyclics in LV79 cultures. These results shed light to the gaps in our knowledge of metacyclogenesis in L. amazonensis, and to proteins that should be studied in the context of infection by this species.

Mechanistic Insight into the Enzymatic Inhibition of β-Amyrin against Mycobacterial Rv1636: In Silico and In Vitro Approaches

Simple Summary

Rv1636 is a mycobacterial universal stress protein whose expression level increases in different type of stress conditions. This protein promotes the growth of Mycobacterium tuberculosis in the host derived stress conditions generated during infection. Therefore in this manuscipt, we are trying to target Rv1636 using natural inhibitor. Targeting essential Mycobacterial protein using natural prodect was hypothesized to generate a molecule with low toxic effects and high inhibitory activity. It was found that Rv1636 contains ATPase activity and its ATPase activity gets disturbed by addition of β-Amyrin in the reaction. β-Amyrin was forund to interfere with the ATP binding site of Rv1636 which was confirmed by molecular docking anad dynamic studies. In addition to the ATPase activity, Rv1636 was also contain the cAMP binding capacity and also involved in balancing the cAMP levels inside cells. So, targeting Rv1636 using β-Amyrin disrupts its ATPase activity and cAMP regulatory activity and these conditions might make Mycobacterium tuberculosis more susceptible to the host derived stress conditions.

Abstract

Mycobacterium tuberculosis has seen tremendous success as it has developed defenses to reside in host alveoli despite various host-related stress circumstances. Rv1636 is a universal stress protein contributing to mycobacterial survival in different host-derived stress conditions. Both ATP and cAMP can be bound with the Rv1636, and their binding actions are independent of one another. β-Amyrin, a triterpenoid compound, is abundant in medicinal plants and has many pharmacological properties and broad therapeutic potential. The current study uses biochemical, biophysical, and computational methods to define the binding of Rv1636 with β-Amyrin. A substantial interaction between β-Amyrin and Rv1636 was discovered by molecular docking studies, which helped decipher the critical residues involved in the binding process. VAL60 is a crucial residue found in the complexes of both Rv1636_β-Amyrin and Rv1636-ATP. Additionally, the Rv1636_β-Amyrin complex was shown to be stable by molecular dynamics simulation studies (MD), with minimal changes observed during the simulation. In silico observations were further complemented by in vitro assays. Successful cloning, expression, and purification of Rv1636 were accomplished using Ni-NTA affinity chromatography. The results of the ATPase activity assay indicated that Rv1636’s ATPase activity was inhibited in the presence of various β-Amyrin concentrations. Additionally, circular dichroism spectroscopy (CD) was used to examine modifications to Rv1636 secondary structure upon binding of β-Amyrin. Finally, isothermal titration calorimetry (ITC) advocated spontaneous binding of β-Amyrin with Rv1636 elucidating the thermodynamics of the Rv1636_β-Amyrin complex. Thus, the study establishes that β-Amyrin binds to Rv1636 with a significant affinity forming a stable complex and inhibiting its ATPase activity. The present study suggests that β-Amyrin might affect the functioning of Rv1636, which makes the bacterium vulnerable to different stress conditions.

Upregulation of protein N-glycosylation plays crucial roles in the response of Camellia sinensis leaves to fluoride

The tea plant (Camellia sinensis) is one of the three major beverage crops in the world with its leaves consumption as tea. However, it can hyperaccumulate fluoride with about 98% fluoride deposition in the leaves. Our previously studies found that cell wall proteins (CWPs) might play a central role in fluoride accumulation/detoxification in C. sinensis. CWP is known to be glycosylated, however the response of CWP N-glycosylation to fluoride remains unknown in C. sinensis. In this study, a comparative N-glycoproteomic analysis was performed through HILIC enrichment coupled with UPLC-MS/MS based on TMT-labeling approach in C. sinensis leaves. Totally, 237 N-glycoproteins containing 326 unique N-glycosites were identified. 73.4%, 18.6%, 6.3% and 1.7% of these proteins possess 1, 2, 3, and ≥4 modification site, respectively. 93.2% of these proteins were predicted to be localized in the secretory pathway and 78.9% of them were targeted to the cell wall and the plasma membrane. 133 differentially accumulated N-glycosites (DNGSs) on 100 N-glycoproteins (DNGPs) were detected and 85.0% of them exhibited upregulated expression after fluoride treatment. 78.0% DNGPs were extracellular DNGPs, which belonged to CWPs, and 53.0% of them were grouped into protein acting on cell wall polysaccharides, proteases and oxido-reductases, whereas the majority of the remaining DNGPs were mainly related to N-glycoprotein biosynthesis, trafficking and quality control. Our study shed new light on the N-glycoproteome study, and revealed that increased N-glycosylation abundance of CWPs might contribute to fluoride accumulation/detoxification in C. sinensis leave.

Growing Maize Root: Lectins Involved in Consecutive Stages of Cell Development

Proteins that carry specific carbohydrate-binding lectin domains have a great variety and are ubiquitous across the plant kingdom. In turn, the plant cell wall has a complex carbohydrate composition, which is subjected to constant changes in the course of plant development. In this regard, proteins with lectin domains are of great interest in the context of studying their contribution to the tuning and monitoring of the cell wall during its modifications in the course of plant organ development. We performed a genome-wide screening of lectin motifs in the Zea mays genome and analyzed the transcriptomic data from five zones of primary maize root with cells at different development stages. This allowed us to obtain 306 gene sequences encoding putative lectins and to relate their expressions to the stages of root cell development and peculiarities of cell wall metabolism. Among the lectins whose expression was high and differentially regulated in growing maize root were the members of the EUL, dirigent–jacalin, malectin, malectin-like, GNA and Nictaba families, many of which are predicted as cell wall proteins or lectin receptor-like kinases that have direct access to the cell wall. Thus, a set of molecular players was identified with high potential to play important roles in the early stages of root morphogenesis.

PPVED: A machine learning tool for predicting the effect of single amino acid substitution on protein function in plants

Single amino acid substitution (SAAS) produces the most common variant of protein function change under physiological conditions. As the number of SAAS events in plants has increased exponentially, an effective prediction tool is required to help identify and distinguish functional SAASs from the whole genome as either potentially causal traits or as variants. Here, we constructed a plant SAAS database that stores 12 865 SAASs in 6172 proteins and developed a tool called Plant Protein Variation Effect Detector (PPVED) that predicts the effect of SAASs on protein function in plants. PPVED achieved an 87% predictive accuracy when applied to plant SAASs, an accuracy that was much higher than those from six human database software: SIFT, PROVEAN, PANTHER-PSEP, PhD-SNP, PolyPhen-2, and MutPred2. The predictive effect of six SAASs from three proteins in Arabidopsis and maize was validated with wet lab experiments, of which five substitution sites were accurately predicted. PPVED could facilitate the identification and characterization of genetic variants that explain observed phenotype variations in plants, contributing to solutions for challenges in functional genomics and systems biology. PPVED can be accessed under a CC-BY (4.0) license via http://www.ppved.org.cn.

merlin, an improved framework for the reconstruction of high-quality genome-scale metabolic models

Genome-scale metabolic models have been recognised as useful tools for better understanding living organisms' metabolism. merlin (https://www.merlin-sysbio.org/) is an open-source and user-friendly resource that hastens the models' reconstruction process, conjugating manual and automatic procedures, while leveraging the user's expertise with a curation-oriented graphical interface. An updated and redesigned version of merlin is herein presented. Since 2015, several features have been implemented in merlin, along with deep changes in the software architecture, operational flow, and graphical interface. The current version (4.0) includes the implementation of novel algorithms and third-party tools for genome functional annotation, draft assembly, model refinement, and curation. Such updates increased the user base, resulting in multiple published works, including genome metabolic (re-)annotations and model reconstructions of multiple (lower and higher) eukaryotes and prokaryotes. merlin version 4.0 is the only tool able to perform template based and de novo draft reconstructions, while achieving competitive performance compared to state-of-the art tools both for well and less-studied organisms.

mebipred: identifying metal binding potential in protein sequence

Motivation

metal-binding proteins have a central role in maintaining life processes. Nearly one-third of known protein structures contain metal ions that are used for a variety of needs, such as catalysis, DNA/RNA binding, protein structure stability, etc. Identifying metal-binding proteins is thus crucial for understanding the mechanisms of cellular activity. However, experimental annotation of protein metal-binding potential is severely lacking, while computational techniques are often imprecise and of limited applicability.

Results

we developed a novel machine learning-based method, mebipred, for identifying metal-binding proteins from sequence-derived features. This method is over 80% accurate in recognizing proteins that bind metal ion-containing ligands; the specific identity of 11 ubiquitously present metal ions can also be annotated. mebipred is reference-free, i.e. no sequence alignments are involved, and is thus faster than alignment-based methods; it is also more accurate than other sequence-based prediction methods. Additionally, mebipred can identify protein metal-binding capabilities from short sequence stretches, e.g. translated sequencing reads, and, thus, may be useful for the annotation of metal requirements of metagenomic samples. We performed an analysis of available microbiome data and found that ocean, hot spring sediments and soil microbiomes use a more diverse set of metals than human host-related ones. For human microbiomes, physiological conditions explain the observed metal preferences. Similarly, subtle changes in ocean sample ion concentration affect the abundance of relevant metal-binding proteins. These results highlight mebipred’s utility in analyzing microbiome metal requirements.

Availability and implementation

mebipred is available as a web server at services.bromberglab.org/mebipred and as a standalone package at https://pypi.org/project/mymetal/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Sirt5 Inhibits BmNPV Replication by Promoting a Relish-Mediated Antiviral Pathway in Bombyx mori

Silent information regulators (Sirtuins) belong to the family of nicotinamide adenine dinucleotide (NAD⁺)-dependent histone deacetylases (HDACs) that have diverse functions in cells. Mammalian Sirtuins have seven isoforms (Sirt1–7) which have been found to play a role in viral replication. However, Sirtuin members of insects are very different from mammals, and the function of insect Sirtuins in regulating virus replication is unclear. The silkworm, Bombyx mori, as a model species of Lepidoptera, is also an important economical insect. B. mori nucleopolyhedrovirus (BmNPV) is a major pathogen that specifically infects silkworms and causes serious losses in the sericulture industry. Here, we used the infection of the silkworm by BmNPV as a model to explore the effect of Sirtuins on virus replication. We initially knocked down all silkworm Sirtuins, and then infected with BmNPV to analyze its replication. Sirt2 and Sirt5 were found to have potential antiviral functions in the silkworm. We further confirmed the antiviral function of silkworm Sirt5 through its effects on viral titers during both knockdown and overexpression experiments. Additionally, Suramin, a Sirt5 inhibitor, was found to promote BmNPV replication. In terms of molecular mechanism, it was found that silkworm Sirt5 might promote the immune pathway mediated by Relish, thereby enhancing the host antiviral response. This study is the first to explore the role of Sirtuins in insect-virus interactions, providing new insights into the functional role of members of the insect Sirtuin family.

The grapevine aspartic protease gene family: characterization and expression modulation in response to Plasmopara viticola

Grapevine aspartic proteases gene family is characterized and five VviAPs appear to be involved in grapevine defense against downy mildew. Grapevine (Vitis vinifera L.) is one of the most important crops worldwide. However, it is highly susceptible to the downy mildew disease caused by Plasmopara viticola (Berk. & Curt.) Berl. & De Toni. To minimize the use of fungicides used to control P. viticola, it is essential to gain a deeper comprehension on this pathosystem and proteases have gained particular interest in the past decade. Proteases were shown to actively participate in plant-pathogen interactions, not only in the processes that lead to plant cell death, stress responses and protein processing/degradation but also as components of the recognition and signalling pathways. The aim of this study was to identify and characterize the aspartic proteases (APs) involvement in grapevine defense against P. viticola. A genome-wide search and bioinformatics characterization of the V. vinifera AP gene family was conducted and a total of 81 APs proteins, coded by 65 genes, were found. VviAPs proteins can be divided into three categories, similar to those previously described for other plants. Twelve APs coding genes were selected, and expression analysis was conducted at several time-points after inoculation in both compatible and incompatible interactions. Five grapevine APs may be involved in grapevine tolerance against P. viticola. Our findings provide an overall understanding of the VviAPs gene family and establish better groundwork to further describe the roles of VviAPs in defense against P. viticola.

Understanding Diversity, Evolution, and Structure of Small Heat Shock Proteins in Annelida Through in Silico Analyses

Small heat shock proteins (sHsps) are oligomeric stress proteins characterized by an α-crystallin domain (ACD). These proteins are localized in different subcellular compartments and play critical roles in the stress physiology of tissues, organs, and whole multicellular eukaryotes. They are ubiquitous proteins found in all living organisms, from bacteria to mammals, but they have never been studied in annelids. Here, a data set of 23 species spanning the annelid tree of life, including mostly transcriptomes but also two genomes, was interrogated and 228 novel putative sHsps were identified and manually curated. The analysis revealed very high protein diversity and showed that a significant number of sHsps have a particular dimeric architecture consisting of two tandemly repeated ACDs. The phylogenetic analysis distinguished three main clusters, two of them containing both monomeric sHsps, and ACDs located downstream in the dimeric sHsps, and the other one comprising the upstream ACDs from those dimeric forms. Our results support an evolutionary history of these proteins based on duplication events prior to the Spiralia split. Monomeric sHsps 76) were further divided into five subclusters. Physicochemical properties, subcellular location predictions, and sequence conservation analyses provided insights into the differentiating elements of these putative functional groups. Strikingly, three of those subclusters included sHsps with features typical of metazoans, while the other two presented characteristics resembling non-metazoan proteins. This study provides a solid background for further research on the diversity, evolution, and function in the family of the sHsps. The characterized annelid sHsps are disclosed as essential for improving our understanding of this important family of proteins and their pleotropic functions. The features and the great diversity of annelid sHsps position them as potential powerful molecular biomarkers of environmental stress for acting as prognostic tool in a diverse range of environments.

Plant Cell Wall Proteomes: The Core of Conserved Protein Families and the Case of Non-Canonical Proteins

Plant cell wall proteins (CWPs) play critical roles during plant development and in response to stresses. Proteomics has revealed their great diversity. With nearly 1000 identified CWPs, the Arabidopsis thaliana cell wall proteome is the best described to date and it covers the main plant organs and cell suspension cultures. Other monocot and dicot plants have been studied as well as bryophytes, such as Physcomitrella patens and Marchantia polymorpha. Although these proteomes were obtained using various flowcharts, they can be searched for the presence of members of a given protein family. Thereby, a core cell wall proteome which does not pretend to be exhaustive, yet could be defined. It comprises: (i) glycoside hydrolases and pectin methyl esterases, (ii) class III peroxidases, (iii) Asp, Ser and Cys proteases, (iv) non-specific lipid transfer proteins, (v) fasciclin arabinogalactan proteins, (vi) purple acid phosphatases and (vii) thaumatins. All the conserved CWP families could represent a set of house-keeping CWPs critical for either the maintenance of the basic cell wall functions, allowing immediate response to environmental stresses or both. Besides, the presence of non-canonical proteins devoid of a predicted signal peptide in cell wall proteomes is discussed in relation to the possible existence of alternative secretion pathways.

Identifying Protein Subcellular Locations With Embeddings-Based node2loc

Identifying protein subcellular locations is an important topic in protein function prediction. Interacting proteins may share similar locations. Thus, it is imperative to infer protein subcellular locations by taking protein-protein interactions (PPIs)into account. In this study, we present a network embedding-based method, node2loc, to identify protein subcellular locations. node2loc first learns distributed embeddings of proteins in a protein-protein interaction (PPI)network using node2vec. Then the learned embeddings are further fed into a recurrent neural network (RNN). To resolve the severe class imbalance of different subcellular locations, Synthetic Minority Over-sampling Technique (SMOTE)is applied to artificially synthesize proteins for minority classes. node2loc is evaluated on our constructed human benchmark dataset with 16 subcellular locations and yields a Matthews correlation coefficient (MCC)value of 0.800, which is superior to baseline methods. In addition, node2loc yields a better performance on a Yeast benchmark dataset with 17 locations. The results demonstrate that the learned representations from a PPI network have certain discriminative ability for classifying protein subcellular locations. However, node2loc is a transductive method, it only works for proteins connected in a PPI network, and it needs to be retrained for new proteins. In addition, the PPI network needs be annotated to some extent with location information. node2loc is freely available at https://github.com/xypan1232/node2loc.

Elucidation of complex respiratory chains: a straightforward strategy to monitor electron transfer between cytochromes

Cytochromes are electron transfer proteins essential in various biological systems, playing crucial roles in the respiratory chains of bacteria. These proteins are particularly abundant in electrogenic microorganisms and are responsible for the efficient delivery of electrons to the cells’ exterior. The capability of sending electron outside the cells open new avenues to be explored for emerging biotechnological applications in bioremediation, microbial electrosynthesis and bioenergy fields. To develop these applications, it is critical to identify the different redox partners and elucidate the stepwise electron transfer along the respiratory paths. However, investigating direct electron transfer events between proteins with identical features in nearly all spectroscopic techniques is extremely challenging. NMR spectroscopy offers the possibility to overcome this difficulty by analysing the alterations of the spectral signatures of each protein caused by electron exchange events. The uncrowded NMR spectral regions containing the heme resonances of the cytochromes display unique and distinct signatures in the reduced and oxidized states, which can be explored to monitor electron transfer within the redox complex. In this study, we present a strategy for a fast and straightforward monitorization of electron transfer between c-type cytochromes, using as model a triheme periplasmic cytochrome (PpcA) and a membrane associated monoheme cytochrome (OmcF) from the electrogenic bacterium Geobacter sulfurreducens. The comparison between the 1D 1H NMR spectra obtained for samples containing the two cytochromes and for samples containing the individual proteins clearly demonstrated a unidirectional electron transfer within the redox complex. This strategy provides a simple and straightforward means to elucidate complex biologic respiratory electron transfer chains.

Subtractive genomic approach toward introduction of novel immunogenic targets against Clostridioides difficile: Thinking out of the box

Clostridioides difficile is one of the major causatives of nosocomial infections worldwide. Antibiotic-associated diarrhea, pseudomembranous colitis, and toxic megacolon are the most common forms of C. difficile infection (CDI). Considering the high antibiotic resistance of C. difficile isolates and the low efficacy of immunization with toxin-related vaccines, we suggested that surface-exposed and secreted proteins could be considered as potential immunogenic targets against CDI. Various immuninformatics databases were used to predict antigenicity, allergenicity, B-cell epitopes, MHC-II binding sites, conserved domains, prevalence and conservation of proteins among the most common sequence types, molecular docking, and immunosimulation of immunogenic targets. Finally, 16 proteins belonging to three functional groups were identified, including proteins involved in the cell wall and peptidoglycan layer (nine proteins), flagellar assembly (five proteins), spore germination (one protein), and a protein with unknown function. Molecular docking results showed that among all the mentioned proteins, WP_009892971.1 (Acd) and WP_009890599.1 (a C40 family peptidase) had the strongest interactions with human Toll-like receptor 2 (TLR-2) and TLR-4. This study proposes a combination of C. difficile toxoid (Tcd) and surface-exposed proteins such as Acd as a promising vaccine formulation for protection against circulating clinical strains of C. difficile.

Spontaneous Phage Resistance in Avian Pathogenic Escherichia coli

Avian pathogenic Escherichia coli (APEC) is one of the most important bacterial pathogens affecting poultry worldwide. The emergence of multidrug-resistant pathogens has renewed the interest in the therapeutic use of bacteriophages (phages). However, a major concern for the successful implementation of phage therapy is the emergence of phage-resistant mutants. The understanding of the phage-host interactions, as well as underlying mechanisms of resistance, have shown to be essential for the development of a successful phage therapy. Here, we demonstrate that the strictly lytic Escherichia phage vB_EcoM-P10 rapidly selected for resistance in the APEC ST95 O1 strain AM621. Whole-genome sequence analysis of 109 spontaneous phage-resistant mutant strains revealed 41 mutants with single-nucleotide polymorphisms (SNPs) in their core genome. In 32 of these, a single SNP was detected while two SNPs were identified in a total of nine strains. In total, 34 unique SNPs were detected. In 42 strains, including 18 strains with SNP(s), gene losses spanning 17 different genes were detected. Affected by genetic changes were genes known to be involved in phage resistance (outer membrane protein A, lipopolysaccharide-, O- antigen-, or cell wall-related genes) as well as genes not previously linked to phage resistance, including two hypothetical genes. In several strains, we did not detect any genetic changes. Infecting phages were not able to overcome the phage resistance in host strains. However, interestingly the initial infection was shown to have a great fitness cost for several mutant strains, with up to ∼65% decrease in overall growth. In conclusion, this study provides valuable insights into the phage-host interaction and phage resistance in APEC. Although acquired resistance to phages is frequently observed in pathogenic E. coli, it may be associated with loss of fitness, which could be exploited in phage therapy.

Protein embeddings and deep learning predict binding residues for various ligand classes

One important aspect of protein function is the binding of proteins to ligands, including small molecules, metal ions, and macromolecules such as DNA or RNA. Despite decades of experimental progress many binding sites remain obscure. Here, we proposed bindEmbed21, a method predicting whether a protein residue binds to metal ions, nucleic acids, or small molecules. The Artificial Intelligence (AI)-based method exclusively uses embeddings from the Transformer-based protein Language Model (pLM) ProtT5 as input. Using only single sequences without creating multiple sequence alignments (MSAs), bindEmbed21DL outperformed MSA-based predictions. Combination with homology-based inference increased performance to F1 = 48 ± 3% (95% CI) and MCC = 0.46 ± 0.04 when merging all three ligand classes into one. All results were confirmed by three independent data sets. Focusing on very reliably predicted residues could complement experimental evidence: For the 25% most strongly predicted binding residues, at least 73% were correctly predicted even when ignoring the problem of missing experimental annotations. The new method bindEmbed21 is fast, simple, and broadly applicable-neither using structure nor MSAs. Thereby, it found binding residues in over 42% of all human proteins not otherwise implied in binding and predicted about 6% of all residues as binding to metal ions, nucleic acids, or small molecules.

Reannotation of Mycoplasma hyopneumoniae hypothetical proteins revealed novel potential virulence factors

Mycoplasma hyopneumoniae is a bacterium that inhabits the swine respiratory tract, causing porcine enzootic pneumonia, which generates significant economic losses to the swine industry worldwide. The knowledge on M. hyopneumoniae biology and virulence have been significantly increased by genomics studies. However, around 30% of the predicted proteins remained of unknown function so far. According to the original annotation, the genome of M. hyopneumoniae 7448, a Brazilian pathogenic strain, had 693 coding DNA sequences, 244 of which were annotated as coding for hypothetical or uncharacterized proteins. Among them, there may be still several genes coding for unknown virulence factors. Therefore, this study aimed to functionally reannotate the whole set of 244 M. hyopneumoniae 7448 proteins of unknown function based on currently available database and bioinformatic tools, in order to predict novel potential virulence factors. Predictions of physicochemical properties, subcellular localization, function, overall association to virulence and antigenicity are provided. With that, 159 out of the set of 244 proteins of unknown function had a putative function associated to them, allowing identification of novel enzymes, membrane transporters, lipoproteins, DNA-binding proteins and adhesins. Furthermore, 139 proteins were generally associated to virulence, 14 of which had a function assigned and were differentially expressed between pathogenic and non-pathogenic strains of M. hyopneumoniae. Moreover, all extracellular or cytoplasmic membrane predicted proteins had putative epitopes identified. Overall, these analyses improved the functional annotation of M. hyopneumoniae 7448 genome from 65% to 87% and allowed the identification of new potential virulence factors.

Investigation of novel putative immunogenic targets against Staphylococcus aureus using a reverse vaccinology strategy

BACKGROUND

The emergence of methicillin-resistant Staphylococcus aureus (MRSA) and vancomycin-resistant S. aureus (VRSA) strains is a significant public health concern. Considering the high morbidity and mortality of invasive S. aureus infections and multi-drug resistant strains, there is an urgent need for non-antibiotic immune-based approaches to cure these infections. Despite all efforts, vaccine candidates targeting S. aureus failed in human clinical trials, and no approved vaccine is available against this pathogen. Therefore, this study aimed to introduce suitable candidates for immunization against S. aureus using a comprehensive reverse vaccinology approach.

METHODS

In this study, we retrieved putative immunogenic targets from three different levels (literature review, automated reverse vaccinology, and manual reverse vaccinology) and evaluated them using several immunoinformatics analyses including antigenicity, allergenicity, PSI-BLAST to human proteome, physiochemical properties, B-cell, and T-cell epitopes. In the next step, the quartile method scoring was used to the shortlisted proteins. Finally, the molecular docking and immune simulation of immunogenic targets were performed.

RESULTS

This study presents 12 vaccine candidates, including three enzymatic proteins (WP_000222271.1, WP_001170274, and WP_000827736.1), three cell wall-associated proteins (WP_001125631.1, WP_000731642, and WP_000751265.1), two hemolysins (WP_000594517.1, and WP_000916697.1), one secretion involved protein (WP_000725226.1), one heme‑iron binding protein (WP_001041573.1), one superantigen like protein (WP_000668994.1) and one hypothetical proteins (WP_000737711.1).

CONCLUSION

Through quartile scoring method, immune simulation and molecular docking, four promising targets including lytic transglycosylase IsaA, HlgA, secretory antigen precursor SsaA, and heme uptake protein IsdB were selected as the shortlisted proteins. It seems that a polarized immunization (Th1/Th17) response is needed for protection against this bacterium. An optimized formulation based on these putative immunogenic proteins and a wisely adjuvant selection may drive the immune system toward a full protection.

Ralstonia solanacearum Effectors Localize to Diverse Organelles in Solanum Hosts

Phytopathogenic bacteria secrete type III effector (T3E) proteins directly into host plant cells. T3Es can interact with plant proteins and frequently manipulate plant host physiological or developmental processes. The proper subcellular localization of T3Es is critical for their ability to interact with plant targets, and knowledge of T3E localization can be informative for studies of effector function. Here we investigated the subcellular localization of 19 T3Es from the phytopathogenic bacteria Ralstonia pseudosolanacearum and Ralstonia solanacearum. Approximately 45% of effectors in our library localize to both the plant cell periphery and the nucleus, 15% exclusively to the cell periphery, 15% exclusively to the nucleus, and 25% to other organelles, including tonoplasts and peroxisomes. Using tomato hairy roots, we show that T3E localization is similar in both leaves and roots and is not impacted by Solanum species. We find that in silico prediction programs are frequently inaccurate, highlighting the value of in planta localization experiments. Our data suggest that Ralstonia targets a wide diversity of cellular organelles and provides a foundation for developing testable hypotheses about Ralstonia effector function.

Subcellular Localization of Seed-Expressed LEA_4 Proteins Reveals Liquid-Liquid Phase Separation for LEA9 and for LEA48 Homo- and LEA42-LEA48 Heterodimers

LEA proteins are involved in plant stress tolerance. In Arabidopsis, the LEA_4 Pfam group is the biggest group with the majority of its members being expressed in dry seeds. To assess subcellular localization in vivo, we investigated 11 seed-expressed LEA_4 proteins in embryos dissected from dry seeds expressing LEA_4 fusion proteins under its native promoters with the Venus fluorescent protein (proLEA_4::LEA_4:Venus). LEA_4 proteins were shown to be localized in the endoplasmic reticulum, nucleus, mitochondria, and plastids. LEA9, in addition to the nucleus, was also found in cytoplasmic condensates in dry seeds dependent on cellular hydration level. Most investigated LEA_4 proteins were detected in 4-d-old seedlings. In addition, we assessed bioinformatic tools for predicting subcellular localization and promoter motifs of 11 seed-expressed LEA_4 proteins. Ratiometric bimolecular fluorescence complementation assays showed that LEA7, LEA29, and LEA48 form homodimers while heterodimers were formed between LEA7-LEA29 and LEA42-LEA48 in tobacco leaves. Interestingly, LEA48 homodimers and LEA42-LEA48 heterodimers formed droplets structures with liquid-like behavior. These structures, along with LEA9 cytoplasmic condensates, may have been formed through liquid-liquid phase separation. These findings suggest possible important roles of LLPS for LEA protein functions.

Light attention predicts protein location from the language of life

Summary

Although knowing where a protein functions in a cell is important to characterize biological processes, this information remains unavailable for most known proteins. Machine learning narrows the gap through predictions from expert-designed input features leveraging information from multiple sequence alignments (MSAs) that is resource expensive to generate. Here, we showcased using embeddings from protein language models for competitive localization prediction without MSAs. Our lightweight deep neural network architecture used a softmax weighted aggregation mechanism with linear complexity in sequence length referred to as light attention. The method significantly outperformed the state-of-the-art (SOTA) for 10 localization classes by about 8 percentage points (Q10). So far, this might be the highest improvement of just embeddings over MSAs. Our new test set highlighted the limits of standard static datasets: while inviting new models, they might not suffice to claim improvements over the SOTA.

Availability and implementation

The novel models are available as a web-service at http://embed.protein.properties. Code needed to reproduce results is provided at https://github.com/HannesStark/protein-localization. Predictions for the human proteome are available at https://zenodo.org/record/5047020.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Proteolytic processing of plant proteins by potyvirus NIa proteases

The NIa protease of potyviruses is a chymotrypsin-like cysteine protease related to the picornavirus 3C protease. It is also a multifunctional protein known to play multiple roles during virus infection. Picornavirus 3C proteases cleave hundreds of host proteins to facilitate virus infection. However, whether or not potyvirus NIa proteases cleave plant proteins has so far not been tested. Regular expression search using the cleavage site consensus sequence [EQN]xVxH[QE]/[SGTA] for the plum pox virus (PPV) protease identified 90 to 94 putative cleavage events in the proteomes of Prunus persica (a crop severely affected by PPV), Arabidopsis thaliana, and Nicotiana benthamiana (two experimental hosts). In vitro processing assays confirmed cleavage of six A. thaliana and five P. persica proteins by the PPV protease. These proteins were also cleaved in vitro by the protease of turnip mosaic virus (TuMV), which has a similar specificity. We confirmed in vivo cleavage of a transiently expressed tagged version of AtEML2, an EMSY-like protein belonging to a family of nuclear histone readers known to be involved in pathogen resistance. Cleavage of AtEML2 was efficient and was observed in plants that coexpressed the PPV or TuMV NIa proteases or in plants that were infected with TuMV. We also showed partial in vivo cleavage of AtDUF707, a membrane protein annotated as lysine ketoglutarate reductase trans-splicing protein. Although cleavage of the corresponding endogenous plant proteins remains to be confirmed, the results show that a plant virus protease can cleave host proteins during virus infection and highlight a new layer of plant-virus interactions.

IMPORTANCE Viruses are highly adaptive and use multiple molecular mechanisms to highjack or modify the cellular resources to their advantage. They must also counteract or evade host defense responses. One well-characterized mechanism used by vertebrate viruses is the proteolytic cleavage of host proteins to inhibit the activities of these proteins and/or to produce cleaved protein fragments that are beneficial to the virus infection cycle. Even though almost half of the known plant viruses encode at least one protease, it was not known whether plant viruses employ this strategy. Using an in silico prediction approach and the well-characterized specificity of potyvirus NIa proteases, we were able to identify hundreds of putative cleavage sites in plant proteins, several of which were validated by downstream experiments. It can be anticipated that many other plant virus proteases also cleave host proteins and that the identification of these cleavage events will lead to novel antiviral strategies.