Next Article in Journal
Plasmodium falciparum Nicotinamidase as A Novel Antimalarial Target
Next Article in Special Issue
Cytochrome P450 Surface Domains Prevent the β-Carotene Monohydroxylase CYP97H1 of Euglena gracilis from Acting as a Dihydroxylase
Previous Article in Journal
Interaction of the Emerging Mycotoxins Beauvericin, Cyclopiazonic Acid, and Sterigmatocystin with Human Serum Albumin
Previous Article in Special Issue
Midazolam as a Probe for Heterotropic Drug-Drug Interactions Mediated by CYP3A4
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biomolecules
Volume 12
Issue 8
10.3390/biom12081107
biomolecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of a Virally Encoded Flavodoxin That Can Drive Bacterial Cytochrome P450 Monooxygenase Activity

by
David C. Lamb
1,*,
Jared V. Goldstone
2,
Bin Zhao
3,
Li Lei
4,
Jonathan G. L. Mullins
1,
Michael J. Allen
5,
Steven L. Kelly
1 and
John J. Stegeman
2
1
Faculty of Medicine, Health and Life Sciences, Swansea University, Swansea SA2 8PP, UK
2
Biology Department, Woods Hole Oceanographic Institution, Woods Hole, MA 02543-1050, USA
3
Cleveland Clinic Lerner Research Institute, 9500 Euclid Avenue, NB21, Cleveland, OH 44195, USA
4
Department of Biochemistry, Vanderbilt University Medical School, Vanderbilt University, Nashville, TN 37232-0146, USA
5
Department of Biosciences, College of Life and Environmental Sciences, University of Exeter, Stocker Road, Exeter EX4 4QD, UK
*
Author to whom correspondence should be addressed.
Biomolecules 2022, 12(8), 1107; https://doi.org/10.3390/biom12081107
Submission received: 31 May 2022 / Revised: 28 July 2022 / Accepted: 4 August 2022 / Published: 11 August 2022
(This article belongs to the Special Issue New Insights into Cytochrome P450s)
Browse Figures
Review Reports Versions Notes

Abstract

:
Flavodoxins are small electron transport proteins that are involved in a myriad of photosynthetic and non-photosynthetic metabolic pathways in Bacteria (including cyanobacteria), Archaea and some algae. The sequenced genome of 0305φ8-36, a large bacteriophage that infects the soil bacterium Bacillus thuringiensis, was predicted to encode a putative flavodoxin redox protein. Here we confirm that 0305φ8-36 phage encodes a FMN-containing flavodoxin polypeptide and we report the expression, purification and enzymatic characterization of the recombinant protein. Purified 0305φ8-36 flavodoxin has near-identical spectral properties to control, purified Escherichia coli flavodoxin. Using in vitro assays we show that 0305φ8-36 flavodoxin can be reconstituted with E. coli flavodoxin reductase and support regio- and stereospecific cytochrome P450 CYP170A1 allyl-oxidation of epi-isozizaene to the sesquiterpene antibiotic product albaflavenone, found in the soil bacterium Streptomyces coelicolor. In vivo, 0305φ8-36 flavodoxin is predicted to mediate the 2-electron reduction of the β subunit of phage-encoded ribonucleotide reductase to catalyse the conversion of ribonucleotides to deoxyribonucleotides during viral replication. Our results demonstrate that this phage flavodoxin has the potential to manipulate and drive bacterial P450 cellular metabolism, which may affect both the host biological fitness and the communal microbiome. Such a scenario may also be applicable in other viral-host symbiotic/parasitic relationships.

1. Importance

Herein, we demonstrate that virally encoded flavodoxin has the potential to manipulate and drive host bacterial P450 cellular metabolism affecting both host biological fitness and the communal microbiome. Viral phenotypic impact was thought limited to merely affecting host mortality. Recent work has cast this aspersion aside, as we have come to realize that close metabolic integration appears to be required for viruses to subdue and coerce their unwitting hosts. Our work herein clearly demonstrates the potential for viruses to not just manipulate ‘virus centric’ replication activities, but also broader host metabolic activities and cellular function(s) as well.

2. Introduction

Flavodoxins (Fld) are small water-soluble electron transfer proteins that contain non-covalently bound flavin mononucleotide (FMN) as the active redox component [1]. Flds were first discovered in the 1960s in both cyanobacteria and Clostridia growing in low iron conditions, where they serve as electron carriers, replacing iron-containing ferredoxin (Fd) in metabolic reactions leading to NADP+ and N2 reduction [2,3]. Since that time, Flds have been isolated from numerous Bacteria, Archaea, cyanobacteria and some eukaryotic algae [1]. Flds are located in the bacterial cytosol and algal chloroplasts and function in electron-consuming metabolic pathways and act as electronic switches between cellular sources of reducing power (e.g., light-driven reactions, pyridine nucleotides, sugars). Biochemically, all characterised flavodoxin proteins to date are highly acidic proteins consisting of ~140–180 amino acid residues (MWs~15–22 kDa) and structurally consist of a five-stranded parallel beta sheet enclosed by alpha helices on either side [1]. Catalytically, the Fld FMN cofactor is responsible for electron transfer, cycling between a neutral semiquinone (sq) and hydroquinone (hq) forms [4,5].
In bacteria the Fld gene is inducible, functioning as an adaptive resource under environmental or nutritional hardships (e.g., iron limitation). Furthermore, Fld has been shown to be an essential gene in some bacteria such as Escherichia coli [6] and Helicobacter pylori [7]. Although not present in higher eukaryotes such as plants and animals, a descendent of the Fld gene appears fused in some eukaryotic genes that encode multidomain redox proteins, e.g., cytochrome P450 reductase (POR) [8]. Flds are divided into short-chain and long-chain types differing by a 20-residue loop of unknown function. Based on phylogenetic analysis it has been proposed that the long-chain Flds may have preceded the shorter ones in evolution [9]. Long chain Flds are found exclusively in cyanobacteria and algae [10,11].
Genes thought to encode flavodoxin have been tentatively annotated in some viral (bacteriophage) genomes. In such phages, the flavodoxin gene is most often associated with genes encoding α and β-ribonucleotide reductase (RNR), and the gene products are thought to form a reversible protein complex that is responsible for converting ribonucleotides to deoxyribonucleotides, a critical reaction for viral DNA synthesis [12,13]. In this reaction, flavodoxin is proposed to mediate the two-electron reduction of the di-ferric cluster of RNR via two successive one-electron transfers. However, this has never been shown experimentally.
Previously, we have reported that genes for multiple and unique monooxygenase cytochromes P450 enzymes are encoded in many genomes of giant viruses (Megavirales) including Mimiviridae and Pandoraviridae [14]. Therein, we also reported that P450 genes are encoded in other viral genomes including a herpesvirus (Ranid herpesvirus 3) and a Mycobacterium phage (Mycobacterium phage Adler). However, genes encoding known P450 redox partners (cytochrome P450 reductase or POR, ferredoxin and ferredoxin reductase, flavodoxin and flavodoxin reductase) were not found in any viral/phage genome thus far described, implying that host redox partners may drive viral P450 activities [14]. Herein, we now show that a bacteriophage 0305φ8-36 [15,16] encodes FMN-bound flavodoxin that can drive bacterial P450 activity, namely CYP170A1 oxidation, to produce the sequiterpene albaflavenone in the soil bacterium Streptomyces coelicolor. This result strongly suggests that phage flavodoxin can influence the biochemical activity of bacterial host enzymes, including P450 enzymes. The environmental relevance of bacteriophages for microbial systems has been widely explored in aquatic environments, but our current understanding of the role of phage(s) in soil ecosystems remains limited. Our results emphasize the importance of accounting for the effect of phage-encoded enzymes on interacting with and manipulating host metabolic pathways and therefore community function, with a novel and exciting example here from the soil microbiome.

3. Materials and Methods

3.1. Bioinformatic and Phylogenetic Analysis

Fld sequences were retrieved from the NCBI sequence databases and BLAST searches employed to identify candidate phage Fld genes. Inferred protein sequences were initially aligned using ClustalOmega [17] and scored with T-Coffee [18]. Alignments were then visualized in BioEdit [19]. Final sequence alignment was constructed as described previously [20]. Virus names are those in the original reports or Genbank. Note that the sequence identities in the text are pairwise amino acid identities based on the length of the shortest protein in the alignment. Identity values used in the tables are based on multiple sequence alignments and are usually not more than 1% different from values stated in the text. Phylogenetic analyses were conducted on a mixture sequences obtained from BLAST searches of the NCBI databases, host sequences, and reference flavodoxin sequences. Sequences were aligned to the PFAM flavodoxin hidden Markov model (FMN_red, PF03358.18) using Clustal-omega v1.2.2 [17]. Maximum likelihood analyses were carried out using RAxML-NG (v1.1.0) [21]. Twenty independent ML searches were carried out followed by 1000 bootstrap replicates. Bootstrap support was assigned using the transfer bootstrap expectation method [22] and visualized using Figtree.

