# Generation of a Library of Carbohydrate-Active Enzymes for Plant Biomass Deconstruction

Vânia Cardoso, Joana Brás, Inês Costa, Luís Ferreira, Luís Gama, Renaud

Vincentelli, Bernard Henrissat, Carlos Fontes

► To cite this version:

Vânia Cardoso, Joana Brás, Inês Costa, Luís Ferreira, Luís Gama, et al.. Generation of a Library of Carbohydrate-Active Enzymes for Plant Biomass Deconstruction. International Journal of Molecular Sciences, 2022, 23 (7), pp.4024. 10.3390/ijms23074024. hal-03781841

HAL Id: hal-03781841

https://hal.inrae.fr/hal-03781841

Submitted on 20 Sep 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Generation of a Library of Carbohydrate-Active Enzymes for Plant Biomass Deconstruction

Vânia Cardoso <fnref n="1" />,<fnref n="2" />,* ID, Joana L. A. Brás <fnref n="2" />, Inês F. Costa <fnref n="2" />, Luís M. A. Ferreira <fnref n="1" />, Luís T. Gama <fnref n="1" /> ID,
Renaud Vincentelli <fnref n="3" />,<fnref n="4" /> ID, Bernard Henrissat <fnref n="3" />,<fnref n="4" />,<fnref n="5" /> ID and Carlos M. G. A. Fontes <fnref n="1" />,<fnref n="2" />,*

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Cardoso, V.; Brás, J.L.A.; Costa, I.F.; Ferreira, L.M.A.; Gama, L.T.; Vincentelli, R.; Henrissat, B.; Fontes, C.M.G.A. Generation of a Library of Carbohydrate-Active Enzymes for Plant Biomass Deconstruction. *Int. J. Mol. Sci.* **2022**, *23*, 4024. https://doi.org/10.3390/ijms23074024

Academic Editors: Marco Moracci,
Beatrice Cobucci-Ponzano and
Andrea Strazzulli

Received: 2 February 2022

Accepted: 3 April 2022

Published: 5 April 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_4.png" /></figure>

Copyright: © 2022 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** In nature, the deconstruction of plant carbohydrates is carried out by carbohydrate-active enzymes (CAZymes). A high-throughput (HTP) strategy was used to isolate and clone 1476 genes obtained from a diverse library of recombinant CAZymes covering a variety of sequence-based families, enzyme classes, and source organisms. All genes were successfully isolated by either PCR (61%) or gene synthesis (GS) (39%) and were subsequently cloned into *Escherichia coli* expression vectors. Most proteins (79%) were obtained at a good yield during recombinant expression. A significantly lower number (p < 0.01) of proteins from eukaryotic (57.7%) and archaeal (53.3%) origin were soluble compared to bacteria (79.7%). Genes obtained by GS gave a significantly lower number (p = 0.04) of soluble proteins while the green fluorescent protein tag improved protein solubility (p = 0.05). Finally, a relationship between the amino acid composition and protein solubility was observed. Thus, a lower percentage of non-polar and higher percentage of negatively charged amino acids in a protein may be a good predictor for higher protein solubility in *E. coli*. The HTP approach presented here is a powerful tool for producing recombinant CAZymes that can be used for future studies of plant cell wall degradation. Successful production and expression of soluble recombinant proteins at a high rate opens new possibilities for the high-throughput production of targets from limitless sources.

**Keywords:** carbohydrate-active enzymes (CAZymes); plant biomass; gene synthesis; PCR; high-throughput (HTP) cloning; HTP expression

## 1. Introduction

Plant biomass is the most abundant renewable source of organic carbon on Earth. The plant cell wall (PCW), its major constituent, is a highly heterogeneous and complex macromolecular structure that surrounds and protects the cell, thus playing a central role in plant survival [1]. PCW is composed of different types of recalcitrant polysaccharides, notably cellulose (40–50%), hemicellulose (25–35%), and lignin (15–20%) [2]. Microbial degradation of plant structural carbohydrates to generate usable sugars is a key step in the carbon cycle. It is of considerable biotechnological importance, particularly for the conversion of lignocellulosic biomass to biofuels, a vital objective for today’s society [3]. Thus, understanding the biochemistry of plant cell wall deconstruction is not only of biological importance but also has growing industrial significance [4–6].

Deconstruction of PCW carbohydrates is mediated by a large repertoire of microbial enzymes, generally termed carbohydrate-active enzymes (CAZymes). CAZymes comprise

<footnote n="1">¹ Centro de Investigação Interdisciplinar em Sanidade Animal—Faculdade de Medicina Veterinária, Universidade de Lisboa, Pólo Universitário do Alto da Ajuda, Avenida da Universidade Técnica, 1300-477 Lisboa, Portugal; luis.ferreira@reitoria.ulisboa.pt (L.M.A.F.); ltgama@fmv.ulisboa.pt (L.T.G.)</footnote>

<footnote n="2">² NZYTech Ltd., Estrada do Paço do Lumiar, Campus do Lumiar, 1649-038 Lisboa, Portugal; joana.bras@nzytech.com (J.L.A.B.); ines.costa@nzytech.com (I.F.C.)</footnote>

<footnote n="3">³ Centre National de la Recherche Scientifique, Unité Mixte de Recherche 7257, Université Aix-Marseille, 13288 Marseille, France; renaud.vincentelli@univ-amu.fr (R.V.); bernard.henrissat@gmail.com (B.H.)</footnote>

<footnote n="4">⁴ Institut National de la Recherche Agronomique, Unité sous Contrat 1408 Architecture et Fonction des Macromolécules Biologiques, 13288 Marseille, France</footnote>

<footnote n="5">⁵ Department of Biological Sciences, King Abdulaziz University, Jeddah 21589, Saudi Arabia</footnote>

<footnote>* Correspondence: vania.cardoso@nzytech.com (V.C.); carlos fontes@nzytech.com (C.M.G.A.F.)</footnote>
a diverse set of enzyme classes, including glycoside hydrolases (GHs), polysaccharide
lyases (PLs), carbohydrate esterases (CEs), and glycosyl transferases (GTs). While PLs,
CEs, and GHs carry out the breakdown of polysaccharides, GTs are mainly involved
in the formation of the glycosidic bond and thus in the biosynthesis of complex car-
bohydrates [7]. The classification of CAZymes in families has been continuously up-
dated in the Carbohydrate-Active EnZymes database (CAZy; www.cazy.org, accessed on
31 January 2022) [8,9]. The CAZy database is an invaluable research tool with respect to the
ever-increasing amount of genomic and metagenomic information relating to carbohydrate
metabolism. Currently, the database lists 172 sequence-based families of GHs, 42 families of
PLs, 19 families of CE s, 114 families of GTs, and, finally, 88 families of CBMs (as of January
2022). Recently, the CAZy database incorporated a new category of proteins, named auxil-
iary activities (AAs), which covers redox enzymes that act cooperatively with CAZymes.
The AA category of CAZymes includes lytic polysaccharide monooxygenases (LPMOs) and
ligninolytic enzymes [10]. CAZymes are classified in families according to their amino acid
sequence similarities, which reflect a common structural fold and catalytic mechanism.

In general, CAZymes have evolved complex molecular architectures containing one
or more catalytic domains generally linked to more than one accessory non-catalytic
carbohydrate-binding module (CBM) [11–13]. CBMs play an important role in CAZyme
function by promoting the hydrolysis of insoluble substrates by the adjoined catalytic
domains. CBMs bind structural carbohydrates, bringing the associated catalytic domains
into close proximity of the recalcitrant substrates and thus targeting the enzyme to their
specific substrates [13,14]. Similarly to CAZymes, CBMs are also classified in families
according to their primary sequence similarity [9].

It is now well established that the functional diversity of CAZymes is enormous and
reflects the wide multiplicity of glycan structures found in nature. Although the CAZy
database provides a solid base for carbohydrate research, the existence of a tangible library
of CAZymes to promote carbohydrate bioengineering research and related applications is
still lacking. Here, a high-throughput (HTP) strategy was used to produce a comprehensive
library of highly diverse, soluble, and pure CAZymes with known biochemical properties
and to cover the maximal number of CAZy families and Enzyme Commission numbers
(EC number). This strategy was also established previously to identify novel activities in a
range of CAZymes with unknown functions [15]. A bioinformatic-informed selection of
characterized CAZymes was initially employed, followed by HTP gene synthesis, cloning,
expression, and purification of recombinant proteins. This approach not only leads to the
generation of a large library of recombinant CAZymes of significant scientific interest to
study PCW deconstruction but also to their applicability in the industrial and bioprocessing
sectors [16].

