Search

The complex structure of plant cell walls resists chemical or biological degradation, challenging the breakdown of lignocellulosic biomass into renewable chemical precursors that could form the basis of future production of green chemicals and transportation fuels. Here, experimental and computational results reveal that the effect of the tetrahydrofuran (THF)-water cosolvents on the structure of lignin and on its interactions with cellulose in the cell wall drives multiple synergistic mechanisms leading to the efficient breakdown and fractionation of biomass into valuable chemical precursors. Molecular simulations show that THF-water is an excellent "theta" solvent, such that lignin dissociates from itself and from cellulose and expands to form a random coil. The expansion of the lignin molecules exposes interunit linkages, rendering them more susceptible to depolymerization by acid-catalyzed cleavage of aryl-ether bonds. Nanoscale infrared sensors confirm cosolvent-mediated molecular rearrangement of lignin in the cell wall of micrometer-thick hardwood slices and track the disappearance of lignin. At bulk scale, adding dilute acid to the cosolvent mixture liberates the majority of the hemicellulose and lignin from biomass, allowing unfettered access of cellulolytic enzymes to the remaining cellulose-rich material, allowing them to sustain high rates of hydrolysis to glucose without enzyme deactivation. Through this multiscale analysis, synergistic mechanisms for biomass deconstruction are identified, portending a paradigm shift toward first-principles design and evaluation of other cosolvent methods to realize low cost fuels and bioproducts.

Pretreatment facilitates more complete deconstruction of plant biomass to enable more economic production of lignocellulosic biofuels and byproducts. Various co-solvent pretreatments have demonstrated advantages relative to aqueous-only methods by enhancing lignin removal to allow unfettered access to cellulose. However, there is a limited mechanistic understanding of the interactions between the co-solvents and cellulose that impedes further improvement of such pretreatment methods. Recently, tetrahydrofuran (THF) has been identified as a highly effective co-solvent for the pretreatment and fractionation of biomass. To elucidate the mechanism of the THF-water interactions with cellulose, we pair simulation and experimental data demonstrating that enhanced solubilization of cellulose can be achieved by the THF-water co-solvent system at equivolume mixtures and moderate temperatures (≤445 K). The simulations show that THF and water spontaneously phase separate on the local surface of a cellulose fiber, owing to hydrogen bonding of water molecules with the hydrophilic cellulose faces and stacking of THF molecules on the hydrophobic faces. Furthermore, a single fully solvated cellulose chain is shown to be preferentially bound by water molecules in the THF-water mixture. In light of these findings, co-solvent reactions were performed on microcrystalline cellulose and maple wood to show that THF significantly enhanced cellulose deconstruction and lignocellulose solubilization at simulation conditions, enabling a highly versatile and efficient biomass pretreatment and fractionation method.

A cellulose synthesis complex with a "rosette" shape is responsible for synthesis of cellulose chains and their assembly into microfibrils within the cell walls of land plants and their charophyte algal progenitors. The number of cellulose synthase proteins in this large multisubunit transmembrane protein complex and the number of cellulose chains in a microfibril have been debated for many years. This work reports a low resolution structure of the catalytic domain of CESA1 from Arabidopsis (Arabidopsis thaliana; AtCESA1CatD) determined by small-angle scattering techniques and provides the first experimental evidence for the self-assembly of CESA into a stable trimer in solution. The catalytic domain was overexpressed in Escherichia coli, and using a two-step procedure, it was possible to isolate monomeric and trimeric forms of AtCESA1CatD. The conformation of monomeric and trimeric AtCESA1CatD proteins were studied using small-angle neutron scattering and small-angle x-ray scattering. A series of AtCESA1CatD trimer computational models were compared with the small-angle x-ray scattering trimer profile to explore the possible arrangement of the monomers in the trimers. Several candidate trimers were identified with monomers oriented such that the newly synthesized cellulose chains project toward the cell membrane. In these models, the class-specific region is found at the periphery of the complex, and the plant-conserved region forms the base of the trimer. This study strongly supports the "hexamer of trimers" model for the rosette cellulose synthesis complex that synthesizes an 18-chain cellulose microfibril as its fundamental product., (© 2016 American Society of Plant Biologists. All Rights Reserved.)

One-dimensional (1D) (spherically averaged) powder diffraction diagrams are commonly used to determine the degree of cellulose crystallinity in biomass samples. Here, it is shown using molecular modeling how disorder in cellulose fibrils can lead to considerable uncertainty in conclusions drawn concerning crystallinity based on 1D powder diffraction data alone. For example, cellulose microfibrils that contain both crystalline and noncrystalline segments can lead to powder diffraction diagrams lacking identifiable peaks, while microfibrils without any crystalline segments can lead to such peaks. This leads to false positives, that is, assigning disordered cellulose as crystalline, and false negatives, that is, categorizing fibrils with crystalline segments as amorphous. The reliable determination of the fraction of crystallinity in any given biomass sample will require a more sophisticated approach combining detailed experiment and simulation., (© 2014 Wiley Periodicals, Inc.)

