Table 1. Details of azo dyes used.

The dyes used in the study consist of Amido Black 10B (AB10B), Reactive Black 5 (RB5), Reactive Blue 160 (RB160). The information concerning the azo dyes, including their molecular formula and structure, is illustrated in Table 1.

The pair-wise sequence alignment of laccase (uniport id: Q5EBY5) with Iron permease FTR1, Cu-oxidase-domain-containing protein, Solute carrier family 40 protein, Glyco_hydro_79C domain-containing protein and several uncharacterized proteins, and heme peroxidase (uniport id: A0A1Q3E2I9) from Lentinula edodes was performed by BLAST analysis using a non-redundant sequence database. The top 10 hits were selected and the multiple-sequence alignment of these sequences was performed using Clustal Omega. The phylogenetic tree from multiple sequence alignment was generated to deduce the evolutionary relationship between the high similarity sequences [9][10].

Protein network analysis was performed for laccase (uniport id: Q5EBY5) and heme peroxidase (uniport id: A0A1Q3E2I9) using STRING protein-protein interaction database by selecting Lentinus tigrinus and Lentinula edodes as organisms respectively. The proteins found to be interacting were evaluated for their functional association with the laccase and heme peroxidase.

The protein structure was predicted by using Alpha Fold (https://colab.research.google.com/github/sokrypton/ColabFold/blob/main/AlphaFold2.ipynb) developed by DeepMind and the quality of the model was evaluated by generating the Ramachandran plot using Discovery Studio software. The quality of the structure was evaluated based on the number of residues in the favorable region of the Ramachandran Plot.

The modeled structure of laccase (Q5EBY5) and heme peroxidase (A0A1Q3E2I9) were used for molecular docking with the dye molecules including Amido Black 10B, Reactive Black 5, and Reactive Blue 160. The 3D structures of the dye molecules were prepared using an online smile translator. Molecular docking was performed using AutoDock tools 1.5.6 version using search parameters of Lamarckian GA algorithm 4.2. The grid box dimensions considered were center X = −13.65, center Y = −4.18, and center Z = −5.79 and the grid size was X = 23.62, Y = 23.55, and Z = 28.73 for laccase (Q5EBY5) and center X = −2.76, center Y = −1.18, and center Z = 26.41 and the grid size was X = 28.6, Y = 29.32, ad Z = 32.49 for heme peroxidase (A0A1Q3E2I9). The genetic algorithm was used as search parameters including the number of GA runs of 10. The clustering of docked conformations was done using default RMSD tolerance of 2.0 Å. The proteins were prepared for docking by polar hydrogen addition, merging of non-polar hydrogens and addition of Kollman charges. The high negative binding energy of docking and the number of hydrogen bond interactions at the active site of the enzyme were considered for the evaluation of the best binding pose [11]-[13].

Table 2. Details of proteins interacting with laccase Q5EBY5.

The complexes of protein-dye molecules having high negative binding energy were considered for MD simulation for a duration of 100 ns. The ligands topologies were made using CHARMM General Force Field online server. The protein topology was made using CHARMM36 force field. The systems were solvated in dodecahedron box using Simple Point Charge water model and ions were added to neutralize the systems. The energy minimization was performed keeping maximum force limit of 10 kJ/mol. The equilibrations at number of atom pressure temperature and number of atom volume temperature at 300 K and 1 bar pressure for a duration of 100 ps each were performed for all the systems. Lastly, the systems were exposed to molecular dynamic run of 100 ns. The analysis of simulation runs was performed by considering the variations in Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), Radius of gyration (Rg) and number of hydrogen bonds formed at each picosecond during the simulation trajectory [14].

The three-dimensional structure of the Lentinus sp. laccase and heme peroxidase were modelled using AlphaFold. Docking analysis serves as a valuable tool for investigating the interaction between proteins and ligands to elucidate the formation of stable docked complexes. This approach is instrumental in exploring both the binding and would also enable in deciphering the mechanistics of enzyme derived catalysis [14][15].

