Molecular docking and dynamics simulation analyses unraveling the differential enzymatic catalysis by plant and fungal laccases with respect to lignin biosynthesis and degradation

Isolation of the Antifungal Compound Alliodorin from the Heartwood of Cordia elaeagnoides A. DC. and the In Silico Analysis of the Laccase

Cordia elaeagnoides A. DC. is an endemic species of Mexico valued for its timber. Renowned for its durability, resistance, and versatile applications in medicine, this tree holds significant commercial importance. Tetrahydrofuran (THF) extract from the heartwood of C. elaeagnoides was studied. Through chromatographic column purification, the compound 8-(2,5-Dihydroxyphenyl)-2,6-dimethylocta-2,6-dienal, also known as alliodorin, was successfully isolated. Identification of alliodorin was confirmed through comprehensive analysis utilizing NMR, IR, and mass spectrometry techniques. Inhibition tests were conducted using both the THF extract and alliodorin against the rotting fungus Trametes versicolor (L.) Lloyd, employing the agar well diffusion assay. Remarkably, alliodorin exhibited 100% inhibition with a median lethal concentration of 0.079 mg/mL and a total lethal concentration of 0.127 mg/mL, in comparison to the commercial fungicide benomyl, which requires a concentration of 1 mg/mL. In silico analysis through molecular docking on the laccase enzyme was proposed in order to explain the inhibitory activity against the fungus T. versicolor, as this enzyme is one of the main sources of nutrients and development for the fungus. Based on these findings, we deduced that alliodorin holds promise as a potent antifungal agent, potentially applicable in a wide array of technological and environmentally friendly initiatives.

In silico and in vitro assays suggests Congo red dye degradation by a Lentinus sp. laccase enzyme

Laccase is a superfamily of ligninolytic enzymes known to degrade a wide variety of xenobiotics, including synthetic dyes. Congo Red (CR) has a diazo dye function, carcinogenic and mutagenic potential, and is currently applied in clinical analysis. The objective of this work was to produce and characterize the crude extract of Lentinus sp. in semi-solid fermentation (FSS) and perform in vitro and in silico studies to assess the potential of the crude extract to discolor the CR dye. Laccase activity was determined using ABTS as substrate and characterized. The in vitro discoloration was carried out using experimental design 22 at room temperature and monitored at 340 nm for 24h. Molecular docking and molecular dynamics simulations were performed between laccase and CR. The maximum laccase activity production was 29.63 U L-1 with six days of FSS. The optimal temperature and pH were 50 °C and 3.0, respectively. Discoloration of the CR dye was obtained only in tests containing CuSO4. Laccase formed stable complexes with the dye, presenting negative binding energy values ranging from -70.94 to -63.16 kcal mol-1 and the occurrence of seven hydrogen bonds. Molecular dynamics results showed the stability of the system (RMSD ranging from 1.0 to 2.5 Ä) and protein-ligand interaction along simulation. RMSF values pointed residues at the end of chains A (residues 300 to 305, 480 to 500) and B (residues 650 to 655 and 950 to 1000) as the most flexible regions of the laccase. This study highlighted the enzymatic action in the bioremediation of CR in vitro in agreement with the in silico simulations that demonstrate the enzyme potential.Communicated by Ramaswamy H. Sarma.

In Silico Insights into the Arsenic Binding Mechanism Deploying Application of Computational Biology-Based Toolsets

An assortment of environmental matrices includes arsenic (As) in its different oxidation states, which is often linked to concerns that pose a threat to public health worldwide. The current difficulty lies in addressing toxicological concerns and achieving sustained detoxification of As. Multiple conventional degradation methods are accessible; however, they are indeed labor-intensive, expensive, and reliant on prolonged laboratory evaluations. Molecular interaction and atomic level degradation mechanisms for enzyme-As exploration are, however, underexplored in those approaches. A feasible approach in this case for tackling this accompanying concern of As might be to cope with undertaking multivalent computational methodologies and tools. This work aimed to provide molecular-level insight into the enzyme-aided As degradation mechanism. AutoDock Vina, CABS-flex 2.0, and Desmond high-performance molecular dynamics simulation (MDS) were utilized in the current investigation to simulate multivalent molecular processes on two protein sets: arsenate reductase (ArsC) and laccase (LAC) corresponding arsenate (ART) and arsenite (AST), which served as model ligands to comprehend binding, conformational, and energy attributes. The structural configurations of both proteins exhibited variability in flexibility and structure framework within the range of 3.5-4.5 Å. The LAC-ART complex exhibited the lowest calculated binding affinity, measuring -5.82 ± 0.01 kcal/mol. Meanwhile, active site residues ILE-200 and HIS-206 were demonstrated to engage in H-bonding with the ART ligand. In contrast to ArsC, the ligand binding affinity of this bound complex was considerably greater. Additional validation of docked complexes was carried out by deploying Desmond MDS of 100 ns to capture protein and ligand conformation behavior. The system achieved stability during the 100 ns simulation run, as confirmed by the average P-L RMSD, which was ∼1 Å. As a preliminary test of the enzyme's ability to catalyze As species, corresponding computational insights might be advantageous for bridging gaps and regulatory consideration.

Oxidoreductase enzymes: Characteristics, applications, and challenges as a biocatalyst

Oxidoreductases are enzymes with distinctive characteristics that favor their use in different areas, such as agriculture, environmental management, medicine, and analytical chemistry. Among these enzymes, oxidases, dehydrogenases, peroxidases, and oxygenases are very interesting. Because their substrate diversity, they can be used in different biocatalytic processes by homogeneous and heterogeneous catalysis. Immobilization of these enzymes has favored their use in the solution of different biotechnological problems, with a notable increase in the study and optimization of this technology in the last years. In this review, the main structural and catalytical features of oxidoreductases, their substrate specificity, immobilization, and usage in biocatalytic processes, such as bioconversion, bioremediation, and biosensors obtainment, are presented.

Comparative homology of Pleurotus ostreatus laccase enzyme: Swiss model or Modeller?

