Understanding Nanocellulose-Water Interactions: Turning a Detriment into an Asset.

Modern technology has enabled the isolation of nanocellulose from plant-based fibers, and the current trend focuses on utilizing nanocellulose in a broad range of sustainable materials applications. Water is generally seen as a detrimental component when in contact with nanocellulose-based materials, just like it is harmful for traditional cellulosic materials such as paper or cardboard. However, water is an integral component in plants, and many applications of nanocellulose already accept the presence of water or make use of it. This review gives a comprehensive account of nanocellulose-water interactions and their repercussions in all key areas of contemporary research: fundamental physical chemistry, chemical modification of nanocellulose, materials applications, and analytical methods to map the water interactions and the effect of water on a nanocellulose matrix.

Entropy Drives Interpolymer Association in Water: Insights into Molecular Mechanisms.

Interpolymer association in aqueous solutions is essential for many industrial processes, new materials design, and the biochemistry of life. However, our understanding of the association mechanism is limited. Classical theories do not provide molecular details, creating a need for detailed mechanistic insights. This work consolidates previous literature with complementary isothermal titration calorimetry (ITC) measurements and molecular dynamics (MD) simulations to investigate molecular mechanisms to provide such insights. The large body of ITC data shows that intermolecular bonds, such as ionic or hydrogen bonds, cannot drive association. Instead, polymer association is entropy-driven due to the reorganization of water and ions. We propose a unifying entropy-driven association mechanism by generalizing previously suggested polyion association principles to include nonionic polymers, here termed polydipoles. In this mechanism, complementary charge densities of the polymers are the common denominators of association, for both polyions and polydipoles. The association of the polymers results mainly from two processes: charge exchange and amphiphilic association. MD simulations indicate that the amphiphilic assembly alone is enough for the initial association. Our proposed mechanism is a step toward a molecular understanding of the formation of complexes between synthetic and biological polymers under ambient or biological conditions.

Capillary forces exerted by a water bridge on cellulose nanocrystals: the effect of an external electric field.

Capillary forces play an important role during the dewatering and drying of nanocellulosic materials. Traditional moisture removal techniques, such as heating, have been proved to be deterimental to the properties of these materials and hence, there is a need to develop novel dewatering techniques without affecting the desired properties of materials. It is, therefore, important to explore novel methods for dewatering these high-added-value materials without negatively influencing their properties. In this context, we explore the effect of electric field on the capillary forces developed by a liquid-water bridge between two cellulosic surfaces, which may be formed during the water removal process following its displacement from the interfibrillar spaces. All-atom molecular dynamics (MD) simulations have been used to study the influence of an externally applied electric field on the capillary force exerted by a water bridge. Our results suggest that the equilibrium contact angle of water and the capillary force exerted by the water bridge between two nanocellulosic surfaces depend on the magnitude and direction of the externally applied electric fields. Hence, an external electric field can be applied to manipulate the capillary forces between two particles. The close agreement between the capillary forces measured through MD simulations and those calculated through classical equations indicates that, within the range of the electric field applied in this study, Young-Laplace equations can be safely employed to predict the capillary forces between two particles. The present study provides insights into the use of electric fields for drying of nanocellulosic materials.

Quantifying the Contribution of the Dispersion Interaction and Hydrogen Bonding to the Anisotropic Elastic Properties of Chitin and Chitosan.

The elastic tensors of chitin and chitosan allomorphs were calculated using density functional theory (DFT) with and without the dispersion correction and compared with experimental values. The longitudinal Young's moduli were 114.9 or 126.9 GPa for α-chitin depending on the hydrogen bond pattern: 129.0 GPa for ß-chitin and 191.5 GPa for chitosan. Furthermore, the moduli were found to vary between 17.0 and 52.8 GPa in the transverse directions and between 2.2 and 15.2 GPa in shear. Switching off the dispersion correction led to a decrease in modulus by up to 63%, depending on the direction. The transverse Young's moduli of α-chitin strongly depended on the hydroxylmethyl group conformation coupled with the dispersion correction, suggesting a synergy between hydrogen bonding and dispersion interactions. The calculated longitudinal Young's moduli were, in general, higher than experimental values obtained in static conditions, and the Poisson's ratios were lower than experimental values obtained in static conditions.

