Search

Lignin is a major component of plant cell walls that is typically underutilized in selective conversion strategies for renewable fuels and chemicals. The mechanisms by which thermal and catalytic treatments deconstruct lignin remain elusive, which is where quantum mechanical calculations can offer fundamental insights. Here, we compute homolytic bond dissociation enthalpies (BDEs) for four prevalent linkages in 69 lignin model compounds, including β-O-4, α-O-4, β-5, and biphenyl bonds, with a large range of natural and oxidized substituents. These calculations include ab initio benchmark values extrapolated to the complete basis set limit and full conformational searches for each compound. The results quantify both the relative BDEs among common lignin bonds and the effect of native and oxidized substituents on the functional groups in lignin. These data yield insights into thermal lignin deconstruction for a large range of prevalent linkages and aid in the identification of targets for catalytic cleavage. © 2011 American Chemical Society.

Methane hydrate nucleation rates are reported from over 200 μs of molecular dynamics simulations across a range of thermodynamic conditions and varying degrees of methane–water interfacial curvatur...

Plants represent a vast, renewable resource and are well suited to provide sustainably for humankind's transportation fuel needs. To produce infrastructure-compatible fuels from biomass, two challenges remain: overcoming plant cell wall recalcitrance to extract sugar and phenolic intermediates, and reduction of oxygenated intermediates to fuel molecules. To compete with fossil-based fuels, two primary routes to deconstruct cell walls are under development, namely biochemical and thermochemical conversion. Here, we focus on overcoming recalcitrance with biochemical conversion, which uses low-severity thermochemical pretreatment followed by enzymatic hydrolysis to produce soluble sugars. Many challenges remain, including understanding how pretreatments affect the physicochemical nature of heterogeneous cell walls; determination of how enzymes deconstruct the cell wall effectively with the aim of designing superior catalysts; and resolution of issues associated with the co-optimization of pretreatment, enzymatic hydrolysis, and fermentation. Here, we highlight some of the scientific challenges and open questions with a particular focus on problems across multiple length scales.

Bundles of cellulose polymers in plant cell walls exhibit a robust network of hydrogen bonds and hydrophobic interactions that must be overcome to decrystallize and hydrolyze the individual polymers to sugars. To investigate the molecular interactions that impart recalcitrance and insolubility to cellulose, we use simulation to determine the decrystallization work, a fundamental, yet experimentally inaccessible measurement, of cello-oligomers from the middle and edge of the hydrophobic face of cellulose Iβ. We demonstrate that cellobiose and cellotetraose decrystallization work does not depend on the position of the oligomer on the crystal surface but that larger oligomers are more difficult to decrystallize depending on the number of intralayer neighbors due to a larger number of stabilizing intralayer hydrogen bonds. The presented results are relevant to mesoscale, morphology-based models of cellulose deconstruction, understanding the molecular details of cellulose decrystallization and insolubility, an...

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.

We determine the mechanism for the liquid-solid phase transition in the Lennard-Jones fluid close to coexistence with aimless shooting and likelihood maximization. The reaction coordinate for this process is a product of a structural descriptor and the size of the nascent solid nucleus and is quantitatively verified with the committor probability histogram test. This study identifies the first accurate scalar reaction coordinate for the liquid-solid nucleation process in Lennard-Jonesium, which will likely extend to nucleation processes in other spherically symmetric fluids. On the basis of our results, we propose a structural correction factor for the commonly cited nucleus size reaction coordinate from classical nucleation theory that enables connection of simulation data to stochastic models of nucleation kinetics. In addition, we show that aimless shooting is able to obtain reasonable acceptance rates for transitions with highly diffusive characteristics, which has been problematic for transition path sampling methods for diffusive processes such as nucleation and macromolecular transitions.

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.

A multiscale simulation model is used to construct potential and free energy surfaces for the carbohydrate-binding module [CBM] from an industrially important cellulase, Trichoderma reesei cellobiohydrolase I, on the hydrophobic face of a coarse-grained cellulose Ibeta polymorph. We predict from computation that the CBM alone exhibits regions of stability on the hydrophobic face of cellulose every 5 and 10 A, corresponding to a glucose unit and a cellobiose unit, respectively. In addition, we predict a new role for the CBM: specifically, that in the presence of hydrolyzed cellulose chain ends, the CBM exerts a thermodynamic driving force to translate away from the free cellulose chain ends. This suggests that the CBM is not only required for binding to cellulose, as has been known for two decades, but also that it has evolved to both assist the enzyme in recognizing a cellulose chain end and exert a driving force on the enzyme during processive hydrolysis of cellulose.

We will synthesize and screen solid catalysts for the depolymerization of lignin to monomeric and oligomeric oxygenated species, which could be fractionated and integrated into refinery intermediate streams for selective upgrading, or catalytically upgraded to fuels and chemicals. This work will primarily focus on the synthesis and application of layered double hydroxides (LDHs) as recyclable, heterogeneous catalysts for depolymerization of lignin model compounds and softwood lignin. LDHs have been shown in our group to offer good supports and catalysts to promote base-catalyzed depolymerization of lignin model compounds and in preliminary experiments for the depolymerization of lignin from an Organosolv process. We will also include additional catalyst supports such as silica, alumina, and carbon as identified in ongoing and past efforts at NREL. This work will consist of two tasks. Overall, this work will be synergistic with ongoing efforts at NREL, funded by the DOE Biomass Program, on the development of catalysts for lignin depolymerization in the context of biochemical and thermochemical conversion of corn stover and other biomass feedstocks to advanced fuels and chemicals.

UNLABELLED Family 1 carbohydrate-binding modules (CBMs) are ubiquitous components of multimodular fungal enzymes that degrade plant cell wall polysaccharides and bind specifically to cellulose. Native glycosylation of family 1 CBMs has been shown to substantially impact multiple physical properties, including thermal and proteolytic stability and cellulose binding affinity. To gain molecular insights into the changes in CBM properties upon glycosylation, solution structures of two glycoforms of a Trichoderma reesei family 1 CBM were studied by NMR spectroscopy: a glycosylated family 1 CBM with a mannose group attached to both Thr1 and Ser3 and a second family 1 CBM with single mannose groups attached to Thr1, Ser3 and Ser14. The structures clearly reveal that monosaccharides at both Ser3 and Ser14 on family 1 CBMs present additional cellulose binding platforms, similar to well-characterized aromatic residues at the binding interface, which align to the cellulose surface. These results are in agreement with previous experimental work demonstrating that glycans at Ser3 and Ser14 impart significant improvements in binding affinity. Additionally, detailed analysis of the NMR structures and molecular simulations indicates that the protein backbone of the CBM is not significantly altered by attachment of monosaccharides, and that the mannose attached to Ser14 may be more flexible than the mannose at Ser3. Overall, the present study reveals how family 1 CBM structures are affected by covalent attachment of monosaccharides, which are likely important post-translational modifications of these common subdomains of fungal plant cell wall degrading enzymes. DATABASE Structural data have been deposited in the RCSB Protein Data Bank (PDB codes: 2MWJ and 2MWK) and the BioMagRes Bank (BMRB codes: 25331 and 25332) for CBM_M2 and CBM_M3, respectively.

