Next Article in Journal
Assessing the Interactions of Statins with Human Adenylate Kinase Isoenzyme 1: Fluorescence and Enzyme Kinetic Studies
Next Article in Special Issue
Ionic Liquids—A Review of Their Toxicity to Living Organisms
Previous Article in Journal
Comparative Antimicrobial Activity of Hp404 Peptide and Its Analogs against Acinetobacter baumannii
Previous Article in Special Issue
Bisphenol A Inhibits the Transporter Function of the Blood-Brain Barrier by Directly Interacting with the ABC Transporter Breast Cancer Resistance Protein (BCRP)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ancient Bacterial Class Alphaproteobacteria Cytochrome P450 Monooxygenases Can Be Found in Other Bacterial Species

1
Department of Biochemistry and Microbiology, Faculty of Science and Agriculture, University of Zululand, KwaDlangezwa 3886, South Africa
2
Department of Pharmaceutical Chemistry, College of Health Sciences, University of KwaZulu-Natal, Durban 4000, South Africa
3
Faculty of Chemistry, Biological and Chemical Research Center, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
4
Department of Molecular Microbiology and Genetics, University of Göttingen, 37077 Göttingen, Germany
5
Department of Microbiology, Immunology and Biochemistry, University of Tennessee Health Science Center, Memphis, TN 38163, USA
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(11), 5542; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22115542
Submission received: 21 April 2021 / Revised: 19 May 2021 / Accepted: 20 May 2021 / Published: 24 May 2021
(This article belongs to the Collection Feature Papers in Molecular Toxicology)

Abstract

:
Cytochrome P450 monooxygenases (CYPs/P450s), heme-thiolate proteins, are well-known players in the generation of chemicals valuable to humans and as a drug target against pathogens. Understanding the evolution of P450s in a bacterial population is gaining momentum. In this study, we report comprehensive analysis of P450s in the ancient group of the bacterial class Alphaproteobacteria. Genome data mining and annotation of P450s in 599 alphaproteobacterial species belonging to 164 genera revealed the presence of P450s in only 241 species belonging to 82 genera that are grouped into 143 P450 families and 214 P450 subfamilies, including 77 new P450 families. Alphaproteobacterial species have the highest average number of P450s compared to Firmicutes species and cyanobacterial species. The lowest percentage of alphaproteobacterial species P450s (2.4%) was found to be part of secondary metabolite biosynthetic gene clusters (BGCs), compared other bacterial species, indicating that during evolution large numbers of P450s became part of BGCs in other bacterial species. Our study identified that some of the P450 families found in alphaproteobacterial species were passed to other bacterial species. This is the first study to report on the identification of CYP125 P450, cholesterol and cholest-4-en-3-one hydroxylase in alphaproteobacterial species (Phenylobacterium zucineum) and to predict cholesterol side-chain oxidation capability (based on homolog proteins) by P. zucineum.

1. Introduction

Cytochrome P450 monooxygenases (CYPs/P450s) are heme-thiolate proteins known to humans for more than five and a half decades [1]. P450s perform enzymatic reactions with stereo- and regio-specific oxidation and because of this capability these enzymes are subjected to various applications in all forms of biology [2,3]. Most living organisms and some non-living entities such as viruses have been found to have P450s in their genomes [4,5], indicating their important role in primary and secondary metabolism.
P450 reactions are critical in determining the drug toxicity of prodrugs and xenobiotic compounds [6]. P450s are well-known drug targets for some of the microbial pathogens [7,8]. Some of the P450s from microbes have been characterized and their application in the production of chemicals that are valuable to humans, such as secondary metabolites (compounds produced by the organisms that have no role in their growth, development or reproduction) has been explored [3,9,10,11,12,13]. Well-known examples of P450s’ involvement are antibiotic production, such as in the biosynthesis of the macrolide antibiotic erythromycin and glycopeptide antibiotics [14,15], the production of the anticancer drugs Taxol and perillyl alcohol [16,17] and the production of pharmaceutical compounds such as pravastatin [13].
P450s’ critical role in the synthesis of secondary metabolites has been thoroughly reviewed [13,18,19]. It was found that the coding sequences (genes) of enzymes involved in the synthesis of different secondary metabolites, including P450s, were part of gene clusters named biosynthetic gene clusters (BGCs) [20], where P450s were found to be critical in contributing to the diversity of the secondary metabolites [13,18,19]. In the pre-genome sequencing era, identification of such gene clusters and P450s that formed part of the clusters required tedious and laborious work. The current genome sequencing era has revolutionized the process and bioinformatics programs are available for the identification of BGCs in organisms [20]. However, to date, identification and annotation of P450s still involve manual sorting and naming as per the International P450 Nomenclature Committee rule, i.e., sequences with >40% identity are assigned to the same family as the named homolog P450 and sequences with >55% identity are assigned to the same subfamily as the named homolog P450 [21,22,23]. Proteins with <40% identity are assigned to a new P450 family. The current genome sequencing era has offered the opportunity to understand P450s’ profiles and thus evolution patterns, particularly in bacterial species where the impact of lifestyle has been found to affect the P450 profiles in an organism [24] profoundly.
The bacterial class Alphaproteobacteria contains the most abundant, extraordinarily diverse and ancient group of bacteria [25,26]. Alphaproteobacteria consist of species with a diversity of lifestyles, geographical distribution and genome size [26]. They inhabit diverse ecological niches, from water to soil, and form both extra- and intracellular associations with eukaryotes, including unicellular organisms, and multicellular plants and mammals [25,26]. Some species are in symbiotic relationships with plants such as species belonging to the genus Agrobacterium; species belonging to Rickettsiaceae, Brucellaceae and Bartonellaceae are human pathogens [27], of which a few are of ecological importance, playing a huge role in carbon, nitrogen and sulfur cycles [28]. These bacteria are found in large amounts in the upper layer of the ocean [29]. In spite of their differences, all alphaproteobacterial species are oligotrophs; therefore, they can survive in an environment with low nutrients [30]. Alphaproteobacteria contain species that have biotechnological significance, especially some of the species found to produce secondary metabolites, compounds that play no role in the growth or reproduction of organisms, but give the organism a selective advantage [31,32]. The biological properties of these secondary metabolites in the field of drug discovery are invaluable, and strenuous efforts to find new secondary metabolites with biotechnological potentials are being made around the world. Some of the secondary metabolites produced by alphaproteobacterial species and their biological function are presented in Table 1.
It is a well-known fact that P450s are involved in the production of secondary metabolites per se; they contribute to the diversity of secondary metabolites produced by an organism [13,18,19]. Recent genome data mining, annotation of P450s and analysis of P450s’ association with secondary metabolism in some bacterial species such as mycobacterial species [39], Streptomyces species [40], cyanobacterial species [41] and Firmicutes species [24] revealed that quite a large number of P450s are involved in the production of secondary metabolites. Apart from these bacterial species, none of the other bacterial species’ P450 repertoire has been analyzed. In addition to this, to date, the P450 repertoire of alphaproteobacterial species and their association with secondary metabolism has not been reported despite the fact that P450s, in general, are involved in the production of secondary metabolites and these species do produce secondary metabolites. Thus, in this study, we address these research gaps by performing genome-wide data mining, annotation and phylogenetic analysis of P450s and identification of P450s involved in secondary metabolism in the bacterial class Alphaproteobacteria.

2. Results and Discussion

2.1. Only 40% of Alphaproteobacterial Species Have P450s

Genome data mining and annotation of P450s in 599 alphaproteobacterial species belonging to 164 genera (Table S1) revealed the presence of P450s in only 241 species belonging to 82 genera (Figure 1). This indicates that only 40% of alphaproteobacterial species and 50% of Alphaproteobacteria genera have P450s (Figure 1). All the species analyzed in the study belonging to genera such as Agrobacterium, Rhodopseudomonas, Hyphomicrobium, Ketogulonicigenium, Filomicrobium, Phenylobacterium, Roseobacter, Gluconacetobacter and Nitrobacter, and most of the species belonging to genera such as Sinorhizobium, Rhizobium, Bradyrhizobium and Sphingobium, have P450s in their genomes (Table S1). A significant number of species belonging to genera such as Rickettsia, Bartonella, Ehrlichia, Wolbachia and Anaplasma were analyzed in this study; no P450s were found in these species, suggesting that species in these genera probably do not have P450s (Table S1). Detailed information on genera, species and P450 information is presented in Table S1. In total, 874 P450s apart from 52 short-P450s were identified in 241 alphaproteobacterial species (Table S2 and Supplementary Dataset 1). On average, four P450s were found in 241 alphaproteobacterial species; 65 species had a single P450 in their genome (Table S1). Of these 65 species, 30 were from the genus Brucella alone (Table S1). Comparative analysis of P450s in alphaproteobacterial species revealed that Bradyrhizobium oligotrophicum has the highest number of P450s (17 P450s) in its genome, followed by 16 P450s in the species Bradyrhizobium japonicum E109, Novosphingobium aromaticivorans and Sphingomonas wittichii and 15 P450s in Bradyrhizobium diazoefficiens USDA 110 (Table S2). Comparative analysis with other bacterial species revealed that alphaproteobacterial species have a higher average number of P450s compared to Firmicutes species [24] and cyanobacterial species [41], but a lower number compared to Streptomyces species [39,40] and mycobacterial species [42] (Table 2). A list of P450s and their sequences along with short-P450s is presented in Supplementary Dataset 1.