Laccases stand out in the industrial context due to their versatile biotechnological applications. Although these enzymes are frequently investigated, currently, Pleurotus ostreatus laccase structural model is unknown. Therefore, this research aims to predict and validate a P. ostreatus laccase theoretical model by means of comparative homology. The laccase target's primary structure (AOM73725.1) was obtained from the NCBI database, the model was predicted from homologous structures obtained from the PDB (PDB-ID: 5A7E, 2HRG, 4JHU, 1GYC) using the Swiss-Model and Modeller, and was refined in GalaxyRefine. The models were validated using PROCHECK, VERIFY 3D, ERRAT, PROVE and QMEAN Z-score servers. Moreover, molecular docking between the laccase model (Lacc4MN) and ABTS was performed on AutoDock Vina. The models that were generated by the Modeller showed superior stereochemical and structural characteristics to those predicted by the Swiss Model. The refinement made it difficult to stabilize the copper atoms which are typical of laccases. The Lacc4MN model showed the interactions between the amino acids in the active site of the laccase and the copper atoms, thereby hinting the stabilization of the metal through electrostatic interactions with histidine and cysteine. The molecular docking between Lacc4MN and ABTS showed negative free energy and the formation of two hydrogen bonds involving the amino acids ASP 208 and GLY 268, and a Pi-sulfur bond between residue HIS 458 and ABTS, which demonstrates the typical catalytic functionality of laccases. Furthermore, the theoretical model Lacc4MN presented stereochemical and structural characteristics that allow its use in silico tests.Communicated by Ramaswamy H. Sarma.

Structural insights, biocatalytic characteristics, and application prospects of lignin-modifying enzymes for sustainable biotechnology-A review

Lignin modifying enzymes (LMEs) have gained widespread recognition in depolymerization of lignin polymers by oxidative cleavage. LMEs are a robust class of biocatalysts that include lignin peroxidase (LiP), manganese peroxidase (MnP), versatile peroxidase (VP), laccase (LAC), and dye-decolorizing peroxidase (DyP). Members of the LMEs family act on phenolic, non-phenolic substrates and have been widely researched for valorization of lignin, oxidative cleavage of xenobiotics and phenolics. LMEs implementation in the biotechnological and industrial sectors has sparked significant attention, although its potential future applications remain underexploited. To understand the mechanism of LMEs in sustainable pollution mitigation, several studies have been undertaken to assess the feasibility of LMEs in correlating to diverse pollutants for binding and intermolecular interactions at the molecular level. However, further investigation is required to fully comprehend the underlying mechanism. In this review we presented the key structural and functional features of LMEs, including the computational aspects, as well as the advanced applications in biotechnology and industrial research. Furthermore, concluding remarks and a look ahead, the use of LMEs coupled with computational frameworks, built upon artificial intelligence (AI) and machine learning (ML), has been emphasized as a recent milestone in environmental research.

A Genome-Wide Identification and Expression Pattern of LMCO Gene Family from Turnip (Brassica rapa L.) under Various Abiotic Stresses

Laccase-like multi-copper oxidases (LMCOs) are a group of enzymes involved in the oxidation of numerous substrates. Recently, these enzymes have become extremely popular due to their practical applications in various fields of biology. LMCOs generally oxidize various substrates by linking four-electron reduction of the final acceptor, molecular oxygen (O₂), to water. Multi-copper oxidases related to laccase are extensively distributed as multi-gene families in the genome sequences of higher plants. The current study thoroughly investigated the LMCO gene family (Br-Lac) and its expression pattern under various abiotic stresses in B. rapa L. A total of 18 Br-Lac gene family members located on five different chromosomes were identified. Phylogenetic analysis classified the documented Br-Lac genes into seven groups: Group-I (four genes), Group-II (nine genes), Group-III (eight genes), Group-IV (four genes), Group-V (six genes), and Group-VI and Group-VII (one gene each). The key features of gene structure and responsive motifs shared the utmost resemblance within the same groups. Additionally, a divergence study also assessed the evolutionary features of Br-Lac genes. The anticipated period of divergence ranged from 12.365 to 39.250 MYA (million years ago). We also identified the pivotal role of the 18 documented members of the LMCO (Br-lac) gene family using quantitative real-time qRT-PCR. Br-Lac-6, Br-Lac-7, Br-Lac-8, Br-Lac-16, Br-Lac-17, and Br-Lac-22 responded positively to abiotic stresses (i.e., drought, heat, and salinity). These findings set the stage for the functional diversity of the LMCO genes in B. rapa.

Assessing chemical hazard and unraveling binding affinity of priority pollutants to lignin modifying enzymes for environmental remediation

Lignin-modifying enzymes (LMEs) are impactful biocatalysts in environmental remediation applications. However, LMEs-assisted experimental degradation neglects the molecular basis of pollutant degradation. Furthermore, throughout the remediation process, the inherent hazards of environmental pollutants remain untapped for in-depth toxicological endpoints. In this investigation, a predictive toxicological framework and a computational framework adopting LMEs were employed to assess the hazards of Priority Pollutants (PP) and its possible LMEs-assisted catalytic screening. The potential hazardous outcomes of PP were assessed using Quantitative structure-activity relationship (QSARs)-based techniques including Toxtree, ECOSAR, and T.E.S.T. tools. Toxicological findings revealed positive outcomes in a multitude of endpoints for all PP. The PP compound 2,3,7,8-TCDD (dioxin) was found to exhibit the lowest concentration of aquatic toxicity implementing aquatic model systems; LC₅₀ as 0.01, 0.01, 0.04 (mg L^-1) for Fish (96 H), Daphnid (48 H), Green algae (96 H) respectively. T.E.S.T. results revealed that chloroform, and 2-chlorophenol both seem to be developmental toxicants. Subsequently, LMEs-assisted docking procedure was employed in predictive mitigation of PP. The docking approach as predicted degradation revealed the far lowest docking energy score for Versatile peroxidase (VP)- 2,3,7,8-TCDD docked complex with a binding energy of -9.2 (kcal mol^-1), involved PHE-46, PRO-139, PRO-141, ILE-148, LEU-165, HIS-169, LEU-228, MET-262, and MET-265 as key interacting amino acid residues. Second most ranked but lesser than VP, Lignin peroxidase (LiP)- 2,3,7,8-TCDD docked complex exhibited a rather lower binding affinity score (-8.8 kcal mol^-1). Predictive degradation screening employing comparative docking revealed varying binding affinities, portraying that each LMEs member has independent feasibility to bind PP as substrate. Predictive findings endorsed the hazardous nature of associated PP in a multitude of endpoints, which could be attenuated by undertaking LMEs as a predictive approach to protect the environment and implement it in regulatory considerations.

Experimental-theoretical study of laccase as a detoxifier of aflatoxins

