Computer simulation studies of microcrystalline cellulose Iβ

Unique aspects of the structure and dynamics of elementary Iβ cellulose microfibrils revealed by computational simulations

The question of how many chains an elementary cellulose microfibril contains is critical to understanding the molecular mechanism(s) of cellulose biosynthesis and regulation. Given the hexagonal nature of the cellulose synthase rosette, it is assumed that the number of chains must be a multiple of six. We present molecular dynamics simulations on three different models of Iβ cellulose microfibrils, 18, 24, and 36 chains, to investigate their structure and dynamics in a hydrated environment. The 36-chain model stays in a conformational space that is very similar to the initial crystalline phase, while the 18- and 24-chain models sample a conformational space different from the crystalline structure yet similar to conformations observed in recent high-temperature molecular dynamics simulations. Major differences in the conformations sampled between the different models result from changes to the tilt of chains in different layers, specifically a second stage of tilt, increased rotation about the O2-C2 dihedral, and a greater sampling of non-TG exocyclic conformations, particularly the GG conformation in center layers and GT conformation in solvent-exposed exocyclic groups. With a reinterpretation of nuclear magnetic resonance data, specifically for contributions made to the C6 peak, data from the simulations suggest that the 18- and 24-chain structures are more viable models for an elementary cellulose microfibril, which also correlates with recent scattering and diffraction experimental data. These data inform biochemical and molecular studies that must explain how a six-particle cellulose synthase complex rosette synthesizes microfibrils likely comprised of either 18 or 24 chains.

Effects of lytic polysaccharide monooxygenase oxidation on cellulose structure and binding of oxidized cellulose oligomers to cellulases

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.

Probing the molecular architecture of Arabidopsis thaliana secondary cell walls using two- and three-dimensional (13)C solid state nuclear magnetic resonance spectroscopy

The plant secondary cell wall is a thickened polysaccharide and phenolic structure, providing mechanical strength to cells, particularly in woody tissues. It is the main feedstock for the developing bioenergy and green chemistry industries. Despite the role that molecular architecture (the arrangement of biopolymers relative to each other, and their conformations) plays in dictating biomass properties, such as recalcitrance to breakdown, it is poorly understood. Here, unprocessed dry (13)C-labeled stems from the model plant Arabidopsis thaliana were analyzed by a variety of (13)C solid state magic angle spinning nuclear magnetic resonance methods, such as one-dimensional cross-polarization and direct polarization, two-dimensional refocused INADEQUATE, RFDR, PDSD, and three-dimensional DARR, demonstrating their viability for the study of native polymer arrangements in intact secondary cell walls. All carbon sites of the two main glucose environments in cellulose (previously assigned to microfibril surface and interior residues) are clearly resolved, as are carbon sites of the other major components of the secondary cell wall: xylan and lignin. The xylan carbon 4 chemical shift is markedly different from that reported previously for solution or primary cell wall xylan, indicating significant changes in the helical conformation in these dried stems. Furthermore, the shift span indicates that xylan adopts a wide range of conformations in this material, with very little in the 31 conformation typical of xylan in solution. Additionally, spatial connections of noncarbohydrate species were observed with both cellulose peaks conventionally assigned as "surface" and as "interior" cellulose environments, raising questions about the origin of these two cellulose signals.

Effect of Moisture on Physicochemical Properties of Cellulose Studied by Molecular Dynamics

In this work, the physicochemical property and its Effect of hydrous amorphous cellulose has been studied using molecular dynamics. Both intramolecular and intermolecular hydrogen bonds in cellulose molecules decrease with increasing water content, directly leading to the decline of cellulose cohesive energy density, solubility parameters, and mechanical parameters. High water content in amorphous cellulose gives bigger interchain distance of cellulose molecules, indicating that the intermolecular interaction of cellulose molecules is weakened greatly by water.

Bubble-surface interactions with graphite in the presence of adsorbed carboxymethylcellulose

