Solvent-Driven Preferential Association of Lignin with Regions of Crystalline Cellulose in Molecular Dynamics Simulation

Advancing plant cell wall modelling: Atomistic insights into cellulose, disordered cellulose, and hemicelluloses - A review

The complexity of plant cell walls on different hierarchical levels still impedes the detailed understanding of biosynthetic pathways, interferes with processing in industry and finally limits applicability of cellulose materials. While there exist many challenges to readily accessing these hierarchies at (sub-) angström resolution, the development of advanced computational methods has the potential to unravel important questions in this field. Here, we summarize the contributions of molecular dynamics simulations in advancing the understanding of the physico-chemical properties of natural fibres. We aim to present a comprehensive view of the advancements and insights gained from molecular dynamics simulations in the field of carbohydrate polymers research. The review holds immense value as a vital reference for researchers seeking to undertake atomistic simulations of plant cell wall constituents. Its significance extends beyond the realm of molecular modeling and chemistry, as it offers a pathway to develop a more profound comprehension of plant cell wall chemistry, interactions, and behavior. By delving into these fundamental aspects, the review provides invaluable insights into future perspectives for exploration. Researchers within the molecular modeling and carbohydrates community can greatly benefit from this resource, enabling them to make significant strides in unraveling the intricacies of plant cell wall dynamics.

Lignin self-assembly phenomena and valorization strategies for pulping, biorefining, and materials development: Part 1. The physical chemistry of lignin self-assembly

Physical chemistry aspects are emphasized in this comprehensive review of self-assembly phenomena involving lignin in various forms. Attention to this topic is justified by the very high availability, low cost, and renewable nature of lignin, together with opportunities to manufacture diverse products, for instance, polymers/resins, bioplastics, carbon fibers, bio-asphalt, sunscreen components, hydrophobic layers, and microcapsules. The colloidal lignin material, nanoparticles, and microstructures that can be formed as a result of changes in solvent properties, pH, or other adjustments to a suspending medium have been shown to depend on many factors. Such factors are examined in this work based on the concepts of self-assembly, which can be defined as an organizing principle dependent on specific attributes of the starting entities themselves. As a means to promote such concepts and to facilitate further development of nano-scale lignin products, this article draws upon evidence from a wide range of studies. These include investigations of many different plant sources of lignin, processes of delignification, solvent systems, anti-solvent systems or other means of achieving phase separation, and diverse means of achieving colloidal stability (if desired) of resulting self-assembled lignin structures. Knowledge of the self-organization behavior of lignin can provide significant structural information to optimize the use of lignin in value-added applications. Examples include chemical conditions and preparation procedures in which lignin-related compounds of particles organize themselves as spheres, hollow spheres, surface-bound layers, and a variety of other structures. Published articles show that such processes can be influenced by the selection of lignin type, pulping or extraction processes, functional groups such as phenolic, carboxyl, and sulfonate, chemical derivatization reactions, solvent applications, aqueous conditions, and physical processes, such as agitation. Precipitation from non-aqueous solutions represents a key focus of lignin self-assembly research. The review also considers stabilization mechanisms of self-assembled lignin-related structures.

Molecular insights into inhibiting effects of lignin on cellulase investigated by molecular dynamics simulation

The hydrolysis of lignocellulose into fermentable monosaccharides using cellulases represents a critical stage in lignocellulosic bioconversion. However, the inactivation of cellulase in the presence of lignin is attributed to the high cost of biofinery. To address this challenge, a comprehensive investigation into the structure-function relationship underlying lignin-driven cellulase inactivation is essential. In this study, molecular docking and molecular dynamics (MD) simulations were employed to explore the impacts of lignin fragments on the catalytic efficiency of cellulase at the atomic level. The findings revealed that soluble lignin fragments and cellulose could spontaneously form stable complexes with cellulase, indicating a competitive binding scenario. The enzyme's structure remained unchanged upon binding to lignin. Furthermore, specific amino acid residues have been identified as involved in interactions with lignin and cellulose. Hydrophobic interactions were found to dominate the binding of lignin to cellulase. Based on the mechanisms underlying the interactions between lignin fragments and cellulase, decreased hydrophobicity and change in the charge of lignin may mitigate the inhibition of cellulase. Furthermore, site mutations and chemical modification are also feasible to improve the efficiency of cellulase. This study may contribute valuable insights into the design of more lignin-resistant enzymes and the optimization of lignocellulosic pretreatment technologies.

