Jakob Wohlert

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.

Free energy of a trans-membrane pore calculated from atomistic molecular dynamics simulations

Atomistic molecular dynamics simulations of a lipid bilayer were performed to calculate the free energy of a trans-membrane pore as a function of its radius. The free energy was calculated as a function of a reaction coordinate using a potential of mean constraint force. The pore radius was then calculated from the reaction coordinate using Monte Carlo particle insertions. The main characteristics of the free energy that comes out of the simulations are a quadratic shape for a radius less than about 0.3 nm, a linear shape for larger radii than this, and a rather abrupt change without local minima or maxima between the two regions. In the outer region, a line tension can be calculated, which is consistent with the experimentally measured values. Further, this line tension can be rationalized and understood in terms of the energetic cost for deforming a part of the lipid bilayer into a hydrophilic pore. The region with small radii can be described and understood in terms of statistical mechanics of density fluctuations. In the region of crossover between a quadratic and linear free energy there was some hysteresis associated with filling and evacuation of the pore with water. The metastable prepore state hypothesized to interpret the experiments was not observed in this region.

Simulation studies of the insolubility of cellulose

Molecular dynamics simulations have been used to calculate the potentials of mean force for separating short cellooligomers in aqueous solution as a means of estimating the contributions of hydrophobic stacking and hydrogen bonding to the insolubility of crystalline cellulose. A series of four potential of mean force (pmf) calculations for glucose, cellobiose, cellotriose, and cellotetraose in aqueous solution were performed for situations in which the molecules were initially placed with their hydrophobic faces stacked against one another, and another for the cases where the molecules were initially placed adjacent to one another in a co-planar, hydrogen-bonded arrangement, as they would be in cellulose Ibeta. From these calculations, it was found that hydrophobic association does indeed favor a crystal-like structure over solution, as might be expected. Somewhat more surprisingly, hydrogen bonding also favored the crystal packing, possibly in part because of the high entropic cost for hydrating glucose hydroxyl groups, which significantly restricts the configurational freedom of the hydrogen-bonded waters. The crystal was also favored by the observation that there was no increase in chain configurational entropy upon dissolution, because the free chain adopts only one conformation, as previously observed, but against intuitive expectations, apparently due to the persistence of the intramolecular O3-O5 hydrogen bond.

Dynamics in atomistic simulations of phospholipid membranes: Nuclear magnetic resonance relaxation rates and lateral diffusion

It is shown that a long, near microsecond, atomistic simulation can shed some light upon the dynamical processes occurring in a lipid bilayer. The analysis focuses on reorientational dynamics of the chains and lateral diffusion of lipids. It is shown that the reorientational correlation functions exhibits an algebraic decay (rather than exponential) for several orders of magnitude in time. The calculated nuclear magnetic resonance relaxation rates agree with experiments for carbons at the C7 position while there are some differences for C3. Lateral diffusion can be divided into two stages. In a first stage occurring at short times, t<5 ns, the center of mass of the lipid moves due to conformational changes of the chains while the headgroup position remains relatively fixed. In this stage, the center of mass can move up to approximately 0.8 nm. The fitted short-time diffusion coefficient is D(1)=13 x 10(-7) cm(2) s(-1) On a longer time scale, the diffusion coefficient becomes D(2)=0.79 x 10(-7) cm(2) s(-1).

Dynamics of cellulose-water interfaces: NMR spin-lattice relaxation times calculated from atomistic computer simulations.

Solid-state nuclear magnetic resonance (CP/MAS 13C NMR) spectroscopy has often been used to study cellulose structure, but some features of the cellulose NMR spectrum are not yet fully understood. One such feature is a doublet around 84 ppm, a signal that has been proposed to originate from C4 atoms at cellulose fibril surfaces. The two peaks yield different T1, differing by approximately a factor of 2 at 75 MHz. In this study, we calculate T1 from C4-H4 vector dynamics obtained from molecular dynamics computer simulations of cellulose I beta-water interfacial systems. Calculated and experimentally obtained T1 values for C4 atoms in surface chains fell within the same order of magnitude, 3-20 s. This means that the applied force field reproduces relevant surface dynamics for the cellulose-water interface sufficiently well. Furthermore, a difference in T1 of about a factor of 2 in the range of Larmor frequencies 25-150 MHz was found for C4 atoms in chains located on top of two different crystallographic planes, namely, (110) and (10). A previously proposed explanation that the C4 peak doublet could derive from surfaces parallel to different crystallographic planes is herewith strengthened by computationally obtained evidence. Another suggested basis for this difference is that the doublet originates from C4 atoms located in surface anhydro-glucose units with hydroxymethyl groups pointing either inward or outward. This was also tested within this study but was found to yield no difference in calculated T1.

Molecular dynamics simulations of membrane-sugar interactions.