The adsorption of carboxymethylcellulose (CMC), and the subsequent effect on bubble-surface interactions, has been studied for a graphite surface. CMC adsorbs on highly oriented pyrolytic graphite (HOPG) in specific patterns: when adsorbed from a solution of low concentration it forms stretched, isolated and sparsely distributed chains, while upon adsorption from a solution of higher concentration, it forms an interconnected network of multilayer features. The amount and topography of the adsorbed CMC affect the electrical properties as well as the wettability of the polymer-modified HOPG surface. Adsorption of CMC onto the HOPG surface causes the zeta potential to be more negative and the modified surface becomes more hydrophilic. This increase in both the absolute value of zeta potential and the surface hydrophilicity originates from the carboxymethyl groups of the CMC polymer. The effect of the adsorbed polymer layer on wetting film drainage and bubble-surface/particle attachment was determined using high speed video microscopy to monitor single bubble-surface collision, and single bubble Hallimond tube flotation experiments. The time of wetting film drainage and the time of three-phase contact line spreading gets significantly longer for polymer-modified HOPG surfaces, indicating that the film rupture and three-phase contact line expansion were inhibited by the presence of polymer. The effect of longer drainage times and slower dewetting correlated with reduced flotation recovery. The molecular kinetic (MK) model was used to quantify the effect of the polymer on dewetting dynamics, and showed an increase in the jump frequency for the polymer adsorbed at the higher concentration.

Water structuring over the hydrophobic surface of cellulose

Many important biological solutes possess not only polar and hydrogen-bonding functionalities but also weakly hydrating, or hydrophobic, surfaces. While the aggregation of these hydrophobic surfaces has been shown to play an important role in the aggregation of individual chains of cellulose, it is not known whether the water structuring imposed by these hydrophobic surfaces more closely resembles that associated with small hydrophobic solutes like methane and fats or more closely resembles that associated with extended hydrophobic surfaces like mica or waxy planes. By using molecular dynamics simulations to characterize the water molecule orientations over different regions of the 100 surface of cellulose in contact with water, it was found that the hydrophobic strips of the cellulose crystal are sufficiently narrow that they hydrate like a fatty acid chain, rather than like a more extended surface, suggesting that their aggregation would be dominated by entropy rather than enthalpy.

Replica-exchange molecular dynamics simulations of cellulose solvated in water and in the ionic liquid 1-butyl-3-methylimidazolium chloride

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.

The impact of alkali pretreatment and post-pretreatment conditioning on the surface properties of rice straw affecting cellulose accessibility to cellulases

Rice straw was pretreated with sodium hydroxide and subsequently conditioned to reduce the pH to 5-6 by either: (1) extensive water washing or (2) acidification with hydrochloric acid then water washing. Alkali pretreatment improved the enzymatic digestibility of rice straw by increasing the cellulose accessibility to cellulases. However, acidification after pretreatment reversed the gains in cellulose accessibility to cellulases and enzymatic digestibility due to precipitation of solubilized compounds. Surface composition analyses by ToF-SIMS confirmed a reduction in surface lignin by pretreatment and water washing, and suggested that acidification precipitated a chemically modified form of lignin on the surfaces of rice straw. The spin-spin relaxation times (T2) of the samples indicated increased porosity in alkali pretreated rice straw. The acidified pretreated rice straw had reduced amounts of water in the longer T2 proton pools associated with water in the pores of the biomass likely due to back-filling by the precipitated components.

Water-polysaccharide interactions in the primary cell wall of Arabidopsis thaliana from polarization transfer solid-state NMR

Polysaccharide-rich plant cell walls are hydrated under functional conditions, but the molecular interactions between water and polysaccharides in the wall have not been investigated. In this work, we employ polarization transfer solid-state NMR techniques to study the hydration of primary-wall polysaccharides of the model plant, Arabidopsis thaliana. By transferring water (1)H polarization to polysaccharides through distance- and mobility-dependent (1)H-(1)H dipolar couplings and detecting it through polysaccharide (13)C signals, we obtain information about water proximity to cellulose, hemicellulose, and pectins as well as water mobility. Both intact and partially extracted cell wall samples are studied. Our results show that water-pectin polarization transfer is much faster than water-cellulose polarization transfer in all samples, but the extent of extraction has a profound impact on the water-polysaccharide spin diffusion. Removal of calcium ions and the consequent extraction of homogalacturonan (HG) significantly slowed down spin diffusion, while further extraction of matrix polysaccharides restored the spin diffusion rate. These trends are observed in cell walls with similar water content, thus they reflect inherent differences in the mobility and spatial distribution of water. Combined with quantitative analysis of the polysaccharide contents, our results indicate that calcium ions and HG gelation increase the amount of bound water, which facilitates spin diffusion, while calcium removal disrupts the gel and gives rise to highly dynamic water, which slows down spin diffusion. The recovery of spin diffusion rates after more extensive extraction is attributed to increased water-exposed surface areas of the polysaccharides. Water-pectin spin diffusion precedes water-cellulose spin diffusion, lending support to the single-network model of plant primary walls in which a substantial fraction of the cellulose surface is surrounded by pectins.

