Free Energy Diagram for the Heterogeneous Enzymatic Hydrolysis of Glycosidic Bonds in Cellulose

Combining thermal model and kinetics: Implications in dynamic simulation of anaerobic digesters

In this study, ADM1-based kinetics were combined with a thermal model that accounts for various heat transfers inside and through the reactor's boundaries. Computing the energy of bioreactions based on kinetic rates prevented an overestimation of approximately 20% in the heat demand of the heat exchanger, compared to calculations using feedstock degradation heat. This also improved the representation of the digester's thermal and reactional inertia. Simulations across different climates demonstrated the necessity to tailor digester construction. In the North West United Kingdom, biogas auto-consumption was 23% higher than that of the same reactor in South West France. Enhancing the thermal insulation of the digester reduced heat losses by 37% in the United Kingdom. Therefore, coupling the kinetic and thermal models expands the insights that can be extracted from simulations, which can be valuable in optimizing the operation and design of biogas digesters.

Integrating electromagnetic cancer stress with immunotherapy: a therapeutic paradigm

An array of published cell-based and small animal studies have demonstrated a variety of exposures of cancer cells or experimental carcinomas to electromagnetic (EM) wave platforms that are non-ionizing and non-thermal. Overall effects appear to be inhibitory, inducing cancer cell stress or death as well as inhibition in tumor growth in experimental models. A variety of physical input variables, including discrete frequencies, amplitudes, and exposure times, have been tested, but drawing methodologic rationale and mechanistic conclusions across studies is challenging. Nevertheless, outputs such as tumor cytotoxicity, apoptosis, tumor membrane electroporation and leak, and reactive oxygen species generation are intriguing. Early EM platforms in humans employ pulsed electric fields applied either externally or using interventional tumor contact to induce tumor cell electroporation with stromal, vascular, and immunologic sparing. It is also possible that direct or external exposures to non-thermal EM waves or pulsed magnetic fields may generate electromotive forces to engage with unique tumor cell properties, including tumor glycocalyx to induce carcinoma membrane disruption and stress, providing novel avenues to augment tumor antigen release, cross-presentation by tumor-resident immune cells, and anti-tumor immunity. Integration with existing checkpoint inhibitor strategies to boost immunotherapeutic effects in carcinomas may also emerge as a broadly effective strategy, but little has been considered or tested in this area. Unlike the use of chemo/radiation and/or targeted therapies in cancer, EM platforms may allow for the survival of tumor-associated immunologic cells, including naïve and sensitized anti-tumor T cells. Moreover, EM-induced cancer cell stress and apoptosis may potentiate endogenous tumor antigen-specific anti-tumor immunity. Clinical studies examining a few of these combined EM-platform approaches are in their infancy, and a greater thrust in research (including basic, clinical, and translational work) in understanding how EM platforms may integrate with immunotherapy will be critical in driving advances in cancer outcomes under this promising combination.

Kinetics and thermodynamics of thermal inactivation for recombinant Escherichia coli cellulases, cel12B, cel8C, and polygalacturonase, peh28; biocatalysts for biofuel precursor production

Lignocellulosic biomass conversion using cellulases/polygalacturonases is a process that can be progressively influenced by several determinants involved in cellulose microfibril degradation. This article focuses on the kinetics and thermodynamics of thermal inactivation of recombinant Escherichia coli cellulases, cel12B, cel8C and a polygalacturonase, peh 28, derived from Pectobacterium carotovorum sub sp. carotovorum. Several consensus motifs conferring the enzymes' thermal stability in both cel12B and peh28 model structures have been detailed earlier, which were confirmed for the three enzymes through the current study of their thermal inactivation profiles over the 20-80°C range using the respective activities on carboxymethylcellulose and polygalacturonic acid. Kinetic constants and half-lives of thermal inactivation, inactivation energy, plus inactivation entropies, enthalpies and Gibbs free energies, revealed high stability, less conformational change and protein unfolding for cel12B and peh28 due to thermal denaturation compared to cel8C. The apparent thermal stability of peh28 and cel12B, along with their hydrolytic efficiency on a lignocellulosic biomass conversion as reported previously, makes these enzymes candidates for various industrial applications. Analysis of the Gibbs free energy values suggests that the thermal stabilities of cel12B and peh28 are entropy-controlled over the tested temperature range.

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.

Enzymatic processing of lignocellulosic biomass: principles, recent advances and perspectives

