dbTEU: a protein database of trace element utilization

deep-Sep: a deep learning-based method for fast and accurate prediction of selenoprotein genes in bacteria

Selenoproteins are a special group of proteins with major roles in cellular antioxidant defense. They contain the 21st amino acid selenocysteine (Sec) in the active sites, which is encoded by an in-frame UGA codon. Compared to eukaryotes, identification of selenoprotein genes in bacteria remains challenging due to the absence of an effective strategy for distinguishing the Sec-encoding UGA codon from a normal stop signal. In this study, we have developed a deep learning-based algorithm, deep-Sep, for quickly and precisely identifying selenoprotein genes in bacterial genomic sequences. This algorithm uses a Transformer-based neural network architecture to construct an optimal model for detecting Sec-encoding UGA codons and a homology search-based strategy to remove additional false positives. During the training and testing stages, deep-Sep has demonstrated commendable performance, including an F1 score of 0.939 and an area under the receiver operating characteristic curve of 0.987. Furthermore, when applied to 20 bacterial genomes as independent test data sets, deep-Sep exhibited remarkable capability in identifying both known and new selenoprotein genes, which significantly outperforms the existing state-of-the-art method. Our algorithm has proved to be a powerful tool for comprehensively characterizing selenoprotein genes in bacterial genomes, which should not only assist in accurate annotation of selenoprotein genes in genome sequencing projects but also provide new insights for a deeper understanding of the roles of selenium in bacteria.IMPORTANCESelenium is an essential micronutrient present in selenoproteins in the form of Sec, which is a rare amino acid encoded by the opal stop codon UGA. Identification of all selenoproteins is of vital importance for investigating the functions of selenium in nature. Previous strategies for predicting selenoprotein genes mainly relied on the identification of a special cis-acting Sec insertion sequence (SECIS) element within mRNAs. However, due to the complexity and variability of SECIS elements, recognition of all selenoprotein genes in bacteria is still a major challenge in the annotation of bacterial genomes. We have developed a deep learning-based algorithm to predict selenoprotein genes in bacterial genomic sequences, which demonstrates superior performance compared to currently available methods. This algorithm can be utilized in either web-based or local (standalone) modes, serving as a promising tool for identifying the complete set of selenoprotein genes in bacteria.

Asgard archaeal selenoproteome reveals a roadmap for the archaea-to-eukaryote transition of selenocysteine incorporation machinery

Selenocysteine (Sec) is encoded by the UGA codon that normally functions as a stop signal and is specifically incorporated into selenoproteins via a unique recoding mechanism. The translational recoding of UGA as Sec is directed by an unusual RNA structure, the SECIS element. Although archaea and eukaryotes adopt similar Sec encoding machinery, the SECIS elements have no similarities to each other with regard to sequence and structure. We analyzed >400 Asgard archaeal genomes to examine the occurrence of both Sec encoding system and selenoproteins in this archaeal superphylum, the closest prokaryotic relatives of eukaryotes. A comprehensive map of Sec utilization trait has been generated, providing the most detailed understanding of the use of this nonstandard amino acid in Asgard archaea so far. By characterizing the selenoproteomes of all organisms, several selenoprotein-rich phyla and species were identified. Most Asgard archaeal selenoprotein genes possess eukaryotic SECIS-like structures with varying degrees of diversity. Moreover, euryarchaeal SECIS elements might originate from Asgard archaeal SECIS elements via lateral gene transfer, indicating a complex and dynamic scenario of the evolution of SECIS element within archaea. Finally, a roadmap for the transition of eukaryotic SECIS elements from archaea was proposed, and selenophosphate synthetase may serve as a potential intermediate for the generation of ancestral eukaryotic SECIS element. Our results offer new insights into a deeper understanding of the evolution of Sec insertion machinery.

MESPEUS: a database of metal coordination groups in proteins