"Like Recycles Like": Selective Ring-Closing Depolymerization of Poly(L-Lactic Acid) to L-Lactide.

Chemical recycling of poly(L-lactic acid) to the cyclic monomer L-lactide is hampered by low selectivity and by epimerization and elimination reactions, impeding its use on a large scale. The high number of side reactions originates from the high ceiling temperature (Tc ) of L-lactide, which necessitates high temperatures or multistep reactions to achieve recycling to L-lactide. To circumvent this issue, we utilized the impact of solvent interactions on the monomer-polymer equilibrium to decrease the Tc of L-lactide. Analyzing the observed Tc in different solvents in relation to their Hildebrand solubility parameter revealed a "like recycles like" relationship. The decreased Tc , obtained by selecting solvents that interact strongly with the monomer (dimethyl formamide or the green solvent γ-valerolactone), allowed chemical recycling of high-molecular-weight poly(L-lactic acid) directly to L-lactide, within 1-4 h at 140 °C, with >95 % conversion and 98-99 % selectivity. Recycled L-lactide was isolated and repolymerized with high control over molecular weight and dispersity, closing the polymer loop.

Regular Motifs in Xylan Modulate Molecular Flexibility and Interactions with Cellulose Surfaces.

Xylan is tightly associated with cellulose and lignin in secondary plant cell walls, contributing to its rigidity and structural integrity in vascular plants. However, the molecular features and the nanoscale forces that control the interactions among cellulose microfibrils, hemicelluloses, and lignin are still not well understood. Here, we combine comprehensive mass spectrometric glycan sequencing and molecular dynamics simulations to elucidate the substitution pattern in softwood xylans and to investigate the effect of distinct intramolecular motifs on xylan conformation and on the interaction with cellulose surfaces in Norway spruce (Picea abies). We confirm the presence of motifs with evenly spaced glycosyl decorations on the xylan backbone, together with minor motifs with consecutive glucuronation. These domains are differently enriched in xylan fractions extracted by alkali and subcritical water, which indicates their preferential positioning in the secondary plant cell wall ultrastructure. The flexibility of the 3-fold screw conformation of xylan in solution is enhanced by the presence of arabinofuranosyl decorations. Additionally, molecular dynamic simulations suggest that the glycosyl substitutions in xylan are not only sterically tolerated by the cellulose surfaces but that they increase the affinity for cellulose and favor the stabilization of the 2-fold screw conformation. This effect is more significant for the hydrophobic surface compared with the hydrophilic ones, which demonstrates the importance of nonpolar driving forces on the structural integrity of secondary plant cell walls. These novel molecular insights contribute to an improved understanding of the supramolecular architecture of plant secondary cell walls and have fundamental implications for overcoming lignocellulose recalcitrance and for the design of advanced wood-based materials.

Toward Improved Understanding of the Interactions between Poorly Soluble Drugs and Cellulose Nanofibers.