Effect of microfibril twisting on theoretical powder diffraction patterns of cellulose Iβ

Previous studies of calculated diffraction patterns for cellulose crystallites suggest that distortions that arise once models have been subjected to MD simulation are the result of both microfibril twisting and changes in unit cell dimensions induced by the empirical force field; to date, it has not been possible to separate the individual contributions of these effects. To provide a better understanding of how twisting manifests in diffraction data, the present study demonstrates a method for generating twisted and linear cellulose structures that can be compared without the bias of dimensional changes, allowing assessment of the impact of twisting alone. Analysis of unit cell dimensions, microfibril volume, hydrogen bond patterns, glycosidic torsion angles, and hydroxymethyl group orientations confirmed that the twisted and linear structures collected with this method were internally consistent, and theoretical powder diffraction patterns for the two were shown to be effectively indistinguishable. These results indicate that differences between calculated patterns for the crystal coordinates and twisted structures from MD simulation can result entirely from changes in unit cell dimensions, and not from microfibril twisting alone. Although powder diffraction patterns for models in the 81-chain size regime were shown to be unaffected by twisting, suggesting that a modest degree of twist is not inconsistent with experimental data, it may be that other diffraction techniques are capable of detecting this structural difference. Until such time as definitive experimental evidence comes to light, the results of this study suggest that both twisted and linear microfibrils may represent an appropriate model for cellulose Iβ.

Supramolecular stacking motifs in the solid state of amide and triazole derivatives of cellobiose

1-Acetamido-1-deoxy-(4-O-β-d-glucopyranosyl-β-d-glucopyranose) (5) and 1-deoxy-1-(4-phenyl-1,2,3-triazolyl)-(4-O-β-d-glucopyranosyl-β-d-glucopyranose) (7) were synthesised from 1-azido-1-deoxy-(4-O-β-d-glucopyranosyl-β-d-glucopyranose) (2) and crystallised as dihydrates. Crystal structural analysis of 5·2H2O displayed an acetamide C(4) chain and stacked cellobiose residues. The structure of 7·2H2O featured π-π stacking and stacking of the cellobiose residues.

Multiscale deconstruction of molecular architecture in corn stover

Lignocellulosic composite in corn stover is a candidate biofuel feedstock of substantial abundance and sustainability. Its utilization is hampered by resistance of constituent cellulose fibrils to deconstruction. Here we use multi-scale studies of pretreated corn stover to elucidate the molecular mechanism of deconstruction and investigate the basis of recalcitrance. Dilute acid pretreatment has modest impact on fibrillar bundles at 0.1 micron length scales while leading to significant disorientation of individual fibrils. It disintegrates many fibrils into monomeric cellulose chains or small side-by-side aggregates. Residual crystalline fibrils lose amorphous surface material, change twist and where still cross-linked, coil around one another. Yields from enzymatic digestion are largely due to hydrolysis of individual cellulose chains and fragments generated during pretreatments. Fibrils that remain intact after pretreatment display substantial resistance to enzymatic digestion. Optimization of yield will require strategies that maximize generation of fragments and minimize preservation of intact cellulosic fibrils.

Correlating the ability of lignocellulosic polymers to constrain water with the potential to inhibit cellulose saccharification

BACKGROUND

Studies in bioconversions have continuously sought the development of processing strategies to overcome the "close physical association" between plant cell wall polymers thought to significantly contribute to biomass recalcitrance [Adv Space Res 18:251-265, 1996],[ Science 315:804-807, 2007]. To a lesser extent, studies have sought to understand biophysical factors responsible for the resistance of lignocelluloses to enzymatic degradation. Provided here are data supporting our hypothesis that the inhibitory potential of different cell wall polymers towards enzymatic cellulose hydrolysis is related to how much these polymers constrain the water surrounding them. We believe the entropy-reducing constraint imparted to polymer associated water plays a negative role by increasing the probability of detrimental interactions such as junction zone formation and the non-productive binding of enzymes.

