An integrative computational model for large-scale identification of metalloproteins in microbial genomes: a focus on iron-sulfur cluster proteins

Proteomic strategies to interrogate the Fe-S proteome

Iron‑sulfur (Fe-S) clusters, inorganic cofactors composed of iron and sulfide, participate in numerous essential redox, non-redox, structural, and regulatory biological processes within the cell. Though structurally and functionally diverse, the list of all proteins in an organism capable of binding one or more Fe-S clusters is referred to as its Fe-S proteome. Importantly, the Fe-S proteome is highly dynamic, with continuous cluster synthesis and delivery by complex Fe-S cluster biogenesis pathways. This cluster delivery is balanced out by processes that can result in loss of Fe-S cluster binding, such as redox state changes, iron availability, and oxygen sensitivity. Despite continued expansion of the Fe-S protein catalogue, it remains a challenge to reliably identify novel Fe-S proteins. As such, high-throughput techniques that can report on native Fe-S cluster binding are required to both identify new Fe-S proteins, as well as characterize the in vivo dynamics of Fe-S cluster binding. Due to the recent rapid growth in mass spectrometry, proteomics, and chemical biology, there has been a host of techniques developed that are applicable to the study of native Fe-S proteins. This review will detail both the current understanding of the Fe-S proteome and Fe-S cluster biology as well as describing state-of-the-art proteomic strategies for the study of Fe-S clusters within the context of a native proteome.

The F plasmid conjutome: the repertoire of E. coli proteins translocated through an F-encoded type IV secretion system

Bacterial conjugation systems pose a major threat to human health through their widespread dissemination of mobile genetic elements (MGEs) carrying cargoes of antibiotic resistance genes. Using the Cre Recombinase Assay for Translocation (CRAfT), we recently reported that the IncFV pED208 conjugation system also translocates at least 16 plasmid-encoded proteins to recipient bacteria. Here, we deployed a high-throughput CRAfT screen to identify the repertoire of chromosomally encoded protein substrates of the pED208 system. We identified 32 substrates encoded by the Escherichia coli W3110 genome with functions associated with (i) DNA/nucleotide metabolism, (ii) stress tolerance/physiology, (iii) transcriptional regulation, or (iv) toxin inhibition. The respective gene deletions did not impact pED208 transfer proficiencies, nor did Group 1 (DNA/nucleotide metabolism) mutations detectably alter the SOS response elicited in new transconjugants upon acquisition of pED208. However, MC4100(pED208) donor cells intrinsically exhibit significantly higher SOS activation than plasmid-free MC4100 cells, and this plasmid carriage-induced stress response is further elevated in donor cells deleted of several Group 1 genes. Among 10 characterized substrates, we gained evidence of C-terminal or internal translocation signals that could function independently or synergistically for optimal protein transfer. Remarkably, nearly all tested proteins were also translocated through the IncN pKM101 and IncP RP4 conjugation systems. This repertoire of E. coli protein substrates, here termed the F plasmid "conjutome," is thus characterized by functions of potential benefit to new transconjugants, diverse TSs, and the capacity for promiscuous transfer through heterologous conjugation systems.

IMPORTANCE

Conjugation systems comprise a major subfamily of the type IV secretion systems (T4SSs) and are the progenitors of a second large T4SS subfamily dedicated to translocation of protein effectors. This study examined the capacity of conjugation machines to function as protein translocators. Using a high-throughput reporter screen, we determined that 32 chromosomally encoded proteins are delivered through an F plasmid conjugation system. The translocated proteins potentially enhance the establishment of the co-transferred F plasmid or mitigate mating-induced stresses. Translocation signals located C-terminally or internally conferred substrate recognition by the F system and, remarkably, many substrates also were translocated through heterologous conjugation systems. Our findings highlight the plasticity of conjugation systems in their capacities to co-translocate DNA and many protein substrates.

Iron-sulfur protein odyssey: exploring their cluster functional versatility and challenging identification