Cellulose nanofibers (CNFs) have interesting physicochemical and colloidal properties that have been recently exploited in novel drug-delivery systems for tailored release of poorly soluble drugs. The morphology and release kinetics of such drug-delivery systems heavily relied on the drug-CNF interactions; however, in-depth understanding of the interactions was lacking. Herein, the interactions between a poorly soluble model drug molecule, furosemide, and cationic cellulose nanofibers with two different degrees of substitution are studied by sorption experiments, Fourier transform infrared spectroscopy, and molecular dynamics (MD) simulation. Both MD simulations and experimental results confirmed the spontaneous sorption of drug onto CNF. Simulations further showed that adsorption occurred by the flat aryl ring of furosemide. The spontaneous sorption was commensurate with large entropy gains as a result of release of surface-bound water. Association between furosemide molecules furthermore enabled surface precipitation as indicated by both simulations and experiments. Finally, sorption was also found not to be driven by charge neutralization, between positive CNF surface charges and the furosemide negative charge, so that surface area is the single most important parameter determining the amount of sorbed drug. An optimized CNF-furosemide drug-delivery vehicle thus needs to have a maximized specific surface area irrespective of the surface charge with which it is achieved. The findings also provide important insights into the design principles of CNF-based filters suitable for removal of poorly soluble drugs from wastewater.

Hydration-Dependent Dynamical Modes in Xyloglucan from Molecular Dynamics Simulation of 13C NMR Relaxation Times and Their Distributions.

Macromolecular dynamics in biological systems, which play a crucial role for biomolecular function and activity at ambient temperature, depend strongly on moisture content. Yet, a generally accepted quantitative model of hydration-dependent phenomena based on local relaxation and diffusive dynamics of both polymer and its adsorbed water is still missing. In this work, atomistic-scale spatial distributions of motional modes are calculated using molecular dynamics simulations of hydrated xyloglucan (XG). These are shown to reproduce experimental hydration-dependent 13C NMR longitudinal relaxation times ( T1) at room temperature, and relevant features of their broad distributions, which are indicative of locally heterogeneous polymer reorientational dynamics. At low hydration, the self-diffusion behavior of water shows that water molecules are confined to particular locations in the randomly aggregated XG network while the average polymer segmental mobility remains low. Upon increasing water content, the hydration network becomes mobile and fully accessible for individual water molecules, and the motion of hydrated XG segments becomes faster. Yet, the polymer network retains a heterogeneous gel-like structure even at the highest level of hydration. We show that the observed distribution of relaxations times arises from the spatial heterogeneity of chain mobility that in turn is a result of heterogeneous distribution of water-chain and chain-chain interactions. Our findings contribute to the picture of hydration-dependent dynamics in other macromolecules such as proteins, DNA, and synthetic polymers, and hold important implications for the mechanical properties of polysaccharide matrixes in plants and plant-based materials.

A molecular dynamics study of the effect of glycosidic linkage type in the hemicellulose backbone on the molecular chain flexibility.

The macromolecular conformation of the constituent polysaccharides in lignocellulosic biomass influences their supramolecular interactions, and therefore their function in plants and their performance in technical products. The flexibility of glycosidic linkages from the backbone of hemicelluloses was studied by evaluating the conformational freedom of the φ and ψ dihedral angles using molecular dynamic simulations, additionally selected molecules were correlated with experimental data by nuclear magnetic resonance spectroscopy. Three types of ß-(1→4) glycosidic linkages involving the monosaccharides (Glcp, Xylp and Manp) present in the backbone of hemicelluloses were defined. Different di- and tetrasaccharides with combinations of such sugar monomers from hemicelluloses were simulated, and free energy maps of the φ - ψ space and hydrogen-bonding patterns were obtained. The glycosidic linkage between Glc-Glc or Glc-Man (C-type) was the stiffest with mainly one probable conformation; the linkage from Man-Man or Man-Glc (M-type) was similar but with an increased probability for an alternative conformation making it more flexible, and the linkage between two Xyl-units (X-type) was the most flexible with two almost equally populated conformations. Glycosidic linkages of the same type showed essentially the same conformational space in both disaccharides and in the central region of tetrasaccharides. Different probabilities of glycosidic linkage conformations in the backbone of hemicelluloses can be directly estimated from the free energy maps, which to a large degree affect the overall macromolecular conformations of these polymers. The information gained contributes to an increased understanding of the function of hemicelluloses both in the cell wall and in technical products.