**2. Results and Discussion**

**2.1. *Identification and Selection of CAZymes***

Initially, a bank of 1955 genes encoding biochemically characterized CAZymes was
constructed. All gene sequences were analyzed using a gene alignment algorithm to exclude
highly similar sequences, explore the CAZy family's diversity, and include a broad range
of EC numbers. As a result, 1476 enzymes of several classes, families, and EC numbers
were selected. The list of proteins included enzymes from 412 different source organisms,
predominantly from bacteria but also from the other domains: Bacteria (1435 targets),
Eukaryote (26 targets), and Archaea (15 targets). The 1476 en-zymes comprised 568 GHs,
58 PLs, 510 CE s, 14 AAs, 304 CBMs, and 22 others (containing different catalytic domains)
(Figure 1). Some sequences were designed to include just a single module while others
were maintained in their native form and all the modules were expressed. All predicted
signal peptides were removed from the target sequences.
<figure><img src="image_5.png" /><figcaption>**Figure 1.** The distribution of the 1476 enzymes produced in this study by the 5 CAZy classes.</figcaption></figure>

## 2.2. *Gene Production and Cloning*

All 1476 genes were successfully produced by amplification from genomic DNA (PCR) or by gene synthesis (GS). In total, 61% of the genes were isolated by PCR amplification using the corresponding genomic DNA as a template and the remaining 39% of the genes was codon-optimized and synthesized. Both PCR and GS products were cloned in pHTP expression vectors, encoding an N-terminal Histidine tag for protein purification. A total of 1350 genes, including all gene-encoding enzymes from the GH, CE, PL, and AA families, and some CBMs, were cloned in pHTP1 expression vector. The remaining 126 genes, encoding for CBMs, were cloned in pHTP9 expression vector, which introduced an N-terminal GFP tag into them.

## 2.3. *Protein Expression and Putative Factors Affecting Solubility*

The 1476 recombinant enzymes were expressed in *E. coli* BL21(DE3) and purified through IMAC using an HTP automated workstation. The molecular integrity of recombinant proteins was evaluated by SDS-PAGE. The data revealed that, in general, the molecular weight (Mw) of the purified enzymes was highly similar to the calculated theoretical value (Figure 2). In addition, 79% of the proteins were obtained in the soluble form (>20 mg of purified protein per liter of media). Since all the 1476 enzymes had previously been biochemically characterized, a higher number of targets at higher levels was expected. However, because only 1 expression condition was tested for all 1476 targets (same expression host, growth medium, induction conditions, and purification method) and that the threshold for solubility (>20 mg/L) was high, it is reasonable to expect that a proportion of ~20% of the enzymes requires further optimization of the expression conditions or a bigger culture scale to allow recovery of sufficient protein for biochemical characterization. Nevertheless, the effect of protein origin, gene production strategy, fusion with solubility tags, protein molecular weight, and amino acid composition on protein expression and solubility was evaluated.
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Protein expression and purity analysis of CAZymes exemplified by SDS-PAGE. Lane M contains NZYTech Low Molecular Weight (LMW) Protein Marker; Lanes 1–8 contain purified recombinant proteins with the following accession numbers: (1) ABD81807.1, (2) AAK06049.1, (3) AAG04556.1, (4) CAA84537.1, (5) CAC83072.1, (6) CAA84537.1, (7) CAB55348.1, and (8) AAP09638.1.</figcaption></figure>

## 2.4. *The Effect of Protein Origin on Recombinant Production*

The data presented in Figure 3A revealed that a significantly lower number ($p < 0.01$) of enzymes originating from eukaryotic or archaeal domains were obtained during this study. Comparison of the recombinant protein solubility across the 3 domains of life showed that a higher number of proteins from bacterial origin, 79.7% (1143 out of 1435), were obtained compared to the eukaryotic and archaeal targets, with this percentage dropping to 57.7% (15 out of 26) and 53.3% (8 out of 15), respectively (Figure 3A). The higher production rate found for bacterial targets is not surprising considering that the expression system selected was *E. coli* [17]. For example, it is widely accepted that the presence of native disulfide bonds in eukaryotic proteins represents an obstacle for proper protein folding and soluble expression in *E. coli* [17,18], resulting in the formation of insoluble inclusion bodies [19]. In addition, archaeal targets are usually poorly expressed in *E. coli* because of their origin in extreme environments and differences in host codon usage. A variety of archaeal targets originated from hyperthermophilic microorganisms, whose gene sequences usually contain a high proportion of rare codons that are poorly processed in *E. coli* [20,21]. Thus, it is expected that the expression of archaeal targets in *E. coli* results in significant levels of insoluble inclusion bodies [22].

## 2.5. *Influence of the Gene Production Strategy, Gene Synthesis, or PCR on Protein Production*

The effect of the gene production strategy on the yield of recombinant proteins was analyzed. The data, presented in Figure 3B, revealed that lower numbers of targets were obtained from genes synthesized artificially ($p = 0.04$). Several studies have demonstrated that codon optimization for heterologous production in *E. coli* can potentially increase recombinant protein expression, especially in sequences of eukaryotic and archaeal origin [23–26]. However, other studies revealed that *E. coli* codon optimization may lead to overexpression of recombinant proteins, thus resulting in a higher tendency to generate insoluble inclusion bodies [27,28]. The formation of inclusion bodies in recombinant expression systems occurs as a result of an imbalance between *in vivo* protein folding and aggregation [29–31]. The data presented here revealed that although the number of obtained proteins was reduced when genes were obtained synthetically, a significant proportion of non-viable targets revealed high levels of expression in the form of inclusion bodies (data not shown). Therefore, we can conclude that *E. coli* codon optimization increased the levels of recombinant protein expression. However, for some proteins, higher levels of expression were beyond *E. coli*'s capacity to provide proper folding and led to the accumulation of the
recombinant proteins in the form of insoluble aggregates. These targets may require further testing of the expression conditions (different strains and simpler media, such as LB).

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Efficacy of the recombinant production of CAZymes. (A) The three domains of life (proportion of recombinant proteins obtained among eukaryotic, bacteria, and archaea targets); (B) production strategy (GS, genes synthetized with codon-optimization for *E. coli*; PCR, genes obtained from genomic DNA); (C) expression vectors: pHTP1, expression vector encoding an N-terminal histidine tag; pHTP9, vector containing an N-terminal GFP tag; (D) protein's molecular mass distribution. SOL: soluble protein, INS: insoluble protein.</figcaption></figure>

## 2.6. Does the GFP Tag Promote Protein Production?

Many studies have focused on optimizing the production processes for recombinant proteins, with the aim of reducing the accumulation of insoluble inclusion bodies. In contrast, solubilization and refolding of protein from inclusion bodies is a common strategy although it requires denaturing conditions and a subsequent renaturing step, usually resulting in poor soluble protein recoveries [32,33]. Different approaches have been developed to prevent the accumulation of inclusion bodies in *E. coli*, such as the optimization of culture conditions, co-expression with molecular chaperones [34,35], lower growth temperature during gene induction [35,36], induction expression in early-log phase culture [37], and induction with lower levels of inducer concentration, such as Isopropyl β-D-1-thiogalactopyranoside (IPTG) [38]. In addition, several fusion tags have been developed to increase the solubility of overexpressed proteins, although with variable degrees of success. Currently, available fusion systems include maltose-binding polypeptide (MBP), glutathione S-transferase (GST), ubiquitin (SUMO), thioredoxin (TrxA), and Green Fluorescent Protein (GFP) [39]. The use of a fusion partner may increase the solubility of the protein and may significantly contribute to an increase in the expression yields under different conditions. The data presented here revealed a higher number (*p* = 0.05) of proteins obtained in the soluble form when fused with GFP tag (pHTP9) (Figure 3C) compared with proteins that only included a histidine tag (pHTP1). Although this suggests higher solubility through the expression with the GFP fusion, only a controlled study in which the same proteins are expressed in the fusion with the two different fusion tags would allow a more insightful elucidation of the impact of GFP on protein solubility.
## 2.7. Influence of the Protein Molecular Weight on Protein Production

