Comparison of Cellulose Iβ Simulations with Three Carbohydrate Force Fields

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.

Advanced Molecular Dynamics Model for Investigating Biological-Origin Microfibril Structures

Understanding the atomic-scale structure of wood microfibrils is essential for establishing fundamental properties in various wood-based research aspects, including moisture impact, wood modification, and pretreatment. In this study, we employed molecular dynamics simulations to investigate the arrangement of wood polymers, including cellulose, hemicellulose, and lignin, with a primary focus on the composition of softwood, specifically Norway Spruce wood. We assessed the accuracy of our molecular dynamics model by comparing it with available experimental data, such as density, Young's modulus, and glass transition temperature, which ensures the reliability of our approach. A key aspect of our study involved modeling the active sorption site for water interaction with wood polymers. Our findings revealed that the interaction between water and hemicellulose, particularly within the hemicellulose-cellulose interphase, was the most prominent binding site. This observation aligns with prior research in this field, further strengthening the validity of our results.

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.

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.

Cellulose Iβ microfibril interaction with pristine graphene in water: Effects of amphiphilicity by molecular simulation

Graphene-cellulose interactions have considerable potential in the development of new materials. In previous computational work (Biomacromolecules2016, 16, 1771), we predicted that the model 100 hydrophobic surface of cellulose interacted favourably with pristine graphene in aqueous solution molecular dynamics simulations; conversely, a model of the hydrophilic 010 surface of cellulose exhibited progressive rearrangement to present a more hydrophobic face with the graphene, with weakened hydrogen bonds between cellulose chains and partial permeation of water. Here, we extend this work by simulating the interaction in aqueous solution of the amphiphilic 110 surface of a cellulose Iβ microfibril model, comprising 36 chains of 40 glucosyl residues, with an infinite sheet of pristine graphene. This face of the microfibril is of intermediate hydrophilicity and progressively associates with graphene over replicate simulations. As cellulose chains adhere to the graphene surface, forming interactions via its CH and OH groups, we observe a degree of local and global untwisting of the microfibril. Complementary rippling of the graphene surface is also observed, as it adapts to interaction with the microfibril. This adsorption process is accompanied by increased exclusion of water between cellulose and graphene although some water localises between chains at the immediate interface. The predicted propensity of a cellulose microfibril to adsorb spontaneously on the graphene surface, with mutual structural accommodation, highlights the amphiphilic nature of cellulose and the types of interactions that can be harnessed to design new graphene-carbohydrate biopolymer materials.

Molecular Dynamics Simulations on the Elastic Properties of Polypropylene Bionanocomposite Reinforced with Cellulose Nanofibrils

Cellulose-reinforced polypropylene bionanocomposites can show improved elastic properties over their pure polypropylene counterparts. We have used equilibrium and non-equilibrium molecular dynamics (MD) simulations to study the elastic properties of polypropylene bionanocomposite systems composed of cellulose nanofibrils (CNF), polypropylene (PP) matrix, and maleic anhydride (MAH) coupling agent. The components of the bionanocomposite were parametrized for compatibility with the AMBER14SB force fields. The elastic properties of pure PP systems converge for the chains with at least 20 monomers. The ratio of cellulose in CNF-PP bionanocomposites strongly affects their elastic properties. The elastic modulus of CNF-PP bionanocomposites shows small improvement when the adhesion between hydrophobic and hydrophilic components is facilitated by a MAH coupling agent. The results demonstrate how fully-atomistic MD simulations can be systematically used to evaluate the elastic properties of CNF-PP bionanocomposites and to make predictions that are in agreement with experiments.

Building an extensible cell wall

This article recounts, from my perspective of four decades in this field, evolving paradigms of primary cell wall structure and the mechanism of surface enlargement of growing cell walls. Updates of the structures, physical interactions, and roles of cellulose, xyloglucan, and pectins are presented. This leads to an example of how a conceptual depiction of wall structure can be translated into an explicit quantitative model based on molecular dynamics methods. Comparison of the model's mechanical behavior with experimental results provides insights into the molecular basis of complex mechanical behaviors of primary cell wall and uncovers the dominant role of cellulose-cellulose interactions in forming a strong yet extensible network.

