Search

Summary This perspective discusses a set of design principles for next-generation kinetically trapped, intrinsically circular polymers (iCPs) that are inherently, selectively, and expediently depolymerizable to their monomer state once their kinetic barriers of deconstruction are overcome, thereby enabling not only the ideal shortest chemical circularity but also tunable performance properties. After describing four elements of the design principles—thermodynamics and kinetics, strategies to overcome trade-offs and unify conflicting properties, predictive modeling, and supply-chain life-cycle assessment and techno-economic analysis, which are illustrated with state-of-the-art examples—it concludes with presenting key challenges and opportunities for sustainable development of iCPs.

Summary Reducing the carbon intensity of plastics production by sourcing sustainable feedstocks while simultaneously enabling effective polymer recycling represents a potential transformation of 21st century manufacturing. To evaluate technologies that could enable such changes, it is imperative to compare the sustainability of bio-based and/or circular plastic flows to those of incumbent manufacturing paradigms. To that end, we estimate the supply chain energy requirements and greenhouse gas (GHG) emissions associated with US-based plastics consumption. Major commodity polymers, each of which has a global consumption of at least 1 MMT per year, account for an estimated annual 3.2 quadrillion Btutus (quads) of energy and 104 MMTCO2e of GHG emissions in the US alone. This study serves as a foundation for comparing the supply chain energy requirements and GHG emissions of today’s plastics manufacturing to tomorrow’s disruptive technologies, to inform the development of bio-based plastics and the circular economy for synthetic polymers.

Valorization of all major lignocellulose components, including lignin, cellulose, and hemicellulose is critical for an economically viable bioeconomy. In most biochemical conversion approaches, the standard process separately upgrades sugar hydrolysates and lignin. Here, we present a new process concept based on an engineered microbe that could enable simultaneous upgrading of all lignocellulose streams, which has the ultimate potential to reduce capital cost and enable new metabolic engineering strategies. Pseudomonas putida is a robust microorganism capable of natively catabolizing aromatics, organic acids, and D-glucose. We engineered this strain to utilize D-xylose by tuning expression of a heterologous D-xylose transporter, catabolic genes xylAB, and pentose phosphate pathway (PPP) genes tal-tkt. We further engineered L-arabinose utilization via the PPP or an oxidative pathway. This resulted in a growth rate on xylose and arabinose of 0.32 h−1 and 0.38 h−1, respectively. Using the oxidative L-arabinose pathway with the PPP xylose pathway enabled D-glucose, D-xylose, and L-arabinose co-utilization in minimal medium using model compounds as well as real corn stover hydrolysate, with a maximum hydrolysate sugar consumption rate of 3.3 g/L/h. After modifying catabolite repression, our engineered P. putida simultaneously co-utilized five representative compounds from cellulose (D-glucose), hemicellulose (D-xylose, L-arabinose, and acetic acid), and lignin-related compounds (p-coumarate), demonstrating the feasibility of simultaneously upgrading total lignocellulosic biomass to value-added chemicals.

Microbial terephthalic acid (TPA) catabolic pathways are conserved among the few bacteria known to turnover this xenobiotic aromatic compound. However, to date there are few reported cases in which this pathway has been successfully expressed in heterologous hosts to impart efficient utilization of TPA as a sole carbon source. In this work, we aimed to engineer TPA conversion in Acinetobacter baylyi ADP1 via the heterologous expression of catabolic and transporter genes from a native TPA-utilizing bacterium. Specifically, we obtained ADP1-derived strains capable of growing on TPA as the sole carbon source using chromosomal insertion and targeted amplification of the tph catabolic operon from Comamonas sp. E6. Adaptive laboratory evolution was then used to improve growth on this substrate. TPA consumption rates of the evolved strains, which retained multiple copies of the tph genes, were ~0.2 g/L/h (or ~1 g TPA/g cells/h), similar to that of Comamonas sp. E6 and almost 2-fold higher than that of Rhodococcus jostii RHA1, another native TPA-utilizing strain. To evaluate TPA transport in the evolved ADP1 strains, we engineered a TPA biosensor consisting of the transcription factor TphR and a fluorescent reporter. In combination with whole-genome sequencing, the TPA biosensor revealed that transport of TPA was not mediated by the heterologous proteins from Comamonas sp. E6. Instead, the endogenous ADP1 muconate transporter MucK, a member of the major facilitator superfamily, was responsible for TPA transport in several evolved strains in which MucK variants were found to enhance TPA uptake. Furthermore, the IclR-type transcriptional regulator DcaS was identified as a repressor of mucK expression. Overall, this work presents an unexpected function of a native protein identified through gene amplification, adaptive laboratory evolution, and a combination of screening methods. This study also provides a TPA biosensor for application in ADP1 and identifies transporter variants for use in metabolic engineering applications focused on plastic upcycling of polyesters.

