Abstract
Introduction
Numerous sectors, including textiles, paint, ink, plastics, and cosmetics, employ various dyes. During their production and usage, some of these dyes are lost, often resulting in environmental concerns. 1 In recent times, many nations have implemented stricter regulations regarding the discharge of colored wastewater that contains dyes. Azo dyes are the most widely manufactured chemically synthesized dyestuffs. 2 These dyes, a class of synthetic compounds primarily used in commercial applications, are aromatic molecules containing one or more azo (N=N) groups. 3 An estimated 1.5 million tons of dyes were produced globally in 2020, with azo dyes comprising more than half of this amount. Azo dyes represent the most significant category of artificial colorants. Azo compounds are extensively utilized across various industries, including paper printing, food, textiles, and cosmetics. 4
Every day, millions of liters of wastewater from textile mills are discharged into public drains, which eventually empty into rivers as untreated effluents. A large portion of this wastewater, particularly that containing azo dyes, demonstrates resistance to degradation. 5 Due to their xenobiotic nature and structural stability, reactive azo dyes are resistant and cannot be fully degraded by conventional wastewater treatment methods, which include chemical treatment, photolysis, and activated sludge. Consequently, these dyes are released into the environment as colored wastewater. The toxicity of these dyes poses serious risks, potentially leading to acute effects on exposed organisms. 6
Numerous approaches have been developed for the reduction of azo dyes. These include physical and chemical methods such as chemical flocculation, activated carbon adsorption, filtration, and specific coagulation techniques. Advanced processes such as reverse osmosis, nanofiltration, and multiple-effect evaporators (MEEs) have also been proven to be effective, although they are often expensive. 7 Moreover, there have also been multiple studies on the application of ionizing radiations to decolorize and degrade the water-contaminating dyes. In research, electron beam radiation has been used to decolorize and degrade various dyes, crystal violet being the representative cationic triphenylmethane dye. The influence of chelating agents, as well as inorganic anions, was explored on the dyes’ degradation, along with total organic carbon elimination. The degradation pathway included primarily N-demethylation, triphenylmethane chromophore cleavage, aromatic product ring breakage and oxidation to carboxylic acid, and mineralization to CO2 and H2O. 8
Biotreatment represents a more cost-effective and environmentally sustainable alternative for azo dye degradation.
9
Various microorganisms, including yeasts, fungi, and bacteria, have been reported to degrade azo dyes. However, despite the structural diversity of dye molecules, only a limited number of enzymes are capable of breaking them down effectively.
10
Among the most promising enzymes for the enzymatic remediation of azo dyes are laccases, azo reductases, and tyrosinase.
11
In this study, we used computational techniques to evaluate the degradation potential of
Recent studies support the role of these enzymes in dye degradation. Therefore, the primary objective of this study was to compare the biodegradation capability of laccase, peroxidase, and tyrosinase in the breakdown of azo dyes using computational approaches.
Methodology
Selection of Enzymes for Comparative Analysis
After identifying the fungal strain, its azo dye-degrading capacity was assessed. Three enzymes—laccase, peroxidase MSP1, and tyrosinase—from
Physicochemical Property Forecasting
The physicochemical properties of the enzymes were analyzed using the ExPASy ProtParam, which calculates parameters such as molecular weight, theoretical pI, amino acid composition, half-life, instability index, aliphatic index, and GRAVY. 15
Homology Modeling of Proteins
Homology modeling is a computational approach used to construct reliable three-dimensional structures of a query protein using the known structures of template proteins.
16
The 3D structures of the three enzymes from
Validation of 3D Structures of Proteins
Model quality was evaluated using PROCHECK and ERRAT via the SAVES server. Ramachandran plots were used to assess residue distribution, and ERRAT was used to identify potential error-prone regions. 17
Selection and Toxicological Analysis of Azo Dyes
Fourteen azo dye structures were selected from the PubChem database, along with their molecular formulas and PubChem identification numbers. Toxicity was predicted using CompTox Chemicals Dashboard, including carcinogenicity, genotoxicity, and the potential for skin or eye irritation. 18
Binding Site Prediction of Enzymes
The active sites of the enzymes were predicted using Discovery Studio 19 to determine potential dye-binding regions and characterize ligand–residue interactions.
Comparative Molecular Docking Analysis
To evaluate the molecular interactions between the selected azo dyes and the three enzymes, docking analysis was performed using PyRx. 20 Dye structures obtained in SDF format were converted to PDB, ligand energies were minimized, and docked to active sites were defined by grid parameters.
Molecular Interaction Analysis of the Docked Complexes
Discovery Studio Visualizer was used to analyze the top protein–ligand complexes based on binding energies. This tool identified the types of bonds and interactions formed between amino acid residues and ligands. 21
Protein–Protein Interaction Analysis
The STRING (Search Tool for the Retrieval of Interacting Genes/Proteins) database was used to perform protein–protein interaction analysis. This analysis was used to identify functional associations of the target enzymes, generating high-confidence interaction networks. 21
Co-expression of Genes
Co-expression analysis of genes encoding laccase, tyrosinase, and dye-degrading peroxidase MSP1 was also performed using the STRING database. The co-expression data for the selected enzymes were analyzed to identify gene clusters potentially involved in azo dye degradation. 23
Flow diagram of the complete process
Results
Selection of Enzymes for Comparative Analysis
The selected enzyme sequences were obtained from the National Center for Biotechnology Information (NCBI) database. The three enzymes chosen for this study were laccase, peroxidase, and tyrosinase. The retrieval process involved querying the NCBI GenBank and Protein databases for sequences corresponding to these enzymes from
Physicochemical Properties of Enzymes
The physicochemical properties of laccase, peroxidase, and tyrosinase were determined via the Expasy ProtParam tool. Their theoretical isoelectric points (pIs) were discovered to be 4.47, 6.84, and 5.83, respectively. Their molecular weights were reported to be 56214.68 Da, 65919.84 Da, and 67691.00 Da, respectively. Subsequently, the stability index was computed and discovered to be less than 40, signifying the enzymes’ stability. When the GRAVY (Grand Average of Hydropathicity) was calculated, a negative value was found, showing that the enzymes are hydrophilic. Other parameters were also predicted, given below in Tables 1–3, respectively, for laccase, peroxidase, and tyrosinase.
Predicted physicochemical properties of laccase from
Predicted physicochemical properties of dye-decolorizing peroxidase MSP1 from
Predicted physicochemical properties of tyrosinase from
Homology Modeling of Proteins
Homology modeling of laccase, peroxidase MSP1, and tyrosinase was performed through SWISS-Model-Expasy, and the model with a high percent identity with the template of each enzyme was downloaded for further analysis. Figures 1–3 show the three-dimensional structure of each enzyme.