2.2. Alphaproteobacterial Species Have the Highest Number of P450 Families Next to Streptomyces

Following the International P450 Nomenclature Committee rules, 874 alphaproteobacterial species’ P450s were grouped into 143 P450 families and 214 P450 subfamilies (grouped in the same family when they share >40% and the same subfamily when they share >55%; with less than 40% identity species are assigned to a new P450 family) [21,22], based on the phylogenetic analysis of P450s (Figure 2) [22,44,45]. A list of P450 families and subfamilies, their count and percentage contribution to the total number of P450s is presented in Table 3.
P450 family and subfamily-level comparative analysis with other bacterial species revealed that alphaproteobacterial species have more P450 families and subfamilies compared to Firmicutes species, mycobacterial species and cyanobacterial species but fewer than Streptomyces species (Table 2). The highest diversity of P450 families and subfamilies observed in alphaproteobacterial species is possibly due to their diverse lifestyle that led to the generation of diverse P450 families and subfamilies, as the lifestyle of an organism is known to have an impact on the P450 repertoire in the genomes [24]. Thus, for this reason, 77 new P450 families were observed in these species, contributing to the P450 family diversity (Table S3). In aligning with high P450 family and subfamily diversity, alphaproteobacterial species also showed the highest P450 diversity percentage compared to Firmicutes species and Streptomyces species, the same as mycobacterial species but lower than cyanobacterial species (Table 2). A point to be noted is that the number of cyanobacterial species analyzed is almost half of alphaproteobacterial species and thus cyanobacterial species have a slightly higher P450 diversity percentage (Table 2). P450 family-level analysis revealed that some P450 families are expanded in alphaproteobacterial species (Table 3). Among the P450 families, the CYP202 family has the highest number of P450s (70 P450s), followed by CYP153 and CYP173 (each 60 P450s), CYP108 and CYP196 (each 37 P450s) (Table 3). Sixty-five P450 families have a single member and 29 P450 families have two members, indicating the diversity of P450 families in alphaproteobacterial species (Table 3). Analysis of P450 families and subfamilies in alphaproteobacterial species revealed that the P450 family with the highest number of subfamilies was CYP108 (nine subfamilies) followed by CYP101, CYP152, CYP199, CYP173 (each with six subfamilies) and CYP206 (fivesubfamilies). Interestingly, some of the subfamilies are expanded in the dominant P450 families such as CYP202, CYP173, CYP153, CYP196, CYP289, CYP201, CYP1101, CYP112, CYP114 and CYP117, where subfamily “A” is dominant (Table 3). A detailed analysis of the subfamilies and their member count is presented in Table 3. Heat map analysis of the presence and absence of P450 families in different alphaproteobacterial species revealed that none of the P450 family is conserved in these species (Figure 3). Non-conservation of P450 families in bacterial species is not common; it was also observed in cyanobacterial species [41]. However, the co-presence of quite a number of P450 families was found in a large number of species, such as CYP173 and CYP202 in 32 species, CYP196 and CYP201 in 16 species, CYP195 and CYP199 in 10 species, CYP195 and CYP196 in 11 species, CY173, CYP117, CYP127, CYP112 and CYP114 in 10 species and CYP173, CYP147, CYP201 and CYP206 in 9 species (Figure 3). A point to be noted is that the co-presence of some P450 families was also observed in other bacterial species [24]. When compared to other bacterial species, the CYP202 family is dominant in alphaproteobacterial species, while CYP110 is dominant in cyanobacterial species, CYP125 in mycobacterial species, and CYP107 in both the Firmicutes species and Streptomyces species (Table 2).

2.3. Only a Few P450s Are Involved in Secondary Metabolism in Alphaproteobacterial Species

Analysis of metabolite BGCs revealed that 504 species of 599 alphaproteobacterial species have secondary metabolite BGCs in their genomes (Supplementary Dataset 3). In total, 2262 secondary metabolite BGCs belonging to 93 types were found in 504 species (Supplementary Dataset 3). Among the 93 cluster types, terpene was the dominant cluster (270 clusters), followed by Homoserine lactone (hserlactone) (180 clusters), bacteriocin (118 clusters), betalactone (82 clusters) and nonribosomal peptides (NRPS) (79 clusters) in these species (Figure 4). Most similar known cluster analysis revealed that of 2262 BGCs, 642 BGCs showed similarity to 132 known clusters and 96 of these clusters showed 100% identity to known clusters, indicating their involvement in the specific secondary metabolite (Supplementary Dataset 3).
Analysis of P450s that are part of BGCs revealed that only 2.4% of alphaproteobacterial species P450s are involved in the production of secondary metabolites (Table 4). The percentage of P450s involved in the production of secondary metabolites in alphaproteobacterial species was found to be lowest compared to other bacterial species such as Cyanobacteria (8%), Firmicutes species (18%), mycobacterial species (11%) and Streptomyces species (23%) (Table 2). Twenty-one P450s from 19 alphaproteobacterial species were found to be part of BGCs (Table 4). Of the 16 P450 families that are part of BGCs, CYP206 is the dominant family (5 P450s-24%), followed by CYP1101 (2 P450s-10%) and the remaining 14 P450 families (CYP195, CYP1101, CYP2334, CYP199, CYP173, CYP153, CYP152, CYP1302, CYP127, CYP1246, CYP1138, CYP1104, CYP108, CYP107, CYP1326) have a single P450 (Table 4). Dominant P450 families such as CYP173, CYP153 and CYP108 have only one representative as part of BGCs, indicating no correlation between the dominant P450 family vs BGCs (Table 4). A point to be noted is that three P450s, CYP1326A2, CYP2334A1 and CYP195A21 from Sulfitobacter sp. AM1-D1, were found to be part of BGCs (Table 4). Most similar known cluster analysis revealed that CYP1101A27 from Celeribacter indicus was certainly involved in the production of ectoine, as the percentage identity with the most similar known cluster was 100% (Table 4).

2.4. Alphaproteobacterial P450 Families Can Be Found in Other Bacterial Species

It is a well-known fact that bacterial species in Alphaproteobacteria are regarded as an ancient group of bacteria [25,26]. A comparison of P450 families with other bacterial species will provide important P450 family evolutionary distribution patterns. To understand this aspect, we performed comprehensive comparative analysis of P450 families from different bacterial species (Figure 5). As shown in Figure 5, four P450 families were commonly found among alphaproteobacterial species, Firmicutes species and cyanobacterial species. Quite a large number of P450 families were found to be common among alphaproteobacterial species and Streptomyces species (22 P450 families), and mycobacterial species (14 P450 families) (Figure 5). Only the CYP107 family was found to be conserved among all bacterial species (Figure 5). CYP152 and CYP197 were commonly found in alphaproteobacterial species, Firmicutes species and Streptomyces species. Neither of these families is present in mycobacterial species. This indicates that the P450 families that are commonly found among alphaproteobacterial species and other bacterial species are indeed passed from alphaproteobacterial species and retained by other bacterial species throughout the speciation, suggesting the important role of these P450 families. One interesting observation is that the CYP125 P450 family, a cholesterol and cholest-4-en-3-one hydroxylase [47,48] and potential drug target against tuberculosis-causing bacteria Mycobacterium tuberculosis H37Rv [49], is commonly found in Alphaproteobacteria and mycobacterial species, and Streptomyces species (Figure 5), indicating that a cholesterol hydroxylation capability already existed in ancient bacteria as described elsewhere [50] and that these bacterial species retained this ability to survive better in the host environment, as mentioned elsewhere [39].
A point to be noted is that only one alphaproteobacterial species, namely Phenylobacterium zucineum, was found to have CYP125 P450 among 599 species used in the study, indicating two scenarios, i.e., loss of this P450 in other alphaproteobacterial species or gain by P. zucineum. This phenomenon needs further investigation. However, the presence of CYP125 P450 indicates cholesterol side-chain degradation and its subsequent utilization as carbon source via the beta-oxidation pathway [51] in P. zucineum. In order to identify the cholesterol side-chain oxidizing enzymes in P. zucineum, if any, we used M. tuberculosis H37Rv cholesterol side-chain oxidizing enzymes as listed in the literature [51]. Based on the homology percentage and annotation of enzymes at KEGG, we identify all possible cholesterol side-chain oxidizing enzymes (31 enzymes) in P. zucineum, indicating that this alphaproteobacterial species is indeed capable of oxidizing the cholesterol side-chain, where the CYP125 reaction will be critical in removing the side-chain from cholesterol (Table 5). As shown in Table 5, two homologs (one named fadD) were found for fadD36 (Rv1193) acyl-CoA synthetase. Our prediction on cholesterol side-chain degradation by P. zucineum is solely based on homolog proteins’ functions, including CYP125, and further validation by experimentation is required.