The sequence of the promoter used to drive expression of glcDEF operon in the engineered strain MFL185 in the original article was incorrect. As described here, we performed additional experiments that indicate expression of this operon was increased in MFL185, as intended. Ultimately, this error is immaterial with respect to the findings and conclusions reported in the original article.

Family 45 glycoside hydrolases (GH45) are endoglucanases that are integral to cellulolytic secretomes, and their ability to break down cellulose has been successfully exploited in textile and detergent industries. In addition to their industrial relevance, understanding the molecular mechanism of GH45-catalyzed hydrolysis is of fundamental importance because of their structural similarity to cell wall–modifying enzymes such as bacterial lytic transglycosylases (LTs) and expansins present in bacteria, plants, and fungi. Our understanding of the catalytic itinerary of GH45s has been incomplete because a crystal structure with substrate spanning the −1 to +1 subsites is currently lacking. Here we constructed and validated a putative Michaelis complex in silico and used it to elucidate the hydrolytic mechanism in a GH45, Cel45A from the fungus Humicola insolens, via unbiased simulation approaches. These molecular simulations revealed that the solvent-exposed active-site architecture results in lack of coordination for the hydroxymethyl group of the substrate at the −1 subsite. This lack of coordination imparted mobility to the hydroxymethyl group and enabled a crucial hydrogen bond with the catalytic acid during and after the reaction. This suggests the possibility of a nonhydrolytic reaction mechanism when the catalytic base aspartic acid is missing, as is the case in some LTs (murein transglycosylase A) and expansins. We calculated reaction free energies and demonstrate the thermodynamic feasibility of the hydrolytic and nonhydrolytic reaction mechanisms. Our results provide molecular insights into the hydrolysis mechanism in HiCel45A, with possible implications for elucidating the elusive catalytic mechanism in LTs and expansins.

Lignin, an abundant aromatic heteropolymer in secondary plant cell walls, is the single largest source of renewable aromatics in the biosphere. Leveraging this resource for renewable bioproducts through targeted microbial action depends on lignin fragment uptake by microbial hosts and subsequent enzymatic action to obtain the desired product. Recent computational work has emphasized that bacterial inner membranes are permeable to many aromatic compounds expected from lignin depolymerization processes. In this study, we expand on these findings through simulations for 42 lignin-related compounds across a gram-negative bacterial outer membrane model. Unbiased simulation trajectories indicate that spontaneous crossing for the full outer membrane is relatively rare at molecular simulation timescales, primarily due to preferential membrane partitioning and slow diffusion within the lipopolysaccharide layer within the outer membrane. Membrane partitioning and permeability coefficients were determined through replica exchange umbrella sampling simulations to overcome sampling limitations. We find that the glycosylated lipopolysaccharides found in the outer membrane increase the permeation barrier to many lignin-related compounds, particularly the most hydrophobic compounds. However, the effect is relatively modest; at industrially relevant concentrations, uncharged lignin-related compounds will readily diffuse across the outer membrane without the need for specific porins. Together, our results provide insight into the permeability of the bacterial outer membrane for assessing lignin fragment uptake and the future production of renewable bioproducts.