Predicted 3D structure of laccase from

Predicted 3D structure of dye-decolorizing peroxidase MSP1 from

Predicted 3D structure of tyrosinase from
Validation of 3D Structures of Proteins
Ramachandran plot results predicted that laccase, peroxidase MSP1, and tyrosinase had 91.7%, 89%, and 88.8% amino acid residues in the favored region. Moreover, the ERRAT graph showed laccase, peroxidase MSP1, and tyrosinase had 91.2%, 88.5%, and 95.5% residues in the non-error region. All these results supported the accuracy of computationally built protein structures. Figures 4 –6 present the results of the RAMACHANDRAN plot for each enzyme, respectively.

Ramachandran plot of the modeled laccase enzyme, used for structural validation, showing the distribution of amino acids in favored, additional allowed, and disallowed regions.

Ramachandran plot of the modeled peroxidase MSP1 enzyme, used for structural validation, showing the distribution of amino acids in favored, additional allowed, and disallowed regions.

Ramachandran plot of the modeled tyrosinase enzyme, used for structural validation, showing the distribution of amino acids in favored, additional allowed, and disallowed regions.
Retrieval of Pollutant Structure and Toxicology Analysis
The PubChem database was used to obtain 14 3D structures of azo dyes as SDF files. Table 4 provides the names, structures, and molecular formulas, along with the PubChem ID of target pollutants and the hazard type of the specific pollutant.
Details of the pollutants of interest, including their chemical names, 3D structures, molecular weights, corresponding PubChem IDs, and toxicological assessment results.
The toxicological assessment of 14 selected azo dyes using the CompTox tool revealed that 12 azo-based compounds were identified as carcinogenic, meaning they have the potential to cause cancer by promoting abnormal cell growth. Additionally, seven were found to be genotoxic, indicating their capacity to damage genetic material such as DNA, which can lead to mutations or other genetic disorders. These findings suggest that a significant proportion of the azo dyes exhibit hazardous properties, underscoring the need for further evaluation and careful consideration of their safety, particularly if they are intended for therapeutic or industrial use. The hazard type of each compound is mentioned in the following table.
Binding Site Prediction of Enzymes
The binding site prediction of an enzyme is crucial in molecular docking analysis because it enables the site-specific docking among proteins and ligands. Binding site pockets of laccase, peroxidase MSP1, and tyrosinase of