2.5. Functional Analysis of Alphaproteobacterial P450s

A literature survey on functional analysis of alphaproteobacterial P450s revealed that several P450s are functionally characterized from these species and are involved in oxidation of xenobiotic compounds (Table 6). These study results revealed that 21 P450s were found to be involved in synthesis, so secondary metabolites were CYP1101A27 from C. indicus, involved in the production of ectoine (Table 4). However, the physiological functional relevance of P450s in these species needs to be investigated.

3. Materials and Methods

3.1. Species and Database

In this study, 599 alphaproteobacterial species’ genomes that are available for public use at Kyoto Encyclopedia of Genes and Genomes (KEGG) database [52] were used (Table S1). Detailed information on species, species codes, genera and GenBank accession codes is presented in Table S1.

3.2. Genome Data Mining and Annotation of P450s

P450 mining in alphaproteobacterial species was carried out using the methods recently described by our laboratory [24,41]. Briefly, the complete proteome of alphaproteobacterial species was downloaded from KEGG and subjected to the NCBI Batch Web CD-Search Tool [68]. Proteins that belong to the P450 superfamily were selected and annotated, as per the International P450 Nomenclature Committee rule, i.e., proteins with >40% identity are grouped under the same family and proteins with >55% identity are grouped under the same subfamily [21,22,44]. Proteins with <40% identity with named P450s are assigned to a new P450 family.

3.3. Phylogenetic Analysis of P450s

Phylogenetic analysis of P450s was carried out following the procedure described recently by our laboratory [24,41]. The phylogenetic tree of P450s was constructed using alphaproteobacterial species P450 protein sequences. Firstly, the MAFFT v6.864 [69] was used to align the protein sequences that are part of the Trex web server [70]. The alignments were then be subjected to interpret the best tree by the Trex web server [70]. Lastly, a web-based tool, iTol, was used to create, visualize and color the tree [71].

3.4. Generation of P450 Profile Heat Maps

P450 profile heat maps were generated following the procedure described recently by our laboratory [24,41]. The heat map was generated using the P450 family data to show the presence or absence of P450s in alphaproteobacterial species. The data were represented as (−3) for family absence (green) and (3) for family presence (red). A tab-delimited file was imported into Mev (Multi-experiment viewer) [72]. Hierarchical clustering using a Euclidean distance metric was used to cluster the data. alphaproteobacterial species with P450s in their genome formed the vertical axis and P450 families formed the horizontal axis.

3.5. Secondary Metabolite BGC Analysis and Identification of P450s That Are Part of BGCs

Secondary metabolite BGC analysis and identification of P450s that are part of BGCs in alphaproteobacterial species was carried out following the procedure described recently by our laboratory [24,41]. Briefly, alphaproteobacterial species’ individual genome ID (Table S1) was submitted to anti-SMASH [46] for identification of secondary metabolite BGCs. Results were downloaded both in the form of Excel spreadsheets representing species-wise cluster information and gene cluster sequences in a Word file. P450s that are part of a specific gene cluster were identified by manually going through the BGCs sequence. Standard gene cluster abbreviation terminology present at the anti-SMASH database [46] was sustained in this study.

3.6. Identification of Cholesterol Side-Chain Oxidizing Genes/Proteins

Thirty-one cholesterol side-chain oxidizing genes/proteins from M. tuberculosis H37Rv [51] were used in the study to identify homologs in P. zucineum. Reference proteins were blasted individually against the P. zucineum genome and based on the percentage identity, homology, coverage and annotation at KEGG, the homolog proteins were identified in P. zucineum. A point to be noted is that M. tuberculosis and P. zucineum are distantly related and thus the lowest percentage identity between these two organisms’ proteins is expected, as previously observed [50].

3.7. Data Analysis

All calculations were carried out following the procedure described previously [24,41,73]. The average number of P450s was calculated using the formula: Average number of P450s = Number of P450s/Number of species. The P450 diversity percentage was calculated using the formula: 100 × Total number of P450 families/(Total number of P450s × Number of species with P450s). The percentage of P450s that formed part of BGCs was calculated using the formula: Percentage of P450s part of BGCs = 100 × Number of P450s part of BGCs/Total number of P450s present in species.

4. Conclusions

Cytochrome P450 monooxygenases (CYPs/P450s) have been well-known proteins for the last six decades. The stereo- and regio-specific oxidation of a variety of compounds by these enzymes led to their applications in quite large areas of biological research. Understanding the evolution of P450s in the bacterial population is now gaining momentum owing to the availability of a large number of bacterial genomes. This study is the first of its kind on the analysis of P450s in an ancient group of bacteria belonging to the class Alphaproteobacteria. Comparative analysis of P450s between different bacterial species revealed that during speciation a large number of P450s became part of secondary metabolite gene clusters (as observed in Streptomyces species and mycobacterial species) and some P450s were passed all the way from Alphaproteobacteria to other bacterial species. This study also reports the first identification of CYP125 P450 in alphaproteobacterial species and predicts that P. zucineum is capable of utilizing the cholesterol side-chain as carbon source. Future study should include thorough profiling of cholesterol-degrading genes/proteins and experimental validation of cholesterol degrading ability. Furthermore, comparative analysis of alphaproteobacterial species and more bacterial species should be carried out to deduce the evolution pattern of P450 families and to identify the loss/gain of new P450 families compared to Alphaproteobacteria with respect to habitat or geographical distribution. The results of this study will serve as reference for future genome data mining and annotation of P450s in species of Alphaproteobacteria.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijms22115542/s1, Table S1: Information on Alphaproteobacterial species and their respective genera used in the study, Table S2: Comparative analysis of P450s and those associated with secondary metabolite biosynthetic gene clusters in the bacterial class Alphaproteobacteria, Table S3: List of new P450 families identified in Alphaproteobacterial species., Dataset 1: P450 sequences identified in Alphaproteobacterial species; Dataset 2: P450 heat-map data; Dataset 3: Secondary metabolite biosynthetic gene clusters (BGCs) analysis and P450s associated with BGCs in Alphaproteobacterial species.

Author Contributions

Conceptualization, K.S.; methodology, N.N., T.P., P.R.S., J.D.K., W.C., D.G., D.R.N. and K.S.; software, N.N., T.P., P.R.S., J.D.K., W.C., D.G., D.R.N. and K.S.; validation, N.N., T.P., P.R.S., J.D.K., W.C., D.G., D.R.N. and K.S.; formal analysis, N.N., T.P., P.R.S., J.D.K., W.C., D.G., D.R.N. and K.S.; investigation, N.N., T.P., P.R.S., J.D.K., W.C., D.G., D.R.N. and K.S.; resources, D.G. and K.S.; data curation, N.N., T.P., P.R.S., J.D.K., W.C., D.G., D.R.N. and K.S.; writing—original draft preparation, N.N., T.P., P.R.S., J.D.K., W.C., D.G., D.R.N. and K.S.; writing—review and editing, N.N., T.P., P.R.S., J.D.K., W.C., D.G., D.R.N. and K.S.; visualization, N.N., T.P., P.R.S., J.D.K., W.C., D.G., D.R.N. and K.S.; supervision, K.S.; project administration, K.S.; funding acquisition, D.G. and K.S. All authors have read and agreed to the published version of the manuscript.

Funding