How lignin sticks to cellulose-insights from atomic force microscopy enhanced by machine-learning analysis and molecular dynamics simulations

Elucidating cellulose-lignin interactions at the molecular and nanometric scales is an important research topic with impacts on several pathways of biomass valorization. Here, the interaction forces between a cellulosic substrate and lignin are investigated. Atomic force microscopy with lignin-coated tips is employed to probe the site-specific adhesion to a cellulose film in liquid water. Over seven thousand force-curves are analyzed by a machine-learning approach to cluster the experimental data into types of cellulose-tip interactions. The molecular mechanisms for distinct types of cellulose-lignin interactions are revealed by molecular dynamics simulations of lignin globules interacting with different cellulose Iβ crystal facets. This unique combination of experimental force-curves, data-driven analysis, and molecular simulations opens a new approach of investigation and updates the understanding of cellulose-lignin interactions at the nanoscale.

Diffusion of volatile organics and water in the epicuticular waxes of petunia petal epidermal cells

Plant cuticles are a mixture of crystalline and amorphous waxes that restrict the exchange of molecules between the plant and the atmosphere. The multicomponent nature of cuticular waxes complicates the study of the relationship between the physical and transport properties. Here, a model cuticle based on the epicuticular waxes of Petunia hybrida flower petals was formulated to test the effect of wax composition on diffusion of water and volatile organic compounds (VOCs). The model cuticle was composed of an n-tetracosane (C₂₄ H₅₀ ), 1-docosanol (C₂₂ H₄₅ OH), and 3-methylbutyl dodecanoate (C₁₇ H₃₄ O₂ ), reflecting the relative chain length, functional groups, molecular arrangements, and crystallinity of the natural waxes. Molecular dynamics simulations were performed to obtain diffusion coefficients for compounds moving through waxes of varying composition. Simulated VOC diffusivities of the model system were found to highly correlate with in vitro measurements in isolated petunia cuticles. VOC diffusivity increased up to 30-fold in completely amorphous waxes, indicating a significant effect of crystallinity on cuticular permeability. The crystallinity of the waxes was highly dependent on the elongation of the lattice length and decrease in gap width between crystalline unit cells. Diffusion of water and higher molecular weight VOCs were significantly affected by alterations in crystalline spacing and lengths, whereas the low molecular weight VOCs were less affected. Comparison of measured diffusion coefficients from atomistic simulations and emissions from petunia flowers indicates that the role of the plant cuticle in the VOC emission network is attributed to the differential control on mass transfer of individual VOCs by controlling the composition, amount, and dynamics of scent emission.

Insight into the dual effect of water on lignin dissolution in ionic liquids

The dual regulation of water on lignin in ionic liquids was studied at the molecular level by molecular dynamics simulation. The simulation results show that a small amount of water will destroy the ion association in ionic liquids, that is, it will produce more free anions and cations. The free ions around lignin are conducive to the dissolution of lignin. On the contrary, excess water will seriously solvate anions and cations. By changing the number of lignin clusters, it is more intuitive to observe that the dissolution of lignin in ILs containing a small amount of water is stronger than that in pure IL, however, the dissolution ability of lignin is reduced after adding a large amount of water in ILs. It is concluded that with the increase of water content, water changes from co-solvent to anti-solvent in the dissolution process. This study provides ideas for the design of IL-water system for economic pretreatment of biomass.

Physico-chemical interactions within lignocellulosic biomass and their importance in developing solvent based deconstruction methods

Fundamental understanding of physico-chemical interactions among the biopolymers in lignocellulosic biomass is crucial to develop atom-efficient deconstruction methods.