Efficient saccharification of lignocellulosic biomass requires concerted development of a pretreatment method, an enzyme cocktail and an enzymatic process, all of which are adapted to the feedstock. Recent years have shown great progress in most aspects of the overall process. In particular, increased insights into the contributions of a wide variety of cellulolytic and hemicellulolytic enzymes have improved the enzymatic processing step and brought down costs. Here, we review major pretreatment technologies and different enzyme process setups and present an in-depth discussion of the various enzyme types that are currently in use. We pay ample attention to the role of the recently discovered lytic polysaccharide monooxygenases (LPMOs), which have led to renewed interest in the role of redox enzyme systems in lignocellulose processing. Better understanding of the interplay between the various enzyme types, as they may occur in a commercial enzyme cocktail, is likely key to further process improvements.

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.

The dissociation mechanism of processive cellulases

Cellulase enzymes deconstruct recalcitrant cellulose into soluble sugars, making them a biocatalyst of biotechnological interest for use in the nascent lignocellulosic bioeconomy. Cellobiohydrolases (CBHs) are cellulases capable of liberating many sugar molecules in a processive manner without dissociating from the substrate. Within the complete processive cycle of CBHs, dissociation from the cellulose substrate is rate limiting, but the molecular mechanism of this step is unknown. Here, we present a direct comparison of potential molecular mechanisms for dissociation via Hamiltonian replica exchange molecular dynamics of the model fungal CBH, Trichoderma reesei Cel7A. Computational rate estimates indicate that stepwise cellulose dethreading from the binding tunnel is 4 orders of magnitude faster than a clamshell mechanism, in which the substrate-enclosing loops open and release the substrate without reversing. We also present the crystal structure of a disulfide variant that covalently links substrate-enclosing loops on either side of the substrate-binding tunnel, which constitutes a CBH that can only dissociate via stepwise dethreading. Biochemical measurements indicate that this variant has a dissociation rate constant essentially equivalent to the wild type, implying that dethreading is likely the predominant mechanism for dissociation.

Predominant Nonproductive Substrate Binding by Fungal Cellobiohydrolase I and Implications for Activity Improvement

Enzymatic conversion of the most abundant renewable source of organic compounds, cellulose to fermentable sugars is attractive for production of green fuels and chemicals. The major component of industrial enzyme systems, cellobiohydrolase I from Hypocrea jecorina (Trichoderma reesei) (HjCel7A) processively splits disaccharide units from the reducing ends of tightly packed cellulose chains. HjCel7A consists of a catalytic domain (CD) and a carbohydrate-binding module (CBM) separated by a linker peptide. A tunnel-shaped substrate-binding site in the CD includes nine subsites for β-d-glucose units, seven of which (-7 to -1) precede the catalytic center. Low catalytic activity of Cel7A is the bottleneck and the primary target for improvement. Here it is shown for the first time that, in spite of much lower apparent k_cat of HjCel7A at the hydrolysis of β-1,4-glucosidic linkages in the fluorogenic cellotetra- and -pentaose compared to the structurally related endoglucanase I (HjCel7B), the specificity constants (catalytic efficiency) k_cat /K_m for both enzymes are almost equal in these reactions. The observed activity difference appears from strong nonproductive substrate binding by HjCel7A, particularly significant for MU-β-cellotetraose (MUG₄ ). Interaction of substrates with the subsites -6 and -5 proximal to the nonconserved Gln101 residue in HjCel7A decreases K_m,ap by >1500 times. HjCel7A can be nonproductively bound onto cellulose surface with K_d ≈2-9 nM via CBM and CD that captures six terminal glucose units of cellulose chain. Decomposition of this nonproductive complex can determine the rate of cellulose conversion. MUG₄ is a promising substrate to select active cellobiohydrolase I variants with reduced nonproductive substrate binding.

Engineering enhanced cellobiohydrolase activity

Glycoside Hydrolase Family 7 cellobiohydrolases (GH7 CBHs) catalyze cellulose depolymerization in cellulolytic eukaryotes, making them key discovery and engineering targets. However, there remains a lack of robust structure–activity relationships for these industrially important cellulases. Here, we compare CBHs from Trichoderma reesei (TrCel7A) and Penicillium funiculosum (PfCel7A), which exhibit a multi-modular architecture consisting of catalytic domain (CD), carbohydrate-binding module, and linker. We show that PfCel7A exhibits 60% greater performance on biomass than TrCel7A. To understand the contribution of each domain to this improvement, we measure enzymatic activity for a library of CBH chimeras with swapped subdomains, demonstrating that the enhancement is mainly caused by PfCel7A CD. We solve the crystal structure of PfCel7A CD and use this information to create a second library of TrCel7A CD mutants, identifying a TrCel7A double mutant with near-equivalent activity to wild-type PfCel7A. Overall, these results reveal CBH regions that enable targeted activity improvements.