The mechanical and dynamical properties of cellulose, the most abundant biomolecule on earth, are essential for its function in plant cell walls and advanced biomaterials. Cellulose is almost always found in a hydrated state, and it is therefore important to understand how hydration influences its dynamics and mechanics. Here, the nanosecond-time scale dynamics of cellulose is characterized using dynamic neutron scattering experiments and molecular dynamics (MD) simulation. The experiments reveal that hydrated samples exhibit a higher average mean-square displacement above ∼240 K. The MD simulation reveals that the fluctuations of the surface hydroxymethyl atoms determine the experimental temperature and hydration dependence. The increase in the conformational disorder of the surface hydroxymethyl groups with temperature follows the cellulose persistence length, suggesting a coupling between structural and mechanical properties of the biopolymer. In the MD simulation, 20% hydrated cellulose is more rigid than the dry form, due to more closely packed cellulose chains and water molecules bridging cellulose monomers with hydrogen bonds. This finding may have implications for understanding the origin of strength and rigidity of secondary plant cell walls. The detailed characterization obtained here describes how hydration-dependent increased fluctuations and hydroxymethyl disorder at the cellulose surface lead to enhancement of the rigidity of this important biomolecule.

Ionic liquids have become a popular solvent for cellulose pretreatment in biorefineries due to their efficiency in dissolution and their reusability. Understanding the interactions between cations, anions, and cellulose is key to the development of better solvents and the improvement of pretreatment conditions. While previous studies described the interactions between ionic liquids and cellulose fibers, shedding light on the initial stages of the cellulose dissolution process, we study the end state of that process by exploring the structure and dynamics of a single cellulose decamer solvated in 1-butyl-3-methyl-imidazolium chloride (BmimCl) and in water using replica-exchange molecular dynamics. In both solvents, global structural features of the cellulose chain are similar. However, analyses of local structural properties show that cellulose explores greater conformational variability in the ionic liquid than in water. For instance, in BmimCl the cellulose intramolecular hydrogen bond O3H'···O5 is disrupted more often resulting in greater flexibility of the solute. Our results indicate that the cellulose chain is more dynamic in BmimCl than in water, which may play a role in the favorable dissolution of cellulose in the ionic liquid. Calculation of the configurational entropy of the cellulose decamer confirms its higher conformational flexibility in BmimCl than in water at elevated temperatures.

Understanding biomass structure and dynamics on multiple time and length scales is important for the development of cellulosic biofuels. To this aim, we have developed a coarse-grain (CG) model for molecular dynamics (MD) simulations of cellulose fibrils in explicit water based on target observables from fully atomistic simulations. This model examines the significance of the presence of explicit solvent and compares results with the previous, implicit solvent CG cellulose models. The present, constraint-free CG model is used to generate a series of noncrystalline fibril structures using a coupling parameter, λ, between fully crystalline and fully amorphous potentials. By exploring various structural parameters, including the root-mean-square deviation, root-mean-square fluctuations, radius of gyration, and persistence length, we find the crystalline-to-amorphous state transition takes place at λ ≈ 0.386. The persistence length of cellulose fibril in the transition region corresponds to that of native cellulose fibrils. The transition between crystalline and amorphous fibrils occurs at larger values of λ in explicit water than in the implicit case. Detailed analysis of individual energetic contribution to the transition reveals that the nonbonded interactions, in particular, that of cellulose-water interaction, plays a significant role in the observed crystalline to amorphous transition of cellulose fibril. The present study thus highlights the importance of solvent presence that cannot be adequately described with the previous implicit solvent model. The present method provides an accurate and constraint-free approach for deriving a variety of structures of cellulose in water, with a wide range of crystallinity, suitable for incorporation into large-scale models of lignocellulosic biomass.

The precipitation of lignin onto cellulose after pretreatment of lignocellulosic biomass is an obstacle to economically viable cellulosic ethanol production. Here, 750 ns nonequilibrium molecular dynamics simulations are reported of a system of lignin and cellulose in aqueous solution. Lignin is found to strongly associate with itself and the cellulose. However, noncrystalline regions of cellulose are observed to have a lower tendency to associate with lignin than crystalline regions, and this is found to arise from stronger hydration of the noncrystalline chains. The results suggest that the recalcitrance of crystalline cellulose to hydrolysis arises not only from the inaccessibility of inner fibers but also due to the promotion of lignin adhesion.

