Biochemical and Structural Properties of a High-Temperature-Active Laccase from Bacillus pumilus and Its Application in the Decolorization of Food Dyes

A novel laccase gene isolated from Bacillus pumilus TCCC 11568 was expressed, and the recombinant laccase (rLAC) displayed maximal activity at 80 °C and at pH 6.0 against ABTS. rLAC maintained its structural integrity at a high temperature (355 K) compared to its tertiary structure at a low temperature (325 K), except for some minor adjustments of certain loops. However, those adjustments were presumed to be responsible for the formation of a more open access aisle that facilitated the binding of ABTS in the active site, resulting in a shorter distance between the catalytic residue and the elevated binding energy. Additionally, rLAC showed good thermostability (≤70 °C) and pH stability over a wide range (3.0–10.0), and displayed high efficiency in decolorizing azo dyes that are applicable to the food industry. This work will improve our knowledge on the relationship of structure–function for thermophilic laccase, and provide a candidate for dye effluent treatment in the food industry.

Laccases can be acquired from different sources, such as bacteria, fungi, actinomycetes, and plants. Fungal laccases show more favorable characteristics for commercial applications in comparison to other laccases from plants and bacteria, because of their high reduction potential [7]. However, fungal laccases also demonstrate some drawbacks, including low fungal growth rate and being operative only under low temperature and pH conditions. They are often deactivated in some industrial sections, such as food and (0.1%, w/v) were added to the bacterial colonies to distinguish the laccase-producing strains by their pink color. Then the obtained colonies that showed a pink color were picked up and purified using single-colony separation.

