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

Laccase, widely distributed in bacteria, fungi, and plants, catalyzes the oxidation of wide range of compounds. With regards to one of the important physiological functions, plant laccases are considered to catalyze lignin biosynthesis while fungal laccases are considered for lignin degradation. The present study was undertaken to explain this dual function of laccases using in-silico molecular docking and dynamics simulation approaches. Modeling and superimposition analyses of one each representative of plant and fungal laccases, namely, Populus trichocarpa and Trametes versicolor, respectively, revealed low level of similarity in the folding of two laccases at 3D levels. Docking analyses revealed significantly higher binding efficiency for lignin model compounds, in proportion to their size, for fungal laccase as compared to that of plant laccase. Residues interacting with the model compounds at the respective enzyme active sites were found to be in conformity with their role in lignin biosynthesis and degradation. Molecular dynamics simulation analyses for the stability of docked complexes of plant and fungal laccases with lignin model compounds revealed that tetrameric lignin model compound remains attached to the active site of fungal laccase throughout the simulation period, while it protrudes outwards from the active site of plant laccase. Stability of these complexes was further analyzed on the basis of binding energy which revealed significantly higher stability of fungal laccase with tetrameric compound than that of plant. The overall data suggested a situation favorable for the degradation of lignin polymer by fungal laccase while its synthesis by plant laccase.