We investigate laccase-mediated detoxification of aflatoxins, fungal carcinogenic food contaminants. Our experimental comparison between two aflatoxins with similar structures (AFB₁ and AFG₂) shows significant differences in laccase-mediated detoxification. A multi-scale modeling approach (Docking, Molecular Dynamics, and Density Functional Theory) identifies the highly substrate-specific changes required to improve laccase detoxifying performance. We employ a large-scale density functional theory-based approach, involving more than 7000 atoms, to identify the amino acid residues that determine the affinity of laccase for aflatoxins. From this study we conclude: (1) AFB₁ is more challenging to degrade, to the point of complete degradation stalling; (2) AFG₂ is easier to degrade by laccase due to its lack of side products and favorable binding dynamics; and (3) ample opportunities to optimize laccase for aflatoxin degradation exist, especially via mutations leading to π-π stacking. This study identifies a way to optimize laccase for aflatoxin bioremediation and, more generally, contributes to the research efforts aimed at rational enzyme optimization.

Comprehensive sequence and structure analysis of algal lipid catabolic enzyme Triacylglycerol lipase: an in silico study to vitalize the development of optimum engineered strains with high lipid productivity

Microalgae as an alternative renewable resource for biofuel production have captured much significance. Nonetheless, its economic viability is a field of major concern for researchers. Unraveling the lipid catabolic pathway and gaining insights into the sequence-structural features of its primary functioning enzyme, Triacylglycerol lipase, will impart valuable information to target microalgae for augmented lipid content. In the present study, a genome-wide comparative study on putative Triacylglycerol lipase (TAGL) enzyme from algal species belonging to varied phylogenetic lineages was performed. The comprehensive sequence analysis revealed that TAGL comprises of three distinct conserved domains, such as, Patatin, Class III Lipase, and Abhydro_lipase, and also confirmed the ubiquitous presence of GXSXG motif in the sequences analyzed. In the absence of a crystal structure of algal TAGL till date, we developed the first 3D model of patatin domain of TAGL from an oleaginous microalga, Phaedactylum tricornutum, employing homology modeling, docking and molecular dynamic simulations methods. The domain-substrate complex having the low-ranking docking score revealed the binding of palmitic acid to the TAGL patatin domain surface with strong hydrogen bond interactions. The simulation results implied that the substrate-complexed patatin domain and the free enzyme adopted a more stable conformation after 40 ns. This is the first ever attempt to provide in-silico insights into the structural and dynamical insights on catalytic mechanism of the TAGL patatin domain. Subsequently, these findings aided our understanding on their structural stability, folding mechanism and protein-substrate interactions, which could be further utilized to design site-specific mutagenic experiments for engineering microalgal strains.Communicated by Ramaswamy H. Sarma.

Investigating NAC Transcription Factor Role in Redox Homeostasis in Solanum lycopersicum L.: Bioinformatics, Physiological and Expression Analysis under Drought Stress

NAC transcription factors regulate stress-defence pathways and developmental processes in crop plants. However, their detailed functional characterization in tomatoes needs to be investigated comprehensively. In the present study, tomato hybrids subjected to 60 and 80 days of drought stress conditions showed a significant increase in membrane damage and reduced relative water, chlorophyll and proline content. However, hybrids viz., VRTH-16-3 and VRTH-17-68 showed superior growth under drought stress, as they were marked with low electrolytic leakage, enhanced relative water content, proline content and an enhanced activity of enzymatic antioxidants, along with the upregulation of NAC and other stress-defence pathway genes. Candidate gene(s) exhibiting maximum expression in all the hybrids under drought stress were subjected to detailed in silico characterization to provide significant insight into its structural and functional classification. The homology modelling and superimposition analysis of predicted tomato NAC protein showed that similar amino acid residues were involved in forming the conserved WKAT domain. DNA docking discovered that the SlNAC1 protein becomes activated and exerts a stress-defence response after the possible interaction of conserved DNA elements using Pro⁷², Asn⁷³, Trp⁸¹, Lys⁸², Ala⁸³, Thr⁸⁴, Gly⁸⁵, Thr⁸⁶ and Asp⁸⁷ residues. A protein-protein interaction analysis identified ten functional partners involved in the induction of stress-defence tolerance.

Genome-Wide Identification of Switchgrass Laccases Involved in Lignin Biosynthesis and Heavy-Metal Responses

Plant laccase genes belong to a multigene family, play key roles in lignin polymerization, and participate in the resistance of plants to biotic and abiotic stresses. Switchgrass is an important resource for forage and bioenergy production, yet information about the switchgrass laccase gene family is scarce. Using bioinformatic approaches, a genome-wide analysis of the laccase multigene family in switchgrass was carried out in this study. In total, 49 laccase genes (PvLac1 to PvLac49) were identified; these can be divided into five subclades, and 20 of them were identified as targets of miR397. The tandem and segmental duplication of laccase genes on Chr05 and Chr08 contributed to the expansion of the laccase family. The laccase proteins shared conserved signature sequences but displayed relatively low sequence similarity, indicating the potential functional diversity of switchgrass laccases. Switchgrass laccases exhibited distinct tissue/organ expression patterns, revealing that some laccases might be involved in the lignification process during stem development. All five of the laccase isoforms selected from different subclades responded to heavy metal. The immediate response of lignin-related laccases, as well as the delayed response of low-abundance laccases, to heavy-metal treatment shed light on the multiple roles of laccase isoforms in response to heavy-metal stress.

In silico analytical toolset for predictive degradation and toxicity of hazardous pollutants in water sources

Different phenolic compounds, including multimeric lignin derivatives in the β-O-4 form, are among the most prevalent compounds in wastewater, often generated from paper industries. Relatively small concentrations of lignin are hazardous to aquatic organisms and can trigger severe environmental hazards. Herein, we present a predictive toolset to insight the induced toxic hazards prediction, and their Lignin peroxidase (LiP)-assisted degradation mechanism of selected multimeric lignin model compounds. T.E.ST and Toxtree toolset were deployed for toxic hazards estimation in different endpoints. To minimize the concerning hazards, we screened multimeric compounds for binding affinity with LiP. The binding affinity was found to be significantly lower than the reference compound. An Extra precision (XP) Glide score of -6.796 kcal/mol was found for dimer (guaiacyl 4-O-5 guaiacyl) complex as lowest compared to reference compound (-4.007 kcal/mol). The active site residues ASP-153, HIP-226, VAL-227, ARG-244, GLU-215, 239, PHE-261 were identified as site-specific key binding AA residues actively involved with corresponding ligands, forming Hydrophobic, H-Bond, π-Stacking, π-π type interactions. The DESMOND-assisted molecular dynamics simulation's (MDS) trajectories of protein-ligand revealed the considerable binding behavior and attained stability and system equilibrium state. Such theoretical and predictive conclusions indicted the feasibility of LiP assisted sustainable mitigation of lignin-based compounds, and such could be used to protect the environment from the potential hazards posed by recognized similar pollutants.