Biomass recalcitrance, the resistance of cellulosic biomass to degradation, is due in part to the stability of the hydrogen bond network and stacking forces between the polysaccharide chains in cellulose microfibers. The fragment molecular orbital (FMO) method at the correlated Møller-Plesset second order perturbation level of theory was used on a model of the crystalline cellulose Iα core with a total of 144 glucose units. These computations show that the intersheet chain interactions are stronger than the intrasheet chain interactions for the crystalline structure, while they are more similar to each other for a relaxed structure. An FMO chain pair interaction energy decomposition analysis for both the crystal and relaxed structures reveals an intricate interplay between electrostatic, dispersion, charge transfer, and exchange repulsion effects. The role of the primary alcohol groups in stabilizing the interchain hydrogen bond network in the inner sheet of the crystal and relaxed structures of cellulose Iα, where edge effects are absent, was analyzed. The maximum attractive intrasheet interaction is observed for the GT-TG residue pair with one intrasheet hydrogen bond, suggesting that the relative orientation of the residues is as important as the hydrogen bond network in strengthening the interaction between the residues.

A molecular level understanding of the structure, dynamics and mechanics of cellulose fibers can aid in understanding the recalcitrance of biomass to hydrolysis in cellulosic biofuel production. Here, a residue-scale REACH (Realistic Extension Algorithm via Covariance Hessian) coarse-grained force field was derived from all-atom molecular dynamics (MD) simulations of the crystalline Iβ cellulose fibril. REACH maps the atomistic covariance matrix onto coarse-grained elastic force constants. The REACH force field was found to reproduce the positional fluctuations and low-frequency vibrational spectra from the all-atom model, allowing elastic properties of the cellulose fibril to be characterized using the coarse-grained force field with a speedup of >20 relative to atomistic MD on systems of the same size. The calculated longitudinal/transversal Young's modulus and the velocity of sound are in agreement with experiment. The persistence length of a 36-chain cellulose microcrystal was estimated to be ~380 μm. Finally, the normal-mode analysis with the REACH force field suggests that intrinsic dynamics might facilitate the deconstruction of the cellulose fibril from the hydrophobic surface.

The use of ionic liquids for non-derivatized cellulose dissolution promises an alternative method for the thermochemical pretreatment of biomass that may be more efficient and environmentally acceptable than more conventional techniques in aqueous solution. Here, we performed equilibrium MD simulations of a cellulose microfibril in the ionic liquid 1-butyl-3-methylimidazolium chloride (BmimCl) and compared the solute structure and the solute-solvent interactions at the interface with those from corresponding simulations in water. The results indicate a higher occurrence of solvent-exposed orientations of cellulose surface hydroxymethyl groups in BmimCl than in water. Moreover, spatial and radial distribution functions indicate that hydrophilic surfaces are a preferred site of interaction between cellulose and the ionic liquid. In particular, hydroxymethyl groups on the hydrophilic fiber surface adopt a different conformation from their counterparts oriented towards the fiber's core. Furthermore, the glucose units with these solvent-oriented hydroxymethyls are surrounded by the heterocyclic organic cation in a preferred parallel orientation, suggesting a direct and distinct interaction scheme between cellulose and BmimCl.

Cellulosic biomass has the potential to serve as a major renewable energy source. However, its strong recalcitrance to degradation hampers its large-scale use in biofuel production. To overcome this problem, a detailed understanding of the origins of the recalcitrance is required. One main biophysical phenomenon leading to the recalcitrance is the high structural ordering of natural cellulose fibrils, that arises largely from an extensive hydrogen-bond network between and within cellulose polymers. Here, we present a lattice-based model of cellulose I(α), one of the two major natural forms, at the resolution of explicit hydrogen bonds. The partition function and thermodynamic properties are evaluated using the transfer matrix method. Two competing hydrogen-bond patterns are found. This plasticity of the hydrogen-bond network leads to an entropic contribution stabilizing the crystalline fibril at intermediate temperatures. At these temperatures, an enhanced probability of bonding between the individual cellulose chains gives rise to increased resistance of the entire cellulose fibril to degradation, before the final disassembly temperature is reached. The results are consistent with the available crystallographic and IR spectroscopic experiments on the thermostability of cellulose I(α)., (© 2011 American Institute of Physics)

The hydrolysis of cellulose is the bottleneck in cellulosic ethanol production. The cellobiohydrolase CelS from Clostridium thermocellum catalyzes the hydrolysis of cello-oligosaccharides via inversion of the anomeric carbon. Here, to examine key features of the CelS-catalyzed reaction, QM/MM (SCCDFTB/MM) simulations are performed. The calculated free energy profile for the reaction possesses a 19 kcal/mol barrier. The results confirm the role of active site residue Glu87 as the general acid catalyst in the cleavage reaction and show that Asp255 may act as the general base. A feasible position in the reactant state of the water molecule responsible for nucleophilic attack is identified. Sugar ring distortion as the reaction progresses is quantified. The results provide a computational approach that may complement the experimental design of more efficient enzymes for biofuel production.