3.2. Cloning, Heterologous Expression and Purification of 0305φ8-36 Flavodoxin

The 0305φ8-36 flavodoxin gene was generated by DNA2.0 (Eurofins) and gene integrity confirmed by DNA sequencing. A unique NdeI restriction site was generated prior to the ATG start codon and a HindIII site following the TGA stop codon to facilitate cloning into the E. coli expression plasmid, pET17b. Additionally, a six residue polyhistidine tag was engineered into the C-terminus of the protein to allow purification of the expressed protein by Ni2+-NTA affinity chromatography. 0305φ8-36 Fld was expressed and purified using similar conditions described for P450 haemoproteins in our laboratory [23,24]. Briefly, transformed BL21(DE3)pLys cells were cultured overnight in Luria Bertani broth containing 100 μg/mL ampicillin. After inoculation (1:100) in 3 litres of Terrific Broth containing 100 μg/mL ampicillin, growth was carried out at 37 °C and 250 rpm for 6 h. Following the addition of riboflavin (1 mM final concentration) for flavin synthesis, flavodoxin expression was induced by the addition of isopropyl-β-D-thiogalactopyranoside (1 mM final concentration). Cell growth was continued for an additional 24 h at 27 °C and 190 rpm. Cells were harvested by centrifugation, resuspended in 250 mM sucrose, 50 mM Tris-HCl, pH 7.4 and disrupted by sonication. The insoluble cell debris was removed by centrifugation at 12,000× g and the supernatant applied to a 5 mL nickel affinity chromatography column (Qiagen) pre-equilibrated with 50 mM Tris-HCl, 100 mM NaCl, 20 mM imidazole, pH 7.4. Unbound protein was removed with further washes, and the His-tagged protein eluted in a single step with 50 mM Tris-HCl, 100 mM NaCl, 500 mM imidazole, pH 7.4. The eluted protein was concentrated to 2 mL using a 3 kDa molecular weight cut off Amicon centrifugal concentrator (Amicon, Merck) and glycerol (20% (v/v)) was added as a cryoprotectant. The purified protein was frozen at −80 °C until used.

3.3. Reconstitution of Streptomyces coelicolor CYP170A1 Activity with 0305φ8-36 Flavodoxin

S. coelicolor CYP170A1 (1 nmol), 0305φ8-36 flavodoxin or E. coli flavodoxin (10 nmol), and E. coli flavodoxin reductase (2 nmol) were reconstituted in 400 μL of 50 mM Tris-HCl buffer (pH 8.2) containing 20% (v/v) glycerol and epiisozizaene (20 nmol) plus 1% (v/v) Me2SO. Following incubation of this mixture for 5 min on ice, the reconstituted enzyme solution was placed in a shaking bath at 35 °C. The reaction was started by the addition of NADPH to a final concentration of 1 mM and 60 μL of an NADPH-regenerating system (glucose-6-phosphate (60 mM), β-nicotinamide adenine dinucleotide phosphate (3 mM), and glucose-6-phosphate dehydrogenase (0.2 unit)). The reaction was carried out for 1.5 h in a 10-mL test tube, at which time it was quenched with 4 μL of concentrated HCl and extracted three times with 400 μL of pentane:methylene chloride (4:1). The extracts were concentrated under a stream of N2, and 2 μL of extract was analyzed by gas chromatography/mass spectrometry (GC/MS). The mixtures were extracted and analyzed by GC/MS. Turnover rates for the epi-isozizaene oxidation was determined from triplicate experiments.

3.4. General Methods

The flavin content of Flds were determined by UV-visible and fluorescence spectrometry [25]. S. coelicolor CYP170A1 was expressed and purified as reported earlier [26]. To improve protein expression, CYP170A1 was coexpressed with molecular chaperones GroES/GroEL. Reduced carbon monoxide (CO) difference spectra for quantification of CYP170A1 content were measured and calculated according to the method described by Omura and Sato [27]. Additionally, recombinant E.coli flavodoxin and flavodoxin reductase were expressed and purified as described previously [28]. Protein quantification was performed by using the bicinchoninic acid assay. Unless otherwise stated, all chemicals were supplied by Sigma Chemical Company (Poole, Dorset, United Kingdom). UV-visible absorption spectra of purified flavodoxins and cytochrome P450 were recorded using a Shimadzu UV-2401 scanning spectrophotometer as described previously [26,29].

4. Results

4.1. Bioinformatic and Phylogenetic Analysis of 0305φ8-36 ORF223

Our earlier work has shown that some giant viruses and a Mycobacterium phage encode cytochrome P450 monooxygenase enzymes, however the associated genes for P450 electron transfer proteins have remained elusive [14]. Hence, we have looked in more detail for the presence of P450 redox partner genes in the virosphere. A BLASTP bioinformatic search of all published viral genomes, using flavodoxin amino acid sequences from a variety of bacterial sources, revealed the presence of a putative flavodoxin gene (ORF223) present in the genome of the bacteriophage 0305φ8-36. 0305φ8-36 has a genome size of 218 kb and contains 247 putative ORFs [15,16]. 0305φ8-36 infects the bacterium Bacillus thuringiensis, a soil bacterium commonly used as a biological pesticide [30]. The complete genome sequence of 0305φ8-36 was published in 2007 (Genbank accession no. EF583821). Examination of the genomic context of 0305φ8-36 ORF223 reveals that this putative flavodoxin gene is directly downstream of 2 ORFs predicted to encode α and β ribonucleotide reductase (ORFs 224 and 225 respectively) and that all three genes are translated from the same mRNA transcript (Figure 1), suggesting gene products that work in a complex.
The 0305φ8-36 ORF223 is predicted to encode a 147 amino acid polypeptide with a predicted molecular weight 16.5 kDa. The isoelectric point of 0305φ8-36 ORF223 was calculated to be approximately 4.9, resulting in an acidic protein similar to other flavodoxin proteins [1]. Amino acid sequence alignment revealed that key bacterial flavodoxin residues essential to binding and stabilising of the FMN cofactor, as determined by protein crystallography analysis of E. coli FldA, are also conserved in the 0305φ8-36 phage flavodoxin protein. These include Thr14 and Thr17 (Thr12 and Thr15 in E. coli flavodoxin), whose side chains form hydrogen bonds to stabilise the FMN phosphate group. Also conserved were phage Trp61 and Tyr97 (Trp57 and Tyr94 in E. coli flavodoxin), which form tight coplanar interactions with the isoalloxazine ring of FMN (Figure 2).
BLASTP analysis of 0305φ8-36 ORF223 against the public databases resulted in the best matches (46–37% sequence identity) to predicted flavodoxin proteins from soil dwelling bacteria including Bacillus massiliosenegalensis (46.5% identity; E-value 2 × 10−35); Bacillus alkalicellulosilyticus (43.2% identity; E-value 5 × 10−33); Rhodococcus qingshengii (43.7% identity; E-value 2 × 10−32) and Mycobacteroides abscessus subsp. Abscessus (42.3% identity; E-value 5 × 10−29). BLASTP analysis of 0305φ8-36 ORF223 against B. thuringiensis genomes revealed matches with a number of predicted flavodoxin proteins but with lower sequence identity (34–31% sequence identity; 95–65% sequence coverage e.g., Bacillus thuringiensis MC28, 33% identity, E-value 7 × 10−16). The closest homologues of most of the 0305φ8-36 virion protein-encoding genes and a few replicative genes were found to reside in a single segment of the chromosome of B. thuringiensis serovar Israelensis. Sequence identity between B. thuringiensis serovar Israelensis flavodoxin and 0305φ8-36 flavodoxin is only 31%. Thus, it is difficult to say with any great certainty whether the 0305φ8-36 phage captured ORF223 by horizontal gene transfer from the host B. thurigiensis or from another soil bacterial species.
We also undertook a BLASTP analysis to search for the presence of flavodoxin genes in other viruses. Putative flavodoxin genes were identified in the following phages: Psychrobacillus phage Perkons (34.7% identity; E-value 2 × 10−9); Enterococcus phage EF1 (31.3% identity; E-value 3 × 10−8); Bacillus phage vB_BcoS-136 (32% identity; E-value 3 × 10−7) and Staphylococcus phage vB_SepS_SEP9 (28.2% identity; E-value 8 × 10−4) and many identical or near identical isolates in Listeria phage (19 isolates; see also Table 1 for a list of phage genomes with encoded flavodoxin genes). We performed maximum likelihood (ML) phylogenetic analyses of 0305phi8-36 along with these additional phage proteins and the respective hosts of most of the phages with flavodoxins. Most of the phage proteins cluster together, in many cases with the host from which they were isolated (e.g., Listeria, Psychrobacillus). In other cases, notably 0305phi8-36, and Staphylococcus, Enterococcus, and Arcanobacterium phages, the phages are not clustered with their hosts, perhaps reflecting an origin in other soil bacteria (Figure 3).