Bioremediation of lignin derivatives and phenolics in wastewater with lignin modifying enzymes: Status, opportunities and challenges

Lignin modifying enzymes from fungi and bacteria are potential biocatalysts for sustainable mitigation of different potentially toxic pollutants in wastewater. Notably, the paper and pulp industry generates enormous amounts of wastewater containing high amounts of complex lignin-derived chlorinated phenolics and sulfonated pollutants. The presence of these compounds in wastewater is a critical issue from environmental and toxicological perspectives. Some chloro-phenols are harmful to the environment and human health, as they exert carcinogenic, mutagenic, cytotoxic, and endocrine-disrupting effects. In order to address these most urgent concerns, the use of oxidative lignin modifying enzymes for bioremediation has come into focus. These enzymes catalyze modification of phenolic and non-phenolic lignin-derived substances, and include laccase and a range of peroxidases, specifically lignin peroxidase (LiP), manganese peroxidase (MnP), versatile peroxidase (VP), and dye-decolorizing peroxidase (DyP). In this review, we explore the key pollutant-generating steps in paper and pulp processing, summarize the most recently reported toxicological effects of industrial lignin-derived phenolic compounds, especially chlorinated phenolic pollutants, and outline bioremediation approaches for pollutant mitigation in wastewater from this industry, emphasizing the oxidative catalytic potential of oxidative lignin modifying enzymes in this regard. We highlight other emerging biotechnical approaches, including phytobioremediation, bioaugmentation, Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-based technology, protein engineering, and degradation pathways prediction, that are currently gathering momentum for the mitigation of wastewater pollutants. Finally, we address current research needs and options for maximizing sustainable biobased and biocatalytic degradation of toxic industrial wastewater pollutants.

A Brief History of Colour, the Environmental Impact of Synthetic Dyes and Removal by Using Laccases

The history of colour is fascinating from a social and artistic viewpoint because it shows the way; use; and importance acquired. The use of colours date back to the Stone Age (the first news of cave paintings); colour has contributed to the social and symbolic development of civilizations. Colour has been associated with hierarchy; power and leadership in some of them. The advent of synthetic dyes has revolutionized the colour industry; and due to their low cost; their use has spread to different industrial sectors. Although the percentage of coloured wastewater discharged by the textile; food; pharmaceutical; cosmetic; and paper industries; among other productive areas; are unknown; the toxic effect and ecological implications of this discharged into water bodies are harmful. This review briefly shows the social and artistic history surrounding the discovery and use of natural and synthetic dyes. We summarise the environmental impact caused by the discharge of untreated or poorly treated coloured wastewater to water bodies; which has led to physical; chemical and biological treatments to reduce the colour units so as important physicochemical parameters. We also focus on laccase utility (EC 1.10.3.2), for discolouration enzymatic treatment of coloured wastewater, before its discharge into water bodies. Laccases (p-diphenol: oxidoreductase dioxide) are multicopper oxidoreductase enzymes widely distributed in plants, insects, bacteria, and fungi. Fungal laccases have employed for wastewater colour removal due to their high redox potential. This review includes an analysis of the stability of laccases, the factors that influence production at high scales to achieve discolouration of high volumes of contaminated wastewater, the biotechnological impact of laccases, and the degradation routes that some dyes may follow when using the laccase for colour removal.

Recombinant laccase rPOXA 1B real-time, accelerated and molecular dynamics stability study

Background

Laccases (EC 1.10.3.2) are multi-copper oxidoreductases with great biotechnological importance due to their high oxidative potential and utility for removing synthetic dyes, oxidizing phenolic compounds, and degrading pesticides, among others.

Methods

A real-time stability study (RTS) was conducted for a year, by using enzyme concentrates from 3 batches (L1, L3, and L4). For which, five temperatures 243.15, 277.15, 298.15, 303.15, 308.15, and 313.15 K were assayed. Using RTS data and the Arrhenius equation, we calculated the rPOXA 1B accelerated stability (AS). Molecular dynamics (MD) computational study results were very close to those obtained experimentally at four different temperatures 241, 278, 298, and 314 K.

Results

In the RTS, 101.16, 115.81, 75.23, 46.09, 5.81, and 4.83% of the relative enzyme activity were recovered, at respective assayed temperatures. AS study, showed that rPOXA 1B is stable at 240.98 ± 5.38, 277.40 ± 1.32 or 297.53 ± 3.88 K; with t_1/2 values of 230.8, 46.2, and 12.6 months, respectively. Kinetic and thermodynamic parameters supported the high stability of rPOXA 1B, with an E_d value of 41.40 KJ mol^− 1, a low variation of K_M and V_max, at 240.98 ± 5.38, and 297.53 ± 3.88 K, and ∆G values showing deactivation reaction does not occur. The MD indicates that fluctuations in loop, coils or loops with hydrophilic or intermediate polarity amino acids as well as in some residues of POXA 1B 3D structure, increases with temperature; changing from three fluctuating residues at 278 K to six residues at 298 K, and nine residues at 314 K.

Conclusions

Laccase rPOXA 1B demonstrated experimentally and computationally to be a stable enzyme, with t_1/2 of 230.8, 46.2 or 12.6 months, if it is preserved impure without preservatives at temperatures of 240.98 ± 5.38, 277.40 ± 1.32 or 297.53 ± 3.88 K respectively; this study could be of great utility for large scale producers.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12896-021-00698-3.

In silico exploration of lignin peroxidase for unraveling the degradation mechanism employing lignin model compounds