Summary Atomic layer deposition (ALD) improves the durability of metal catalysts using nanoscale metal oxide coatings. However, targeted coating strategies and economic models are lacking for process-specific deactivation challenges that account for implications at scale. Herein, we apply Al2O3 ALD to Pd/TiO2 to increase durability during hydrogenation of muconic acid, a bio-based platform chemical, to adipic acid. Initial coating development and characterization are performed on the milligram scale using stop-flow ALD. Subsequently, ALD coating scale is increased by 3 orders of magnitude using fluidized bed ALD. Activity, leaching resistance, and thermal stability are evaluated at each synthesis scale. ALD-coated catalysts retain up to 2-fold greater muconic acid hydrogenation activity and undergo significantly less physical restructuring than uncoated Pd/TiO2 after high-temperature treatments, while reducing Pd leaching by over 4-fold. Techno-economic analysis for an adipic acid biorefinery supports increased ALD material costs through catalyst lifetime extension, underscoring the potential viability of this technology.

Summary To drive innovation in chemical and material applications beyond what has been afforded by the mature petrochemical industry, new molecules that possess diverse chemical functionality are needed. One source of such molecules lies in the varied metabolic pathways that soil microbes utilize to catabolize aromatic compounds generated during plant decomposition. Here, we have engineered Pseudomonas putida KT2440 to convert these aromatic compounds to 15 catabolic intermediates that exhibit substantial chemical diversity. Bioreactor cultivations, analytical methods, and bench-scale separations were developed to enable production (up to 58 g/L), detection, and purification of each target molecule. We further engineered strains for production of a subset of these molecules from glucose, achieving a 41% molar yield of muconic acid. Finally, we produce materials from three compounds to illustrate the potential for realizing performance-advantaged properties relative to petroleum-derived analogs.

Summary Polyethylene terephthalate (PET) is the largest produced polyester globally with an annual production exceeding 26 million tons for use in carpet, clothing, and single-use beverage bottles, among others. Today, only PET bottles are reclaimed for recycling, albeit at a low reclamation rate, with most of the waste PET accumulating in landfills or the environment. In this study, PET is upcycled to higher-value, long-lifetime materials, namely two types of fiber-reinforced plastics (FRPs), via combination with renewably sourceable monomers. By harnessing the embodied energy in reclaimed PET (rPET) and implementing renewably sourceable monomers with distinct chemical functionality relative to petroleum building blocks, the resultant rPET-FRPs exhibit considerably improved material properties and are predicted to save 57% in the total supply chain energy and reduce greenhouse gas emissions by 40% over standard petroleum-based FRPs. Overall, this study enables a route to PET upcycling via bio-based monomers that could incentivize both improved plastics reclamation and acceleration of the bioeconomy.

Robust fluorescence-based biosensors are emerging as critical tools for high-throughput strain improvement in synthetic biology. Many biosensors are developed in model organisms where sophisticated synthetic biology tools are also well established. However, industrial biochemical production often employs microbes with phenotypes that are advantageous for a target process, and biosensors may fail to directly transition outside the host in which they are developed. In particular, losses in sensitivity and dynamic range of sensing often occur, limiting the application of a biosensor across hosts. Here we demonstrate the optimization of an Escherichia coli-based biosensor in a robust microbial strain for the catabolism of aromatic compounds, Pseudomonas putida KT2440, through a generalizable approach of modulating interactions at the protein-DNA interface in the promoter and the protein-protein dimer interface. The high-throughput biosensor optimization approach demonstrated here is readily applicable towards other allosteric regulators. Keywords: Whole cell biosensor, Aromatic catabolism, Transcription factor, PcaU, Shikimate

Summary Reductive catalytic fractionation (RCF) has emerged as a leading biomass fractionation and lignin valorization strategy. Here, flowthrough reactors were used to investigate RCF of poplar. Most RCF studies to date have been conducted in batch, but a flow-based process enables the acquisition of intrinsic kinetic and mechanistic data essential to accelerate the design, optimization, and scale-up of RCF processes. Time-resolved product distributions and yields obtained from experiments with different catalyst loadings were used to identify and deconvolute events during solvolysis and hydrogenolysis. Multi-bed RCF experiments provided unique insights into catalyst deactivation, showing that leaching, sintering, and surface poisoning are causes for decreased catalyst performance. The onset of catalyst deactivation resulted in higher concentrations of unsaturated lignin intermediates and increased occurrence of repolymerization reactions, producing high-molecular-weight species. Overall, this study demonstrates the concept of flowthrough RCF, which will be vital for realistic scale-up of this promising approach.