Unlocking the secret of lignin-enzyme interactions: Recent advances in developing state-of-the-art analytical techniques

Bioconversion of renewable lignocellulosics to produce liquid fuels and chemicals is one of the most effective ways to solve the problem of fossil resource shortage, energy security, and environmental challenges. Among the many biorefinery pathways, hydrolysis of lignocellulosics to fermentable monosaccharides by cellulase is arguably the most critical step of lignocellulose bioconversion. In the process of enzymatic hydrolysis, the direct physical contact between enzymes and cellulose is an essential prerequisite for the hydrolysis to occur. However, lignin is considered one of the most recalcitrant factors hindering the accessibility of cellulose by binding to cellulase unproductively, which reduces the saccharification rate and yield of sugars. This results in high costs for the saccharification of carbohydrates. The various interactions between enzymes and lignin have been explored from different perspectives in literature, and a basic lignin inhibition mechanism has been proposed. However, the exact interaction between lignin and enzyme as well as the recently reported promotion of some types of lignin on enzymatic hydrolysis is still unclear at the molecular level. Multiple analytical techniques have been developed, and fully unlocking the secret of lignin-enzyme interactions would require a continuous improvement of the currently available analytical techniques. This review summarizes the current commonly used advanced research analytical techniques for investigating the interaction between lignin and enzyme, including quartz crystal microbalance with dissipation (QCM-D), surface plasmon resonance (SPR), attenuated total reflectance-Fourier transform infrared (ATR-FTIR) spectroscopy, atomic force microscopy (AFM), nuclear magnetic resonance (NMR) spectroscopy, fluorescence spectroscopy (FLS), and molecular dynamics (MD) simulations. Interdisciplinary integration of these analytical methods is pursued to provide new insight into the interactions between lignin and enzymes. This review will serve as a resource for future research seeking to develop new methodologies for a better understanding of the basic mechanism of lignin-enzyme binding during the critical hydrolysis process.

Interaction of lignin dimers with model cell membranes: A quartz crystal microbalance and molecular dynamics simulation study

