key: cord-0077064-fgynhz35
authors: Paul, Parag Kumar; Al Azad, Salauddin; Rahman, Mohammad Habibur; Farjana, Mithila; Uddin, Muhammad Ramiz; Dey, Dipta; Mahmud, Shafi; Ema, Tanzila Ismail; Biswas, Partha; Anjum, Maliha; Akhi, Ozifatun Jannat; Ahmed, Shahlaa Zernaz
title: Catabolic profiling of selective enzymes in the saccharification of non-food lignocellulose parts of biomass into functional edible sugars and bioenergy: An in silico bioprospecting
date: 2022-01-14
journal: J Adv Vet Anim Res
DOI: 10.5455/javar.2022.i565
sha: 93f5ac6dc2f2e303e867d4608b653b27fe6ed0fd
doc_id: 77064
cord_uid: fgynhz35

OBJECTIVES: The research aims to analyze the catabolic strength of different hydrolytic enzymes in assessing the biological conversion potential of lignocellulose parts of agricultural biomass wastes into functional edible sugars and biofuels. MATERIALS AND METHODS: The enzymes’ hydrolytic properties—versatile peroxidase, manganese peroxidase, and lignin peroxidase were used to identify their complexing strength with the lignin substrate, whereas endoglucanase cel12A, acidocaldarius cellulase, and Melanocarpus albomyces endoglucanase were tested on the cellulose gel substrate. Because the biodegradation properties are heavily influenced by the “enzyme-substrate complexing energy level,” proper molecular optimization and energy minimization of the enzymes and substrates were carried out, as well as the identification of the enzyme’s active sites prior to complexing.comprehensive molecular dynamic simulation was run to study their—alpha carbon, root-mean-square deviation (Å), molecular surface area (Å2), root-mean-square fluctuation (Å), radius of gyration (nm), hydrogen bonds with hydrophobic interactions, and solvent accessible surface area (Å2) values for 50 ns. The simulated data mining was conducted using advanced programming algorithms to establish the final enzyme-substrate complexing strength in binding and catalysis. RESULTS: Among the lignin-degrading enzymes, versatile peroxidase shows promising catalytic activity with the best docking pose and significant values in all the dynamic simulation parameters. Similarly, Melanocarpus albomyces endoglucanase shows the best activity in all aspects of molecular docking and dynamics among the cellulose-degrading enzymes. CONCLUSION: The lignin content of biomass wastes can be degraded into cellulose and hemicellulose using lignin-degrading enzymes. The cellulose can be further degraded into glucose and xylose sugars following the cellulose-degrading enzyme activity. These sugars can be further degraded into biofuel through anaerobic fermentation. Systematic bioconversion of the lignocellulosic components can ensure sustainable biomass management, creating an alternative food and energy source for human beings to face the challenges of global hunger where the enzymes can pave the way.

The world's population has quadrupled over the last century [1] . The world population was 1.8 billion in 1915, according to the US census [2] , and it is now 7.3 billion, according to the US census [3] .This number will continue to grow to 10 billion as soon as 2050 [1] . The global food crisis has long been a significant threat, and this increasing population further threatens food insurance for the near future [3] . Even now, 8.9% of the world population is hungry and malnourished, and more than half of the African population will suffer from food shortages by 2030 [3] . Therefore, the traditional means of food security for the ever-increasing population is almost questionable, indicating the need for alternative food and nutrient production sources. Lignocellulosic biomass degradation can be an alternative and renewable energy source of nutrients and a food production process [4] . It is a plant-based material that is not used for food or energy sources. However, because of its chemical nature, mainly composed of carbohydrate polymers (lignin) and polysaccharides (celluloses and hemicelluloses) [4] , functional, simple sugars like glucose and xylose can be biosynthesized from lignocellulosic biomass. Glucose can be directly consumed as a nutrient and as energy. It could also be used as a starting or intermediate material for various industrial and medical needs, such as candy, beer, glucose syrup, sweetener, thickener, moisture-retaining agent, fondant, and certain canned and premade baked goods [5] .