MESPEUS is a freely accessible database which uses carefully selected metal coordination groups found in metalloprotein structures taken from the Protein Data Bank. The database contains geometrical information of metal sites within proteins, including 40 metal types. In order to completely determine the metal coordination, the symmetry-related units of a given protein structure are taken into account and are generated using the appropriate space group symmetry operations. This permits a more complete description of the metal coordination geometry by including all coordinating atoms. The user-friendly web interface allows users to directly search for a metal site of interest using several useful options, including searching for metal elements, metal-donor distances, coordination number, donor residue group, and structural resolution. These searches can be carried out singly or in combination. The details of a metal site and the metal site(s) in the whole structure can be graphically displayed using the interactive web interface. MESPEUS is automatically updated monthly by synchronizing with the PDB database. An investigation for the Mg-ATP interaction is given to demonstrate how MESPEUS can be used to extract information about metal sites by selecting structure and coordination features. MESPEUS is available at http://mespeus.nchu.edu.tw/.

The Mechanism of Metal Homeostasis in Plants: A New View on the Synergistic Regulation Pathway of Membrane Proteins, Lipids and Metal Ions

Heavy metal stress (HMS) is one of the most destructive abiotic stresses which seriously affects the growth and development of plants. Recent studies have shown significant progress in understanding the molecular mechanisms underlying plant tolerance to HMS. In general, three core signals are involved in plants' responses to HMS; these are mitogen-activated protein kinase (MAPK), calcium, and hormonal (abscisic acid) signals. In addition to these signal components, other regulatory factors, such as microRNAs and membrane proteins, also play an important role in regulating HMS responses in plants. Membrane proteins interact with the highly complex and heterogeneous lipids in the plant cell environment. The function of membrane proteins is affected by the interactions between lipids and lipid-membrane proteins. Our review findings also indicate the possibility of membrane protein-lipid-metal ion interactions in regulating metal homeostasis in plant cells. In this review, we investigated the role of membrane proteins with specific substrate recognition in regulating cell metal homeostasis. The understanding of the possible interaction networks and upstream and downstream pathways is developed. In addition, possible interactions between membrane proteins, metal ions, and lipids are discussed to provide new ideas for studying metal homeostasis in plant cells.

Macrophage Selenoproteins Restrict Intracellular Replication of Francisella tularensis and Are Essential for Host Immunity

The essential micronutrient Selenium (Se) is co-translationally incorporated as selenocysteine into proteins. Selenoproteins contain one or more selenocysteines and are vital for optimum immunity. Interestingly, many pathogenic bacteria utilize Se for various biological processes suggesting that Se may play a role in bacterial pathogenesis. A previous study had speculated that Francisella tularensis, a facultative intracellular bacterium and the causative agent of tularemia, sequesters Se by upregulating Se-metabolism genes in type II alveolar epithelial cells. Therefore, we investigated the contribution of host vs. pathogen-associated selenoproteins in bacterial disease using F. tularensis as a model organism. We found that F. tularensis was devoid of any Se utilization traits, neither incorporated elemental Se, nor exhibited Se-dependent growth. However, 100% of Se-deficient mice (0.01 ppm Se), which express low levels of selenoproteins, succumbed to F. tularensis-live vaccine strain pulmonary challenge, whereas 50% of mice on Se-supplemented (0.4 ppm Se) and 25% of mice on Se-adequate (0.1 ppm Se) diet succumbed to infection. Median survival time for Se-deficient mice was 8 days post-infection while Se-supplemented and -adequate mice was 11.5 and >14 days post-infection, respectively. Se-deficient macrophages permitted significantly higher intracellular bacterial replication than Se-supplemented macrophages ex vivo, corroborating in vivo observations. Since Francisella replicates in alveolar macrophages during the acute phase of pneumonic infection, we hypothesized that macrophage-specific host selenoproteins may restrict replication and systemic spread of bacteria. F. tularensis infection led to an increased expression of several macrophage selenoproteins, suggesting their key role in limiting bacterial replication. Upon challenge with F. tularensis, mice lacking selenoproteins in macrophages (TrspM) displayed lower survival and increased bacterial burden in the lung and systemic tissues in comparison to WT littermate controls. Furthermore, macrophages from TrspM mice were unable to restrict bacterial replication ex vivo in comparison to macrophages from littermate controls. We herein describe a novel function of host macrophage-specific selenoproteins in restriction of intracellular bacterial replication. These data suggest that host selenoproteins may be considered as novel targets for modulating immune response to control a bacterial infection.

