Defining the Frontiers of Synergism between Cellulolytic Enzymes for Improved Hydrolysis of Lignocellulosic Feedstocks

: Lignocellulose has economic potential as a bio-resource for the production of value-added products (VAPs) and biofuels. The commercialization of biofuels and VAPs requires efﬁcient enzyme cocktail activities that can lower their costs. However, the basis of the synergism between enzymes that compose cellulolytic enzyme cocktails for depolymerizing lignocellulose is not understood. This review aims to address the degree of synergism (DS) thresholds between the cellulolytic enzymes and how this can be used in the formulation of effective cellulolytic enzyme cocktails. DS is a powerful tool that distinguishes between enzymes’ synergism and anti-synergism during the hydrolysis of biomass. It has been established that cellulases, or cellulases and lytic polysaccharide monooxygenases (LPMOs), always synergize during cellulose hydrolysis. However, recent evidence suggests that this is not always the case, as synergism depends on the speciﬁc mechanism of action of each enzyme in the combination. Additionally, expansins, nonenzymatic proteins responsible for loosening cell wall ﬁbers, seem to also synergize with cellulases during biomass depolymerization. This review highlighted the following four key factors linked to DS: (1) a DS threshold at which the enzymes synergize and produce a higher product yield than their theoretical sum, (2) a DS threshold at which the enzymes display synergism, but not a higher product yield, (3) a DS threshold at which enzymes do not synergize, and (4) a DS threshold that displays anti-synergy. This review deconvolutes the DS concept for cellulolytic enzymes, to postulate an experimental design approach for achieving higher synergism and cellulose conversion yields. to prevent product inhibition. c The enzymes were immobilized on synthetized Fe 3 O 4 nanoparticles.