Lignocellulosic biomass has typically been utilized to produce bioethanol, a novel and renewable energy resource for bioethanol production as a sustainable substitute for fossil fuels and toxic chemicals [6] . However, by utilizing the catabolic activities of several enzymes, this biomass can be converted into glucose and xylose [7] as the edible energy sources of humans. First, lignocellulosic biomass can be enzymatically hydrolyzed to expose the unstable cellulose and hemicellulose portions. Second, cellulose can be further catabolized to produce glucose and xylose [8] . The total process is termed an enzyme-induced saccharification reaction [9] .

Manganese peroxidase is a well-known lignin oxidizing and depolymerizing peroxidase, produced by soil fungi and basidiomycetes, causing white rot [10] , which oxidizes the lignocellulose biomass containing manganese (II) ions (Mn 2+ ) into reactive Mn 3+ . This Mn 3+ is stabilized by chelators such as oxalic acid, which is found in fungus and bacteria. Stabilized Mn 3+ attacks phenolic lignin to disintegrate lignocellulosic biomass spontaneously [11] . Other enzymes of the lignin-degrading family are versatile peroxidase and lignin peroxidase, which possess the unique characteristics of lignin degradation without the use of redox mediators [12] . As reported by established research, the hemicellulose is degraded by cellulases such as acidocaldarius cellulase, Melanocarpus albomyces endoglucanase, and standard endoglucanase cel12A. These cellulases catalyze the disintegration of hemicellulose by hydrolyzing the β-1,4-glycosidic bonds [13] .

Considering these factors, the present study aims to analyze the activity of selective industrial enzymes in biosynthesizing edible glucose and xylose to mitigate the global food and energy demand in the near future. At the same time, it is predicted that the catabolic properties of enzymes can enable human beings to use lignocellulose biomass as a novel alternative source of edible sugars as well as biofuel, depending on their molecular docking and dynamic simulation profiles.