Lignin peroxidase is a heme-containing biocatalyst, well-known for its diverse applications in the fields from environmental chemistry to biotechnology. LiP-mediated oxidative catalysis is H₂O₂-dependent, and can oxidize phenolic, and non-phenolic substrates by oxidative cleavage of the C-C and C-O bonds of lignin. In contrast to fungi-derived LiP, the binding affinity of bacterial-derived LiP to lignin at the molecular level is poorly known to date. Tremendous wet-lab studies have been unveiled that provide degradation and biotransformation information on kraft lignin, whilst studies on the completely transformed compounds and the degradation of each transformed compounds simultaneously during degradation are scarce. To gain an understanding of the degradation process using docking, and MDS based studies, we assessed the binding affinity of selected lignin model compounds with bacterial origin LiP and validated such docked complexes exploiting 30 ns molecular dynamics simulations. We selected and picked a total of 12 lignin model compounds for molecular modeling analysis, namely two chlorinated lignin model compounds (monomer) (2-chlorosyringaldehyde and 5-chlorovanillin), eight standard lignin model compounds (veratryl alcohol, syringyl alcohol, sinapyl alcohol, methyl hydroquinone, guaiacol, coniferyl alcohol, catechol, and 4-methoxy phenol), while, two 4-O-5, and β-O-4 linkage-based multimeric model compounds (dimer: 2-methoxy-6-(2-methoxy-4-methylphenoxy)-4-methylphenol; trimer: syringyl β-O-4 syringyl β-O-4 sinapyl alcohol). Far more specific binding residues were observed from XP-Glide docking, as TYR, HIP (protonated histidine), PHE, VAL, ASP, THR, LYS and GLN. The binding affinity was confirmed by the Gibbs free energy or binding energy (ΔG) score; furthermore, it is found that the maximum binding energy seems to be observed for 4-methoxyphenol with a Glide score of -3.438 with Pi-Pi stacking and H-bond type bonding interactions, whilst the lowest XP Gscore as -8.136 with Pi-Pi stacking and H-bond (side chain) type bonding interactions were found for the trimer model compound. The docked complexes were further evaluated for deep rigorous structural and functional fluctuation analyses through high-performance molecular dynamics simulations-DESMOND, after a post simulation run of 30 ns. The RMSD trajectory analyses of the protein-ligands were found to be in the equilibrium state at the end of simulation run for multimeric lignin model compounds. In addition, ionic ligand-protein interaction occurs among chlorinated compounds, while hydrophobic and H-bond contacts have frequently been observed in all lignin-model compounds. The findings herein demonstrate that bacterial LiP can effectively catalyze multiple lignin model compounds, and it might further be used as an effective tool for sustainable mitigation of diverse environmental contaminants.

An effective computational-screening strategy for simultaneously improving both catalytic activity and thermostability of α-l-rhamnosidase

Catalytic efficiency and thermostability are the two most important characteristics of enzymes. However, it is always tough to improve both catalytic efficiency and thermostability of enzymes simultaneously. In the present study, a computational strategy with double-screening steps was proposed to simultaneously improve both catalysis efficiency and thermostability of enzymes; and a fungal α-l-rhamnosidase was used to validate the strategy. As the result, by molecular docking and sequence alignment analysis within the binding pocket, seven mutant candidates were predicted with better catalytic efficiency. By energy variety analysis, A355N, S356Y, and D525N among the seven mutant candidates were predicted with better thermostability. The expression and characterization results showed the mutant D525N had significant improvements in both enzyme activity and thermostability. Molecular dynamics simulations indicated that the mutations located within the 5 Å range of the catalytic domain, which could improve root mean squared deviation, electrostatic, Van der Waal interaction, and polar salvation values, and formed water bridge between the substrate and the enzyme. The study indicated that the computational strategy based on the binding energy, conservation degree and mutation energy analyses was effective to develop enzymes with better catalysis and thermostability, providing practical approach for developing industrial enzymes.

A gene that underwent adaptive evolution, LAC2 (LACCASE), in Populus euphratica improves drought tolerance by improving water transport capacity

Drought severely limits plant development and growth; accordingly, plants have evolved strategies to prevent water loss and adapt to water deficit conditions. However, experimental cases that corroborate these evolutionary processes are limited. The LACCASEs (LACs) family is involved in various plant development and growth processes. Here, we performed an evolutionary analysis of LACs from Populus euphratica and characterized the functions of LACs in Arabidopsis and poplar. The results showed that in PeuLACs, multiple gene duplications led to apparent functional redundancy as the result of various selective pressures. Among them, PeuLAC2 underwent strong positive selection. Heterologous expression analyses showed that the overexpression of PeuLAC2 alters the xylem structure of plants, including thickening the secondary cell wall (SCW) and increasing the fiber cell length and stem tensile strength. Altogether, these changes improve the water transport capacity of plants. The analysis of the physiological experimental results showed that PeuLAC2-OE lines exhibited a stronger antioxidant response and greater drought tolerance than WT. Three genes screened by transcriptome analysis, NAC025, BG1, and UGT, that are associated with SCW synthesis and drought stress were all upregulated in the PeuLAC2-OE lines, implying that the overexpression of PeuLAC2 thickened the SCW, improved the water transport capacity of the plant, and further enhanced its drought tolerance. Our study highlights that genes that have undergone adaptive evolution may participate in the development of adaptive traits in P. euphratica and that PeuLAC2 could be a candidate gene for molecular genetic breeding in trees.

The transport of a charged peptide through carbon nanotubes under an external electric field: a molecular dynamics simulation

The study of interactions between biomolecules and carbon nanotubes (CNTs) is of great importance in CNT-based drug delivery systems and biomedical devices. In this work, the transport of polyarginine (R8) peptide through CNTs under an external electric field was investigated via all-atom molecular dynamics (AAMD) simulation. It was found that the electric field can assist the R8 peptide to overcome the resistance and make the transport smooth. Moreover, the efficiency of transport was improved with the increasing intensity of the electric field in a suitable range. In addition, we also investigated the effects of different types of CNTs on the transport of the R8 peptide and found that the single-walled carbon nanotube (SWCNT) was more suitable for transporting the R8 peptide than the double-walled carbon nanotube (DWCNT) due to its lower energy barrier to the R8 peptide. All these findings shed light on the role of the electric field on the transport of the R8 peptide through CNTs and also gave some valuable insights into the effects of CNT types on the transport process of the peptide.

The role of electric field and types of carbon nanotube influencing the delivery process of peptide through CNTs were studied via all-atom molecular dynamics simulation.

Genes Identification, Molecular Docking and Dynamics Simulation Analysis of Laccases from Amylostereum areolatum Provides Molecular Basis of Laccase Bound to Lignin