A particularly promising approach to deconstructing and fractionating lignocellulosic biomass to produce green renewable fuels and high-value chemicals pretreats the biomass with organic solvents in aqueous solution. Here, neutron scattering and moleculardynamics simulations reveal the temperature-dependent morphological changes in poplar wood biomass during tetrahydrofuran (THF):water pretreatment and provide a mechanism by which the solvent components drive efficient biomass breakdown. Whereas lignin dissociates over a wide temperature range (>25 °C) cellulose disruption occurs only above 150 °C. Neutron scattering with contrast variation provides direct evidence for the formation of THF-rich nanoclusters (Rg ~ 0.5 nm) on the nonpolar cellulose surfaces and on hydrophobic lignin, and equivalent waterrich nanoclusters on polar cellulose surfaces. The disassembly of the amphiphilic biomass is thus enabled through the local demixing of highly functional cosolvents, THF and water, which preferentially solvate specific biomass surfaces so as to match the local solute polarity. A multiscale description of the efficiency of THF:water pretreatment is provided: matching polarity at the atomic scale prevents lignin aggregation and disrupts cellulose, leading to improvements in deconstruction at the macroscopic scale. [ABSTRACT FROM AUTHOR]

It has been previously shown that cellulose-lignin droplets’ strong interactions, resulting from lignin coalescence and redisposition on cellulose surface during thermochemical pretreatments, increase cellulose recalcitrance to biological conversion, especially at commercially viable low enzyme loadings. However, information on the impact of cellulose–hemicellulose interactions on cellulose recalcitrance following relevant pretreatment conditions are scarce. Here, to investigate the effects of plausible hemicellulose precipitation and re-association with cellulose on cellulose conversion, different pretreatments were applied to pure Avicel® PH101 cellulose alone and Avicel mixed with model hemicellulose compounds followed by enzymatic hydrolysis of resulting solids at both low and high enzyme loadings. Solids produced by pretreatment of Avicel mixed with hemicelluloses (AMH) were found to contain about 2 to 14.6% of exogenous, precipitated hemicelluloses and showed a remarkably much lower digestibility (up to 60%) than their respective controls. However, the exogenous hemicellulosic residues that associated with Avicel following high temperature pretreatments resulted in greater losses in cellulose conversion than those formed at low temperatures, suggesting that temperature plays a strong role in the strength of cellulose–hemicellulose association. Molecular dynamics simulations of hemicellulosic xylan and cellulose were found to further support this temperature effect as the xylan–cellulose interactions were found to substantially increase at elevated temperatures. Furthermore, exogenous, precipitated hemicelluloses in pretreated AMH solids resulted in a larger drop in cellulose conversion than the delignified lignocellulosic biomass containing comparably much higher natural hemicellulose amounts. Increased cellulase loadings or supplementation of cellulase with xylanases enhanced cellulose conversion for most pretreated AMH solids; however, this approach was less effective for solids containing mannan polysaccharides, suggesting stronger association of cellulose with (hetero) mannans or lack of enzymes in the mixture required to hydrolyze such polysaccharides. [ABSTRACT FROM AUTHOR]

Background: The conversion of plant biomass to ethanol via enzymatic cellulose hydrolysis offers a potentially sustainable route to biofuel production. However, the inhibition of enzymatic activity in pretreated biomass by lignin severely limits the efficiency of this process. Results: By performing atomic-detail molecular dynamics simulation of a biomass model containing cellulose, lignin, and cellulases (TrCel7A), we elucidate detailed lignin inhibition mechanisms. We find that lignin binds preferentially both to the elements of cellulose to which the cellulases also preferentially bind (the hydrophobic faces) and also to the specific residues on the cellulose-binding module of the cellulase that are critical for cellulose binding of TrCel7A (Y466, Y492, and Y493). Conclusions: Lignin thus binds exactly where for industrial purposes it is least desired, providing a simple explanation of why hydrolysis yields increase with lignin removal. [ABSTRACT FROM AUTHOR]

Lignocellulosic biomass, a potentially important renewable organic source of energy and chemical feedstock, resists degradation to glucose in industrial hydrolysis processes and thus requires expensive thermochemical pretreatments. Understanding the mechanism of biomass breakdown during these pretreatments will lead to more efficient use of biomass. By combining multiple probes of structure, sensitive to different length scales, with molecular dynamics simulations, we reveal two fundamental processes responsible for the morphological changes in biomass during steam explosion pretreatment: cellulose dehydration and ligninhemicellulose phase separation. We further show that the basic driving forces are the same in other leading thermochemical pretreatments, such as dilute acid pretreatment and ammonia fiber expansion. [ABSTRACT FROM AUTHOR]