InterMetalDB: A Database and Browser of Intermolecular Metal Binding Sites in Macromolecules with Structural Information

InterMetalDB is a free-of-charge database and browser of intermolecular metal binding sites that are present on the interfaces of macromolecules forming larger assemblies based on structural information deposited in Protein Data Bank (PDB). It can be found and freely used at https://intermetaldb.biotech.uni.wroc.pl/. InterMetalDB collects the interfacial binding sites with involvement of metal ions and clusters them on the basis of 50% sequence similarity and the nearest metal environment (5 Å radius). The data are available through the web interface where they can be queried, viewed, and downloaded. Complexity of the query depends on the user, because the questions in the query are connected with each other by a logical AND. InterMetalDB offers several useful options for filtering records including searching for structures by particular parameters such as structure resolution, structure description, and date of deposition. Records can be filtered by coordinated metal ion, number of bound amino acid residues, coordination sphere, and other features. InterMetalDB is regularly updated and will continue to be regularly updated with new content in the future. InterMetalDB is a useful tool for all researchers interested in metalloproteins, protein engineering, and metal-driven oligomerization.

Bioinformatics of Selenoproteins

Significance: Bioinformatics has brought important insights into the field of selenium research. The progress made in the development of computational tools in the last two decades, coordinated with growing genome resources, provided new opportunities to study selenoproteins. The present review discusses existing tools for selenoprotein gene finding and other bioinformatic approaches to study the biology of selenium. Recent Advances: The availability of complete selenoproteomes allowed assessing a global distribution of the use of selenocysteine (Sec) across the tree of life, as well as studying the evolution of selenoproteins and their biosynthetic pathway. Beyond gene identification and characterization, human genetic variants in selenoprotein genes were used to examine adaptations to selenium levels in diverse human populations and to estimate selective constraints against gene loss. Critical Issues: The synthesis of selenoproteins is essential for development in mice. In humans, several mutations in selenoprotein genes have been linked to rare congenital disorders. And yet, the mechanism of Sec insertion and the regulation of selenoprotein synthesis in mammalian cells are not completely understood. Future Directions: Omics technologies offer new possibilities to study selenoproteins and mechanisms of Sec incorporation in cells, tissues, and organisms.

Bioinformatics of Metalloproteins and Metalloproteomes

Trace metals are inorganic elements that are required for all organisms in very low quantities. They serve as cofactors and activators of metalloproteins involved in a variety of key cellular processes. While substantial effort has been made in experimental characterization of metalloproteins and their functions, the application of bioinformatics in the research of metalloproteins and metalloproteomes is still limited. In the last few years, computational prediction and comparative genomics of metalloprotein genes have arisen, which provide significant insights into their distribution, function, and evolution in nature. This review aims to offer an overview of recent advances in bioinformatic analysis of metalloproteins, mainly focusing on metalloprotein prediction and the use of different metals across the tree of life. We describe current computational approaches for the identification of metalloprotein genes and metal-binding sites/patterns in proteins, and then introduce a set of related databases. Furthermore, we discuss the latest research progress in comparative genomics of several important metals in both prokaryotes and eukaryotes, which demonstrates divergent and dynamic evolutionary patterns of different metalloprotein families and metalloproteomes. Overall, bioinformatic studies of metalloproteins provide a foundation for systematic understanding of trace metal utilization in all three domains of life.

Abundance and taxonomic affiliation of molybdenum transport and utilization genes in Tengchong hot springs, China

The nitrogen, sulfur and carbon cycles all rely on critical microbial transformations that are carried out by enzymes that require molybdenum (Mo) as a cofactor. Despite Mo importance in these biogeochemical cycles, little information exists about microbial Mo utilization in extreme environments where, due to geochemical conditions, bioavailable Mo may be limited. Using metagenomic data from nine hot springs in Tengchong, Yunnan Province, China, which range in temperature from 42°C to 96°C and pH from 2.3 to 9, the effects of pH, temperature and spring geochemistry on the abundance and taxonomic affiliation of genes related to Mo were studied. Dissolved Mo was only detected at sites with circumneutral pH. However, processes and organisms that require Mo were detected at all sites across all temperature and pH gradients. All sites contained xanthine dehydrogenase, formate dehydrogenase, carbon-monoxide dehydrogenase, nitrate reductase, sulfite oxidase and methionine-sulfoxide reductase despite different community compositions. This suggests that different microbial communities, resulting from different physicochemical conditions, may be performing similar metabolic functions. Furthermore, the abundance and taxonomic diversity of Mo-related annotations increased with higher concentrations of Mo. This study shows that despite geochemical conditions that can limit Mo bioavailability, microbes require Mo for a variety of processes.