No relationship was observed between the solubility of the recombinant proteins and protein molecular weight (Figure 3D). The absence of a relationship between these two factors is not surprising. Usually, proteins with a molecular mass below 100 kDa are well tolerated in *E. coli* and are therefore expressed to considerable levels. In contrast, proteins above 100 kDa are usually not properly processed in the bacterium and hence are degraded or subjected to premature termination [39]. Other authors also suggested that protein size can affect secretion performance as large cytoplasmic proteins are physically impossible to translocate [40,41]. Since only 10 proteins in this study (Figure 3D) had a molecular mass higher than 105 kDa, it is not surprising that the protein molecular weight had no effect on production. In contrast, smaller proteins, below 10 kDa, are difficult to express stably in *E. coli* because of their improper folding and higher tendency for proteolytic degradation. These proteins can be stabilized when expressed as fusions to large proteins, such as MBP, GST, and GFP, which allow them to fold properly [39]. For the smaller CBMs that were directly cloned into pHTP9 expression vector (121 out of 304 CBMs), which contains an N-terminal GFP tag, this problem was attenuated.

## 2.8. Influence of the Primary Sequence Composition on Protein Production

Previous studies have shown a relationship between amino acid composition and recombinant protein production. Therefore, the amino acid frequency for proteins expressed in *E. coli* using the described HTP approach was analyzed. Comparison of the amino acid frequency for proteins obtained or not obtained at significant levels revealed differences in the composition of non-polar amino acids (*p* = 0.002), negatively charged amino acids (*p* < 0.001), and the content of specific amino acids (see Table 1) The non-polar amino acid group comprises glycine, cysteine, alanine, leucine, isoleucine, valine, methionine, proline, phenylalanine, and tryptophan. The negatively charged amino acid group, also named acidic amino acids, contains aspartic acid and glutamic acid. The data suggest that, in *E. coli*, a lower frequency of non-polar amino-acids and higher percentage of acidic amino acids increases the probability of a protein being soluble. This result is in agreement with Bertone et al. [42] and Christendat et al. [43], who verified that protein solubility is significantly influenced by the frequency of acidic amino acids, basic amino acids, and non-polar amino acids. Niu et al. [44] confirmed the importance of negatively charged residues. Moreover, the authors observed that dipeptides comprising acidic amino acids and basic amino acids, especially dipeptides comprising acidic amino acids, were the major determinants of protein solubility [44]. The relative content of negatively charged residues seems to be the strongest determinant of protein's solubility as it was selected as a key attribute in several studies [42,43,45].

Analysis of the amino acid frequencies in proteins obtained with and without success revealed that lysine, tryptophan, arginine, alanine, aspartic acid, glutamic acid, glycine, and tyrosine (Table 1) are important in discriminating between insoluble and soluble proteins. In addition, the effect of aspartic acid, glutamic acid, and glutamine residues on protein solubility was previously described by Bertone et al. [42]. However, the present analysis revealed that other amino acids, in particular the relative content of lysine, tyrosine, and arginine, seem to be related to the propensity of a protein to be obtained, with a confirmed statistically significant effect (*p* < 0.002).

Over the past two decades, different studies have explored the relationship between protein solubility and amino acid sequence composition. Solubility prediction can increase the overall success rate of experiments by avoiding potentially insoluble targets and exclusively choosing promising candidates. Therefore, further analysis of the dipeptide content should be performed to develop a bioinformatics tool for CAZymes solubility prediction, based on primary sequence analysis. An accurate theoretical prediction of the solubility from a sequence is instrumental for setting priorities on targets in large-scale proteomics projects. The data presented here suggest that the predominance of specific amino acids is a good predictor of increased solubility of specific recombinant proteins. Considering the
large number of CAZyme variants available in the CAZy database, we can envision that the selection of the best variants for recombinant production should consider those that present a more adequate amino acid composition for the soluble expression in *E. coli*.

<table><caption>**Table 1.** Analysis of the amino acid composition in primary sequences of soluble (SOL) and insoluble (INS) CAZymes produced in this study.</caption><thead><tr><th></th><th>INS</th><th>SOL</th><th>SEM</th><th>p-Value</th></tr></thead><tbody><tr><td>(n)</td><td>310</td><td>1166</td><td></td><td></td></tr><tr><td>Amino acid groups (%)</td><td></td><td></td><td></td><td></td></tr><tr><td>Non-polar amino acids <sup>1</sup></td><td>51.4 <sup>a</sup></td><td>50.2 <sup>b</sup></td><td>0.31</td><td>0.002</td></tr><tr><td>Polar neutral amino acids <sup>2</sup></td><td>26.0</td><td>26.0</td><td>0.38</td><td>0.976</td></tr><tr><td>Negatively charged amino acids <sup>3</sup></td><td>11.2 <sup>b</sup></td><td>12.2 <sup>a</sup></td><td>0.17</td><td>0.001</td></tr><tr><td>Positively charged amino acids <sup>4</sup></td><td>11.4</td><td>11.5</td><td>0.17</td><td>0.420</td></tr><tr><td>Amino acid contents (%)</td><td></td><td></td><td></td><td></td></tr><tr><td>Isoleucine (I)</td><td>5.0</td><td>5.2</td><td>0.11</td><td>0.051</td></tr><tr><td>Leucine (L)</td><td>7.6</td><td>7.4</td><td>0.15</td><td>0.213</td></tr><tr><td>Lysine (K)</td><td>4.4 <sup>b</sup></td><td>5.0 <sup>a</sup></td><td>0.15</td><td>0.001</td></tr><tr><td>Methionine (M)</td><td>1.9</td><td>2.0</td><td>0.06</td><td>0.138</td></tr><tr><td>Phenylalanine (F)</td><td>4.1</td><td>4.1</td><td>0.08</td><td>0.828</td></tr><tr><td>Threonine (T)</td><td>6.2</td><td>6.0</td><td>0.14</td><td>0.231</td></tr><tr><td>Tryptophan (W)</td><td>2.3 <sup>a</sup></td><td>2.2 <sup>b</sup></td><td>0.07</td><td>0.042</td></tr><tr><td>Valine (V)</td><td>6.5</td><td>6.5</td><td>0.10</td><td>0.610</td></tr><tr><td>Arginine (R)</td><td>4.7 <sup>a</sup></td><td>4.3 <sup>b</sup></td><td>0.12</td><td>0.005</td></tr><tr><td>Histidine (H)</td><td>2.3</td><td>2.3</td><td>0.07</td><td>0.659</td></tr><tr><td>Alanine (A)</td><td>9.0 <sup>a</sup></td><td>8.5 <sup>b</sup></td><td>0.20</td><td>0.014</td></tr><tr><td>Asparagine (N)</td><td>5.1</td><td>5.2</td><td>0.15</td><td>0.511</td></tr><tr><td>Aspartic acid (D)</td><td>6.0 <sup>b</sup></td><td>6.5 <sup>a</sup></td><td>0.10</td><td>0.001</td></tr><tr><td>Cysteine (C)</td><td>1.0</td><td>1.0</td><td>0.07</td><td>0.784</td></tr><tr><td>Glutamic acid (E)</td><td>5.2 <sup>b</sup></td><td>5.6 <sup>a</sup></td><td>0.12</td><td>0.002</td></tr><tr><td>Glutamine (Q)</td><td>3.6</td><td>3.4</td><td>0.09</td><td>0.287</td></tr><tr><td>Glycine (G)</td><td>8.9 <sup>a</sup></td><td>8.6 <sup>b</sup></td><td>0.13</td><td>0.016</td></tr><tr><td>Proline (P)</td><td>5.0</td><td>4.9</td><td>0.10</td><td>0.418</td></tr><tr><td>Serine (S)</td><td>6.6</td><td>6.5</td><td>0.15</td><td>0.246</td></tr><tr><td>Tyrosine (Y)</td><td>4.5 <sup>b</sup></td><td>4.9 <sup>a</sup></td><td>0.11</td><td>0.001</td></tr></tbody></table>