Cellulose Nanocrystals: Tensile Strength and Failure Mechanisms Revealed Using Reactive Molecular Dynamics

Cellulose nanocrystals (CNCs) offer excellent mechanical properties. However, measuring the strength by performing reliable experiments at the nanoscale is challenging. In this paper, we model Iβ crystalline cellulose using reactive molecular dynamics simulations. Taking the fibril twist into account, structural changes and hydrogen-bonding characteristics of CNCs during the tensile test are inspected and the failure mechanism of CNCs is analyzed down to the scale of individual bonds. The C4-O4 glycosidic bond is found to be responsible for the failure of CNCs. Finally, the effect of strain rate on ultimate properties is analyzed and a nonlinear model is used to predict the ultimate strength of 9.2 GPa and ultimate strain of 8.5% at a 1 s^-1 strain rate. This study sheds light on the applications of cellulose in nanocomposites and further modeling of cellulose nanofibres.

Dietary Fiber in Plant Cell Walls—The Healthy Carbohydrates

Dietary fiber (DF) is one of the major classes of nutrients for humans. It is widely distributed in the edible parts of natural plants, with the cell wall being the main DF-containing structure. The DF content varies significantly in different plant species and organs, and the processing procedure can have a dramatic effect on the DF composition of plant-based foods. Given the considerable nutritional value of DF, a deeper understanding of DF in food plants, including its composition and biosynthesis, is fundamental to the establishment of a daily intake reference of DF and is also critical to molecular breeding programs for modifying DF content. In the past decades, plant cell wall biology has seen dramatic progress, and such knowledge is of great potential to be translated into DF-related food science research and may provide future research directions for improving the health benefits of food crops. In this review, to spark interdisciplinary discussions between food science researchers and plant cell wall biologists, we focus on a specific category of DF—cell wall carbohydrates. We first summarize the content and composition of carbohydrate DF in various plant-based foods, and then discuss the structure and biosynthesis mechanism of each carbohydrate DF category, in particular the respective biosynthetic enzymes. Health impacts of DF are highlighted, and finally, future directions of DF research are also briefly outlined.

Carbohydrate Force Fields: The Role of Small Partial Atomic Charges in Preventing Conformational Collapse

Although the quality of current additive all-atom force fields for carbohydrates has been demonstrated in many applications, occasional significant differences reported for the hydrodynamic behavior of specific polysaccharides modeled with different force fields is a cause for concern. In particular, irreversible conformational collapse has been noted for some polysaccharide simulations with the GLYCAM06j force field. Here, we investigate the cause of this phenomenon through comparative simulations of a range of saccharides with both the GLYCAM06j and the CHARMM36 carbohydrate force fields. We find that conformational collapse in GLYCAM06j occurs for saccharide chains containing the deoxy sugar α-l-rhamnose after relatively long simulation intervals. Further, we explore the mechanism of conformational collapse and show that this phenomenon arises because of the anomalous low energy in GLYCAM06j (as compared to quantum mechanical calculations) of a specific orientation of α-l-Rha to α-l-Rha glycosidic linkages, which are subsequently sustained by intramolecular interactions in the saccharide chain. We identify the lack of partial charges on aliphatic hydrogens in GLYCAM as the source of this anomaly, demonstrating that addition of small partial atomic charges on the aliphatic protons in rhamnose removes the conformational collapse phenomenon. This work reveals the large cumulative impact that small partial charges may have on the dynamic behavior of polysaccharides and indicates that future reparameterization of the GLYCAM06j force field should investigate the addition of partial charges on all aliphatic hydrogens.

Combining Computational Chemistry and Crystallography for a Better Understanding of the Structure of Cellulose