Phylogenetic Analysis of the Laccase-Producing Strain
The 16S rDNA was amplified with the genomic DNA of the laccase-producing strain according to the method described previously [17]. The homologous sequences of the resulting 16S rDNA sequence were searched with GenBank BLAST (http://www.ncbi.nlm. nih.gov/BLAST/ accessed on 12 March 2020). Subsequently, a bootstrap consensus tree was built using MEGA 6.0 software with the neighbor-joining method for phylogenetic analysis [18]. To further identify this strain, taxonomic analysis was conducted as noted in Bergey's Manual of Determinative Bacteriology.

Heterologous Expression of Laccase
The gene encoding laccase (lac) was amplified with two primers, Lac-F (5 -CGCGG ATCCGATGGCACTGGAAAAATTTG-3 ; the underlined bases encode the BamHI site) and Lac-R (5 -ACGCGTCGACCTGCTTATCCGTGACGTCC-3 ; the underlined bases encode the SalI site) from the genomic DNA of B. pumilus. After double digestion with BamHI and SalI, it was integrated into the BamHI-SalI-linearized expression plasmid pET-22b (+) to construct the plasmid pET-lac. It was then transformed into E. coli BL21 (DE3). A positive colony was initially selected to grow at 37 • C in LB medium containing ampicillin (100 µg/mL) for 12 h. Subsequently, the preculture (1 mL) was diluted with 50 mL of fresh LB medium with ampicillin (100 µg/mL). When OD 600 reached 0.6-0.8, 1 mM isopropyl-β-D-1-thiogalactopyranoside (IPTG) was used to induce the recombinant laccase (rLAC) at 16 • C for 20 h. E. coli BL21 cells harboring the empty plasmid pET-22b (+) were used as the control.

Molecular Docking and Molecular Dynamics Simulation Analysis of rLAC
The software AutoDock Vina was used to perform molecular docking between the rLAC and the ligand. The 3D structure of rLAC was BLAST-searched in the UniProt database using SWISS-model (https://swissmodel.expasy.org/ accessed on 15 May 2021). A protein structure (Protein Data Bank entry 1GSK) with more than 70% homology was used as the template for homology modeling. The 3D structure of ABTS was obtained from PubChem (https://pubchem.ncbi.nlm.nih.gov/ accessed on 15 May 2021). A method of semi-flexible docking was utilized to allow the chemical bonds of the ligand to freely rotate, while fixing the coordinates of the atoms in rLAC. The grid parameters were set to build a sphere grid near the active sites of rLAC1 (H419, C492, H497, and M502) with a radius of 5.0 Å in the docking area.
The docking conformation with the lowest protein-ligand binding energy was selected for molecular dynamics simulation by Gromacs 5.1.4, combined with the AMBER99SB forcefield. The parameters for the molecular dynamics simulations were set based on the previous study [19,20], in which the SPC216 model was set as the water molecule. The simulation system was neutralized by a total of 11 Na + , which were added into a cubic box with a distance of 1.5 nm in all three dimensions. The incorrect geometries and collisions among the atoms in the simulation system were removed by the steepest descent minimization method, using a maximum number of 50,000 steps. Equilibration simulations for the solvent around the protein were performed under constant volume-temperature (NVT) and constant pressure-temperature (NPT) ensembles under harmonic restraints. Short-range van der Waals interactions and short-range electrostatic interactions were truncated smoothly at 1 nm. The leap-frog integrator algorithm with a step size of 2 fs was utilized to integrate the equation of motion for equilibrium dynamics. The LINCS algorithm was used to constrain the hydrogen bonds. The Particle Mesh Ewald (PME) method was used to calculate the long-range electrostatic interactions with a grid size of 0.16 Å. The pressure (1 bar) and temperature (300 K) under the isothermal-isostatic (NPT) ensemble were maintained using Parrinello-Rahman and the V-rescale method, respectively. Both the NVT and NPT simulation times were set as 100 ps and the final MD simulations were carried out for 100 ns. The root-mean-square fluctuation (RMSF) and the root-mean-square deviation (RMSD) were analyzed using the gmx rmsf and gmx rms commands respectively. The snapshots were performed using Visual Molecular Dynamics (VMD) 1.9.4.
The approach of MM/PBSA (Molecular Mechanics Poisson Boltzmann Surface Area) was utilized for energy decomposition according to the trajectory of the molecular dynamics simulations. The binding energy within a time period of 0~10 ns was calculated using the tool g_mmpbsa.py and the software Numpy (Numerical Python). The binding free energy is calculated using Equation (1): Among them, ∆G bind stands for the binding free-energy, ∆E ele stands for the electrostatic interaction, ∆E vdw stands for the van der Waals force, ∆G SASA stands for the non-polar contribution energy calculated using the solvent accessible surface (SASA), and ∆G solv stands for the solvent free-energy.

Purification of rLAC
Cells expressing rLAC were centrifugally harvested via centrifugation at 8000× g and 4 • C for 15 min, and then resuspended in 20 mM Tris-HCl buffer (pH 7.0) with 500 mM NaCl and 20 mM imidazole. Subsequently, the cells were disrupted via sonication at 320 W with 4 s strokes and 3 s intervals. After centrifugation (12,000 × g, 30 min), the supernatant was injected into a nickel-nitrilotriacetic acid (Ni-NTA) agarose gel column (Shenggong, Shanghai, China). rLAC was eluted using Tris-buffer containing 500 mM imidazole and 500 mM NaCl (pH 7.0) after removing the impurities. After dialysis, the purity and molecular mass of rLAC were determined using sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) analysis.

Enzyme Assay
The activity of rLAC was determined with ABTS as the substrate at 80 • C using the previously described method [21,22]. Concisely, the enhanced absorbance at 420 nm (ε 420 = 36,000 M −1 cm −1 ) was monitored for ABTS (6 mM) oxidation in 0.1 M citratephosphate buffer (pH 6.0). One unit of the activity for laccase was designated as the amount of laccase that oxidized 1 µmol of substrate per minute.

Characterization of rLAC
The optimal temperature of rLAC was determined with ABTS as the substrate at different temperatures from 30 • C to 90 • C. To identify the optimal pH of rLAC, the enzymatic reaction was conducted over a pH range of 3.0 to 9.0 with citrate-phosphate buffer (0.1 M, pH 3.0 to 8.0) or glycine-sodium hydroxide buffer (0.1 M, pH 9.0). The relative enzymatic activity was determined with the maximal activity as 100%.
To identify the thermostability of rLAC, it was incubated at three temperatures (60 • C, 70 • C, and 80 • C) for 0 h to 2 h without substrate. The pH stability was identified by keeping the purified rLAC at pH 3.0, 7.0, and 10.0 (4 • C) for various periods (0-10 days). After incubation, the residual activities of the samples were measured with the substrate ABTS under optimal conditions (80 • C, pH 6.0). The activity of the initial rLAC was taken as 100%.

Identification of the Strain with Laccase Activity
Numerous studies have been focused upon bacterial laccases over the past decade, especially for those derived from the Bacillus genus, due to their advantages over fungal laccases, such as a higher pH stability, wider pH adaptation, and thermostability [12,23]. In this work, a bacterial strain demonstrated a high degree of laccase activity towards syringaldazine among 216 bacteria screened from a soil sample gathered from Guangdong Province of China. After aligning the 16S rDNA sequence of this strain with BLAST included in the NCBI database, the result suggested this strain was a member of the Bacillus family. Phylogenetic tree analysis suggested that this strain was most closely associated with B. pumilus ( Figure S1 in Supplementary Materials).
The morphological properties of this strain are as follows: rod-shaped, Gram-positive, rough-surfaced, and with spore-producing colonies (data not shown). Therefore, it was finally characterized and named B. pumilus TCCC 11568 on the basis of its biochemical tests and morphological characteristics, along with the 16S rDNA analysis. In summary, we isolated a new Bacillus strain showing laccase activity from a soil sample, and our study enriched the family of laccase-producing bacteria.
In this work, a significant band of approximately 57 kDa was observed on the SDS-PAGE gel of the cell extracts, while no band was found in the cell extract of the control ( Figure S3a). A target band with a molecular mass of approximately 57 kDa was detected on the SDS-PAGE gel after purification ( Figure S3b). The obtained molecular mass was similar to some other reported bacterial laccases, including B. amyloliquefaciens [27] and B. halodurans C-125 (56 kDa) [32]. Its specific activity towards ABTS was 190 U/mg, which was notably higher than other Bacillus laccases to ABTS, such as B. amyloliquefaciens LC02 (20.7 ± 1.2 U/mg for ABTS) [33]. According to the sequence alignment, protein modelling, and enzyme activity analysis, it can be concluded that a new laccase with a high level of activity and a conserved catalytic motif that has been observed in other reported laccases was obtained in our study.

Effect of Temperature and pH on the Activity and Stability of rLAC
The maximum activity of rLAC for ABTS oxidation was detected at 80 • C, and a relatively high activity level was sustained at 60-90 • C (Figure 2a). To investigate its thermostability, the residual activities were determined after incubation at differing temperatures (50-80 • C) for 0 min to 120 min (Figure 2c). It suggested that rLAC was quite stable at 50 • C or 60 • C, and retained 87.3% or 70% of the original activity at these temperatures after a 120-min incubation, respectively. Although rLAC showed an evidently decreased stability after incubation at a high temperature (70 • C or 80 • C), it still retained 50.3% and 24.6% of the original activity after a 120-min incubation, respectively (Figure 2c). The high optimum temperature was one of the most significant properties demonstrated by the bacterial laccases, especially for those derived from the Bacillus genus, such as 80 • C for CotA laccase from B. amyloliquefaciens TCCC 111018 [27], B. velezensis [29], and B. licheniformis [34]. By contrast, a similar phenomenon was not observed for fungal laccases whose optimum temperatures were usually below 60 • C. The optimal temperatures for laccases from Cerrena unicolorstrain GSM-01, Trametes sp. F1635, and Aureobasidium melanogenum strain 11-1 were 45 • C, 50 • C, and 40 • C, respectively [35][36][37]. Compared to fungal laccases, the high thermostability is another impressive feature of bacterial laccases. The thermostability of rLAC from B. pumilus TCCC 11568 was similar to laccases isolated from other Bacillus strains, such as B. amyloliquefaciens TCCC 111018 [27], and B. tequilensis SN4 [38]. Moreover, its thermostability was even higher than B. licheniformis DSM 13 laccase, which was deprived of approximately 92% of its initial activity after 1 h of incubation at 80 • C [33]. Thus, rLAC might be directly applied in treating hot effluents from the dyeing process, which is generally performed under high-temperature conditions [39]. In addition, using thermostable laccase to treat hot textile effluents is economically attractive because the recycling of hot water saves a significant amount of energy in the process.
rLAC was capable of oxidizing ABTS over a wide pH range and demonstrated its maximum activity at pH 6.0 (Figure 2b), which was higher than most Bacillus laccases, such as B. velezensis TCCC 111904 [29], B. amyloliquefaciens TCCC 111018 [27], B. vallismortis fmb-103 [26], B. subtilis X1 [28], and B. clausii KSM-K16 [13]. Generally, fungal laccases are merely stable under acidic and neutral pH conditions, but bacterial laccases are stable under alkaline conditions [40]. Additionally, rLAC exhibited a higher degree of stability (Figure 2d) than other Bacillus laccases over a wide pH range, such as the B. velezensis TCCC 111904 laccase, which retained 67.6% of its initial activity after incubation at 4 • C and pH 9.0 for 24 h [29]. Additionally, rLAC also showed quite a good degree of stability in an acidic environment, retaining 77.1% of its original activity after incubation at pH 3.0 for 10 d. The above results demonstrate that rLAC from B. pumilus TCCC 11568 has a high potential for applications in a wide range of pH environments, especially under alkaline conditions.
In summary, the B. pumilus TCCC 11568 laccase demonstrated thermophilic features and a higher degree of thermostability than some reported Bacillus laccases. Therefore, this laccase has excellent potential for becoming a new industrial enzyme formulation.

Influence of Metal Ions and Inhibitors on the Activity of rLAC
rLAC activity was seriously affected by 5 mM Co 2+ and Mn 2+ , retaining 24% and 9.4% of the control, respectively, but other metal ions did not lead to critical activity loss, indicating the high tolerance of rLAC towards most metal ions (Table 1). A comparable phenomenon was also found with the B. safensis sp. S31 laccase, which retained only 13.2% activity when incubated with 1 mM Mn 2+ [41]. Electrolytes are required to obtain a high degree of efficiency in the industrial dyeing process. NaCl is the most extensively used electrolyte, and its concentration is approximately 25-30 g/L (around 0.5 M). Most fungal laccases are deactivated in solutions containing over 100 mM NaCl because they are intrinsically sensitive to halides [42]. This is a major barrier for applying fungal laccases to treat wastewater containing chloride ions at a high concentration. By contrast, some Bacillus laccases were highly tolerant to chloride, and their activities were even promoted by increased NaCl concentrations [32]. Different tolerances to halide inhibition among fungal and bacterial laccases may result from their different target localizations after protein translation. Laccases from fungal species are primarily secreted for degrading lignin, but Bacillus laccases anchor in the spore coat. This protein produced a high degree of stability as an indirect consequence of its evolutionary constraints, since it must adapt to a tightly packed protein coat structure [43]. Here, we found that rLAC from B. pumilus TCCC 11568 also exhibited a high tolerance to NaCl, and it even retained 50.7% activity of the control in a solution containing 500 mM NaCl. Similar phenomena were also reported for laccases from B. pumilus W3 and B. vallismortis fmb-103 [26,44]. Therefore, rLAC with its high salinity tolerance would demonstrate more advantages in dealing with dye effluents that are high in salinity [45]. Additionally, no vital activity loss was found during incubation with 5 mM EDTA (62.1% activity retained, Table 1). On the contrary, a complete loss was found with laccase isoforms of Aspergillus ochraceus and Thermus thermophilus in 1 mM EDTA [46,47].
The effects of some reported inhibitors, such as L-cysteine, sodium dodecyl sulfate (SDS), β-mercaptoethanol, and dithiothreitol (DTT) on their rLAC activity are also exhibited in Table 1. rLAC activity was greatly decreased by 0.5 mM DTT, L-cysteine, and β-mercaptoethanol; this phenomenon was also noted in other fungal and bacterial laccases [26,44,48]. Additionally, laccases from diverse Bacillus strains also demonstrated differing amounts of tolerance against the same inhibitor. SDS (0.1 mM) significantly inhibited B. safensis sp. S31 laccase activity, retaining only 1.3% the activity of the control [41]. On the contrary, the activity of rLAC was not even affected by SDS at a high concentration (5 mM). A similar phenomenon was also noted for the B. vallismortis fmb-103 laccase [26].

Investigation of the Thermophilic Features of rLAC via Molecular Docking and MD Simulations
To present a reasonable explanation for the thermo-stable features of rLAC, MD simulations were performed at 325 K, 355 K, and 365 K, respectively. The root mean square deviation (RMSD), displayed as a function of simulation time, is a key indicator for assessing the stability of protein structure upon the binding of a ligand. In general, the lower the RMSD value displayed by the protein complex during MD simulations, the higher the stability shown [49]. Our results suggested that the complex systems reached equilibrium in the 75 ns simulation (Figure 3a). The RMSD for the rLAC-ABTS complex at high temperatures (355 K and 365 K) was lower in comparison to that at a low temperature (325 K), indicating an enhanced degree of stability for the complex upon its increase in temperature. Rg, calculated as the root mean square distance from the center of mass, indicates the level of the secondary structure backbone rigidity and the compactness of the protein complex system. The Rg values of the rLAC-ABTS complex at a high temperature were larger than those at low temperature, suggesting a reduction in the compactness and structural backbone rigidity (Figure 3b). The above results suggested that the stability of the rLAC-ABTS complex is not merely determined by the protein compactness, but by a synergic conformation change. In fact, the subsequent RMSF analysis that is commonly used to characterize the structural fluctuation and protein structural integrity also displayed expansion in some of the loop regions of rLAC, indicating a decrease in protein compactness in some regions, along with an increase in temperature from 325 K to 365 K. Actually, the RMSF values for all rLAC residues under high-temperature conditions (355 K and 365 K) did not demonstrate remarkable differences to those of rLAC under mild temperature conditions (325 K), except for a few residues that were all within the loop region with the largest enhancement of RMSF values, such as Glu90, Leu219, Gly323, and Asp379 (Figures 3c and 4). The distance between residues Glu90 and Asp379, Glu90 and Leu219, Asp379, and Leu210 gradually increased along with the elevated temperature. As a result, the loop regions of I, II, and IV were far away from each other at high temperatures (Figure 4b, 355 K; Figure 4c, 365 K) at the end of the MD simulations. On the contrary, the distance between loop region II and loop region III became shorter. In conclusion, the overall conformation of rLAC did not show significant alterations at high temperatures (355 K, 365 K), but minor location adjustments in some loop regions occurred at high temperatures. As can be supposed, the good flexibility of the loop region may potentially make it simple to expose the active center, and may therefore facilitate the binding of rLAC with the substrate (ABTS) and the subsequent catalytic reaction. In fact, this speculation was also supported by observations of the surface representation and binding modes of ABTS with rLAC ( Figure 5). As can be seen from the trajectories recorded after 100 ns of MD simulations, a different spatial conformation of the binding pocket in rLAC at different temperatures was observed ( Figure 5). A narrow aisle for substrate access to the active site formed under low temperatures, which obviously did not facilitate effective binding of ABTS with the rLAC in the right manner, nor subsequent catalysis (Figure 5a). By contrast, the conformation changed at high temperatures (355 K, 365 K), and a more open space was formed. As a result, ABTS could enter the binding pocket easily and it formed a tighter complex with rLAC in the active site, which was one of the most important prerequisites for an effective catalytic reaction by laccase (Figure 5b,c).
The distance between the catalytic residue and ABTS was also analyzed to explore the influence of temperature on the catalysis performance of rLAC. The distance between the C-atom of H497 located in the active site and the C-atom of ABTS decreases along with an increase in temperature, indicating that a closer interaction between ABTS and rLAC occurred at the active center ( Figure 6). In fact, the released binding free-energy analysis for rLAC-ABTS verified this phenomenon further (Table 2).  The electrostatic interaction and the van der Waals interaction made a major contribution to the production of the binding energy during the binding process of ABTS in the active center of rLAC. Moreover, the released binding energy was much higher for rLAC-ABTS at high temperatures, especially at 355 K. These results also matched the RMSD analysis in Figure 3, in which the RMSD value obtained at 355 K was the lowest, indicating a high thermostability of the rLAC-substrate complex and a much higher substrate-binding affinity at high temperature. Mollania et al. (2017) investigated the variation of the protein structure of the laccase obtained from a local Bacillus species (HR03) after thermal activation, using gas-phase electrophoretic mobility macromolecule analysis, far-UV CD-spectra, common biochemical methods, fluorescence analysis, etc. [50]. They speculated that the improved activity of Bacillus sp. HR03 laccase was not the result of oligomerization, but of the generation of a more active conformation due to the thermal treatment at a high temperature (70 • C). In other words, the incubation of Bacillus laccase at high temperatures did not severely affect the overall structure in an irreversible way, but instead led to a more beneficial conformation of the Cu sites and of the active site for substrate binding. According to our literature investigation, the mechanism for the high thermophilic features of Bacillus laccase is still not very clear. We herein presented another reasonable and vivid explanation for interpreting the thermophilic features of rLAC via docking and molecular dynamics simulations analyses. We hypothesized that the formation of a correctly matching substrate access aisle to the active site, through the minor adjustment of protein spatial conformation, was responsible for the improved catalytic activity of Bacillus laccase when incubated at high temperatures. As a result, the released binding energy greatly increased, due to the tighter binding of ABTS at the active site driving the accelerated catalytic reaction.

Dye Decolorization
Color is a vital quality attribute of food in stimulating customers' appetites and is, therefore, one of the most important concerns for food manufacturers when considering market acceptance. Synthetic dyes as colorants are widely applied in the food industry to dye foodstuffs such as mustard, sweets, jams, cakes, and juices beverages, and they are also used in the production of drugs and cosmetics, because they have certain fascinating advantages compared to natural pigments, including brighter colors, low production costs, high water solubility, and high fading resistance against exposure to chemicals, light, and water [51].
However, large amounts of effluent with intense colors are produced after the dyeing process, as approximately 10-20% of the dyes are lost during this process [52]. Therefore, without treatment, the discharged effluents are an important source of pollution for water bodies, and bring in severe damage to aquatic ecosystems through impeding light penetration to the water bodies, which consequently reduces the photosynthetic rate of aquatic plants and the dissolved oxygen levels in the water [23,53]. In recent years, more and more bioremediation methods with favorable properties, such as higher efficiencies, lower costs, and eco-friendly properties, have been developed and applied to decolorizing dye effluents [54][55][56].
Laccases can degrade numerous environmental pollutants, including micropollutants, personal-care products, and textile dyes [6,57]. Since the typical dye effluents are generally identified via their high temperatures and high pH values (e.g., over 40 • C and around pH 9.0) [39], most fungal laccases will be deactivated under these harsh conditions; e.g., laccases from Sclerotium rolfsii and T. hirsuta can only work best with degrading indigo carmine under an acidic pH environment [45,58]. Thus, for bacterial laccases, conditions of much higher pH stability and thermostability show a higher potential for treating dye effluents.
In this study, a new laccase with a high thermostability and pH stability was obtained from the newly identified bacterium B. pumilus TCCC 11568 ( Figure 2). Therefore, we assessed its decolorization ability against several azo dyes at 60 • C and over a broad pH range, starting from pH 5.0 to 9.0 ( Figure 7). The decolorization rates for each azo dye were all above 70% under optimum conditions (an appropriate redox mediator and pH environment) (Figure 7). Taking the dye Sunset Yellow for example, the decolorization rate was below 20% or 40% when using ABTS or syringaldehyde as a redox mediator, respectively. However, the decolorization rate soared up to 84.4% when incubated with syringaldehyde at pH 9.0 ( Figure 7). Additionally, we also found that rLAC displayed a higher decolorizing ability against all of the dyes in a neutral (pH 7.0) or alkaline environment (pH 9.0), except for the dye etythrosine, which is more easily decolorized in an acid environment. This differs from B. velezensis TCCC 111904 laccase, which was more effective for decolorizing the dyes when ABTS was used as the mediator at pH 5.5 and pH 7.0 [25]. Therefore, it was inferred that synthetic dyes with diverse structures could be efficiently decolorized by rLAC isolated from B. pumilus TCCC 11568 over a broad pH range (5.0-9.0) at a high temperature (60 • C) with an appropriate mediator. However, depending on the extensive specificity of the CotA laccase, no tendency of laccase could be summarized for the decolorizing activity towards synthetic dyes with different structures [59].
As is known, most dye effluents are released at high temperatures [39]. Thus, it would be helpful to decrease the unnecessary cost of cooling the dye effluents to meet the required reaction conditions due to the alkaline stability and high thermostability of rLAC. Moreover, recycling hot water after decolorization could save a lot of energy. Therefore, the B. pumilus TCCC 11568 laccase presented promising applications because of its excellent efficiency in decolorizing dyes under high temperatures and alkaline conditions.

Conclusions
The recombinant laccase (rLAC) of Bacillus pumilus TCCC 11568 demonstrated a high thermostability and pH stability over a wide pH range. The overall conformation of rLAC did not show significant alterations at a high temperature (355 K), but minor location adjustments in some loop regions occurred at a high temperature. However, those adjustments were presumed to be responsible for the formation of a more open access aisle that facilitated ABTS binding in the active site, resulting in a shorter distance from the catalytic residue and elevated binding energy. Due to its thermophilic feature, rLAC could effectively decolorize azo dyes at high temperatures over an extensive pH range in the presence of the appropriate mediator. These aforementioned features make rLAC a potential candidate for industrial applications in dye decolorization.

Supplementary Materials:
The following supporting information can be downloaded at: https://www. mdpi.com/article/10.3390/foods11101387/s1, Figure S1: 16S rDNA-based phylogenetic comparison of B. pumilus TCCC 11568 and other Bacillus species. The phylogenetic tree was engineered using MEGA 6.0 software with the neighbor-joining method. The numbers at branch points suggested the percentage of bootstrap sampling from 1000 replications; Figure S2

Conflicts of Interest:
The authors declare no conflict of interest.