Khajamohiddin Syed expresses sincere gratitude to the NRF, South Africa for a research grant (Grant Number 114159) and the University of Zululand (Grant number C686). Dominik Gront and Justyna D. Kryś were supported by the Polish National Science Centre (NCN) 2018/29/B/ST6/01989. Master’s students Tiara Padayachee and Nomfundo Nzuza thank the NRF, South Africa for master’s scholarships (Grant numbers MND190619448759 and MND190626451135). Puleng Rosinah Syed also thanks the NRF, South Africa for a doctoral scholarship (Grant number MND190606443406).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors want to thank Barbara Bradley, Pretoria, South Africa for English language editing.

Conflicts of Interest

The authors declare no conflict of interest and the funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Yamazaki, H. Fifty Years of Cytochrome P450 Research; Springer Japan: Tokyo, Japan, 2014; pp. 1–409. [Google Scholar]
  2. Li, Z.; Jiang, Y.; Guengerich, F.P.; Ma, L.; Li, S.; Zhang, W. Engineering cytochrome P450 enzyme systems for biomedical and biotechnological applications. J. Biol. Chem. 2020, 295, 833–849. [Google Scholar] [CrossRef]
  3. Urlacher, V.B.; Girhard, M. Cytochrome P450 monooxygenases in biotechnology and synthetic biology. Trends Biotechnol. 2019, 37, 882–897. [Google Scholar] [CrossRef] [PubMed]
  4. Nelson, D.R. Cytochrome P450 diversity in the tree of life. Biochim. Biophys. Acta (BBA) Proteins Proteom. 2018, 1866, 141–154. [Google Scholar] [CrossRef] [PubMed]
  5. 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] [Green Version]
  6. Guengerich, F.P. A history of the roles of cytochrome P450 enzymes in the toxicity of drugs. Toxicol. Res. 2021, 37, 1–23. [Google Scholar] [CrossRef]
  7. Lepesheva, G.I.; Friggeri, L.; Waterman, M.R. CYP51 as drug targets for fungi and protozoan parasites: Past, present and future. Parasitology 2018, 145, 1820–1836. [Google Scholar] [CrossRef]
  8. Jawallapersand, P.; Mashele, S.S.; Kovačič, L.; Stojan, J.; Komel, R.; Pakala, S.B.; Kraševec, N.; Syed, K. Cytochrome P450 monooxygenase CYP53 family in fungi: Comparative structural and evolutionary analysis and its role as a common alternative anti-fungal drug target. PLoS ONE 2014, 9, e107209. [Google Scholar] [CrossRef] [Green Version]
  9. Girvan, H.M.; Munro, A.W. Applications of microbial cytochrome P450 enzymes in biotechnology and synthetic biology. Curr. Opin. Chem. Biol. 2016, 31, 136–145. [Google Scholar] [CrossRef]
  10. McLean, K.J.; Leys, D.; Munro, A.W. Microbial cytochrome P450s In Cytochrome P450: Structure, Mechanism, and Biochemistry, 4th ed.; Oritz de Montellano, P.R., Ed.; Springer International Publishing: Berlin/Heidelberg, Germany, 2015; pp. 261–407. [Google Scholar]
  11. Kelly, S.L.; Kelly, D.E. Microbial cytochromes P450: Biodiversity and biotechnology. Where do cytochromes P450 come from, what do they do and what can they do for us? Philos. Trans. R. Soc. B Biol. Sci. 2013, 368, 20120476. [Google Scholar] [CrossRef] [Green Version]
  12. Girhard, M.; Bakkes, P.J.; Mahmoud, O.; Urlacher, V.B. P450 Biotechnology. In Cytochrome P450; Springer: Berlin/Heidelberg, Germany, 2015; pp. 451–520. [Google Scholar]
  13. Di Nardo, G.; Gilardi, G. Natural compounds as pharmaceuticals: The key role of cytochromes P450 reactivity. Trends Biochem. Sci. 2020, 45, 511–525. [Google Scholar] [CrossRef]
  14. Andersen, J.F.; Tatsuta, K.; Gunji, H.; Ishiyama, T.; Hutchinson, C.R. Substrate specificity of 6-deoxyerythronolide B hydroxylase, a bacterial cytochrome P450 of erythromycin A biosynthesis. Biochemistry 1993, 32, 1905–1913. [Google Scholar] [CrossRef]
  15. Bischoff, D.; Bister, B.; Bertazzo, M.; Pfeifer, V.; Stegmann, E.; Nicholson, G.J.; Keller, S.; Pelzer, S.; Wohlleben, W.; Süssmuth, R. The biosynthesis of vancomycin-type glycopeptide antibiotics-a model for oxidative side-chain cross-linking by oxygenases coupled to the action of peptide synthetases. ChemBioChem 2005, 6, 267–272. [Google Scholar] [CrossRef]
  16. Jennewein, S.; Park, H.; DeJong, J.M.; Long, R.M.; Bollon, A.P.; Croteau, R.B. Coexpression in yeast of Taxus cytochrome P450 reductase with cytochrome P450 oxygenases involved in Taxol biosynthesis. Biotechnol. Bioeng. 2005, 89, 588–598. [Google Scholar] [CrossRef]
  17. Van Beilen, J.B.; Holtackers, R.; Lüscher, D.; Bauer, U.; Witholt, B.; Duetz, W.A. Biocatalytic production of perillyl alcohol from limonene by using a novel Mycobacterium sp. cytochrome P450 alkane hydroxylase expressed in Pseudomonas putida. Appl. Environ. Microbiol. 2005, 71, 1737–1744. [Google Scholar] [CrossRef] [Green Version]
  18. Podust, L.M.; Sherman, D.H. Diversity of P450 enzymes in the biosynthesis of natural products. Nat. Prod. Rep. 2012, 29, 1251–1266. [Google Scholar] [CrossRef] [Green Version]
  19. Greule, A.; Stok, J.E.; De Voss, J.J.; Cryle, M.J. Unrivalled diversity: The many roles and reactions of bacterial cytochromes P450 in secondary metabolism. Nat. Prod. Rep. 2018, 35, 757–791. [Google Scholar] [CrossRef] [Green Version]
  20. Cimermancic, P.; Medema, M.H.; Claesen, J.; Kurita, K.; Brown, L.C.W.; Mavrommatis, K.; Pati, A.; Godfrey, P.A.; Koehrsen, M.; Clardy, J.; et al. Insights into secondary metabolism from a global analysis of prokaryotic biosynthetic gene clusters. Cell 2014, 158, 412–421. [Google Scholar] [CrossRef] [Green Version]
  21. Nelson, D.R. Cytochrome P450 nomenclature. Methods Mol. Biol. 1998, 107, 15–24. [Google Scholar]
  22. Nelson, D.R. Cytochrome P450 nomenclature, 2004. Methods Mol. Biol. 2006, 320, 1–10. [Google Scholar]
  23. Nelson, D.R. The cytochrome p450 homepage. Hum. Genom. 2009, 4, 59–65. [Google Scholar] [CrossRef] [Green Version]
  24. 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, 1–12. [Google Scholar] [CrossRef] [PubMed]
  25. Muñoz-Gómez, S.; Hess, S.; Burger, G.; Lang, B.F.; Susko, E.; Slamovits, C.H.; Roger, A.J. An updated phylogeny of the Alphaproteobacteria reveals that the parasitic Rickettsiales and Holosporales have independent origins. eLife 2019, 8, 42535. [Google Scholar] [CrossRef] [PubMed]
  26. Ettema, T.J.; Andersson, S.G. The α-proteobacteria: The Darwin finches of the bacterial world. Biol. Lett. 2009, 5, 429–432. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Mukhopadhya, I.; Hansen, R.; El-Omar, E.M.; Hold, G.L. IBD—what role do Proteobacteria play? Nat. Rev. Gastroenterol. Hepatol. 2012, 9, 219–230. [Google Scholar] [CrossRef]
  28. Kersters, K.; De Vos, P.; Gillis, M.; Swings, J.; Vandamme, P.; Stackebrandt, E. Introduction to the Proteobacteria. In The Prokaryotes; Springer Science and Business Media LLC: Berlin/Heidelberg, Germany, 2006; pp. 3–37. [Google Scholar]
  29. Gupta, R.S. Protein signatures distinctive of Alphaproteobacteria and its subgroups and a model for α –Proteobacterial evolution. Crit. Rev. Microbiol. 2005, 31, 101–135. [Google Scholar] [CrossRef]
  30. Phung, N.T.; Lee, J.; Kang, K.H.; Chang, I.S.; Gadd, G.M.; Kim, B.H. Analysis of microbial diversity in oligotrophic microbial fuel cells using 16S rDNA sequences. FEMS Microbiol. Lett. 2004, 233, 77–82. [Google Scholar] [CrossRef] [Green Version]
  31. Timmermans, M.L.; Paudel, Y.P.; Ross, A.C. Investigating the biosynthesis of natural products from marine Proteobacteria: A survey of molecules and strategies. Mar. Drugs 2017, 15, 235. [Google Scholar] [CrossRef] [Green Version]
  32. O’Brien, J.; Wright, G.D. An ecological perspective of microbial secondary metabolism. Curr. Opin. Biotechnol. 2011, 22, 552–558. [Google Scholar] [CrossRef]
  33. Tsukimoto, M.; Nagaoka, M.; Shishido, Y.; Fujimoto, J.; Nishisaka, F.; Matsumoto, S.; Harunari, E.; Imada, C.; Matsuzaki, T. Bacterial production of the tunicate-derived antitumor cyclic depsipeptide didemnin B. J. Nat. Prod. 2011, 74, 2329–2331. [Google Scholar] [CrossRef]
  34. Oh, D.-C.; Strangman, W.K.; Kauffman, C.A.; Jensen, P.R.; Fenical, W. Thalassospiramides A (I) and B (II), Immunosuppressive peptides from the marine bacterium Thalassospira sp. Org. Lett. 2007, 9, 1525–1528. [Google Scholar] [CrossRef]
  35. Meena, K.K.; Kumar, M.; Kalyuzhnaya, M.G.; Yandigeri, M.S.; Singh, D.P.; Saxena, A.K.; Arora, D.K. Epiphytic pink-pigmented methylotrophic bacteria enhance germination and seedling growth of wheat (Triticum aestivum) by producing phytohormone. Antonie van Leeuwenhoek 2011, 101, 777–786. [Google Scholar] [CrossRef]
  36. Siqueira, A.F.; Ormeño-Orrillo, E.; Souza, R.C.; Rodrigues, E.P.; Almeida, L.G.P.; Barcellos, F.G.; Batista, J.S.S.; Nakatani, A.S.; Martínez-Romero, E.; Vasconcelos, A.T.R.; et al. Comparative genomics of Bradyrhizobium japonicum CPAC 15 and Bradyrhizobium diazoefficiens CPAC 7: Elite model strains for understanding symbiotic performance with soybean. BMC Genom. 2014, 15, 1–21. [Google Scholar] [CrossRef] [Green Version]
  37. Khan, A.L.; Waqas, M.; Kang, S.-M.; Al-Harrasi, A.; Hussain, J.; Al-Rawahi, A.; Al-Khiziri, S.; Ullah, I.; Ali, L.; Jung, H.-Y.; et al. Bacterial endophyte Sphingomonas sp. LK11 produces gibberellins and IAA and promotes tomato plant growth. J. Microbiol. 2014, 52, 689–695. [Google Scholar] [CrossRef]
  38. Long, C.; Lu, X.-L.; Gao, Y.; Jiao, B.-H.; Liu, X.-Y. Description of a sulfitobacter strain and its extracellular cyclodipeptides. Evidence Based Complement. Altern. Med. 2011, 2011, 1–6. [Google Scholar]
  39. Senate, L.M.; Tjatji, M.P.; Pillay, K.; Chen, W.; Zondo, N.M.; Syed, P.R.; Mnguni, F.C.; Chiliza, Z.E.; Bamal, H.D.; Karpoormath, R.; et al. Similarities, variations, and evolution of cytochrome P450s in Streptomyces versus Mycobacterium. Sci. Rep. 2019, 9, 1–12. [Google Scholar] [CrossRef] [Green Version]
  40. Mnguni, F.C.; Padayachee, T.; Chen, W.; Gront, D.; Yu, J.-H.; Nelson, D.R.; Syed, K. More P450s are involved in secondary metabolite biosynthesis in Streptomyces compared to Bacillus, Cyanobacteria and Mycobacterium. Int. J. Mol. Sci. 2020, 21, 4814. [Google Scholar] [CrossRef]
  41. Khumalo, M.J.; Nzuza, N.; Padayachee, T.; Chen, W.; Yu, J.-H.; Nelson, D.R.; Syed, K. Comprehensive analyses of cytochrome P450 monooxygenases and secondary metabolite biosynthetic gene clusters in Cyanobacteria. Int. J. Mol. Sci. 2020, 21, 656. [Google Scholar] [CrossRef] [Green Version]
  42. Parvez, M.; Qhanya, L.B.; Mthakathi, N.T.; Kgosiemang, I.K.R.; Bamal, H.D.; Pagadala, N.S.; Xie, T.; Yang, H.; Chen, H.; Theron, C.W.; et al. Molecular evolutionary dynamics of cytochrome P450 monooxygenases across kingdoms: Special focus on mycobacterial P450s. Sci. Rep. 2016, 6, 33099. [Google Scholar] [CrossRef] [Green Version]
  43. 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] [Green Version]
  44. Nelson, D.R.; Kamataki, T.; Waxman, D.J.; Guengerich, F.P.; Estabrook, R.W.; Feyereisen, R.; Gonzalez, F.J.; Coon, M.J.; Gunsalus, I.C.; Gotoh, O.; et al. The P450 Superfamily: Update on new sequences, gene mapping, accession numbers, early trivial names of enzymes, and nomenclature. DNA Cell Biol. 1993, 12, 1–51. [Google Scholar] [CrossRef]
  45. Nelson, D.R. Cytochrome P450 nomenclature. In Cytochrome P450 Protocols; Springer: Berlin/Heidelberg, Germany, 1998; pp. 15–24. [Google Scholar]
  46. Blin, K.; Shaw, S.; Steinke, K.; Villebro, R.; Ziemert, N.; Lee, S.Y.; Medema, M.H.; Weber, T. antiSMASH 5.0: Updates to the secondary metabolite genome mining pipeline. Nucleic Acids Res. 2019, 47, W81–W87. [Google Scholar] [CrossRef] [Green Version]
  47. McLean, K.J.; Lafite, P.; Levy, C.; Cheesman, M.R.; Mast, N.; Pikuleva, I.A.; Leys, D.; Munro, A.W. The structure of Mycobacterium tuberculosis CYP125: Molecular basis for cholesterol binding in a P450 needed for host infection. J Biol Chem 2009, 284, 35524–35533. [Google Scholar] [CrossRef] [Green Version]
  48. Ouellet, H.; Guan, S.; Johnston, J.B.; Chow, E.D.; Kells, P.M.; Burlingame, A.L.; Cox, J.S.; Podust, L.M.; De Montellano, P.R.O. Mycobacterium tuberculosis CYP125A1, a steroid C27 monooxygenase that detoxifies intracellularly generated cholest-4-en-3-one. Mol. Microbiol. 2010, 77, 730–742. [Google Scholar] [CrossRef] [Green Version]
  49. Ortiz de Montellano, P.R. Potential drug targets in the Mycobacterium tuberculosis cytochrome P450 system. J. Inorg. Biochem. 2018, 180, 235–245. [Google Scholar] [CrossRef]
  50. Holert, J.; Cardenas, E.; Bergstrand, L.H.; Zaikova, E.; Hahn, A.S.; Hallam, S.J.; Mohn, W.W. Metagenomes reveal global distribution of bacterial steroid catabolism in natural, engineered, and host environments. mBio 2018, 9, e02345-17. [Google Scholar] [CrossRef] [Green Version]
  51. Van Wyk, R.; Van Wyk, M.; Mashele, S.S.; Nelson, D.R.; Syed, K. Comprehensive comparative analysis of cholesterol catabolic genes/proteins in mycobacterial species. Int. J. Mol. Sci. 2019, 20, 1032. [Google Scholar] [CrossRef] [Green Version]
  52. Kanehisa, M.; Sato, Y.; Furumichi, M.; Morishima, K.; Tanabe, M. New approach for understanding genome variations in KEGG. Nucleic Acids Res. 2019, 47, D590–D595. [Google Scholar] [CrossRef] [Green Version]
  53. Hall, E.A.; Bell, S.G. The efficient and selective biocatalytic oxidation of norisoprenoid and aromatic substrates by CYP101B1 from Novosphingobium aromaticivorans DSM12444. RSC Adv. 2014, 5, 5762–5773. [Google Scholar] [CrossRef] [Green Version]
  54. Ma, M.; Bell, S.G.; Yang, W.; Hao, Y.; Rees, N.H.; Bartlam, M.; Zhou, W.; Wong, L.-L.; Rao, Z. Structural Analysis of CYP101C1 from Novosphingobium aromaticivorans DSM12444. ChemBioChem 2010, 12, 88–99. [Google Scholar] [CrossRef]
  55. Yang, W.; Bell, S.G.; Wang, H.; Zhou, W.; Hoskins, N.; Dale, A.; Bartlam, M.; Wong, L.-L.; Rao, Z. Molecular characterization of a class I P450 electron transfer system from Novosphingobium aromaticivorans DSM12444. J. Biol. Chem. 2010, 285, 27372–27384. [Google Scholar] [CrossRef] [Green Version]
  56. Batabyal, D.; Poulos, T.L. Crystal structures and functional characterization of wild-type CYP101D1 and its active site mutants. Biochem. 2013, 52, 8898–8906. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Yang, W.; Bell, S.G.; Wang, H.; Zhou, W.; Bartlam, M.; Wong, L.-L.; Rao, Z. The structure of CYP101D2 unveils a potential path for substrate entry into the active site. Biochem. J. 2010, 433, 85–93. [Google Scholar] [CrossRef] [PubMed]
  58. Bell, S.G.; Yang, W.; Yorke, J.A.; Zhou, W.; Wang, H.; Harmer, J.; Copley, R.; Zhang, A.; Zhou, R.; Bartlam, M.; et al. Structure and function of CYP108D1 from Novosphingobium aromaticivorans DSM12444: An aromatic hydrocarbon-binding P450 enzyme. Acta Crystallogr. Sect. D Biol. Crystallogr. 2012, 68, 277–291. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Zhou, R.; Huang, C.; Zhang, A.; Bell, S.G.; Zhou, W.; Wong, L.-L. Crystallization and preliminary X-ray analysis of CYP153C1 from Novosphingobium aromaticivorans DSM12444. Acta Crystallogr. Sect. F Struct. Biol. Cryst. Commun. 2011, 67, 964–967. [Google Scholar] [CrossRef]
  60. Krings, U.; Berger, R.G. Terpene bioconversion–how does its future look? Nat. Prod. Commun. 2010, 5, 1934578X1000500927. [Google Scholar] [CrossRef] [Green Version]
  61. Cho, M.-A.; Han, S.; Lim, Y.-R.; Kim, V.; Kim, H.; Kim, D. Streptomyces cytochrome P450 enzymes and their roles in the biosynthesis of macrolide therapeutic agents. Biomol. Ther. 2019, 27, 127–133. [Google Scholar] [CrossRef]
  62. Agematu, H.; Matsumoto, N.; Fujii, Y.; Kabumoto, H.; Doi, S.; Machida, K.; Ishikawa, J.; Arisawa, A. Hydroxylation of testosterone by bacterial cytochromes P450 using the Escherichia coli expression system. Biosci. Biotechnol. Biochem. 2006, 70, 307–311. [Google Scholar] [CrossRef] [Green Version]
  63. Guo, D.; Xu, F.; Bell, S.G.; Pang, X.; Bartlam, M.; Wong, L.-L. Purification, crystallization and preliminary crystallographic analysis of CYP 195A2, a P450 enzyme from Rhodopseudomonas palustris. Protein Pept. Lett. 2008, 15, 423–426. [Google Scholar] [CrossRef]
  64. Furuya, T.; Kino, K. Discovery of 2-naphthoic acid monooxygenases by genome mining and their use as biocatalysts. ChemSusChem 2009, 2, 645–649. [Google Scholar] [CrossRef]
  65. Bell, S.G.; Yang, W.; Tan, A.B.H.; Zhou, R.; Johnson, E.O.D.; Zhang, A.; Zhou, W.; Rao, Z.; Wong, L.-L. The crystal structures of 4-methoxybenzoate bound CYP199A2 and CYP199A4: Structural changes on substrate binding and the identification of an anion binding site. Dalton Trans. 2012, 41, 8703–8714. [Google Scholar] [CrossRef]
  66. Coleman, T.; Wong, S.H.; Podgorski, M.N.; Bruning, J.B.; De Voss, J.J.; Bell, S.G. Cytochrome P450 CYP199A4 from Rhodopseudomonas palustris catalyzes heteroatom dealkylations, sulfoxidation, and amide and cyclic hemiacetal formation. ACS Catal. 2018, 8, 5915–5927. [Google Scholar] [CrossRef] [Green Version]
  67. Berne, C.; Pignol, D.; Lavergne, J.; Garcia, D. CYP201A2, a cytochrome P450 from Rhodopseudomonas palustris, plays a key role in the biodegradation of tributyl phosphate. Appl. Microbiol. Biotechnol. 2007, 77, 135–144. [Google Scholar] [CrossRef]
  68. Lu, S.; Wang, J.; Chitsaz, F.; Derbyshire, M.K.; Geer, R.C.; Gonzales, N.R.; Gwadz, M.; I Hurwitz, D.; Marchler, G.H.; Song, J.S.; et al. CDD/SPARCLE: The conserved domain database in 2020. Nucleic Acids Res. 2020, 48, D265–D268. [Google Scholar] [CrossRef] [Green Version]
  69. Katoh, K.; Kuma, K.; Toh, H.; Miyata, T. MAFFT version 5: Improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005, 33, 511–518. [Google Scholar] [CrossRef]
  70. Boc, A.; Diallo, A.B.; Makarenkov, V. T-REX: A web server for inferring, validating and visualizing phylogenetic trees and networks. Nucleic Acids Res. 2012, 40, W573–W579. [Google Scholar] [CrossRef] [Green Version]
  71. Letunic, I.; Bork, P. Interactive tree of life (iTOL) v3: An online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016, 44, W242–W245. [Google Scholar] [CrossRef]
  72. Howe, E.A.; Sinha, R.; Schlauch, D.; Quackenbush, J. RNA-Seq analysis in MeV. Bioinformatics 2011, 27, 3209–3210. [Google Scholar] [CrossRef] [Green Version]
  73. Msomi, N.K.; Padayachee, R.; Nzuza, N.; Syed, P.R.; Krys, J.D.; Chen, W.; Gront, D.; Nelson, R.D.; Syed, K. In silico analysis of P450s and their role in secondary metabolism in the bacterial class Gammaproteobacteria. Int. J. Mol. Sci. 2021, 26, 1538. [Google Scholar]