The approaches in this article seek to enhance understanding of cellulose at the molecular level, independent of the source and the particular crystalline form of cellulose. Four main areas of structure research are reviewed. Initially, the molecular shape is inferred from the crystal structures of many small molecules that have β-(1→4) linkages. Then, conformational analyses with potential energy calculations of cellobiose are covered, followed by the use of Atoms-In-Molecules theory to learn about interactions in experimental and theoretical structures. The last section covers models of cellulose nanoparticles. Controversies addressed include the stability of twofold screw-axis conformations, the influence of different computational methods, the predictability of crystalline conformations by studies of isolated molecules, and the twisting of model cellulose crystals.

Cellulose-hemicellulose interactions - A nanoscale view

In this work, we study interactions of five different hemicellulose models, i.e. Galactoglucomannan, O-Acetyl-Galactoglucomannan, Fuco-Galacto-Xyloglucan, 4-O-Methylglucuronoxylan, and 4-O-Methylglucuronoarabinoxylan, and their respective binding strength to cellulose nanocrystals by molecular dynamics simulations. Glucuronoarabinoxylan showed the highest free energy of binding, whereas Xyloglucan had the lowest interaction energies amongst the five models. We further performed simulated shear tests and concluded that failure mostly happens at the inter-molecular interaction level within the hemicellulose fraction, rather than at the interface with cellulose. The presence of water molecules seems to have a weakening effect on the interactions of hemicellulose and cellulose, taking up the available hydroxyl groups on the surface of the cellulose for hydrogen bonding. We believe that these studies can shed light on better understanding of plant cell walls, as well as providing evidence on variability of the structures of different plant sources for extractions, purification, and operation of biorefineries.

Strengthening and Toughening Hierarchical Nanocellulose via Humidity-Mediated Interface

Undoubtedly humidity is a non-negligible and sensitive problem for cellulose, which is usually regarded as one disadvantage to cellulose-based materials because of the uncontrolled deformation and mechanical decline. But the lack of an in-depth understanding of the interfacial behavior of nanocellulose in particular makes it challenging to maintain anticipated performance for cellulose-based materials under varied relative humidity (RH). Starting from multiscale mechanics, we herein carry out first-principles calculations and large-scale molecular dynamics simulations to demonstrate the humidity-mediated interface in hierarchical cellulose nanocrystals (CNCs) and associated deformation modes. More intriguingly, the simulations and subsequent experiments reveal that water molecules (moisture) as the interfacial media can strengthen and toughen nanocellulose simultaneously within a suitable range of RH. From the perspective of interfacial design in materials, the anomalous mechanical behavior of nanocellulose with humidity-mediated interfaces indicates that flexible hydrogen bonds (HBs) play a pivotal role in the interfacial sliding. The difference between CNC-CNC HBs and CNC-water-CNC HBs triggers the humidity-mediated interfacial slipping in nanocellulose, resulting in the arising of a pronounced strain hardening stage and the suppression of strain localization during uniaxial tension. This inelastic deformation of nanocellulose with humidity-mediated interfaces is similar to the Velcro-like behavior of a wet wood cell wall. Our investigations give evidence that the humidity-mediated interface can promote the mechanical enhancement of nanocellulose, which would provide a promising strategy for the bottom-up design of cellulose-based materials with tailored mechanical properties.

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.

Benchmark assessment of molecular geometries and energies from small molecule force fields

Background: Force fields are used in a wide variety of contexts for classical molecular simulation, including studies on protein-ligand binding, membrane permeation, and thermophysical property prediction. The quality of these studies relies on the quality of the force fields used to represent the systems. Methods: Focusing on small molecules of fewer than 50 heavy atoms, our aim in this work is to compare nine force fields: GAFF, GAFF2, MMFF94, MMFF94S, OPLS3e, SMIRNOFF99Frosst, and the Open Force Field Parsley, versions 1.0, 1.1, and 1.2. On a dataset comprising 22,675 molecular structures of 3,271 molecules, we analyzed force field-optimized geometries and conformer energies compared to reference quantum mechanical (QM) data. Results: We show that while OPLS3e performs best, the latest Open Force Field Parsley release is approaching a comparable level of accuracy in reproducing QM geometries and energetics for this set of molecules. Meanwhile, the performance of established force fields such as MMFF94S and GAFF2 is generally somewhat worse. We also find that the series of recent Open Force Field versions provide significant increases in accuracy. Conclusions: This study provides an extensive test of the performance of different molecular mechanics force fields on a diverse molecule set, and highlights two (OPLS3e and OpenFF 1.2) that perform better than the others tested on the present comparison. Our molecule set and results are available for other researchers to use in testing.

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.