Family 7 glycoside hydrolases (GH7) are among the principal enzymes for cellulose degradation in nature and industrially. These important enzymes are often bimodular, comprised of a catalytic domain attached to a carbohydrate binding module (CBM) via a flexible linker, and exhibit a long active site that binds cello-oligomers of up to ten glucosyl moieties. GH7 cellulases consist of two major subtypes: cellobiohydrolases (CBH) and endoglucanases (EG). Despite the critical biological and industrial importance of GH7 enzymes, there remain gaps in our understanding of how GH7 sequence and structure relate to function. Here, we employed machine learning to gain insights into relationships between sequence, structure, and function across the GH7 family. Machine-learning models, using the number of residues in the active-site loops as features, were able discriminate GH7 CBHs and EGs with up to 99% accuracy. The lengths of the A4, B2, B3, and B4 loops were strongly correlated with functional subtype across the GH7 family. Position-specific classification rules were derived such that specific amino acids at 42 different sequence positions predicted the functional subtype with accuracies greater than 87%. A random forest model trained on residues at 19 positions in the catalytic domain predicted the presence of a CBM with 89.5% accuracy. We propose these positions play vital roles in the functional variation of GH7 cellulases. Taken together, our results complement numerous experimental findings and present functional relationships that can be applied when prospecting GH7 cellulases from nature, for sequence annotation, and to understand or manipulate function.

The sugar dehydration products, furfural and 5-(hydroxymethyl)furfural (HMF), are commonly formed during high-temperature processing of lignocellulose, most often in thermochemical pretreatment, liquefaction, or pyrolysis. Typically, these two aldehydes are considered major inhibitors in microbial conversion processes. Many microbes can convert these compounds to their less toxic, dead-end alcohol counterparts, furfuryl alcohol and 5-(hydroxymethyl)furfuryl alcohol. Recently, the genes responsible for aerobic catabolism of furfural and HMF were discovered in Cupriavidus basilensis HMF14 to enable complete conversion of these compounds to the TCA cycle intermediate, 2-oxo-glutarate. In this work, we engineer the robust soil microbe, Pseudomonas putida KT2440, to utilize furfural and HMF as sole carbon and energy sources via complete genomic integration of the 12 kB hmf gene cluster previously reported from Burkholderia phytofirmans. The common intermediate, 2-furoic acid, is shown to be a bottleneck for both furfural and HMF metabolism. When cultured on biomass hydrolysate containing representative amounts of furfural and HMF from dilute-acid pretreatment, the engineered strain outperforms the wild type microbe in terms of reduced lag time and enhanced growth rates due to catabolism of furfural and HMF. Overall, this study demonstrates that an approach for biological conversion of furfural and HMF, relative to the typical production of dead-end alcohols, enables both enhanced carbon conversion and substantially improves tolerance to hydrolysate inhibitors. This approach should find general utility both in emerging aerobic processes for the production of fuels and chemicals from biomass-derived sugars and in the biological conversion of high-temperature biomass streams from liquefaction or pyrolysis where furfural and HMF are much more abundant than in biomass hydrolysates from pretreatment., Highlights • HMF and furfural are common microbial inhibitors in biomass conversion. • HMF and furfural gene cluster was isolated from Burkholderia phytofirmans.. • We heterologously express the HMF/furfural gene cluster in Pseudomonas putida.. • Expression enables cultivation on HMF and furfural as a sole carbon source. • Expression also enables enhanced conversion on lignocellulosic hydrolysate.