Sequence based alignment of laccase of Lentinus(Uniprot ID: Q5EBY5) showed high sequence identity of 87.93% with laccase of Polyporus brumalis 82.54% with laccase of Ganoderma lucidum and 82.72% with laccase of Ganoderma weberianum. Further multiple sequence alignment followed by phylogenetic analysis revealed close evolutionary affiliation of laccase of Lentinus tigrinus with Polyporus brumalis (ABN13591.1), Cerioporus squamosus (KAI0699344.1), Ganoderma lucidum (AHA83584.1) while the considered laccase is relatively less evolutionarily related with Ganoderma weberianum (ANA53145.1) and multiple sequence alignment analysis of laccase showed 100% similarity with Lentinus (Figure 1(a), Figure 1(b)).

(a)

(b)

Figure 1. Phylogenetic tree (a) and multiple sequence alignment (b) analysis of laccase (Q5EBY5).

Network analysis of Laccase of Lentinus(Uniprot ID: Q5EBY5) displayed close interaction of laccase (Figure 2) with Iron permease FTR1 (A0A5C2S0T8), Cu-oxidase-domain-containing protein (A0A5C2S1A8), Solute carrier family 40 protein (A0A5C2S2R5), Glyco_hydro_79C domain-containing protein (A0A5C2SSW2) and several uncharacterized proteins (A0A5C2SIH4, A0A5C2SWJ8) (Table 2).

The Ramachandran plot analysis denoted most of the residues in favorable region explaining the high accuracy of the structure (Figure 3(a)).In silico analyses illustrate the enzymatic capability of Lentinus sp. laccase in mitigating diazo dyes through the formation of hydrogen bonds. The 3D structure of laccase (Uniprot ID: Q5EBY5) was modelled using Alpha Fold which is an AI based high accuracy protein structure prediction tool (Figure 3(b)).

Figure 2. String network of proteins interacting with laccase (Red).

Table 3. Structural and catalytic features of laccase and heme peroxidase used for molecular docking.

(a)

(b)

Figure 3. Ramachandran Plot (a) and modeled structure (b) of Q5EBY5 laccase of Lentinus.

Molecular docking analysis revealed significant binding of three azo dyes, Amido black 10B (AB10B), Reactive Black 5 (RB5) and Reactive blue 160 (RB160), at the catalytic site of laccase (Q5EBY5) with significant binding energies of −7.6 kcal/mol, −6.9 kcal/mol and −7.4 kcal/mol respectively. Further interaction analysis revealed that Amido black 10B formed 3 hydrogen bonds with Gly101, Val402, and Gln441 (Figure 4(a), Figure 4(b)) [16]-[18]. Additionally, these analyses indicate a slight shift in the enzyme’s conformational state due to laccase-azo dyes interaction. These findings encourage the refinement of technologies tailored for synthetic dye treatment.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 4. Docked complex of Amido black 10B with laccase (a), representation of molecular interactions in the organic residue of AB10B dye from the enzyme-ligand complex (b). Gly101, Val402, and Gln441—conventional hydrogen bonding; PHE69—Pi-sulfur; MET 310—Pi-Sigma; PRO313—Pi-Alkyl; THR306—Pi-donor hydrogen bond. Docked complex of Reactive Blue 160 with laccase (c), representation of molecular interactions in the organic residue of RB160 dye from the enzyme-ligand complex (d), LYS71, ASP443, HIS401, VAL402, THR308, AND SER227—ASP128 conventional hydrogen bonding; GLY101—carbon-hydrogen bond; MET310—Pi-Sigma; PRO132—Pi-Alkyl. Docked complex of Reactive Black 5 with laccase (e), representation of molecular interactions in the organic residue of RB5 dye from the enzyme-ligand complex (f). ASP443, PHE440, THR308, SER227—conventional hydrogen bonding; PHE69, MET310—Pi-sulfur; PRO132—Pi-Alkyl.

The interactions of Reactive black 160 had demarcated 7 hydrogen bonds with Lys71, Asp443, His401, Val402, Thr308, Ser227 and Asp128 (Figure 4(c), Figure 4(d)).