Molecular adhesion at clay nanocomposite interfaces depends on counterion hydration-molecular dynamics simulation of montmorillonite/xyloglucan.

Nacre-mimetic clay/polymer nanocomposites with clay platelet orientation parallel to the film surface show interesting gas barrier and mechanical properties. In moist conditions, interfacial adhesion is lowered and mechanical properties are reduced. Molecular dynamic simulations (MD) have been performed to investigate the effects of counterions on molecular adhesion at montmorillonite clay (Mnt)-xyloglucan (XG) interfaces. We focus on the role of monovalent cations K(+), Na(+), and Li(+) and the divalent cation Ca(2+) for mediating and stabilizing the Mnt/XG complex formation. The conformation of adsorbed XG is strongly influenced by the choice of counterion and so is the simulated work of adhesion. Free energy profiles that are used to estimate molecular adhesion show stronger interaction between XG and clay in the monovalent cation system than in divalent cation system, following a decreasing order of K-Mnt, Na-Mnt, Li-Mnt, and Ca-Mnt. The Mnt clay hydrates differently in the presence of different counterions, leading to a chemical potential of water that is highest in the case of K-Mnt, followed by Na-Mnt and Li-Mnt, and lowest in the case of Ca-Mnt. This means that water is most easily displaced from the interface in the case of K-Mnt, which contributes to the relatively high work of adhesion. In all systems, the penalty of replacing polymer with water at the interface gives a positive contribution to the work of adhesion of between 19 and 35%. Our work confirms the important role of counterions in mediating the adsorption of biopolymer XG to Mnt clays and predicts potassium or sodium as the best choice of counterions for a Mnt-based biocomposite design.

Molecular dynamics simulations and NMR spectroscopy studies of trehalose-lipid bilayer systems.

The disaccharide trehalose (TRH) strongly affects the physical properties of lipid bilayers. We investigate interactions between lipid membranes formed by 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and TRH using NMR spectroscopy and molecular dynamics (MD) computer simulations. We compare dipolar couplings derived from DMPC/TRH trajectories with those determined (i) experimentally in TRH using conventional high-resolution NMR in a weakly ordered solvent (bicelles), and (ii) by solid-state NMR in multilamellar vesicles (MLV) formed by DMPC. Analysis of the experimental and MD-derived couplings in DMPC indicated that the force field used in the simulations reasonably well describes the experimental results with the exception for the glycerol fragment that exhibits significant deviations. The signs of dipolar couplings, not available from the experiments on highly ordered systems, were determined from the trajectory analysis. The crucial step in the analysis of residual dipolar couplings (RDCs) in TRH determined in a bicelle-environment was access to the conformational distributions derived from the MD trajectory. Furthermore, the conformational behavior of TRH, investigated by J-couplings, in the ordered and isotropic phases is essentially identical, indicating that the general assumptions in the analyses of RDCs are well founded.

Acetylation of Nanocellulose: Miscibility and Reinforcement Mechanisms in Polymer Nanocomposites.

The improvement of properties in nanocomposites obtained by topochemical surface modification, e.g., acetylation, of the nanoparticles is often ascribed to improved compatibility between the nanoparticle and the matrix. It is not always clear however what is intended: specific interactions at the interface leading to increased adhesion or the miscibility between the nanoparticle and the polymer. In this work, it is demonstrated that acetylation of cellulose nanocrystals greatly improves mechanical properties of their nanocomposites with polycaprolactone. In addition, molecular dynamics simulations with a combination of potential of mean force calculations and computational alchemy are employed to analyze the surface energies between the two components. The work of adhesion between the two phases decreases with acetylation. It is discussed how acetylation can still contribute to the miscibility, which leads to a stricter use of the concept of compatibility. The integrated experimental-modeling toolbox used has wide applicability for assessing changes in the miscibility of polymer nanocomposites.

An ultrastrong nanofibrillar biomaterial: the strength of single cellulose nanofibrils revealed via sonication-induced fragmentation.