Insight into Cellulose Dissolution with the Tetrabutylphosphonium Chloride-Water Mixture using Molecular Dynamics Simulations

All-atom molecular dynamics simulations are utilized to determine the properties and mechanisms of cellulose dissolution using the ionic liquid tetrabutylphosphonium chloride (TBPCl)–water mixture, from 63.1 to 100 mol % water. The hydrogen bonding between small and large cellulose bundles with 18 and 88 strands, respectively, is compared for all concentrations. The Cl, TBP, and water enable cellulose dissolution by working together to form a cooperative mechanism capable of separating the cellulose strands from the bundle. The chloride anions initiate the cellulose breakup, and water assists in delaying the cellulose strand reformation; the TBP cation then more permanently separates the cellulose strands from the bundle. The chloride anion provides a net negative pairwise energy, offsetting the net positive pairwise energy of the peeling cellulose strand. The TBP–peeling cellulose strand has a uniquely favorable and potentially net negative pairwise energy contribution in the TBPCl–water solution, which may partially explain why it is capable of dissolving cellulose at moderate temperatures and high water concentrations. The cellulose dissolution declines rapidly with increasing water concentration as hydrogen bond lifetimes of the chloride–cellulose hydroxyl hydrogens fall below the cellulose’s largest intra-strand hydrogen bonding lifetime.

Wood–Moisture Relationships Studied with Molecular Simulations: Methodological Guidelines

This paper aims at providing a methodological framework for investigating wood polymers using atomistic modeling, namely, molecular dynamics (MD) and grand canonical Monte Carlo (GCMC) simulations. Atomistic simulations are used to mimic water adsorption and desorption in amorphous polymers, make observations on swelling, mechanical softening, and on hysteresis. This hygromechanical behavior, as observed in particular from the breaking and reforming of hydrogen bonds, is related to the behavior of more complex polymeric composites. Wood is a hierarchical material, where the origin of wood-moisture relationships lies at the nanoporous material scale. As water molecules are adsorbed into the hydrophilic matrix in the cell walls, the induced fluid–solid interaction forces result in swelling of these cell walls. The interaction of the composite polymeric material, that is the layer S2 of the wood cell wall, with water is known to rearrange its internal material structure, which makes it moisture sensitive, influencing its physical properties. In-depth studies of the coupled effects of water sorption on hygric and mechanical properties of different polymeric components can be performed with atomistic modeling. The paper covers the main components of knowledge and good practice for such simulations.

Molecular insight into the wetting behavior and amphiphilic character of cellulose nanocrystals

The study of nanocellulose is a field of growing interest due to its many applications and its use in the development of biocompatible and eco-friendly materials. In spite of the vast number of studies in the field, many questions about the role of the molecular structure in the properties of cellulose are still subject of debate. One of these fundamental questions is the possible amphiphilic nature of cellulose and the relative role of hydrogen bonding and hydrophobic effect on the interactions of cellulose. In this work we present an extensive molecular dynamics simulation study of this question by analyzing the wetting of cellulose with water and organic solvent, its interaction with hydrophilic and hydrophobic ions and its interaction with a protein (human epidermal growth factor, hEGF). We consider two characteristic cellulose crystal planes of Iβ cellulose with very different roughness, different hydrogen bonding capability and different exposure of cellulose hydrophobic groups (the (010) plane which has exposed -OH groups and the (100) plane with buried -OH groups). Our results show that both surfaces are simultaneously hydrophilic and lipophilic, with both surfaces having very similar contact angles. In spite of the global similarity of wetting of both surfaces, the molecular details of wetting are very different and substantial local wetting heterogeneities (which strongly depend on the surface) appear for both solvents. We also observe a weak interaction of both surfaces with hydrophobic and hydrophilic solutes. These weak interactions are attributed to the simultaneous lipophilic and hydrophilic character of both (100) and (010) cellulose surfaces. Interestingly, we found a substantial interaction of both cellulose planes with polar and apolar residues of the hEGF protein.