Iron-sulfur (Fe-S) clusters are an essential and ubiquitous class of protein-bound prosthetic centers that are involved in a broad range of biological processes (e.g. respiration, photosynthesis, DNA replication and repair and gene regulation) performing a wide range of functions including electron transfer, enzyme catalysis, and sensing. In a general manner, Fe-S clusters can gain or lose electrons through redox reactions, and are highly sensitive to oxidation, notably by small molecules such as oxygen and nitric oxide. The [2Fe-2S] and [4Fe-4S] clusters, the most common Fe-S cofactors, are typically coordinated by four amino acid side chains from the protein, usually cysteine thiolates, but other residues (e.g. histidine, aspartic acid) can also be found. While diversity in cluster coordination ensures the functional variety of the Fe-S clusters, the lack of conserved motifs makes new Fe-S protein identification challenging especially when the Fe-S cluster is also shared between two proteins as observed in several dimeric transcriptional regulators and in the mitoribosome. Thanks to the recent development of in cellulo, in vitro, and in silico approaches, new Fe-S proteins are still regularly identified, highlighting the functional diversity of this class of proteins. In this review, we will present three main functions of the Fe-S clusters and explain the difficulties encountered to identify Fe-S proteins and methods that have been employed to overcome these issues.

Monitoring Fe-S cluster occupancy across the E. coli proteome using chemoproteomics

Iron-sulfur (Fe-S) clusters are ubiquitous metallocofactors involved in redox chemistry, radical generation and gene regulation. Common methods to monitor Fe-S clusters include spectroscopic analysis of purified proteins and autoradiographic visualization of radiolabeled iron distribution in proteomes. Here, we report a chemoproteomic strategy that monitors changes in the reactivity of Fe-S cysteine ligands to inform on Fe-S cluster occupancy. We highlight the utility of this platform in Escherichia coli by (1) demonstrating global disruptions in Fe-S incorporation in cells cultured under iron-depleted conditions, (2) determining Fe-S client proteins reliant on five scaffold, carrier and chaperone proteins within the Isc Fe-S biogenesis pathway and (3) identifying two previously unannotated Fe-S proteins, TrhP and DppD. In summary, the chemoproteomic strategy described herein is a powerful tool that reports on Fe-S cluster incorporation directly within a native proteome, enabling the interrogation of Fe-S biogenesis pathways and the identification of previously uncharacterized Fe-S proteins.

Revealing a New Family of D-2-Hydroxyglutarate Dehydrogenases in Escherichia coli and Pantoea ananatis Encoded by ydiJ

In E. coli and P. ananatis, L-serine biosynthesis is initiated by the action of D-3-phosphoglycerate dehydrogenase (SerA), which converts D-3-phosphoglycerate into 3-phosphohydroxypyruvate. SerA can concomitantly catalyze the production of D-2-hydroxyglutarate (D-2-HGA) from 2-ketoglutarate by oxidizing NADH to NAD⁺. Several bacterial D-2-hydroxyglutarate dehydrogenases (D2HGDHs) have recently been identified, which convert D-2-HGA back to 2-ketoglutarate. However, knowledge about the enzymes that can metabolize D-2-HGA is lacking in bacteria belonging to the Enterobacteriaceae family. We found that ydiJ encodes novel D2HGDHs in P. ananatis and E. coli, which were assigned as D2HGDHPa and D2HGDHEc, respectively. Inactivation of ydiJ in P. ananatis and E. coli led to the significant accumulation of D-2-HGA. Recombinant D2HGDHEc and D2HGDHPa were purified to homogeneity and characterized. D2HGDHEc and D2HGDHPa are homotetrameric with a subunit molecular mass of 110 kDa. The pH optimum was 7.5 for D2HGDHPa and 8.0 for D2HGDHEc. The Km for D-2-HGA was 208 μM for D2HGDHPa and 83 μM for D2HGDHEc. The enzymes have strict substrate specificity towards D-2-HGA and displayed maximal activity at 45 °C. Their activity was completely inhibited by 0.5 mM Mn²⁺, Ni²⁺ or Co²⁺. The discovery of a novel family of D2HGDHs may provide fundamental information for the metabolic engineering of microbial chassis with desired properties.