Isopropanol is an important target molecule for sustainable production of fuels and chemicals. Increases in DNA synthesis and synthetic biology capabilities have resulted in the development of a range of new strategies for the more rapid design, construction, and testing of production strains. Here, we report on the use of such capabilities to construct and test 903 different variants of the isopropanol production pathway in Escherichia coli. We first constructed variants to explore the effect of codon optimization, copy number, and translation initiation rates on isopropanol production. The best strain (PA06) produced isopropanol at titers of 17.5 g/L, with a yield of 0.62 (mol/mol), and maximum productivity of 0.40 g/L/h. We next integrated the isopropanol synthetic pathway into the genome and then used the CRISPR EnAbled Trackable genome Engineering (CREATE) strategy to generate an additional 640 individual RBS library variants for further evaluation. After testing each of these variants, we constructed a combinatorial library containing 256 total variants from four different RBS levels for each gene. The best producing variant, PA14, produced isopropanol at titers of 7.1 g/L at 24 h, with a yield of 0.75 (mol/mol), and maximum productivity of 0.62 g/L/h (which was 0.22 g/L/h above the parent strain PA07). We demonstrate the ability to rapidly construct and test close to ~1000 designer strains and identify superior performers.

The production of chemicals alongside fuels will be essential to enhance the feasibility of lignocellulosic biorefineries. Succinic acid (SA), a naturally occurring C4-diacid, is a primary intermediate of the tricarboxylic acid cycle and a promising building block chemical that has received significant industrial attention. Basfia succiniciproducens is a relatively unexplored SA-producing bacterium with advantageous features such as broad substrate utilization, genetic tractability, and facultative anaerobic metabolism. Here B. succiniciproducens is evaluated in high xylose-content hydrolysates from corn stover and different synthetic media in batch fermentation. SA titers in hydrolysate at an initial sugar concentration of 60 g/L reached up to 30 g/L, with metabolic yields of 0.69 g/g, and an overall productivity of 0.43 g/L/h. These results demonstrate that B. succiniciproducens may be an attractive platform organism for bio-SA production from biomass hydrolysates.

Lucas D. Ellis received his BS from California Polytechnic State University, MS from Dartmouth College, and PhD focusing on heterogeneous catalysis from the University of Colorado. He is currently a Director’s Postdoctoral Fellow at the National Renewable Energy Laboratory (NREL). His research focuses on developing heterogeneous catalysts and chemical processes capable of reducing the impact of climate change. Gregg T. Beckham is a senior research fellow at NREL. He received a PhD and MSCEP from the Massachusetts Institute of Technology and a BS from Oklahoma State University. He currently leads an interdisciplinary team on the conversion of biomass to chemicals and materials and in the area of plastics upcycling.

Lignin represents an untapped feedstock for the production of fuels and chemicals, but its intrinsic heterogeneity makes lignin valorization a significant challenge. In nature, many aerobic organisms degrade lignin-derived aromatic molecules through conserved central intermediates including catechol and protocatechuate. Harnessing this microbial approach offers potential for lignin upgrading in modern biorefineries, but significant technical development is needed to achieve this end. Catechol and protocatechuate are subjected to aromatic ring cleavage by dioxygenase enzymes that, depending on the position, ortho or meta relative to adjacent hydroxyl groups, result in different products that are metabolized through parallel pathways for entry into the TCA cycle. These degradation pathways differ in the combination of succinate, acetyl-CoA, and pyruvate produced, the reducing equivalents regenerated, and the amount of carbon emitted as CO2-factors that will ultimately impact the yield of the targeted product. As shown here, the ring-cleavage pathways can be interchanged with one another, and such substitutions have a predictable and substantial impact on product yield. We demonstrate that replacement of the catechol ortho degradation pathway endogenous to Pseudomonas putida KT2440 with an exogenous meta-cleavage pathway from P. putida mt-2 increases yields of pyruvate produced from aromatic molecules in engineered strains. Even more dramatically, replacing the endogenous protocatechuate ortho pathway with a meta-cleavage pathway from Sphingobium sp. SYK-6 results in a nearly five-fold increase in pyruvate production. We further demonstrate the aerobic conversion of pyruvate to l-lactate with a yield of 41.1 ± 2.6% (wt/wt). Overall, this study illustrates how aromatic degradation pathways can be tuned to optimize the yield of a desired product in biological lignin upgrading.