Predicting the Structures of Glycans, Glycoproteins, and Their Complexes

Complex carbohydrates are ubiquitous in nature, and together with proteins and nucleic acids they comprise the building blocks of life. But unlike proteins and nucleic acids, carbohydrates form nonlinear polymers, and they are not characterized by robust secondary or tertiary structures but rather by distributions of well-defined conformational states. Their molecular flexibility means that oligosaccharides are often refractory to crystallization, and nuclear magnetic resonance (NMR) spectroscopy augmented by molecular dynamics (MD) simulation is the leading method for their characterization in solution. The biological importance of carbohydrate-protein interactions, in organismal development as well as in disease, places urgency on the creation of innovative experimental and theoretical methods that can predict the specificity of such interactions and quantify their strengths. Additionally, the emerging realization that protein glycosylation impacts protein function and immunogenicity places the ability to define the mechanisms by which glycosylation impacts these features at the forefront of carbohydrate modeling. This review will discuss the relevant theoretical approaches to studying the three-dimensional structures of this fascinating class of molecules and interactions, with reference to the relevant experimental data and techniques that are key for validation of the theoretical predictions.

Relative Contributions of Solubility and Mobility to the Stability of Amorphous Solid Dispersions of Poorly Soluble Drugs: A Molecular Dynamics Simulation Study

Amorphous solid dispersions are considered a promising formulation strategy for the oral delivery of poorly soluble drugs. The limiting factor for the applicability of this approach is the physical (in)stability of the amorphous phase in solid samples. Minimizing the risk of reduced shelf life for a new drug by establishing a suitable excipient/polymer-type from first principles would be desirable to accelerate formulation development. Here, we perform Molecular Dynamics simulations to determine properties of blends of eight different polymer–small molecule drug combinations for which stability data are available from a consistent set of literature data. We calculate thermodynamic factors (mixing energies) as well as mobilities (diffusion rates and roto-vibrational fluctuations). We find that either of the two factors, mobility and energetics, can determine the relative stability of the amorphous form for a given drug. Which factor is rate limiting depends on physico-chemical properties of the drug and the excipients/polymers. The methods outlined here can be readily employed for an in silico pre-screening of different excipients for a given drug to establish a qualitative ranking of the expected relative stabilities, thereby accelerating and streamlining formulation development.

Validation and Comparison of Force Fields for Native Cyclodextrins in Aqueous Solution

Molecular dynamics simulations of native α-, β-, and γ-cyclodextrin in aqueous solution have been conducted with the goal to investigate the performance of the CHARMM36 force field, the AMBER-compatible q4md-CD force field, and five variants of the GROMOS force field. The properties analyzed are structural parameters derived from X-ray diffraction and NMR experiments as well as hydrogen bonds and hydration patterns, including hydration free enthalpies. Recent revisions of the torsional-angle parameters for carbohydrate systems within the GROMOS family of force fields lead to a significant improvement of the agreement between simulated and experimental NMR data. Therefore, we recommend using the variant 53A6_GLYC instead of 53A6 and 56A6_{CARBO_R} or 2016H66 instead of 56A6_CARBO to simulate cyclodextrins in solution. The CHARMM36 and q4md-CD force fields show a similar performance as the three recommended GROMOS parameter sets. A significant difference is the more flexible nature of the cyclodextrins modeled with the CHARMM36 and q4md-CD force fields compared to the three recommended GROMOS parameter sets.

Hydrogen bonds and twist in cellulose microfibrils