An obligate mutualistic relationship exists between the fungus Amylostereum areolatum and woodwasp Sirex noctilio. The fungus digests lignin in the host pine, providing essential nutrients for the growing woodwasp larvae. However, the functional properties of this symbiosis are poorly described. In this study, we identified, cloned, and characterized 14 laccase genes from A. areolatum. These genes encoded proteins of 508 to 529 amino acids and contained three typical copper-oxidase domains, necessary to confer laccase activity. Besides, we performed molecular docking and dynamics simulation of the laccase proteins in complex with lignin compounds (monomers, dimers, trimers, and tetramers). AaLac2, AaLac3, AaLac6, AaLac8, and AaLac10 were found that had low binding energies with all lignin model compounds tested and three of them could maintain stability when binding to these compounds. Among these complexes, amino acid residues ALA, GLN, LEU, PHE, PRO, and SER were commonly present. Our study reveals the molecular basis of A. areolatum laccases interacting with lignin, which is essential for understanding how the fungus provides nutrients to S. noctilio. These findings might also provide guidance for the control of S. noctilio by informing the design of enzyme mutants that could reduce the efficiency of lignin degradation.

Role of laccase gene in wheat NILs differing at QTL-Fhb1 for resistance against Fusarium head blight

Fusarium head blight (FHB), caused mainly by Fusarium graminearum (Fg), is one of the most severe diseases of wheat. It affects grain yield and quality due to mycotoxin contamination, which is harmful for both human and livestock consumption. Cell wall lignification, following pathogen invasion, is one of the innate defense responses. Plant laccases are known to lignify the secondary cell walls. A metabolo-genomics study identified laccase as one of the candidate genes in QTL-Fhb1 of wheat NILs derived from Sumai 3*5/Thatcher cross. Based on phylogenetics, it was named as TaLAC4. Real-time qPCR revealed a strongly induced expression of TaLAC4 in NIL-R. The VIGS based transient silencing of TaLAC4 in NIL-R resulted in an increased susceptibility leading to Fg spread within the entire spike in 15dpi, contrasting to non-silenced where the infection was limited to inoculated spikelets. Histopathology revealed thickened cell walls, mainly due to G-lignin, in non-silenced NIL-R, relative to silenced, in conjunction with higher total lignin content. Metabolic profiling of TaLAC4 silenced NILs identified the accumulation of several precursor metabolites higher in abundances upstream TaLAC4. These results confirm that the resistance function of TaLAC4 in NIL-R is due to pathogen-induced lignification of secondary cell walls in the rachis.

Modeling to Understand Plant Protein Structure-Function Relationships-Implications for Seed Storage Proteins

Proteins are among the most important molecules on Earth. Their structure and aggregation behavior are key to their functionality in living organisms and in protein-rich products. Innovations, such as increased computer size and power, together with novel simulation tools have improved our understanding of protein structure-function relationships. This review focuses on various proteins present in plants and modeling tools that can be applied to better understand protein structures and their relationship to functionality, with particular emphasis on plant storage proteins. Modeling of plant proteins is increasing, but less than 9% of deposits in the Research Collaboratory for Structural Bioinformatics Protein Data Bank come from plant proteins. Although, similar tools are applied as in other proteins, modeling of plant proteins is lagging behind and innovative methods are rarely used. Molecular dynamics and molecular docking are commonly used to evaluate differences in forms or mutants, and the impact on functionality. Modeling tools have also been used to describe the photosynthetic machinery and its electron transfer reactions. Storage proteins, especially in large and intrinsically disordered prolamins and glutelins, have been significantly less well-described using modeling. These proteins aggregate during processing and form large polymers that correlate with functionality. The resulting structure-function relationships are important for processed storage proteins, so modeling and simulation studies, using up-to-date models, algorithms, and computer tools are essential for obtaining a better understanding of these relationships.

Laccases: structure, function, and potential application in water bioremediation

The global rise in urbanization and industrial activity has led to the production and incorporation of foreign contaminant molecules into ecosystems, distorting them and impacting human and animal health. Physical, chemical, and biological strategies have been adopted to eliminate these contaminants from water bodies under anthropogenic stress. Biotechnological processes involving microorganisms and enzymes have been used for this purpose; specifically, laccases, which are broad spectrum biocatalysts, have been used to degrade several compounds, such as those that can be found in the effluents from industries and hospitals. Laccases have shown high potential in the biotransformation of diverse pollutants using crude enzyme extracts or free enzymes. However, their application in bioremediation and water treatment at a large scale is limited by the complex composition and high salt concentration and pH values of contaminated media that affect protein stability, recovery and recycling. These issues are also associated with operational problems and the necessity of large-scale production of laccase. Hence, more knowledge on the molecular characteristics of water bodies is required to identify and develop new laccases that can be used under complex conditions and to develop novel strategies and processes to achieve their efficient application in treating contaminated water. Recently, stability, efficiency, separation and reuse issues have been overcome by the immobilization of enzymes and development of novel biocatalytic materials. This review provides recent information on laccases from different sources, their structures and biochemical properties, mechanisms of action, and application in the bioremediation and biotransformation of contaminant molecules in water. Moreover, we discuss a series of improvements that have been attempted for better organic solvent tolerance, thermo-tolerance, and operational stability of laccases, as per process requirements.

Effect of surfactants on the interaction of phenol with laccase: Molecular docking and molecular dynamics simulation studies

Some surfactants can enhance the removal of phenol by laccase (Lac) in various industrial effluents. Their behavior and function in the biodegradation of phenolic wastewater have been experimentally reported by many researchers, but the underlying molecular mechanism is still unclear. Therefore, the interaction mechanisms of phenol with Lac from Trametes versicolor were investigated in the presence or absence of Triton X-100 (TX100) or rhamnolipid (RL) by molecular docking and molecular dynamics (MD) simulations. The results indicate that phenol contacts with an active site of Lac by hydrogen bonds (HBs) and van der Waals (vdW) interactions in aqueous solution for maintaining its stability. The presence of TX100 or RL results in the significant changes of enzymatic conformations. Meanwhile, the hydrophobic parts of surfactants contact with the outside surface of Lac. These changes lead to the decrease of binding energy between phenol and Lac. The migration behavior of water molecules within hydration shell is also inevitably affected. Therefore, the amphipathic TX100 or RL may influence the phenol degradation ability of Lac by modulating their interactions and water environment. This study offers molecular level of understanding on the function of surfactants in biosystem.