4.2. Expression, Purification and Spectral Analysis of Recombinant 0305φ8-36 ORF223

The heterologous expression of histidine-tagged 0305φ8-36 ORF223 in E. coli and purification using Ni2+-NTA affinity chromatography yielded an average of 10 mg of purified flavodoxin protein per liter of culture. Cell fractionation by ultracentrifugation at 100,000× g revealed 0305φ8-36 ORF223 to be a soluble protein being located in the cytosolic fraction. No FMN-bound protein was detected in the isolated E. coli membranes. After two consecutive Ni2+-affinity chromatography purification steps a homogeneous band was observed on sodium dodecyl sulphate-polyacrylamide gel electrophoresis at approximately ~20 kDa, compared with the predicted molecular mass of the protein (16.5 kDa). This is probably due to inclusion of the engineered His tag into the recombinant protein, which results in mobility retardation in the PAGE gel. The UV-visible absorption spectrum of purified 0305φ8-36 ORF223 showed maxima at 360 and 460 nm with a shoulder at 490 nm typical of oxidized Flds (Figure 4). The spectra generated for 0305φ8-36 ORF223 protein was similar with those produced for control purified E. coli flavodoxin (Figure 4; [23,26]) and these characteristics are typical of other purified and spectrally characterized flavodoxin proteins described in the literature [11,28].

4.3. Characterization of the Catalytic Activity of 0305φ8-36 Flavodoxin in P450-Mediated Monooxygenation

In the soil bacterium S. coelicolor, cytochrome P450 170A1 (CYP170A1) catalyzes the sequential allylic oxidation of epi-isozizaene to the sequiterpene, albaflavenone (Figure 5A). Because the endogenous reductase complex is not known for CYP170A1, we previously used the E. coli flavodoxin and flavodoxin reductase system to drive this monooxygenase activity [26]. To investigate if 0305φ8-36 flavodoxin is functionally active we reconstituted S. coelicolor CYP170A1 enzymatic activity, replacing E. coli flavodoxin in the E. coli flavodoxin/flavodoxin system with the viral flavodoxin protein. Analysis of the CYP170A1-catalyzed reaction products by GC/MS revealed the time-dependent formation of two isomeric albaflavenol products with m/z 220 (tR 11.8 min; tR 12 min) and the ketone albaflaveone sequiterpene product with m/z 218 (tR 12.7 min) (Figure 5B). Using the E. coli flavodoxin/flavodoxin reductase couple, the turnover number was calculated to be ~0.36 ± 0.04 min−1 for the conversion of epi-isozizaene to albaflavenone. When the E. coli flavodoxin was replaced with 0305φ8-36 flavodoxin the turnover number was ~0.29 ± 0.03 min−1. Negative control incubations lacking CYP170A1, flavodoxins, or NADPH did not generate any products.

5. Discussion