There is increasing experimental and computational evidence that cellulose microfibrils can exist in a stable twisted form. In this study, atomistic molecular dynamics (MD) simulations are performed to investigate the importance of intrachain hydrogen bonds on the twist in cellulose microfibrils. We systematically enforce or block the formation of these intrachain hydrogen bonds by either constraining dihedral angles or manipulating charges. For the majority of simulations a consistent right handed twist is observed. The exceptions are two sets of simulations that block the O2-O6' intrachain hydrogen bond, where no consistent twist is observed in multiple independent simulations suggesting that the O2-O6' hydrogen bond can drive twist. However, in a further simulation where exocyclic group rotation is also blocked, right-handed twist still develops suggesting that intrachain hydrogen bonds are not necessary to drive twist in cellulose microfibrils.

Computing converged free energy differences between levels of theory via nonequilibrium work methods: Challenges and opportunities

We demonstrate that Jarzynski's equation can be used to reliably compute free energy differences between low and high level representations of systems. The need for such a calculation arises when employing the so-called "indirect" approach to free energy simulations with mixed quantum mechanical/molecular mechanical (QM/MM) Hamiltonians; a popular technique for circumventing extensive simulations involving quantum chemical computations. We have applied this methodology to several small and medium sized organic molecules, both in the gas phase and explicit solvent. Test cases include several systems for which the standard approach; that is, free energy perturbation between low and high level description, fails to converge. Finally, we identify three major areas in which the difference between low and high level representations make the calculation of ΔAlow→high difficult: bond stretching and angle bending, different preferred conformations, and the response of the MM region to the charge distribution of the QM region. © 2016 Wiley Periodicals, Inc.

Calculations of the absolute binding free energies for Ralstonia solanacearum lectins bound with methyl-α-l-fucoside at molecular mechanical and quantum mechanical/molecular mechanical levels

In this work, both a molecular mechanical (MM) method and a hybrid quantum mechanical/molecular mechanical (QM/MM) method have been applied in the study of the binding affinities of methyl-α-l-fucoside to Ralstonia solanacearum lectins.

Diffraction pattern simulation of cellulose fibrils using distributed and quantized pair distances

Intensity simulation of X-ray scattering from large twisted cellulose molecular fibrils is important in understanding the impact of chemical or physical treatments on structural properties such as twisting or coiling. This paper describes a highly efficient method for the simulation of X-ray diffraction patterns from complex fibrils using atom-type-specific pair-distance quantization. Pair distances are sorted into arrays which are labelled by atom type. Histograms of pair distances in each array are computed and binned and the resulting population distributions are used to represent the whole pair-distance data set. These quantized pair-distance arrays are used with a modified and vectorized Debye formula to simulate diffraction patterns. This approach utilizes fewer pair distances in each iteration, and atomic scattering factors are moved outside the iteration since the arrays are labelled by atom type. This algorithm significantly reduces the computation time while maintaining the accuracy of diffraction pattern simulation, making possible the simulation of diffraction patterns from large twisted fibrils in a relatively short period of time, as is required for model testing and refinement.

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.

Unraveling cellulose microfibrils: a twisted tale

Molecular dynamics (MD) simulations of cellulose microfibrils are pertinent to the paper, textile, and biofuels industries for their unique capacity to characterize dynamic behavior and atomic-level interactions with solvent molecules and cellulase enzymes. While high-resolution crystallographic data have established a solid basis for computational analysis of cellulose, previous work has demonstrated a tendency for modeled microfibrils to diverge from the linear experimental structure and adopt a twisted conformation. Here, we investigate the dependence of this twisting behavior on computational approximations and establish the theoretical basis for its occurrence. We examine the role of solvent, the effect of nonbonded force field parameters [partial charges and van der Waals (vdW) contributions], and the use of explicitly modeled oxygen lone pairs in both the solute and solvent. Findings suggest that microfibril twisting is favored by vdW interactions, and counteracted by both intrachain hydrogen bonds and solvent effects at the microfibril surface.