The Fe-S proteome of Escherichia coli: prediction, function and fate

Iron-sulfur (Fe-S) clusters are inorganic ubiquitous and ancient cofactors. Fe-S bound proteins contribute to most cellular processes, including DNA replication and integrity, genetic expression and regulation, metabolism, biosynthesis and most bioenergetics systems. Also, Fe-S proteins hold a great biotechnological potential in metabolite and chemical production, including antibiotics. From classic biophysics and spectroscopy methodologies to recent development in bioinformatics, including structural modeling and chemoproteomics, our capacity to predict and identify Fe-S proteins has spectacularly increased over the recent years. Here, these developments are presented and collectively used to update the composition of Escherichia coli Fe-S proteome, for which we predict 181 occurrences, i.e. 40 more candidates than in our last catalog (Py and Barras, 2010), and equivalent to 4% of its total proteome. Besides, Fe-S clusters can be targeted by redox active compounds or reactive oxygen and nitrosative species, and even be destabilized by contaminant metals. Accordingly, we discuss how cells handle damaged Fe-S proteins, i.e. degradation, recycling or repair.

Machine Learning Approaches for Metalloproteins

Metalloproteins are a family of proteins characterized by metal ion binding, whereby the presence of these ions confers key catalytic and ligand-binding properties. Due to their ubiquity among biological systems, researchers have made immense efforts to predict the structural and functional roles of metalloproteins. Ultimately, having a comprehensive understanding of metalloproteins will lead to tangible applications, such as designing potent inhibitors in drug discovery. Recently, there has been an acceleration in the number of studies applying machine learning to predict metalloprotein properties, primarily driven by the advent of more sophisticated machine learning algorithms. This review covers how machine learning tools have consolidated and expanded our comprehension of various aspects of metalloproteins (structure, function, stability, ligand-binding interactions, and inhibitors). Future avenues of exploration are also discussed.

Identification of Iron-Sulfur (Fe-S) Cluster and Zinc (Zn) Binding Sites Within Proteomes Predicted by DeepMind's AlphaFold2 Program Dramatically Expands the Metalloproteome

DeepMind's AlphaFold2 software has ushered in a revolution in high quality, 3D protein structure prediction. In very recent work by the DeepMind team, structure predictions have been made for entire proteomes of twenty-one organisms, with >360,000 structures made available for download. Here we show that thousands of novel binding sites for iron-sulfur (Fe-S) clusters and zinc (Zn) ions can be identified within these predicted structures by exhaustive enumeration of all potential ligand-binding orientations. We demonstrate that AlphaFold2 routinely makes highly specific predictions of ligand binding sites: for example, binding sites that are comprised exclusively of four cysteine sidechains fall into three clusters, representing binding sites for 4Fe-4S clusters, 2Fe-2S clusters, or individual Zn ions. We show further: (a) that the majority of known Fe-S cluster and Zn binding sites documented in UniProt are recovered by the AlphaFold2 structures, (b) that there are occasional disputes between AlphaFold2 and UniProt with AlphaFold2 predicting highly plausible alternative binding sites, (c) that the Fe-S cluster binding sites that we identify in E. coli agree well with previous bioinformatics predictions, (d) that cysteines predicted here to be part of ligand binding sites show little overlap with those shown via chemoproteomics techniques to be highly reactive, and (e) that AlphaFold2 occasionally appears to build erroneous disulfide bonds between cysteines that should instead coordinate a ligand. These results suggest that AlphaFold2 could be an important tool for the functional annotation of proteomes, and the methodology presented here is likely to be useful for predicting other ligand-binding sites.

Prediction of the Iron–Sulfur Binding Sites in Proteins Using the Highly Accurate Three-Dimensional Models Calculated by AlphaFold and RoseTTAFold