International audience; Lytic polysaccharide monooxygenases (LPMOs) are a recently discovered class of enzymes that employ a copper-mediated, oxidative mechanism to cleave glycosidic bonds. The LPMO catalytic mechanism likely requires that molecular oxygen first binds to Cu(I), but the oxidation state in many reported LPMO structures is ambiguous, and the changes in the LPMO active site required to accommodate both oxidation states of copper have not been fully elucidated. Here, a diffraction data collection strategy minimizing the deposited x-ray dose was used to solve the crystal structure of a chitin-specific LPMO from Enterococcus faecalis (EfaCBM33A) in the Cu(II)-bound form. Subsequently, the crystalline protein was photoreduced in the x-ray beam, which revealed structural changes associated with the conversion from the initial Cu(II)-oxidized form with two coordinated water molecules, which adopts a trigonal bipyramidal geometry, to a reduced Cu(I) form in a T-shaped geometry with no coordinated water molecules. A comprehensive survey of Cu(II) and Cu(I) structures in the Cambridge Structural Database unambiguously shows that the geometries observed in the least and most reduced structures reflect binding of Cu(II) and Cu(I), respectively. Quantum mechanical calculations of the oxidized and reduced active sites reveal little change in the electronic structure of the active site measured by the active site partial charges. Together with a previous theoretical investigation of a fungal LPMO, this suggests significant functional plasticity in LPMO active sites. Overall, this study provides molecular snapshots along the reduction process to activate the LPMO catalytic machinery and provides a general method for solving LPMO structures in both copper oxidation states.

Enzymes acting on β-linked arabinofuranosides have been unknown until recently, in spite of wide distribution of β- l -arabinofuranosyl oligosaccharides in plant cells. Recently, a β- l -arabinofuranosidase from the glycoside hydrolase family 127 (HypBA1) was discovered in the newly characterized degradation system of hydroxyproline-linked β- l -arabinooligosaccharides in the bacterium Bifidobacterium longum . Here, we report the crystal structure of HypBA1 in the ligand-free and β- l -arabinofuranose complex forms. The structure of HypBA1 consists of a catalytic barrel domain and two additional β-sandwich domains, with one β-sandwich domain involved in the formation of a dimer. Interestingly, there is an unprecedented metal-binding motif with Zn 2+ coordinated by glutamate and three cysteines in the active site. The glutamate residue is located far from the anomeric carbon of the β- l -arabinofuranose ligand, but one cysteine residue is appropriately located for nucleophilic attack for glycosidic bond cleavage. The residues around the active site are highly conserved among GH127 members. Based on biochemical experiments and quantum mechanical calculations, a possible reaction mechanism involving cysteine as the nucleophile is proposed.

Carbohydrate structures are modified and degraded in the biosphere by a myriad of mostly hydrolytic enzymes. Recently, lytic polysaccharide mono-oxygenases (LPMOs) were discovered as a new class of enzymes for cleavage of recalcitrant polysaccharides that instead employ an oxidative mechanism. LPMOs employ copper as the catalytic metal and are dependent on oxygen and reducing agents for activity. LPMOs are found in many fungi and bacteria, but to date no basidiomycete LPMO has been structurally characterized. Here we present the three-dimensional crystal structure of the basidiomycete Phanerochaete chrysosporium GH61D LPMO, and, for the first time, measure the product distribution of LPMO action on a lignocellulosic substrate. The structure reveals a copper-bound active site common to LPMOs, a collection of aromatic and polar residues near the binding surface that may be responsible for regio-selectivity, and substantial differences in loop structures near the binding face compared with other LPMO structures. The activity assays indicate that this LPMO primarily produces aldonic acids. Last, molecular simulations reveal conformational changes, including the binding of several regions to the cellulose surface, leading to alignment of three tyrosine residues on the binding face of the enzyme with individual cellulose chains, similar to what has been observed for family 1 carbohydrate-binding modules. A calculated potential energy surface for surface translation indicates that P. chrysosporium GH61D exhibits energy wells whose spacing seems adapted to the spacing of cellobiose units along a cellulose chain.