Optimizing Solute-Solute Interactions in the GLYCAM06 and CHARMM36 Carbohydrate Force Fields Using Osmotic Pressure Measurements

GLYCAM06 and CHARMM36 are successful force fields for modeling carbohydrates. To correct recently identified deficiencies with both force fields, we adjusted intersolute nonbonded parameters to reproduce the experimental osmotic coefficient of glucose at 1 M. The modified parameters improve behavior of glucose and sucrose up to 4 M and improve modeling of a dextran 55-mer. While the modified parameters may not be applicable to all carbohydrates, they highlight the use of osmotic simulations to optimize force fields.

Cellulose Nanofibrils and Mechanism of their Mineralization in Biomimetic Synthesis of Hydroxyapatite/Native Bacterial Cellulose Nanocomposites: Molecular Dynamics Simulations

Molecular dynamics (MD) simulation of a nanofibril of native bacterial cellulose (BC) in solutions of mineral ions is presented. The supersaturated calcium-phosphate (CP) solution with the ionic composition of hydroxyapatite and CaCl2 solutions with the concentrations below, equal to, and above the solubility limits are simulated. The influence of solvation models (TIP3P and TIP4P-ew water models) on structural characteristics of the simulated nanofibril and on the crystal nucleation process is assessed. The structural characteristics of cellulose nanofibrils (in particular, of the surface layer) are found to be nearly independent of the solvation models used in the simulation and on the presence of ions in the solutions. It is shown that ionic clusters are formed in the solution rather than on the fibril surface. The cluster sizes are slightly different for the two water models. The effect of the ion-ion interaction parameters on the results is discussed. The main conclusion is that the activity of hydroxyl groups on the BC fibril surface is not high enough to cause adsorption of Ca(2+) ions from the solution. Therefore, the nucleation of CP crystals takes place initially in solution, and then the crystallites formed can be adsorbed on BC nanofibril surfaces.

The structure and dynamics of chitin nanofibrils in an aqueous environment revealed by molecular dynamics simulations

The structure–property relations reveal the typical size of chitin nanofibrils observed in natural systems.

Electric Interfacial Layer of Modified Cellulose Nanocrystals in Aqueous Electrolyte Solution: Predictions by the Molecular Theory of Solvation

The X-ray crystal structure-based models of Iα cellulose nanocrystals (CNC), both pristine and containing surface sulfate groups with negative charge 0-0.34 e/nm(2) produced by sulfuric acid hydrolysis of softwood pulp, feature a highly polarized "crystal-like" charge distribution. We perform sampling using molecular dynamics (MD) of the structural relaxation of neutral pristine and negatively charged sulfated CNC of various lengths in explicit water solvent and then employ the statistical mechanical 3D-RISM-KH molecular theory of solvation to evaluate the solvation structure and thermodynamics of the relaxed CNC in ambient aqueous NaCl solution at a concentration of 0.0-0.25 mol/kg. The MD sampling induces a right-hand twist in CNC and rearranges its initially ordered structure with a macrodipole of high-density charges at the opposite faces into small local spots of alternating charge at each face. This surface charge rearrangement observed for both neutral and charged CNC significantly affects the distribution of ions around CNC in aqueous electrolyte solution. The solvation free energy (SFE) of charged sulfated CNC has a minimum at a particular electrolyte concentration depending on the surface charge density, whereas the SFE of neutral CNC increases linearly with NaCl concentration. The SFE contribution from Na(+) counterions exhibits behavior similar to the NaCl concentration dependence of the whole SFE. An analysis of the 3D maps of Na(+) density distributions shows that these model CNC particles exhibit the behavior of charged nanocolloids in aqueous electrolyte solution: an increase in electrolyte concentration shrinks the electric interfacial layer and weakens the effective repulsion between charged CNC particles. The 3D-RISM-KH method readily treats solvent and electrolyte of a given nature and concentration to predict effective interactions between CNC particles in electrolyte solution. We provide CNC structural models and a modeling procedure for studies of effective interactions and the formation of ordered phases of CNC suspensions in electrolyte solution.