AlphaFold and RoseTTAFold are deep learning-based approaches that predict the structure of proteins from their amino acid sequences. Remarkable success has recently been achieved in the prediction accuracy of not only the fold of the target protein but also the position of its amino acid side chains. In this article, I question the accuracy of these methods to predict iron–sulfur binding sites. I analyze three-dimensional models calculated by AlphaFold and RoseTTAFold of Fe–S–dependent enzymes, for which no structure of a homologous protein has been solved experimentally. In all cases, the amino acids that presumably coordinate the cluster were gathered together and facing each other, which led to a quite accurate model of the Fe–S cluster binding site. Yet, cysteine candidates were often involved in intramolecular disulfide bonds, and the number and identity of the protein amino acids that should ligate the cluster were not always clear. The experimental structure determination of the protein with its Fe–S cluster and in complex with substrate/inhibitor/product is still needed to unambiguously visualize the coordination state of the cluster and understand the conformational changes occurring during catalysis.

Mechanism of Iron–Sulfur Cluster Assembly: In the Intimacy of Iron and Sulfur Encounter

Iron–sulfur (Fe–S) clusters are protein cofactors of a multitude of enzymes performing essential biological functions. Specialized multi-protein machineries present in all types of organisms support their biosynthesis. These machineries encompass a scaffold protein on which Fe–S clusters are assembled and a cysteine desulfurase that provides sulfur in the form of a persulfide. The sulfide ions are produced by reductive cleavage of the persulfide, which involves specific reductase systems. Several other components are required for Fe–S biosynthesis, including frataxin, a key protein of controversial function and accessory components for insertion of Fe–S clusters in client proteins. Fe–S cluster biosynthesis is thought to rely on concerted and carefully orchestrated processes. However, the elucidation of the mechanisms of their assembly has remained a challenging task due to the biochemical versatility of iron and sulfur and the relative instability of Fe–S clusters. Nonetheless, significant progresses have been achieved in the past years, using biochemical, spectroscopic and structural approaches with reconstituted system in vitro. In this paper, we review the most recent advances on the mechanism of assembly for the founding member of the Fe–S cluster family, the [2Fe2S] cluster that is the building block of all other Fe–S clusters. The aim is to provide a survey of the mechanisms of iron and sulfur insertion in the scaffold proteins by examining how these processes are coordinated, how sulfide is produced and how the dinuclear [2Fe2S] cluster is formed, keeping in mind the question of the physiological relevance of the reconstituted systems. We also cover the latest outcomes on the functional role of the controversial frataxin protein in Fe–S cluster biosynthesis.

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.

Ubiquinone Biosynthesis over the Entire O2 Range: Characterization of a Conserved O2-Independent Pathway

In order to colonize environments with large O₂ gradients or fluctuating O₂ levels, bacteria have developed metabolic responses that remain incompletely understood. Such adaptations have been recently linked to antibiotic resistance, virulence, and the capacity to develop in complex ecosystems like the microbiota. Here, we identify a novel pathway for the biosynthesis of ubiquinone, a molecule with a key role in cellular bioenergetics. We link three uncharacterized genes of Escherichia coli to this pathway and show that the pathway functions independently from O₂. In contrast, the long-described pathway for ubiquinone biosynthesis requires O₂ as a substrate. In fact, we find that many proteobacteria are equipped with the O₂-dependent and O₂-independent pathways, supporting that they are able to synthesize ubiquinone over the entire O₂ range. Overall, we propose that the novel O₂-independent pathway is part of the metabolic plasticity developed by proteobacteria to face various environmental O₂ levels.

Most bacteria can generate ATP by respiratory metabolism, in which electrons are shuttled from reduced substrates to terminal electron acceptors, via quinone molecules like ubiquinone. Dioxygen (O₂) is the terminal electron acceptor of aerobic respiration and serves as a co-substrate in the biosynthesis of ubiquinone. Here, we characterize a novel, O₂-independent pathway for the biosynthesis of ubiquinone. This pathway relies on three proteins, UbiT (YhbT), UbiU (YhbU), and UbiV (YhbV). UbiT contains an SCP2 lipid-binding domain and is likely an accessory factor of the biosynthetic pathway, while UbiU and UbiV (UbiU-UbiV) are involved in hydroxylation reactions and represent a novel class of O₂-independent hydroxylases. We demonstrate that UbiU-UbiV form a heterodimer, wherein each protein binds a 4Fe-4S cluster via conserved cysteines that are essential for activity. The UbiT, -U, and -V proteins are found in alpha-, beta-, and gammaproteobacterial clades, including several human pathogens, supporting the widespread distribution of a previously unrecognized capacity to synthesize ubiquinone in the absence of O₂. Together, the O₂-dependent and O₂-independent ubiquinone biosynthesis pathways contribute to optimizing bacterial metabolism over the entire O₂ range.