We report the mechanical strength of native cellulose nanofibrils. Native cellulose nanofibrils, purified from wood and sea tunicate, were fully dispersed in water via a topochemical modification of cellulose nanofibrils using 2,2,6,6-tetramethylpiperidinyl-1-oxyl (TEMPO) as a catalyst. The strength of individual nanofibrils was estimated based on a model for the sonication-induced fragmentation of filamentous nanostructures. The resulting strength parameters were then analyzed based on fracture statistics. The mean strength of the wood cellulose nanofibrils ranged from 1.6 to 3 GPa, depending on the method used to measure the nanofibril width. The highly crystalline, thick tunicate cellulose nanofibrils exhibited higher mean strength of 3-6 GPa. The strength values estimated for the cellulose nanofibrils in the present study are comparable with those of commercially available multiwalled carbon nanotubes.

The thermodynamics of enhanced dope stability of cellulose solution in NaOH solution by urea.

The addition of urea in pre-cooled alkali aqueous solution is known to improve the dope stability of cellulose solution. However, its thermodynamic mechanism at a molecular level is not fully understood yet. By using molecular dynamics simulation of an aqueous NaOH/urea/cellulose system using an empirical force field, we found that urea was concentrated in the first solvation shell of the cellulose chain stabilized mainly by dispersion interaction. When adding a glucan chain into the solution, the total solvent entropy reduction is smaller if urea is present. Each urea molecule expelled an average of 2.3 water molecules away from the cellulose surface, releasing water entropy that over-compensates the entropy loss of urea and thus maximizing the total entropy. Scaling the Lennard-Jones parameter and atomistic partial charge of urea revealed that direct urea/cellulose interaction was also driven by dispersion energy. The mixing of urea solution and cellulose solution in the presence or absence of NaOH are both exothermic even after correcting for the contribution from dilution.

Water as an Intrinsic Structural Element in Cellulose Fibril Aggregates.

While strong water association with cellulose in plant cell walls and man-made materials is well-established, its molecular scale aspects are not fully understood. The thermodynamic consequences of having water molecules located at the microfibril-microfibril interfaces in cellulose fibril aggregates are therefore analyzed by molecular dynamics simulations. We find that a thin layer of water molecules at those interfaces can be in a state of thermal equilibrium with water surrounding the fibril aggregates because such an arrangement lowers the free energy of the total system. The main reason is enthalpic: water at the microfibril-microfibril interfaces enables the cellulose surface hydroxyls to experience a more favorable electrostatic environment. This enthalpic gain overcomes the entropic penalty from strong immobilization of water molecules. Hence, those particular water molecules stabilize the cellulose fibril aggregates, akin to the role of water in some proteins. Structural and functional hypotheses related to this finding are presented.

Timescales for convergence in all-atom molecular dynamics simulations of hydrated amorphous xylan.

Atomistic molecular dynamics simulation is an important complement to experimental techniques for investigating properties of hydrated carbohydrate polymers at the molecular scale. A critical problem is to determine whether or not a simulation has converged to thermal equilibrium before data collection can begin. In this work, simulations of xylan oligomers starting from random configurations at different levels of hydration are performed. The simulations show clear evidence of phase separation into water-rich and polymer-rich phases at higher hydration, in spite of standard indicators of equilibrium, such as density and energy, remaining constant. Using instead a set of parameters that are coupled to the structural and dynamical heterogeneity of the system, it is shown that simulation times on the order of one microsecond are needed to reach an equilibrated state. Moreover, qualitative similarities in the temporal evolution of these parameters suggest significant interplay between the structure and both polymer and water dynamics.

Xylan adsorption on cellulose: Preferred alignment and local surface immobilizing effect.