The importance of the glycan structure and size, amino acid residues near the glycosylation site, and glycosidic linkage in controlling the effects of CBM O-glycosylation is shown., Protein glycosylation is a ubiquitous post-translational modification in all kingdoms of life. Despite its importance in molecular and cellular biology, the molecular-level ramifications of O-glycosylation on biomolecular structure and function remain elusive. Here, we took a small model glycoprotein and changed the glycan structure and size, amino acid residues near the glycosylation site, and glycosidic linkage while monitoring any corresponding changes to physical stability and cellulose binding affinity. The results of this study reveal the collective importance of all the studied features in controlling the most pronounced effects of O-glycosylation in this system. Going forward, this study suggests the possibility of designing proteins with multiple improved properties by simultaneously varying the structures of O-glycans and amino acids local to the glycosylation site.

Clathrate hydrates are solid crystalline structures most commonly formed from solutions that have nucleated to form a mixed solid composed of water and gas. Understanding the mechanism of clathrate hydrate nucleation is essential to grasp the fundamental chemistry of these complex structures and their applications. Molecular dynamics (MD) simulation is an ideal method to study nucleation at the molecular level because the size of the critical nucleus and formation rate occur on the nano scale. Various analysis methods for nucleation have been developed through MD to analyze nucleation. In particular, the mean first-passage time (MFPT) and survival probability (SP) methods have proven to be effective in procuring the nucleation rate and critical nucleus size for monatomic systems. This study assesses the MFPT and SP methods, previously used for monatomic systems, when applied to analyzing clathrate hydrate nucleation. Because clathrate hydrate nucleation is relatively difficult to observe in MD simulations (due to its high free energy barrier), these methods have yet to be applied to clathrate hydrate systems. In this study, we have analyzed the nucleation rate and critical nucleus size of methane hydrate using MFPT and SP methods from data generated by MD simulations at 255 K and 50 MPa. MFPT was modified for clathrate hydrate from the original version by adding the maximum likelihood estimate and growth effect term. The nucleation rates calculated by MFPT and SP methods are within 5%, and the critical nucleus size estimated by the MFPT method was 50% higher, than values obtained through other more rigorous but computationally expensive estimates. These methods can also be extended to the analysis of other clathrate hydrates.

Dehydration over acidic zeolites is an important reaction class for the upgrading of biomass pyrolysis vapors to hydrocarbon fuels or to precursors for myriad chemical products. Here, we examine the dehydration of ethanol at a Brønsted acid site, T12, found in HZSM-5 using density functional theory (DFT). The geometries of both cluster and mixed quantum mechanics/molecular mechanics (QM:MM) models are prepared from the ZSM-5 crystal structure. Comparisons between these models and different DFT methods are conducted to show similar results among the models and methods used. Inclusion of the full catalyst cavity through a QM:MM approach is found to be important, since activation barriers are computed on average as 7 kcal mol(-1) lower than those obtained with a smaller cluster model. Two different pathways, concerted and stepwise, have been considered when examining dehydration and deprotonation steps. The current study shows that a concerted dehydration process is possible with a lower (4-5 kcal mol(-1)) activation barrier while previous literature studies have focused on a stepwise mechanism. Overall, this work demonstrates that fairly high activation energies (∼50 kcal mol(-1)) are required for ethanol dehydration. A concerted mechanism is favored over a stepwise mechanism because charge separation in the transition state is minimized. QM:MM approaches appear to provide superior results to cluster calculations due to a more accurate representation of charges on framework oxygen atoms.

In nature, polysaccharide glycosidic bonds are cleaved by hydrolytic enzymes for a vast array of biological functions. Recently, a new class of enzymes that utilize an oxidative mechanism to cleave glycosidic linkages was discovered; these enzymes are called lytic polysaccharide monooxygenases (LPMO). These oxidative enzymes are synergistic with cocktails of hydrolytic enzymes and are thought to act primarily on crystalline regions, in turn providing new sites of productive attachment and detachment for processive hydrolytic enzymes. In the case of cellulose, the homopolymer of β-1,4-d-glucose, enzymatic oxidation occurs at either the reducing end or the nonreducing end of glucose, depending on enzymatic specificity, and results in the generation of oxidized chemical substituents at polymer chain ends. LPMO oxidation of cellulose is thought to produce either a lactone at the reducing end of glucose that can spontaneously or enzymatically convert to aldonic acid or 4-keto-aldose at the nonreducing end that may further oxidize to a geminal diol. Here, we use molecular simulation to examine the effect of oxidation on the structure of crystalline cellulose. The simulations highlight variations in behaviors depending on the chemical identity of the oxidized species and its location within the cellulose fibril, as different oxidized species introduce steric effects that disrupt local crystallinity and in some cases reduce the work needed for polymer decrystallization. Reducing-end oxidations are easiest to decrystallize when located at the end of the fibril, whereas nonreducing end oxidations readily decrystallize from internal cleavage sites despite their lower solvent accessibility. The differential in decrystallization free energy suggests a molecular mechanism consistent with experimentally observed LPMO/cellobiohydrolase synergy. Additionally, the soluble oxidized cellobiose products released by hydrolytic cellulases may bind to the active sites of cellulases with different affinities relative to cellobiose itself, which potentially affects hydrolytic turnover through product inhibition. To examine the effect of oxidation on cello-oligomer binding, we use thermodynamic integration to compute the relative change in binding free energy between the hydrolyzed and oxidized products in the active site of Family 7 and Family 6 processive glycoside hydrolases, Trichoderma reesei Cel7A and Cel6A, which are key industrial cellulases and commonly used model systems for fungal cellulases. Our results suggest that the equilibrium between the two reducing end oxidized products, favoring the linear aldonic acid, may increase product inhibition, which would in turn reduce processive substrate turnover. In the case of LMPO action at the nonreducing end, oxidation appears to lower affinity with the nonreducing end specific cellulase, reducing product inhibition and potentially promoting processive cellulose turnover. Overall, this suggests that oxidation of recalcitrant polysaccharides by LPMOs accelerates degradation not only by increasing the concentration of chain termini but also by reducing decrystallization work, and that product inhibition may be somewhat reduced as a result.

