Cytochrome P450 monooxygenases (P450s) found in all domains of life are known for their catalytic versatility and stereo- and regio-specific activity. While the impact of lifestyle on P450 evolution was reported in many eukaryotes, this remains to be addressed in bacteria. In this report, Streptomyces and Mycobacterium , belonging to the phylum Actinobacteria , were studied owing to their contrasting lifestyles and impacts on human. Analyses of all P450s and those predicted to be associated with secondary metabolism have revealed that different lifestyles have affected the evolution of P450s in these bacterial genera. We have found that while species in both genera have essentially the same number of P450s in the genome, Streptomyces P450s are much more diverse than those of Mycobacterium . Moreover, despite both belonging to Actinobacteria , only 21 P450 families were common, and 123 and 56 families were found to be unique to Streptomyces and Mycobacterium , respectively. The presence of a large and diverse number of P450s in Streptomyces secondary metabolism contributes to antibiotic diversity, helping to secure the niche. Conversely, based on the currently available functional data, types of secondary metabolic pathways and associated P450s, mycobacterial P450s seem to play a role in utilization or synthesis of lipids.

Unraveling the role of cytochrome P450 monooxygenases (CYPs/P450s), heme-thiolate proteins present in living and non-living entities, in secondary metabolite synthesis is gaining momentum. In this direction, in this study, we analyzed the genomes of 203 Streptomyces species for P450s and unraveled their association with secondary metabolism. Our analyses revealed the presence of 5460 P450s, grouped into 253 families and 698 subfamilies. The CYP107 family was found to be conserved and highly populated in Streptomyces and Bacillus species, indicating its key role in the synthesis of secondary metabolites. Streptomyces species had a higher number of P450s than Bacillus and cyanobacterial species. The average number of secondary metabolite biosynthetic gene clusters (BGCs) and the number of P450s located in BGCs were higher in Streptomyces species than in Bacillus, mycobacterial, and cyanobacterial species, corroborating the superior capacity of Streptomyces species for generating diverse secondary metabolites. Functional analysis via data mining confirmed that many Streptomyces P450s are involved in the biosynthesis of secondary metabolites. This study was the first of its kind to conduct a comparative analysis of P450s in such a large number (203) of Streptomyces species, revealing the P450s’ association with secondary metabolite synthesis in Streptomyces species. Future studies should include the selection of Streptomyces species with a higher number of P450s and BGCs and explore the biotechnological value of secondary metabolites they produce.

Cytochrome P450 monooxygenases (CYPs/P450s) are among the most catalytically-diverse enzymes, capable of performing enzymatic reactions with chemo-, regio-, and stereo-selectivity. Our understanding of P450s’ role in secondary metabolite biosynthesis is becoming broader. Among bacteria, Bacillus species are known to produce secondary metabolites, and recent studies have revealed the presence of secondary metabolite biosynthetic gene clusters (BGCs) in these species. However, a comprehensive comparative analysis of P450s and P450s involved in the synthesis of secondary metabolites in Bacillus species has not been reported. This study intends to address these two research gaps. In silico analysis of P450s in 128 Bacillus species revealed the presence of 507 P450s that can be grouped into 13 P450 families and 28 subfamilies. No P450 family was found to be conserved in Bacillus species. Bacillus species were found to have lower numbers of P450s, P450 families and subfamilies, and a lower P450 diversity percentage compared to mycobacterial species. This study revealed that a large number of P450s (112 P450s) are part of different secondary metabolite BGCs, and also identified an association between a specific P450 family and secondary metabolite BGCs in Bacillus species. This study opened new vistas for further characterization of secondary metabolite BGCs, especially P450s in Bacillus species.

Tremellomycetes, a fungal class in the subphylum Agaricomycotina, contain well-known opportunistic and emerging human pathogens. The azole drug fluconazole, used in the treatment of diseases caused by some species of Tremellomycetes, inhibits cytochrome P450 monooxygenase CYP51, an enzyme that converts lanosterol into an essential component of the fungal cell membrane ergosterol. Studies indicate that mutations and over-expression of CYP51 in species of Tremellomycetes are one of the reasons for fluconazole resistance. Moreover, the novel drug, VT-1129, that is in the pipeline is reported to exert its effect by binding and inhibiting CYP51. Despite the importance of CYPs, the CYP repertoire in species of Tremellomycetes has not been reported to date. This study intends to address this research gap. Comprehensive genome-wide CYP analysis revealed the presence of 203 CYPs (excluding 16 pseudo-CYPs) in 23 species of Tremellomycetes that can be grouped into 38 CYP families and 72 CYP subfamilies. Twenty-three CYP families are new and three CYP families (CYP5139, CYP51 and CYP61) were conserved across 23 species of Tremellomycetes. Pathogenic cryptococcal species have 50% fewer CYP genes than non-pathogenic species. The results of this study will serve as reference for future annotation and characterization of CYPs in species of Tremellomycetes.