Cellulase enzymes cleave glycosidic bonds in cellulose to produce cellobiose via either retaining or inverting hydrolysis mechanisms, which are significantly pH-dependent. Many fungal cellulases function optimally at pH ~5, and their activities decrease dramatically at higher or lower pH. To understand the molecular-level implications of pH in cellulase structure, we use a hybrid, solvent-based, constant pH molecular dynamics method combined with pH-based replica exchange to determine the pK(a) values of titratable residues of a glycoside hydrolase (GH) family 6 cellobiohydrolase (Cel6A) and a GH family 7 cellobiohydrolase (Cel7A) from the fungus Hypocrea jecorina. For both enzymes, we demonstrate that a bound substrate significantly affects the pKa values of the acid residues at the catalytic center. The calculated pK(a) values of catalytic residues confirm their proposed roles from structural studies and are consistent with the experimentally measured apparent pKa values. Additionally, GHs are known to impart a strained pucker conformation in carbohydrate substrates in active sites for catalysis, and results from free energy calculations combined with constant pH molecular dynamics suggest that the correct ring pucker is stable near the optimal pH for both Cel6A and Cel7A. Much longer molecular dynamics simulations of Cel6A and Cel7A with fixed protonation states based on the calculated pK(a) values suggest that pH affects the flexibility of tunnel loops, which likely affects processivity and substrate complexation. Taken together, this work demonstrates several molecular-level effects of pH on GH enzymes important for cellulose turnover in the biosphere and relevant to biomass conversion processes.

Cellulases and hemicellulases are responsible for the turnover of plant cell wall polysaccharides in the biosphere, and thus form the foundation of enzyme engineering efforts in biofuels research. Many of these carbohydrate-active enzymes from filamentous fungi contain both N-linked and O-linked glycosylation, the extent and heterogeneity of which depends on growth conditions, expression host, and the presence of glycan trimming enzymes in the secretome, all of which in turn impact enzyme activity. As the roles of glycosylation in enzyme function have not been fully elucidated, here we discuss the potential roles of glycosylation on glycoside hydrolase enzyme structure and function after secretion. We posit that glycosylation, instead of hindering cellulase engineering, can be used as an additional tool to enhance enzyme activity, given deeper understanding of its molecular-level role in biomass deconstruction.

Background: Family 1 carbohydrate-binding modules (CBMs) bind selectively to the hydrophobic surfaces of cellulose. Results: Simulations have shown that the planar face of the CBM binds preferentially to the hydrophobic face. Conclusion: Thermodynamic driving forces enable transfer of the CBM from the hydrophilic to hydrophobic surfaces. Significance: Selectivity of CBM provides access of cellulases to active surfaces of cellulose., Cellulase enzymes often contain carbohydrate-binding modules (CBMs) for binding to cellulose. The mechanisms by which CBMs recognize specific surfaces of cellulose and aid in deconstruction are essential to understand cellulase action. The Family 1 CBM from the Trichoderma reesei Family 7 cellobiohydrolase, Cel7A, is known to selectively bind to hydrophobic surfaces of native cellulose. It is most commonly suggested that three aromatic residues identify the planar binding face of this CBM, but several recent studies have challenged this hypothesis. Here, we use molecular simulation to study the CBM binding orientation and affinity on hydrophilic and hydrophobic cellulose surfaces. Roughly 43 μs of molecular dynamics simulations were conducted, which enables statistically significant observations. We quantify the fractions of the CBMs that detach from crystal surfaces or diffuse to other surfaces, the diffusivity along the hydrophobic surface, and the overall orientation of the CBM on both hydrophobic and hydrophilic faces. The simulations demonstrate that there is a thermodynamic driving force for the Cel7A CBM to bind preferentially to the hydrophobic surface of cellulose relative to hydrophilic surfaces. In addition, the simulations demonstrate that the CBM can diffuse from hydrophilic surfaces to the hydrophobic surface, whereas the reverse transition is not observed. Lastly, our simulations suggest that the flat faces of Family 1 CBMs are the preferred binding surfaces. These results enhance our understanding of how Family 1 CBMs interact with and recognize specific cellulose surfaces and provide insights into the initial events of cellulase adsorption and diffusion on cellulose.