Cellulose in plant cell walls is the largest reservoir of renewable carbon on Earth. The saccharification of cellulose from plant biomass into soluble sugars can be achieved using fungal and bacterial cellulolytic enzymes, cellulases, and further converted into fuels and chemicals. Most fungal cellulases are both N- and O-glycosylated in their native form, yet the consequences of glycosylation on activity and structure are not fully understood. Studying protein glycosylation is challenging as glycans are extremely heterogeneous, stereochemically complex, and glycosylation is not under direct genetic control. Despite these limitations, many studies have begun to unveil the role of cellulase glycosylation, especially in the industrially relevant cellobiohydrolase from Trichoderma reesei, Cel7A. Glycosylation confers many beneficial properties to cellulases including enhanced activity, thermal and proteolytic stability, and structural stabilization. However, glycosylation must be controlled carefully as such positive effects can be dampened or reversed. Encouragingly, methods for the manipulation of glycan structures have been recently reported that employ genetic tuning of glycan-active enzymes expressed from homogeneous and heterologous fungal hosts. Taken together, these studies have enabled new strategies for the exploitation of protein glycosylation for the production of enhanced cellulases for biofuel production.

The goal for catalyst development, both homogeneous and heterogeneous, using bio-based feedstocks is designing new and robust materials for the dehydration, decarboxylation, decarbonylation, hydrogenolysis, esterification, and ketonization reactions that are required for converting these renewable, generally oxygenated feedstocks to desirable and renewable products. This chapter summarizes the work done on developing bio-based products in the last decade using the seminal US Department of Energy report on the top twelve bio-based chemicals as a starting point to assess catalyst improvement, novel process options, commercialization potential, and market penetration when possible.

The complete mitochondrial genome of Limnoria quadripunctata, a marine wood-eating isopod crustacean, was determined from whole genome sequence data. The mitogenome is 16,503 bp in length and contains 39 genes: 13 protein-coding, 2 ribosomal RNA, 22 tRNA, two of which are repeated and a control region. The start codon most commonly used by the Limnoria protein-coding genes is ATN, as is the case in the two other available complete isopod mitogenomes. The gene arrangement differs among these complete isopod mitogenomes, as does the AT-content of H-strand protein-coding genes. The latter observations, coupled with the considerable nucleotide diversity observed between the isopod mitogenomes, support the idea that each isopod species belongs to a distinct lineage as implied by their current placement in separate suborders.

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.

Translocation of carbohydrate polymers through protein tunnels and clefts is a ubiquitous biochemical phenomenon in proteins such as polysaccharide synthases, glycoside hydrolases, and carbohydrate-binding modules. Although static snapshots of carbohydrate polymer binding in proteins have long been studied via crystallography and spectroscopy, the molecular details of polysaccharide chain processivity have not been elucidated. Here, we employ simulation to examine how a cellulose chain translocates by a disaccharide unit during the processive cycle of a glycoside hydrolase family 7 cellobiohydrolase. Our results demonstrate that these biologically and industrially important enzymes employ a two-step mechanism for chain threading to form a Michaelis complex and that the free energy barrier to chain threading is significantly lower than the hydrolysis barrier. Taken with previous studies, our findings suggest that the rate-limiting step in enzymatic cellulose degradation is the glycosylation reaction, not chain processivity. Based on the simulations, we find that strong electrostatic interactions with polar residues that are conserved in GH7 cellobiohydrolases, but not in GH7 endoglucanases, at the leading glucosyl ring provide the thermodynamic driving force for polysaccharide chain translocation. Also, we consider the role of aromatic-carbohydrate interactions, which are widespread in carbohydrate-active enzymes and have long been associated with processivity. Our analysis suggests that the primary role for these aromatic residues is to provide tunnel shape and guide the carbohydrate chain to the active site. More broadly, this work elucidates the role of common protein motifs found in carbohydrate-active enzymes that synthesize or depolymerize polysaccharides by chain translocation mechanisms coupled to catalysis.

Background Lignin, nature’s dominant aromatic polymer, is found in most terrestrial plants in the approximate range of 15 to 40% dry weight and provides structural integrity. Traditionally, most large-scale industrial processes that use plant polysaccharides have burned lignin to generate the power needed to productively transform biomass. The advent of biorefineries that convert cellulosic biomass into liquid transportation fuels will generate substantially more lignin than necessary to power the operation, and therefore efforts are underway to transform it to value-added products. Production of biofuels from cellulosic biomass requires separation of large quantities of the aromatic polymer lignin. In planta genetic engineering, enhanced extraction methods, and a deeper understanding of the structure of lignin are yielding promising opportunities for efficient conversion of this renewable resource to carbon fibers, polymers, commodity chemicals, and fuels. [Credit: Oak Ridge National Laboratory, U.S. Department of Energy] Advances Bioengineering to modify lignin structure and/or incorporate atypical components has shown promise toward facilitating recovery and chemical transformation of lignin under biorefinery conditions. The flexibility in lignin monomer composition has proven useful for enhancing extraction efficiency. Both the mining of genetic variants in native populations of bioenergy crops and direct genetic manipulation of biosynthesis pathways have produced lignin feedstocks with unique properties for coproduct development. Advances in analytical chemistry and computational modeling detail the structure of the modified lignin and direct bioengineering strategies for targeted properties. Refinement of biomass pretreatment technologies has further facilitated lignin recovery and enables catalytic modifications for desired chemical and physical properties. Outlook Potential high-value products from isolated lignin include low-cost carbon fiber, engineering plastics and thermoplastic elastomers, polymeric foams and membranes, and a variety of fuels and chemicals all currently sourced from petroleum. These lignin coproducts must be low cost and perform as well as petroleum-derived counterparts. Each product stream has its own distinct challenges. Development of renewable lignin-based polymers requires improved processing technologies coupled to tailored bioenergy crops incorporating lignin with the desired chemical and physical properties. For fuels and chemicals, multiple strategies have emerged for lignin depolymerization and upgrading, including thermochemical treatments and homogeneous and heterogeneous catalysis. The multifunctional nature of lignin has historically yielded multiple product streams, which require extensive separation and purification procedures, but engineering plant feedstocks for greater structural homogeneity and tailored functionality reduces this challenge.

