REACH Coarse-Grained Simulation of a Cellulose Fiber

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.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

A Review of Cellulose Coarse-Grained Models and Their Applications

Cellulose is the most common biopolymer and widely used in our daily life. Due to its unique properties and biodegradability, it has been attracting increased attention in the recent years and various new applications of cellulose and its derivatives are constantly being found. The development of new materials with improved properties, however, is not always an easy task, and theoretical models and computer simulations can often help in this process. In this review, we give an overview of different coarse-grained models of cellulose and their applications to various systems. Various coarse-grained models with different mapping schemes are presented, which can efficiently simulate systems from the single cellulose fibril/crystal to the assembly of many fibrils/crystals. We also discuss relevant applications of these models with a focus on the mechanical properties, self-assembly, chiral nematic phases, conversion between cellulose allomorphs, composite materials and interactions with other molecules.

Extending the Martini Coarse-Grained Force Field to N-Glycans

Glycans play a vital role in a large number of cellular processes. Their complex and flexible nature hampers structure-function studies using experimental techniques. Molecular dynamics (MD) simulations can help in understanding dynamic aspects of glycans if the force field parameters used can reproduce key experimentally observed properties. Here, we present optimized coarse-grained (CG) Martini force field parameters for N-glycans, calibrated against experimentally derived binding affinities for lectins. The CG bonded parameters were obtained from atomistic (ATM) simulations for different glycan topologies including high mannose and complex glycans with various branching patterns. In the CG model, additional elastic networks are shown to improve maintenance of the overall conformational distribution. Solvation free energies and octanol-water partition coefficients were also calculated for various N-glycan disaccharide combinations. When using standard Martini nonbonded parameters, we observed that glycans spontaneously aggregated in the solution and required down-scaling of their interactions for reproduction of ATM model radial distribution functions. We also optimized the nonbonded interactions for glycans interacting with seven lectin candidates and show that a relatively modest scaling down of the glycan-protein interactions can reproduce free energies obtained from experimental studies. These parameters should be of use in studying the role of glycans in various glycoproteins and carbohydrate binding proteins as well as their complexes, while benefiting from the efficiency of CG sampling.

Recent Advances in Coarse-Grained Models for Biomolecules and Their Applications

Molecular dynamics simulations have emerged as a powerful tool to study biological systems at varied length and timescales. The conventional all-atom molecular dynamics simulations are being used by the wider scientific community in routine to capture the conformational dynamics and local motions. In addition, recent developments in coarse-grained models have opened the way to study the macromolecular complexes for time scales up to milliseconds. In this review, we have discussed the principle, applicability and recent development in coarse-grained models for biological systems. The potential of coarse-grained simulation has been reviewed through state-of-the-art examples of protein folding and structure prediction, self-assembly of complexes, membrane systems and carbohydrates fiber models. The multiscale simulation approaches have also been discussed in the context of their emerging role in unravelling hierarchical level information of biosystems. We conclude this review with the future scope of coarse-grained simulations as a constantly evolving tool to capture the dynamics of biosystems.

Coarse-grained molecular dynamics simulations of α-1,3-glucan

In this paper we present a computational study of aggregation in aqueous solutions of α-1,3-glucan captured using a coarse-grained (CG) model that can be extended to other polysaccharides. This CG model captures atomistic geometry (i.e., relative placement of the hydrogen bonding donors and acceptors within the monomer) of the α-1,3-glucan monomer, the directional interactions due to the donor-acceptor hydrogen bonds, and their effect on aggregation of multiple α-1,3-glucan chains without the extensive computational resources needed for simulations with atomistic models. Using this CG model, we conduct molecular dynamics simulations to assess the effect of varying α-1,3-glucan chain length and hydrogen bond interaction strengths on the aggregation of multiple chains at finite concentrations in implicit solvent. We quantify the hydrogen bonding strength needed for multiple chains to aggregate, the distribution of inter- and intra-chain hydrogen bonds within the aggregate and in some cases, the shapes of the aggregate. We also explore the effect of substitution/silencing of some randomly selected or specific hydrogen bonding sites in the chain on the aggregation and aggregate structure. In the unmodified α-1,3-glucan solution, the inter-chain hydrogen bonds cause the chains to aggregate into sheets. Random silencing of hydrogen bonding donor sites only increases the hydrogen bond strength needed for aggregation but retains the same aggregate structure as the unmodified chains. Specific silencing of the hydrogen-bonding site on the C6 carbon leads to the chains aggregating into planar sheets that then fold over to form hollow cylinders at intermediate hydrogen bond strength - 4.7 to 5.3 kcal mol-1. These cylindrical aggregates assemble end-to-end to form larger aggregates at higher hydrogen bond strengths.

Mechanical Properties of Cellulose Nanocrystal (CNC) Bundles: Coarse-Grained Molecular Dynamic Simulation