Trace Elements and Healthcare: A Bioinformatics Perspective

Biological trace elements are essential for human health. Imbalance in trace element metabolism and homeostasis may play an important role in a variety of diseases and disorders. While the majority of previous researches focused on experimental verification of genes involved in trace element metabolism and those encoding trace element-dependent proteins, bioinformatics study on trace elements is relatively rare and still at the starting stage. This chapter offers an overview of recent progress in bioinformatics analyses of trace element utilization, metabolism, and function, especially comparative genomics of several important metals. The relationship between individual elements and several diseases based on recent large-scale systematic studies such as genome-wide association studies and case-control studies is discussed. Lastly, developments of ionomics and its recent application in human health are also introduced.

Comparative genomics reveals new evolutionary and ecological patterns of selenium utilization in bacteria

Selenium (Se) is an important micronutrient for many organisms, which is required for the biosynthesis of selenocysteine, selenouridine and Se-containing cofactor. Several key genes involved in different Se utilization traits have been characterized; however, systematic studies on the evolution and ecological niches of Se utilization are very limited. Here, we analyzed more than 5200 sequenced organisms to examine the occurrence patterns of all Se traits in bacteria. A global species map of all Se utilization pathways has been generated, which demonstrates the most detailed understanding of Se utilization in bacteria so far. In addition, the selenophosphate synthetase gene, which is used to define the overall Se utilization, was also detected in some organisms that do not have any of the known Se traits, implying the presence of a novel Se form in this domain. Phylogenetic analyses of components of different Se utilization traits revealed new horizontal gene transfer events for each of them. Moreover, by characterizing the selenoproteomes of all organisms, we found a new selenoprotein-rich phylum and additional selenoprotein-rich species. Finally, the relationship between ecological environments and Se utilization was investigated and further verified by metagenomic analysis of environmental samples, which indicates new macroevolutionary trends of each Se utilization trait in bacteria. Our data provide insights into the general features of Se utilization in bacteria and should be useful for a further understanding of the evolutionary dynamics of Se utilization in nature.

Evolution of the Selenoproteome in Helicobacter pylori and Epsilonproteobacteria

By competing for the acquisition of essential nutrients, Helicobacter pylori has the unique ability to persist in the human stomach, also causing nutritional insufficiencies in the host. Although the H. pylori genome apparently encodes selenocysteine synthase (SelA, HP1513), a key pyridoxal phosphate (PLP)-dependent enzyme for the incorporation of selenium into bacterial proteins, nothing is known about the use of this essential element in protein synthesis by this pathogen. We analyzed the evolution of the complete machinery for incorporation of selenium into proteins and the selenoproteome of several H. pylori strains and related Epsilonproteobacteria. Our searches identified the presence of selenoproteins-including the previously unknown DUF466 family-in various Epsilonproteobacteria, but not in H. pylori. We found that a complete system for selenocysteine incorporation was present in the Helicobacteriaceae ancestor and has been recently lost before the split of Helicobacter acinonychis and H. pylori. Our results indicate that H. pylori, at variance with other gastric and enterohepatic Helicobacter, does not use selenocysteine in protein synthesis and does not use selenium for tRNA wobble base modification. However, selA has survived as a functional gene, having lost the domain for the binding of selenocysteine tRNA, but maintaining the ability to bind the PLP cofactor. The evolutionary modifications described for the SelA protein of H. pylori find parallels in other bacterial and archaeal species, suggesting that an alternative enzymatic function is hidden in many proteins annotated as selenocysteinyl-tRNA synthase.