The majority of biological turnover of lignocellulosic biomass in nature is conducted by fungi, which commonly use Family 1 carbohydrate-binding modules (CBMs) for targeting enzymes to cellulose. Family 1 CBMs are glycosylated, but the effects of glycosylation on CBM function remain unknown. Here, the effects of O-mannosylation are examined on the Family 1 CBM from the Trichoderma reesei Family 7 cellobiohydrolase at three glycosylation sites. To enable this work, a procedure to synthesize glycosylated Family 1 CBMs was developed. Subsequently, a library of 20 CBMs was synthesized with mono-, di-, or trisaccharides at each site for comparison of binding affinity, proteolytic stability, and thermostability. The results show that, although CBM mannosylation does not induce major conformational changes, it can increase the thermolysin cleavage resistance up to 50-fold depending on the number of mannose units on the CBM and the attachment site. O-Mannosylation also increases the thermostability of CBM glycoforms up to 16 °C, and a mannose disaccharide at Ser3 seems to have the largest themostabilizing effect. Interestingly, the glycoforms with small glycans at each site displayed higher binding affinities for crystalline cellulose, and the glycoform with a single mannose at each of three positions conferred the highest affinity enhancement of 7.4-fold. Overall, by combining chemical glycoprotein synthesis and functional studies, we show that specific glycosylation events confer multiple beneficial properties on Family 1 CBMs.

Polysaccharide depolymerization in nature is primarily accomplished by processive glycoside hydrolases (GHs), which abstract single carbohydrate chains from polymer crystals and cleave glycosidic linkages without dissociating after each catalytic event. Understanding the molecular-level features and structural aspects of processivity is of importance due to the prevalence of processive GHs in biomass-degrading enzyme cocktails. Here, we describe recent advances towards the development of a molecular-level theory of processivity for cellulolytic and chitinolytic enzymes, including the development of novel methods for measuring rates of key steps in processive action and insights gained from structural and computational studies. Overall, we present a framework for developing structure-function relationships in processive GHs and outline additional progress towards developing a fundamental understanding of these industrially important enzymes.

Lytic polysaccharide monooxygenases (LPMOs) exhibit a mononuclear copper-containing active site and use dioxygen and a reducing agent to oxidatively cleave glycosidic linkages in polysaccharides. LPMOs represent a unique paradigm in carbohydrate turnover and exhibit synergy with hydrolytic enzymes in biomass depolymerization. To date, several features of copper binding to LPMOs have been elucidated, but the identity of the reactive oxygen species and the key steps in the oxidative mechanism have not been elucidated. Here, density functional theory calculations are used with an enzyme active site model to identify the reactive oxygen species and compare two hypothesized reaction pathways in LPMOs for hydrogen abstraction and polysaccharide hydroxylation; namely, a mechanism that employs a η(1)-superoxo intermediate, which abstracts a substrate hydrogen and a hydroperoxo species is responsible for substrate hydroxylation, and a mechanism wherein a copper-oxyl radical abstracts a hydrogen and subsequently hydroxylates the substrate via an oxygen-rebound mechanism. The results predict that oxygen binds end-on (η(1)) to copper, and that a copper-oxyl-mediated, oxygen-rebound mechanism is energetically preferred. The N-terminal histidine methylation is also examined, which is thought to modify the structure and reactivity of the enzyme. Density functional theory calculations suggest that this posttranslational modification has only a minor effect on the LPMO active site structure or reactivity for the examined steps. Overall, this study suggests the steps in the LPMO mechanism for oxidative cleavage of glycosidic bonds.

We report a novel approach to concurrently improve the tolerance to ionic liquids (ILs) as well as reduce lignin inhibition of Trichoderma reesei cellulase via engineering enzyme charge. Succinylation of the cellulase enzymes led to a nearly twofold enhancement in cellulose conversion in 15% (v/v) 1-butyl-3-methylimidazolium chloride ([BMIM][Cl]). The improvement in activity upon succinylation correlated with the apparent preferential exclusion of the [Cl] anion in fluorescence quenching assays. Additionally, modeling analysis of progress curves of Avicel hydrolysis in buffer indicated that succinylation had a negligible impact on the apparent KM of cellulase. As evidence of reducing lignin inhibition of T. reesei cellulase, succinylation resulted in a greater than twofold increase in Avicel conversion after 170 h in buffer with 1 wt% lignin. The impact of succinylation on lignin inhibition of cellulase further led to the reduction in apparent KM of the enzyme cocktail for Avicel by 2.7-fold. These results provide evidence that naturally evolved cellulases with highly negative surface charge densities may similarly repel lignin, resulting in improved cellulase activity. Ultimately, these results underscore the potential of rational charge engineering as a means of enhancing cellulase function and thus conversion of whole biomass in ILs.

In the last several decades, significant efforts have been conducted to understand the fundamental reactivity of glucose derived from plant biomass in various chemical environments for conversion to renewable fuels and chemicals. For reactions of glucose in water, it is known that inorganic salts naturally present in biomass alter the product distribution in various deconstruction processes. However, the molecular-level interactions of alkali metal ions and glucose are unknown. These interactions are of physiological interest as well, for example, as they relate to cation-glucose cotransport. Here, we employ quantum mechanics (QM) to understand the interaction of a prevalent alkali metal, sodium, with glucose from a structural and thermodynamic perspective. The effect on β-glucose is subtle: a sodium ion perturbs bond lengths and atomic partial charges less than rotating a hydroxymethyl group. In contrast, the presence of a sodium ion significantly perturbs the partial charges of α-glucose anomeric and ring oxygens. Molecular dynamics (MD) simulations provide dynamic sampling in explicit water, and both the QM and the MD results show that sodium ions associate at many positions with respect to glucose with reasonably equivalent propensity. This promiscuous binding nature of Na(+) suggests that computational studies of glucose reactions in the presence of inorganic salts need to ensure thorough sampling of the cation positions, in addition to sampling glucose rotamers. The effect of NaCl on the relative populations of the anomers is experimentally quantified with light polarimetry. These results support the computational findings that Na(+) interacts similarly with α- and β-glucose.