Interaction between xylan and cellulose microfibrils is required to maintain the integrity of secondary cell walls. However, the mechanisms governing their assembly and the effects on cellulose surface polymers are not fully clear. Here, molecular dynamics simulations are used to study xylan adsorption onto hydrated cellulose fibrils. Based on multiple spontaneous adsorption simulations it is shown that an antiparallel orientation is thermodynamically preferred over a parallel one, and that adsorption is accompanied by the formation of regular but orientation-dependent hydrogen bond patterns. Furthermore, xylan adsorption restricts the local dynamics of the adjacent glucose residues in the surface layer to a level of the crystalline core, which is manifested as a three-fold increase in their 13C NMR T1 relaxation time. These results suggest that xylan forms a rigid and ordered layer around the cellulose fibril that functions as a transition phase to more flexible and disordered polysaccharide and lignin domains.

Entropy drives the adsorption of xyloglucan to cellulose surfaces - A molecular dynamics study.

The adsorption of nonionic polymers to cellulose is of large importance both in the plant cell wall during synthesis and for the development of sustainable materials from wood. Here, the thermodynamics of adsorption of the polysaccharide xyloglucan (XG) to both native and chemically modified cellulose with carboxyl groups was investigated using molecular dynamics simulations. The free energy of adsorption was calculated as the potential of mean force between an XG oligomer and model cellulose surfaces in a range of temperatures from 298 K to 360 K. It was found that the adsorption near room temperature is an endothermic process dominated by the entropy of released interfacial water molecules. This was corroborated by quantitative assessment of the absolute entropy per water molecule both at the interface and in the bulk. In the case of native cellulose, the adsorption became exothermic at higher temperatures, while the relatively strong interactions between water and the charged groups of the oxidized cellulose impede such a transition. The results also indicate that the extraction of strongly associated hemicelluloses would be facilitated by low temperature.

O-Methylation in Carbohydrates: An NMR and MD Simulation Study with Application to Methylcellulose.

Methylated carbohydrates are important from both biological and technical perspectives. Specifically, methylcellulose is an interesting cellulose derivative that has applications in foods, materials, cosmetics, and many other fields. While the molecular dynamics simulation technique has the potential for both advancing the fundamental understanding of this polymer and aiding in the development of specific applications, a general drawback is the lack of experimentally validated interaction potentials for the methylated moieties. In the present study, simulations using the GROMOS 56 carbohydrate force field are compared to NMR spin-spin coupling constants related to the conformation of the exocyclic torsion angle ω in d-glucopyranose and derivatives containing a 6-O-methyl substituent and a 13C-isotopologue thereof. A 3JCC Karplus-type relationship is proposed for the C5-C6-O6-CMe torsion angle. Moreover, solvation free energies are compared to experimental data for small model compounds. Alkylation in the form of 6-O-methylation affects exocyclic torsion only marginally. Computed solvation free energies between nonmethylated and methylated molecules were internally consistent, which validates the application of these interaction potentials for more specialized purposes.

Wettability of cellulose surfaces under the influence of an external electric field.

HYPOTHESIS: Interfacial tensions play an important role in dewatering of hydrophilic materials like nanofibrillated cellulose, and are affected by the molecular organization of water at the interface. Application of an electric field influences the orientation of water molecules along the field direction. Hence, it should be possible to alter the interfacial free energies to tune the wettability of cellulose surface through application of an external electric field thus, aiding the dewatering process. SIMULATIONS: Molecular dynamics simulations of cellulose surface in contact with water under the influence of an external electric field have been conducted with GLYCAM-06 forcefield. The effect of variation in electric field intensity and directions on the spreading coefficient has been addressed via orientational preference of water molecules and interfacial free energy analyses. FINDINGS: The application of electric field influences the interfacial free energy difference at the cellulose-water interface. The spreading coefficient increases with the electric field directed parallel to the cellulose-water interface while it decreases in the perpendicular electric field. Variation in interfacial free energies seems to explain the change in contact angle adequately in presence of an electric field. The wettability of cellulose surface can be tuned by the application of an external electric field.