SEM: Standard error of the mean; Means in the same line with different letter superscripts (a, b) are significantly different (p < 0.05) or tend to be significantly different (p < 0.1). ¹ The non-polar amino acid group is composed of glycine (G), cysteine (C), alanine (A), leucine (L), isoleucine (I), valine (V), methionine (M), proline (P), phenylalanine (F), and tryptophan (W); ² The polar neutral amino acid group is composed of asparagine (N), glutamine (Q), serine (S), threonine (T), and tyrosine (Y); ³ The negatively charged amino acid group, also named acidic amino acids, comprises aspartic acid (D) and glutamic acid (E); ⁴ Positively charged amino acids, also named basic amino acids, comprise arginine (R), histidine (H), and lysine (K).

## 2.9. CAZymes Activity

To confirm that the library of soluble recombinant CAZymes displayed the expected activities, a subset of the CAZymes produced was functionally tested. Thus, all 16 xylanases (EC3.2.1.8) and 11 glucuronoyl esterases (EC 3.1.1.-) of the bank of 1166 soluble proteins were functionally tested using their target substrates, under standard conditions. The data revealed that all 27 recombinant CAZymes displayed the expected activity, thus suggesting that all enzymes of the bank are active when expressed in their recombinant form and reflecting their biochemical properties described in the literature.

The influence of enzyme activity (as expressed by the EC number) on protein solubility was evaluated. As described in Table 2, the biochemical properties of the 1476 CAZymes produced in this study covered 153 different EC number activities. The activity profile of the soluble proteins was compared with the original bank, and it was observed that 17 enzyme activities (EC numbers 1.1.3.-, 1.14.99.54, 2.4.1.140, 2.4.1.20, 2.4.1.216, 2.4.1.8, 3.2.1.109, 3.2.1.111, 3.2.1.120, 3.2.1.136, 3.2.1.157, 3.2.1.164, 3.2.1.177, 3.2.1.154, 3.2.1.61, 4.2.2.16, 4.2.2.24) were absent in the bank of the 1166 soluble enzymes. Most of these activities
were represented only by one member (only two with three members) and thus the limited number of sequences per EC number limits the statistical analysis. The number of enzyme activities covered by the bank is depicted in Table 2.