Cellobiohydrolases (CBHs) are typically major components of natural enzyme cocktails for biomass degradation. Their active sites are enclosed in a tunnel, enabling processive hydrolysis of cellulose chains. Glycoside hydrolase Family 6 (GH6) CBHs act from nonreducing ends by an inverting mechanism and are present in many cellulolytic fungi and bacteria. The bacterial Thermobifida fusca Cel6B (TfuCel6B) exhibits a longer and more enclosed active site tunnel than its fungal counterparts. Here, we determine the structures of two TfuCel6B mutants co-crystallized with cellobiose, D274A (catalytic acid), and the double mutant D226A/S232A, which targets the putative catalytic base and a conserved serine that binds the nucleophilic water. The ligand binding and the structure of the active site are retained when compared with the wild type structure, supporting the hypothesis that these residues are directly involved in catalysis. One structure exhibits crystallographic waters that enable construction of a model of the α-anomer product after hydrolysis. Interestingly, the product sites of TfuCel6B are completely enclosed by an “exit loop” not present in fungal GH6 CBHs and by an extended “bottom loop”. From the structures, we hypothesize that either of the loops enclosing the product subsites in the TfuCel6B active site tunnel must open substantially for product release. With simulation, we demonstrate that both loops can readily open to allow product release with equal probability in solution or when the enzyme is engaged on cellulose. Overall, this study reveals new structural details of GH6 CBHs likely important for functional differences among enzymes from this important family.

We examine how the catalytic domain of a glycoside hydrolase family 7 endoglucanase catalytic domain (Cel7B CD) facilitates complexation of cellulose chains from a crystal surface. With direct relevance to the science of biofuel production, this problem also represents a model system of biopolymer processing by proteins in Nature. Interactions of Cel7B CD with a cellulose microfibril along different paths of complexation are characterized by mapping the atomistic fluctuations recorded in free-energy simulations onto the parameters of a coarse-grain model. The resulting patterns of protein-biopolymer couplings also uncover the sequence signatures of the enzyme in peeling off glucan chains from the microfibril substrate. We show that the semiopen active site of Cel7B CD exhibits similar barriers and free energies of complexation over two distinct routes; namely, scooping of a chain into the active-site cleft and threading from the chain end into the channel. On the other hand, the complexation energetics strongly depends on the surface packing of the targeted chain and the resulting interaction sites with the enzyme. A revealed principle is that Cel7B CD facilitates cellulose deconstruction via adaptive coupling to the emergent substrate. The flexible, peripheral segments of the protein outside of the active-site cleft are able to accommodate the varying features of cellulose along the simulated paths of complexation. The general strategy of linking physics-based molecular interactions to protein sequence could also be helpful in elucidating how other protein machines process biopolymers.

Plant cell-wall polysaccharides represent a vast source of food in nature. To depolymerize polysaccharides to soluble sugars, many organisms use multifunctional enzyme mixtures consisting of glycoside hydrolases, lytic polysaccharide mono-oxygenases, polysaccharide lyases, and carbohydrate esterases, as well as accessory, redox-active enzymes for lignin depolymerization. Many of these enzymes that degrade lignocellulose are multimodular with carbohydrate-binding modules (CBMs) and catalytic domains connected by flexible, glycosylated linkers. These linkers have long been thought to simply serve as a tether between structured domains or to act in an inchworm-like fashion during catalytic action. To examine linker function, we performed molecular dynamics (MD) simulations of the Trichoderma reesei Family 6 and Family 7 cellobiohydrolases ( Tr Cel6A and Tr Cel7A, respectively) bound to cellulose. During these simulations, the glycosylated linkers bind directly to cellulose, suggesting a previously unknown role in enzyme action. The prediction from the MD simulations was examined experimentally by measuring the binding affinity of the Cel7A CBM and the natively glycosylated Cel7A CBM-linker. On crystalline cellulose, the glycosylated linker enhances the binding affinity over the CBM alone by an order of magnitude. The MD simulations before and after binding of the linker also suggest that the bound linker may affect enzyme action due to significant damping in the enzyme fluctuations. Together, these results suggest that glycosylated linkers in carbohydrate-active enzymes, which are intrinsically disordered proteins in solution, aid in dynamic binding during the enzymatic deconstruction of plant cell walls.

Dolichyl-P-Man:Man(5)GlcNAc(2)-PP-dolichyl α-1,3-mannosyltransferase (also known as “asparagine-linked glycosylation 3”, or ALG3) is involved in early N -linked glycan synthesis and thus is essential for formation of N -linked protein glycosylation. In this study, we examined the effects of alg3 gene deletion ( alg3Δ ) on growth, development, pigment production, protein secretion and recombinant Trichoderma reesei cellobiohydrolase (rCel7A) expressed in Aspergillus niger . The alg3Δ delayed spore germination in liquid cultures of complete medium (CM), potato dextrose (PD), minimal medium (MM) and CM with addition of cAMP (CM + cAMP), and resulted in significant reduction of hyphal growth on CM, potato dextrose agar (PDA), and CM + cAMP and spore production on CM. The alg3Δ also led to a significant accumulation of red pigment on both liquid and solid CM cultures. The relative abundances of 54 of the total 215 proteins identified in the secretome were significantly altered as a result of alg3Δ , 63% of which were secreted at higher levels in alg3Δ strain than the parent. The rCel7A expressed in the alg3Δ mutant was smaller in size than that expressed in both wild-type and parental strains, but still larger than T. reesei Cel7A. The circular dichroism (CD)-melt scans indicated that change in glycosylation of rCel7A does not appear to impact the secondary structure or folding. Enzyme assays of Cel7A and rCel7A on nanocrystalline cellulose and bleached kraft pulp demonstrated that the rCel7As have improved activities on hydrolyzing the nanocrystalline cellulose. Overall, the results suggest that alg3 is critical for growth, sporulation, pigment production, and protein secretion in A. niger , and demonstrate the feasibility of this alternative approach to evaluate the roles of N -linked glycosylation in glycoprotein secretion and function.