It is well documented that disaccharides in general and trehalose (TRH) in particular strongly affect physical properties and functionality of lipid bilayers. We investigate interactions between lipid membranes formed by 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) and TRH by means of molecular dynamics (MD) computer simulations. Ten different TRH concentrations were studied in the range wTRH = 0-0.20 (w/w). The potential of mean force (PMF) for DMPC bilayer-TRH interactions was determined using two different force fields, and was subsequently used in a simple analytical model for description of sugar binding at the membrane interface. The MD results were in good agreement with the predictions of the model. The net affinities of TRH for the DMPC bilayer derived from the model and MD simulations were compared with experimental results. The area per lipid increases and the membrane becomes thinner with increased TRH concentration, which is interpreted as an intercalation effect of the TRH molecules into the polar part of the lipids, resulting in conformational changes in the chains. These results are consistent with recent experimental observations. The compressibility modulus related to the fluctuations of the membrane increases dramatically with increased TRH concentration, which indicates higher order and rigidity of the bilayer. This is also reflected in a decrease (by a factor of 15) of the lateral diffusion of the lipids. We interpret these observations as a formation of a glassy state at the interface of the membrane, which has been suggested in the literature as a hypothesis for the membrane-sugar interactions.

Free energy surfaces for the interaction of D-glucose with planar aromatic groups in aqueous solution

Multidimensional potentials of mean force for the interactions in aqueous solution of both anomers of D-glucopyranose with two planar aromatic molecules, indole and para-methyl-phenol, have been calculated using molecular dynamics simulations with umbrella sampling and were subsequently used to estimate binding free energies. Indole and para-methyl-phenol serve as models for the side chains of the amino acids tryptophan and tyrosine, respectively. In all cases, a weak affinity between the glucose molecules and the flat aromatic surfaces was found. The global minimum for these interactions was found to be for the case when the pseudoplanar face of β-D-glucopyranose is stacked against the planar surfaces of the aromatic residues. The calculated binding free energies are in good agreement with both experiment and previous simulations. The multidimensional free energy maps suggest a mechanism that could lend kinetic stability to the complexes formed by sugars bound to sugar-binding proteins.

Molecular Dynamics Simulations of the Ionic Liquid 1-n-Butyl-3-Methylimidazolium Chloride and Its Binary Mixtures with Ethanol.

Room temperature ionic liquids (ILs) of the imidazolium family have attracted much attention during the past decade for their capability to dissolve biomass. Besides experimental work, numerous compuational studies have been concerned with the physical properties of both neat ILs and their interactions with different solutes, in particular, carbohydrates. Many classical force fields designed specifically for ILs have been found to yield viscosities that are too high for the liquid state, which has been attributed to the fact that the effective charge densities are too high due to the lack of electronic polarizability. One solution to this problem has been uniform scaling of the partial charges by a scale factor in the range 0.6-0.9, depending on model. This procedure has been shown to improve the viscosity of the models, and also to positively affect other properties, such as diffusion constants and ionic conductivity. However, less attention has been paid to how this affects the overall thermodynamics of the system, and the problems it might create when the IL models are combined with other force fields (e.g., for solutes). In the present work, we employ three widely used IL force fields to simulate 1-n-butyl-3-methyl-imidazolium chloride in both the crystal and the liquid state, as well as its binary mixture with ethanol. Two approaches are used: one in which the ionic charge is retained at its full integer value and one in which the partial charges are uniformly reduced to 85%. We investigate and calculate crystal and liquid structures, molar heat capacities, heats of fusion, self-diffusion constants, ionic conductivity, and viscosity for the neat IL, and ethanol activity as a function of ethanol concentration for the binary mixture. We show that properties of the crystal are less affected by charge scaling compared to the liquid. In the liquid state, transport properties of the neat IL are generally improved by scaling, whereas values for the heat of fusion are unaffected, and results for the heat capacity are ambiguous. Neither full nor reduced charges could reproduce experimental ethanol activities for the whole range of compositions.

Molecular dynamics simulation of strong interaction mechanisms at wet interfaces in clay–polysaccharide nanocomposites

Bio-composites comprised of the polysaccharide xyloglucan (XG) and montmorillonite (MTM) clay has potential as a ‘green’ replacement of conventional petroleum-derived polymers in the packaging industry. These materials have been shown to possess excellent material properties, even in high relative humidity. Although interfacial interaction between XG and MTM, and the molecular structure of XG can be identified as key parameters for the complex formation process and the resulting tensile properties, these properties are usually difficult to address using experimental methods. Here we use molecular dynamics (MD) simulations to study the adsorption of fully atomistic models of both native and chemically modified XG to MTM clay surfaces in explicit water. We show that the driving force for adsorption is a favorable change in enthalpy, and furthermore that native XG adsorbs stronger than modified XG. This highlights the importance of molecular structure details to molecular adhesion. The present study provides insights into the molecular scale adsorption mechanisms and can therefore help in designing routes for further improvements of bio-composite materials.

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.