A study of the interaction between cell membranes and small molecules derived from lignin, a protective phenolic biopolymer found in vascular plants, is crucial for identifying their potential as pharmacological and toxicological agents. In this work, the interactions of model cell membranes [supported 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) lipid bilayers] are compared for three βO4 dimers of coniferyl alcohol (G lignin monomer): guaiacylglycerol guaiacol ester with a hydroxypropenyl (HOC₃H₄-) tail (G-βO4'-G), a truncated GG dimer without HOC₃H₄- (G-βO4'-truncG), and a benzylated GG dimer (benzG-βO4'-G). The uptake of the lignin dimers (per mass of lipid) and the energy dissipation (a measure of bilayer disorder) are higher for benzG-βO4'-G and G-βO4'-truncG than those for G-βO4'-G in the gel-phase DPPC bilayer, as measured using quartz crystal microbalance with dissipation (QCM-D). A similar uptake of G-βO4'-truncG is observed for a fluid-phase bilayer of 1,2-dioleoyl-sn-glycero-3-phosphocholine, suggesting that the effect of the bilayer phase on dimer uptake is minimal. The effects of increasing lignin dimer concentration are examined through an analysis of density profiles, potential of mean force curves, lipid order parameters, and bilayer area compressibilities (disorder) in the lipid bilayers obtained from molecular dynamics simulations. Dimer distributions and potentials of mean force indicate that the penetration into bilayers is higher for benzG-βO4'-G and G-βO4'-truncG than that for G-βO4'-G, consistent with the QCM-D results. Increased lipid tail disorder due to dimer penetration leads to a thinning and softening of the bilayers. Minor differences in the structure of lignin derivatives (such as truncating the hydroxypropenyl tail) have significant impacts on their ability to penetrate lipid bilayers.

Molecular Insight into the Cosolvent Effect on Lignin-Cellulose Adhesion

Understanding and controlling the physical adsorption of lignin compounds on cellulose pulp are key parameters in the successful optimization of organosolv processes. The effect of binary organic-aqueous solvents on the coordination of lignin to cellulose was studied with molecular dynamics simulations, considering ethanol and acetonitrile to be organic cosolvents in aqueous solutions in comparison to their monocomponent counterparts. The structures of the solvation shells around cellulose and lignin and the energetics of lignin-cellulose adhesion indicate a more effective disruption of lignin-cellulose binding by binary solvents. The synergic effect between solvent components is explained by their preferential interactions with lignin-cellulose complexes. In the presence of pure water, long-lasting H-bonds in the lignin-cellulose complex are observed, promoted by the nonfavorable interactions of lignin with water. Ethanol and acetonitrile compete with water and lignin for cellulose oxygen binding sites, causing a nonlinear decrease in the lignin-cellulose interactions with the amount of the organic component. This effect is modulated by the water exclusion from the cellulose solvation shell by the organic solvent component. The amount and rate of water exclusion depend on the type of organic cosolvent and its concentration.

Enzymes to unravel bioproducts architecture

Enzymes are essential and ubiquitous biocatalysts involved in various metabolic pathways and used in many industrial processes. Here, we reframe enzymes not just as biocatalysts transforming bioproducts but also as sensitive probes for exploring the structure and composition of complex bioproducts, like meat tissue, dairy products and plant materials, in both food and non-food bioprocesses. This review details the global strategy and presents the most recent investigations to prepare and use enzymes as relevant probes, with a focus on glycoside-hydrolases involved in plant deconstruction and proteases and lipases involved in food digestion. First, to expand the enzyme repertoire to fit bioproduct complexity, novel enzymes are mined from biodiversity and can be artificially engineered. Enzymes are further characterized by exploring sequence/structure/dynamics/function relationships together with the environmental factors influencing enzyme interactions with their substrates. Then, the most advanced experimental and theoretical approaches developed for exploring bioproducts at various scales (from nanometer to millimeter) using active and inactive enzymes as probes are illustrated. Overall, combining multimodal and multiscale approaches brings a better understanding of native-form or transformed bioproduct architecture and composition, and paves the way to mainstream the use of enzymes as probes.

A Multifunctional Cosolvent Pair Reveals Molecular Principles of Biomass Deconstruction

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.

Effects of the Surface Morphology and Conformations of Lignocellulosic Biomass Biopolymers on Their Nanoscale Interactions with Hydrophobic Self-Assembled Monolayers

The effects of the morphology and conformations of the surface biopolymers present on lignocellulosic biomass as well as their steric hindrance on enzymatic adsorption to biomass surfaces remain elusive. In a step to better understand these effects, nanoscale steric forces between a model surface that represents the hydrophobic residues of a cellulase enzyme and a set of reference lignocellulosic substrates were measured using atomic force microscopy (AFM) in liquid media. The reference substrates investigated were prepared by kraft, sulfite, and organosolv pulping pretreatment methods and varied in their surface lignin, xylan, and acetone extractives' contents. Measured steric forces were quantified through fitting to a model developed to describe polyelectrolytes brushes in terms of a brush thickness and a brush grafting density. Our data indicated that cellulose microfibrils extend from the microfibril matrix leading to a long-range steric repulsion and low attractive forces to the hydrophobic model of the enzyme, suggesting that steric hindering can be a possible mechanism for nonproductive binding of enzymes to cellulose. When the amount of xylan increased in the absence of lignin, steric repulsions between the hydrophobic model of the enzyme, and biomass biopolymers decreased as a result of collapsed cellulose microfibrils and adhesion forces increased. This suggests that leaving a small amount of xylan after biomass pretreatment can help improve enzymatic binding to cellulose. Irrespective of the type of lignin present on biomass, grafting densities increased and brush thicknesses decreased compared to those of lignin-free substrates. When compared to lignin-free substrates, lignin-containing substrates had higher attractive forces and lower steric repulsive forces. In addition, AFM images of the reference substrates in the wet and dry states showed that lignin precipitates on the biomass surface where kraft lignin had the highest particle size leading to a limited accessibility of the enzyme to the cellulose in biomass. When the effects of lignin precipitate size, the adhesion force, and steric forces on nonproductive enzymatic binding were all considered, our results indicate that organosolv pretreatment should be the treatment of choice to minimize enzymatic nonproductive binding to lignin.

Nanoscale Engineering of Designer Cellulosomes

Biocatalysts showcase the upper limit obtainable for high-speed molecular processing and transformation. Efforts to engineer functionality in synthetic nanostructured materials are guided by the increasing knowledge of evolving architectures, which enable controlled molecular motion and precise molecular recognition. The cellulosome is a biological nanomachine, which, as a fundamental component of the plant-digestion machinery from bacterial cells, has a key potential role in the successful development of environmentally-friendly processes to produce biofuels and fine chemicals from the breakdown of biomass waste. Here, the progress toward so-called "designer cellulosomes", which provide an elegant alternative to enzyme cocktails for lignocellulose breakdown, is reviewed. Particular attention is paid to rational design via computational modeling coupled with nanoscale characterization and engineering tools. Remaining challenges and potential routes to industrial application are put forward.

Conformations of Low-Molecular-Weight Lignin Polymers in Water

Low-molecular-weight lignin binds to cellulose during the thermochemical pretreatment of biomass for biofuel production, which prevents the efficient hydrolysis of the cellulose to sugars. The binding properties of lignin are influenced strongly by the conformations it adopts. Here, we use molecular dynamics simulations in aqueous solution to investigate the dependence of the shape of lignin polymers on chain length and temperature. Lignin is found to adopt collapsed conformations in water at 300 and 500 K. However, at 300 K, a discontinuous transition is found in the shape of the polymer as a function of the chain length. Below a critical degree of polymerization, Nc =15, the polymer adopts less spherical conformations than above Nc. The transition disappears at high temperatures (500 K) at which only spherical shapes are adopted. An implication relevant to cellulosic biofuel production is that lignin will self-aggregate even at high pretreatment temperatures.

Investigation of the binding properties of a multi-modular GH45 cellulase using bioinspired model assemblies

BACKGROUND

Enzymes degrading plant biomass polymers are widely used in biotechnological applications. Their efficiency can be limited by non-specific interactions occurring with some chemical motifs. In particular, the lignin component is known to bind enzymes irreversibly. In order to determine interactions of enzymes with their substrates, experiments are usually performed on isolated simple polymers which are not representative of plant cell wall complexity. But when using natural plant substrates, the role of individual chemical and structural features affecting enzyme-binding properties is also difficult to decipher.

RESULTS

We have designed and used lignified model assemblies of plant cell walls as templates to characterize binding properties of multi-modular cellulases. These three-dimensional assemblies are modulated in their composition using the three principal polymers found in secondary plant cell walls (cellulose, hemicellulose, and lignin). Binding properties of enzymes are obtained from the measurement of their mobility that depends on their interactions with the polymers and chemical motifs of the assemblies. The affinity of the multi-modular GH45 cellulase was characterized using a statistical analysis to determine the role played by each assembly polymer. Presence of hemicellulose had much less impact on affinity than cellulose and model lignin. Depending on the number of CBMs appended to the cellulase catalytic core, binding properties toward cellulose and lignin were highly contrasted.

CONCLUSIONS

Model assemblies bring new insights into the molecular determinants that are responsible for interactions between enzymes and substrate without the need of complex analysis. Consequently, we believe that model bioinspired assemblies will provide relevant information for the design and optimization of enzyme cocktails in the context of biorefineries.

Molecular Origin of Strength and Stiffness in Bamboo Fibrils

Bamboo, a fast-growing grass, has a higher strength-to-weight ratio than steel and concrete. The unique properties of bamboo come from the natural composite structure of fibers that consists mainly of cellulose microfibrils in a matrix of intertwined hemicellulose and lignin called lignin-carbohydrate complex (LCC). Here, we have used atomistic simulations to study the mechanical properties of and adhesive interactions between the materials in bamboo fibers. With this aim, we have developed molecular models of lignin, hemicellulose and LCC structures to study the elastic moduli and the adhesion energies between these materials and cellulose microfibril faces. Good agreement was observed between the simulation results and experimental data. It was also shown that the hemicellulose model has stronger mechanical properties than lignin while lignin exhibits greater tendency to adhere to cellulose microfibrils. The study suggests that the abundance of hydrogen bonds in hemicellulose chains is responsible for improving the mechanical behavior of LCC. The strong van der Waals forces between lignin molecules and cellulose microfibril is responsible for higher adhesion energy between LCC and cellulose microfibrils. We also found out that the amorphous regions of cellulose microfibrils are the weakest interfaces in bamboo fibrils. Hence, they determine the fibril strength.

Molecular dynamics simulations of large macromolecular complexes

Connecting dynamics to structural data from diverse experimental sources, molecular dynamics simulations permit the exploration of biological phenomena in unparalleled detail. Advances in simulations are moving the atomic resolution descriptions of biological systems into the million-to-billion atom regime, in which numerous cell functions reside. In this opinion, we review the progress, driven by large-scale molecular dynamics simulations, in the study of viruses, ribosomes, bioenergetic systems, and other diverse applications. These examples highlight the utility of molecular dynamics simulations in the critical task of relating atomic detail to the function of supramolecular complexes, a task that cannot be achieved by smaller-scale simulations or existing experimental approaches alone.

Supramolecular Interactions in Secondary Plant Cell Walls: Effect of Lignin Chemical Composition Revealed with the Molecular Theory of Solvation

Plant biomass recalcitrance, a major obstacle to achieving sustainable production of second generation biofuels, arises mainly from the amorphous cell-wall matrix containing lignin and hemicellulose assembled into a complex supramolecular network that coats the cellulose fibrils. We employed the statistical-mechanical, 3D reference interaction site model with the Kovalenko-Hirata closure approximation (or 3D-RISM-KH molecular theory of solvation) to reveal the supramolecular interactions in this network and provide molecular-level insight into the effective lignin-lignin and lignin-hemicellulose thermodynamic interactions. We found that such interactions are hydrophobic and entropy-driven, and arise from the expelling of water from the mutual interaction surfaces. The molecular origin of these interactions is carbohydrate-π and π-π stacking forces, whose strengths are dependent on the lignin chemical composition. Methoxy substituents in the phenyl groups of lignin promote substantial entropic stabilization of the ligno-hemicellulosic matrix. Our results provide a detailed molecular view of the fundamental interactions within the secondary plant cell walls that lead to recalcitrance.

Mechanism of lignin inhibition of enzymatic biomass deconstruction

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.

Why genetic modification of lignin leads to low-recalcitrance biomass

Molecular dynamics simulations show genetically modified lignin to associate less with hemicellulose than does wild type.

Plant biomass recalcitrance: effect of hemicellulose composition on nanoscale forces that control cell wall strength

Efficient conversion of lignocellulosic biomass to second-generation biofuels and valuable chemicals requires decomposition of resilient plant cell wall structure. Cell wall recalcitrance varies among plant species and even phenotypes, depending on the chemical composition of the noncellulosic matrix. Changing the amount and composition of branches attached to the hemicellulose backbone can significantly alter the cell wall strength and microstructure. We address the effect of hemicellulose composition on primary cell wall assembly forces by using the 3D-RISM-KH molecular theory of solvation, which provides statistical-mechanical sampling and molecular picture of hemicellulose arrangement around cellulose. We show that hemicellulose branches of arabinose, glucuronic acid, and especially glucuronate strengthen the primary cell wall by strongly coordinating to hydrogen bond donor sites on the cellulose surface. We reveal molecular forces maintaining the cell wall structure and provide directions for genetic modulation of plants and pretreatment design to render biomass more amenable to processing.