Nature uses a diversity of glycoside hydrolase (GH) enzymes to convert polysaccharides to sugars. As lignocellulosic biomass deconstruction for biofuel production remains costly, natural GH diversity offers a starting point for developing industrial enzymes, and fungal GH family 7 (GH7) cellobiohydrolases, in particular, provide significant hydrolytic potential in industrial mixtures. Recently, GH7 enzymes have been found in other kingdoms of life besides fungi, including in animals and protists. Here, we describe the in vivo spatial expression distribution, properties, and structure of a unique endogenous GH7 cellulase from an animal, the marine wood borer Limnoria quadripunctata (LqCel7B). RT-quantitative PCR and Western blot studies show that LqCel7B is expressed in the hepatopancreas and secreted into the gut for wood degradation. We produced recombinant LqCel7B, with which we demonstrate that LqCel7B is a cellobiohydrolase and obtained four high-resolution crystal structures. Based on a crystallographic and computational comparison of LqCel7B to the well-characterized Hypocrea jecorina GH7 cellobiohydrolase, LqCel7B exhibits an extended substrate-binding motif at the tunnel entrance, which may aid in substrate acquisition and processivity. Interestingly, LqCel7B exhibits striking surface charges relative to fungal GH7 enzymes, which likely results from evolution in marine environments. We demonstrate that LqCel7B stability and activity remain unchanged, or increase at high salt concentration, and that the L. quadripunctata GH mixture generally contains cellulolytic enzymes with highly acidic surface charge compared with enzymes derived from terrestrial microbes. Overall, this study suggests that marine cellulases offer significant potential for utilization in high-solids industrial biomass conversion processes.

Root rot fungi of the Heterobasidion annosum complex are the most damaging pathogens in temperate forests, and the recently sequenced Heterobasidion irregulare genome revealed over 280 carbohydrate-active enzymes. Here, H. irregulare was grown on biomass, and the most abundant protein in the culture filtrate was identified as the only family 7 glycoside hydrolase in the genome, which consists of a single catalytic domain, lacking a linker and carbohydrate-binding module. The enzyme, HirCel7A, was characterized biochemically to determine the optimal conditions for activity. HirCel7A was crystallized and the structure, refined at 1.7 Å resolution, confirms that HirCel7A is a cellobiohydrolase rather than an endoglucanase, with a cellulose-binding tunnel that is more closed than Phanerochaete chrysosporium Cel7D and more open than Hypocrea jecorina Cel7A, suggesting intermediate enzyme properties. Molecular simulations were conducted to ascertain differences in enzyme-ligand interactions, ligand solvation, and loop flexibility between the family 7 glycoside hydrolase cellobiohydrolases from H. irregulare, H. jecorina, and P. chrysosporium. The structural comparisons and simulations suggest significant differences in enzyme-ligand interactions at the tunnel entrance in the −7 to −4 binding sites and suggest that a tyrosine residue at the tunnel entrance of HirCel7A may serve as an additional ligand-binding site. Additionally, the loops over the active site in H. jecorina Cel7A are more closed than loops in the other two enzymes, which has implications for the degree of processivity, endo-initiation, and substrate dissociation. Overall, this study highlights molecular level features important to understanding this biologically and industrially important family of glycoside hydrolases.

Cellulase enzymes deconstruct cellulose to glucose, and are often comprised of glycosylated linkers connecting glycoside hydrolases (GHs) to carbohydrate-binding modules (CBMs). Although linker modifications can alter cellulase activity, the functional role of linkers beyond domain connectivity remains unknown. Here we investigate cellulase linkers connecting GH Family 6 or 7 catalytic domains to Family 1 or 2 CBMs, from both bacterial and eukaryotic cellulases to identify conserved characteristics potentially related to function. Sequence analysis suggests that the linker lengths between structured domains are optimized based on the GH domain and CBM type, such that linker length may be important for activity. Longer linkers are observed in eukaryotic GH Family 6 cellulases compared to GH Family 7 cellulases. Bacterial GH Family 6 cellulases are found with structured domains in either N to C terminal order, and similar linker lengths suggest there is no effect of domain order on length. O-glycosylation is uniformly distributed across linkers, suggesting that glycans are required along entire linker lengths for proteolysis protection and, as suggested by simulation, for extension. Sequence comparisons show that proline content for bacterial linkers is more than double that observed in eukaryotic linkers, but with fewer putative O-glycan sites, suggesting alternative methods for extension. Conversely, near linker termini where linkers connect to structured domains, O-glycosylation sites are observed less frequently, whereas glycines are more prevalent, suggesting the need for flexibility to achieve proper domain orientations. Putative N-glycosylation sites are quite rare in cellulase linkers, while an N-P motif, which strongly disfavors the attachment of N-glycans, is commonly observed. These results suggest that linkers exhibit features that are likely tailored for optimal function, despite possessing low sequence identity. This study suggests that cellulase linkers may exhibit function in enzyme action, and highlights the need for additional studies to elucidate cellulase linker functions.

Cellulases hydrolyze β-1,4 glycosidic linkages in cellulose, which are among the most prevalent and stable bonds in Nature. Cellulases comprise many glycoside hydrolase families and exist as processive or nonprocessive enzymes. Product inhibition negatively impacts cellulase action, but experimental measurements of product-binding constants vary significantly, and there is little consensus on the importance of this phenomenon. To provide molecular level insights into cellulase product inhibition, we examine the impact of product binding on processive and nonprocessive cellulases by calculating the binding free energy of cellobiose to the product sites of catalytic domains of processive and nonprocessive enzymes from glycoside hydrolase families 6 and 7. The results suggest that cellobiose binds to processive cellulases much more strongly than nonprocessive cellulases. We also predict that the presence of a cellodextrin bound in the reactant site of the catalytic domain, which is present during enzymatic catalysis, has no effect on product binding in nonprocessive cellulases, whereas it significantly increases product binding to processive cellulases. This difference in product binding correlates with hydrogen bonding between the substrate-side ligand and the cellobiose product in processive cellulase tunnels and the additional stabilization from the longer tunnel-forming loops. The hydrogen bonds between the substrate- and product-side ligands are disrupted by water in nonprocessive cellulase clefts, and the lack of long tunnel-forming loops results in lower affinity of the product ligand. These findings provide new insights into the large discrepancies reported for binding constants for cellulases and suggest that product inhibition will vary significantly based on the amount of productive binding for processive cellulases on cellulose.

Myriad scientific domains concern themselves with biological electron transfer (ET) events that span across vast scales of rate and efficiency through a remarkably fine-tuned integration of amino acid (AA) sequences, electronic structure, dynamics, and environment interactions. Within this intricate scheme, many questions persist as to how proteins modulate electron-tunneling properties. To help elucidate these principles, we develop a model set of peptides representing the common α-helix and β-strand motifs including all natural AAs within implicit protein-environment solvation. Using an effective Hamiltonian strategy with density functional theory, we characterize the electronic coupling through these peptides, furthermore considering side-chain dynamics. For both motifs, predictions consistently show that backbone-mediated electronic coupling is distinctly sensitive to AA type (aliphatic, polar, aromatic, negatively charged and positively charged), and to side-chain orientation. The unique properties of these residues may be employed to design activated, deactivated, or switch-like superexchange pathways. Electronic structure calculations and Green's function analyses indicate that localized shifts in the electron density along the peptide play a role in modulating these pathways, and further substantiate the experimentally observed behavior of proline residues as superbridges. The distinct sensitivities of tunneling pathways to sequence and conformation revealed in this electronic coupling database help improve our fundamental understanding of the broad diversity of ET reactivity and provide guiding principles for peptide design.