Cellobiohydrolases (CBHs) are critical for natural and industrial biomass degradation but their structure–activity relationships are not fully understood. Here, the authors present the biochemical and structural characterization of two CBHs, identifying protein regions that confer enhanced CBH activity.

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.

Cellular automata modeling depicts degradation of cellulosic material by a cellulase system with single-molecule resolution

BACKGROUND

Enzymatic hydrolysis of cellulose involves the spatiotemporally correlated action of distinct polysaccharide chain cleaving activities confined to the surface of an insoluble substrate. Because cellulases differ in preference for attacking crystalline compared to amorphous cellulose, the spatial distribution of structural order across the cellulose surface imposes additional constraints on the dynamic interplay between the enzymes. Reconstruction of total system behavior from single-molecule activity parameters is a longstanding key goal in the field.

RESULTS

We have developed a stochastic, cellular automata-based modeling approach to describe degradation of cellulosic material by a cellulase system at single-molecule resolution. Substrate morphology was modeled to represent the amorphous and crystalline phases as well as the different spatial orientations of the polysaccharide chains. The enzyme system model consisted of an internally chain-cleaving endoglucanase (EG) as well as two processively acting, reducing and non-reducing chain end-cleaving cellobiohydrolases (CBHs). Substrate preference (amorphous: EG, CBH II; crystalline: CBH I) and characteristic frequencies for chain cleavage, processive movement, and dissociation were assigned from biochemical data. Once adsorbed, enzymes were allowed to reach surface-exposed substrate sites through "random-walk" lateral diffusion or processive motion. Simulations revealed that slow dissociation of processive enzymes at obstacles obstructing further movement resulted in local jamming of the cellulases, with consequent delay in the degradation of the surface area affected. Exploiting validation against evidence from atomic force microscopy imaging as a unique opportunity opened up by the modeling approach, we show that spatiotemporal characteristics of cellulose surface degradation by the system of synergizing cellulases were reproduced quantitatively at the nanometer resolution of the experimental data. This in turn gave useful prediction of the soluble sugar release rate.

CONCLUSIONS

Salient dynamic features of cellulose surface degradation by different cellulases acting in synergy were reproduced in simulations in good agreement with evidence from high-resolution visualization experiments. Due to the single-molecule resolution of the modeling approach, the utility of the presented model lies not only in predicting system behavior but also in elucidating inherently complex (e.g., stochastic) phenomena involved in enzymatic cellulose degradation. Thus, it creates synergy with experiment to advance the mechanistic understanding for improved application.

Temperature Effects on Kinetic Parameters and Substrate Affinity of Cel7A Cellobiohydrolases

We measured hydrolytic rates of four purified cellulases in small increments of temperature (10-50 °C) and substrate loads (0-100 g/liter) and analyzed the data by a steady state kinetic model that accounts for the processive mechanism. We used wild type cellobiohydrolases (Cel7A) from mesophilic Hypocrea jecorina and thermophilic Rasamsonia emersonii and two variants of these enzymes designed to elucidate the role of the carbohydrate binding module (CBM). We consistently found that the maximal rate increased strongly with temperature, whereas the affinity for the insoluble substrate decreased, and as a result, the effect of temperature depended strongly on the substrate load. Thus, temperature had little or no effect on the hydrolytic rate in dilute substrate suspensions, whereas strong temperature activation (Q10 values up to 2.6) was observed at saturating substrate loads. The CBM had a dual effect on the activity. On one hand, it diminished the tendency of heat-induced desorption, but on the other hand, it had a pronounced negative effect on the maximal rate, which was 2-fold larger in variants without CBM throughout the investigated temperature range. We conclude that although the CBM is beneficial for affinity it slows down the catalytic process. Cel7A from the thermophilic organism was moderately more activated by temperature than the mesophilic analog. This is in accord with general theories on enzyme temperature adaptation and possibly relevant information for the selection of technical cellulases.