RESULTS

Selected commercial lignocellulose-derived polymers, including hemicelluloses, pectins, and lignin, showed varied potential to inhibit 24-h cellulose conversion by a mix of purified cellobiohydrolase I and β-glucosidase. At low dry matter loadings (0.5% w/w), insoluble hemicelluloses were most inhibitory (reducing conversion relative to cellulose-only controls by about 80%) followed by soluble xyloglucan and wheat arabinoxylan (reductions of about 70% and 55%, respectively), while the lignin and pectins tested were the least inhibitory (approximately 20% reduction). Low field nuclear magnetic resonance (LF-NMR) relaxometry used to observe water-related proton relaxation in saturated polymer suspensions (10% dry solids, w/w) showed spin-spin, T2, relaxation time curves generally approached zero faster for the most inhibitory polymer preparations. The manner of this decline varied between polymers, indicating different biophysical aspects may differentially contribute to overall water constraint in each case. To better compare the LF-NMR data to inhibitory potential, T2 values from monocomponent exponential fits of relaxation curves were used as a measure of overall water constraint. These values generally correlated faster relaxation times (greater water constraint) with greater inhibition of the model cellulase system by the polymers.

CONCLUSIONS

The presented correlation of cellulase inhibition and proton relaxation data provides support for our water constraint-biomass recalcitrance hypothesis. Deeper investigation into polymer-cellulose-cellulase interactions should help elucidate the types of interactions that may be propagating this correlation. If these observations can be verified to be more than correlative, the hypothesis and data presented suggest that a focus on water-polymer interactions and ways to alter them may help resolve key biological lignocellulose processing bottlenecks.

Levoglucosan formation from crystalline cellulose: importance of a hydrogen bonding network in the reaction

Levoglucosan (1,6-anhydro-β-D-glucopyranose) formation by the thermal degradation of native cellulose was investigated by MP4(SDQ)//DFT(B3LYP) and DFT(M06-2X)//DFT(B3LYP) level computations. The computational results of dimer models lead to the conclusion that the degradation occurs by a concerted mechanism similar to the degradation of methyl β-D-glucoside reported in our previous study. One-chain models of glucose hexamer, in which the interchain hydrogen bonds of real cellulose crystals are absent, do not exhibit the correct reaction behavior of levoglucosan formation; for instance, the activation enthalpy (Ea =≈38 kcal mol(-1) ) is considerably underestimated compared to the experimental value (48-60 kcal mol(-1) ). This problem is solved with the use of two-chain models that contain interchain hydrogen bonds. The theoretical study of this model clearly shows that the degradation of the internal glucosyl residue leads to the formation of a levoglucosan precursor at the chain end and levoglucosan is selectively formed from this levoglucosan end. The calculated Ea (56-62 kcal mol(-1) ) agrees well with the experimental value. The computational results of three-chain models indicate that this degradation occurs selectively on the crystalline surface. All these computational results provide a comprehensive understanding of several experimental facts, the mechanisms of which have not yet been elucidated.

3D electron tomography of pretreated biomass informs atomic modeling of cellulose microfibrils

Fundamental insights into the macromolecular architecture of plant cell walls will elucidate new structure-property relationships and facilitate optimization of catalytic processes that produce fuels and chemicals from biomass. Here we introduce computational methodology to extract nanoscale geometry of cellulose microfibrils within thermochemically treated biomass directly from electron tomographic data sets. We quantitatively compare the cell wall nanostructure in corn stover following two leading pretreatment strategies: dilute acid with iron sulfate co-catalyst and ammonia fiber expansion (AFEX). Computational analysis of the tomographic data is used to extract mathematical descriptions for longitudinal axes of cellulose microfibrils from which we calculate their nanoscale curvature. These nanostructural measurements are used to inform the construction of atomistic models that exhibit features of cellulose within real, process-relevant biomass. By computational evaluation of these atomic models, we propose relationships between the crystal structure of cellulose Iβ and the nanoscale geometry of cellulose microfibrils.

CRAFS: a model to analyze two-dimensional X-ray diffraction patterns of plant cellulose