The penetration of a charged peptide across a membrane under an external electric field: a coarse-grained molecular dynamics simulation

The processes of single polyarginine (R8) peptide penetration through planar and vesicle membranes under an external electric field are simulated via a coarse-grained molecular dynamics (CGMD) simulation.

Comparative modeling and molecular docking analysis of white, brown and soft rot fungal laccases using lignin model compounds for understanding the structural and functional properties of laccases

Extrinsic catalytic properties of laccase enable it to oxidize a wide range of aromatic (phenolic and non-phenolic) compounds which makes it commercially an important enzyme. In this study, we have extensively compared and analyzed the physico-chemical, structural and functional properties of white, brown and soft rot fungal laccases using standard protein analysis software. We have computationally predicted the three-dimensional comparative models of these laccases and later performed the molecular docking studies using the lignin model compounds. We also report a customizable rapid and reliable protein modelling and docking pipeline for developing structurally and functionally stable protein structures. We have observed that soft rot fungal laccases exhibited comparatively higher structural variation (higher random coil) when compared to brown and white rot fungal laccases. White and brown rot fungal laccase sequences exhibited higher similarity for conserved domains of Trametes versicolor laccase, whereas soft rot fungal laccases shared higher similarity towards conserved domains of Melanocarpus albomyces laccase. Results obtained from molecular docking studies showed that aminoacids PRO, PHE, LEU, LYS and GLN were commonly found to interact with the ligands. We have also observed that white and brown rot fungal laccases showed similar docking patterns (topologically monomer, dimer and trimer bind at same pocket location and tetramer binds at another pocket location) when compared to soft rot fungal laccases. Finally, the binding efficiencies of white and brown rot fungal laccases with lignin model compounds were higher compared to the soft rot fungi. These findings can be further applied in developing genetically efficient laccases which can be applied in growing biofuel and bioremediation industries.

Interactions of carbon nanotubes and/or graphene with manganese peroxidase during biodegradation of endocrine disruptors and triclosan

Molecular-level biodegradation processes of bisphenol A (BPA), nonylphenol (NP) and triclosan (TCS) mediated by manganese peroxidase (MnP) were investigated with and without single-walled carbon nanotube (SWCNT) and/or graphene (GRA). Although the incorporation of SWCNT, GRA or their combination (SWCNT+GRA) did not break up the complexes composed of manganese peroxidase (MnP) and these substrates, they had different effects on the native contacts between the substrates and MnP. GRA tended to decrease the overall stability of the binding between MnP and its substrates. SWCNT or SWCNT+GRA generally had a minor impact on the mean binding energy between MnP and its substrates. We detected some sensitive residues from MnP that were dramatically disturbed by the GRA, SWCNT or SWCNT+GRA. Nanomaterials changed the number and behavior of water molecules adjacent to both MnP and its substrates, which was not due to the destruction of H-bond network formed by sensitive regions and water molecules. The present results are useful for understanding the molecular basis of pollutant biodegradation affected by the nanomaterials in the environment, and are also helpful in assessing the risks of these materials to the environment.

A Kepler Workflow Tool for Reproducible AMBER GPU Molecular Dynamics

With the drive toward high throughput molecular dynamics (MD) simulations involving ever-greater numbers of simulation replicates run for longer, biologically relevant timescales (microseconds), the need for improved computational methods that facilitate fully automated MD workflows gains more importance. Here we report the development of an automated workflow tool to perform AMBER GPU MD simulations. Our workflow tool capitalizes on the capabilities of the Kepler platform to deliver a flexible, intuitive, and user-friendly environment and the AMBER GPU code for a robust and high-performance simulation engine. Additionally, the workflow tool reduces user input time by automating repetitive processes and facilitates access to GPU clusters, whose high-performance processing power makes simulations of large numerical scale possible. The presented workflow tool facilitates the management and deployment of large sets of MD simulations on heterogeneous computing resources. The workflow tool also performs systematic analysis on the simulation outputs and enhances simulation reproducibility, execution scalability, and MD method development including benchmarking and validation.

Gel-Based Purification and Biochemical Study of Laccase Isozymes from Ganoderma sp. and Its Role in Enhanced Cotton Callogenesis

Basidiomycetous fungi, Ganoderma lucidum MDU-7 and Ganoderma sp. kk-02 secreted multiple laccase isozymes under diverse growth condition. Aromatic compounds and metal salts were also found to regulate the differential expression of laccase isozymes from both the Ganoderma sp. Laccase isozymes induced in the presence of copper from G. lucidum MDU-7 were purified by gel-based (native-PAGE) purification method. The purity of laccase isozymes was checked by zymogram and SDS-PAGE. The SDS-PAGE of purified proteins confirmed the multimeric nature of laccase isozymes. The molecular mass of isozymes was found to be in the range of 40–66 kDa. Further, the purified laccase isozymes and their peptides were confirmed with the help of MALDI-TOF peptide fingerprinting. The biochemical characterization of laccase isozymes viz. Glac L2, Glac L3, Glac L4, and Glac L5 have shown the optimum temperature in the range of 30°–45°C and pH 3.0. The K_m values of all the laccase isozymes determined for guaiacol were (96–281 μM), ABTS (15–83 μM) and O-tolidine (78–724 μM). Further, laccase isozymes from G. lucidum whole genome were studied using bioinformatics tools. The molecular modeling and docking of laccase isozymes with different substrates showed a significant binding affinity, which further validates our experimental results. Interestingly, copper induced laccase of 40 U/ml in culture medium was found to significantly induce cotton callogenesis. Interestingly, all the laccase isozymes were found to have an antioxidative role and therefore capable in free radicals scavenging during callogenesis. This is the first detailed study on the biochemical characterization of all the laccase isozymes purified by a gel-based novel method.

Laccase-13 Regulates Seed Setting Rate by Affecting Hydrogen Peroxide Dynamics and Mitochondrial Integrity in Rice