Cellulose nanocrystals (CNCs) is a promising biodegradable nanomaterial with outstanding physical, chemical, and mechanical properties for many applications. Although aligned CNCs can self-assemble into bundles, their mechanical performance is reduced by interfacial strength between CNCs and a twisted structure. In this paper, we employ developed coarse-grained (CG) molecular dynamic (MD) simulations to investigate the influence of twist and interface energy on the tensile performance of CNC bundles. CNC bundles of different sizes (number of particles) are tested to also include the effect of size on mechanical performance. The effect of interfacial energy and twist on the mechanical performance shows that elastic modulus, strength, and toughness are more sensitive to twisted angle than interfacial energy. In addition, the effect of size on the bundle and twist on their mechanical performance revealed that both size and twist have a significant effect on the results and can reduce the strength and elastic modulus by 75% as a results of covalent bond dissociation. In addition, a comparison of the broken regions for different values of twist shows that by increasing the twist angle the crack propagates in multiple locations with a twisted shape.

Generalization of the elastic network model for the study of large conformational changes in biomolecules

The elastic network (EN) is a prime model that describes the long-time dynamics of biomolecules. However, the use of harmonic potentials renders this model insufficient for studying large conformational changes of proteins (e.g. stretching of proteins, folding and thermal unfolding). Here, we extend the capabilities of the EN model by using a harmonic approximation described by Lennard-Jones (LJ) interactions for far contacts and native contacts obtained from the standard overlap criterion as in the case of Gō-like models. While our model is validated against the EN model by reproducing the equilibrium properties for a number of proteins, we also show that the model is suitable for the study of large conformation changes by providing various examples. In particular, this is illustrated on the basis of pulling simulations that predict with high accuracy the experimental data on the rupture force of the studied proteins. Furthermore, in the case of DDFLN4 protein, our pulling simulations highlight the advantages of our model with respect to Gō-like approaches, where the latter fail to reproduce previous results obtained by all-atom simulations that predict an additional characteristic peak for this protein. In addition, folding simulations of small peptides yield different folding times for α-helix and β-hairpin, in agreement with experiment, in this way providing further opportunities for the application of our model in studying large conformational changes of proteins. In contrast to the EN model, our model is suitable for both normal mode analysis and molecular dynamics simulation. We anticipate that the proposed model will find applications in a broad range of problems in biology, including, among others, protein folding and thermal unfolding.

Elastic moduli of biological fibers in a coarse-grained model: crystalline cellulose and β-amyloids

We study the mechanical response of cellulose and β-amyloid microfibrils to three types of deformation: tensile, indentational, and shear. The cellulose microfibrils correspond to the allomorphs Iα or Iβ whereas the β-amyloid microfibrils correspond to the polymorphs of either two- or three-fold symmetry. This response can be characterized by three elastic moduli, namely, Y_L, Y_T, and S. We use a structure-based coarse-grained model to analyze the deformations in a unified manner. We find that each of the moduli is almost the same for the two allomorphs of cellulose but Y_L is about 20 times larger than Y_T (140 GPa vs. 7 GPa), indicating the existence of significant anisotropy. For cellulose we note that the anisotropy results from the involvement of covalent bonds in stretching. For β-amyloid, the sense of anisotropy is opposite to that of cellulose. In the three-fold symmetry case, Y_L is about half of Y_T (3 vs. 7) whereas for two-fold symmetry the anisotropy is much larger (1.6 vs. 21 GPa). The S modulus is derived to be 1.2 GPa for three-fold symmetry and one half of it for the other symmetry and 3.0 GPa for cellulose. The values of the moduli reflect deformations in the hydrogen-bond network. Unlike in our theoretical approach, no experiment can measure all three elastic moduli with the same apparatus. However, our theoretical results are consistent with various measured values: typical Y_L for cellulose Iβ ranges from 133 to 155 GPa, Y_T from 2 to 25 GPa, and S from 1.8 to 3.8 GPa. For β-amyloid, the experimental values of S and Y_T are about 0.3 GPa and 3.3 GPa respectively, while the value of Y_L has not been reported.

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.

A coarse-grain force-field for xylan and its interaction with cellulose

We have built a coarse-grain (CG) model describing xylan and its interaction with crystalline cellulose surfaces. Each xylosyl or glucosyl unit was represented by a single grain. Our calculations rely on force-field parameters adapted from the atomistic description of short xylan fragments and their adsorption on cellulose. This CG model was first validated for xylan chains both isolated and in the bulk where a good match was found with its atomistic counterpart as well as with experimental measurements. A similar agreement was also found when short xylan fragments were adsorbed on the (110) surface of crystalline cellulose. The CG model, which was extended to the (100) and (1-10) surfaces, revealed that the adsorbed xylan, which was essentially extended in the atomistic situation, could also adopt coiled structures, especially when laying on the hydrophobic cellulose surfaces.

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β.

Biological materials by design

In this topical review we discuss recent advances in the use of physical insight into the way biological materials function, to design novel engineered materials 'from scratch', or from the level of fundamental building blocks upwards and by using computational multiscale methods that link chemistry to material function. We present studies that connect advances in multiscale hierarchical material structuring with material synthesis and testing, review case studies of wood and other biological materials, and illustrate how engineered fiber composites and bulk materials are designed, modeled, and then synthesized and tested experimentally. The integration of experiment and simulation in multiscale design opens new avenues to explore the physics of materials from a fundamental perspective, and using complementary strengths from models and empirical techniques. Recent developments in this field illustrate a new paradigm by which complex material functionality is achieved through hierarchical structuring in spite of simple material constituents.

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.