Reactive Black 5 upon docking with laccase had resulted in 4 hydrogen bonds with Asp443, Phe440, Thr308, and Ser227 (Figure 4(e), Figure 4(f)). Thus, these indicate a potential interaction of the azo dyes with the active site of the laccase and hence further accelerate the laccase in the azo dye degradation (as shown Table 3) [14].

Figure 5. Analysis of MD simulation parameters for Q5EBY5-Amido black 10B complex stability, (a) Root mean square deviation (RMSD), (b) Root mean square fluctuation (RMSF), (c) Radius of gyration, (d) Hydrogen bonds analysis for Q5EBY5-Amido black 10B complex during the MD simulation trajectory of 100 ns.

Molecular dynamic simulation analysis was performed using GROMACS software for the high binding energy best docked complex of Q5EBY5 with AB10B to evaluate the stability of protein with the azo dye. Molecular dynamic simulation of apo-protein was also performed for comparative analysis. The RMSD analysis was done to assess the structural stability of the obtained complexes after virtual screening, and also to evaluate structural agreement with the crystal structure. The RMSF for each protein-ligand complex and native protein were calculated for fluctuation of the atom coordinates of the Cα to assess the flexibility of the structures, these values were evaluated with respect to the residues in the protein to estimate the motion in various residues. The RMSD values tend to converge after 20 ns and deviate from 0.12 to 0.32 nm with an average value of 0.22 nm. The RMSD values of apo and complex forms were closer, delineating stable interaction of the ligand at the binding site (Figure 5(a)). The RMSF values of apo and complex showed the small fluctuation within 0.5nm in turn demarcating high stability of ligand interaction (Figure 5(b)). The Rg values indicated high compactness of the structure in both apo and complex forms providing evidence that the ligand binding did not impact into any major change in the structure of enzyme (Figure 5(c)). This study was used to assess the changes brought about by ligand binding to the protein. The consistent hydrogen bonds formed during the trajectory of 100ns apparently prove effective interaction of the AB10B at the catalytic site of the enzyme (Figure 5(d)). Hence, these denote Q5EBY5 laccase as effective an competent enzyme that can lead into degradation of AB10B with higher efficiency [14][16][19].

(a)

(b)

Figure 6. Phylogenetic (a) and multiple sequence alignment analysis of A0A1Q3E2I9 heme peroxidase from Lentinula edodes (b) and with heme peroxidase of different species.

The multiple sequence alignment followed by phylogenetic analysis of heme peroxidase from Lentinula edodes(A0A1Q3E2I9) revealed the effective evolutionary relationship with heme peroxidase of Lentinula lateritia(KAJ4494436.1). These two closely related heme peroxidases in turn showed phylogenetic relationship with heme peroxidases of Lentinula novae-zelandiae(KAJ3865441.1) and Lentinula lateritia(KAJ3854868.1). Moreover, heme peroxidases of Lentinula aff. lateritia(KAJ3812978.1), Lentinula novae zelandiae, Lentinula aciculospora(KAJ4490906.1), Lentinula boryana(KAJ3993518.1), and Lentinula detonsa(KAJ3749620.1) were found to have significant sequence identity and phylogenetic association with the heme peroxidase of Lentinula edodes.Multiple sequence alignment analysis of heme peroxidase showed 100% similarity with Lentinus (Figure 6(a), Figure 6(b)).

The protein-protein network analysis revealed the closely associated proteins with that of heme peroxidase in Lentinula edodes (A0A1Q3E2I9) (Figure 7). The considered heme peroxidase closely interacts with two aspartate aminotransferases (A0A1Q3EHG3; A0A1Q3EG41) and one lipoxygenase (A0A1Q3ERW1). The other interacting enzymes with that of the considered heme peroxidase are given in Table 4 which include two lysophospholipase (A0A1Q3E3F4), glycoside hydrolase (A0A1Q3EQ81), protein transport protein (A0A1Q3ET19), and transmembrane protein (A0A1Q3E2D7).