Understanding the enzymatic mechanism that cellulases employ to degrade cellulose is critical to efforts to efficiently utilize plant biomass as a sustainable energy resource. A key component of cellulase action on cellulose is product inhibition from monosaccharide and disaccharides in the product site of cellulase tunnel. The absolute binding free energy of cellobiose and glucose to the product site of the catalytic tunnel of the Family 7 cellobiohydrolase (Cel7A) of Trichoderma reesei (Hypocrea jecorina) was calculated using two different approaches: steered molecular dynamics (SMD) simulations and alchemical free energy perturbation molecular dynamics (FEP/MD) simulations. For the SMD approach, three methods based on Jarzynski's equality were used to construct the potential of mean force from multiple pulling trajectories. The calculated binding free energies, -14.4 kcal/mol using SMD and -11.2 kcal/mol using FEP/MD, are in good qualitative agreement. Analysis of the SMD pulling trajectories suggests that several protein residues (Arg-251, Asp-259, Asp-262, Trp-376, and Tyr-381) play key roles in cellobiose and glucose binding to the catalytic tunnel. Five mutations (R251A, D259A, D262A, W376A, and Y381A) were made computationally to measure the changes in free energy during the product expulsion process. The absolute binding free energies of cellobiose to the catalytic tunnel of these five mutants are -13.1, -6.0, -11.5, -7.5, and -8.8 kcal/mol, respectively. The results demonstrated that all of the mutants tested can lower the binding free energy of cellobiose, which provides potential applications in engineering the enzyme to accelerate the product expulsion process and improve the efficiency of biomass conversion.

Fungi and bacteria secrete glycoprotein cocktails to deconstruct cellulose. Cellulose-degrading enzymes (cellulases) are often modular, with catalytic domains for cellulose hydrolysis and carbohydrate-binding modules connected by linkers rich in serine and threonine with O-glycosylation. Few studies have probed the role that the linker and O-glycans play in catalysis. Since different expression and growth conditions produce different glycosylation patterns that affect enzyme activity, the structure-function relationships that glycosylation imparts to linkers are relevant for understanding cellulase mechanisms. Here, the linker of the Trichoderma reesei Family 7 cellobiohydrolase (Cel7A) is examined by simulation. Our results suggest that the Cel7A linker is an intrinsically disordered protein with and without glycosylation. Contrary to the predominant view, the O-glycosylation does not change the stiffness of the linker, as measured by the relative fluctuations in the end-to-end distance; rather, it provides a 16 Å extension, thus expanding the operating range of Cel7A. We explain observations from previous biochemical experiments in the light of results obtained here, and compare the Cel7A linker with linkers from other cellulases with sequence-based tools to predict disorder. This preliminary screen indicates that linkers from Family 7 enzymes from other genera and other cellulases within T. reesei may not be as disordered, warranting further study.

Lignin is an abundant aromatic biopolymer within plant cell walls formed through radical coupling chemistry, whose composition and topology can vary greatly depending on the biomass source. Computational modeling provides a complementary approach to traditional experimental techniques to probe lignin interactions, lignin structure, and lignin material properties. However, current modeling approaches are limited based on the subset of lignin chemistries covered by existing lignin force fields. To fill the gap, we developed a comprehensive lignin force field that accounts for more lignin–lignin and lignin–carbohydrate interlinkages than existing lignin force fields, and also greatly expands the lignin monomer chemistries that can be modeled beyond simple alcohols and into the rich mixture of natural lignin varieties. The development of this force field utilizes recent developments in parameterization methodology, and synthesizes them into a workflow that combines target data from multiple molecules simultaneously into a single consistent and comprehensive parameter set. The parameter set represents a significant improvement to alternatives for atomic modeling of diverse lignin topologies, more accurately reproducing experimental observables while also significantly reducing the error relative to quantum calculations. The improved energetics, as well as the rigid adherence to CHARMM parameterization philosophy, enables simulation of lignin within its biological context with greater accuracy than was previously possible. The lignin force field presented here is therefore a crucial first step towards modeling lignin structure across a broad range of environments, including within plant cell walls where lignin is complexed with carbohydrates and deconstructed by bacterial or fungal enzymes, or as it exists within industrial solvent mixtures. Future simulations enabled by this updated lignin force field will thus lead to better chemical and structural understanding of lignin, providing new insight into its role in biomass recalcitrance or probing the potential for lignin to be used within industrial processes.