Transition Metal Homeostasis in Streptococcus pyogenes and Streptococcus pneumoniae

Trace metals such as Fe, Mn, Zn and Cu are essential for various biological functions including proper innate immune function. The host immune system has complicated and coordinated mechanisms in place to either starve and/or overload invading pathogens with various metals to combat the infection. Here, we discuss the roles of Fe, Mn and Zn in terms of nutritional immunity, and also the roles of Cu and Zn in metal overload in relation to the physiology and pathogenesis of two human streptococcal species, Streptococcus pneumoniae and Streptococcus pyogenes. S. pneumoniae is a major human pathogen that is carried asymptomatically in the nasopharynx by up to 70% of the population; however, transition to internal sites can cause a range of diseases such as pneumonia, otitis media, meningitis and bacteraemia. S. pyogenes is a human pathogen responsible for diseases ranging from pharyngitis and impetigo, to severe invasive infections. Both species have overlapping capacity with respect to metal acquisition, export and regulation and how metal homeostasis relates to their virulence and ability to invade and survive within the host. It is becoming more apparent that metals have an important role to play in the control of infection, and with further investigations, it could lead to the potential use of metals in novel antimicrobial therapies.

Biological and Chemical Adaptation to Endogenous Hydrogen Peroxide Production in Streptococcus pneumoniae D39

Adaptation to endogenous oxidative stress is an integral aspect of Streptococcus pneumoniae colonization and virulence. In this work, we identify key transcriptomic and proteomic features of the pneumococcal endogenous oxidative stress response. The thiol peroxidase TpxD plays a critical role in adaptation to endogenous H₂O₂ and serves to limit protein sulfenylation of glycolytic, capsule, and nucleotide biosynthesis enzymes in S. pneumoniae.

The catalase-negative, facultative anaerobe Streptococcus pneumoniae D39 is naturally resistant to hydrogen peroxide (H₂O₂) produced endogenously by pyruvate oxidase (SpxB). Here, we investigate the adaptive response to endogenously produced H₂O₂. We show that lactate oxidase, which converts lactate to pyruvate, positively impacts pyruvate flux through SpxB and that ΔlctO mutants produce significantly lower H₂O₂. In addition, both the SpxB pathway and a candidate pyruvate dehydrogenase complex (PDHC) pathway contribute to acetyl coenzyme A (acetyl-CoA) production during aerobic growth, and the pyruvate format lyase (PFL) pathway is the major acetyl-CoA pathway during anaerobic growth. Microarray analysis of the D39 strain cultured under aerobic versus strict anaerobic conditions shows upregulation of spxB, a gene encoding a rhodanese-like protein (locus tag spd0091), tpxD, sodA, piuB, piuD, and an Fe-S protein biogenesis operon under H₂O₂-producing conditions. Proteome profiling of H₂O₂-induced sulfenylation reveals that sulfenylation levels correlate with cellular H₂O₂ production, with endogenous sulfenylation of ≈50 proteins. Deletion of tpxD increases cellular sulfenylation 5-fold and has an inhibitory effect on ATP generation. Two major targets of protein sulfenylation are glyceraldehyde-3-phosphate dehydrogenase (GapA) and SpxB itself, but targets also include pyruvate kinase, LctO, AdhE, and acetate kinase (AckA). Sulfenylation of GapA is inhibitory, while the effect on SpxB activity is negligible. Strikingly, four enzymes of capsular polysaccharide biosynthesis are sulfenylated, as are enzymes associated with nucleotide biosynthesis via ribulose-5-phosphate. We propose that LctO/SpxB-generated H₂O₂ functions as a signaling molecule to downregulate capsule production and drive altered flux through sugar utilization pathways.