Carbohydrate-binding modules (CBMs) are ubiquitous components of glycoside hydrolases, which degrade polysaccharides in nature. CBMs target specific polysaccharides, and CBM binding affinity to cellulose is known to be proportional to cellulase activity, such that increasing binding affinity is an important component of performance improvement. To ascertain the impact of protein and glycan engineering on CBM binding, we use molecular simulation to quantify cellulose binding of a natively glycosylated Family 1 CBM. To validate our approach, we first examine aromatic-carbohydrate interactions on binding, and our predictions are consistent with previous experiments, showing that a tyrosine to tryptophan mutation yields a 2-fold improvement in binding affinity. We then demonstrate that enhanced binding of 3–6-fold over a nonglycosylated CBM is achieved by the addition of a single, native mannose or a mannose dimer, respectively, which has not been considered previously. Furthermore, we show that the addition of a single, artificial glycan on the anterior of the CBM, with the native, posterior glycans also present, can have a dramatic impact on binding affinity in our model, increasing it up to 140-fold relative to the nonglycosylated CBM. These results suggest new directions in protein engineering, in that modifying glycosylation patterns via heterologous expression, manipulation of culture conditions, or introduction of artificial glycosylation sites, can alter CBM binding affinity to carbohydrates and may thus be a general strategy to enhance cellulase performance. Our results also suggest that CBM binding studies should consider the effects of glycosylation on binding and function.

Proteins employ aromatic residues for carbohydrate binding in a wide range of biological functions. Glycoside hydrolases, which are ubiquitous in nature, typically exhibit tunnels, clefts, or pockets lined with aromatic residues for processing carbohydrates. Mutation of these aromatic residues often results in significant activity differences on insoluble and soluble substrates. However, the thermodynamic basis and molecular level role of these aromatic residues remain unknown. Here, we calculate the relative ligand binding free energy by mutating tryptophans in the Trichoderma reesei family 6 cellulase (Cel6A) to alanine. Removal of aromatic residues near the catalytic site has little impact on the ligand binding free energy, suggesting that aromatic residues immediately upstream of the active site are not directly involved in binding, but play a role in the glucopyranose ring distortion necessary for catalysis. Removal of aromatic residues at the entrance and exit of the Cel6A tunnel, however, dramatically impacts the binding affinity, suggesting that these residues play a role in chain acquisition and product stabilization, respectively. The roles suggested from differences in binding affinity are confirmed by molecular dynamics and normal mode analysis. Surprisingly, our results illustrate that aromatic-carbohydrate interactions vary dramatically depending on the position in the enzyme tunnel. As aromatic-carbohydrate interactions are present in all carbohydrate-active enzymes, these results have implications for understanding protein structure-function relationships in carbohydrate metabolism and recognition, carbon turnover in nature, and protein engineering strategies for biomass utilization. Generally, these results suggest that nature employs aromatic-carbohydrate interactions with a wide range of binding affinities for diverse functions.

At phase boundaries, physical activities of enzymes such as substrate complexation play critical roles in driving biocatalysis. A prominent example is the cellulase cocktails secreted by fungi and bacteria for deconstructing crystalline cellulose in biomass into soluble sugars. At interfaces, molecular mechanisms of the physical steps in biocatalysis remain elusive due to the difficulties of characterizing protein action with high temporal and spatial resolution. Here, we focus on endoglucanase I (Cel7B) from the fungus Trichoderma reesei that hydrolyzes glycosidic bonds on cellulose randomly. We employ all-atom molecular dynamics (MD) simulations to elucidate the interactions of the catalytic domain (CD) of Cel7B with a cellulose microfibril before and after complexing a glucan chain in the binding cleft. The calculated mechanical coupling networks in Cel7B-glucan and Cel7B-microfibril complexes reveal a previously unresolved allosteric coupling at the solid-liquid interface: attachment of the Cel7B CD to the cellulose surface affects glucan chain clenching in the binding cleft. Alternative loop segments of the Cel7B CD were found to affix to intact or defective surface structures on the microfibril, depending on the complexation state. From a multiple sequence alignment, residues in surface-affixing segments show strong conservation, highlighting the functional importance of the physical activities that they facilitate. Surface-affixing residues also demonstrate significant sequence correlation with active-site residues, revealing the functional connection between complexation and hydrolysis. Analysis of the Cel7B CD exemplifies that the mechanical coupling networks calculated from atomistic MD simulations can be used to capture the conservation and correlation in sequence alignment.

Chitin is the primary structural material of insect and crustacean exoskeletons and fungal and algal cell walls, and as such it is the one of the most abundant biological materials on Earth. Chitin forms linear polymers of β1,4-linked-N-acetyl-D-glucosamine (GlcNAc), and in Nature, enzyme cocktails deconstruct chitin to GlcNAc. The mechanism of chitin deconstruction, like that of cellulose deconstruction, has been under investigation due to its importance in the global carbon cycle and in production of renewable and sustainable products from biological matter. To further understand the nanoscale properties of chitin, here we simulate crystals of α-chitin, which is the most prevalent form in Nature. We find excellent agreement with the recently reported crystal structure and we report the salient features of the simulations related to crystalline stability. We also compute the thermodynamic work required to peel individual chains from α-chitin surfaces, which a chitinase enzyme must conduct to deconstruct chitin. Compared with previous simulations of native plant cellulose Iβ, α-chitin exhibits higher decrystallization work for chains in the middle of surfaces and similar work for chains on the edges of crystals. Unlike cellulose, the free energy profile is dominated by a single bifurcated hydrogen bond between chains formed by the GlcNAc side chains and the O6 atoms on the primary alcohol group. This study highlights the molecular features of chitin that make it such a tough, recalcitrant material, and provides a key thermodynamic parameter in our quantitative understanding of how enzymes contribute to the turnover of carbohydrates in the biosphere.