Cellulose from higher plants is a vast renewable resource organized as crystals. Analysis of these crystals by X-ray diffraction poses very specific challenges, including ubiquitous crystallite texture and substantial overlapping of diffraction peaks. In this article, a tailor-made model named Cellulose Rietveld Analysis for Fine Structure (CRAFS) is developed to analyze two-dimensional X-ray diffraction patterns from raw and processed plant cellulose. One-dimensional powder diffractograms are analyzable as a particular case. The CRAFS model considers cellulose Iβ crystal structure, fibrillar crystal shape, paracrystalline peak broadening, pseudo-Voigt peak profiles, harmonic crystallite orientation distribution function and diffraction in fiber geometry. Formulated on the basis of the Rietveld method, CRAFS is presently written in the MATLAB computing language. A set of meaningful coefficients are output from each analyzed pattern. To exemplify model applicability, representative samples are analyzed, bringing some general insights and evidencing the model's potential for systematic parameterization of the fine structure of raw and processed plant celluloses.

Structure of cellulose microfibrils in primary cell walls from collenchyma

In the primary walls of growing plant cells, the glucose polymer cellulose is assembled into long microfibrils a few nanometers in diameter. The rigidity and orientation of these microfibrils control cell expansion; therefore, cellulose synthesis is a key factor in the growth and morphogenesis of plants. Celery (Apium graveolens) collenchyma is a useful model system for the study of primary wall microfibril structure because its microfibrils are oriented with unusual uniformity, facilitating spectroscopic and diffraction experiments. Using a combination of x-ray and neutron scattering methods with vibrational and nuclear magnetic resonance spectroscopy, we show that celery collenchyma microfibrils were 2.9 to 3.0 nm in mean diameter, with a most probable structure containing 24 chains in cross section, arranged in eight hydrogen-bonded sheets of three chains, with extensive disorder in lateral packing, conformation, and hydrogen bonding. A similar 18-chain structure, and 24-chain structures of different shape, fitted the data less well. Conformational disorder was largely restricted to the surface chains, but disorder in chain packing was not. That is, in position and orientation, the surface chains conformed to the disordered lattice constituting the core of each microfibril. There was evidence that adjacent microfibrils were noncovalently aggregated together over part of their length, suggesting that the need to disrupt these aggregates might be a constraining factor in growth and in the hydrolysis of cellulose for biofuel production.

Interactions binding mineral and organic phases in nanocomposites based on bacterial cellulose and calcium phosphates

The interactions responsible for the adhesion of calcium phosphate (CP) nanocrystals and bacterial cellulose (BC) nanofibrils in the composite material obtained by mixing aqueous suspensions of presynthesized CP and BC and the dependence of these interactions on the CP morphology and chemical structure have been elucidated by molecular mechanics calculations of the CP-BC interfacial structures. The interactions between the superficial CP and BC crystal layers have been simulated. Two crystalline CP structures (i.e., hydroxyapatite (HAP) and whitlockite) with two morphologies (plate-shaped and rod-shaped) were considered. Electrostatics has been found to be the major contributor to the adhesion of the CP crystallites and BC nanofibers, and the formation of interfacial hydrogen bonds makes a minor contribution to the interaction energy. It has also been found that, in general, the energy gain resulting from whitlockite-BC binding is greater than that for HAP-BC binding, and the binding of the rod-shaped crystallites of whitlockite with BC is the most profitable. The energy loss and entropy gain upon replacement of the BC-water and CP-water contacts by the BC-CP contacts have been estimated.

Structural basis for substrate targeting and catalysis by fungal polysaccharide monooxygenases

The use of cellulases remains a major cost in the production of renewable fuels and chemicals from lignocellulosic biomass. Fungi secrete copper-dependent polysaccharide monooxygenases (PMOs) that oxidatively cleave crystalline cellulose and improve the effectiveness of cellulases. However, the means by which PMOs recognize and cleave their substrates in the plant cell wall remain unclear. Here, we present structures of Neurospora crassa PMO-2 and PMO-3 at 1.10 and 1.37 Å resolution, respectively. In the structures, dioxygen species are found in the active sites, consistent with the proposed cleavage mechanism. Structural and sequence comparisons between PMOs also reveal that the enzyme substrate-binding surfaces contain highly varied aromatic amino acid and glycosylation positions. The structures reported here provide evidence for a wide range of PMO substrate recognition patterns in the plant cell wall, including binding modes that traverse multiple glucan chains.

Binding preferences, surface attachment, diffusivity, and orientation of a family 1 carbohydrate-binding module on cellulose

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