Figure 1. Comparative analysis of P450 statistics in the bacterial class Alphaproteobacteria. Detailed information is presented in Table S1.
Figure 1. Comparative analysis of P450 statistics in the bacterial class Alphaproteobacteria. Detailed information is presented in Table S1.
Ijms 22 05542 g001
Figure 2. Phylogenetic analysis of alphaproteobacterial P450s. P450 families that are expanded in these species were highlighted in different colors and indicated in the figure. Alphaproteobacterial P450 protein sequences used to construct the phylogenetic tree are presented in Supplementary Dataset 1.
Figure 2. Phylogenetic analysis of alphaproteobacterial P450s. P450 families that are expanded in these species were highlighted in different colors and indicated in the figure. Alphaproteobacterial P450 protein sequences used to construct the phylogenetic tree are presented in Supplementary Dataset 1.
Ijms 22 05542 g002
Figure 3. Heat map figure representing the presence or absence of cytochrome P450 families in 599 alphaproteobacterial species. The data have been represented as −3 for family absence (green) and 3 for family presence (red). One hundred and forty-three alphaproteobacterial species form the vertical axis and 214 P450 families form the horizontal axis. The respective data used in the generation of this figure are presented in Supplementary Dataset 2.
Figure 3. Heat map figure representing the presence or absence of cytochrome P450 families in 599 alphaproteobacterial species. The data have been represented as −3 for family absence (green) and 3 for family presence (red). One hundred and forty-three alphaproteobacterial species form the vertical axis and 214 P450 families form the horizontal axis. The respective data used in the generation of this figure are presented in Supplementary Dataset 2.
Ijms 22 05542 g003
Figure 4. Comparative analysis of types of secondary metabolite BGCs in alphaproteobacterial species. BGCs that are populated in alphaproteobacterial species were presented in the figure. The number at the top of each bar represents the total number of clusters. Detailed information is presented in Supplementary Dataset 3.
Figure 4. Comparative analysis of types of secondary metabolite BGCs in alphaproteobacterial species. BGCs that are populated in alphaproteobacterial species were presented in the figure. The number at the top of each bar represents the total number of clusters. Detailed information is presented in Supplementary Dataset 3.
Ijms 22 05542 g004
Figure 5. Comparative analysis of P450 families among different bacterial species. The number of P450 families indicated with a number (bold case) and P450 families that are common between alphaproteobacterial species and other bacterial species were listed in the figure.
Figure 5. Comparative analysis of P450 families among different bacterial species. The number of P450 families indicated with a number (bold case) and P450 families that are common between alphaproteobacterial species and other bacterial species were listed in the figure.
Ijms 22 05542 g005
Table 1. Some of the secondary metabolites produced by alphaproteobacterial species and their biological functions.
Table 1. Some of the secondary metabolites produced by alphaproteobacterial species and their biological functions.
Secondary MetaboliteSpecies (Source)Biological FunctionReference
Didemnin BTistrella mobilis YIT 12409Anti-viral and anti-tumor properties[33]
ThalassospiramideThalassospira sp. CNJ-328Immunosuppression[34]
Tropodithietic acidPhaeobacter gallaeciensis and P. porticolaAntibiotic, biocontrol and anti-cancer[31]
PhytohormoneMethylobacterium sp. NC4Improve plant growth[35]
6-gingerol and curcuminBradyrhizobium japonicum CPAC 15 and Bradyrhizobium diazoefficiens CPAC 7Antitumor activity[36]
Gibberellins, salicylic acid, indole-3-acetic acid, zeatin, and abscisic acidSphingomonas speciesPromote plant germination and growth[37]
CyclodipeptidesSulfitobacter speciesBioactive agents[38]
Table 2. Comparative analysis of main characteristics of P450s in different bacterial species.
Table 2. Comparative analysis of main characteristics of P450s in different bacterial species.
CategoryAlphaproteobacterial SpeciesFirmicutes SpeciesStreptomyces SpeciesMycobacterial SpeciesCyanobacterial Species
Total no of species analyzed59997220360114
No of species that have P450s24122920360114
No of P450s87471254601784341
No of families143142537736
No of subfamilies2145369813279
Dominant P450 familyCYP202CYP107CYP107CYP125CYP110
Average no. of P450s4127303
P450 diversity percentage0.070.0080.020.070.09
No of P450s part of BGCs21126123120427
Percentage of P450s part of BGCs2.41823118
Reference(s)This work[24,43][39,40][39,42][41]
Table 3. Comparative analysis of P450 families and subfamilies in alphaproteobacterial species.
Table 3. Comparative analysis of P450 families and subfamilies in alphaproteobacterial species.
P450 FamilyP450 CountPercentage CountSubfamilyP450 Count
CYP101010.11%A1
CYP101710.11%A1
CYP101101.14%B1
C3
D2
E2
Q1
R1
CYP102121.37%A10
K1
AC1
CYP103620.23%B1
C1
CYP10340.46%A4
CYP104131.49%A13
CYP10540.46%BF3
CX1
CYP106120.23%B2
CYP106840.46%A4
CYP107510.11%A1
CYP107610.11%A1
CYP107710.11%A1
CYP107880.92%A7
C1
CYP107101.14%AN8
_1
JE1
CYP108120.23%B2
CYP108210.11%A1
CYP108350.57%A3
C2
CYP108610.11%B1
CYP108820.23%B2
CYP108910.11%A1
CYP108374.23%A2
D6
G18
L1
P2
U1
V1
W1
X5
CYP109020.23%C2
CYP109130.34%A3
CYP109440.46%A3
B1
CYP109860.69%A6
CYP1101293.32%A29
CYP1104101.14%A4
C3
E3
CYP1107161.83%A12
B3
C1
CYP111170.80%A5
B1
C1
CYP111630.34%B3
CYP111810.11%B1
CYP11120.23%A1
B1
CYP112232.63%A23
CYP113510.11%A1
CYP113710.11%A1
CYP113820.23%B2
CYP114110.11%A1
CYP114510.11%A1
CYP114222.52%A22
CYP115510.11%B1
CYP11510.11%A1
CYP11610.11%B1
CYP117070.80%A7
CYP117110.11%A1
CYP117310.11%A1
CYP117420.23%A2
CYP117510.11%A1
CYP117222.52%A22
CYP118110.11%A1
CYP118210.11%A1
CYP118620.23%A2
CYP118710.11%A1
CYP120420.23%A2
CYP120610.11%B1
CYP122120.23%B2
CYP122940.46%A1
B2
C1
CYP124620.23%A2
CYP125010.11%A1
CYP125810.11%A1
CYP12510.11%P1
CYP127510.11%B1
CYP127192.17%A18
C1
CYP128110.11%A1
CYP129150.57%B3
C2
CYP130220.23%A2
CYP131120.23%A2
CYP131210.11%A1
CYP132620.23%B2
CYP133010.11%A1
CYP133710.11%A1
CYP13310.11%F1
CYP134910.11%A1
CYP135010.11%A1
CYP137110.11%B1
CYP137620.23%B1
C1
CYP138410.11%A1
CYP139610.11%A1
CYP140510.11%B1
CYP140630.34%B2
C1
CYP146010.11%C1
CYP147161.83%D16
CYP151520.23%A2
CYP15120.23%C2
CYP152131.49%B4
C3
E2
P1
AA1
AB2
CYP153606.86%A44
C2
D14
CYP15540.46%J3
K1
CYP159110.11%A1
CYP159710.11%A1
CYP173230.34%A3
CYP173320.23%A1
B1
CYP173410.11%A1
CYP173510.11%A1
CYP173610.11%A1
CYP173720.23%A2
CYP173810.11%A1
CYP173920.23%A2
CYP173606.86%A45
B8
C3
G1
H1
J2
CYP174010.11%A1
CYP174110.11%A1
CYP174210.11%A1
CYP174310.11%A1
CYP174410.11%A1
CYP174510.11%A1
CYP174610.11%A1
CYP174710.11%A1
CYP174810.11%A1
CYP174920.23%A1
B1
CYP175010.11%A1
CYP175110.11%A1
CYP17524 A4
CYP175310.11%A1
CYP175410.11%A1
CYP175510.11%A1
CYP18010.11%D1
CYP18670.80%K7
CYP19110.11%B1
CYP19230.34%A3
CYP193171.95%A17
CYP194101.14%A10
CYP195212.40%A18
D1
E2
CYP196374.23%A34
B2
C1
CYP19720.23%R2
CYP199161.83%A11
B1
J1
K1
L1
M1
CYP20080.92%A5
B2
C1
CYP201313.55%A25
B3
C3
CYP202708.01%A46
B24
CYP203131.49%A12
B1
CYP20420.23%A2
CYP206192.17%A9
B1
C3
D6
CYP21920.23%A2
CYP22330.34%A2
E1
CYP22410.11%A1
CYP22570.80%A7
CYP22620.23%C1
D1
CYP27820.23%A1
C1
CYP28610.11%C1
CYP28820.23%B2
CYP289333.78%A33
CYP29060.69%A5
B1
CYP214010.11%A1
CYP123610.11%A1
CYP233420.23%A2
Table 4. Identification of P450s that are part of the secondary metabolite BGCs in alphaproteobacterial species. Reference cluster information was obtained by performing BLAST at the anti-SMASH database [46], as indicated in the materials and methods section. The cluster type and most similar known cluster names available at the anti-SMASH database [46] were listed in the table.
Table 4. Identification of P450s that are part of the secondary metabolite BGCs in alphaproteobacterial species. Reference cluster information was obtained by performing BLAST at the anti-SMASH database [46], as indicated in the materials and methods section. The cluster type and most similar known cluster names available at the anti-SMASH database [46] were listed in the table.
Species NameP450Reference Cluster Information
Cluster TypeMost Similar Known ClusterSimilarity
Novosphingobium aromaticivoransCYP153C1TerpeneAstaxanthin dideoxyglycoside75%
Mesorhizobium japonicum MAFF 303099CYP127A3v1Hserlactone--
Agrobacterium fabrumCYP206A1Terpene--
Agrobacterium sp. H13-3CYP206A4Terpene--
Agrobacterium tumefaciensCYP206A4Terpene--
Agrobacterium rhizogenesCYP206A2Terpene--
Rhizobium sp. IRBG74CYP206A3Terpene--
Rhizobium sp. NT-26CYP107JE1NAGGN--
Bradyrhizobium sp. S23321CYP199A26Hserlactone, t2pksColabomycin4%
Bradyrhizobium oligotrophicumCYP108L2NRPS--
Bosea vaviloviaeCYP1101A30T1PKSS56-p1, NRPS3%
Beijerinckia indicaCYP173J1NRPS, T1PKS--
Pseudorhodoplanes sinuspersiciCYP1104E2T1PKSSphingan polysaccharide, saccharide13%
Rhodobacter sphaeroides ATCC 17025CYP152C2hserlactoneConglobatin, NRPS10%
Celeribacter indicusCYP1101A27EctoineEctoine, other100%
Hyphomonadaceae bacterium UKL13-1CYP1246A4Terpene--
Sphingopyxis macrogoltabida 203CYP1302A1NRPS--
Xanthobacter autotrophicusCYP1138B1T1PKS--
Sulfitobacter sp.AM1-D1CYP1326A2acyl_amino_acids--
CYP195A21bacteriocin--
CYP2334A1NRPS-like--
Table 5. Identification of homolog proteins involved in cholesterol side-chain oxidation in Phenylobacterium zucineum. Mycobacterium tuberculosis H37Rv cholesterol side-chain degrading proteins from the published literature [51] were used to identify homolog proteins in P. zucineum. Gene and protein IDs and enzyme annotations were from KEGG [52].
Table 5. Identification of homolog proteins involved in cholesterol side-chain oxidation in Phenylobacterium zucineum. Mycobacterium tuberculosis H37Rv cholesterol side-chain degrading proteins from the published literature [51] were used to identify homolog proteins in P. zucineum. Gene and protein IDs and enzyme annotations were from KEGG [52].
Proteins Involved in Cholesterol Side-Chain Oxidation in Mycobacterium tuberculosis H37RvHomolog Proteins in Phenylobacterium zucineum
Gene NameGene IDEnzymeProtein ID% Identity% HomologyEnzyme
fadD36Rv1193Acyl-CoA synthetasePHZ_c1345 (fadD)2941Long-chain acyl-CoA synthetase
pzu:PHZ_c11553852Malonyl-CoA/methylmalonyl-CoA synthetase
fadD19Rv3515cProbable fatty-acid-CoA ligasepzu:PHZ_c20653957Fatty-acid-CoA ligase
fadD3Rv3561Acyl-CoA synthetase (AMP forming)pzu:PHZ_c19093547Long-chain-fatty-acid—CoA ligase
fadD17Rv3506Possible fatty-acid-CoA ligasepzu:PHZ_c05973046Fatty-acyl-CoA synthase
fadD19Rv3515cProbable fatty-acid-CoA ligasepzu:PHZ_c20653957Long-chain fatty acid:CoA ligase
fadD10Rv0099Fatty acid-CoA synthasepzu:PHZ_c25122943Long-chain-fatty-acid--CoA ligase
fadD9Rv2590Fatty acid-CoA synthasepzu:PHZ_c21232538Long-chain acyl-CoA synthetase
fadD18Rv3513cPossible fatty-acid-CoA ligasepzu:PHZ_c20655573CoA-synthetase, long-chain fatty acid:CoA ligase
fadE22Rv3061cAcyl-CoA dehydrogenasepzu:PHZ_c26784054Dehydrogenase family protein
fadE23Rv3140Acyl-CoA dehydrogenasepzu:PHZ_c08903351Acyl-CoA dehydrogenase
fadE24Rv3139Acyl-CoA dehydrogenasepzu:PHZ_c23653248Acyl-CoA dehydrogenase
fadE25Rv3274cAcyl-CoA dehydrogenase FADE25pzu:PHZ_c16804060Acyl-CoA dehydrogenase
fadE26Rv3504Probable acyl-CoA dehydrogenasepzu:PHZ_c23363350Acyl-CoA dehydrogenase
fadE27Rv3505Probable acyl-CoA dehydrogenasepzu:PHZ_c24062845Isovaleryl CoA dehydrogenase
fadE28Rv3544cProbable acyl-CoA dehydrogenasepzu:PHZ_c24982643Acyl-CoA dehydrogenase
fadE29Rv3543cProbable acyl-CoA dehydrogenasepzu:PHZ_c25413251Acyl-CoA dehydrogenase
fadE30Rv3560cProbable acyl-CoA dehydrogenasepzu:PHZ_c25414157Acyl-CoA dehydrogenase
fadE31Rv3562Probable acyl-CoA dehydrogenasepzu:PHZ_c26783653Dehydrogenase family protein
fadE32Rv3563Probable acyl-CoA dehydrogenasepzu:PHZ_c24983245Acyl-CoA dehydrogenase
fadE33Rv3564Probable acyl-CoA dehydrogenasepzu:PHZ_c26793348Acyl-CoA dehydrogenase family protein
fadE34Rv3573cProbable acyl-CoA dehydrogenasepzu:PHZ_c26783856Dehydrogenase family protein
fadE5Rv0244cAcyl-CoA dehydrogenasepzu:PHZ_c33543651Acyl-CoA dehydrogenase
mbtN (fadE14)Rv1346Acyl-CoA dehydrogenasepzu:PHZ_c03883048Acyl-CoA dehydrogenase
echA9Rv1071c3-Hydroxyisobutyryl-CoA hydrolasepzu:PHZ_c16793957Enoyl-CoA hydratase/isomerase family protein
echA19Rv3516Possible enoyl-CoA hydratasepzu:PHZ_c25353854Enoyl-CoA hydratase/carnithine racemase
echA20Rv3550Possible enoyl-CoA hydratasepzu:PHZ_c34293351Enoyl-CoA hydratase/isomerase family protein
fadB2Rv0468hydroxybutyryl-CoA dehydrogenasepzu:PHZ_c315243603-Hydroxyacyl-CoA dehydrogenase
fadB3Rv1715Hydroxybutyryl-CoA dehydrogenasepzu:PHZ_c315235493-Hydroxyacyl-CoA dehydrogenase
fadA5Rv3546Acetoacetyl-CoA thiolasepzu:PHZ_c25044060Acetyl-CoA C-acetyltransferase
hsd4ARv3502c17β-Hydroxysteroid dehydrogenase (17β-HSD)pzu:PHZ_c20083649Short-chain dehydrogenase/reductase SDR
ltp2Rv3540cProbable ketoacyl-CoA thiolasepzu:PHZ_c32453649Acetyl-CoA C-acetyltransferase
Table 6. Functional analysis of alphaproteobacterial P450s.
Table 6. Functional analysis of alphaproteobacterial P450s.
P450FunctionReference
CYP101B1α-ionone, β-damascone, phenylcyclohexane and p-cymene hydroxylation[53]
CYP101C1Ionone derivative hydroxylation[54]
CYP101D1Terpenoid (camphor) hydroxylase[55,56]
CYP101D2Camphor 5-exo hydroxylase[57]
CYP108D1Aromatic hydrocarbon hydroxylase[55,58]
CYP153C1, CYP153D1Alkane hydroxylase[59]
CYP111A2Oxidizes linalool to 8-hydroxylinalool[55]
CYP111B1β-Ionone hydroxylation[60]
CYP112A2Oxidation of Rapamycin[61]
CYP127A3Hydroxylation of testosterone[62]
CYP195A2Degradation of 4-fluoro-, 4-chloro- and 4-methylsalicylic acid, and 3-chloro- and 3-methylsalicylic acid[63]
CYP199A1Hydroxylation of 2-naphthoic acid[64]
CYP199A2Hydroxylation of 2-naphthoic acid, para-substituted benzoic acids and involved in the degradation of ligninolic compounds[64,65]
CYP199A4Catalyzes heteroatom dealkylations, sulfoxidation, and amide and cyclic hemiacetal formation[66]
CYP200A1Hydroxylation of testosterone[62]
CYP201A2Involved in the biodegradation of tributyl phosphate[67]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nzuza, N.; Padayachee, T.; Syed, P.R.; Kryś, J.D.; Chen, W.; Gront, D.; Nelson, D.R.; Syed, K. Ancient Bacterial Class Alphaproteobacteria Cytochrome P450 Monooxygenases Can Be Found in Other Bacterial Species. Int. J. Mol. Sci. 2021, 22, 5542. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22115542

AMA Style

Nzuza N, Padayachee T, Syed PR, Kryś JD, Chen W, Gront D, Nelson DR, Syed K. Ancient Bacterial Class Alphaproteobacteria Cytochrome P450 Monooxygenases Can Be Found in Other Bacterial Species. International Journal of Molecular Sciences. 2021; 22(11):5542. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22115542

Chicago/Turabian Style

Nzuza, Nomfundo, Tiara Padayachee, Puleng Rosinah Syed, Justyna Dorota Kryś, Wanping Chen, Dominik Gront, David R. Nelson, and Khajamohiddin Syed. 2021. "Ancient Bacterial Class Alphaproteobacteria Cytochrome P450 Monooxygenases Can Be Found in Other Bacterial Species" International Journal of Molecular Sciences 22, no. 11: 5542. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms22115542

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