We use molecular simulation to elucidate the structural behavior of small hydrated cellulose Iβ microfibrils heated to 227 °C (500 K) with two carbohydrate force fields. In contrast to the characteristic two-dimensional hydrogen-bonded layer sheets present in the cellulose Iβ crystal structure, we show that at high temperature a three-dimensional hydrogen bond network forms, made possible by hydroxymethyl groups changing conformation from trans-gauche (TG) to gauche-gauche (GG) in every second layer corresponding to "center" chains in cellulose Iβ and from TG to gauche-trans (GT) in the "origin" layer. The presence of a regular three-dimensional hydrogen bond network between neighboring sheets eliminates the possibility of twist, whereas two-dimensional hydrogen bonding allows for microfibril twist to occur. Structural features of this high-temperature phase as determined by molecular simulation may explain several experimental observations for which no detailed structural basis has been offered. This includes an explanation for the observed temperature and crystal size dependence for the extent of hydrogen/deuterium exchange, and diffraction patterns of cellulose at high temperature.

Most bacteria use free enzymes to degrade plant cell walls in nature. However, some bacteria have adopted a different strategy wherein enzymes can either be free or tethered on a protein scaffold forming a complex called a cellulosome. The study of the structure and mechanism of these large macromolecular complexes is an active and ongoing research topic, with the goal of finding ways to improve biomass conversion using cellulosomes. Several mechanisms involved in cellulosome formation remain unknown, including how cellulosomal enzymes assemble on the scaffoldin and what governs the population of cellulosomes created during self-assembly. Here, we present a coarse-grained model to study the self-assembly of cellulosomes. The model captures most of the physical characteristics of three cellulosomal enzymes (Cel5B, CelS, and CbhA) and the scaffoldin (CipA) from Clostridium thermocellum. The protein structures are represented by beads connected by restraints to mimic the flexibility and shapes of these proteins. From a large simulation set, the assembly of cellulosomal enzyme complexes is shown to be dominated by their shape and modularity. The multimodular enzyme, CbhA, binds statistically more frequently to the scaffoldin than CelS or Cel5B. The enhanced binding is attributed to the flexible nature and multimodularity of this enzyme, providing a longer residence time around the scaffoldin. The characterization of the factors influencing the cellulosome assembly process may enable new strategies to create designers cellulosomes.

We probe the molecular-level behavior of the Family 1 carbohydrate-binding module (CBM) from a commonly studied fungal cellulase, the Family 7 cellobiohydrolase (Cel7A) from Trichoderma reesei, on the hydrophobic face of crystalline cellulose. With a fully atomistic model, we predict that the CBM alone exhibits regions of thermodynamic stability along a cellulose chain corresponding to a cellobiose unit, which is the catalytic product of the entire Cel7A enzyme. In addition, we determine which residues and the types of interactions that are responsible for the observed processivity length scale of the CBM: Y5, Q7, N29, and Y32. These results imply that the CBM can anchor the Cel7A enzyme at discrete points along a cellulose chain and thus aid in both recognizing cellulose chain ends for initial attachment to cellulose as well as aid in enzymatic catalysis by diffusing between stable wells on a length scale commensurate with the catalytic, processive cycle of Cel7A during cellulose hydrolysis. Comparison of other Family 1 CBMs show high functional homology to the four amino acids responsible for the processivity length scale on the surface of crystalline cellulose, which suggests that Family 1 CBMs may generally employ this type of approach for translation on the cellulose surface. Overall, this work provides further insight into the molecular-level mechanisms by which a CBM recognizes and interacts with cellulose.

This study applies aimless shooting and likelihood maximization to determine the molecular mechanism in the solid state polymorph transformation in terephthalic acid from over 500 candidate order parameters. The crystals examined here extend the range of crystal sizes considered in our previous work (J. Amer. Chem. Soc. 2007, 129, 4714) and reveal a change in the mechanism with increasing system size. As the crystal size increases beyond that studied in our previous work, the polymorph transformation mechanism changes from a global distortion of the crystal to a local corner nucleation mechanism. In the corner nucleation mechanism, the interfacial area between the two polymorphs is minimized for a given nucleus size. However, this mechanism differs from classical nucleation theory in that the molecular level details are essential to describe the nucleation process, which involves nonspherical domains at the corner of the crystal. These new findings suggest that there is a range of sizes for which corner nucleation is the dominant mechanism of polymorph transitions, thus implying that different mechanistic regimes exist for nucleation based on crystal size. From a computational standpoint, this study demonstrates the utility of aimless shooting and likelihood maximization to identify nonintuitive reaction coordinates in complex systems.

This paper extends our previous work on obtaining reaction coordinates from aimless shooting and likelihood maximization. We introduce a simplified version of aimless shooting and a half-trajectory likelihood score based on the committor probability. Additionally, we analyze and compare the absolute log-likelihood score for perfect and approximate reaction coordinates. We also compare the aimless shooting and likelihood maximization approach to the earlier genetic neural network (GNN) approach of Ma and Dinner [J. Phys. Chem. B 109, 6769 (2005)]. For a fixed number of total trajectories in the GNN approach, the accuracy of the transition state ensemble decreases as the number of trajectories per committor probability estimate increases. This quantitatively demonstrates the benefit of individual committor probability realizations over committor probability estimates. Furthermore, when the least squares score of the GNN approach is applied to individual committor probability realizations, the likelihood score still provides a better approximation to the true transition state surface. Finally, the polymorph transition in terephthalic acid demonstrates that the new half-trajectory likelihood scheme estimates the transition state location more accurately than likelihood schemes based on the probability of being on a transition path.

A molecular mechanism for nucleation for the solid-state polymorph transformation of terephthalic acid is presented. New methods recently developed in our group, aimless shooting and likelihood maximization, are employed to construct a model for the reaction coordinate for the two system sizes studied. The reaction coordinate approximation is validated using the committor probability analysis. The transformation proceeds via a localized, elongated nucleus along the crystal edge formed by fluctuations in the supramolecular synthons, suggesting a nucleation and growth mechanism in the macroscopic system.

Methane clathrate hydrate nucleation and growth is investigated via analysis of molecular dynamics simulations using a new order parameter. This order parameter (OP), named the Mutually Coordinated Guest (MCG) OP, quantifies the appearance and connectivity of molecular clusters composed of guests separated by water clusters. It is the first two-component OP used for quantifying hydrate nucleation and growth. The algorithm for calculating the MCG OP is described in detail. Its physical motivation and advantages compared to existing methods are discussed.