Active sites predicted by Discovery Studio for (1) laccase, (2) peroxidase MSP1, and (3) tyrosinase.
Binding site predictions for each enzyme (laccase, peroxidase MSP1, and tyrosinase) using Discovery Studio, with their respective XYZ dimensions and point counts.
Molecular docking analysis
Molecular docking was performed between 14 selected azo-based compounds and each enzyme on PyRx. The grid box was adjusted on the proteins according to the active sites for better docking poses. The binding affinities were shown to be above −6 kcal/mol for all 14 ligands, predicting the good catalyzing ability of enzymes on azo dyes. This indicates that either one enzyme or a combination of all can be beneficial in removing textile dyes from the environment. The values of binding affinities (kcal/mol) of the protein–ligand complexes are shown in Table 6.
Binding affinities (in kcal/mol) of laccase, peroxidase MSP1, and tyrosinase with 14 target ligands.
Molecular Interaction Analysis of the Docked Complexes
The PDB files (protein data bank) were downloaded of the top three docked complexes of laccase, peroxidase MSP1, and tyrosinase enzymes, and Discovery Studio was employed to analyze the interactions between pollutants and each enzyme. The distance between laccase and Congo red, Mordant black 11, and Indigocarmine ranges from 2.18 to 5.22 Å. The distance between peroxidase MSP1 and Congo red, Mordant black 11, and Acid blue 93 ranges from 2.03 to 5.33 Å. The distance between tyrosinase and Congo red, Acid yellow 99, and Indigocarmine ranges from 1.75 to 5.42 Å. As a result, these distances among the two enzymes and three substrates suggest the existence of hydrogen bonds (2.8–3.4 Å) and van der Waals forces (3.8–4.2 Å). The results are shown in Table 7 on laccase, peroxidase, and tyrosinase. Interaction studies of each enzyme and azo dye at the molecular level are shown in Figures 8–10.
Amino acid residues, distance range, and types of bonds present between each enzyme (laccase, peroxidase MSP1, and tyrosinase) and the pollutants.

Molecular interaction analysis of laccase with (a) Congo red, (b) Mordant black 11, and (c) Indigocarmine.

Molecular interaction analysis of peroxidase MSP1 with (a) Congo red, (b) Mordant black 11, and (c) Acid blue 93.

Molecular interaction analysis of tyrosinase with (a) Congo red, (b) acid yellow 99, and (c) Indigocarmine.
Protein–Protein Interaction Analysis
The construction of protein–protein interaction (PPI) networks using the STRING database provides a visual and analytical framework to understand the relationships between proteins. These networks reveal functionally coherent modules—groups of proteins that collaborate within specific biological processes. In our analysis of enzymes from