Tuberculosis (TB) is one of the top infectious diseases causing numerous human deaths in the world. Despite enormous efforts, the physiology of the causative agent, Mycobacterium tuberculosis, is poorly understood. To contribute to better understanding the physiological capacity of these microbes, we have carried out extensive in silico analyses of the 1111 mycobacterial species genomes focusing on revealing the role of the orphan cytochrome P450 monooxygenase (CYP) CYP139 family. We have found that CYP139 members are present in 894 species belonging to three mycobacterial groups: M. tuberculosis complex (850-species), Mycobacterium avium complex (34-species), and non-tuberculosis mycobacteria (10-species), with all CYP139 members belonging to the subfamily “A”. CYP139 members have unique amino acid patterns at the CXG motif. Amino acid conservation analysis placed this family in the 8th among CYP families belonging to different biological domains and kingdoms. Biosynthetic gene cluster analyses have revealed that 92% of CYP139As might be associated with producing different secondary metabolites. Such enhanced secondary metabolic potentials with the involvement of CYP139A members might have provided mycobacterial species with advantageous traits in diverse niches competing with other microbial or viral agents, and might help these microbes infect hosts by interfering with the hosts’ metabolism and immune system.

Phylogenetic and structural analysis of P450 proteins fused to peroxidase/dioxygenase has not been reported yet. We present phylogenetic and in silico structural analysis of the novel P450 fusion family CYP5619 from the deadliest fish pathogenic oomycete, Saprolegnia diclina. Data-mining and annotation of CYP5619 members revealed their unique presence in oomycetes. CYP5619 members have the highest number of conserved amino acids among eukaryotic P450s. The highest number of conserved amino acids (78%) occurred in the peroxidase/dioxygenase domain compared to the P450 domain (22%). In silico structural analysis using a high-quality CYP5619A1 model revealed that CYP5619A1 has characteristic P450 structural motifs including EXXR and CXG. However, the heme-binding domain (CXG) in CYP5619 members was found to be highly degenerated. The in silico substrate binding pattern revealed that CYP5619A1 have a high affinity to medium chain fatty acids. Interestingly, the controlling agent of S. diclina malachite green was predicted to have the highest binding affinity, along with linoleic acid. However, unlike fatty acids, none of the active site amino acids formed hydrogen bonds with malachite green. The study’s results will pave the way for assessing CYP5619A1’s role in S. diclina physiology, including the nature of malachite green binding.

Cytochrome P450 monooxygenases (CYPs/P450s) and their redox partners, ferredoxins, are ubiquitous in organisms. P450s have been studied in biology for over six decades owing to their distinct catalytic activities, including their role in drug metabolism. Ferredoxins are ancient proteins involved in oxidation-reduction reactions, such as transferring electrons to P450s. The evolution and diversification of P450s in various organisms have received little attention and no information is available for archaea. This study is aimed at addressing this research gap. Genome-wide analysis revealed 1204 P450s belonging to 34 P450 families and 112 P450 subfamilies, where some families and subfamilies are expanded in archaea. We also identified 353 ferredoxins belonging to the four types 2Fe-2S, 3Fe-4S, 7Fe-4S and 2[4Fe-4S] in 40 archaeal species. We found that bacteria and archaea shared the CYP109, CYP147 and CYP197 families, as well as several ferredoxin subtypes, and that these genes are co-present on archaeal plasmids and chromosomes, implying the plasmid-mediated lateral transfer of these genes from bacteria to archaea. The absence of ferredoxins and ferredoxin reductases in the P450 operons suggests that the lateral transfer of these genes is independent. We present different scenarios for the evolution and diversification of P450s and ferredoxins in archaea. Based on the phylogenetic analysis and high affinity to diverged P450s, we propose that archaeal P450s could have diverged from CYP109, CYP147 and CYP197. Based on this study’s results, we propose that all archaeal P450s are bacterial in origin and that the original archaea had no P450s.

Cytochrome P450 monooxygenases (CYPs/P450s) are well known for their role in organisms’ primary and secondary metabolism. Among 20 P450s of the tuberculosis-causing Mycobacterium tuberculosis H37Rv, CYP128A1 is particularly important owing to its involvement in synthesizing electron transport molecules such as menaquinone-9 (MK9). This study employs different in silico approaches to understand CYP128 P450 family’s distribution and structural aspects. Genome data-mining of 4250 mycobacterial species has revealed the presence of 2674 CYP128 P450s in 2646 mycobacterial species belonging to six different categories. Contrast features were observed in the CYP128 gene distribution, subfamily patterns, and characteristics of the secondary metabolite biosynthetic gene cluster (BGCs) between M. tuberculosis complex (MTBC) and other mycobacterial category species. In all MTBC species (except one) CYP128 P450s belong to subfamily A, whereas subfamily B is predominant in another four mycobacterial category species. Of CYP128 P450s, 78% was a part of BGCs with CYP124A1, or together with CYP124A1 and CYP121A1. The CYP128 family ranked fifth in the conservation ranking. Unique amino acid patterns are present at the EXXR and CXG motifs. Molecular dynamic simulation studies indicate that the CYP128A1 bind to MK9 with the highest affinity compared to the azole drugs analyzed. This study provides comprehensive comparative analysis and structural insights of CYP128A1 in M. tuberculosis.