Figure 7. Protein-protein interaction network of heme peroxidase A0A1Q3E2I9 (LENED_ 003065) in Lentinula edodes.

Table 4. Details of proteins interacting with heme peroxidase A0A1Q3E2I9 (LENED_003065).

The 3D structure of heme peroxidase (A0A1Q3E2I9) was predicted by Alpha Fold and the Ramachandran analysis was performed to evaluate the modeled structure quality (Figure 8(a), Figure 8(b)). Most of the residues are in favorable regions of the Ramachandran plot proving high quality of the predicted structure [20].

(a)

(b)

Figure 8. Ramachandran Plot (a) and modelled structure (b) of A0A1Q3E2I9 heme peroxidase of Lentinula edodes.

Molecular docking analysis of A0A1Q3E2I9 heme peroxidase of Lentinula edodes resulted in high negative binding energies of −7.7, −6.4 kcal/mol for Amido Black 10B form Arg980, Ile803 in conventional hydrogen bonding and Gly688-carbon-hydrogen bond (Figure 9(a), Figure 9(b)). The interaction of RB5 shows multiple hydrogen bonds involving Tyr622, Gln685, Ile803, Thr804, Thr818, Tyr955, Asp1061, in addition Gly688, Ala821, and Lys1060 involved in carbon hydrogen bond (Figure 9(c), Figure 9(d)). The interaction analysis showed several hydrogen bonds between heme peroxidase and AB10B and RB5. The interaction analysis of heme peroxidase with RB160 displayed unfavorable interactions of the ligand with the protein (as shown in Table 3) [21].

(a)

(b)

(c)

(d)

Figure 9. Docking interactions between heme peroxidase and Amido Black 10B (a), representation of molecular interactions in the organic residue of AB10B dye from the enzyme-ligand complex (b). ARG980, ILE803 conventional hydrogen bonding; GLY688—carbon-hydrogen bond; LEU814, THR818—Pi-Sigma. Docking interactions between heme peroxidase and Reactive Black 5 (c), representation of molecular interactions in the organic residue of RB5 dye from the enzyme-ligand complex (d). TYR622, GLN685, ILE803, THR804, THR818, TYR955, ASP1061—conventional hydrogen bonding; GLY688, ALA821, LYS1060—carbon-hydrogen bond; LEU814—Pi-Sigma; PHE822—Pi-Sulfur.

Molecular dynamic simulation of the apo-heme peroxidase and complexes of heme peroxidase with Amido Black 10B and Reactive Black 5 were performed. The RMSD values with respect to time of apo-protein as well as that of complexes were in the range of 0.15 to 0.5 nm with an average value of 0.3 nm during the trajectory denoting high stability of the complexes due to the least deviations in complexes (Figure 10(a)). Stable interaction have been observed earlier [22]. The RMSF values showed similar fluctuations in AB10B complex to that of the apo-heme peroxidase denoting relatively higher stability of the interactions compared to the other dye RB5 (Figure 10(b)). The overall fluctuations were less than 1nm thus both the ligands did not pose any major structural change in the protein and demarcated higher stability of complexes. The Rg values of both the complexes were similar to that of the apo-protein and varied in the range of 2.80 to 3.00 nm during the simulation trajectory thus attributing high compactness of the protein during simulation period of 100 ns (Figure 10(c)). The number of hydrogen bonds calculated at each picosecond in both the complexes between the protein and the ligand were consistent during the simulation and revealed effective interactions between protein and ligand during the 100 ns simulation (Figure 10(d)).

Several structural characteristics contributing to novel structural, catalytic, and stability features have been unveiled through bioinformatics driven genome-wide analysis [23][24].

Figure 10. Analysis of MD simulation parameters for A0A1Q3E2I9-ligand complex stability, (a) RMSD, (b) RMSF, (c) Rg, (d) Number of hydrogen bonds with respect to time in A0A1Q3E2I9-ligand complexes during the trajectory period of 100 ns.