Based on a comprehensive literature review, two substrates, lignin (PubChem CID-73555271) and cellulose microcrystal (PubChem CID-14055602), were selected in the current research. PubChem (https://pubchem. ncbi.nlm.nih.gov/), an National Center for Biotechnology Information authorized open chemistry repository, has been utilized to screen substrates based on in-depth analysis of their physicochemical properties for interacting with enzymes.

On the other hand, lignin peroxidase Protein Data Bank (PDB ID-1LLP), manganese peroxidase (PDB ID-1YYD), and versatile peroxidase (PDB ID-2BOQ) were considered as the main lignin-degrading enzymes, whereas endoglucanase cel12A (PDB ID-1H0B), Acidocaldarius cellulase (PDB ID-3GZK), and Melanocarpus albomyces endoglucanase (PDB ID-1OA9) were selected for cellulose gel (microcrystal) bio-degradation analysis. The set of enzymes was prepared using PDB (http://www.rcsb. org/pdb/).

The active sites of the selective enzymes were identified to ensure supramolecular docking with their corresponding substrates (lignin and cellulose gel) rather than conventional docking [14] . To ensure the best-fit regions of the enzymes to be docked with the substrates, core algorithm (COACH-D) (https://yanglab.nankai.edu.cn/COACH-D/) was used, following the suggested algorithm [15] , so that the binding energy (kcal/mol), residue numbers, and the binding pose can be determined simultaneously.

The 3D structures of the selected substrates (lignin and cellulose gel) were procured in structure data file format from PubChem [16] . Energy minimization was done by reducing the accumulative charge on substrates to zero based on the Gasteiger and Marsili method (1978) and by using the University of California San Francisco Chimera (UCSF) Chimera Software Package (version-1.14) (https:// www.cgl.ucsf.edu/chimera/) [17] . Following optimization, the substrates of interest were transformed into a "mol2 file" for further assessment and molecular docking. All the optimized structures of the substrates have been shown in ( Fig. 1A and 1B) .

The enzymes were optimized sequentially, considering the specific chain selection of the enzymes, the removal of peripheral metal ions, co-factors, water molecules, heteroatoms, additional chains, and non-bonded substrate subunits that were interconnecting with the amino acid residues were made. Besides, the addition of hydrogen atoms and the minimization of energy were conducted due to point-specific docking. The optimization process was also conducted with the UCSF Chimera Software (version-1.14) (https://www.cgl.ucsf.edu/chimera/) [18] . All optimized enzymes were conserved (.pdb format) for the next steps [19] . The lignin and cellulose gel hydrolytic enzymes are presented in Figure 1 , as sequenced-1LLP (C), 1YYD (D), 2B0Q (E), 1H0B (F), 3GZK (G), and 1OA9 (H).

Depending on the positions of the point-specific residues of the enzymes, point-specific super docking was conducted between the optimized substrates and their corresponding enzymes to examine the probable enzyme-substrate binding interactions using PyRx 0.8 docking software [20] . Here, the selected enzyme-substrate complex was changed into the "pdbqt format". Following the supramolecular docking operation, the binding affinity of the substrates was calculated, and root mean square deviation (RMSD) was saved in the "CSV file" for each of the enzyme-substrate complexes.

The visualization software DS Visualizer (https://accelrys-discovery-studio-visualizer.software.informer.com/3.0/) was used for the initial visualization of enzyme-substrate interactions. Afterward, PyMOL version 2.4.1 (https://pymol. org/) was used to analyze and visualize the interactions developed. The enzyme-substrate complexes were saved as "PDB files" from PyMOL for further analysis. Finally, the saved PDB files were analyzed in Ligplot+ (version 2.2) (https:// www.ebi.ac.uk/thornton-srv/software/LigPlus/) for secondary and quantitative networking of interactions (based on the Java Runtime Environment interface) in aspects of marking the potential hydrogen bonding interactions and scopes for non-covalent bond formation [21] .

First, the MDS of the enzymes were run on CABS-flex 2.0 (http://biocomp.chem.uw.edu.pl/CABSflex2) for 10 ns to observe its natural changes in structural orientation and the interactions with adjacent water molecules and ions [22] . Afterward, MDS for individual enzyme-substrate complexes were run on the Substrate and Receptor Molecular Dynamics tools (http://chemyang.ccnu.edu.cn/ccb/server/ LARMD/index.php) up to 3.1 ns for an initial assessment of RMSD, root mean square fluctuation (RMSF), solvent-accessible surface area (SASA), principle component analysis, and Debye-Waller factor for thermostability (B-Factor) [23] .

Finally, the comprehensive MDS of the biological complexes was executed in yet another scientific artificial reality application (YASARA) (version 11.9.18) to confirm the docking pattern and stability of various interactions [24] . The assisted model building with energy refinement-14 force field was used in this system, and the complexes were cleaned at the start of the experiment, and the network of hydrogen bonds was enhanced, and the water molecules were added. Besides, 0.9% NaCl was supplemented at a 310 K temperature, maintaining a pH level of 7.4 [25] . To evaluate the electrostatic-interactions, Particle Mesh Ewald's method was applied [26] . The artificial cell was set at 20Å, which was bigger than all the biological complexes combined. This ensured free movement. The normal simulation time step of 1.25 fs was applied, and each simulation trajectory was conserved after a 100 ps interval. The simulation was run for 50 ns to analyze Cα-RMSD, Å, RMSF, Å, SASA, (Å 2 ), the radius of gyration (Rg, nm), hydrogen bond, and molecular surface area (MolSA, Å 2 ) [19, 21] .

The data generated from the MDS for RMSD, RMSF, SASA, MolSA, and Rg were statistically analyzed using the "R programming language" (version R-4.0.2, for Linux) [27, 28] and GraphPad Prism (version 8.0.1) [29] [30] [31] . The enzyme-substrate complexes were refined and graphically visualized using the aforementioned software packages.

In the enzymes' active site selection process, a diversified number of amino acid residues were found in the target points (ranged between 9 and 27), depending on the enzyme types. Besides, distinguishable binding energy (kcal/mol) levels were predicted, considering the number of residues on the binding domains of the enzymes ( Table 1) .

The lignin binding energy of the optimized enzymes 1LLP, 1YYD, and 2BOQ were −7.0; −6.5, and −8.0 (kcal/mol), respectively ( Table 2) . "1LLP-lignin" complex ( Fig. 2A) (Table 2) . Similarly, "1YYD-lignin" complex ( Fig. 2B) is stabilized 190, 191, 193, 194, 196, 197, 198, 199, 200 Table 3 ). The MDS of 50 ns run time revealed that the RMSD values of the "Cα atom stability" of all complexes ranged between 0.8 and 1.25 Å. The "2BOQ-lignin" complex showed a better fluctuation rate between 30 and 50 ns time-period than the resting complexes (Fig. 4A) . The MolSA (Å 2 ) of the enzyme-substrate complexes was scored between 15,202.69 and 17,497.95 Å 2 (Fig. 4C) for the lignin hydrolytic enzymes. Analysis of each enzyme's flexibility resulted in RMSF values of less than 2 Å, indicating the enzymes' more stable nature and a lower number of alpha-helix and beta-sheet entities (Fig. 5A) at the structural interaction points. The Rg (nm) values represent the tightness of the protein structure, where a lower degree of fluctuation stands for the significant uniformity and rigidity of the system. Among the three complexes, the 2BOQ-lignin complex showed the more compact structure following the value of 19.50, while the other complexes showed values as −19.80 (1YYD-lignin) and 19.60 (1LLPlignin), respectively, which are less significant than the 2BOQ-lignin complexes (Fig. 5C) . Besides, SASA confirmed the expansion of the enzymatic volume during the 50ns dynamics simulation process, where the "2BOQ-lignin" complex showed an effective output greater than the other complexes and represented the significant SASA value at 12,600 Å 2 . On the other hand, enzyme complexes such as 1LLP-lignin and 1YYD-lignin demonstrated higher SASA values at 140,002 and 145,002, respectively. Notably, the lowest SASA value of the system displayed more significant results that were found from the 2BOQ-lignin interaction (Fig. 5D) . Among the lignin-degrading enzymes, 1LLP displayed polar and apolar energy per area at 4,592.05 and 9,123.59, respectively, whereas 1YYD produces polar (Table 4 ). Enzymesubstrate structural strength considering the number of intermolecular hydrogen bonds are 238, 240, and 240 for the 2BOQ-lignin, 1YYD-lignin, and 1LLP-lignin sequentially. As the lowest number of hydrogen bonds, which showed maximum stability, 2BOQ was the significant one.

The binding affinity of the complex between "1OA9cellulose" was found to be more significant (−6.7 kcal/mol) than the "3GZK-cellulose" and "1H0B-cellulose" complexes, with energy scales of −5.3 and −6.0 kcal/mol, respectively ( Table 2 ). Tyr163 (3.05), Glu124 (2.71), Trp26 (3.14), and Glu207 (3.19, 3.04, and 2.93) hydrogen bond interactions form complex 1H0B-cellulose (Fig. 3A) . Simultaneously, complex hydrophobic networking is developed by cellulose with the neighboring amino acid residues (Pro137, Gly138, Trp159, Trp161, Trp108, Asp106, Tyr59, Trp128, Met126) of the enzyme (Table 3) . Similarly, the complex between "3GZK-cellulose" (Fig. 3B) is mainly a network of hydrogen bonds with the residues-Arg410 (3.17 and 3.32 Å), Asp398 (3.29 Å), Asp345 (2.84 Å), Trp401 (2.94 Å), Gly342 (2.97 Å), and hydrophobic interactions with Glu340, IIe400, Ala344, Met4005, Trp343, and Val406 (Table 3 ). Asp10 (3.18), Leu130 (2.98), Val128 (3.05), and other residues (Tyr146, Gly127, Gly129, Tyr8, Tyr110, Gly112, Thr6, and Asn178) ( Table 3) .

Representative RMSD values of cellulose gel complexes complexed with 1HOB, 3GZK, and 1OA9 were between 0.6 and 1.20, with the 1OA9-cellulose complex showing a higher fluctuation rate within the 30 to 45 ns time period than the other complexes (Fig. 4B) . The enzyme-substrate complexes' MolSA resulted in values ranging from 9,138.552 to 22,534.41 Å 2 (Fig. 4C) . The RMSF values for the enzymes' flexibility assessment did not exceed 2.4 Å (except 3GZK), referring to a more stable enzyme structure with fewer beta-sheet and alpha-helix regions. The enzyme-3GZK possesses an RMSF value of 3.7 Å that the protein structure is less stable and bears more α-helix and/or β-sheet regions resulting from the dynamic simulation (Fig. 5B) . Among the three complexes, the 1OA9cellulose gel complex showed a more compact and rigid structure following the value of 15.80 collected from the Rg output data. Besides, the other complexes showed values of 23.20 (3GZK-Cellulose gel) and 16.50 (1H0B-Cellulose gel) that are less significant than the former complex (Fig. 5C ). The SASA profiling conformed to the expansion of the "1OA9-cellulose gel" complex surface area during the 50 ns dynamics simulation process as more significant than the other complexes, representing an area size of 84,002, whereas the enzyme complexes' "1H0B-cellulose" and "3GZK-Cellulose" documented comparatively higher SASA values, respectively. It is important to recommend that the lowest SASA value of the system display more significant results with the 1OA9-cellulose gel complex (Fig. 5D) . Amongst all the enzymes, 1OA9 had the lowest polar energy per area (3,340.94). On the contrary, 1H0B had the lowest apolar energy generated per 2 area (3,340.94). 3GZK demonstrated stark contrast, as it had the largest difference between the polar and apolar values, 7,632.39 and 11,981.65, respectively ( Table 4 ). The intermolecular hydrogen bonds are 160, 180, and 400 for the 1OA9-Cellulose, 1H0B-Cellulose, and 3GZK-Cellulose complexes, respectively. The lowest number of hydrogen bonds showed a more stable complex formation of cellulose with 1OA9 than with all other enzymes. Considering all the MDS profiles like RMSD, RMSF, SASA, MolSA, Rg, and hydrogen bonds, the catabolic significance of substrate-degrading enzymes in the current research has resulted in 2BOQ <1LLP <1YYD. Similarly, the catabolic profiles of the cellulose-degrading enzymes resulted in 1OA9 <1H0B < 3GZK. The consequence of cellulose degradation is the production of glucose and xylose, which are simple forms of sugar and are edible for human beings as an energy source.

Food production, preservation, and supply for the ever-increasing global population have become a matter of great challenge, where alternative edible sources of food and energy are desperately needed. The use of non-food lignocellulose biomass can be a good option for producing edible sugars for energy if they are properly manipulated [32, 33] . The catalytic activity of the enzymes can play a significant role in the lignocellulose manipulation procedure, where the lignin-degrading enzymes first degrade the lignin wall to expose cellulose and hemicellulose [34, 35] . Afterward, the cellulose microcrystals are hydrolyzed with the cellulose-degrading enzymes to produce glucose and xylose, which are completely edible for human beings [36] . Considering all those things, in this current research, lignin (PubChem CID-73555271) and cellulose gel (PubChem CID-14055602) were used as the main substrates on which the catabolic activity of enzymes was studied depending on their enzyme-substrate complexing characteristics, inter-molecular interactions, and MDS properties. It's very pivotal to investigate the active sites of enzymes to conduct point specific super docking instead of conventional docking, which was undertaken by this researcher using the high-throughput tool COACH-D and its algorithm [14, 15] . This resulted in a sufficient number of amino acid residues with enormous binding energies (kcal/ mol) ( Table 1) . To study lignin hydrolysis activity-manganese peroxidase (PDB ID-1YYD), versatile peroxidase (PDB ID-2BOQ), and lignin peroxidase (PDB ID-1LLP) enzymes were used [37, 38] , whereas Acidocaldarius cellulase (PDB ID-3GZK), Melanocarpus albomyces endoglucanase (PDB ID-1OA9) [39, 40] , and standard endoglucanase cel12A (PDB ID-1H0B) [41] were used for cellulose degradation.

To understand the catabolic activity of the enzymes toward lignin and cellulose substrates, depending on their binding affinity and RMSD values, molecular docking was conducted individually [42] following the optimization of all the substrates and enzymes. The "2BOQ-lignin" and "1OA9-cellulose gel" resulted in the best fitting scores of −8.00 and −6.7 (kcal/mol), respectively ( Table 2 ). The docking score detects comparatively selective, potent, and efficient candidates for surface modeling [43] . Moreover, for the analysis of the enzyme-substrate interactions, to understand the hydrophobic and hydrogen bonding interactions ( Table 3 ) between each of the enzyme-substrate complexes, their quantitative networking model was developed using "Java interface" of computer programming ( Figs. 2 and 3 ). MDS has proven to be the best authentic method for investigating biomolecular interactions, evaluating the atomic level's stability and the resulting output data of the dynamic trajectories, which paves the way for the quantitative relationship between enzyme and substrate structure and functions [44] . In this current study, the YASARA (version 11.9.18) dynamics simulator was conducted, maintaining all the physio-chemical and physiological parameters (Temperature-310K, pH-7.4; 0.9% salt concentration; addition of ions like Na+ and Cl−, etc.) during the 50 ns simulation run time [45] , and the output data resulted from the dynamics state of the enzyme's capability to catabolize the lead substrates (Lignin, Cellulose gel The RMSD values of the lignin-degrading 1LLP, 1YYD, and 2BOQ enzymes were 1.25, 1.30, and 1.10 Å sequentially (Fig. 4A) . On the other hand, the RMSD values were 1.20, 0.80, and 0.75 Å for 1H0B, 3GZK, and 1OA9, respectively (Fig. 4B) . The finding relates the catabolic properties of the enzymes with the other established findings [46, 47] .

To analyze the RMSD values (for frame "x"), the following algorithm was used-

Here, "N" is the number of atoms in the atom selection; "t ref " is the reference time (typically the first frame is used as the reference, and it is regarded as time, "t = 0"); and "r" is the position of the selected atoms in frame "x" after superimposing on the reference frame, where frame x is recorded at the time : "t x ." The procedure is repeated for every frame in the simulation trajectory.

There are two types of available MolSA calculations, namely-vdW and solvent accessible [48] . More commonly, the enzyme-substrate complex stability should be increased with the decrease of the surface area (Å 2 ) [49] . Considering the criteria, 2BOQ-lignin and 1OA9-cellulose gel complexes were the more stable ones, with 15,202.688 and 9,138.552 Å 2 , respectively (Fig. 4C ). In the same way, the fluctuation of all the complexes resulted below 2 Å (Fig. 5A and B) except the "3GZK-cellulose gel" complex ( Fig. 5B) , from the RMSF values of the 50 ns MDS. RMSF data represent the flexibility and strength of enzyme-substrate complexes (Fig. 5) . The algorithm used for the characterization of local changes along the enzyme chain-

Here, "T" is the trajectory time over which the RMSF is calculated, "t ref " is the reference time (usually for the first frame, and is regarded as the zero of time); "r" is the position of atom "i" in the reference at the time "t ref ", and "r" is the position of the atom "i" at the time "t" after superposition on the reference frame.

The Rg represents the tightness of the protein structure, where a decreased fluctuation level indicates the increased firmness and rigidity of the system [21] . Among the six complexes, "2BOQ-lignin" and "1OA9-cellulose gel" complexes showed the more compact values means 19.50 and 15.80 nm, respectively (Fig. 5C) . The SASA uses parameters to check the changes in the arrangement of the enzyme-substrate complexes retrieved from the molecular docking and MDS to figure out the efficacy of the interaction between the enzyme and substrate. It detects the interactions occurring between the complex surface and the water molecules based on the total energy per area of the substrate and the compound [50] . Higher SASA values indicate unstable structures due to increased contact between the hydrophobic amino acids with the water [51] . Contact between the hydrophobic part and the water is highly unfavorable, destabilizing the interaction between the enzyme and substrate. The interaction between lignin and 1YYD was highly unfavorable due to their higher total energy per area. At the same time, 2BOQ-lignin scored the lowest total energy per area (Table 4) , causing them to have the strongest interaction [52] . Similarly, the "1OA9-cellulose gel" complex possesses the lowest total energy per area, ensuring their highest interaction profile than the others (Fig. 5D ). It is important to mention that enzymes in both the lignin as well as Cellulose Gel had zero value under unknown. This is because the dynamic simulation conducted was accurate, as no unknown residue was found [53] . The enzyme-substrate complexes' conformational stability was determined via the number of hydrogen bonds retrieved from the 2BOQlignin (238) and 1OA9-Cellulose gel (160) complexes.

This in silico study focuses on the catabolic activity of numerous industrial enzymes in the degradation of nonfood lignocellulose biomass into simple sugars (glucose, xylose) based on enzyme-substrate interactions and structure-based molecular networking. This cutting-edge method for biosynthesizing food energy from completely non-food sources by modulating enzyme bioconversion activity could be a revolutionary option for securing a hunger-free world for future generations, which is the primary goal of this research.

The catabolic activity of six enzymes was examined in silico to configure their efficacy in degrading lignocellulose biomass into functional sugars as food energy, as a part of a study of alternate food supplies for the future to alleviate hunger. The designed enzymes were optimized using energy minimization and applied to molecular docking on lignin and cellulosic components, with a 50 ns runtime of the MDS, resulting in RMSD, MolSA, RMSF, SASA, Rg, and Hydrogen bond amount values in the favorable ranges, indicating the catabolic strength of the enzymes of interest. Using all parameters of docking and dynamics simulations, the researchers discovered that Versatile Peroxidase and Melanocarpus albomyces endoglucanase enzymes have the most efficient features for bioconverting lignin and cellulose microcrystals, respectively. Researchers have had a lot of trouble figuring out which biocatalysts are appropriate for the saccharification process before fermentation. This research gives an authentic means of monitoring and predicting prospective enzyme activity for a successful hydrolytic biological conversion process before doing so directly in the fermenter to produce edible sugars at that time. The destroyed lignin structure exposes the cellulosic contents, which are eventually converted into glucose and xylose as edible sugars for humans, who are in desperate need of ways to alleviate the global food crisis, produce biofuels, and manage agricultural waste. This study will aid the next generation of scientists, particularly bioinformatics and bioengineering experts, in pre-examining the catalytic potentialities of their enzymes of interest for a successful biological reaction before conducting it inside any bioreactor system. Agricultural waste materials can be anticipated to produce a variety of value-added goods, ranging from basic food components to carbon-neutral biofuel generation, using this type of in silico technique. 

World population growth. Our World in Data

A national profile of family and unpaid caregivers who assist older adults with health care activities

Youth food insecurity in Ghana and South Africa: prevalence, socioeconomic correlates, and moderation effect of gender

Lignocellulosic biomass: understanding recalcitrance and predicting hydrolysis

Uptake and trophic transfer of nitrogen and carbon in a temperate forested headwater stream

Consolidated briefing of biochemical ethanol production from lignocellulosic biomass

β-Glucosidase immobilisation on synthetic superparamagnetic magnetite nanoparticles and their application in saccharification of wheat straw and Eucalyptus globulus pulps

Pretreatment of lignocellulosic materials as substrates for fermentation processes

Effective one-step saccharification of lignocellulosic biomass using magnetite-biocatalysts containing saccharifying enzymes

Induction, purification and characterization of a novel manganese peroxidase from Irpex lacteus CD2 and its application in the decolorization of different types of dye

Mechanisms of manganese (ii) oxidation by filamentous ascomycete fungi vary with species and time as a function of secretome composition

Membrane-bound class III peroxidases: unexpected enzymes with exciting functions

Cellulase-lignin interactions-the role of carbohydrate-binding module and pH in non-productive binding

COACH-D: improved protein-ligand binding sites prediction with refined ligand-binding poses through molecular docking

Protein-ligand binding site recognition using complementary binding-specific substructure comparison and sequence profile alignment

PubChem in 2021: new data content and improved web interfaces

Molecular docking using chimera and autodock vina software for nonbioinformaticians

UCSF ChimeraX: meeting modern challenges in visualization and analysis

Molecular optimization, docking, and dynamic simulation profiling of selective aromatic phytochemical ligands in blocking the SARS-CoV-2 S protein attachment to ACE2 receptor: an in-silico approach of targeted drug designing

Small-molecule library screening by docking with PyRx

Molecular optimization, docking and dynamic simulation study of selective natural aromatic components to block E2-CD81 complex formation in predating protease inhibitor-resistant HCV influx

CABS-flex 2.0: a web server for fast simulations of flexibility of protein structures

LARMD: integration of bioinformatic resources to profile ligand-driven protein dynamics with a case on the activation of estrogen receptor

New ways to boost molecular dynamics simulations

Identification and in silico molecular modelling study of newly isolated Bacillus subtilis SI-18 strain against S9 protein of Rhizoctonia solani

Screening of phytochemicals as potent inhibitor of 3-chymotrypsin and papain-like proteases of SARS-CoV2: an in-silico approach to combat COVID-19

RT-PCR mediated identification of SARS-CoV-2 patients from particular regions of Bangladesh and the multi-factorial analysis considering their pre and post infection health conditions

The serological basis of the correlation between iron deficiency anemia and thyroid disorders in women: A community based study

Molecular identification of a Bacillus cereus strain from Murrah buffalo milk showed in vitro bioremediation properties on selective heavy metals

In ovo inoculation of duck embryos with different strains of Bacillus cereus to analyse their synergistic post-hatch anti-allergic potentialities

An experimental analysis of different point specific musculoskeletal pain among selected adolescent-club cricketers in Dhaka city

Utilization of lignocellulosic waste by the edible mushroom, Pleurotus

Characterization of bionanocomposite films prepared with agar and paper-mulberry pulp nanocellulose

Bioavailability of carbohydrate content in natural and transgenic switchgrasses for the extreme thermophile Caldicellulosiruptor bescii

Lignin-degrading enzymes

Enzymatic degradation of (ligno) cellulose

Lignin degradation: a microbial approach

Redesign of a new manganese peroxidase highly expressed in Pichia pastoris towards a lignin-degrading versatile peroxidase

Evaluation of temperature, pH and nutrient conditions in bacterial growth and extracellular hydrolytic activities of two Alicyclobacillus spp. strains

Cellulolytic potential of thermophilic species from four fungal orders

Effect of endoglucanases from different glycoside hydrolase families on enzymatic preparation of cellulose nanocrystal

Exploration and evaluation of bioactive phytocompounds against BRCA proteins by in silico approach

Structure-based docking studies towards exploring potential anti-androgen activity of selected phytochemicals against

Molecular docking and dynamic simulation of AZD3293 and solanezumab effects against BACE1 to treat Alzheimer's disease

Putative inhibitors of SARS-CoV-2 main protease from a library of marine natural products: A virtual screening and molecular modeling study

Lignin peroxidase functionalities and prospective applications

Heterologous expression of the thermostable cellulase CEL45A in Trichoderma reesei (Doctoral dissertation

Improving the economy of lignocellulose-based biorefineries with organosolv pretreatment

Natural compounds as inhibitors of SARS-COV-2 main Protease (3CLpro): a molecular docking and simulation approach to combat COVID-19

Molecular dynamics simulation revealed binding of nucleotide inhibitors to ZIKV polymerase over 444 nanoseconds

Computational discovery of plant-based inhibitors against human carbonic anhydrase IX and molecular dynamics simulation

Identification of conformational B-cell epitopes in diphtheria toxin at varying temperatures using molecular dynamics simulations

Improving prediction of protein secondary structure, backbone angles, solvent accessibility and contact numbers by using predicted contact maps and an ensemble of recurrent and residual convolutional neural networks

The authors are grateful to the RPG (Govt. License ID: 05-060-06021) for providing all types of unconditional supports in the technical issues (under the Project of Category: D1; ID. #04-2020/21). The authors are also grateful to Dr. Sharmin Ahmed (Chulalongkorn University, Thailand)

The authors have no competing interest at all with the others.