Laccase interaction network predicted by String.

Tyrosinase interaction network predicted by String.

Dye-decolorizing peroxidase MSP1 Interaction network predicted by String.
Co-expression of Gene
The co-expression analysis of the laccase, tyrosinase, and dye-degrading peroxidase MSP1 genes in

Co-expression analysis of the gene encoding laccase in

Co-expression analysis of the gene encoding tyrosinase in

Co-expression analysis of the gene encoding dye-decolorizing peroxidase MSP1 in
Discussion
The biodegradation of azo dyes presents a significant environmental challenge due to the stability of these synthetic compounds and their potential toxicity, particularly when discharged into aquatic systems from textile industry effluents. 24 Azo dyes, widely used in the textile industry, are notorious for their persistent nature due to the azo bonds linking aromatic rings, which render them resistant to natural degradation processes. 25 This persistence can lead to the accumulation of dyes in water bodies, posing severe risks to aquatic life and human health, including toxic, mutagenic, and carcinogenic effects. 26 Bioremediation provides an eco-friendly solution to azo dye pollution by using enzymes from microorganisms like fungi and bacteria. It provides a sustainable alternative to conventional chemical treatments, which often involve hazardous substances and generate harmful by-products. 27 Enzymatic degradation, by contrast, typically operates under mild conditions and produces fewer by-products, making it a greener option. Additionally, the ability to engineer and optimize these enzymes for specific dye degradation tasks enhances their efficiency and applicability in various industrial and environmental settings. 28 The development of enzyme-based treatments represents a promising direction for addressing pollution from the textile industry and other sectors that utilize azo dyes.29,30
In this study, the retrieved structures of laccase, tyrosinase, and peroxidase from NCBI were analyzed for their physicochemical properties, indicating that all three proteins are stable, as evidenced by their instability indices being less than 40. An instability index below 40 suggests that these proteins have a stable structure, making them suitable for functional applications in environmental bioremediation processes, such as the degradation of azo dyes.
31
This stability is crucial for the effectiveness of these enzymes in industrial applications, where they are exposed to various environmental conditions. The findings align with previous studies where protein stability was directly linked to enzyme efficiency and longevity during catalytic processes.
32
In contrast to previous studies where bacterial enzymes were used, our study focuses on fungal-derived enzymes, which demonstrated greater predicted stability, suggesting their potential suitability for long-term and large-scale applications. Compared to enzymes from
Moreover, molecular docking analysis was performed against 14 azo dyes with laccase, peroxidase (MSP1), and tyrosinase to check their potential against dyes. The results show that the enzymatic degradation of azo dyes by laccase, peroxidase (MSP1), and tyrosinase is characterized by specific molecular interactions that play a crucial role in the breakdown process. The docking analysis revealed that all the proteins—laccase, peroxidase, and tyrosinase—exhibit the highest binding affinity against Congo red dye, with values of −10.1, −10.6, and −11.3, respectively. Toxicology analysis also shows that Congo red dye is carcinogenic and causes toxicological effects. Suryamathi et al. examined the use of tyrosinase from
The interaction analysis of laccase, tyrosinase, and peroxidase MSP1 with azo dyes reveals distinct roles and effectiveness in the degradation process.
35
Laccase shows moderate interaction distances (2.18–5.22 Å) with substrates like Congo red, indicating its ability to form essential hydrogen bonds and van der Waals interactions for substrate binding and oxidative cleavage of azo bonds. In a study by Junior et al., the application of laccase from
The co-expression analysis of laccase, tyrosinase, and peroxidase genes in
The comparative analysis in this study indicates that while all three enzymes—laccase, peroxidase, and tyrosinase—can degrade azo dyes, peroxidase from
Conclusion
This study highlights the significant potential of laccase, tyrosinase, and peroxidase MSP1 from