<table><caption>**Table 2.** Activity of the 1166 soluble recombinant enzymes organized by CAZy family and expressed through the EC number.</caption><thead><tr><th>CAZy Family</th><th>EC Number</th></tr></thead><tbody><tr><td>AA1</td><td>1.10.3.2</td></tr><tr><td>AA3</td><td>1.1.3.12</td></tr><tr><td>AA7</td><td>1.1.3.-</td></tr><tr><td>AA10</td><td>1.-.-.-/1.14.99.54</td></tr><tr><td>AA NC</td><td>1.10.3.-/1.3.3.5</td></tr><tr><td>CE1</td><td>3.1.1.73</td></tr><tr><td>CE2</td><td>3.1.1.72</td></tr><tr><td>CE3</td><td>3.1.1.72</td></tr><tr><td>CE4</td><td>3.2.1.8/3.5.1.-</td></tr><tr><td>CE6</td><td>3.1.1.72</td></tr><tr><td>CE7</td><td>3.1.1.41/3.1.1.72</td></tr><tr><td>CE8</td><td>3.1.1.11</td></tr><tr><td>CE9</td><td>3.5.1.25</td></tr><tr><td>CE11</td><td>3.5.1.-</td></tr><tr><td>CE12</td><td>3.1.1.-/3.1.1.72</td></tr><tr><td>CE14</td><td>3.5.1.-/3.5.1.89</td></tr><tr><td>CE15</td><td>3.1.1.-/3.1.1.72</td></tr><tr><td>GH1</td><td>3.2.1.-/3.2.1.21/3.2.1.23/3.2.1.25/3.2.1.37/3.2.1.74/3.2.1.85/3.2.1.86</td></tr><tr><td>GH2</td><td>3.2.1.23/3.2.1.25/3.2.1.31/3.2.1.165</td></tr><tr><td>GH3</td><td>3.2.1.21/3.2.1.37/3.2.1.45/3.2.1.52/3.2.1.74/3.2.1.120</td></tr><tr><td>GH4</td><td>3.2.1.20/3.2.1.22/3.2.1.67/3.2.1.86/3.2.1.122/3.2.1.139</td></tr><tr><td>GH5</td><td>3.2.1.4/3.2.1.8/3.2.1.73/3.2.1.74/3.2.1.78/3.2.1.91/3.2.1.123/3.2.1.132/3.2.1.151</td></tr><tr><td>GH6</td><td>3.2.1.4</td></tr><tr><td>GH8</td><td>3.2.1.4/3.2.1.73/3.2.1.132/3.2.1.156</td></tr><tr><td>GH9</td><td>3.2.1.-/3.2.1.4/3.2.1.91/3.2.1.151/3.2.1.165</td></tr><tr><td>GH10</td><td>3.2.1.4/3.2.1.8</td></tr><tr><td>GH11</td><td>3.2.1.8</td></tr><tr><td>GH12</td><td>3.2.1.4/3.2.1.151</td></tr><tr><td>GH13</td><td>2.4.1.4/2.4.1.7/2.4.1.18/2.4.1.19/2.4.1.25/3.2.1.1/3.2.1.4/3.2.1.10/3.2.1.20/3.2.1.41/3.2.1.68/3.2.1.70/3.2.1.93/3.2.1.98/3.2.1.133/3.2.1.135/3.2.1.141/5.4.99.11/5.4.99.15/5.4.99.16</td></tr><tr><td>GH14</td><td>3.2.1.2</td></tr><tr><td>GH15</td><td>3.2.1.3/3.2.1.28/3.2.1.70</td></tr><tr><td>GH16</td><td>3.2.1.-/3.2.1.4/3.2.1.6/3.2.1.39/3.2.1.73/3.2.1.81/3.2.1.83/3.2.1.103/3.2.1.178</td></tr><tr><td>GH17</td><td>2.4.1.-</td></tr><tr><td>GH18</td><td>3.2.1.-/3.2.1.14/3.2.1.96</td></tr><tr><td>GH19</td><td>3.2.1.14</td></tr><tr><td>GH20</td><td>3.2.1.52</td></tr><tr><td>GH23</td><td>3.2.1.17/4.2.2.n1</td></tr><tr><td>GH24</td><td>3.2.1.17</td></tr><tr><td>GH25</td><td>3.2.1.17</td></tr><tr><td>GH26</td><td>3.2.1.78/3.2.1.100</td></tr><tr><td>GH27</td><td>3.2.1.88/3.2.1.94</td></tr><tr><td>GH28</td><td>3.2.1.15/3.2.1.67/3.2.1.82</td></tr><tr><td>GH29</td><td>3.2.1.51/3.2.1.111</td></tr><tr><td>GH30</td><td>3.2.1.8/3.2.1.31/3.2.1.38/3.2.1.136/3.2.1.164</td></tr><tr><td>GH31</td><td>3.2.1.-/3.2.1.20/3.2.1.84/2.4.1.161/3.2.1.177</td></tr><tr><td>GH32</td><td>3.2.1.26/3.2.1.64/3.2.1.65/3.2.1.80/3.2.1.153/4.2.2.16</td></tr><tr><td>GH33</td><td>3.2.1.-/3.2.1.18/2.4.1.-</td></tr><tr><td>GH35</td><td>3.2.1.23/3.2.1.165</td></tr><tr><td>GH36</td><td>3.2.1.22/3.2.1.49</td></tr><tr><td>GH37</td><td>3.2.1.28</td></tr><tr><td>GH38</td><td>3.2.1.-/3.2.1.24</td></tr><tr><td>GH39</td><td>3.2.1.37</td></tr></tbody></table>
<table><caption>**Table 2. Cont.**</caption><thead><tr><th>CAZy Family</th><th>EC Number</th></tr></thead><tbody><tr><td>GH42</td><td>3.2.1.-/3.2.1.23</td></tr><tr><td>GH43</td><td>3.2.1.-/3.2.1.37/3.2.1.55/3.2.1.99</td></tr><tr><td>GH44</td><td>3.2.1.4</td></tr><tr><td>GH45</td><td>3.2.1.4</td></tr><tr><td>GH46</td><td>3.2.1.132</td></tr><tr><td>GH47</td><td>3.2.1.113</td></tr><tr><td>GH48</td><td>3.2.1.4/3.2.1.176</td></tr><tr><td>GH49</td><td>3.2.1.11/3.2.1.95</td></tr><tr><td>GH50</td><td>3.2.1.23/3.2.1.81</td></tr><tr><td>GH51</td><td>3.2.1.55</td></tr><tr><td>GH52</td><td>3.2.1.37</td></tr><tr><td>GH53</td><td>3.2.1.89</td></tr><tr><td>GH55</td><td>3.2.1.39</td></tr><tr><td>GH57</td><td>2.4.1.18/2.4.1.25/3.2.1.1/3.2.1.41/3.2.1.54</td></tr><tr><td>GH62</td><td>3.2.1.55</td></tr><tr><td>GH63</td><td>3.2.1.20/3.2.1.84/3.2.1.170</td></tr><tr><td>GH64</td><td>3.2.1.39</td></tr><tr><td>GH65</td><td>2.4.1.8/2.4.1.64/2.4.1.216/2.4.1.230/2.4.1.279/2.4.1.282</td></tr><tr><td>GH66</td><td>3.2.1.11</td></tr><tr><td>GH67</td><td>3.2.1.139</td></tr><tr><td>GH68</td><td>2.4.1.9/2.4.1.10/3.2.1.26</td></tr><tr><td>GH70</td><td>2.4.1.5/2.4.1.140/2.4.4.-</td></tr><tr><td>GH73</td><td>3.2.1.-</td></tr><tr><td>GH74</td><td>3.2.1.-/3.2.1.4</td></tr><tr><td>GH75</td><td>3.2.1.132</td></tr><tr><td>GH76</td><td>3.2.1.101</td></tr><tr><td>GH77</td><td>2.4.1.25</td></tr><tr><td>GH78</td><td>3.2.1.40</td></tr><tr><td>GH79</td><td>3.2.1.31</td></tr><tr><td>GH80</td><td>3.2.1.132</td></tr><tr><td>GH81</td><td>3.2.1.39</td></tr><tr><td>GH82</td><td>3.2.1.157</td></tr><tr><td>GH84</td><td>3.2.1.35/3.2.1.52/3.2.1.169</td></tr><tr><td>GH85</td><td>3.2.1.96</td></tr><tr><td>GH86</td><td>3.2.1.81/3.2.1.178</td></tr><tr><td>GH87</td><td>3.2.1.59/3.2.1.61</td></tr><tr><td>GH88</td><td>3.2.1.-</td></tr><tr><td>GH91</td><td>4.2.2.18</td></tr><tr><td>GH92</td><td>3.2.1.-/3.2.1.24/3.2.1.113</td></tr><tr><td>GH94</td><td>2.4.1.-/2.4.1.20/2.4.1.49</td></tr><tr><td>GH95</td><td>3.2.1.51/3.2.1.63</td></tr><tr><td>GH97</td><td>3.2.1.3/3.2.1.20/3.2.1.22</td></tr><tr><td>GH98</td><td>3.2.1.102</td></tr><tr><td>GH99</td><td>3.2.1.130</td></tr><tr><td>GH100</td><td>3.2.1.26</td></tr><tr><td>GH101</td><td>3.2.1.97</td></tr><tr><td>GH102</td><td>4.2.2.n1</td></tr><tr><td>GH103</td><td>4.2.2.n1</td></tr><tr><td>GH104</td><td>4.2.2.n1</td></tr><tr><td>GH105</td><td>3.2.1.-/3.2.1.172</td></tr><tr><td>GH106</td><td>3.2.1.40</td></tr><tr><td>GH107</td><td>3.2.1.-</td></tr><tr><td>GH108</td><td>3.2.1.17</td></tr><tr><td>GH109</td><td>3.2.1.49</td></tr><tr><td>GH110</td><td>3.2.1.-/3.2.1.22</td></tr><tr><td>GH111</td><td>3.2.1.-</td></tr><tr><td>GH112</td><td>2.4.1.211/2.4.1.247</td></tr><tr><td>GH113</td><td>3.2.1.78</td></tr><tr><td>GH114</td><td>3.2.1.109</td></tr><tr><td>GH115</td><td>3.2.1.-</td></tr></tbody></table>
<table><caption>**Table 2.** Cont.</caption><thead><tr><th>CAZy Family</th><th>EC Number</th></tr></thead><tbody><tr><td>GH116</td><td>3.2.1.21/3.2.1.37/3.2.1.52</td></tr><tr><td>GH117</td><td>3.2.1.-</td></tr><tr><td>GH118</td><td>3.2.1.81</td></tr><tr><td>GH119</td><td>3.2.1.1</td></tr><tr><td>GH120</td><td>3.2.1.37</td></tr><tr><td>GH121</td><td>3.2.1.-</td></tr><tr><td>GH122</td><td>3.2.1.20</td></tr><tr><td>GH123</td><td>3.2.1.53</td></tr><tr><td>GH125</td><td>3.2.1.-</td></tr><tr><td>GH126</td><td>3.2.1.-</td></tr><tr><td>GH127</td><td>3.2.1.185</td></tr><tr><td>GH129</td><td>3.2.1.49</td></tr><tr><td>GH130</td><td>2.4.1.281/2.4.1.319</td></tr><tr><td>GH134</td><td>3.2.1.78</td></tr><tr><td>GH137</td><td>3.2.1.31</td></tr><tr><td>GH142</td><td>3.2.1.185</td></tr><tr><td>GH143</td><td>3.2.1.185</td></tr><tr><td>PL1</td><td>4.2.2.2/4.2.2.10</td></tr><tr><td>PL2</td><td>4.2.2.2/4.2.2.9</td></tr><tr><td>PL3</td><td>4.2.2.2</td></tr><tr><td>PL4</td><td>4.2.2.-</td></tr><tr><td>PL5</td><td>4.2.2.3</td></tr><tr><td>PL6</td><td>4.2.2.-/4.2.2.3</td></tr><tr><td>PL7</td><td>4.2.2.-/4.2.2.3/4.2.2.11</td></tr><tr><td>PL8</td><td>4.2.2.1/4.2.2.5/4.2.2.12/4.2.2.20</td></tr><tr><td>PL9</td><td>4.2.2.2/4.2.2.9</td></tr><tr><td>PL10</td><td>4.2.2.2</td></tr><tr><td>PL11</td><td>4.2.2.23/4.2.2.24</td></tr><tr><td>PL12</td><td>4.2.2.8</td></tr><tr><td>PL13</td><td>4.2.2.7</td></tr><tr><td>PL15</td><td>4.2.2.-/4.2.2.3</td></tr><tr><td>PL17</td><td>4.2.2.-</td></tr><tr><td>PL18</td><td>4.2.2.3/4.2.2.11</td></tr><tr><td>PL21</td><td>4.2.2.7/4.2.2.8</td></tr><tr><td>PL22</td><td>4.2.2.6</td></tr><tr><td>PL24</td><td>4.2.2.-</td></tr></tbody></table>

## 2.10. Final Outputs

In the present study, we started with 1955 putative candidates from 486 different organisms, covering 59% and 62% of the families and EC numbers described in the CAZy database, respectively (Figure 4A,B) (www.cazy.org, accessed on 31 January 2022). Following sequence analysis, 1476 enzymes (see Table S1), comprising approximately 55% of the families and EC numbers described in the CAZy database, were selected for production (Figure 4C). From 1166 CAZymes obtained in the soluble form (Figure 4D), 732 proteins were selected to build a large and diverse CAZyme library to develop fundamental and applied research (Figure 4E). Together with a previously implemented bank, the present library includes 1023 CAZymes (Figure 4E), belonging to 218 different organisms with different origins: 993 bacterial, 21 eukaryotic, and 9 archaeal. These enzymes are distributed in 5 CAZyme classes, comprising many families, and comprising 638 GHs, 77 PLs, 137 CEs, 9 AAs, and 162 CBMs (Figure 5). Therefore, approximately 57% of all families described in the CAZy database are covered by the present recombinant library, which constitutes a valuable resource for exploiting the biological and biotechnological potential of carbohydrates.
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Number of CAZy families and ECs obtained in this study. (A) CAZy database numbers (www.cazy.org, accessed on 31 January 2022), (B) selected biochemically characterized CAZymes, (C) CAZymes produced in this study, (D) enzymes obtained in the soluble form, and (E) final Nzytech' CAZy commercial library.</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 5.** Number of CAZy classes and families represented in the current library.</figcaption></figure>

### 3. Materials and Methods

#### 3.1. *Identification and Selection of CAZymes*

CAZymes were drawn from the continuously updated CAZy database, from which 1955 biochemically characterized enzymes from distinct classes, families, and activities described in the literature were identified. The diversity of EC numbers was also explored. To tune the selection, a gene alignment algorithm (NZYTech Genes & Enzymes, Lisbon, Portugal) was used to exclude highly similar sequences, namely primary sequences with homology higher than 90%. Signal peptides were predicted using SignalP 4.1 (https://services.healthtech.dtu.dk/service.php?SignalP-5.0, accessed on 31 January 2022; [46]) and removed from all candidate sequences. Enzyme domains were predicted using dbCAN,
which provides pre-computed CAZyme sequence and annotation data for several bacterial genomes [47].

## 3.2. *Polymerase Chain Reaction (PCR) and Gene Synthesis (GS)*

Genes encoding 1476 CAZymes were obtained through PCR or GS (Table S1). When genomic DNA was commercially available (DSMZ-German Collection of Microorganisms and Cell Cultures GmbH., Braunschweig, Germany), genes were obtained by PCR. DNA amplifications were performed with Supreme NZYProof DNA Polymerase (NZYTech Genes & Enzymes, Lisbon, Portugal) using genomic DNA as the template and gene-specific primers. The PCR conditions followed the manufacturer's protocols. Following amplification, the assembled PCR products were purified using an NZYDNA Clean-up 96 well plate (NZYTech Genes & Enzymes, Lisbon, Portugal) in a Liquid Handling Robot workstation (TECAN, Freedom EVO series, Mannedorf, Switzerland). On the other hand, the primary sequences of CAZymes, for which no template was available, were obtained by de novo gene synthesis. Gene sequences were designed by back translating the protein sequences and optimizing codon usage for high levels of expression in *E. coli*, using the ATGenium codon optimization algorithm. In addition, genes were designed to ensure a Codon Adaptation Index (CAI) value higher than 0.8. Synthetic genes were produced in an HTP pipeline using optimized procedures, as described in Sequeira et al. [48,49]. The databank describing the 1476 genes and enzymes is presented in Table S1.

## 3.3. *High-Throughput Cloning, Transformation, and Sequencing*

Purified PCR and GS products were cloned into pHTP1 or pHTP9 expression vectors (NZYTech Genes & Enzymes, Lisbon, Portugal) using NZYEasy Cloning & Expression kits I and IX, respectively (NZYTech Genes & Enzymes, Lisbon, Portugal), according to the protocol reported in Turchetto et al. [50] and Duhoo et al. [51]. Following the cloning reaction, recombinant plasmids were transformed using an HTP method into *E. coli* NZY5α competent cells (NZYTech Genes & Enzymes, Lisbon, Portugal). The transformed bacteria were spread on LB-agar kanamycin 24-deep-well plates (24-DW). After an overnight incubation at 37 °C, 1 colony per transformation was picked and grown in liquid LB media supplemented with 50 µg/mL of kanamycin in 24-DW plates (5 mL) sealed with gas-permeable adhesive seals. The plasmids were purified from the bacterial pellets using an NZYMiniprep 96-well plate kit (NZYTech Genes & Enzymes, Lisbon, Portugal) on a Tecan workstation (TECAN, Freedom EVO series, Mannedorf, Switzerland). All constructs were completely sequenced in both directions to ensure 100% consistency with the gene sequences.

## 3.4. *High-Throughput Protein Expression, Purification, and Quantification*

The protein expression and purification steps were based on the protocol described by Duhoo et al. [51] with few modifications. *E. coli* BL21 (DE3) cells (NZYTech Genes & Enzymes, Lisbon, Portugal) were used to transform the recombinant pHTP clones. Recombinant strains were cultured in 24-DW with 5 mL of Auto-Induction LB medium supplemented with 50 µg/mL of kanamycin and grown at 37 °C for 4 h followed by 18 h at 25 °C. At the end of the culture, cells were harvested by centrifugation at 1500× g for 15 min at 4 °C. The cell pellets were resuspended in 1 mL of NZY Bacterial Cell Lysis Buffer (NZYTech Genes & Enzymes, Lisbon, Portugal) with 4 µg/mL of DNAse and 100 µg/mL of lysozyme. The proteins were then purified by immobilized metal ion-affinity chromatography (IMAC) and all the steps were automated on a Tecan workstation (TECAN, Freedom EVO series, Mannedorf, Switzerland) containing a vacuum manifold. Briefly, 1 mL of crude cell lysate was incubated with 200 µL of Ni²⁺ Sepharose 6 Fast flow resin (GE Healthcare, 17-5318-02) with bound nickel and then transferred into 96-well filter plates (20 µm) (Macherey-Nagel). The wells were washed twice with 1000 µL of buffer A (50 mM NaHepes, 1 M NaCl, 10 mM imidazole, 5 mM CaCl₂, pH 7.5) followed by 1 wash with 1000 µL of buffer B (NaHepes, 1 M NaCl, 35 mM imidazole, 5 mM CaCl₂,
pH 7.5), and, finally, the proteins were eluted with 300 µL of elution buffer (NaHepes, 1 M NaCl, 300 mM imidazole, 5 mM CaCl₂, pH 7.5). The integrity of purified recombinant proteins, in terms of the solubility, purity, and molecular weight, was assessed by sodium dodecyl sulphate 14% polyacrylamide gel electrophoresis (SDS-PAGE) and the protein concentration was measured using the Bradford protein assay (NZYTech, Lisbon, Portugal) against a standard curve of bovine serum albumin (BSA). Recombinant proteins' molecular weight and extinction coefficient were calculated using the ExPASy-ProtParam tool (http://www.expasy.ch/tools/protparam.html, accessed on 31 January 2022).

### 3.5. Enzyme Activity

Determination of xylanase (EC3.2.1.8) and glucuronoyl esterase (EC 3.1.1.-) activity was performed as described by Correia et al. [52] and Hüttner et al. [53], respectively.

### 3.6. Statistical Analyses

The effects of different organisms with origin, gene production strategies, vector with GFP tag, and molecular masses were analyzed with the chi-squared test. Data related to the primary sequence analysis were subjected to ANOVA according to the general linear model procedures. All statistical analyses were performed with a 95% confidence interval and conducted with SAS [54].

## 4. Conclusions

The main drawback of biomass use is related to the complexity of macromolecular polysaccharides' composition, which requires a wide plethora of CAZymes presenting a broad spectrum of specificities. The strategy used in this study combined bioinformatic tools for enzyme selection and an optimized HTP platform for the production, cloning, expression, and purification of recombinant proteins with a wide diversity of enzyme specificities. This approach resulted in a library comprising more than 1000 CAZymes involved in carbohydrate degradation, making its application in both industrial and bioprocessing sectors possible. The data presented here reinforce the observations that different factors, such as production strategies and protein primary sequence, play a major role in determining the expression and solubility of a protein heterologously produced in *E. coli*. Further analysis of the amino acid composition of the proteins selected for this study must be performed to develop innovative strategies to predict protein solubility. Nevertheless, the effectiveness of the HTP approach described here is illustrated by the high number of soluble targets obtained, most of them with high solubility rates. Taken together, these studies reveal that the implemented strategy described here is a powerful tool not only for CAZymes production but also for the generation of large libraries of recombinant proteins that will allow exploration of the biological functions of the extensive genomic and meta-genomic information available from various sources.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23074024/s1.

**Author Contributions:** Conceptualization, V.C. and C.M.G.A.F.; methodology, J.L.A.B. and R.V.; formal analysis, L.T.G.; investigation, V.C. and I.F.C.; writing—original draft preparation, V.C.; writing—review and editing, L.M.A.F., R.V., B.H. and C.M.G.A.F.; supervision, C.M.G.A.F.; project administration, J.L.A.B. and C.M.G.A.F. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by Portugal 2020, grant number 47033, project name, Glycomed.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Not applicable.
**Acknowledgments:** V. Cardoso was supported by FCT, Fundação para a Ciência e a Tecnologia, with a PhD fellowship (SFRH/BDE/101026/2014). We thank Shabir Najmudin for insightful discussions and the careful revision of the manuscript.

**Conflicts of Interest:** Luis M. A. Ferreira and Carlos M.G.A. Fontes are beneficiaries of the company that commercializes some of the enzymes described in this work.

## References

1. Caffall, K.H.; Mohnen, D. The structure, function, and biosynthesis of plant cell wall pectic polysaccharides. *Carbohydr. Res.* **2009**, *344*, 1879–1900. [CrossRef] [PubMed]

2. Gray, K.A.; Zhao, L.; Emptage, M. Bioethanol. *Curr. Opin. Chem. Biol.* **2006**, *10*, 141–146. [CrossRef] [PubMed]

3. Lynd, L.R.; Laser, M.S.; Bransby, D.; Dale, B.E.; Davison, B.; Hamilton, R.; Himmel, M.; Keller, M.; McMillan, J.D.; Sheehan, J.; et al. How biotech can transform biofuels. *Nat. Biotechnol.* **2008**, *26*, 169–172. [CrossRef] [PubMed]

4. Sticklen, M.B. Plant genetic engineering for biofuel production: Towards affordable cellulosic ethanol. *Nat. Rev. Genet.* **2008**, *9*, 433–443. [CrossRef] [PubMed]

5. Himmel, M.E.; Bayer, E.A. Lignocellulose conversion to biofuels: Current challenges, global perspectives. *Curr. Opin. Biotechnol.* **2009**, *20*, 316–317. [CrossRef]

6. Hervé, C.; Rogowski, A.; Blake, A.W.; Marcus, S.E.; Gilbert, H.J.; Knox, J.P. Carbohydrate-binding modules promote the enzymatic deconstruction of intact plant cell walls by targeting and proximity effects. *Proc. Natl. Acad. Sci. USA* **2010**, *107*, 15293–15298. [CrossRef]

7. Chakraborty, S.; Rani, A.; Dhillon, A.; Goyal, A. Polysaccharide lyases. In *Current Developments in Biotechnology and Bioengineering: Production, Isolation and Purification of Industrial Products*; Elsevier B.V.: Gurgaon, India, 2017; pp. 527–539.

8. Cantarel, B.L.; Coutinho, P.M.; Rancurel, C.; Bernard, T.; Lombard, V.; Henrissat, B. The carbohydrate-active enzymes database (CAZy): An expert resource for glycogenomics. *Nucleic Acids Res.* **2009**, *37*, 233–238. [CrossRef]

9. Lombard, V.; Golaconda Ramulu, H.; Drula, E.; Coutinho, P.M.; Henrissat, B. The carbohydrate-active enzymes database (CAZy) in 2013. *Nucleic Acids Res.* **2014**, *42*, 490–495. [CrossRef]

10. Levasseur, A.; Drula, E.; Lombard, V.; Coutinho, P.M.; Henrissat, B. Expansion of the enzymatic repertoire of the CAZy database to integrate auxiliary redox enzymes. *Biotechnol. Biofuels* **2013**, *6*, 1–14. [CrossRef]

11. Boraston, A.B.; Bolam, D.N.; Gilbert, H.J.; Davies, G.J. Carbohydrate-binding modules: Fine-tuning polysaccharide recognition. *Biochem. J.* **2004**, *382*, 769–781. [CrossRef]

12. Hashimoto, H. Recent structural studies of carbohydrate-binding modules. *Cell. Mol. Life Sci.* **2006**, *63*, 2954–2967. [CrossRef] [PubMed]

13. Guillén, D.; Sánchez, S.; Rodríguez-Sanoja, R. Carbohydrate-binding domains: Multiplicity of biological roles. *Appl. Microbiol. Biotechnol.* **2010**, *85*, 1241–1249. [CrossRef] [PubMed]

14. Gilbert, H.J. The biochemistry and structural biology of plant cell wall deconstruction. *Plant Physiol.* **2010**, *153*, 444–455. [CrossRef] [PubMed]

15. Helbert, W.; Poulet, L.; Drouillard, S.; Mathieu, S.; Loiodice, M.; Couturier, M.; Lombard, V.; Terrapon, N.; Turchetto, J.; Vincentelli, R.; et al. Discovery of novel carbohydrate-active enzymes through the rational exploration of the protein sequences space. *Proc. Natl. Acad. Sci. USA* **2019**, *116*, 6063–6068. [CrossRef] [PubMed]

16. Zeng, W.; Guo, L.; Xu, S.; Chen, J.; Zhou, J. High-throughput screening technology in industrial biotechnology. *Trends Biotechnol.* **2020**, *38*, 888–906. [CrossRef] [PubMed]

17. Camilo, C.M.; Polikarpov, I. High-throughput cloning, expression and purification of glycoside hydrolases using ligation-independent cloning (LIC). *Protein Expr. Purif.* **2014**, *99*, 35–42. [CrossRef] [PubMed]

18. Nguyen, V.D.; Hatahet, F.; Salo, K.E.H.; Enlund, E.; Zhang, C.; Ruddock, L.W. Pre-expression of a sulfhydryl oxidase significantly increases the yields of eukaryotic disulfide bond containing proteins expressed in the cytoplasm of *E.coli*. *Microb. Cell Factories* **2011**, *10*, 1–10. [CrossRef]

19. Khow, O.; Suntrarachun, S. Strategies for production of active eukaryotic proteins in bacterial expression system. *Asian Pac. J. Trop. Biomed.* **2012**, *2*, 159–162. [CrossRef]

20. Kim, S.; Lee, S.B. Rare codon clusters at 5'-end influence heterologous expression of archaeal gene in *Escherichia coli*. *Protein Expr. Purif.* **2006**, *50*, 49–57. [CrossRef]

21. Lipinszki, Z.; Vernyik, V.; Farago, N.; Sari, T.; Puskas, L.G.; Blattner, F.R.; Posfai, G.; Gyorfy, Z. Enhancing the translational capacity of *E. coli* by resolving the codon bias. *ACS Synth. Biol.* **2018**, *7*, 2656–2664. [CrossRef]

22. Kim, S.; Lee, S.B. Soluble expression of archaeal proteins in *Escherichia coli* by using fusion-partners. *Protein Expr. Purif.* **2008**, *62*, 116–119. [CrossRef]

23. Piubelli, L.; Campa, M.; Temporini, C.; Binda, E.; Mangione, F.; Amicosante, M.; Terreni, M.; Marinelli, F.; Pollegioni, L. Optimizing *Escherichia coli* as a protein expression platform to produce *Mycobacterium tuberculosis* immunogenic proteins. *Microb. Cell Factories* **2013**, *12*, 1–14. [CrossRef] [PubMed]
24. Lai, G.H.; Lin, Y.C.; Tsai, Y.L.; Lien, Y.Y.; Lin, M.K.; Chen, H.J.; Chang, W.T.; Tzen, J.T.; Lee, M.S. High yield production of pigeon circovirus capsid protein in the *E. coli* by evaluating the key parameters needed for protein expression. *BMC Vet. Res.* **2014**, *10*, 1–11. [CrossRef] [PubMed]

25. Devi, N.; Adivitiya, Khasa, Y.P. A combinatorial approach of N-terminus blocking and codon optimization strategies to enhance the soluble expression of recombinant hIL-7 in *E. coli* fed-batch culture. *Appl. Microbiol. Biotechnol.* **2016**, *100*, 9979–9994. [CrossRef]

26. Mirzaei, M.; Saffar, B.; Shareghi, B. Cloning, codon optimization, and expression of *Yersinia intermedia* phytase gene in *E. coli*. *Iran. J. Biotechnol.* **2016**, *14*, 63–69. [CrossRef] [PubMed]

27. Fahnert, B.; Lilie, H.; Neubauer, P. Inclusion bodies: Formation and utilisation. *Adv. Biochem. Eng. Biotechnol.* **2004**, *89*, 93–142. [PubMed]

28. Zhao, M.; Wu, F.; Xu, P. Development of a rapid high-efficiency scalable process for acetylated *Sus scrofa* cationic trypsin production from *Escherichia coli* inclusion bodies. *Protein Expr. Purif.* **2015**, *116*, 120–126. [CrossRef] [PubMed]

29. Mannall, G.J.; Titchener-Hooker, N.J.; Dalby, P.A. Factors affecting protein refolding yields in a fed-batch and batch-refolding system. *Biotechnol. Bioeng.* **2007**, *97*, 1523–1534. [CrossRef]

30. de Groot, N.S.; Espargaró, A.; Morell, M.; Ventura, S. Studies on bacterial inclusion bodies. *Future Microbiol.* **2008**, *3*, 423–435. [CrossRef]

31. Eiberle, M.K.; Jungbauer, A. Technical refolding of proteins: Do we have freedom to operate? *Biotechnol. J.* **2010**, *5*, 547–559. [CrossRef]

32. Sahdev, S.; Khattar, S.K.; Saini, K.S. Production of active eukaryotic proteins through bacterial expression systems: A review of the existing biotechnology strategies. *Mol. Cell. Biochem.* **2008**, *307*, 249–264. [CrossRef] [PubMed]

33. Singh, A.; Upadhyay, V.; Upadhyay, A.K.; Singh, S.M.; Panda, A.K. Protein recovery from inclusion bodies of *Escherichia coli* using mild solubilization process. *Microb. Cell Factories* **2015**, *14*, 1–10. [CrossRef] [PubMed]

34. Voulgaridou, G.P.; Mantso, T.; Chlichlia, K.; Panayiotidis, M.I.; Pappa, A. Efficient *E. coli* expression strategies for production of soluble human crystallin ALDH3A1. *PLoS ONE* **2013**, *8*, 15–17. [CrossRef] [PubMed]

35. Arya, R.; Sabir, J.S.M.; Bora, R.S.; Saini, K.S. Optimization of culture parameters and novel strategies to improve protein solubility. In *Methods in Molecular Biology (Methods and Protocols)*; Humana Press: New York, NY, USA, 2015; pp. 45–63.

36. Vera, A.; González-Montalbán, N.; Arís, A.; Villaverde, A. The conformational quality of insoluble recombinant proteins is enhanced at low growth temperatures. *Biotechnol. Bioeng.* **2007**, *96*, 1101–1106. [CrossRef]

37. San-Miguel, T.; Pérez-Bermúdez, P.; Gavidia, I. Production of soluble eukaryotic recombinant proteins in *E. coli* is favoured in early log-phase cultures induced at low temperature. *Springerplus* **2013**, *2*, 89. [CrossRef]

38. Huyen, D.T.; Giang, L.Q.; Hai, T.N. Expression of flagellin FLjB derived from *Salmonella enterica* serovar typhimurium in *Escherichia coli* BL21. *Tap Chi Sinh Hoc* **2014**, *36*, 506–514. [CrossRef]

39. Kaur, J.J.; Kumar, A.; Kaur, J.J. Strategies for optimization of heterologous protein expression in *E. coli*: Roadblocks and reinforcements. *Int. J. Biol. Macromol.* **2018**, *106*, 803–822. [CrossRef]

40. Feilmeier, B.J.; Iseminger, G.; Schroeder, D.; Webber, H.; Phillips, G.J. Green fluorescent protein functions as a reporter for protein localization in *Escherichia coli*. *J. Bacteriol.* **2000**, *182*, 4068–4076. [CrossRef]

41. Palacios, J.L.; Zaror, I.; Martínez, P.; Uribe, F.; Opazo, P.; Socías, T.; Gidekel, M.; Venegas, A. Subset of hybrid eukaryotic proteins is exported by the type I secretion system of *Erwinia chrysanthemi*. *J. Bacteriol.* **2001**, *183*, 1346–1358. [CrossRef]

42. Bertone, P.; Kluger, Y.; Lan, N.; Zheng, D.; Christendat, D.; Yee, A.; Edwards, A.M.; Arrowsmith, C.H.; Montelione, G.T.; Gerstein, M. SPINE: An integrated tracking database and data mining approach for identifying feasible targets in high-throughput structural proteomics. *Nucleic Acids Res.* **2001**, *29*, 2884–2898. [CrossRef]

43. Christendat, D.; Yee, A.; Dharamsi, A.; Kluger, Y.; Savchenko, A.; Cort, J.R.; Booth, V.; Mackereth, C.D.; Saridakis, V.; Ekiel, I.; et al. Structural proteomics of an archaeon. *Nat. Struct. Biol.* **2000**, *7*, 903–909.

44. Niu, X.; Li, N.; Chen, D.; Wang, Z. Interconnection between the protein solubility and amino acid and dipeptide compositions. *Protein Pept. Lett.* **2013**, *20*, 88–95. [CrossRef] [PubMed]

45. Goh, C.S.; Lan, N.; Douglas, S.M.; Wu, B.; Echols, N.; Smith, A.; Milburn, D.; Montelione, G.T.; Zhao, H.; Gerstein, M. Mining the structural genomics pipeline: Identification of protein properties that affect high-throughput experimental analysis. *J. Mol. Biol.* **2004**, *336*, 115–130. [CrossRef] [PubMed]

46. Petersen, T.N.; Brunak, S.; Heijne, G.V.; Nielsen, H. SignalP 4.0: Discriminating signal peptides from transmembrane regions. *Nat. Methods* **2011**, *8*, 785–786. [CrossRef]

47. Yin, Y.; Mao, X.; Yang, J.; Chen, X.; Mao, F.; Xu, Y. dbCAN: A web resource for automated carbohydrate-active enzyme annotation. *Nucleic Acids Res.* **2012**, *40*, 445–451. [CrossRef] [PubMed]

48. Sequeira, A.F.; Brás, J.L.A.; Fernandes, V.O.; Guerreiro, C.I.P.D.; Vincentelli, R.; Fontes, C.M.G.A. A novel platform for high-throughput gene synthesis to maximize recombinant expression in *Escherichia coli*. In *Methods in Molecular Biology*, 1620; Domingues, L., Ed.; Springer Science+Business Media: Berlin/Heidelberg, Germany, 2017; pp. 113–128.

49. Sequeira, A.F.; Turchetto, J.; Saez, N.J.; Peysson, F.; Ramond, L.; Duhoo, Y.; Blémont, M.; Fernandes, V.O.; Gama, L.T.; Ferreira, L.; et al. Gene design, fusion technology and TEV cleavage conditions influence the purification of oxidized disulphide-rich venom peptides in *Escherichia coli*. *Microb. Cell Factories* **2017**, *16*, 4. [CrossRef]
50. Turchetto, J.; Sequeira, A.F.; Ramond, L.; Peysson, F.; Brás, J.L.; Saez, N.J.; Duhoo, Y.; Blémont, M.; Guerreiro, C.I.; Quinton, L.; et al. High-throughput expression of animal venom toxins in *Escherichia coli* to generate a large library of oxidized disulphide-reticulated peptides for drug discovery. *Microb. Cell Factories* **2017**, *16*, 6. [CrossRef] [PubMed]

51. Duhoo, Y.; Girault, J.; Vourtsos, J.; Ramond, L.; Durbé, F.; Fourquet, P.; Nominé, Y.; Cardoso, V.; Sequeira, A.F.; Brás, J.L.A.; et al. High-throughput production of a new library of human single and tandem PDZ domains allows quantitative PDZ-peptide interaction screening through high-throughput holdup assay. In *Methods in Molecular Biology*; Springer: Berlin/Heidelberg, Germany, 2019.

52. Correia, M.A.; Mazumder, K.; Brás, J.L.; Firbank, S.J.; Zhu, Y.; Lewis, R.J.; York, W.S.; Fontes, C.M.; Gilbert, H.J. Structure and function of an arabinoxylan-specific xylanase. *J. Biol Chem.* **2011**, *286*, 22510–22520. [CrossRef]

53. Hüttner, S.; Klaubauf, S.; de Vries, R.P.; Olsson, L. Characterisation of three fungal glucuronoyl esterases on glucuronic acid ester model compounds. *Appl. Microbiol. Biotechnol.* **2017**, *101*, 5301–5311. [CrossRef]

54. SAS. SAS User's Guide: Statistics; SAS Institute: Cary, NC, USA, 2004.