Introduction
Lignocellulosic feedstocks have huge economical potential as a source of value-added products (VAPs) and biofuel [1][2][3]. Cellulose is the major component of lignocellulose, and it can be used as a source of fermentable sugars that can serve as precursors for VAPs synthesis. Cellulose is sourced from plant material (wood pulp, cotton, or cereals, such as wheat, sugarcane, and rice bagasse), and is a linear polymer that consists of β-Dglucose molecules linked by glycosidic bonds [1,[4][5][6]. The linear structure of a cellulose chain is directional, as it consists of a reducing-end glucose that contains an anomeric carbon (C1) and a non-reducing-end glucose, consisting of hydrogen and a hydroxyl group on the C4 carbon ( Figure 1). Many cellulose chains bundle via hydrogen bonding to constitute cellulose microfibrils, which consist of crystalline regions intersected by amorphous regions. The β-D-glucose residues in the crystalline regions of the cellulose microfibril are anhydrous; this is because the orientation of the three, free hydroxyl groups of the sugar create strong molecular bonds with glucose molecules in adjacent chains, which make cellulose insoluble in water [2]. Thus, crystalline cellulose is recalcitrant to Figure 1. A representation of a cellulose chain, with reducing and non-reducing ends indicated. The blue brackets with subscript "n" indicate that the degree of polymerization (DP) of the cellulose chain is longer than a DP of 10. It is well understood that several glycoside hydrolases (GHs) from different families act synergistically to achieve complete hydrolysis of the cellulosic fraction of the biomass feedstock into fermentable sugars [3,7,[10][11][12]. There are the following two types of mechanisms of synergy that have been documented to date: (1) simultaneous synergism, whereby all enzymes are added into the reaction at the same time and (2) sequential synergism, whereby enzymes are added into a reaction sequentially, based on their substrate cleavage preferences [11,13]. The simultaneous application of synergistic enzymes is the most preferred technique for degrading biomass feedstocks; for instance, the cellulolytic enzyme homeo-synergism is based on the simultaneous use of the following cellobiohydrolases: CBHI (EC 3.2.1.176) and CBHII (EC 3.2.1.91), endoglucanase (EG, EC 3.2.1.4) and β-glucosidase (β-gl, EC 3.2. 1.21). Several research groups have demonstrated synergism between CBHI, CBHII, EG, and β-gl sourced from fungi or bacteria, e.g., Trichoderma reesei and Clostridium thermocellum, during the degradation of cellulose [5,[12][13][14][15]. Figure 2 shows the well-established mode of action of cellulases during the synergistic hydrolysis of cellulose. Further analysis of the literature reveals that cellulolytic enzyme synergism can be grouped into the following two classes: (1) synergism between non-catalytic and catalytic cellulolytic enzymes (e.g., expansins and carbohydrate-binding modules (CBMs)), and (2) synergism between catalytic cellulolytic enzymes, e.g., exo-exo, exoendo, exo-endo-β-gl, or cellulase and lytic polysaccharide monooxygenase (LPMO) synergy [16][17][18][19][20][21]. There are only a few studies that have documented endo-endo synergy, even though EGs constitute the most diverse cellulase group, classified under GH families 5,6,7,8,9,10,12,44,45,48,51,74, and not classified sequences (NC) (http://www.cazy.org/Glycoside-Hydrolases.html; accessed on 30 June 2021). It is also well documented that CBHs from different GH families that hydrolyze cellulose on the same cleavage sites do not synergize together, but compete for occupancy on the same binding sites [15,22].
Two other enzyme systems play a significant role in cellulose saccharification, these being LPMOs and non-hydrolytic enzymes (expansins and swollenins). Forsberg et al. [23] demonstrated that the LPMOs sourced from bacteria and fungi were originally classified into CBM family 33 and GH61, respectively. According to the CAZy database, the enzymes classified under CBM33 were later assigned to auxiliary activities family 10 (AA10), while those in GH61 were assigned to AA9 (http://www.cazy.org/Auxiliary-Activities.html; accessed on 30 June 2021). The LPMOs in family AA9 and AA10 are copper-dependent enzymes that cleave the glycosidic bonds of crystalline cellulose regions by oxidizing the C1 (EC 1.14.99.54) or C4 (EC 1.14.99.56) and C1/C4 of the glucose residues [23,24] (see Figure 3). The LPMOs in family AA16 are newly classified cellulose-specific enzymes with C1 regiospecificity [25]. Kadowaki et al. [24] demonstrated that some LPMOs sourced from Myceliophthora thermophila, belonging to AA9D, have regiospecificity for C1, and those sourced from Aspergillus nidulans, belonging to AA9 (also known as 22), and from M. thermophila, belonging to AA9J, have regiospecificity for C4. Additionally, the LPMOs sourced from M. thermophila (belonging to AA9I), Thielavia terrestris (belonging to AA9E), and Thermoascus aurantiacus (belonging to AA9A) displayed regiospecificity for both C1 and C4 (referred to as C1/C4). Both bacterial and fungal LPMOs have similar structures and mechanisms of action [23]. The catalytic domain of LPMOs can also be attached to CBMs by a flexible linker [24,26]. Furthermore, Isaksen et al. [27] demonstrated that the fungal LPMO-AA9C, sourced from Neurospora crassa, was capable of cleaving polymeric cellulose and cello-oligosaccharides with a degree of polymerization of about four to six. The LPMOs cleaved the C1 of cellulose or cello-oligosaccharides to produce aldouronic acids; however, when they acted on the C4, they produced gem-diol [26,27]. The LPMOs from families AA9, AA10, and AA16 have been reported to synergize with cellulases to enhance the saccharification of cellulose [25,28,29]. LPMOs sourced from Myceliophthora thermophila, belonging to AA9D, have regiospecificity for C1, and those sourced from Aspergillus nidulans, belonging to AA9 (also known as 22), and from M. thermophila, belonging to AA9J, have regiospecificity for C4. Additionally, the LPMOs sourced from M. thermophila (belonging to AA9I), Thielavia terrestris (belonging to AA9E), and Thermoascus aurantiacus (belonging to AA9A) displayed regiospecificity for both C1 and C4 (referred to as C1/C4). Both bacterial and fungal LPMOs have similar structures and mechanisms of action [23]. The catalytic domain of LPMOs can also be attached to CBMs by a flexible linker [24,26]. Furthermore, Isaksen et al. [27] demonstrated that the fungal LPMO-AA9C, sourced from Neurospora crassa, was capable of cleaving polymeric cellulose and cello-oligosaccharides with a degree of polymerization of about four to six. The LPMOs cleaved the C1 of cellulose or cello-oligosaccharides to produce aldouronic acids; however, when they acted on the C4, they produced gem-diol [26,27]. The LPMOs from families AA9, AA10, and AA16 have been reported to synergize with cellulases to enhance the saccharification of cellulose [25,28,29]. Expansins and expansin-like proteins are non-hydrolytic or non-oxidative enzymes that can bind to the crystalline cellulose and weaken the degree of crystallinity of the polymer by breaking the inter-or intra-molecular hydrogen bonding network [30,31]. Kim et al. [32] demonstrated that Bacillus subtilis-derived expansin sequences were homologous to maize plant expansins, and their 3D structures were also similar. The role of expansins in plants is well established, and plant-derived expansins are classified into two groups, Comment shown co Figure 2. A schematic representation of cellulose and corresponding cellulolytic enzymes that degrade it. The enzyme cleavage sites on the cellulose are represented by red arrows. Microcrystalline and amorphous (green color) regions of cellulose consist of β-D-glucose linked by β-1,4-glycosidic bonds. Where β-gl, CBHI, CBHII and EG represent β-glucosidase, cellobiohydrolase I, cellobiohydrolase II and endo-glucanase.

Comment
Expansins and expansin-like proteins are non-hydrolytic or non-oxidative enzymes that can bind to the crystalline cellulose and weaken the degree of crystallinity of the polymer by breaking the inter-or intra-molecular hydrogen bonding network [30,31]. Kim et al. [32] demonstrated that Bacillus subtilis-derived expansin sequences were homologous to maize plant expansins, and their 3D structures were also similar. The role of expansins in plants is well established, and plant-derived expansins are classified into two groups, αand β-expansins [4]. The α-and β-expansins only share 20% amino acid sequence similarity. Interestingly, expansins generally have about 30% sequence similarity to endoglucanases classified under GH45 [4]. Several studies have shown that the expansins synergize with cellulases that hydrolyze the crystalline microfibrils of cellulose [4,30,31,33]. Fungi also possess another non-hydrolytic protein called swollenin, which is similar to the expansins [33]. Swollenins are reported to modify the chemistry and structure of microcrystalline cellulose by reducing its degree of crystallinity, allowing the cellulases to hydrolyze cellulose effectively. Additionally, studies have demonstrated that cellulases and swollenins synergize to improve the saccharification of cellulose [17,33]. In this literature review, we used the DS, a powerful parameter that describes molecular interactions between the enzymes during synergy, to explain hydrolytic cellulolytic enzyme to non-hydrolytic enzyme synergism thresholds, to develop a model to assist in the improvement of cellulose hydrolysis by cocktails. gize with cellulases that hydrolyze the crystalline microfibrils of cellulose [4,30,31,33]. Fungi also possess another non-hydrolytic protein called swollenin, which is similar to the expansins [33]. Swollenins are reported to modify the chemistry and structure of microcrystalline cellulose by reducing its degree of crystallinity, allowing the cellulases to hydrolyze cellulose effectively. Additionally, studies have demonstrated that cellulases and swollenins synergize to improve the saccharification of cellulose [17,33]. In this literature review, we used the DS, a powerful parameter that describes molecular interactions between the enzymes during synergy, to explain hydrolytic cellulolytic enzyme to nonhydrolytic enzyme synergism thresholds, to develop a model to assist in the improvement of cellulose hydrolysis by cocktails. Figure 3. A schematic representation of cellulose chain and corresponding lytic polysaccharide monooxygenases (LPMOs) that degrade it. LPMOs can cleave glycosidic bonds joining the glucose units of cellulose at C1 (green block), and at either C1 or C4 (grey triangles). The C1 regiospecific LPMOs produce glucono-δ-lactone and gluconic acid (brown arrows), which is stable in solution. In contrast, LPMOs with C1/C4 regiospecificity can produce a mixture of gluconic acid (green and brown arrow) and 4-gem-diol-glucose (blue arrows).
The DS between enzymes is defined as the observed activity of two or more enzyme combinations divided by the theoretical sum of the activities of the individual enzymes acting on the same substrate. The enzyme combinations are formulated using three protein loading systems, which are protein mass ratio, molar ratio, and activity ratio [11,12,15,20,34]. According to van Dyk and Pletschke [11], the DS can provide important information about the chemical structure of the substrate and enzyme interactions during the hydrolysis of the substrate. Two or more enzymes cooperating during the hydrolysis of the substrate lead to higher soluble sugar production ( Figure 2). Additionally, the mass enzyme loading in the reaction may be reduced due to synergism between enzymes, and this could lead to lowered enzyme costs during biomass valorization [5,32,35]. Van Dyk and Pletschke [11] further demonstrated that the DS is a particularly important component for synergy studies, yet not all experiments performed to investigate the GHs or GHto-LPMO synergism factor in the quantification of the DS.

Commente
shown com Figure 3. A schematic representation of cellulose chain and corresponding lytic polysaccharide monooxygenases (LPMOs) that degrade it. LPMOs can cleave glycosidic bonds joining the glucose units of cellulose at C1 (green block), and at either C1 or C4 (grey triangles). The C1 regiospecific LPMOs produce glucono-δ-lactone and gluconic acid (brown arrows), which is stable in solution.
The DS between enzymes is defined as the observed activity of two or more enzyme combinations divided by the theoretical sum of the activities of the individual enzymes acting on the same substrate. The enzyme combinations are formulated using three protein loading systems, which are protein mass ratio, molar ratio, and activity ratio [11,12,15,20,34]. According to van Dyk and Pletschke [11], the DS can provide important information about the chemical structure of the substrate and enzyme interactions during the hydrolysis of the substrate. Two or more enzymes cooperating during the hydrolysis of the substrate lead to higher soluble sugar production ( Figure 2). Additionally, the mass enzyme loading in the reaction may be reduced due to synergism between enzymes, and this could lead to lowered enzyme costs during biomass valorization [5,32,35]. Van Dyk and Pletschke [11] further demonstrated that the DS is a particularly important component for synergy studies, yet not all experiments performed to investigate the GHs or GH-to-LPMO synergism factor in the quantification of the DS.
The lack of investigation of the DS during enzymatic cellulose saccharification is attributed to a focus on hydrolysis yields by researchers. It appears that some enzyme combinations do not display a positive DS, yet these combinations demonstrate higher yields, in terms of soluble sugar production, compared to when the individual enzymes are used to hydrolyze the substrate. However, we propose that investigating the synergism between enzymes should not only be limited to yield improvement, but should also include the DS. This would enable researchers and companies to design better enzyme cocktails for specific feedstocks. This short review will address how we can achieve higher yields and DS by re-evaluating the established synergisms and setting the DS threshold benchmarks, which can lead to higher yields.

Defining DS Threshold and Its Impact on Substrate Hydrolysis
It is well accepted that a DS of 1 arbitrary unit (AU) demonstrates no synergism between the enzymes used to formulate a said enzyme cocktail [11,12,15,20,36,37]. These studies suggest that the enzymes do not cooperate during substrate hydrolysis, leading to no improvement in the rate of substrate deconstruction compared to when the monocomponents act individually; for example, endoglucanase and exoglucanase enzymes do not synergize, and generally result in a DS of 1 AU, during the hydrolysis of the model amorphous cellulose carboxyl methylcellulose (CMC), because of their uncomplimentary modes of action on the substrate. A value below 1 AU shows that there is no synergism between the enzymes during the hydrolysis of the substrate. This observation suggests that one or more of the enzymes cannot hydrolyze the substrate efficiently, due to competitive binding of two or more enzymes on the same substrate cleavage site (same regiospecific enzymes). Instances that resulted in a DS below 1 AU have been documented between CBHII and EG, processive EGs from GH9 and CBHI or CBHI (from GH48), and CBHI (from GH7) or CBHI (GH7), and an LPMO (AA9E). In all these cases, substrate hydrolysis was impeded, and product formation significantly reduced. Generally, DS values above one are achieved when the enzymes in the cocktails hydrolyze different cleavage sites of the substrate, e.g., CBHI, CBHII, EG, and β-gl degrade cellulose at different cleaving sites, resulting in significantly higher product formation ( Figure 2).
It is reported that enzyme synergism and its DS are influenced not only by the substrate, but by other factors covered in many other studies [5,11,37]. The quantification of DS levels is generally based on the enzyme activity or specific activity, the conversion rate and/or the yield of product in relation to substrate quantity [11,12,15,38]. There are no thresholds that distinguish between the patterns of enzymes competing for the same cleavage sites on the substrate and those that do not cooperate at all during substrate hydrolysis. Studies have only focused on reporting that a DS value < 1 results in no synergism. We, therefore, hypothesize that 1 ≥ DS values > 0.5 demonstrates that the enzymes are not cooperating on the substrate, while 0.5 ≥ DS values > 0 could signify competitive binding to the same substrate's cleavage site (same regiospecificity), and DS values > 1 could signify that there is synergism between non-hydrolytic and hydrolytic cellulolytic enzymes. The following sections of this review attempt to explain and formulate the DS threshold, whereby cellulolytic enzyme synergism can result in a higher yield (DS > 1), or inhibited yield due to non-cooperating (1 ≥ DS value > 0.5) and competing enzymes (0.5 ≥ DS value > 0).

Unravelling Enzyme Synergism
We hypothesize that, to formulate effective enzyme cocktails, researchers should not only focus on the yield of soluble sugars, as conventional synergy studies have done for the past several decades, but should focus on both the yield of soluble sugars and the DS between enzymes. This approach will generate holistic information about the reaction, with regards to the interactions between the enzymes in the cocktails, as well as the products formed. Hence, we attempted to establish the DS threshold that results in high yields of soluble sugars. A detailed discussion of various cellulolytic enzyme synergy studies that were conducted with different enzyme combinations are presented below and summarized in Table 1.

Non-Hydrolytic and Catalytic Cellulolytic Enzyme Synergy
The crystalline nature of cellulose requires a concerted effort of non-hydrolytic enzymes that weaken the hydrogen bonds between the cellulose fibers, or assist with threading and transferring the cellulose fibers to the catalytic cleft of hydrolytic enzymes [4,5]. These non-hydrolytic cellulolytic enzymes are grouped under CBMs, expansins, expansinlike proteins, and swollenins [4,31,[39][40][41]. Lee et al. [30] demonstrated that expansins reduce the crystallinity of the cellulose, enhancing the activity of the commercial cellulase preparation Celluclast ® 1.5 L 4.8-fold compared to when the cellulase preparation was used alone. A DS value of 4.6 AU between Celluclast ® 1.5 L and an expansin was predicted by Lee et al. [30]. Pech-Cervantes et al. [41] demonstrated that the synergy between an expansin-like protein and cellulase did not have a significant impact on the hydrolysis of amorphous cellulose substrate (CMC). However, the expansin-like protein modified the structure of the crystalline cellulose substrate, filter paper, leading to synergy between it and the cellulase, whereby a four-fold improvement in sugar release was recorded. Again, the calculated DS value, using the published data, computed a DS value of 4.2 AU. A similar impact of synergism was reported by Kim et al. [42], where expansin-to-cellulase synergism increased the release of reducing sugars by 240% (corresponding to a calculated DS value of 2.4 AU), compared to the use of cellulase alone. These studies show that there can only be the following two thresholds with regards to expansin-to-cellulase synergism: (1) DS greater than or equal to 1 AU (DS ≥ 1) on crystalline cellulose, and (2) DS equal to one arbitrary unit (DS = 1) on amorphous cellulose. It is important to note that the currently reported synergistic studies were conducted using commercial cellulase preparations, suggesting that there is a need to use the expansins or expansin-like proteins with pure cellulases (CBHs or EGs), so that the mechanism of enzyme interaction can be more closely elucidated.
Another group of enzymes with a similar function to expansins are swollenins (SWO), which are mostly produced by fungi [35,43]. In contrast to expansins, SWOs have crossactivity on amorphous cellulose (CMC), lichenan, mixed-linked glucan (MLG), and laminarin, while on microcrystalline celluloses (e.g., Avicel, cotton, and filter paper), they display non-hydrolytic activity by remodifying its structure [35,39,43]. A SWO displayed synergy with commercial cellulase on Avicel, leading to the release of 30% more reducing sugars [35]. Zhang et al. [39] purified and used CBHI (Cel7A), or endoglucanase (EG) supplemented with β-glucosidase (β-gl) and SWO, to hydrolyze phosphoric acid swollen cellulose (PASC) and cellulose nanocrystals (CNC). Interestingly, the EG, β-gl, and SWO combination displayed higher activity than the Cel7A, β-gl, and SWO combination on PASC, while on CNC, the EG-containing mixture did not show any significant activity, while the Cel7A mixture showed higher activity. Our highest estimated DS value for the EG-containing mixture was 1.2 AU, and the Cel7A-containing mixture displayed 2.0 AU on PASC, suggesting that the Cel7A and SWO synergistic interaction was more productive. Also, the synergism between Cel7A and SWO was effective because both enzymes have proven to be effective on the crystalline substrate (CNC), while EG cleaves amorphous regions of PASC. Our estimated DS values from the published data demonstrate that the EG mixture had an estimated DS value of 1.0 AU (with a 20% conversion rate of CNC), while the Cel7A mixture had the highest DS value of 1.4 AU (with a 60% conversion rate of CNC). These observations agree with a study by Santos and co-workers, which demonstrated that SWOs reduce the crystallinity of cellulose and modify the microfibril structure to create more binding and cleavage sites for cellulolytic enzymes, in particular, CBHs [43].
It is important to note that both expansins and SWOs facilitate the activity of cellulase cocktails or pure mono-component cellulases (CBHI or EG) on regenerated amorphous cellulose (RAC) or PASC, or filter paper and/or Avicel [30,35,39,42,43]. However, further laoding high concentrations of expansins and SWO, beyond a certain threshold, on the substrate does not translate to a linear increment in overall cellulose hydrolysis. The concentration that resulted in higher cellulase and non-catalytic synergism before expansin/SWO binding saturation was 10 mg protein/g biomass [17,35]. At higher loadings of expansin (30 mg protein/g biomass), a reduction in biomass conversion, from a 15% (for expansin at 10 mg/g biomass) to 4.5% conversion rate, was observed [35]. In contrast, at a similar enzyme loading (30 mg protein/g biomass), SWO did not show any significantly different activity compared to the activity displayed by 10 mg protein/mg biomass [17]. These findings support our hypothesis that expansin and SWO synergism with cellulases cannot result in DS values that are below 1 AU, and in instances where there is no synergism, the DS values will be equal to 1 AU (Table 1).

Synergism between Processive Cellulolytic Enzymes; Exo-Exo or Exo-Endo Processive Synergy
There is strong evidence that supports the hypothesis that some backbone cleaving enzymes always interact synergistically, resulting in a DS that is higher than 1 AU and a soluble sugar yield that is higher than that produced by the sum of the individual enzymes during the hydrolysis of biomass [11]; for instance, cellulase enzymes that cleave microcrystalline cellulose, such as Cel6A and Cel7A (exo-exo synergy), or Cel6A, Cel7A, and EG (exo-endo synergy; see Section 3.5), generally synergize with a DS that is higher than 1 AU [11,15,20]. Exo-exo synergy between cellulases has been demonstrated by Boisset et al. [44] using bacterial microcrystalline cellulose (BMC), while other researchers from our laboratory (Enzyme Science Programme (ESP) at Rhodes University) have used Avicel, PASC, or NaOH-treated Avicel, also called regenerated amorphous cellulose (RAC), as suitable substrates [12,15,34]. Badino et al. [20] showed that the DS between GH7 (Cel7A) and GH6 (Cel6A) enzymes sourced from Hypocrea jecorina was above 2.2 AU, while the glucose yield produced by synergistic actions was four-fold higher than that produced by individual enzymes. A study by Boisset et al. [45] also found similar results using Humicola insolens GH7 (Cel7A) and GH6 (Cel6A), which showed a DS of about 4 AU between enzyme combinations, and converted more than 50% of BMC compared to 15% by Cel6A and 28% by Cel7A alone. Based on the literature and our experiences, the fungalsourced Cel7A (GH7) and Cel6A (GH6) synergism generally displays a DS above 2 AU, and will always result in a higher yield of soluble sugars than individual enzymes. In some bacterial systems, such as Clostridium phytofermentans ISDg, a model microorganism for consolidated bioprocessing (CBP), there is a different cellulase system that degrades the crystalline regions of the cellulose [47]. The CBH enzymes employed by C. phytofermentans ISDg are Cel48 and Cel9. Zhang and co-workers suggested that Cel48 attached to the substrate from the reducing end, while Cel9 (processive EG) attached from the non-reducing end. In addition, CBHs that belong to GH 48, such as Cel48, usually contain a family 3 CBM [47]. Anaerocellum thermophilum is another bacterium that has demonstrated that the GH48 and GH9 enzyme system is generally employed to hydrolyze microcrystalline cellulose via exo-exo synergy [58]. Additionally, Irwin and co-workers showed that there was higher synergism between Cel48A and Cel9A sourced from Thermobifida fusca compared to other enzymes from different GH families, such as GH5 and GH6 [59]. Cel9 from Bacillus licheniformis was shown to have a CBM3c linked to the catalytic domain, and displayed the highest activity on PASC, BMC, and filter paper. Cel9 (GH 9) from Clostridium cellulosi possesses five CBMs, namely, CBM3c, three CBMX2s, and CBM3b [60]. Zhang and co-workers proposed that the multi-modular Cel9A cleaves crystalline cellulose fibers via a "wire-walking mode", whereby the CBMX2s and CBM3b bind to the cellulose substrate, while CBM3c, linked adjacent to the catalytic domain (CD), extends the cleft up to 17 subsites. Thus, it appears that the union of Cel9 s CBM3c and its CD assists in the processive nature of this enzyme on crystalline cellulose. Therefore, Cel9A qualifies as a pseudo CBH, based on its hydrolysis mechanism on microcrystalline cellulose. The exo-exo synergy between Cel48 and Cel9A on filter paper displayed a DS value of 1.5 AU [49]. Another study similarly reported a DS of about 1.5 AU for Cel48 (GH48) and Cel9 (GH9), while the soluble sugars produced from Avicel were four-fold higher than those produced by the individual enzymes [47]. These observations support our hypothesis that there are the following two or more distinct exo-cellulase systems that result in higher DS and yields of soluble sugars: (1) one from the fungal cellulolytic systems (mostly Cel7A and Cel6A), and (2) another from the bacterial cellulolytic systems (Cel 48 and Cel 9A).
Fungal exo-exo synergy is evolutionally designed to depend on the CBHs from GH family 6 and 7, because their fusion always results in a higher DS and yield of soluble sugars [11,16]. We propose that a DS greater than 2 AU should be considered as the threshold that always results in a higher yield of soluble sugars during GH6 and GH7 CBHs synergism (Table 1). Thus, the DS threshold of 2 AU must be used as a benchmark for Cel7A and Cel6A during the formulation of enzyme cocktails. However, a lower DS threshold for the bacterial exo-exo synergy systems, of about 1.5 AU, must be achieved because Cel48A and Cel9A displayed a very high yield of soluble sugars from the crystalline cellulosic substrate at low DS values, compared to the fungal system [47,49]. Lastly, we recommend that the fungal or bacterial exo-exo synergy systems should not be mixed in one formulation of an enzyme cocktail, since this could lead to competitive inhibition (see Section 4 for details on this inhibition).

Endo-Endo Synergy
Endoglucanases of fungal origin are diverse and complex, as some strains from certain species secrete more than three endoglucanases, which are classified under GH families 5, 6, 7, 9, 12, 44, 45, and 74 [61]. Foreman et al. [61] demonstrated that the filamentous fungus Trichoderma reesei secretes large amounts of Cel7B (EGI), Cel5A (EGII), Cel12A (EGIII), Cel61A (EGIV), and Cel45A (EGV). Even though Foreman and colleagues assigned EGIV or Cel61A to GH family 61, all the enzymes that belong to this family are now classified under auxiliary activity (AA) family 9 (http://www.cazy.org/AA9.html; accessed on 30 June 2021) (see Section 3.6 for details). In addition, Wood et al. [62] purified five distinct endoglucases (EGI to EG-V, with molecular weights of between 25 kDa and 62.5 kDa) from Penicillium pinophilum and hypothesized that each endoglucanase plays a specific role during biomass degradation. Mode of hydrolysis deciphering assays of these EGs revealed that EGI only hydrolyzed cello-oligosaccharides with a degree of polymerization (DP) higher than six, while EGII, EGIII, EGIV, and EGV hydrolyzed oligosaccharides with a DP higher than two [63]. Interestingly, EGII and EGIV only produced cellobiose, even when hydrolyzing cellotriose or cellopentaose; this suggests that the two EGs also performed transglycosylation [62][63][64]. In addition, EG from Aspegillus terreus, which belongs to GH12, could only hydrolyze cello-oligosaccharides with a DP > 4 [65]. The observations regarding the different DPs required for various classes of EGs, by Bhat et al. [63] and Segato et al. [65], suggest that it is vital to test the synergism between the different classes of EGs from different GH families and determine their degrees of synergy to develop a DS threshold for the fungal EG system. The fact that EGVI, or Cel61A, is currently classified into AA family 9 is another reason why the synergism between endoglucases from different families must be considered.
Five endoglucanases of the ruminant bacterial symbionts were detected, and their molecular weights were 42, 50, 52, 53, and 101 kDa [66]. Even though the 52 and 53 kDa EGs were more active compared to three other proteins, the authors did not classify these EGs into different GH families. Some of the bacterial EGs belong to GH family 8, and are specific for the β-1,4-glycosidic bonds of cellulose or β-glucans [67]. The GH8 endoglucanase produces mostly cellotetraose, cellotriose, and cellobiose from polymeric substrates, such as CMC and β-glucan. However, on oligosaccharides such as cellopentaose, the enzyme produced mostly cellotetraose and cellobiose, while cellohexaose hydrolysis resulted in cellotetraose, cellotriose, and cellobiose. Scapin and co-workers [67] showed that the protein structure of EG from GH8 resembles an (α/α) 6 barrel and possesses a catalytic cleft with six subsites, which results in a unique mode of action that is specific for substrates with a DP > 5. Hakamada et al. [68] also reported on a 43 kDa EG from GH8 that was sourced from Bacillus circulans KSM-N257, which had a similar mode of action and structural features to the EG reported by Scapin et al. [67], except that it had higher activity towards lichenin, followed by CMC. The Thermobifida fusca strain UPMC 901 was shown to have multiple EGs using zymography, with about four bands equaling 82, 60, 35, and 30 kDa [69]. Some Bacillus sp. also produced diverse EGs with high molecular weights of 80 to 100 kDa. It is apparent that the bacterial endoglucanase system is not well classified into old categories as fungal systems are (which are grouped into EGI to EGV); however, some studies have classified them into new GH families 5, 6, 7, 8, 9, 10, 12, 44, 51, and 74 [70] (http://www.cazy.org/; accessed on 30 June 2021).
The literature shows that EGs from both fungal and bacterial systems are biochemically characterized, but there is a lack of evidence supporting their endo-endo synergy. We propose that endo-endo synergism could be vital for the effective conversion of cellulose or complex oligosaccharides from β-glucan to soluble sugars. Also, endo-endo synergy could reveal more details on the mechanism of action of EGs classified as being members of GH family 51 and 74, which are not well characterized. Lastly, there is a lack of information regarding the role of cellulolytic enzymes from GH 61 (now classified in AA family 9) with EGs in cellulase synergy during lignocellulose saccharification.

Endoglucanase (EG)-β-Glucosidase (β-gl) Synergy
Even though there is a lack of information regarding the synergism between endoglucanases, there are several studies that demonstrate that EGs and β-gl do synergize during cellulose degradation. This information is vital because it can be used to define the threshold of DS between EG and β-gl interactions, which can result in the generation of higher yields of soluble sugars from biomass hydrolysis. Endoglucanase sourced from Scytalidium thermophilum synergized with β-gl from Humicola insolens and released more than a two-fold higher amount of reducing sugars from DMSO-pretreated sugarcane bagasse [71]. The DS of about 1.6 AU demonstrated that the EG and β-gl synergized; however, the authors only used one enzyme combination (3 EG: 1 β-gl) for their synergy studies. An EG-to-β-gl cocktail was also used to produce cellulose nanocrystals (CNC); in this case, their combinations were not well defined [72]. Additionally, Teixeira et al. [72] used the synergy of Pyrococcus horikoshii endoglucanase and Pyrococcus furiosus β-glucosidase to achieve conversion yields of 60% from wet-disk-milled Celish ® KY and 40% from unbleached Kraft wood pulp. The EG and β-gl synergism was superior to the commercial β-gl (referred to as OptimashTMBG) because it converted more than two-fold more of the amorphous regions of the biomass to soluble sugars, compared to commercial β-gl, resulting in shorter CNC. Cairo et al. [51] also demonstrated that the synergism between EG and β-gl resulted in a DS of about 9.79, and about a 10-fold higher conversion rate of CMC to soluble sugars when compared to the use of individual enzymes. Theoretically, true EG and β-gl synergy always results in DS values higher than 1 AU, and yields of soluble sugars that are higher compared to when the enzymes are acting individually. This theoretical observation is supported by the fact that the EGs produce various oligosaccharides that have a DP > 2, and β-gls cleave the EG products to glucose [12,51,73,74]. Thus, we suggest that for as long as a true EG (non-processive EG) hydrolyzes amorphous cellulose, the β-gl will produce glucose from these products, resulting in high DS and yields of soluble sugars. Several studies have demonstrated the co-evolution of the EG and β-gl interaction by co-expression of their genes in engineered fungi or bacteria to produce a protein that displays a dual function of EG and β-gl [73,74]. Thus, we propose that the DS threshold of 1.3 AU will always result in high yields of soluble sugars between true EG and β-gl.

Exo-Endo and/or β-Glucosidase Synergy (Cellulolytic Enzyme Cocktail)
Based on both theoretical and empirical evidence, CBHs, EGs, and β-gls can synergize during the hydrolysis of biomass (Figure 2), due to their varied substrate specificities. The synergism occurs in cases where the reaction contains the following: (1) CBHI and CBHII, which attack the biomass from the reducing and non-reducing ends of cellulose chains to produce cellobiose, respectively; (2) true EGs, which randomly cleave the amorphous regions of the biomass to produce oligosaccharides with a DP between 2 and 10, or cleave the amorphous regions of the microcrystalline cellulose to produce new chain ends for CBHs; (3) β-gls, which degrade the CBHs and EG-generated products to monosaccharides (glucose) [75]. Boisset et al. [44] demonstrated that CBHI (Cel7A), CBHII (Cel6A), and EG (Cel45A) sourced from Humicola insolens displayed synergism during the hydrolysis of bacterial cellulose ribbons (microcrystalline cellulose). The compositions of Cel6A and Cel7A were varied between 0 and 100%, with dosage variation steps of 20%, while Cel45A was fixed at 1.2%. The cocktail composed of the three enzymes converted more than 60% of the bacterial cellulose ribbons and displayed a DS of about 4.5 AU, while the combination of Cel6A and Cel7A converted only 30% and had a DS of 2 AU [45]. Some of the commercial cocktails, such as Celluclast ® 1.5L, consist of Trichoderma reesei cellulases in the following combinations: 55% CBH-7A, 10% CBH-6A, 10% EG-7B, 10% EG-5A, and 1% β-gl [76]. The highest biomass conversion that the authors could achieve was about a 50% yield, only when they added an extra 10% β-gl, but their experimental design did not factor in the DS, which shows the molecular-level interactions between the enzymes. The addition of 50 U/g substrate of a GH12 EG sourced from Gloeophyllum trabeum to 10 U/g substrate of Celluclast ® 1.5L resulted in a synergistic effect of about 14.5%, 16.1%, 29.0%, and 13.4% on filter paper, hydrogen peroxide-acetic acid-pretreated pine, corn stover, and rice straw, respectively [77].
The synergy between CBHI (Cel7A) and EG (Cel7B) from T. reesei effectively hydrolyzed steam-pretreated spruce and displayed a DS of about 2 AU throughout the time of the reaction [55]. Woodward et al. [53] demonstrated synergism between cellulases with a final concentration of 20 µg/mL, consisting of 5 µg of EG II/mL, 10 µg of CBH I/mL, and 5 µg of CBH II/mL, and a fixed β-gl concentration of 1.08 units. The DS reached a peak at 2 AU and, at the same point, the soluble sugars reached the highest yield [53]. Bhat and Bhat [64] also suggested that both EGIII and EGV, sourced from Penicillium pinophilum, were required for the complete hydrolysis of crystalline cellulose because they increased the DS values and yields of soluble sugars during synergy with CBHI and CBHII. The DS appears to be time dependent, as was revealed in the study between CBHI (Cel7A) and EG (Cel7B), i.e., at 0-12 h, the enzyme combinations displayed the highest DS of about 2.35 AU, followed by a decrease to 1.6 AU at 12-24 h, and the DS values of 1.6 AU did not change significantly between 24 and 48 h [14]. In the same study, Pellegrini and collaborators demonstrated that the size of the biomass affects the DS between cellulases, because they only hydrolyzed the finer fibers of the filter paper and left the thick fibers untouched. Zhang and Zhang [5] also echoed the same views that the DS between enzymes, which results in higher yields of soluble sugars, is affected by the reaction conditions and the size of the fibers they hydrolyzed. Based on the above information, we propose that the fungal cellulolytic enzyme cocktails have various DS threshold values based on the different model substrates, i.e., 4-10 AU for bacterial cellulose and cotton (DP ≥ 2000) [10,44,78], 1.5-2.5 AU for Avicel/filter paper (DP~300), and 1.3-15 AU for PASC/RAC (DP~60) [10,11,14,15].
The above information relates to the fungal cellulolytic enzyme cocktails; however, it is clear that exoglucanase systems of bacterial organisms are different to those of fungal systems (as we proposed in exo-exo synergy; Section 3.2); for instance, Cel6A, Cel48A, Cel6B, and Cel9A sourced from Cellulomonas fimi were all processive on the different types of cellulose biomasses utilized, namely, PASC, Avicel, and crystalline celluloses Iα or III from green algae [22]. It is important to note that Cel6B and Cel9A are classified as endoglucanases, but Uchiyama et al. [22] demonstrated that these enzymes are arguably highly processive, and they used HS-AFM to show that the enzymes were able to move on crystalline cellulose III. Furthermore, combinations of TrCel7A and CfCel6B resulted in high yields of soluble sugars and synergism, while TrCel6A and CfCel6B resulted in antisynergy. This observation suggests that CfCel6B is a CBH and not an EG, as it competes with CBHII (TrCel6A) for substrate binding sites. Interestingly, CfCel6B and CfCel5B demonstrated a DS of 2.16 AU and an improved yield of soluble sugar at a ratio of 3:1 when the enzyme loading was at a cellulose concentration of 2.5 mg/g [52]. However, when the enzyme concentration of CfCel6B and CfCel5B was increased to 10 mg/g cellulose, the DS values were decreased by 1 AU for all of the following combinations: 3:1, 1:1, and 1:3. The processive and non-processive EGs from GH family 9, and a CBHI from GH family 48, sourced from C. cellulosi CS-4-4, displayed synergy on filter paper when they were supplemented with β-gl from Caldicellulosiruptor sp. [49]. Various molar ratios were tested in this study, and the best combination for the formulated cellulase cocktail was 25:25:10:18 for CcCel9A: CcCel9B: CcCel48A: BlgA, which resulted in a DS of 2.6 and 1.86 mg/mL of glucose release. The literature contains more information regarding fungal cellulolytic enzyme cocktails, but, so far, only a few bacterial cellulolytic enzyme cocktails have been recorded. This suggests that there is still a lot that needs to be done with regard to bacterial cellulolytic enzymes, to establish if their synergy can be superior to that of fungi, or even commercial cocktails. We propose that a DS threshold above 1.5 AU can result in higher yields of soluble sugars for bacterial cellulolytic cocktails, but more studies are required to establish this claim.

Cellulase and LPMO Synergism
We have pointed out that some of the putative EGs that were initially deposited in GH family 61 were later moved to auxiliary activities (AA) family 9 in the CAZy database (http://www.cazy.org/; accessed on 30 July 2021). Thus, we decided to also investigate the role of AAs, particularly those from AA family 9, 10, and 16, in improving the DS during cellulose hydrolysis and the yield of soluble sugars. LPMOs (also known as AA9, AA10, or AA16) are auxiliary enzymes to cellulases, which cleave the cellulose chains via an oxidative mechanism [21,23,25]. Thus, the mechanism of substrate cleavage by LPMOs is different from cellulases classified under GH families, because they cleave the biomass via hydrolytic mechanisms [5,10]. The DS between these two classes of enzymes can reveal whether they cooperate or act individually during the hydrolysis of the cellulosic biomass.
Synergism between an LPMO and CBHI was demonstrated by core expression of the genes that encode for these proteins in the fungus (Penicillium funiculosum) that produced the LPMO naturally and five engineered strains [21]. It was evident that the secretome from the engineered strains (referred to as Pf OAO1 and Pf OA3) contained higher Avicelase activity, filter paper activity, and production of hydrogen peroxide, which suggest the presence of cellulases and LPMOs. In addition, the use of PfOAO1 and Pf OA3 for the saccharification of acid-pretreated wheat straw increased the conversion to 80% and 75%, respectively. This suggested that the overexpression of CBHI and LPMO in Pf OAO1 and Pf OA3 increased the synergy between these enzymes. Interestingly, the individual application of three LPMOs from Thermoascus aurantiacus (TaAA9A), Lentinus similis (LsAA9A), and Thielavia terrestris (TtAA9E) on Avicel or PASC resulted in no release of detectable reducing sugars, according to Tokin et al. [37]. In contrast, the LPMOs synergized with TrCel7A and TrCel6A on Avicel or PASC, except for TrCel7A and TtAA9E, which displayed anti-synergy (more details regarding anti-synergy will be provided in the next section). The binary synergy between TrCel6A and TtAA9E resulted in the highest DS of about 2.5 and 2 AU during PASC and Avicel degradation, respectively. The high synergism between these enzymes can be explained by the fact that TtAA9E cleaves cellulose substrates by oxidizing the C1 of sugar moieties. TrCel6A also synergized with TaAA9A or LsAA9A, as their DS values were 2.0 and 1.8 AU on PASC, and 1.5 and 1.3 AU on Avicel, respectively [37]. The CBHI and TaAA9A or LsAA9A displayed a DS of 2.25 and 1.4 AU on PASC, or 1.8 and 1.5 AU, respectively [37]. EG (Cel7B) and MtLPMO9A synergy improved the hydrolysis of Avicel, bacterial cellulose, and sugarcane bagasse (SCB), which resulted in higher glucose production in the presence of Anβ-gl. The DS values between cellulases (TrCel7A, TrCel6A, and TrCel7B) and MtLPMO9A during the degradation of Avicel, BC, and SCB were 2.8, 2.5, and 2.6 AU, respectively. We propose that DS values greater than 1.5 AU should be considered as a threshold that results in higher yields of soluble sugars during the synergism between cellulases and LPMOs from AA9 (LPMO9A) during the degradation of PASC. However, DS values greater than 2 AU between cellulases and LPMO9A should be considered as a threshold that results in the improved degradation of crystalline substrates, such as Avicel. Unfortunately, to date (30 July 2021), the LPMOs from family 16 (AA16) are not yet well characterized, and only a few studies have been conducted; hence, we cannot conclusively provide predictions of the DS threshold [25].
Bacterial LPMOs that possess cellulolytic activity are generally classified under AA family 10 (http://www.cazy.org/AA10.html; accessed on 30 July 2021). Similarly, to the fungal LPMO system, which can oxidize C1 or C4, bacterial LPMOs can cleave glycosidic bonds by oxidizing C1 and C4 of the cellulose chain [23,79] (http://www.cazy.org/AA10. html; accessed on 30 July 2021). Forsberg et al. [79] demonstrated that Streptomyces coelicolor was equipped with LPMOs, which strictly oxidize C1 (ScLPMO10C, also known as CelS2), and a C1/C4 oxidizer (ScLPMO10B). Interestingly, ScLPMO10C and ScLPMO10B synergy improved the degradation of PASC. The enzyme combinations produced 140 µM glucose/gluconic acid compared to 20 µM glucose/gluconic acid produced by the individual enzymes, or 400 µM cellobiose compared to 100 µM produced by the individual enzymes. Even though Forsberg et al. [79] did not measure the DS values for synergism during these studies, our calculations from their glucose/gluconic acid data suggest that a DS value of about 2.3 AU was achieved. In addition, Tuveng et al. [26] also demonstrated that LPMO10 synergized with Cel48A and Cel6B, displaying a DS value of about 1.2 AU during Avicel hydrolysis. It is worth noting that, to date, most studies have focused on fungal cellulase-LPMO synergism-but few have focused on bacterial cellulase-LPMO synergism. Therefore, we cannot speculate on the DS threshold regarding bacterial cellulase-LPMO synergism. These observations suggest that more work still needs to be conducted to establish molecular interactions between bacterial cellulases and LPMOs, or fungal cellulases and the bacterial LPMO system.

Unraveling Cellulolytic Enzyme Anti-Synergism
There are the following three cases whereby cellulolytic enzymes do not synergize and the DS values are equal to, or less than, 1 AU (DS = 1 AU): (1) the first case is where the DS values are equal to 1 AU because the cellulases, or cellulases and LPMOs do not cooperate during the reaction, meaning that each enzyme is working on its own; (2) the DS values can be less than 1 AU because of cellulases, or cellulases and non-hydrolytic enzymes that compete for the same substrate binding sites, in order for enzymes to cleave the substrate; (3) the DS values can also be lower than, or equal to, 1 AU because the rate-limiting cellulases/LPMOs and the non-rate-limiting cellulases/LPMOs in the reaction display different rates of activity on their respective components of the substrate (e.g., binary synergism between processive EG and the CBHII). A detailed discussion on the factors that lead to reduced DS thresholds or synergism between cellulosic enzymes continues below, and a summary of this can be found in Table 2. Table 2. Factors resulting in anti-synergy between cellulolytic enzymes.

Factors Effects on Enzyme Synergy Enzyme Combination DS Reference
Substrate surface EG generated uneven layers of chains on cellulose that inhibit processive CBHs.
Sometimes the DS is lower than 1 AU, but the yield of soluble sugars produced by the combined cellulases, or cellulases and LPMOs, is higher compared to when the enzymes act individually. This type of anti-synergy occurs between cellulase enzymes such as CBHII (non-rate-limiting enzyme) and EG (rate-limiting enzymes) during the hydrolysis of the crystalline substrate. Anderson et al. [79] used H. insolens EG (Cel45A) and CBHII (Cel6) to show that the rate-limiting Cel45A converted about 1.7% of Avicel, while the non-ratelimiting Cel6A converted 0.1% of Avicel after 52 h. The combination of Cel6A and Cel45A displayed the highest DS of about 0.9 AU after 3 h and 0.5 AU after 52 h [79]. This study demonstrated that the endoglucanase was producing oligosaccharides at a higher rate than CBHII was able to convert the dextrans to cellobiose. The addition of P. brasillianum IBT 20888 β-gl (Cel3A) to a Cel6A/Cel45A mixture did not improve the conversion rate of Avicel (0.4%) or DS (0.5 AU), suggesting that the enzymes were not cooperating during Avicel hydrolysis. Liu et al. [52] demonstrated that the synergy between CBHII (Cel6B) and EGs is largely "endo-driven" because they are involved in creating shortened cellulose chains that can be hydrolyzed by the processive Cel6B. This study corroborates the work conducted by Anderson et al. [79], suggesting that EG is rate limiting, while CBHII is non-rate limiting during their synergic action on cellulose. Lastly, the EG created more sites that could be cleaved by CBHII on the cellulose III allomorph (similar to PASC) compared to Avicel, explaining why the EG and CBHII synergized on PASC, and did not synergize on Avicel [7,79].
Jalak et al. [16] demonstrated that during the hydrolysis of bacterial cellulose by the EG (TrCel5A) and CBHI (TrCel7A), the EG was the rate-limiting enzyme because it removed the amorphous regions on the bacterial cellulose substrate. The amorphous regions of the biomass led to CBHI stalling and becoming redundant; however, the removal of the amorphous regions by EG resulted in higher efficiency of CBHI [16]. A simulation study by Shang and Chu [80] also revealed that the cellulose hydrolyzed by EGs can contain uneven microfiber chains, which can stall CBHI activity, which results in anti-synergy. We argue that EG has the following two distinct functions as a cellulose hydrolysis rate-limiting enzyme during its synergism with CBHs: (1) EG produces the cleavage sites for CBHII on the microcrystalline substrates, and (2) EG removes the amorphous regions for CBHI. These observations suggest that EG is very important during cellulose hydrolysis because, as a rate-limiting enzyme, it determines or improves the efficiency of CBHs (Figure 2; Malgas et al. [3]). It is worth noting that EGI and EGII compete for the same cleavage sites on the cellulose biomass, which leads to inhibition (anti-synergy) towards Avicel hydrolysis [53,81]. Bhat and Bhat [64] reported that the enzyme EGII possesses a CBM, while EGI consists of only a CD. Therefore, the processive nature of EGII leads to antisynergy due to a "traffic jam" at the site of cleavage between the two EGs. Additionally, Thoresen et al. [12] showed that enzyme cocktails formulated with CsCel48A, microbial Cel6A, and various endoglucanases from different GH families displayed anti-synergy, which resulted in reduced soluble sugars, glucose, and DS values below 1 AU.
In addition, Tokin et al. [37] and Mafa et al. [15] demonstrated that CBHI from GH7 could be inhibited/hindered by cellulolytic enzymes with the same regiospecificity, leading to DS values below 1 AU. The observed DS values below 1 AU were not only due to the same substrate binding sites, but because the enzyme mixture of the CBHI (GH7) and LPMO (AA9E), or CBHI (GH7) and exoglucanase (Exg-D: GH5_38) did not effectively hydrolyze the crystalline substrate (Avicel), and regenerated amorphous cellulose (RAC). To investigate this enzyme obstruction phenomenon, Mafa et al. [15] used a different substrate, β-(1,3-1,4)-mixed linked glucan (MLG), which could be hydrolyzed by CBHI (GH7), CBHII (GH6), and Exg-D (GH5_8). The enzyme mixture of the CBHI and Exg-D continued to show anti-synergy, with DS values below 1 AU. Even though the CBHII and Exg-D mixture did not show synergism on Avicel and RAC, the enzyme mixture showed synergism on MLG, suggesting that substrate chemistry and cleavage site availability were indeed the cause of the reduced synergism between CBHII and Exg-D. Therefore, it is apparent that the enzyme obstruction/inhibition on the same substrate cleavage site could lead to enzymes displaying anti-synergy, with DS values below 1 AU. Zhou et al. [56] also demonstrated that there is some competition between CBHI and LPMO sourced from Myceliophthora thermophila (MtLPMO9L), which belongs to C1-oxidizing LPMOs ( Figure  3 and Table 2). Zhou and co-workers demonstrated this competition by performing sequential synergy, where the MtLPMO9L was applied on the substrate first, and, after 5 h, it was removed by centrifugation at 13,000× g; the remaining substrate was boiled to deactivate LPMO, and then 5 mg/mL CBHI was added to the reactions [56]. This approach did not only prevent competition between two enzymes, but it also resulted in a 1.6-fold higher yield during the hydrolysis of Avicel. The findings by Keller et al. [29] also echo the same anti-synergistic patterns between CBHI and LPMO (MtLPMO9A) that were reported by other studies [37,56]. Thus, we propose that 0.6 ≤ DS values < 1.0 demonstrates non-cooperation between enzymes (meaning that each cellulolytic enzyme works independently), while 0.5 ≥ DS values > 0 demonstrates competition of cellulolytic enzymes for the same binding or cleavage sites on the cellulose substrates.

Perspectives and Future Recommendations
Cellulolytic enzyme synergism is essential for the deconstruction or modification of cellulose substrates in the biorefinery sector. Hence, research on synergistic interactions has attracted extensive interest from all over the world over the past few decades. Cellulases were shown to be the main enzymes that degrade or modify cellulose substrates; however, an accumulation of information on the modes of action and mechanisms of expansins, SWOs and LPMOs has improved our understanding of cellulolytic enzyme synergism. Hence, we attempted to classify the various synergisms between non-cellulolytic active enzymes and cellulolytic active enzymes, and to theorize their respective DS thresholds. The LPMOs and cellulases, or expansins/SWOs and cellulases synergisms result in enhanced conversion rates of cellulose to soluble sugars, and aldoronic and gem-diol acids. Interestingly, to date, there have been no studies that have investigated the synergy between expansins/SWOs and cellulolytic LPMOs. This suggests that a knowledge gap exists in how expansins, or SWOs and LPMOs (AA9, AA10, and AA16) interact during biomass deconstruction. The enzyme cocktail(s) containing expansins/SWOs, cellulases, and LPMOs could perhaps be superior to the currently available commercial cellulolytic enzyme cocktails. It is vital to understand the different forms of synergy between the enzymes before formulating enzyme cocktails containing expansins/SWOs, cellulases, and LPMOs. This literature review explored the possible types of synergy and grouped them into two major classes. The first class of synergism resulted from non-hydrolytic and catalytically active cellulolytic enzyme interactions, and the second class resulted from the interaction of hydrolytic enzymes. It was clear that the enzymes displayed different interactions during exo-exo, endo-glucosidase, exo-endo, and cellulase-LPMO synergy, and that the varied DS values validated the proposal of different enzyme interactions. As a result of these various DS values, we have theorized DS value thresholds, which indicated the synergism between the cellulolytic enzymes that resulted in higher yields, as a benchmark for each type of synergy in Sections 3 and 4. Lastly, based on the published data to date, we have also theorized the DS value thresholds for anti-synergism due to non-cooperating enzymes (1 ≥ DS values > 0.6), or competition for the binding/cleavage site (0.5 ≥ DS values > 0) of the cellulose substrates. This literature review can assist studies focusing on producing nano-microcrystalline cellulose through expansins/SWOs and endoglucanases, endo-endo, and endo-β-gl synergism. In addition, the theorized DS value thresholds can be applied to understand the level of synergism between cellulolytic enzymes during the formulation of enzyme cocktails.
We propose a new way of formulating cellulolytic enzyme cocktails effectively, as described below. Firstly, develop the fungal CBH core enzyme set (CES) between CBHI from GH family 7 and CBHII from GH families 5, 6, or 9; or develop the bacterial CBH-CES between CBHI from GH family 48 and CBHII from GH families 5, 6, or 9. It is important to note that CBHIs from family 7 and 48 must be kept separate at all times because their combination in the CES formulations leads to non-productive competition. Perhaps cues may also be found in the CAZy database, which showed that GH family 7 consists mostly of fungal enzymes, about 5220 entries, while the bacterial enzymes only make up 10 entries, and only four entries are unclassified (http://www.cazy.org/GH7 _eukaryota.html; accessed on 31 July 2021). GH family 48 consists mostly of bacterial enzymes, constituting 1181 entries, followed by 34 fungal enzyme entries and 4 viral enzyme entries (http://www.cazy.org/GH48_eukaryota.html; accessed on 31 July 2021).
The best performing combinations of fungal or bacterial CBH-CES (supplemented with 10% constant β-gl, relative to CBHs concentrations) must be identified and used in combination with LPMOs to formulate a quaternary CES. This means that C1-and C4oxidizing LPMOs (from AA9, AA10, or AA16) binary synergism must be established first via differential % mass protein loading or molar ratio. Therefore, the best CBH-CES and LPMO-CES can be mixed using different enzyme molar ratios to formulate a quaternary CES. Similarly, the binary or ternary synergism between EGs that are not processive from GH family 5,6,7,9,12,44,45,51, and 74 must be formulated. We caution against the use of the processive EGs in cellulolytic enzyme cocktails intended to produce a high glucose yield because they compete with non-processive EGs for the cleavage sites on cellulose, resulting in "traffic jamming". The best EG combinations for binary or ternary synergism must be used to formulate a quaternary or quinary synergy combination with LPMOs (C1-or/and C2-oxidizing enzymes). Finally, the CBH and LPMO-CES, and the EG and LPMOs can be used to formulate a cellulolytic enzyme cocktail that is effective, based on the theoretic concept that the cocktail will contain most of the enzyme required for effective cellulose hydrolysis. Expansin or SWO addition to the proposed cellulolytic enzyme cocktail must be at 1% and 10% of the total concentration because expansins, in particular, inhibit the current commercial cellulase cocktails at a higher concentrations. Institutional Review Board Statement: The current study did not use humans or animals.

Informed Consent Statement:
The current study did not use humans or animals.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to reasons of privacy.