Seed setting rate is one of the most important components of rice grain yield. To date, only several genes regulating setting rate have been identified in plant. In this study, we showed that laccase-13 (OsLAC13), a member of laccase family genes which are known for their roles in modulating phenylpropanoid pathway and secondary lignification in cell wall, exerts a regulatory function in rice seed setting rate. OsLAC13 expressed in anthers and promotes hydrogen peroxide production both in vitro and in the filaments and anther connectives. Knock-out of OsLAC13 showed significantly increased seed setting rate, while overexpression of this gene exhibited induced mitochondrial damage and suppressed sugar transportation in anthers, which in turn affected seed setting rate. OsLAC13 also induced H₂O₂ production and mitochondrial damage in the root tip cells which caused the lethal phenotype. We also showed that high abundant of OsmiR397, the suppressor of OsLAC13 mRNA, increased the seed setting rate of rice plants, and restrains H₂O₂ accumulation in roots during oxidative stress. Our results suggested a novel regulatory role of OsLAC13 gene in regulating seed setting rate by affecting H₂O₂ dynamics and mitochondrial integrity in rice.

Molecular docking and molecular dynamics simulation analyses of urea with ammoniated and ammoxidized lignin

Ammoniated lignin, prepared through the Mannich reaction of lignin, has more advantages as a slow-release carrier of urea molecules than ammoxidized lignin and lignin. The advantages of the ammoniated lignin include its amine groups added and its high molecular mass kept as similar as that of lignin. Three organic molecules including guaiacyl, 2-hydroxybenzylamine and 5-carbamoylpentanoic acid are monomers respectively in lignin, ammoniated lignin and ammoxidized lignin. We studied the difference between the interactions of lignin, ammoniated lignin and ammoxidized lignin with respect to urea, based on radial distribution functions (RDFs) results from molecular dynamics (MD) simulations. Glass transition temperature (T_g) and solubility parameter (δ) of ammoniated and ammoxidized lignin have been calculated by MD simulations in the constant-temperature and constant-pressure ensemble (NPT). Molecular docking results showed the interaction sites of the urea onto the ammoniated and ammoxidized lignin and three different interaction modes were identified. Root mean square deviation (RMSD) values could indicate the mobilities of the urea molecule affected by the three different interaction modes. A series of MD simulations in the constant-temperature and constant-volume ensemble (NVT) helped us to calculate the diffusivity of urea which was affected by the content of urea in ammoniated and ammoxidized lignin.

Knockdown of a laccase in Populus deltoides confers altered cell wall chemistry and increased sugar release

Plant laccases are thought to function in the oxidation of monolignols which leads to higher order lignin formation. Only a hand-full of laccases in plants have been functionally evaluated, and as such little is known about the breadth of their impact on cell wall chemistry or structure. Here, we describe a previously uncharacterized laccase from Populus, encoded by locus Potri.008G064000, whose reduced expression resulted in transgenic Populus trees with changes in syringyl/guaiacyl ratios as well as altered sugar release phenotypes. These phenotypes are consistent with plant biomass exhibiting reduced recalcitrance. Interestingly, the transgene effect on recalcitrance is dependent on a mild pretreatment prior to chemical extraction of sugars. Metabolite profiling suggests the transgene modulates phenolics that are associated with the cell wall structure. We propose that this particular laccase has a range of functions related to oxidation of phenolics and conjugation of flavonoids that interact with lignin in the cell wall.

Conversion of lignin into value-added materials and chemicals via laccase-assisted copolymerization

With today's environmental concerns and the diminishing supply of the world's petroleum-based chemicals and materials, much focus has been directed toward alternative sources. Woody biomass presents a promising option due to its sheer abundance, renewability, and biodegradability. Lignin, a highly irregular polyphenolic compound, is one of the major chemical constituents of woody biomass and is the second most abundant biopolymer on Earth, surpassed only by cellulose. The pulp and paper and cellulosic ethanol industries produce lignin on the scale of millions of tons each year as a by-product. Traditionally, lignin has been viewed as a waste material and burned as an inefficient fuel. However, in recent decades, research has focused on more economical ways to convert lignin into value-added commodities, such as biofuels, biomaterials, and biochemicals, thus developing and strengthening the concept of fully integrated biorefineries. Owing to the phenolic structure of lignin, it is possible to enzymatically graft molecules onto its surface using laccases (benzenediol:oxygen oxidoreductases, EC 1.10.3.2) to create exciting novel biomaterials. These environmentally friendly enzymes use oxygen as their only co-substrate and produce water as their sole by-product, and have thus found great industrial application. This mini-review highlights recent advances in the field of laccase-facilitated functionalization of lignin as well as promising future directions for lignin-based polymers.

Two Decades of Laccases: Advancing Sustainability in the Chemical Industry

Given the current state of environmental affairs and that our future on this planet as we know it is in jeopardy, research and development into greener and more sustainable technologies within the chemical and forest products industries is at its peak. Given the global scale of these industries, the need for environmentally benign practices is propelling new green processes. These challenges are also impacting academic research and our reagents of interest are laccases. These enzymes are employed in a variety of biotechnological applications due to their native function as catalytic oxidants. They are about as green as it gets when it comes to chemical processes, requiring O₂ as their only co-substrate and producing H₂ O as the sole by-product. The following account will review our twenty year journey on the use of these enzymes within our research group, from their initial use in biobleaching of kraft pulps and for fiber modification within the pulp and paper industry, to their current application as green catalytic oxidants in the field of synthetic organic chemistry.

Lignin engineering through laccase modification: a promising field for energy plant improvement

Laccase (p-diphenol:dioxygen oxidoreductase, EC 1.10.3.2) is a member of the multicopper oxidases and catalyzes the one-electron oxidation of a wide range of substrates, coupled with the reduction of oxygen to water. It is widely distributed in bacteria, fungi, plants and insects. Laccases are encoded by multigene family, and have been characterized mostly from fungi till now, with abundant industrial applications in pulp and paper, textile, food industries, organic synthesis, bioremediation and nanobiotechnology, while limited researches have been performed in plants, and no application has been reported. Plant laccases share the common molecular architecture and reaction mechanism with fungal ones, despite of difference in redox potential and pH optima. Plant laccases are implicated in lignin biosynthesis since genetic evidence was derived from the Arabidopsis LAC4 and LAC17. Manipulation of plant laccases has been considered as a promising and innovative strategy in plant biomass engineering for desirable lignin content and/or composition, since lignin is the major recalcitrant component to saccharification in biofuel production from lignocellulose, and therefore directly limits the fermentation yields. Moreover, plant laccases have been reported to be involved in wound healing, maintenance of cell wall structure and integrity, and plant responses to environmental stresses. Here, we summarize the properties and functions of plant laccase, and discuss the potential of biotechnological application, thus providing a new insight into plant laccase, an old enzyme with a promising beginning in lignocellulose biofuel production.