Molybdenum cofactor: a key component of Mycobacterium tuberculosis pathogenesis?

Mycobacterium tuberculosis (Mtb) and other members of the Mtb complex possess an expanded complement of genes for the biosynthesis of molybdenum cofactor (MoCo), a tricyclic pterin molecule that is covalently attached to molybdate. This cofactor allows the redox properties of molybdenum to be harnessed by enzymes in order to catalyze redox reactions in carbon, nitrogen and sulfur metabolism. In this article, we summarize recent advances in elucidating the MoCo biosynthetic pathway in Mtb and highlight the evidence implicating the biosynthesis of this cofactor, as well as the enzymes that depend upon it for activity, in Mtb pathogenesis.

Comparative genomics analysis of the metallomes

Biological trace metals are needed in small quantities, but used by all living organisms. They are employed in key cellular functions in a variety of biological processes, resulting in the various degree of dependence of organisms on metals. Most effort in the field has been placed on experimental studies of metal utilization pathways and metal-dependent proteins. On the other hand, systemic level analyses of metalloproteomes (or metallomes) have been limited for most metals. In this chapter, we focus on the recent advances in comparative genomics, which provides many insights into evolution and function of metal utilization. These studies suggested that iron and zinc are widely used in biology (presumably by all organisms), whereas some other metals such as copper, molybdenum, nickel, and cobalt, show scattered occurrence in various groups of organisms. For these metals, most user proteins are well characterized and their dependence on a specific element is evolutionarily conserved. We also discuss evolutionary dynamics of the dependence of user proteins on different metals. Overall, comparative genomics analysis of metallomes provides a foundation for the systemic level understanding of metal utilization as well as for investigating the general features, functions, and evolutionary dynamics of metal use in the three domains of life.

Selenium-dependent biogenesis of formate dehydrogenase in Campylobacter jejuni is controlled by the fdhTU accessory genes

The food-borne bacterial pathogen Campylobacter jejuni efficiently utilizes organic acids such as lactate and formate for energy production. Formate is rapidly metabolized via the activity of the multisubunit formate dehydrogenase (FDH) enzyme, of which the FdhA subunit is predicted to contain a selenocysteine (SeC) amino acid. In this study we investigated the function of the cj1500 and cj1501 genes of C. jejuni, demonstrate that they are involved in selenium-controlled production of FDH, and propose the names fdhT and fdhU, respectively. Insertional inactivation of fdhT or fdhU in C. jejuni resulted in the absence of FdhA and FdhB protein expression, reduced fdhABC RNA levels, the absence of FDH enzyme activity, and the lack of formate utilization, as assessed by (1)H nuclear magnetic resonance. The fdhABC genes are transcribed from a single promoter located two genes upstream of fdhA, and the decrease in fdhABC RNA levels in the fdhU mutant is mediated at the posttranscriptional level. FDH activity and the ability to utilize formate were restored by genetic complementation with fdhU and by supplementation of the growth media with selenium dioxide. Disruption of SeC synthesis by inactivation of the selA and selB genes also resulted in the absence of FDH activity, which could not be restored by selenium supplementation. Comparative genomic analysis suggests a link between the presence of selA and fdhTU orthologs and the predicted presence of SeC in FdhA. The fdhTU genes encode accessory proteins required for FDH expression and activity in C. jejuni, possibly by contributing to acquisition or utilization of selenium.

Redox biology: computational approaches to the investigation of functional cysteine residues

Cysteine (Cys) residues serve many functions, such as catalysis, stabilization of protein structure through disulfides, metal binding, and regulation of protein function. Cys residues are also subject to numerous post-translational modifications. In recent years, various computational tools aiming at classifying and predicting different functional categories of Cys have been developed, particularly for structural and catalytic Cys. On the other hand, given complexity of the subject, bioinformatics approaches have been less successful for the investigation of regulatory Cys sites. In this review, we introduce different functional categories of Cys residues. For each category, an overview of state-of-the-art bioinformatics methods and tools is provided, along with examples of successful applications and potential limitations associated with each approach. Finally, we discuss Cys-based redox switches, which modify the view of distinct functional categories of Cys in proteins.