When first discovered, flavodoxin was characterized as a key electron transfer protein functioning in photosynthetic cyanobacteria growing under limiting iron conditions [2,3]. Since that time flavodoxins have been discovered in numerous non-photosynthetic organisms, mainly Bacteria and in Archaea, where they function in diverse metabolic pathways. For example, in E.coli flavodoxin is involved in the biosynthesis of methionine through activation of cobalamine-dependent methionine synthase [60]; in B. subtilis flavodoxin plays a role in electron transfer to drive the cytochrome P450 enzyme (CYP107H1) involved in the biosynthesis of the water-soluble vitamin biotin [61]; and in the diazotroph Azotobacter vinelandii flavodoxin transfers reducing equivalents to the nitrogenase enzyme involved in nitrate reduction [62]. In the archaeon Methanosarcina acetivorans a flavodoxin plays a role in electron transfer in the metabolic pathway converting acetate to methane [63]. In bacterial cytochrome P450 systems flavodoxin, while driving P450 as a stand-alone protein, can be structurally linked to P450 in three ways: (i) fused to a P450 domain e.g., Rhodococcus rhodochrous XplA which is used in biotechnology to degrade the military explosive hexahydro-1,3,5-trinitro-1,3,5-triazine [64] (ii) fused to both the P450 domain and a flavodoxin reductase domain e.g., Rhodococcus ruber CYP116B3 a protein which can metabolize an array of hydrocarbons including the environmental pollutant xylene [65] and (iii) genetically organised in an operon with P450 and flavodoxin reductase e.g., Citrobacter braakii CYP176A1 catalyzes the conversion of cineole to 6β-hydroxycineole in energy metabolism [66]. In all cases genetic organisation is thought to favour optimal flavodoxin electron transfer during P450 catalysis.
Microsynteny analysis of phage flavodoxin genes reveals that in all cases that we examined, the flavodoxin gene is linked to or in the close neighbourhood of genes encoding α and β ribonucleotide reductase (Figure 6). This strongly suggests that the endogenous function of phage flavodoxin is to provide reducing equivalents to these enzymes during the catalytic removal of the 2′-hydroxyl group of the ribose ring of nucleoside diphosphates [67]. This reaction occurs during the formation of deoxyribonucleotides from ribonucleotides, which are necessary for DNA synthesis during phage replication.
0305φ8-36 is a bacteriophage that infects the soil bacterium B. thuringiensis [15,16]. Examining the metabolic pathways of this bacterium, our analysis of the sequenced genome of the B. thuringiensis reference strain Bt407 [68] revealed the presence of three cytochrome P450 genes: CYP102A5, CYP106B4 and CYP107J3. Furthermore, other strains of B. thuringiensis, infected by phage encoding flavodoxin [31], have been shown to encode P450 family members: CYP102 (CYP102A8); CYP106 (CYP106B1, B5); CYP107DY4; CYP107J (CYP107J2); CYP109T5 and CYP2467A1 [69] (D. Nelson personal communication). It is feasible that the infected B. thuringiensis strains can have their P450 enzyme activities enhanced, modified or inhibited by the phage flavodoxin. The phenotypic results will be dependent on the host P450 activity. CYP102A enzymes are generally self-sufficient P450s with a fused C-terminal reductase domain that are known to metabolize fatty acids and lipids, although no endogenous function for a CYP102 has been determined [24]. Similarly, the endogenous functions of the CYP106 and CYP107 enzymes are unknown. However, individual recombinant CYP106 and CYP107 enzymes have been shown to selectively hydroxylate a variety of steroid molecules including testosterone, progesterone and their derivatives [70]. Phages that encode flavodoxin can also infect other Bacillus species including B. cereus, B. subtilis and B. anthracis (Figure 6; [71]. Previous in silico analysis of P450s in 128 Bacillus species revealed the presence of 507 P450s grouped into 13 P450 families and 28 subfamilies with no P450 family found to be conserved across Bacillus species [72]. Bioinformatic analysis of those Bacillus P450 genes revealed that a limited number of P450 families dominate - CYP102, CYP106, CYP107, CYP109, CYP134, CYP152, and CYP197. Conversely, no P450 genes have been found to be encoded in any sequenced genome of a Listeria species [69]. However, of the more than 90 Staphylococcus species that have been sequenced, nine have been shown to encode P450. Staphylococcus P450s are limited to members of the CYP134 and CYP152 families [73]. Of note is P450 P450 OleTSA, a cytochrome P450 enzyme from Staphylococcus aureus that catalyzes the oxidative decarboxylation and hydroxylation of fatty acids to generate terminal alkenes and fatty alcohols [74]. It is possible that Staphylococcus phage flavodoxin can interact with S. aureus OleT and impact this catalytic activity.
Data are now revealing a remarkable diversity of metabolic genes in viruses, including many involved in sphingolipid biosynthesis, fermentation, and nitrogen metabolism [75,76]. For example, numerous giant viruses are predicted to encode the components of glycolysis and the tricarboxylic acid cycle, suggesting that they can re-program fundamental aspects of their host’s central carbon metabolism [77,78]. However, many of these viral enzymes also require cofactors to function. For example, aconitase and fumarase require Fe-S clusters—and other enzymes show other prosthetic group requirements (FAD/FMN/heme). This has led to a re-evaluation of the complex metabolic capabilities of viruses, particularly given their roles as important drivers of global bio-geochemical cycles [79].
In this study, we have found that a phage encoded FMN-containing flavodoxin can drive bacterial P450 activity. Similar phages infect different species of soil-inhabiting bacteria, mainly of the Bacillus genus and this suggests that the impacts of this metabolic activity may be widespread. Our results clearly show that viruses can engage with and possibly hijack host oxidative enzyme capacity. Historically, viruses were considered as metabolically inactive, non-living entities in comparison to the metabolic complexity displayed by cellular life. Their phenotypic impact was thought limited to merely affecting host mortality. Recent work has cast this aspersion aside, as we have come to realize that close metabolic integration appears to be required for viruses to subdue and coerce their unwitting hosts. This work clearly demonstrates the potential for viruses to not just manipulate ‘virus centric’ replication activities, but broader host metabolic activities and cellular function(s) as well. The large number of viral metabolic genes, including the different classes of redox proteins, requires further research into how these virus-specific enzymes affect both host physiology and the global microbiome and biogeochemical cycling.

Author Contributions

Conceptualization, D.C.L., S.L.K., J.V.G. and J.J.S.; methodology, D.C.L., J.V.G., B.Z., L.L., M.J.A., J.G.L.M., S.L.K. and J.J.S.; validation, D.C.L., J.V.G., B.Z., L.L., M.J.A., J.G.L.M., S.L.K. and J.J.S.; data analysis, D.C.L., J.V.G., B.Z., L.L., M.J.A., J.G.L.M., S.L.K. and J.J.S.; writing—original draft preparation, D.C.L., S.L.K., J.V.G. and J.J.S.; writing—review and editing, D.C.L., J.V.G., B.Z., L.L., M.J.A., J.G.L.M., S.L.K. and J.J.S.; project administration, D.C.L. and J.J.S. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by the National Institutes of Health grant 5U41HG003345 (J.V.G.), by the Woods Hole Center for Oceans and Human Health, NIH P01 ES021923 and NSF OCE-1314642 (J.J.S.), and by a Fulbright Scholarship (to D.C.L.). Funding at Swansea University supported by the European Regional Development Fund/Welsh European Funding Office via the BEACON project (S.L.K).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data used to support the findings of this research are available upon request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sancho, J. Flavodoxins: Sequence, folding, binding, function and beyond. Cell Mol. Life Sci. 2006, 63, 855–864. [Google Scholar] [CrossRef] [PubMed]
  2. Smillie, R.M. Isolation of two proteins with chloroplast ferredoxin activity from a blue-green alga. Biochem. Biophys. Res. Commun. 1965, 20, 621–629. [Google Scholar] [CrossRef]
  3. Knight, E.; D’Eustachio, A.J.; Hardy, R.W. Flavodoxin: A flavoprotein with ferredoxin activity from Clostrium pasteurianum. Biochim. Biophys. Acta 1966, 113, 626–628. [Google Scholar] [CrossRef]
  4. Anderson, R.F. Energetics of the one-electron reduction steps of riboflavin, FMN and FAD to their fully reduced forms. Biochim. Biophys. Acta 1983, 722, 158–162. [Google Scholar] [CrossRef]
  5. Mayhew, S.G. The effects of pH and semiquinone formation on the oxidation-reduction potentials of flavin mononucleotide. A reappraisal. Eur. J. Biochem. 1999, 265, 698–702. [Google Scholar] [CrossRef] [PubMed]
  6. Gaudu, P.; Weiss, B. Flavodoxin mutants of Escherichia coli K-12. J. Bacteriol. 2000, 182, 1788–1793. [Google Scholar] [CrossRef] [PubMed]
  7. Freigang, J.; Diederichs, K.; Schafer, K.P.; Welte, W.; Paul, R. Crystal structure of oxidized flavodoxin, an essential protein in Helicobacter pylori. Protein Sci. 2002, 11, 253–261. [Google Scholar] [CrossRef] [PubMed]
  8. Zhao, Q.; Modi, S.; Smith, G.; Paine, M.; McDonagh, P.D.; Wolf, C.R.; Tew, D.; Lian, L.Y.; Roberts, G.C.; Driessen, H.P. Crystal structure of the FMN-binding domain of human cytochrome P450 reductase at 1.93 A resolution. Protein Sci. 1999, 8, 298–306. [Google Scholar] [CrossRef]
  9. Lopez-Llano, J.; Maldonado, S.; Jain, S.; Lostao, A.; Godoy-Ruiz, R.; Sanchez-Ruiz, J.M.; Cortijo, M.; Fernandez-Recio, J.; Sancho, J. The long and short flavodoxins: II. The role of the differentiating loop in apoflavodoxin stability and folding mechanism. J. Biol. Chem. 2004, 279, 47184–47191. [Google Scholar] [CrossRef] [PubMed]
  10. Fukuyama, K.; Matsubara, H.; Rogers, L.J. Crystal structure of oxidized flavodoxin from a red alga Chondrus crispus refined at 1.8 A resolution. Description of the flavin mononucleotide binding site. J. Mol. Biol. 1992, 225, 775–789. [Google Scholar] [CrossRef]
  11. Peleato, M.L.; Ayora, S.; Inda, L.A.; Gomez-Moreno, C. Isolation and characterization of two different flavodoxins from the eukaryote Chlorella fusca. Biochem. J. 1994, 302, 807–811. [Google Scholar] [CrossRef] [PubMed]
  12. Dwivedi, B.; Xue, B.; Lundin, D.; Edwards, R.A.; Breitbart, M. A bioinformatic analysis of ribonucleotide reductase genes in phage genomes and metagenomes. BMC Evol. Biol. 2013, 13, 33. [Google Scholar] [CrossRef] [PubMed]
  13. Harrison, A.O.; Moore, R.M.; Polson, S.W.; Wommack, K.E. Reannotation of the Ribonucleotide Reductase in a Cyanophage Reveals Life History Strategies Within the Virioplankton. Front. Microbiol. 2019, 10, 134. [Google Scholar] [CrossRef]
  14. Lamb, D.C.; Follmer, A.H.; Goldstone, J.V.; Nelson, D.R.; Warrilow, A.G.; Price, C.L.; True, M.Y.; Kelly, S.L.; Poulos, T.L.; Stegeman, J.J. On the occurrence of cytochrome P450 in viruses. Proc. Natl. Acad. Sci. USA 2019, 116, 12343–12352. [Google Scholar] [CrossRef] [PubMed]
  15. Hardies, S.C.; Thomas, J.A.; Serwer, P. Comparative genomics of Bacillus thuringiensis phage 0305phi8-36: Defining patterns of descent in a novel ancient phage lineage. Virol. J. 2007, 4, 97. [Google Scholar] [CrossRef] [PubMed]
  16. Thomas, J.A.; Hardies, S.C.; Rolando, M.; Hayes, S.J.; Lieman, K.; Carroll, C.A.; Weintraub, S.T.; Serwer, P. Complete genomic sequence and mass spectrometric analysis of highly diverse, atypical Bacillus thuringiensis phage 0305phi8-36. Virology 2007, 368, 405–421. [Google Scholar] [CrossRef] [PubMed]
  17. Sievers, F.; Wilm, A.; Dineen, D.; Gibson, T.J.; Karplus, K.; Li, W.; Lopez, R.; McWilliam, H.; Remmert, M.; Soding, J.; et al. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol. Syst. Biol. 2011, 7, 539. [Google Scholar] [CrossRef]
  18. Chang, J.M.; Di Tommaso, P.; Notredame, C. TCS: A new multiple sequence alignment reliability measure to estimate alignment accuracy and improve phylogenetic tree reconstruction. Mol. Biol. Evol. 2014, 31, 1625–1637. [Google Scholar] [CrossRef]
  19. Hall, T.A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl. Acids Symp. Series 1999, 41, 95–98. [Google Scholar]
  20. Robert, X.; Gouet, P. Deciphering key features in protein structures with the new ENDscript server. Nucleic Acids Res. 2014, 42, 320–324. [Google Scholar] [CrossRef]
  21. Kozlov, A.; Darriba, D.; Flouri, T.; Morel, B.; Stamatakis, A. RAxML-NG: A fast, scalable, and user-friendly tool for maximum likelihood phylogenetic inference. Bioinformatics 2019, 35, 4453–4455. [Google Scholar] [CrossRef] [PubMed]
  22. Lemoine, F.; Domelevo Entfellner, J.B.; Wilkinson, E.; Correia, D.; Felipe, M.D.; De Oliveira, T.; Gascuel, O. Renewing Felsenstein’s phylogenetic bootstrap in the era of big data. Nature 2018, 556, 452–456. [Google Scholar] [CrossRef] [PubMed]
  23. Lamb, D.C.; Lei, L.; Warrilow, A.G.; Lepesheva, G.I.; Mullins, J.G.; Waterman, M.R.; Kelly, S.L. The first virally encoded cytochrome P450. J. Virol. 2009, 83, 8266–8269. [Google Scholar] [CrossRef] [PubMed]
  24. Lamb, D.C.; Lei, L.; Zhao, B.; Yuan, H.; Jackson, C.J.; Warrilow, A.G.; Skaug, T.; Dyson, P.J.; Dawson, E.S.; Kelly, S.L.; et al. Streptomyces coelicolor A3(2) CYP102 protein, a novel fatty acid hydroxylase encoded as a heme domain without an N-terminal redox partner. Appl. Environ. Microbiol. 2010, 76, 1975–1980. [Google Scholar] [CrossRef] [PubMed]
  25. Jenkins, C.M.; Waterman, M.R. Flavodoxin and NADPH-flavodoxin reductase from Escherichia coli support bovine cytochrome P450c17 hydroxylase activities. J. Biol. Chem. 1994, 269, 27401–27408. [Google Scholar] [CrossRef]
  26. Zhao, B.; Lin, X.; Lei, L.; Lamb, D.C.; Kelly, S.L.; Waterman, M.R.; Cane, D.E. Biosynthesis of the sesquiterpene antibiotic albaflavenone in Streptomyces coelicolor A3(2). J. Biol. Chem. 2008, 283, 8183–8189. [Google Scholar] [CrossRef]
  27. Omura, T.; Sato, R. The Carbon Monoxide-Binding Pigment of Liver Microsomes. I. Evidence for Its Hemoprotein Nature. J. Biol. Chem. 1964, 239, 2370–2378. [Google Scholar] [CrossRef]
  28. Jenkins, C.M.; Waterman, M.R. NADPH-flavodoxin reductase and flavodoxin from Escherichia coli: Characteristics as a soluble microsomal P450 reductase. Biochemistry 1998, 37, 6106–6113. [Google Scholar] [CrossRef]
  29. Zhao, B.; Guengerich, F.P.; Bellamine, A.; Lamb, D.C.; Izumikawa, M.; Lei, L.; Podust, L.M.; Sundaramoorthy, M.; Kalaitzis, J.A.; Reddy, L.M.; et al. Binding of two flaviolin substrate molecules, oxidative coupling, and crystal structure of Streptomyces coelicolor A3(2) cytochrome P450 158A2. J. Biol. Chem. 2005, 280, 11599–11607. [Google Scholar] [CrossRef]
  30. Bravo, A.; Likitvivatanavong, S.; Gill, S.S.; Soberon, M. Bacillus thuringiensis: A story of a successful bioinsecticide. Insect Biochem. Mol. 2011, 41, 423–431. [Google Scholar] [CrossRef]
  31. Sauder, A.B.; Quinn, M.R.; Brouillette, A.; Caruso, S.; Cresawn, S.; Erill, I.; Lewis, L.; Loesser-Casey, K.; Pate, M.; Scott, C.; et al. Genomic characterization and comparison of seven Myoviridae bacteriophage infecting Bacillus thuringiensis. Virology 2016, 489, 243–251. [Google Scholar] [CrossRef] [PubMed]
  32. Erill, I.; Caruso, S.M.; 2015 UMBC Phage Hunters. Genome sequence of Bacillus cereus group phage SalinJah. Genome Announc. 2016, 4, e00953-16. [Google Scholar] [CrossRef] [PubMed]
  33. Erill, I.; Caruso, S.M. Genome sequences of two Bacillus cereus group Bacteriophages, Eyuki and AvesoBmore. Genome Announc. 2015, 3, e01199-15. [Google Scholar] [CrossRef] [PubMed]
  34. Ting, J.H.; Smyth, T.B.; Chamakura, K.R.; Kuty Everett, G.F. Complete genome of Bacillus thuringiensis myophage BigBertha. Genome Announc. 2013, 1, e00853-13. [Google Scholar] [CrossRef]
  35. Maroun, J.W.; Whitcher, K.J.; Chamakura, K.R.; Kuty Everett, G.F. Complete genome of Bacillus thuringiensis myophage Spock. Genome Announc. 2013, 1, e00863-13. [Google Scholar] [CrossRef]
  36. Breslin, E.F.; Cornell, J.; Schuhmacher, Z.; Himelright, M.; Andos, A.; Childs, A.; Clem, A.; Gerber, M.; Gordillo, A.; Harb, L.; et al. Complete genome sequence of Bacillus phage Belinda from Grand Cayman Island. Genome Announc. 2016, 4, e00571-16. [Google Scholar] [CrossRef]
  37. Denes, T.; Vongkamjan, K.; Ackermann, H.W.; Moreno Switt, A.I.; Wiedmann, M.; den Bakker, H.C. Comparative genomic and morphological analyses of Listeria phages isolated from farm environments. Appl. Environ. Microbiol. 2014, 80, 4616–4625. [Google Scholar] [CrossRef]
  38. Peters, T.L.; Song, Y.; Bryan, D.W.; Hudson, L.K.; Denes, T.G. Mutant and recombinant phages selected from in vitro coevolution conditions overcome phage-resistant Listeria monocytogenes. Appl. Environ. Microbiol. 2020, 86, e02138-20. [Google Scholar] [CrossRef]
  39. Duperier, J.; Bulpitt, M.; Bispo, F.; Greguske, E. Genome annotations of two Bacillus phages, Tomato and Baseball Field. Microbiol. Resour. Announc. 2021, 10, e01196-20. [Google Scholar] [CrossRef]
  40. Peters, T.L.; Hudson, L.K.; Song, Y.; Denes, T.G. Complete genome sequences of two Listeria phages of the genus Pecentumvirus. Microbiol. Resour. Announc. 2019, 8, e01229-19. [Google Scholar] [CrossRef]
  41. Kostyk, N.; Chigbu, O.; Cochran, E.; Davis, J.; Essig, J.; Do, L.; Farooque, N.; Gbadamosi, Z.; Gnanodayan, A.; Hale, A.; et al. Complete genome sequences of Bacillus cereus group phages AaronPhadgers, ALPS, Beyonphe, Bubs, KamFam, OmnioDeoPrimus, Phireball, PPIsBest, YungSlug, and Zainny. Microbiol. Resour. Announc. 2021, 10, e0030021. [Google Scholar] [CrossRef] [PubMed]
  42. Nordmann, B.; Schilling, T.; Hoppert, M.; Hertel, R. Complete genome sequence of the virus isolate vB_BthM-Goe5 infecting Bacillus thuringiensis. Arch Virol. 2019, 164, 1485–1488. [Google Scholar] [CrossRef] [PubMed]
  43. Foltz, S.; Johnson, A.A. Complete genome sequences of nine Bacillus cereus group phages. Genome Announc. 2016, 4, e00473-16. [Google Scholar] [CrossRef] [PubMed]
  44. Flounlacker, K.; Miller, R.; Marquez, D.; Johnson, A.; the 2015–2016 VCU Phage Hunters. complete genome sequences of Bacillus phages DirtyBetty and Kida. Genome Announc. 2017, 5, e01385-16. [Google Scholar] [CrossRef]
  45. Erill, I.; Caruso, S.M. Complete genome sequence of Bacillus cereus group phage TsarBomba. Genome Announc. 2015, 3, e01178-15. [Google Scholar] [CrossRef]
  46. Lee, M.; Puglisi, K.M.; UMBC STEM-BUILD Cohort 1; Erill, I.; Caruso, S.M. Complete genome sequences of HonestAbe, Anthony, and Taffo16, three cluster C Bacillus cereus group bacteriophages. Genome Announc. 2018, 6, e00493-18. [Google Scholar] [CrossRef]
  47. Greguske, E.; Nadeau, A.; Fitzmeyer, E.; Fucikova, K. Genome sequence of Bacillus phage Saddex. Microbiol. Resour. Announc. 2018, 7, e01044-18. [Google Scholar] [CrossRef]
  48. Kent, B.; Raymond, T.; Mosier, P.D.; Johnson, A.A.; the 2016–2017 VCU Phage Hunters. Complete genome sequences of Bacillus phages Janet and OTooleKemple52. Genome Announc. 2018, 6, e00083-18. [Google Scholar] [CrossRef]
  49. Moreland, R.; Korn, A.; Newkirk, H.; Liu, M.; Gill, J.J.; Cahill, J.; Ramsey, J. Complete genome sequence of Staphylococcus aureus myophage Maine. Microbiol. Resour. Announc. 2019, 8, e01050-19. [Google Scholar] [CrossRef]
  50. Korn, A.M.; Hillhouse, A.E.; Sun, L.; Gill, J.J. Comparative genomics of three novel jumbo bacteriophages infecting Staphylococcus aureus. J. Virol. 2021, 95, e0239120. [Google Scholar] [CrossRef]
  51. Gutiérrez, D.; Vandenheuvel, D.; Martínez, B.; Rodríguez, A.; Lavigne, R.; García, P. Two phages, phiIPLA-RODI and phiIPLA-C1C, lyse mono- and dual-species Staphylococcal biofilms. Appl. Environ. Microbiol. 2015, 81, 3336–3348. [Google Scholar] [CrossRef] [PubMed]
  52. El Haddad, L.; Ben Abdallah, N.; Plante, P.L.; Dumaresq, J.; Katsarava, R.; Labrie, S.; Corbeil, J.; St-Gelais, D.; Moineau, S. Improving the safety of Staphylococcus aureus polyvalent phages by their production on a Staphylococcus xylosus strain. PLoS ONE 2014, 9, e102600. [Google Scholar]
  53. Alves, D.R.; Gaudion, A.; Bean, J.E.; Perez Esteban, P.; Arnot, T.C.; Harper, D.R.; Kot, W.; Hansen, L.H.; Enright, M.C.; Jenkins, A.T. Combined use of bacteriophage K and a novel bacteriophage to reduce Staphylococcus aureus biofilm formation. Appl. Environ. Microbiol. 2014, 80, 6694–6703. [Google Scholar] [CrossRef] [PubMed]
  54. Gu, J.; Liu, X.; Lu, R.; Li, Y.; Song, J.; Lei, L.; Sun, C.; Feng, X.; Du, C.; Yu, H.; et al. Complete genome sequence of Staphylococcus aureus bacteriophage GH15. J. Virol. 2012, 86, 8914–8915. [Google Scholar] [CrossRef]
  55. Crane, A.; Abaidoo, J.; Beltran, G.; Fry, D.; Furey, C.; Green, N.; Johal, R.; La Rosa, B.; Jimenez, C.L.; Luong, L.; et al. The complete genome sequence of the Staphylococcus bacteriophage Metroid. G3 2020, 10, 2975–2979. [Google Scholar] [CrossRef]
  56. Jun, J.W.; Giri, S.S.; Kim, H.J.; Chi, C.; Yun, S.; Kim, S.G.; Kim, S.W.; Kang, J.W.; Park, S.C. Complete genome sequence of the novel bacteriophage pSco-10 Infecting Staphylococcus cohnii. Genome Announc. 2017, 5, e01032-17. [Google Scholar] [CrossRef]
  57. Sáez Moreno, D.; Visram, Z.; Mutti, M.; Restrepo-Córdoba, M.; Hartmann, S.; Kremers, A.I.; Tišáková, L.; Schertler, S.; Wittmann, J.; Kalali, B.; et al. ε2-phages are naturally bred and have a vastly improved host range in Staphylococcus aureus over wild type phages. Pharmaceuticals 2021, 14, 325. [Google Scholar] [CrossRef]
  58. Oduor, J.M.O.; Kadija, E.; Nyachieo, A.; Mureithi, M.W.; Skurnik, M. Bioprospecting Staphylococcus phages with therapeutic and bio-control potential. Viruses 2020, 12, 133. [Google Scholar] [CrossRef]
  59. Abatángelo, V.; Peressutti Bacci, N.; Boncompain, C.A.; Amadio, A.F.; Carrasco, S.; Suárez, C.A.; Morbidoni, H.R. Broad-range lytic bacteriophages that kill Staphylococcus aureus local field strains. PLoS ONE 2017, 12, e0181671. [Google Scholar]
  60. Hoover, D.M.; Jarrett, J.T.; Sands, R.H.; Dunham, W.R.; Ludwig, M.L.; Matthews, R.G. Interaction of Escherichia coli cobalamin-dependent methionine synthase and its physiological partner flavodoxin: Binding of flavodoxin leads to axial ligand dissociation from the cobalamin cofactor. Biochemistry 1997, 36, 127–138. [Google Scholar] [CrossRef]
  61. Green, A.J.; Rivers, S.L.; Cheeseman, M.; Reid, G.A.; Quaroni, L.G.; Macdonald, I.D.; Chapman, S.K.; Munro, A.W. Expression, purification and characterization of cytochrome P450 Biol: A novel P450 involved in biotin synthesis in Bacillus subtilis. J. Biol. Inorg. Chem. 2001, 6, 523–533. [Google Scholar] [CrossRef]
  62. Gangeswaran, R.; Eady, R.R. Flavodoxin 1 of Azotobacter vinelandii: Characterization and role in electron donation to purified assimilatory nitrate reductase. Biochem. J. 1996, 317, 103–108. [Google Scholar] [CrossRef]
  63. Prakash, D.; Iyer, P.R.; Suharti, S.; Walters, K.A.; Santiago-Martinez, M.G.; Golbeck, J.H.; Murakami, K.S.; Ferry, J.G. Structure and function of an unusual flavodoxin from the domain Archaea. Proc. Natl. Acad. Sci. USA 2019, 116, 25917–25922. [Google Scholar] [CrossRef]
  64. Jackson, R.G.; Rylott, E.L.; Fournier, D.; Hawari, J.; Bruce, N.C. Exploring the biochemical properties and remediation applications of the unusual explosive-degrading P450 system XplA/B. Proc. Natl. Acad. Sci. USA 2007, 104, 16822–16827. [Google Scholar] [CrossRef]
  65. Liu, L.; Schmid, R.D.; Urlacher, V.B. Cloning, expression, and characterization of a self-sufficient cytochrome P450 monooxygenase from Rhodococcus ruber DSM 44319. Appl. Microbiol. Biotechnol. 2006, 72, 876–882. [Google Scholar] [CrossRef]
  66. Hawkes, D.B.; Adams, G.W.; Burlingame, A.L.; Ortiz de Montellano, P.R.; De Voss, J.J. Cytochrome P450(cin) (CYP176A), isolation, expression, and characterization. J. Biol. Chem. 2002, 277, 27725–27732. [Google Scholar] [CrossRef]
  67. Jordan, A.; Reichard, P. Ribonucleotide reductases. Annu. Rev. Biochem. 1998, 67, 71–98. [Google Scholar] [CrossRef]
  68. Sheppard, A.E.; Poehlein, A.; Rosenstiel, P.; Liesegang, H.; Schulenburg, H. Complete Genome Sequence of Bacillus thuringiensis Strain 407 Cry. Genome Announc. 2013, 1, e00158-12. [Google Scholar] [CrossRef]
  69. Mthethwa, B.C.; Chen, W.; Ngwenya, M.L.; Kappo, A.P.; Syed, P.R.; Karpoormath, R.; Yu, J.H.; Nelson, D.R.; Syed, K. Comparative Analyses of Cytochrome P450s and Those Associated with Secondary Metabolism in Bacillus Species. Int. J. Mol. Sci. 2018, 19, 3623. [Google Scholar] [CrossRef]
  70. Szaleniec, M.; Wojtkiewicz, A.M.; Bernhardt, R.; Borowski, T.; Donova, M. Bacterial steroid hydroxylases: Enzyme classes, their functions and comparison of their catalytic mechanisms. Appl. Microbiol. Biotechnol. 2018, 102, 8153–8171. [Google Scholar] [CrossRef]
  71. Grose, J.H.; Jensen, G.L.; Burnett, S.H.; Breakwell, D.P. Genomic comparison of 93 Bacillus phages reveals 12 clusters, 14 singletons and remarkable diversity. BMC Genom. 2014, 15, 855. [Google Scholar]
  72. Furuya, T.; Shibata, D.; Kino, K. Phylogenetic analysis of Bacillus P450 monooxygenases and evaluation of their activity towards steroids. Steroids 2009, 74, 906–912. [Google Scholar] [CrossRef]
  73. Padayachee, T.; Nzuza, N.; Chen, W.; Nelson, D.R.; Syed, K. Impact of lifestyle on cytochrome P450 monooxygenase repertoire is clearly evident in the bacterial phylum Firmicutes. Sci. Rep. 2020, 10, 13982. [Google Scholar] [CrossRef]
  74. Jiang, Y.; Li, Z.; Wang, C.; Zhou, Y.J.; Xu, H.; Li, S. Biochemical characterization of three new α-olefin-producing P450 fatty acid decarboxylases with a halophilic property. Biotechno. Biofuels 2019, 12, 79. [Google Scholar] [CrossRef]
  75. Vardi, A.; Haramaty, L.; Van Mooy, B.A.; Fredricks, H.F.; Kimmance, S.A.; Larsen ABidle, K.D. Host-virus dynamics and subcellular controls of cell fate in a natural coccolithophore population. Proc. Natl. Acad. Sci. USA 2012, 109, 19327–19332. [Google Scholar] [CrossRef]
  76. Monier, A.; Chambouvet, A.; Milner, D.S.; Attah, V.; Terrado, R.; Lovejoy, C.; Moreau, H.; Santoro, A.E.; Derelle, E.; Richards, T.A. Host-derived viral transporter protein for nitrogen uptake in infected marine phytoplankton. Proc. Natl. Acad. Sci. USA 2017, 114, 489–498. [Google Scholar] [CrossRef]
  77. Moniruzzaman, M.; Martinez-Gutierrez, C.A.; Weinheimer, A.R.; Aylward, F.O. Dynamic genome evolution and complex virocell metabolism of globally-distributed giant viruses. Nat. Commun. 2020, 11, 1710. [Google Scholar] [CrossRef]
  78. Blanc-Mathieu, R.; Dahle, H.; Hofgaard, A.; Brandt, D.; Ban, H.; Kalinowski, J.; Ogata, H.; Sandaa, R.A. A persistent giant algal virus, with a unique morphology, encodes an unprecedented number of genes involved in energy metabolism. J. Virol. 2021, 95, e02446-20. [Google Scholar] [CrossRef]
  79. Brahim Belhaouari, D.; Pires De Souza, G.A.; Lamb, D.C.; Kelly, S.L.; Goldstone, J.V.; Stegeman, J.J.; Colson, P.; La Scola, B.; Aherfi, S. Metabolic arsenal of giant viruses: Host hijack or self-use? eLife 2022, 11, e78674. [Google Scholar] [CrossRef]
Biomolecules 12 01107 g001 550
Figure 1. Schematic diagram showing the genomic region of bacteriophage 0305φ8-36 where the flavodoxin gene is located. The scale is in kilobase pairs. Each arrow represents an individual orf. The flavodoxin gene, colored green, is 0305φ8-36 ORF number 223. ORFs colored yellow and blue have been tentatively assigned dCMP deaminase and α and β-ribonucleotide reductase (nrd), respectively. Arrows shown in grey represent genes predicted to encode virion proteins of unknown function.
Figure 1. Schematic diagram showing the genomic region of bacteriophage 0305φ8-36 where the flavodoxin gene is located. The scale is in kilobase pairs. Each arrow represents an individual orf. The flavodoxin gene, colored green, is 0305φ8-36 ORF number 223. ORFs colored yellow and blue have been tentatively assigned dCMP deaminase and α and β-ribonucleotide reductase (nrd), respectively. Arrows shown in grey represent genes predicted to encode virion proteins of unknown function.
Biomolecules 12 01107 g001
Biomolecules 12 01107 g002 550
Figure 2. Amino acid sequence alignment of 0305φ8-36 phage flavodoxin with selected cellular flavodoxins. Alignment was generated using the Clustal Omega program via the EBI (https://www.ebi.ac.uk/Tools/msa/clustalo/, accessed on 22 July 2022). The flavodoxin sequences used were: long chain flavodoxins, E. coli (strain K12) FldA (accession number: P61949), E. coli (strain K12) FplB (P0ABY4) and short chain flavodoxins, B. subtilis (strain 168) ykuN (O34737), B. subtilis ykuP (strain 168) (O34589), Methanosarcina acetivorans (strain ATCC 35395) (Q8TPV5), B. thuringiensis serovar israelensis (strain ATCC 35646) (Q3ESQ3), Azotobacter vinelandii (strain DJ/ATCC BAA-1303) (P52964), Clostridium beijerinckii (Clostridium MP) (P00322) and Synechocystis sp. (strain PCC 6803/Kazusa) (P27319). Residues indicated by asterisks are conserved in all sequences. Residues indicated by single and double dots are of similar and highly similar chemical character, respectively, according to criteria set in the program. The residues comprising the three FMN-binding loops are underlined in red.
Figure 2. Amino acid sequence alignment of 0305φ8-36 phage flavodoxin with selected cellular flavodoxins. Alignment was generated using the Clustal Omega program via the EBI (https://www.ebi.ac.uk/Tools/msa/clustalo/, accessed on 22 July 2022). The flavodoxin sequences used were: long chain flavodoxins, E. coli (strain K12) FldA (accession number: P61949), E. coli (strain K12) FplB (P0ABY4) and short chain flavodoxins, B. subtilis (strain 168) ykuN (O34737), B. subtilis ykuP (strain 168) (O34589), Methanosarcina acetivorans (strain ATCC 35395) (Q8TPV5), B. thuringiensis serovar israelensis (strain ATCC 35646) (Q3ESQ3), Azotobacter vinelandii (strain DJ/ATCC BAA-1303) (P52964), Clostridium beijerinckii (Clostridium MP) (P00322) and Synechocystis sp. (strain PCC 6803/Kazusa) (P27319). Residues indicated by asterisks are conserved in all sequences. Residues indicated by single and double dots are of similar and highly similar chemical character, respectively, according to criteria set in the program. The residues comprising the three FMN-binding loops are underlined in red.
Biomolecules 12 01107 g002
Biomolecules 12 01107 g003 550
Figure 3. Relationship of 0305φ8-36 flavodoxin with other phage and bacterial host flavodoxins. Phylogenetic tree of 0305φ8-36 flavodoxin with additional phage flavodoxins and the respective hosts of most of the phages with flavodoxins. Numbers indicate the bootstrap probability values for the branch topology shown.
Figure 3. Relationship of 0305φ8-36 flavodoxin with other phage and bacterial host flavodoxins. Phylogenetic tree of 0305φ8-36 flavodoxin with additional phage flavodoxins and the respective hosts of most of the phages with flavodoxins. Numbers indicate the bootstrap probability values for the branch topology shown.
Biomolecules 12 01107 g003
Biomolecules 12 01107 g004 550
Figure 4. UV-visible absorbance spectra of purified E. coli and 0305φ8-36 phage flavodoxins. (A) Oxidized spectrum obtained for purified E. coli flavodoxin (15 μM). (B) Oxidized spectrum obtained for purified phage 0305φ8-36 flavodoxin (15 μM).
Figure 4. UV-visible absorbance spectra of purified E. coli and 0305φ8-36 phage flavodoxins. (A) Oxidized spectrum obtained for purified E. coli flavodoxin (15 μM). (B) Oxidized spectrum obtained for purified phage 0305φ8-36 flavodoxin (15 μM).
Biomolecules 12 01107 g004
Biomolecules 12 01107 g005 550
Figure 5. Epi-isozizaene allylic oxidation by CYP170A1 supported by 0305φ8-36 flavodoxin. (A) Metabolic pathway for the formation of albaflavenone catalysed by CYP170A1. (B) In vitro catalytic activity of CYP170A1 supported by 0305φ8-36 and E. coli flavodoxin reductase. Product profile was obtained using individual component concentrations as follows: 1 nmol CYP170A1, 10 nmol 0305φ8-36 flavodoxin or 10 nmol E. coli flavodoxin, 2.0 nmol E. coli flavodoxin reductase and 40 nmol epi-isozizaene. Three products are noted, 2 and 3 (4(R)-albaflavenol and 4(S)-albaflavenol epimers respectively) and 4 (albaflavenone), respectively; substrate epi-isozizaene as 1. Experiments were performed in triplicate. In negative control incubations (without CYP170A1 or without flavodoxin/flavodoxin reductase), no product formation was observed (data not shown).
Figure 5. Epi-isozizaene allylic oxidation by CYP170A1 supported by 0305φ8-36 flavodoxin. (A) Metabolic pathway for the formation of albaflavenone catalysed by CYP170A1. (B) In vitro catalytic activity of CYP170A1 supported by 0305φ8-36 and E. coli flavodoxin reductase. Product profile was obtained using individual component concentrations as follows: 1 nmol CYP170A1, 10 nmol 0305φ8-36 flavodoxin or 10 nmol E. coli flavodoxin, 2.0 nmol E. coli flavodoxin reductase and 40 nmol epi-isozizaene. Three products are noted, 2 and 3 (4(R)-albaflavenol and 4(S)-albaflavenol epimers respectively) and 4 (albaflavenone), respectively; substrate epi-isozizaene as 1. Experiments were performed in triplicate. In negative control incubations (without CYP170A1 or without flavodoxin/flavodoxin reductase), no product formation was observed (data not shown).
Biomolecules 12 01107 g005
Biomolecules 12 01107 g006 550
Figure 6. Microsynteny around phage flavodoxin genes. Comparison of the genetic loci and surrounding genes for flavodoxins in Bacillus phages, Listeria phages and Staphylococcus phages. Each phage flavodoxin gene (fld) is colored green and α and β ribonucleotide reductase genes (α and β nrb) are colored blue. In each Bacillus, Listeria and Staphylococcus phage genomes examined the fld gene was next to or close to α and β nrb genes, strongly suggesting functional linkage.
Figure 6. Microsynteny around phage flavodoxin genes. Comparison of the genetic loci and surrounding genes for flavodoxins in Bacillus phages, Listeria phages and Staphylococcus phages. Each phage flavodoxin gene (fld) is colored green and α and β ribonucleotide reductase genes (α and β nrb) are colored blue. In each Bacillus, Listeria and Staphylococcus phage genomes examined the fld gene was next to or close to α and β nrb genes, strongly suggesting functional linkage.
Biomolecules 12 01107 g006
Table
Table 1. Flavodoxin genes found in Bacillus, Listeria and Staphylococcus phage genomes.
Table 1. Flavodoxin genes found in Bacillus, Listeria and Staphylococcus phage genomes.
ClassPhagePhage HostSample LocationPredicted Phage ORFsFlavodoxin
aa 		Genbank
Identifier		Ref.
Bacillus0305φ8-36 Bacillus thuringiensissoil246148ABS83781.1[16]
Hakuna Bacillus thuringiensissoil294159YP_009036585.1[31]
Megatron Bacillus thuringiensissoil290159YP_009036208.1[31]
Riley Bacillus thuringiensissoil290149YP_009055891.1[31]
CAM003 Bacillus thuringiensissoil287150YP_009037037.1[31]
EvoliBacillus thuringiensissoil294150YP_009035659.1[31]
Hoody T Bacillus thuringiensissoil270153YP_009035330.1[31]
TrollBacillus thuringiensissoil289149YP_008430917.1[31]
SalinJahBacillus cereussoil292149ANH50605.1[32]
EyukiBacillus cereussoil300159YP_009212082.1[33]
AvesoBmoreBacillus thuringiensissoil301149YP_009206486.1[33]
BigBerthaBacillus thuringiensissoil291153YP_008771155.1[34]
SpockBacillus thuringiensissoil283149YP_008770352.1[35]
BelindaBacillus thuringiensissoil295159ANM46069.1[36]
TomatoBacillus thuringiensissoil200159QLF85935.1[37]
AaronPhadgersBacillus thuringiensissoil301159ASR78796.1[38]
BeyonpheBacillus cereuswater300149QDH49824.1[38]
BubsBacillus thuringiensissoil302159ASR78601.1[38]
KamfamBacillus thuringiensissoil293150AXQ67322.1[38]
PIPsBestBacillus thuringiensissoil301159ASR78379.1[38]
ALPSBacillus thuringiensissoil295159QDH50122.1[38]
PhireballBacillus thuringiensissoil299159QDH49415.1[38]
OmnioDeoPrimusBacillus thuringiensissoil299159AXQ67477.1[38]
ZainnyBacillus thuringiensissoil303159ASR79375.1[38]
vB_BthM-Goe5Bacillus thuringiensissoil272150AZF89241.1[39]
SageFayge Bacillus thuringiensissoil300159YP_009280944.1[40]
NemoBacillus thuringiensissoil301159YP_009287019.1[40]
NigalanaBacillus thuringiensissoil302159YP_009282535.1[40]
DIGNKCBacillus thuringiensissoil291159AMW62869.1[40]
ZukoBacillus thuringiensissoil294159AMW62552.1[40]
PhrodoBacillus thuringiensissoil288149YP_009290006.1[40]
NotTheCreekBacillus thuringiensissoil296159YP_009284467.1[40]
JugloneBacillus thuringiensissoil293149AMW61744.1[40]
VinnyBacillus thuringiensissoil297150AMW61891.1[40]
DirtyBettyBacillus thuringiensissoil302159YP_009285086.1[41]
KidaBacillus thuringiensissoil304159YP_009279309.1[41]
TsarBombaBacillus thuringiensissoil247151YP_009206942.1[42]
HonestAbeBacillus thuringiensissoil286159AUV57777.1[43]
Taffo16Bacillus thuringiensissoil284149ASZ75860.1[43]
AnthonyBacillus thuringiensissoil279150ASU00988.1[43]
SaddexBacillus thuringiensissoil208159AXF41931.1[44]
JanetBacillus thuringiensissoil285150ASR79938.1[45]
OTooleKemple52Bacillus thuringiensissoil291153ASR79614.1[45]
SBP8aBacillus cereussoil298146AOZ62389.1DS *
QCM8 Bacillus cereussoil288151AOZ62044.1DS *
BJ4Bacillus cereussoil298146AOZ61763.1DS *
BastilleBacillus cereussoil273150YP_006907364.1DS *
ChotacabrasBacillus thuringiensissoil285149QEM43180.1DS *
vB_BcoS-136Bacillus cohniilake sediment238146AYP68319.1DS *
FlapjackBacillus thuringiensissoil288149ARQ95040.1DS *
vB_BanH_
RonSwanson		Bacillus anthracissoil-202UGO50429.1DS *
vB_BanH_
Emiliahah		Bacillus anthracissoil-174UGO49205.1DS *
vB_BanH_JarJarBacillus anthracissoil-202UGO48939.1DS *
vB_BanH_
McCartney		Bacillus anthracissoil-174UGO47691.1DS *
vB_BanH_
Abinadi		Bacillus anthracissoil-159UGO46389.1DS *
AnthosBacillus thuringiensissoil290149QPY77365.1DS *
SmudgeBacillus thuringiensissoil292159ANI24755.1DS *
BC-T25Bacillus sp.soil236151QEG04194.1DS *
ListeriavB_LmoM_AG20Listeria monocytogenesabbatoir178141YP_007676786.1DS *
LP-124Listeria monocytogenessilage188150YP_009784573.1[46]
LP-064Listeria monocytogenessilage188150YP_009592673.1[46]
LP-083-2Listeria monocytogenessilage189150YP_009044596.1[46]
LP-048Listeria monocytogenessilage177150YP_009042933.1[46]
LP-125Listeria monocytogenessilage189150YP_008240106.1[46]
LP-Mix_6.2Listeria monocytogeneslaboratory197150QNL32088.1[47]
LP-Mix_6.1Listeria monocytogeneslaboratory195150QNL31890.1[47]
LP-039Listeria monocytogenessilage201150QEP53123.2[48]
LP-066Listeria monocytogenessilage189150QDK04972.2[48]
LMTA-94Listeria monocytogeneslaboratory189150AID17150.1DS
StaphylococcusMaineStaphylococcus aureuspig barn swab219130QEM41386.1[49]
MarsHillStaphylococcus aureuspig barn swab262134QQM14610.1[50]
MadawaskaStaphylococcus aureuspig barn swab264134QQO92731.1[50]
MachiasStaphylococcus aureuspig barn swab263140QQO92468.1[50]
phiIPLA-RODIStaphylococcus aureussewage213139AJA42083.1[51]
phiIPLA-C1CStaphylococcus aureussewage203132YP_009214514.1[51]
Team1Staphylococcus aureushospital217143YP_009098273.1[52]
MCE-2014Staphylococcus aureussewage204130P_009098058.1[53]
GH15Staphylococcus aureussewage214130YP_007002257.1[54]
MetroidStaphylococcus sp.soil254143KE56191.1[55]
pSco-10Staphylococcus cohniiduck feces131130ANH50479.1[56]
BT3Staphylococcus aureusanimal232139QVD58098.1[57]
Stab23Staphylococcus xylosussewage247130VEV88569.1[58]
vB_Sau_S24Staphylococcus aureussoil209130ARM69410.1[59]
vB_Sau_Clo6Staphylococcus aureussewage213130ARM69197.1[59]
vB_SsapH-Golestan-100Staphylococcus saprophyticus-192133BDA81541.1DS *
phiSA_BS1Staphylococcus sp.dairy farm200129YP_009799552.1DS *
vB_SsapH-Golestan-105-MStaphylococcus saprophyticus-203149BDA82285.1DS *
SA3Staphylococcus aureussewage223130ASZ78055.1DS *
PALS_1Staphylococcus aureusanimal191130QDJ97591.1DS *
* DS direct submission to Genbank.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lamb, D.C.; Goldstone, J.V.; Zhao, B.; Lei, L.; Mullins, J.G.L.; Allen, M.J.; Kelly, S.L.; Stegeman, J.J. Characterization of a Virally Encoded Flavodoxin That Can Drive Bacterial Cytochrome P450 Monooxygenase Activity. Biomolecules 2022, 12, 1107. https://doi.org/10.3390/biom12081107

AMA Style

Lamb DC, Goldstone JV, Zhao B, Lei L, Mullins JGL, Allen MJ, Kelly SL, Stegeman JJ. Characterization of a Virally Encoded Flavodoxin That Can Drive Bacterial Cytochrome P450 Monooxygenase Activity. Biomolecules. 2022; 12(8):1107. https://doi.org/10.3390/biom12081107

Chicago/Turabian Style

Lamb, David C., Jared V. Goldstone, Bin Zhao, Li Lei, Jonathan G. L. Mullins, Michael J. Allen, Steven L. Kelly, and John J. Stegeman. 2022. "Characterization of a Virally Encoded Flavodoxin That Can Drive Bacterial Cytochrome P450 Monooxygenase Activity" Biomolecules 12, no. 8: 1107. https://doi.org/10.3390/biom12081107

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop