Relationships of crystallinity and reaction rates for enzymatic degradation of poly (ethylene terephthalate), PET

Biocatalytic degradation of plastic waste is anticipated to play an important role in future recycling systems. However, enzymatic degradation of crystalline poly (ethylene terephthalate) (PET) remains consistently poor. Herein, we employed functional assays to elucidate the molecular underpinnings of this limitation. This included utilizing complementary activity assays to monitor the degradation of PET disks with varying crystallinity (XC), as well as determining enzymatic kinetic parameters for soluble PET fragments. The results indicate that an efficient PET-hydrolase, LCCICCG, operates through an endolytic mode of action, and that its activity is limited by conformational constraints in the PET polymer. Such constraints become more pronounced at high XC values, and this limits the density of productive sites on the PET surface. Endolytic chain-scissions are the dominant reaction type in the initial stage, and this means that little or no soluble organic product are released. However, endolytic cuts gradually and locally promote chain mobility and hence the density of attack sites on the surface. This leads to an upward concave progress curve; a behavior sometimes termed lag-phase kinetics.

ITC analysis of polydisperse systems: Unravelling the impact of sample heterogeneity

Binding interactions often involve heterogeneous samples displaying a distribution of binding sites that vary in affinity and binding enthalpy. Examples include biological samples like proteins and chemically produced samples like modified cyclodextrins. Experimental studies often ignore sample heterogeneity and treat the system as an interaction of two homogeneous species, i.e. a chemically well-defined ligand binding to one type of site. The present study explores, by simulations and experiments, the impact of heterogeneity in isothermal titration calorimetry (ITC) setups where one of the binding components is heterogeneous. It is found that the standard single-site model, based on the assumption of two homogeneous binding components, provides excellent fits to simulated ITC data when the binding free energy is normally distributed and all sites have similar binding enthalpies. In such cases, heterogeneity can easily go undetected but leads to underestimated binding constants. Heterogeneity in the binding enthalpy is a bigger problem and may result in enthalpograms of increased complexity that are likely to be misinterpreted as two-site binding or other complex binding models. Finally, it is shown that heterogeneity can account for previously observed experimental anomalies. All simulations are accessible in Google Colab for readers to experiment with the simulation parameters.

Unveiling PET Hydrolase Surface Dynamics through Fluorescence Microscopy

PET hydrolases are an emerging class of enzymes that are being heavily researched for their use in bioprocessing polyethylene terephthalate (PET). While work has been done in studying the binding of PET oligomers to the active site of these enzymes, the dynamics of PET hydrolases binding to a bulk PET surface is an unexplored area. Here, methods were developed for total internal reflection fluorescence (TIRF) microscopy and fluorescence recovery after photobleaching (FRAP) microscopy to study the adsorption and desorption dynamics of these proteins onto a PET surface. TIRF microscopy was employed to measure both on and off rates of two of the most commonly studied PET hydrolases, PHL7 and LCC, on a PET surface. It was found that these proteins have a much slower off rates on the order of 10-3  s-1 , comparable to non-productive binding in enzymes such as cellulose. In combination with FRAP microscopy, a dynamic model is proposed in which adsorption and desorption dominates over lateral diffusion over the surface. The results of this study could have implications for the future engineering of PET hydrolases, either to target them to a PET surface or to modulate interaction with their substrate.

Impact of Synergy Partner Cel7B on Cel7A Binding Rates: Insights from Single-Molecule Data

Enzymatic degradation of cellulosic biomass is a well-established route for the sustainable production of biofuels, chemicals, and materials. A strategy employed by nature and industry to achieve an efficient degradation of cellulose is that cellobiohydrolases (or exocellulases), such as Cel7A, work synergistically with endoglucanases, such as Cel7B, to achieve the complete degradation of cellulose. However, a complete mechanistic understanding of this exo-endo synergy is still lacking. Here, we used single-molecule fluorescence microscopy to quantify the binding kinetics of Cel7A on cellulose when it is acting alone on the cellulose fibrils and in the presence of its synergy partner, the endoglucanase Cel7B. To this end, we used a fluorescently tagged Cel7A and studied its binding in the presence of the unlabeled Cel7B. This provided the single-molecule data necessary for the estimation of the rate constants of association kON and dissociation kOFF of Cel7A for the substrate. We show that the presence of Cel7B does not impact the dissociation rate constant, kOFF. But, the association rate of Cel7A decreases by a factor of 2 when Cel7B is present at a molar proportion of 10:1. This ratio has previously been shown to lead to synergy. This decrease in association rate is observed in a wide range of total enzyme concentrations, from sub nM to μM concentrations. This decrease in kON is consistent with the formation of cellulase clusters recently observed by others using atomic force microscopy.

Improved Hydrolysis of Granular Starches by a Psychrophilic α-Amylase Starch Binding Domain-Fusion

Degradation of starch granules by a psychrophilic α-amylase, AHA, from the Antarctic bacterium Pseudoalteromonas haloplanktis TAB23 was facilitated by C-terminal fusion to a starch-binding domain (SBD) from either Aspergillus niger glucoamylase (SBD_GA) or Arabidopsis thaliana glucan, water dikinase 3 (SBD_GWD3) via a decapeptide linker. Depending on the waxy, normal or high-amylose starch type and the botanical source, the AHA-SBD fusion enzymes showed up to 3 times higher activity than AHA wild-type. The SBD-fusion thus increased the density of enzyme attack-sites and binding-sites on the starch granules by up to 5- and 7-fold, respectively, as measured using an interfacial catalysis approach that combined conventional Michaelis-Menten kinetics, with the substrate in excess, and inverse kinetics, having enzyme in excess, with enzyme-starch granule adsorption isotherms. Higher substrate affinity of the SBD_GA compared to SBD_GWD3 was accompanied by the superior activity of AHA-SBD_GA in agreement with the Sabatier principle of adsorption limited heterogenous catalysis.

Purification and biochemical characterization of SM14est, a PET-hydrolyzing enzyme from the marine sponge-derived Streptomyces sp. SM14

The successful enzymatic degradation of polyester substrates has fueled worldwide investigation into the treatment of plastic waste using bio-based processes. Within this realm, marine-associated microorganisms have emerged as a promising source of polyester-degrading enzymes. In this work, we describe the hydrolysis of the synthetic polymer PET by SM14est, a polyesterase which was previously identified from Streptomyces sp. SM14, an isolate of the marine sponge Haliclona simulans. The PET hydrolase activity of purified SM14est was assessed using a suspension-based assay and subsequent analysis of reaction products by UV-spectrophotometry and RP-HPLC. SM14est displayed a preference for high salt conditions, with activity significantly increasing at sodium chloride concentrations from 100 mM up to 1,000 mM. The initial rate of PET hydrolysis by SM14est was determined to be 0.004 s^-1 at 45°C, which was increased by 5-fold to 0.02 s^-1 upon addition of 500 mM sodium chloride. Sequence alignment and structural comparison with known PET hydrolases, including the marine halophile PET6, and the highly efficient, thermophilic PHL7, revealed conserved features of interest. Based on this work, SM14est emerges as a useful enzyme that is more similar to key players in the area of PET hydrolysis, like PHL7 and IsPETase, than it is to its marine counterparts. Salt-tolerant polyesterases such as SM14est are potentially valuable in the biological degradation of plastic particles that readily contaminate marine ecosystems and industrial wastewaters.

Structural and functional characterization of a multi-domain GH92 α-1,2-mannosidase from Neobacillus novalis

Many secreted eukaryotic proteins are N-glycosylated with oligosaccharides composed of a high-mannose N-glycan core and, in the specific case of yeast cell-wall proteins, an extended α-1,6-mannan backbone carrying a number of α-1,2- and α-1,3-mannose substituents of varying lengths. α-Mannosidases from CAZy family GH92 release terminal mannose residues from these N-glycans, providing access for the α-endomannanases, which then degrade the α-mannan backbone. Most characterized GH92 α-mannosidases consist of a single catalytic domain, while a few have extra domains including putative carbohydrate-binding modules (CBMs). To date, neither the function nor the structure of a multi-domain GH92 α-mannosidase CBM has been characterized. Here, the biochemical investigation and crystal structure of the full-length five-domain GH92 α-1,2-mannosidase from Neobacillus novalis (NnGH92) with mannoimidazole bound in the active site and an additional mannoimidazole bound to the N-terminal CBM32 are reported. The structure of the catalytic domain is very similar to that reported for the GH92 α-mannosidase Bt3990 from Bacteroides thetaiotaomicron, with the substrate-binding site being highly conserved. The function of the CBM32s and other NnGH92 domains was investigated by their sequential deletion and suggested that whilst their binding to the catalytic domain was crucial for the overall structural integrity of the enzyme, they appear to have little impact on the binding affinity to the yeast α-mannan substrate. These new findings provide a better understanding of how to select and optimize other multi-domain bacterial GH92 α-mannosidases for the degradation of yeast α-mannan or mannose-rich glycans.

Interfacial Catalysis during Amylolytic Degradation of Starch Granules: Current Understanding and Kinetic Approaches

Enzymatic hydrolysis of starch granules forms the fundamental basis of how nature degrades starch in plant cells, how starch is utilized as an energy resource in foods, and develops efficient, low-cost saccharification of starch, such as bioethanol and sweeteners. However, most investigations on starch hydrolysis have focused on its rates of degradation, either in its gelatinized or soluble state. These systems are inherently more well-defined, and kinetic parameters can be readily derived for different hydrolytic enzymes and starch molecular structures. Conversely, hydrolysis is notably slower for solid substrates, such as starch granules, and the kinetics are more complex. The main problems include that the surface of the substrate is multifaceted, its chemical and physical properties are ill-defined, and it also continuously changes as the hydrolysis proceeds. Hence, methods need to be developed for analyzing such heterogeneous catalytic systems. Most data on starch granule degradation are obtained on a long-term enzyme-action basis from which initial rates cannot be derived. In this review, we discuss these various aspects and future possibilities for developing experimental procedures to describe and understand interfacial enzyme hydrolysis of native starch granules more accurately.

Protein-plastic interactions: The driving forces behind the high affinity of a carbohydrate-binding module for polyethylene terephthalate

Polyethylene terephthalate (PET) waste is a common pollutant in the environment, mainly due to resistance of the plastic to bio-degradation. Nevertheless, hydrolytic enzymes have been identified with activity on this substrate, which are continually being engineered to increase activity. Some insoluble biological polymers are degraded by enzymes with a multi-domain architecture, comprising of a catalytic domain, and a substrate-binding domain, such as a carbohydrate-binding module (CBM). Enzymes that degrade PET have been shown to have a higher activity when fused with these CBMs, indicating a promising route for engineering better enzymes for plastic bioprocessing. However, no detailed study of the affinity and binding mechanism of these domains on PET has yet been made. Here, we perform an in depth analysis of a binding domain from CBM family 2 on PET, showing that the affinity of the protein for the plastic is highly dependent on temperature and crystallinity of the plastic. We also investigate the mechanism of the interaction, and show how affinity may be engineered in both directions. CBM affinity for other synthetic polymers is also demonstrated for the first time. Our results demonstrate that the substrate affinity of fusion enzymes with binding modules can be tuned to the desired level.

Rate Response of Poly(Ethylene Terephthalate) (PET)-Hydrolases to Substrate Crystallinity: Basis for Understanding the Lag Phase

The rate response of PET-hydrolases to increased substrate crystallinity (XC) of poly(ethylene terephthalate) (PET) manifests as a rate-lowering effect that varies significantly for different enzymes. Herein, we report the influence of XC on the product release rate of six thermostable PET-hydrolases. All enzyme reactions displayed a distinctive lag phase until measurable product formation occurred. The duration of the lag phase increased with XC. The recently discovered PET-hydrolase PHL7 worked efficiently on "amorphous" PET disks (XC ∼10%), but this enzyme was extremely sensitive to increased XC, whereas the enzymes LCCICCG, LCC, and DuraPETase had higher tolerance to increases in XC and had activity on PET disks having XC of 24.4%. Microscopy revealed that the XC tolerant hydrolases generated smooth and more uniform substrate surface erosion than PHL7 during reaction. Structural and molecular dynamics analysis of the PET-hydrolyzing enzymes disclosed that surface electrostatics and enzyme flexibility may account for the observed differences.

Tuning alginate β-lactoglobulin complex coacervation by modulating pH and temperature

The use of biomolecules in food matrices and encapsulation systems is, as in other areas, moving towards greener solutions and a center piece here is the complex coacervation between natural anionic polysaccharides and proteins. Both alginate and β-lactoglobulin (β-Lg) are used in different sectors and have been shown to coacervate at pH < 5.2. Albeit with increased interest, complex coacervation has almost exclusively been studied from a macromolecular perspective, and described as an interaction based on charge-charge attraction. Here, we show that through changes in pH and temperature, alginate β-Lg complex coacervation can be tuned to purpose. By detailed biophysical and chemical characterization of coacervation and coacervate particles, insights into the molecular interaction and effect of external factors are obtained. We find that carboxylate resonance stabilization causes a release of protons at pH < pK_a,alginate and an uptake of protons at pH > pK_a,alginate upon coacervation. Proton release and uptake were quantified at pH 2.65 and 4.00 by isothermal titration calorimetry to be 4 and 2 protons per β-Lg molecule, respectively. By increasing the temperature to 65 °C, we discovered a secondary β-Lg concentration dependent coacervation step, where the formed particles change into large assemblies driven by entropy. These findings bring new insights to complex coacervation and its applicability in microencapsulation and drug delivery.

Molecular insights into alginate β-lactoglobulin A multivalencies-The foundation for their amorphous aggregates and coacervation

For improved control of biomaterial property design, a better understanding of complex coacervation involving anionic polysaccharides and proteins is needed. Here, we address the initial steps in condensate formation of β-lactoglobulin A (β-LgA) with nine defined alginate oligosaccharides (AOSs) and describe their multivalent interactions in structural detail. Binding of AOSs containing four, five, or six uronic acid residues (UARs), either all mannuronate (M), all guluronate (G), or alternating M and G embodying the block structural components of alginates, was characterized by isothermal titration calorimetry, nuclear magnetic resonance spectroscopy (NMR), and molecular docking. β-LgA was highly multivalent exhibiting binding stoichiometries decreasing from five to two AOSs with increasing degree of polymerization (DP) and similar affinities in the mid micromolar range. The different AOS binding sites on β-LgA were identified by NMR chemical shift perturbation analyses and showed diverse compositions of charged, polar and hydrophobic residues. Distinct sites for the shorter AOSs merged to accommodate longer AOSs. The AOSs bound dynamically to β-LgA, as concluded from saturation transfer difference and ¹ H-ligand-targeted NMR analyses. Molecular docking using Glide within the Schrödinger suite 2016-1 revealed the orientation of AOSs to only vary slightly at the preferred β-LgA binding site resulting in similar XP glide scores. The multivalency coupled with highly dynamic AOS binding with lack of confined conformations in the β-LgA complexes may help explain the first steps toward disordered β-LgA alginate coacervate structures.

Reaction Pathways for the Enzymatic Degradation of Poly(Ethylene Terephthalate): What Characterizes an Efficient PET-Hydrolase?

Bioprocessing of polyester waste has emerged as a promising tool in the quest for a cyclic plastic economy. One key step is the enzymatic breakdown of the polymer, and this entails a complicated pathway with substrates, intermediates, and products of variable size and solubility. We have elucidated this pathway for poly(ethylene terephthalate) (PET) and four enzymes. Specifically, we combined different kinetic measurements and a novel stochastic model and found that the ability to hydrolyze internal bonds in the polymer (endo-lytic activity) was a key parameter for overall enzyme performance. Endo-lytic activity promoted the release of soluble PET fragments with two or three aromatic rings, which, in turn, were broken down with remarkable efficiency (k_cat /K_M values of about 10⁵  M^-1 s^-1 ) in the aqueous bulk. This meant that approximatly 70 % of the final, monoaromatic products were formed via soluble di- or tri-aromatic intermediates.

The Sabatier principle as a tool for discovery and engineering of industrial enzymes

The recent breakthrough in all-atom, protein structure prediction opens new avenues for a range of computational approaches in enzyme design. These new approaches could become instrumental for the development of technical biocatalysts, and hence our transition toward more sustainable industries. Here, we discuss one approach, which is well-known within inorganic catalysis, but essentially unexploited in biotechnology. Specifically, we review examples of linear free-energy relationships (LFERs) for enzyme reactions and discuss how LFERs and the associated Sabatier Principle may be implemented in algorithms that estimate kinetic parameters and enzyme performance based on model structures.

Unaided efficient transglutaminase cross-linking of whey proteins strongly impacts the formation and structure of protein alginate particles

•

Microbial transglutaminase (MTG) cross-linked >70% β-lactoglobulin (β-Lg) at pH 8.5.

•

Initial MTG catalyzed isopeptide bond formation caused partial unfolding of β-Lg.

•

>75% of whey protein cross-linked, forming hetero-polymers containing β-Lg.

•

50% less alginate is needed to form particles with cross-linked than with native β-Lg.

•

Cross-linked β-Lg and alginate formed suspendable hydrophobically driven particles.

There is a dogma within whey protein modification, which dictates the necessity of pretreatment to enzymatic cross-linking of β-lactoglobulin (β-Lg). Here microbial transglutaminase (MTG) cross-linked whey proteins and β-Lg effectively in 50 mM NaHCO₃, pH 8.5, without pretreatment. Cross-linked β-Lg spanned 18 to >240 kDa, where 6 of 9 glutamines reacted with 8 of 15 lysines. The initial isopeptide bond formation caused loss of β-Lg native structure with t_1/2 = 3 h, while the polymerization occurred with t_1/2 = 10 h. Further, cross-linking effects on protein carbohydrate interaction have been overlooked, leaving a gap in understanding of these complex food matrices. Complexation with alginate showed that β-Lg cross-linking decreased onset of particle formation, hydrodynamic diameter, stoichiometry (β-Lg/alginate) and dissociation constant. The complexation was favored at higher temperatures (40 °C), suggesting that hydrophobic interactions were important. Thus, β-Lg was cross-linked without pretreatment and the resulting polymers gave rise to altered complexation with alginate.

Sabatier Principle for Rationalizing Enzymatic Hydrolysis of a Synthetic Polyester

Interfacial enzyme reactions are common in Nature and in industrial settings, including the enzymatic deconstruction of poly(ethylene terephthalate) (PET) waste. Kinetic descriptions of PET hydrolases are necessary for both comparative analyses, discussions of structure-function relations and rational optimization of technical processes. We investigated whether the Sabatier principle could be used for this purpose. Specifically, we compared the kinetics of two well-known PET hydrolases, leaf-branch compost cutinase (LCC) and a cutinase from the bacterium Thermobifida fusca (TfC), when adding different concentrations of the surfactant cetyltrimethylammonium bromide (CTAB). We found that CTAB consistently lowered the strength of enzyme-PET interactions, while its effect on enzymatic turnover was strongly biphasic. Thus, at gradually increasing CTAB concentrations, turnover was initially promoted and subsequently suppressed. This correlation with maximal turnover at an intermediate binding strength was in accordance with the Sabatier principle. One consequence of these results was that both enzymes had too strong intrinsic interaction with PET for optimal turnover, especially TfC, which showed a 20-fold improvement of k _cat at the maximum. LCC on the other hand had an intrinsic substrate affinity closer to the Sabatier optimum, and the turnover rate was 5-fold improved at weakened substrate binding. Our results showed that the Sabatier principle may indeed rationalize enzymatic PET degradation and support process optimization. Finally, we suggest that future discovery efforts should consider enzymes with weakened substrate binding because strong adsorption seems to limit their catalytic performance.

Analysis of fungal high-mannose structures using CAZymes

Glycoengineering ultimately allows control over glycosylation patterns to generate new glycoprotein variants with desired properties. A common challenge is glycan heterogeneity, which may affect protein function and limit the use of key techniques such as mass spectrometry. Moreover, heterologous protein expression can introduce nonnative glycan chains that may not fulfill the requirement for therapeutic proteins. One strategy to address these challenges is partial trimming or complete removal of glycan chains, which can be obtained through selective application of exoglycosidases. Here, we demonstrate an enzymatic O-deglycosylation toolbox of a GH92 α-1,2-mannosidase from Neobacillus novalis, a GH2 β-galactofuranosidase from Amesia atrobrunnea and the jack bean α-mannosidase. The extent of enzymatic O-deglycosylation was mapped against a full glycosyl linkage analysis of the O-glycosylated linker of cellobiohydrolase I from Trichoderma reesei (TrCel7A). Furthermore, the influence of deglycosylation on TrCel7A functionality was evaluated by kinetic characterization of native and O-deglycosylated forms of TrCel7A. This study expands structural knowledge on fungal O-glycosylation and presents a ready-to-use enzymatic approach for controlled O-glycan engineering in glycoproteins expressed in filamentous fungi.

Structure-function analysis of two closely related cutinases from Thermobifida cellulosilytica

Cutinases can play a significant role in a biotechnology-based circular economy. However, relatively little is known about the structure-function relationship of these enzymes, knowledge that is vital to advance optimized, engineered enzyme candidates. Here, two almost identical cutinases from Thermobifida cellulosilytica DSM44535 (Thc_Cut1 and Thc_Cut2) with only 18 amino acids difference were used for a rigorous biochemical characterization of their ability to hydrolyze poly(ethylene terephthalate) (PET), PET-model substrates, and cutin-model substrates. Kinetic parameters were compared with detailed in silico docking studies of enzyme-ligand interactions. The two enzymes interacted with, and hydrolyzed PET differently, with Thc_Cut1 generating smaller PET-degradation products. Thc_Cut1 also showed higher catalytic efficiency on long-chain aliphatic substrates, an effect likely caused by small changes in the binding architecture. Thc_Cut2, in contrast, showed improved binding and catalytic efficiency when approaching the glass transition temperature of PET, an effect likely caused by longer amino acid residues in one area at the enzyme's surface. Finally, the position of the single residue Q93 close to the active site, rotated out in Thc_Cut2, influenced the ligand position of a trimeric PET-model substrate. In conclusion, we illustrate that even minor sequence differences in cutinases can affect their substrate binding, substrate specificity, and catalytic efficiency drastically.

Semi-empirical Analysis of Complex ITC Data from Protein-Surfactant Interactions

Isothermal titration calorimetry (ITC) is a widely used method to determine binding affinities and thermodynamics in ligand-receptor interactions, but it also has the capability of providing detailed information on much more complex events. However, the lack of available methods to analyze ITC data is limiting the use of the technique in such multifaceted cases. Here, we present the software ANISPROU. Through a semi-empirical approach that allows for extraction of quantitative information from complex ITC data, ANISPROU solves an inverse problem where three parameters describing a set of predefined functions must be found. In analogy to strategies adopted in other scientific fields, such as geophysics, imaging, and many others, it employs an optimization algorithm which minimizes the difference between calculated and experimental data. In contrast to the existing methods, ANISPROU provides automated and objective analysis of ITC data on sodium dodecyl sulfate (SDS)-induced protein unfolding, and in addition, more information can be extracted from the data. Here, data series on SDS-mediated protein unfolding is analyzed, and binding isotherms and thermodynamic information on the unfolding events are extracted. The obtained binding isotherms as well as the enthalpy of different events are similar to those obtained using the existing manual methods, but our methodology ensures a more robust result, as the entire data set is used instead of single data points. We foresee that ANISPROU will be useful in other cases with complex enthalpograms, for example, in cases with coupled interactions in biomolecular, polymeric, and amphiphilic systems including cases where both structural changes and interactions occur simultaneously.

Two different regimes in alcohol-induced coil-helix transition: effects of 2,2,2-trifluoroethanol on proteins being either independent of or enhanced by solvent structural fluctuations

Inhomogeneous distribution of constituent molecules in a mixed solvent has been known to give remarkable effects on the solute, e.g., conformational changes of biomolecules in an alcohol-water mixture. We investigated the general effects of 2,2,2-trifluoroethanol (TFE) on proteins/peptides in a mixture of water and TFE using melittin as a model protein. Fluctuations and Kirkwood-Buff integrals (KBIs) in the TFE-H2O mixture, quantitative descriptions of inhomogeneity, were determined by small-angle X-ray scattering investigation and compared with those in the aqueous solutions of other alcohols. The concentration fluctuation for the mixtures ranks as methanol < ethanol ≪ TFE < tert-butanol < 1-propanol, indicating that the inhomogeneity of molecular distribution in the TFE-H2O mixture is unexpectedly comparable to those in the series of mono-ols. On the basis of the concentration dependence of KBIs between the TFE molecules, it was found that a strong attraction between the TFE molecules is not necessarily important to induce helix conformation, which is inconsistent with the previously proposed mechanism. To address this issue, by combining the KBIs and the helix contents reported by the experimental spectroscopic studies, we quantitatively evaluated the change in the preferential binding parameter of TFE to melittin attributed to the coil-helix transition. As a result, we found two different regimes on TFE-induced helix formation. In the dilute concentration region of TFE below ∼2 M, where the TFE molecules are not aggregated among themselves, the excess preferential binding of TFE to the helix occurs due to the direct interaction between them, namely independent of the solvent fluctuation. In the higher concentration region above ∼2 M, in addition to the former effect, the excess preferential binding is significantly enhanced by the solvent fluctuation. This scheme should be held as general cosolvent effects of TFE on proteins/peptides.

Physical constraints and functional plasticity of cellulases

Enzyme reactions, both in Nature and technical applications, commonly occur at the interface of immiscible phases. Nevertheless, stringent descriptions of interfacial enzyme catalysis remain sparse, and this is partly due to a shortage of coherent experimental data to guide and assess such work. In this work, we produced and kinetically characterized 83 cellulases, which revealed a conspicuous linear free energy relationship (LFER) between the substrate binding strength and the activation barrier. The scaling occurred despite the investigated enzymes being structurally and mechanistically diverse. We suggest that the scaling reflects basic physical restrictions of the hydrolytic process and that evolutionary selection has condensed cellulase phenotypes near the line. One consequence of the LFER is that the activity of a cellulase can be estimated from its substrate binding strength, irrespectively of structural and mechanistic details, and this appears promising for in silico selection and design within this industrially important group of enzymes.

Surface display as a functional screening platform for detecting enzymes active on PET

Poly(ethylene terephthalate) (PET) is the world’s most abundant polyester plastic, and its ongoing accumulation in nature is causing a global environmental problem. Currently, the main recycling processes utilize thermomechanical or chemical means, resulting in the deterioration of the mechanical properties of PET. Consequently, polluting de novo synthesis remains preferred, creating the need for more efficient and bio-sustainable ways to hydrolyze the polymer. Recently, a PETase enzyme from the bacterium Ideonella sakaiensis was shown to facilitate PET biodegradation, albeit at slow rate. Engineering of more efficient PETases is required for industrial relevance, but progress is currently hampered by the dependency on intracellular expression in Escherichia coli. To create a more efficient screening platform in E. coli, we explore different surface display anchors for fast and easy assaying of PETase activity. We show that PETases can be functionally displayed on the bacterial cell surface, enabling screening of enzyme activity on PET microparticles – both while anchored to the cell and following solubilization of the enzymes.

Thermodynamic and structural study of DMPC-alkanol systems

The thermodynamic and structural behaviors of lamellar dimyristoylphosphatidylcholine-alkanol (abbreviation DMPC-CnOH, n = 8-18 is the even number of carbons in the alkyl chain) systems were studied by using DSC and SAXD/WAXD methods at a 0-0.8 CnOH : DMPC molar ratio range. Up to n≤ 10 a significant biphasic effect depending on the main transition temperature t_m on the CnOH concentration was observed. Two breakpoints were revealed: turning point (TP), corresponding to the minimum, and threshold concentration (c_T), corresponding to the end of the biphasic tendency. These breakpoints were also observed in the alkanol concentration dependent change in the enthalpy of the main transition ΔH_m. In the case of CnOHs with n > 10 we propose a marked shift of TP and c_T to very low concentrations; consequently, only increase of t_m is observed. A partial phase diagram was constructed for a pseudo-binary DMPC-C12OH system. We suggest a fluid-fluid immiscibility of the DMPC-C12OH system above c_T with a consequent formation of domains with different C12OH contents. At a constant CnOH concentration, the effects of CnOHs on ΔH_m and bilayer repeat distance were found to depend predominantly on the mismatch between CnOH and lipid chain lengths. Observed effects are suggested to be underlined by a counterbalancing effect of interchain van der Waals interactions and headgroup repulsion.

Comparative Biochemistry of Four Polyester (PET) Hydrolases*

The potential of bioprocessing in a circular plastic economy has strongly stimulated research into the enzymatic degradation of different synthetic polymers. Particular interest has been devoted to the commonly used polyester, poly(ethylene terephthalate) (PET), and a number of PET hydrolases have been described. However, a kinetic framework for comparisons of PET hydrolases (or other plastic-degrading enzymes) acting on the insoluble substrate has not been established. Herein, we propose such a framework, which we have tested against kinetic measurements for four PET hydrolases. The analysis provided values of k_cat and K_M , as well as an apparent specificity constant in the conventional units of M^-1 s^-1 . These parameters, together with experimental values for the number of enzyme attack sites on the PET surface, enabled comparative analyses. A variant of the PET hydrolase from Ideonella sakaiensis was the most efficient enzyme at ambient conditions; it relied on a high k_cat rather than a low K_M . Moreover, both soluble and insoluble PET fragments were consistently hydrolyzed much faster than intact PET. This suggests that interactions between polymer strands slow down PET degradation, whereas the chemical steps of catalysis and the low accessibility associated with solid substrate were less important for the overall rate. Finally, the investigated enzymes showed a remarkable substrate affinity, and reached half the saturation rate on PET when the concentration of attack sites in the suspension was only about 50 nM. We propose that this is linked to nonspecific adsorption, which promotes the nearness of enzyme and attack sites.

Computing Cellulase Kinetics with a Two-Domain Linear Interaction Energy Approach

While heterogeneous enzyme reactions play an essential role in both nature and green industries, computational predictions of their catalytic properties remain scarce. Recent experimental work demonstrated the applicability of the Sabatier principle for heterogeneous biocatalysis. This provides a simple relationship between binding strength and the catalytic rate and potentially opens a new way for inexpensive computational determination of kinetic parameters. However, broader implementation of this approach will require fast and reliable prediction of binding free energies of complex two-phase systems, and computational procedures for this are still elusive. Here, we propose a new framework for the assessment of the binding strengths of multidomain proteins, in general, and interfacial enzymes, in particular, based on an extended linear interaction energy (LIE) method. This two-domain LIE (2D-LIE) approach was successfully applied to predict binding and activation free energies of a diverse set of cellulases and resulted in robust models with high accuracy. Overall, our method provides a fast computational screening tool for cellulases that have not been experimentally characterized, and we posit that it may also be applicable to other heterogeneously acting biocatalysts.

Impact of Alginate Mannuronic-Guluronic Acid Contents and pH on Protein Binding Capacity and Complex Size

Alginates, serving as hydrocolloids in the food and pharma industries, form particles at pH < 4.5 with positively charged proteins, such as β-lactoglobulin (β-Lg). Alginates are linear anionic polysaccharides composed of 1,4-linked β-d-mannuronate (M) and α-l-guluronate (G) residues. The impact of M and G contents and pH is investigated to correlate with the formation and size of β-Lg alginate complexes under relevant ionic strength. It is concluded, using three alginates of M/G ratios 0.6, 1.1, and 1.8 and similar molecular mass, that β-Lg binding capacity is higher at pH 4.0 than at pH 2.65 and for high M content. By contrast, the largest particles are obtained at pH 2.65 and with high G content. At pH 4.0 and 2.65, the stoichiometry was 28-48 and 3-10 β-Lg molecules bound per alginate, respectively, increasing with higher M content. The findings will contribute to the design of formation of the desired alginate-protein particles in the acidic pH range.

Activity of fungal β-glucosidases on cellulose

BACKGROUND

Fungal beta-glucosidases (BGs) from glucoside hydrolase family 3 (GH3) are industrially important enzymes, which convert cellooligosaccharides into glucose; the end product of the cellulolytic process. They are highly active against the β-1,4 glycosidic bond in soluble substrates but typically reported to be inactive against insoluble cellulose.

RESULTS

We studied the activity of four fungal GH3 BGs on cellulose and found significant activity. At low temperatures (10 ℃), we derived the approximate kinetic parameters k _cat = 0.3 ± 0.1 s^-1 and K _M  = 80 ± 30 g/l for a BG from Aspergillus fumigatus (AfBG) acting on Avicel. Interestingly, this maximal turnover is higher than reported values for typical cellobiohydrolases (CBH) at this temperature and comparable to those of endoglucanases (EG). The specificity constant of AfGB on Avicel was only moderately lowered compared to values for EGs and CBHs.

CONCLUSIONS

Overall these observations suggest a significant promiscuous side activity of the investigated GH3 BGs on insoluble cellulose. This challenges the traditional definition of a BG and supports suggestions that functional classes of cellulolytic enzymes may represent a continuum of overlapping modes of action.

Functional analysis of chimeric TrCel6A enzymes with different carbohydrate binding modules

The glycoside hydrolase (GH) family 6 is an important group of enzymes that constitute an essential part of industrial enzyme cocktails used to convert lignocellulose into fermentable sugars. In nature, enzymes from this family often have a carbohydrate binding module (CBM) from the CBM family 1. These modules are known to promote adsorption to the cellulose surface and influence enzymatic activity. Here, we have investigated the functional diversity of CBMs found within the GH6 family. This was done by constructing five chimeric enzymes based on the model enzyme, TrCel6A, from the soft-rot fungus Trichoderma reesei. The natural CBM of this enzyme was exchanged with CBMs from other GH6 enzymes originating from different cellulose degrading fungi. The chimeric enzymes were expressed in the same host and investigated in adsorption and quasi-steady-state kinetic experiments. Our results quantified functional differences of these phylogenetically distant binding modules. Thus, the partitioning coefficient for substrate binding varied 4-fold, while the maximal turnover (kcat) showed a 2-fold difference. The wild-type enzyme showed the highest cellulose affinity on all tested substrates and the highest catalytic turnover. The CBM from Serendipita indica strongly promoted the enzyme's ability to form productive complexes with sites on the substrate surface but showed lower turnover of the complex. We conclude that the CBM plays an important role for the functional differences between GH6 wild-type enzymes.

The structural basis of fungal glucuronoyl esterase activity on natural substrates

Structural and functional studies were conducted of the glucuronoyl esterase (GE) from Cerrena unicolor (CuGE), an enzyme catalyzing cleavage of lignin-carbohydrate ester bonds. CuGE is an α/β-hydrolase belonging to carbohydrate esterase family 15 (CE15). The enzyme is modular, comprised of a catalytic and a carbohydrate-binding domain. SAXS data show CuGE as an elongated rigid molecule where the two domains are connected by a rigid linker. Detailed structural information of the catalytic domain in its apo- and inactivated form and complexes with aldouronic acids reveal well-defined binding of the 4-O-methyl-a-D-glucuronoyl moiety, not influenced by the nature of the attached xylo-oligosaccharide. Structural and sequence comparisons within CE15 enzymes reveal two distinct structural subgroups. CuGE belongs to the group of fungal CE15-B enzymes with an open and flat substrate-binding site. The interactions between CuGE and its natural substrates are explained and rationalized by the structural results, microscale thermophoresis and isothermal calorimetry.

Substrate binding in the processive cellulase Cel7A: Transition state of complexation and roles of conserved tryptophan residues

Cellobiohydrolases effectively degrade cellulose and are of biotechnological interest because they can convert lignocellulosic biomass to fermentable sugars. Here, we implemented a fluorescence-based method for real-time measurements of complexation and decomplexation of the processive cellulase Cel7A and its insoluble substrate, cellulose. The method enabled detailed kinetic and thermodynamic analyses of ligand binding in a heterogeneous system. We studied WT Cel7A and several variants in which one or two of four highly conserved Trp residues in the binding tunnel had been replaced with Ala. WT Cel7A had on/off-rate constants of 1 × 105 m-1 s-1 and 5 × 10-3 s-1, respectively, reflecting the slow dynamics of a solid, polymeric ligand. Especially the off-rate constant was many orders of magnitude lower than typical values for small, soluble ligands. Binding rate and strength both were typically lower for the Trp variants, but effects of the substitutions were moderate and sometimes negligible. Hence, we propose that lowering the activation barrier for complexation is not a major driving force for the high conservation of the Trp residues. Using so-called Φ-factor analysis, we analyzed the kinetic and thermodynamic results for the variants. The results of this analysis suggested a transition state for complexation and decomplexation in which the reducing end of the ligand is close to the tunnel entrance (near Trp-40), whereas the rest of the binding tunnel is empty. We propose that this structure defines the highest free-energy barrier of the overall catalytic cycle and hence governs the turnover rate of this industrially important enzyme.

Removal of N-linked glycans in cellobiohydrolase Cel7A from Trichoderma reesei reveals higher activity and binding affinity on crystalline cellulose

BACKGROUND

Cellobiohydrolase from glycoside hydrolase family 7 is a major component of commercial enzymatic mixtures for lignocellulosic biomass degradation. For many years, Trichoderma reesei Cel7A (TrCel7A) has served as a model to understand structure-function relationships of processive cellobiohydrolases. The architecture of TrCel7A includes an N-glycosylated catalytic domain, which is connected to a carbohydrate-binding module through a flexible, O-glycosylated linker. Depending on the fungal expression host, glycosylation can vary not only in glycoforms, but also in site occupancy, leading to a complex pattern of glycans, which can affect the enzyme's stability and kinetics.

RESULTS

Two expression hosts, Aspergillus oryzae and Trichoderma reesei, were utilized to successfully express wild-types TrCel7A (WT _Ao and WT _Tr ) and the triple N-glycosylation site deficient mutants TrCel7A N45Q, N270Q, N384Q (ΔN-glyc _Ao and ΔN-glyc _Tr ). Also, we expressed single N-glycosylation site deficient mutants TrCel7A (N45Q _Ao , N270Q _Ao , N384Q _Ao ). The TrCel7A enzymes were studied by steady-state kinetics under both substrate- and enzyme-saturating conditions using different cellulosic substrates. The Michaelis constant (K _M ) was consistently found to be lowered for the variants with reduced N-glycosylation content, and for the triple deficient mutants, it was less than half of the WTs' value on some substrates. The ability of the enzyme to combine productively with sites on the cellulose surface followed a similar pattern on all tested substrates. Thus, site density (number of sites per gram cellulose) was 30-60% higher for the single deficient variants compared to the WT, and about twofold larger for the triple deficient enzyme. Molecular dynamic simulation of the N-glycan mutants TrCel7A revealed higher number of contacts between CD and cellulose crystal upon removal of glycans at position N45 and N384.

CONCLUSIONS

The kinetic changes of TrCel7A imposed by removal of N-linked glycans reflected modifications of substrate accessibility. The presence of N-glycans with extended structures increased K _M and decreased attack site density of TrCel7A likely due to steric hindrance effect and distance between the enzyme and the cellulose surface, preventing the enzyme from achieving optimal conformation. This knowledge could be applied to modify enzyme glycosylation to engineer enzyme with higher activity on the insoluble substrates.

pH profiles of cellulases depend on the substrate and architecture of the binding region

Understanding the pH effect of cellulolytic enzymes is of great technological importance. In this study, we have examined the influence of pH on activity and stability for central cellulases (Cel7A, Cel7B, Cel6A from Trichoderma reesei, and Cel7A from Rasamsonia emersonii). We systematically changed pH from 2 to 7, temperature from 20°C to 70°C, and used both soluble (4-nitrophenyl β- d-lactopyranoside [pNPL]) and insoluble (Avicel) substrates at different concentrations. Collective interpretation of these data provided new insights. An unusual tolerance to acidic conditions was observed for both investigated Cel7As, but only on real insoluble cellulose. In contrast, pH profiles on pNPL were bell-shaped with a strong loss of activity both above and below the optimal pH for all four enzymes. On a practical level, these observations call for the caution of the common practice of using soluble substrates for the general characterization of pH effects on cellulase activity. Kinetic modeling of the experimental data suggested that the nucleophile of Cel7A experiences a strong downward shift in pK_a upon complexation with an insoluble substrate. This shift was less pronounced for Cel7B, Cel6A, and for Cel7A acting on the soluble substrate, and we hypothesize that these differences are related to the accessibility of water to the binding region of the Michaelis complex.

A practical approach to steady-state kinetic analysis of cellulases acting on their natural insoluble substrate

Measurement of steady-state rates (v_SS) is straightforward in standard enzymology with soluble substrate, and it has been instrumental for comparative biochemical analyses within this area. For insoluble substrate, however, experimental values of v_ss remain controversial, and this has strongly limited the amount and quality of comparative analyses for cellulases and other enzymes that act on the surface of an insoluble substrate. In the current work, we have measured progress curves over a wide range of conditions for two cellulases, TrCel6A and TrCel7A from Trichoderma reesei, acting on their natural, insoluble substrate, cellulose. Based on this, we consider practical compromises for the determination of experimental v_SS values, and propose a basic protocol that provides representative reaction rates and is experimentally simple so that larger groups of enzymes and conditions can be readily assayed with standard laboratory equipment. We surmise that the suggested experimental approach can be useful in comparative biochemical studies of cellulases; an area that remains poorly developed.

Systematic deletions in the cellobiohydrolase (CBH) Cel7A from the fungus Trichoderma reesei reveal flexible loops critical for CBH activity

Glycoside hydrolase family 7 (GH7) cellulases are some of the most efficient degraders of cellulose, making them particularly relevant for industries seeking to produce renewable fuels from lignocellulosic biomass. The secretome of the cellulolytic model fungus Trichoderma reesei contains two GH7s, termed TrCel7A and TrCel7B. Despite having high structural and sequence similarities, the two enzymes are functionally quite different. TrCel7A is an exolytic, processive cellobiohydrolase (CBH), with high activity on crystalline cellulose, whereas TrCel7B is an endoglucanase (EG) with a preference for more amorphous cellulose. At the structural level, these functional differences are usually ascribed to the flexible loops that cover the substrate-binding areas. TrCel7A has an extensive tunnel created by eight peripheral loops, and the absence of four of these loops in TrCel7B makes its catalytic domain a more open cleft. To investigate the structure-function relationships of these loops, here we produced and kinetically characterized several variants in which four loops unique to TrCel7A were individually deleted to resemble the arrangement in the TrCel7B structure. Analysis of a range of kinetic parameters consistently indicated that the B2 loop, covering the substrate-binding subsites -3 and -4 in TrCel7A, was a key determinant for the difference in CBH- or EG-like behavior between TrCel7A and TrCel7B. Conversely, the B3 and B4 loops, located closer to the catalytic site in TrCel7A, were less important for these activities. We surmise that these results could be useful both in further mechanistic investigations and for guiding engineering efforts of this industrially important enzyme family.

Thermoactivation of a cellobiohydrolase

We have measured activity and substrate affinity of the thermostable cellobiohydrolase, Cel7A, from Rasamsonia emersonii over a broad range of temperatures. For the wild type enzyme, which does not have a Carbohydrate Binding Module (CBM), higher temperature only led to moderately increased activity against cellulose, and we ascribed this to a pronounced, temperature induced desorption of enzyme from the substrate surface. We also tested a "high affinity" variant of R. emersonii Cel7A with a linker and CBM from a related enzyme. At room temperature, the activity of the variant was similar to the wild type, but the variant was more accelerated by temperature and about two-fold faster around 70 °C. This better thermoactivation of the high-affinity variant could not be linked to differences in stability or the catalytic process, but coincided with less desorption as temperature increased. Based on these observations and earlier reports on moderate thermoactivation of cellulases, we suggest that better cellulolytic activity at industrially relevant temperatures may be attained by engineering improved substrate affinity into enzymes that already possess good thermostability.

Correlation of structure, function and protein dynamics in GH7 cellobiohydrolases from Trichoderma atroviride, T. reesei and T. harzianum

BACKGROUND

The ascomycete fungus Trichoderma reesei is the predominant source of enzymes for industrial conversion of lignocellulose. Its glycoside hydrolase family 7 cellobiohydrolase (GH7 CBH) TreCel7A constitutes nearly half of the enzyme cocktail by weight and is the major workhorse in the cellulose hydrolysis process. The orthologs from Trichoderma atroviride (TatCel7A) and Trichoderma harzianum (ThaCel7A) show high sequence identity with TreCel7A, ~ 80%, and represent naturally evolved combinations of cellulose-binding tunnel-enclosing loop motifs, which have been suggested to influence intrinsic cellobiohydrolase properties, such as endo-initiation, processivity, and off-rate.

RESULTS

The TatCel7A, ThaCel7A, and TreCel7A enzymes were characterized for comparison of function. The catalytic domain of TatCel7A was crystallized, and two structures were determined: without ligand and with thio-cellotriose in the active site. Initial hydrolysis of bacterial cellulose was faster with TatCel7A than either ThaCel7A or TreCel7A. In synergistic saccharification of pretreated corn stover, both TatCel7A and ThaCel7A were more efficient than TreCel7A, although TatCel7A was more sensitive to thermal inactivation. Structural analyses and molecular dynamics (MD) simulations were performed to elucidate important structure/function correlations. Moreover, reverse conservation analysis (RCA) of sequence diversity revealed divergent regions of interest located outside the cellulose-binding tunnel of Trichoderma spp. GH7 CBHs.

CONCLUSIONS

We hypothesize that the combination of loop motifs is the main determinant for the observed differences in Cel7A activity on cellulosic substrates. Fine-tuning of the loop flexibility appears to be an important evolutionary target in Trichoderma spp., a conclusion supported by the RCA data. Our results indicate that, for industrial use, it would be beneficial to combine loop motifs from TatCel7A with the thermostability features of TreCel7A. Furthermore, one region implicated in thermal unfolding is suggested as a primary target for protein engineering.

The influence of different linker modifications on the catalytic activity and cellulose affinity of cellobiohydrolase Cel7A from Hypocrea jecorina

Various cellulases consist of a catalytic domain connected to a carbohydrate-binding module (CBM) by a flexible linker peptide. The linker if often strongly O-glycosylated and typically has a length of 20-50 amino acid residues. Functional roles, other than connecting the two folded domains, of the linker and its glycans, have been widely discussed, but experimental evidence remains sparse. One of the most studied cellulose degrading enzymes is the multi-domain cellobiohydrolase Cel7A from Hypocrea jecorina. Here, we designed variants of Cel7A with mutations in the linker region to elucidate the role of the linker. We found that moderate modification of the linker could result in significant changes in substrate affinity and catalytic efficacy. These changes were quite different for different linker variants. Thus, deletion of six residues near the catalytic domain had essentially no effects on enzyme function. Conversely, a substitution of four glycosylation sites near the middle of the linker reduced substrate affinity and increased maximal turnover. The observation of weaker binding provides some support of recent suggestions that linker glycans may be directly involved in substrate interactions. However, a variant with several inserted glycosylation sites near the CBM also showed lower affinity for the substrate compared to the wild-type, and we suggest that substrate interactions of the glycans depend on their exact location as well as other factors such as changes in structure and dynamics of the linker peptide.

Anomeric Selectivity and Product Profile of a Processive Cellulase

Cellobiohydrolases (CBHs) make up an important group of enzymes for both natural carbon cycling and industrial deconstruction of lignocellulosic biomass. The consecutive hydrolysis of one cellulose strand relies on an intricate pattern of enzyme-substrate interactions in the long, tunnel-shaped binding site of the CBH. In this work, we have investigated the initial complexation mode with cellulose of the most thoroughly studied CBH, Cel7A from Hypocrea jecorina (HjCel7A). We found that HjCel7A predominantly produces glucose when it initiates a processive run on insoluble microcrystalline cellulose, confirming the validity of an even and odd product ratio as an estimate of processivity. Moreover, the glucose released from cellulose was predominantly α-glucose. A link between the initial binding mode of the enzyme and the reducing end configuration was investigated by inhibition studies with the two anomers of cellobiose. A clear preference for β-cellobiose in product binding site +2 was observed for HjCel7A, but not the homologous endoglucanase, HjCe7B. Possible relationships between this anomeric preference in the product site and the prevalence of odd-numbered initial-cut products are discussed, and a correlation between processivity and anomer selectivity is proposed.