IMPORTANCE Adaptation to endogenous oxidative stress is an integral aspect of Streptococcus pneumoniae colonization and virulence. In this work, we identify key transcriptomic and proteomic features of the pneumococcal endogenous oxidative stress response. The thiol peroxidase TpxD plays a critical role in adaptation to endogenous H₂O₂ and serves to limit protein sulfenylation of glycolytic, capsule, and nucleotide biosynthesis enzymes in S. pneumoniae.

Contemplating a role for titanium in organisms

Titanium is the ninth most abundant element in the Earth's crust and some organisms sequester it avidly, though no essential biological role has yet been recognized. This Minireview addresses how the properties of titanium, especially in an oxic aqueous environment, might make a biological role difficult to recognize. It further considers how new -omic technologies might overcome the limitations of the past and help to reveal a specific role for this metal. While studies with well established model organisms have their rightful place, organisms that are known avid binders or sequesterers of titanium should be promising places to investigate a biological role.

MetalPredator: a web server to predict iron-sulfur cluster binding proteomes

MOTIVATION

The prediction of the iron-sulfur proteome is highly desirable for biomedical and biological research but a freely available tool to predict iron-sulfur proteins has not been developed yet.

RESULTS

We developed a web server to predict iron-sulfur proteins from protein sequence(s). This tool, called MetalPredator, is able to process complete proteomes rapidly with high recall and precision.

AVAILABILITY AND IMPLEMENTATION

The web server is freely available at: http://metalweb.cerm.unifi.it/tools/metalpredator/

CONTACT

andreini@cerm.unifi.it

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Microbial Metalloproteomics

Metalloproteomics is a rapidly developing field of science that involves the comprehensive analysis of all metal-containing or metal-binding proteins in a biological sample. The purpose of this review is to offer a comprehensive overview of the research involving approaches that can be categorized as inductively coupled plasma (ICP)-MS based methods, X-ray absorption/fluorescence, radionuclide based methods and bioinformatics. Important discoveries in microbial proteomics will be reviewed, as well as the outlook to new emerging approaches and research areas.

Chemical and Biological Reduction of the Radical SAM Enzyme 7-Carboxy-7-deazaguanine [corrected] Synthase

The radical S-adenosyl-L-methionine (SAM) superfamily is a large and growing group of enzymes that conduct complex radical-mediated transformations. A one-electron reduction of SAM via the +1 state of the cubane [4Fe-4S] cluster generates a 5'-deoxyadenosyl radical, which initiates turnover. The [4Fe-4S] cluster must be reduced from its resting +2 state to the catalytically active +1 oxidation state by an electron. In practice, dithionite or the Escherichia coli flavodoxin (EcFldA)/ferredoxin (flavodoxin):NADP(+) oxidoreductase (Fpr)/NADPH system is used. Herein, we present a systematic investigation of the reductive activation of the radical SAM enzyme CDG synthase (BsQueE) from Bacillus subtilis comparing biological and chemical reductants. These data show that either of the flavodoxin homologues encoded by the B. subtilis genome, BsYkuN or BsYkuP, as well as a series of small molecule redox mediators, supports BsQueE activity. With dithionite as a reductant, the activity of BsQueE is ~75-fold greater in the presence of BsYkuN and BsYkuP compared to that in the presence of dithionite alone. By contrast, EcFldA supports turnover to ~10-fold greater levels than dithionite alone under the same conditions. Comparing the ratio of the rate of turnover to the apparent binding constant for the flavodoxin homologues reveals 10- and 240-fold preferences for BsYkuN over BsYkuP and EcFldA, respectively. The differential activation of the enzyme cannot be explained by the abortive cleavage of SAM. We conclude from these observations that the differential activation of BsQueE by Fld homologues may reside in the details of the interaction between the flavodoxin and the radical SAM enzyme.