Estimation of the Persistence Length of Polymers by MD Simulations on Small Fragments in Solution. Application to Cellulose

The Conformation of Glycosidic Linkages According to Various Force Fields: Monte Carlo Modeling of Polysaccharides Based on Extrapolation of Short-Chain Properties

The conformational features of the glycosidic linkage are the most important variable to consider when studying di-, oligo-, and polysaccharide molecules using molecular dynamics (MD) simulations. The accuracy of the theoretical model describing this degree of freedom influences the quality of the results obtained from MD calculations based on this model. This article focuses on the following two issues related to the conformation of the glycosidic linkage. First, we describe the results of a comparative analysis of the predictions of three carbohydrate-dedicated classical force fields for MD simulations, namely, CHARMM, GLYCAM, and GROMOS, in the context of different parameters of structural and energetic nature related to the conformation of selected types of glycosidic linkages, α(1 → 4), β(1 → 3), and β(1 → 4), connecting glucopyranose units. This analysis revealed several differences, mainly concerning the energy levels of the secondary and tertiary conformers and the linkage flexibility within the dominant exo-syn conformation for α(1 → 4) and β(1 → 3) linkages. Some aspects of the comparative analysis also included the newly developed, carbohydrate-dedicated Martini 3 coarse-grained force field. Second, to overcome the time-scale problem associated with sampling slow degrees of freedom in polysaccharide chains during MD simulations, we developed a coarse-grained (CG) model based on the data from MD simulations and designed for Monte Carlo modeling. This model (CG MC) is based on information from simulations of short saccharide chains, effectively sampled in atomistic MD simulations, and is capable of extrapolating local conformational properties to the case of polysaccharides of arbitrary length. The CG MC model has the potential to estimate the conformations of very long polysaccharide chains, taking into account the influence of secondary and tertiary conformations of glycosidic linkages. With respect to the comparative analysis of force fields, the application of CG MC modeling showed that relatively small differences in the predictions of individual force fields with respect to a single glycosidic linkage accumulate when considering their effect on the structure of longer chains, leading to drastically different predictions with respect to parameters describing the polymer conformation, such as the persistence length.

Prediction and Control in DNA Nanotechnology

DNA nanotechnology is a rapidly developing field that uses DNA as a building material for nanoscale structures. Key to the field's development has been the ability to accurately describe the behavior of DNA nanostructures using simulations and other modeling techniques. In this Review, we present various aspects of prediction and control in DNA nanotechnology, including the various scales of molecular simulation, statistical mechanics, kinetic modeling, continuum mechanics, and other prediction methods. We also address the current uses of artificial intelligence and machine learning in DNA nanotechnology. We discuss how experiments and modeling are synergistically combined to provide control over device behavior, allowing scientists to design molecular structures and dynamic devices with confidence that they will function as intended. Finally, we identify processes and scenarios where DNA nanotechnology lacks sufficient prediction ability and suggest possible solutions to these weak areas.

Penetration of Cell Surface Glycocalyx by Enveloped Viruses Is Aided by Weak Multivalent Adhesive Interaction

Viral infection usually begins with adhesion between the viral particle and viral receptors displayed on the cell membrane. The exterior surface of the cell membrane is typically coated with a brush-like layer of molecules, the glycocalyx, that the viruses need to penetrate. Although there is extensive literature on the biomechanics of virus-cell adhesion, much of it is based on continuum-level models that do not address the question of how virus/cell-membrane adhesion occurs through the glycocalyx. In this work, we present a simulation study of the penetration mechanism. Using a coarse-grained molecular model, we study the force-driven and diffusive penetration of a brush-like glycocalyx by viral particles. For force-driven penetration, we find that viral particles smaller than the spacing of molecules in the brush reach the membrane surface readily. For a given maximum force, viral particles larger than the minimum spacing of brush molecules arrest at some distance from the membrane, governed by the balance of elastic and applied forces. For the diffusive case, we find that weak but multivalent attraction between the glycocalyx molecules and the virus effectively leads to its engulfment by the glycocalyx. Our finding provides potential guidance for developing glycocalyx-targeting drugs and therapies by understanding how virus-cell adhesion works.

Single-molecule investigations of single-chain cellulose biosynthesis

Cellulose biosynthesis in sessile bacterial colonies originates in the membrane-integrated bacterial cellulose synthase (Bcs) AB complex. We utilize optical tweezers to measure single-strand cellulose biosynthesis by BcsAB from Rhodobacter sphaeroides. Synthesis depends on uridine diphosphate glucose, Mg²⁺, and cyclic diguanosine monophosphate, with the last displaying a retention time of ∼80 min. Below a stall force of 12.7 pN, biosynthesis is relatively insensitive to force and proceeds at a rate of one glucose addition every 2.5 s at room temperature, increasing to two additions per second at 37°. At low forces, conformational hopping is observed. Single-strand cellulose stretching unveiled a persistence length of 6.2 nm, an axial stiffness of 40.7 pN, and an ability for complexes to maintain a tight grip, with forces nearing 100 pN. Stretching experiments exhibited hysteresis, suggesting that cellulose microstructure underpinning robust biofilms begins to form during synthesis. Cellohexaose spontaneously binds to nascent single cellulose strands, impacting polymer mechanical properties and increasing BcsAB activity.

Hydrodynamic and conformational characterization of aqueous sodium alginate solutions with varying salinity

Insight into the role of electrostatic interactions on the hydrodynamics and conformation of aqueous sodium alginate was gained through viscometry. Alginate chains are found to shrink in salt-free solutions more rapidly with increasing polymer concentration compared to salt-solutions. For salt-free solutions, a reduced polymer concentration of less than 1 suffices to make the alginate coil volume half of that at infinite dilution which becomes invariant when the reduced concentration exceeds 8. In saline media having salt concentration greater than 0.1 mol·L^-1, the chains become more flexible, caused by the shielding of intra-chain repulsions. The chains effectively reached unperturbed state when the added salt concentration becomes ≥0.5 mol·L^-1. Alginate chains are shown to remain stiff up to about 8-10 monomers within the investigated temperature range. This study explores the possible modification of the individual chain behavior induced by the neighboring chains or by the variation of temperature.

Sodium ion transport across the endothelial glycocalyx layer under electric field conditions: A molecular dynamics study

In the present research, the sodium ion transport across the endothelial glycocalyx layer (EGL) under an imposed electric field is investigated, for the first time, using a series of molecular dynamics simulations. The electric field is perpendicularly imposed on the EGL with varying strengths. The sodium ion molarity difference between the inner and outer layers of EGL, Δc, is used to quantify the sodium transport in the presence of the negatively charged glycocalyx sugar chains. Results suggest that a weak electric field increases Δc, regardless of whether the electric field is imposed perpendicularly inward or outward. By contrast, a strong electric field drives sodium ions to travel in the same orientation as the electric field. Scrutiny of the charge distribution of the glycocalyx sugar chains suggests that the electric field modifies the spatial layouts of glycocalyx atoms as it drives the transport of sodium ions. The modification in glycocalyx layouts further changes the inter-molecular interactions between glycocalyx sugar chains and sodium ions, thereby limiting the electric field control of ion transport. The sodium ions, in turn, alter the apparent bending stiffness of glycocalyx. Moreover, the negative charges of the glycocalyx sugar chains play an important role in maintaining structural stability of endothelial glycocalyx. Based on the findings, a hypothesis is proposed regarding the existence of a strength threshold of the electric field in controlling charged particles in the endothelium, which offers an alternative explanation for contrasting results in previous experimental observations.

Revisiting the Dissolution of Cellulose in NaOH as "Seen" by X-rays

Cotton production is reaching a global limit, leading to a growing demand for bio-based textile fibers produced by other means. Textile fibers based on regenerated cellulose from wood holds great potential, but in order to produce fibers, the components need to be dissolved in suitable solvents. Furthermore, the dissolution process of cellulose is not yet fully understood. In this study, we investigated the dissolution state of microcrystalline cellulose in aqueous NaOH by using primarily scattering methods. Contrary to previous findings, this study indicated that cellulose concentrations of up to 2 wt % are completely molecularly dissolved in 8 wt % NaOH. Scattering data furthermore revealed the presence of semi-flexible cylinders with stiff segments. In order to improve the dissolution capability of NaOH, the effects of different additives have been of interest. In this study, scattering data indicated that the addition of ZnO decreased the formation of aggregates, while the addition of PEG did not improve the dissolution properties significantly, although preliminary NMR data did suggest a weak attraction between PEG and cellulose. Overall, this study sheds further light on the dissolution of cellulose in NaOH and highlights the use of scattering methods to assess solvent quality.

Enhanced Antimicrobial Activity and Structural Transitions of a Nanofibrillated Cellulose-Nisin Biocomposite Suspension

Resistance to antibiotics has posed a high demand for novel strategies to fight bacterial infections. Antimicrobial peptides (AMPs) are a promising alternative to conventional antibiotics. However, their poor solubility in water and sensitivity to degradation has limited their application. Here, we report the design of a smart, pH-responsive antimicrobial nanobiocomposite material based on the AMP nisin and 2,2,6,6-tetramethyl-1-piperidinyloxyl-oxidized nanofibrillated cellulose (TONFC). Morphological transformations of the nanoscale structure of nisin functionalized-TONFC fibrils were discovered at pH values between 5.8 and 8.0 using small-angle X-ray scattering. Complementary ζ potential measurements indicate that electrostatic attractions between the negatively charged TONFC surface and the positively charged nisin molecules are responsible for the integration of nisin. Modification of the pH level or increasing the ionic strength reduces the nisin binding capacity of TONFC. Biological evaluation studies using a bioluminescence-based reporter strain of Bacillus subtilis and a clinically relevant strain of Staphylococcus aureus indicated a significantly higher antimicrobial activity of the TONFC-nisin biocomposite compared to the pure nisin against both strains under physiological pH and ionic strength conditions. The in-depth characterization of this new class of antimicrobial biocomposite material based on nanocellulose and nisin may guide the rational design of sustainable antimicrobial materials.

Molecular dynamics study of the LCST transition in aqueous poly(N-n-propylacrylamide)

The breadth of technological applications of smart polymers relies on the possibility of tuning their molecular structure to respond to external stimuli. In this context, N-substituted acrylamide-based polymers are widely studied thermoresponsive polymers. Poly(N-n-propylacrylamide) (PNnPAm), which is a structural isomer of the poly(N-isopropylacrylamide) (PNIPAm) exhibits however, a lower phase transition in aqueous solution. In this work, we use all-atom molecular dynamics simulations of PNnPAm in aqueous solutions to study, from a microscopic point-of-view, the influence of chain size and concentration on the LCST of PNnPAm. Our analysis shows that the collapse of a single oligomer of PNnPAm upon heating is dependent on the chain length and corresponds to a complex interplay between hydration and intermolecular interactions. Analysis of systems with multiple chains shows an aggregation of PNnPAm chains above the LCST.

A Facile Method to Fabricate Anisotropic Hydrogels with Perfectly Aligned Hierarchical Fibrous Structures

Natural structural materials (such as tendons and ligaments) are comprised of multiscale hierarchical architectures, with dimensions ranging from nano- to macroscale, which are difficult to mimic synthetically. Here a bioinspired, facile method to fabricate anisotropic hydrogels with perfectly aligned multiscale hierarchical fibrous structures similar to those of tendons and ligaments is reported. The method includes drying a diluted physical hydrogel in air by confining its length direction. During this process, sufficiently high tensile stress is built along the length direction to align the polymer chains and multiscale fibrous structures (from nano- to submicro- to microscale) are spontaneously formed in the bulk material, which are well-retained in the reswollen gel. The method is useful for relatively rigid polymers (such as alginate and cellulose), which are susceptible to mechanical signal. By controlling the drying with or without prestretching, the degree of alignment, size of superstructures, and the strength of supramolecular interactions can be tuned, which sensitively influence the strength and toughness of the hydrogels. The mechanical properties are comparable with those of natural ligaments. This study provides a general strategy for designing hydrogels with highly ordered hierarchical structures, which opens routes for the development of many functional biomimetic materials for biomedical applications.

Modeling of the morphological change of cellulose microfibrils caused with aqueous NaOH solution: the longitudinal contraction and laterally swelling during decrystallization

The conformation of cellulose microfibrils treated with aqueous NaOH was modeled as partially decrystallized cellulose chains before completing conversion to cellulose II, in order to elucidate the change in morphology of ramie fiber caused by NaOH treatment. Equations for the relative length and width of the microfibrils were derived on the basis of partially decrystallized microfibrils modeling. Each equation contains four parameters, n, β, w _c , and c _r , which correspond to the number of glucose residues between periodic defects along the untreated ramie cellulose microfibrils, the extension ratio of amorphous cellulose chain along length, the cross-section crystallinity, and the correction term of crystallinity, respectively. The validity of the derived equations was confirmed by two types of simulations. One is performed using experimental data L/L ₀ and W/W ₀ as a function of crystallinity, while the other is done using the relationship between the relative length and width obtained from the experimental data, which is independent of crystallinity, was performed. The best-fit simulation was obtained under n = 277, β = 2.813, and c _r w _c  = 0.671 for the former and under n = 301 and β = 2.792 for the latter. These values of n and β correspond closely to the values reported in references for ramie microfibrils. Both simulation results show that macroscopic changes in the morphology of ramie fibers is attributable to the changes in cellulose chain conformation in the decrystallized regions created along the microfibrils upon NaOH treatment.

Nanoscale Concentration Quantification of Pharmaceutical Actives in Amorphous Polymer Matrices by Electron Energy-Loss Spectroscopy

We demonstrated the use of electron energy-loss spectroscopy (EELS) to evaluate the composition of phenytoin:hydroxypropyl methylcellulose acetate succinate (HPMCAS) spin-coated solid dispersions (SDs). To overcome the inability of bright-field and high-angle annular dark-field TEM imaging to distinguish between glassy drug and polymer, we used the π-π* transition peak in the EELS spectrum to detect phenytoin within the HPMCAS matrix of the SD. The concentration of phenytoin within SDs of 10, 25, and 50 wt % drug loading was quantified by a multiple least-squares analysis. Evaluating the concentration of 50 different regions in each SD, we determined that phenytoin and HPMCAS are intimately mixed at a length scale of 200 nm, even for drug loadings up to 50 wt %. At length scales below 100 nm, the variance of the measured phenytoin concentration increases; we speculate that this increase is due to statistical fluctuations in local concentration and chemical changes induced by electron irradiation. We also performed EELS analysis of an annealed 25 wt % phenytoin SD and showed that the technique can resolve concentration differences between regions that are less than 50 nm apart. Our findings indicate that EELS is a useful tool for quantifying, with high accuracy and sub-100 nm spatial resolution, the composition of many pharmaceutical and soft matter systems.

Understanding nanocellulose chirality and structure-properties relationship at the single fibril level

Nanocellulose fibrils are ubiquitous in nature and nanotechnologies but their mesoscopic structural assembly is not yet fully understood. Here we study the structural features of rod-like cellulose nanoparticles on a single particle level, by applying statistical polymer physics concepts on electron and atomic force microscopy images, and we assess their physical properties via quantitative nanomechanical mapping. We show evidence of right-handed chirality, observed on both bundles and on single fibrils. Statistical analysis of contours from microscopy images shows a non-Gaussian kink angle distribution. This is inconsistent with a structure consisting of alternating amorphous and crystalline domains along the contour and supports process-induced kink formation. The intrinsic mechanical properties of nanocellulose are extracted from nanoindentation and persistence length method for transversal and longitudinal directions, respectively. The structural analysis is pushed to the level of single cellulose polymer chains, and their smallest associated unit with a proposed 2 × 2 chain-packing arrangement.

Cellulose is a material found in many different biological systems, but the fine structure at the single-molecule level is still being assessed. Here, the authors present high-resolution imaging of cellulose structures at the single particle level, finding evidence of chirality in bundles and fibrils.

Conformational analysis of oligosaccharides and polysaccharides using molecular dynamics simulations

Complex carbohydrates usually have a large number of rotatable bonds and consequently a large number of theoretically possible conformations can be generated (combinatorial explosion). The application of systematic search methods for conformational analysis of carbohydrates is therefore limited to disaccharides and trisaccharides in a routine analysis. An alternative approach is to use Monte-Carlo methods or (high-temperature) molecular dynamics (MD) simulations to explore the conformational space of complex carbohydrates. This chapter describes how to use MD simulation data to perform a conformational analysis (conformational maps, hydrogen bonds) of oligosaccharides and how to build realistic 3D structures of large polysaccharides using Conformational Analysis Tools (CAT).

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

Ionic liquids have become a popular solvent for cellulose pretreatment in biorefineries due to their efficiency in dissolution and their reusability. Understanding the interactions between cations, anions, and cellulose is key to the development of better solvents and the improvement of pretreatment conditions. While previous studies described the interactions between ionic liquids and cellulose fibers, shedding light on the initial stages of the cellulose dissolution process, we study the end state of that process by exploring the structure and dynamics of a single cellulose decamer solvated in 1-butyl-3-methyl-imidazolium chloride (BmimCl) and in water using replica-exchange molecular dynamics. In both solvents, global structural features of the cellulose chain are similar. However, analyses of local structural properties show that cellulose explores greater conformational variability in the ionic liquid than in water. For instance, in BmimCl the cellulose intramolecular hydrogen bond O3H'···O5 is disrupted more often resulting in greater flexibility of the solute. Our results indicate that the cellulose chain is more dynamic in BmimCl than in water, which may play a role in the favorable dissolution of cellulose in the ionic liquid. Calculation of the configurational entropy of the cellulose decamer confirms its higher conformational flexibility in BmimCl than in water at elevated temperatures.

Structure and response to flow of the glycocalyx layer

The glycocalyx is a sugar-rich layer located at the luminal part of the endothelial cells. It is involved in key metabolic processes and its malfunction is related to several diseases. To understand the function of the glycocalyx, a molecular level characterization is necessary. In this article, we present large-scale molecular-dynamics simulations that provide a comprehensive description of the structure and dynamics of the glycocalyx. We introduce the most detailed, to-date, all-atom glycocalyx model, composed of lipid bilayer, proteoglycan dimers, and heparan sulfate chains with realistic sequences. Our results reveal the folding of proteoglycan ectodomain and the extended conformation of heparan sulfate chains. Furthermore, we study the glycocalyx response under shear flow and its role as a flypaper for binding fibroblast growth factors (FGFs), which are involved in diverse functions related to cellular differentiation, including angiogenesis, morphogenesis, and wound healing. The simulations show that the glycocalyx increases the effective concentration of FGFs, leading to FGF oligomerization, and acts as a lever to transfer mechanical stimulus into the cytoplasmic side of endothelial cells.

Escherichia coli peptidoglycan structure and mechanics as predicted by atomic-scale simulations

Bacteria face the challenging requirement to maintain their shape and avoid rupture due to the high internal turgor pressure, but simultaneously permit the import and export of nutrients, chemical signals, and virulence factors. The bacterial cell wall, a mesh-like structure composed of cross-linked strands of peptidoglycan, fulfills both needs by being semi-rigid, yet sufficiently porous to allow diffusion through it. How the mechanical properties of the cell wall are determined by the molecular features and the spatial arrangement of the relatively thin strands in the larger cellular-scale structure is not known. To examine this issue, we have developed and simulated atomic-scale models of Escherichia coli cell walls in a disordered circumferential arrangement. The cell-wall models are found to possess an anisotropic elasticity, as known experimentally, arising from the orthogonal orientation of the glycan strands and of the peptide cross-links. Other features such as thickness, pore size, and disorder are also found to generally agree with experiments, further supporting the disordered circumferential model of peptidoglycan. The validated constructs illustrate how mesoscopic structure and behavior emerge naturally from the underlying atomic-scale properties and, furthermore, demonstrate the ability of all-atom simulations to reproduce a range of macroscopic observables for extended polymer meshes.

The structure of the bacterial cell wall has been a point of controversy and contention since it was first discovered. Although the basic chemical composition of peptidoglycan, the key constituent of the cell wall, is now well established, its long-range organization is not. This dearth of information at the mesoscopic scale is a result of the inability of experimental imaging techniques to simultaneously visualize both the atomic-level detail of the peptidoglycan network and its macroscopic arrangement around the bacterium. Now, using molecular dynamics (MD) simulations, we have carefully constructed and validated models of sections of the Escherichia coli cell wall in full atomic detail. By comparing various properties of these models, including elasticity, pore size, and thickness with experiments, we can discriminate between them, resolving which best represents the native wall structure. In doing so, our study provides approaches for connecting measurements made in atomic-scale MD simulations with large-scale and even macroscopic properties.

Marine polysaccharides in microencapsulation and application to aquaculture: "from sea to sea"

This review's main objective is to discuss some physico-chemical features of polysaccharides as intrinsic determinants for the supramolecular structures that can efficiently provide encapsulation of drugs and other biological entities. Thus, the general characteristics of some basic polysaccharides are outlined in terms of their conformational, dynamic and thermodynamic properties. The analysis of some polysaccharide gelling properties is also provided, including the peculiarity of the charged polysaccharides. Then, the way the basic physical chemistry of polymer self-assembly is made in practice through the laboratory methods is highlighted. A description of the several literature procedures used to influence molecular interactions into the macroscopic goal of the encapsulation is given with an attempt at classification. Finally, a practical case study of specific interest, the use of marine polysaccharide matrices for encapsulation of vaccines in aquaculture, is reported.

Cellulose nanomaterials review: structure, properties and nanocomposites

This critical review provides a processing-structure-property perspective on recent advances in cellulose nanoparticles and composites produced from them. It summarizes cellulose nanoparticles in terms of particle morphology, crystal structure, and properties. Also described are the self-assembly and rheological properties of cellulose nanoparticle suspensions. The methodology of composite processing and resulting properties are fully covered, with an emphasis on neat and high fraction cellulose composites. Additionally, advances in predictive modeling from molecular dynamic simulations of crystalline cellulose to the continuum modeling of composites made with such particles are reviewed (392 references).

A reoptimized GROMOS force field for hexopyranose-based carbohydrates accounting for the relative free energies of ring conformers, anomers, epimers, hydroxymethyl rotamers, and glycosidic linkage conformers

This article presents a reoptimization of the GROMOS 53A6 force field for hexopyranose-based carbohydrates (nearly equivalent to 45A4 for pure carbohydrate systems) into a new version 56A(CARBO) (nearly equivalent to 53A6 for non-carbohydrate systems). This reoptimization was found necessary to repair a number of shortcomings of the 53A6 (45A4) parameter set and to extend the scope of the force field to properties that had not been included previously into the parameterization procedure. The new 56A(CARBO) force field is characterized by: (i) the formulation of systematic build-up rules for the automatic generation of force-field topologies over a large class of compounds including (but not restricted to) unfunctionalized polyhexopyranoses with arbritrary connectivities; (ii) the systematic use of enhanced sampling methods for inclusion of experimental thermodynamic data concerning slow or unphysical processes into the parameterization procedure; and (iii) an extensive validation against available experimental data in solution and, to a limited extent, theoretical (quantum-mechanical) data in the gas phase. At present, the 56A(CARBO) force field is restricted to compounds of the elements C, O, and H presenting single bonds only, no oxygen functions other than alcohol, ether, hemiacetal, or acetal, and no cyclic segments other than six-membered rings (separated by at least one intermediate atom). After calibration, this force field is shown to reproduce well the relative free energies of ring conformers, anomers, epimers, hydroxymethyl rotamers, and glycosidic linkage conformers. As a result, the 56A(CARBO) force field should be suitable for: (i) the characterization of the dynamics of pyranose ring conformational transitions (in simulations on the microsecond timescale); (ii) the investigation of systems where alternative ring conformations become significantly populated; (iii) the investigation of anomerization or epimerization in terms of free-energy differences; and (iv) the design of simulation approaches accelerating the anomerization process along an unphysical pathway.

Conformational analysis in aqueous solution and estimation of the persistence length of exopolysaccharides produced by Lactobacillus helveticus Lh59 and Streptococcus macedonicus Sc136

The conformations of two exopolysaccharides, one produced by Lactobacillus helveticus Lh59 and one by Streptococcus macedonicus Sc136, were investigated by molecular modelling to get insight into their physical properties. In particular, the influence of side chains and glycosidic linkage types on the overall shapes and persistence lengths of the polysaccharides were studied. It appeared that the side chains had only minor effects on the persistence lengths of the polysaccharides, with the exception of the monosaccharide residue directly attached to the backbone. The extensiveness of the backbone is enhanced by β-D-Hexp-(1→4) structural elements, whereas α-D-Hexp-(1→3) and α-D-Hexp-(1→4) elements create structural variability by introducing bends. The occurrence of either (1→5)- or (1→6)-linked monosaccharides, resulting in one extra bond between monosaccharide units, causes the overall flexibility of the polysaccharide chain to be enhanced and extensiveness to be reduced. The presence of such flexible linkages in the backbone leads to a shorter persistence length than when present in the side chain: a value of only 4.5 nm was found for the L. helveticus Lh59 EPS with a (1→5) linkage in the backbone, and a value of 8.5 nm for the S. macedonicus Sc136 EPS with two (1→6) linkages in the side chain.Key words: conformational analysis, exopolysaccharide, Lactobacillus helveticus, molecular modelling, persistence length, Streptococcus macedonicus.

Ultrasonic wave spectroscopy study of sugar oligomers and polysaccharides in aqueous solutions: The hydration length concept

The determination of apparent persistence length and radius of gyration of maltodextrins in water is achievable through high-resolution ultrasonic spectroscopy measurements. Classical hydration number for those carbohydrates is characteristic of an apparent persistence degree of polymerisation of the polymer. A force-field based molecular modeling of a 10DP malto-oligomer allows measurement of the corresponding length for the lowest energetic conformation in solution. A good agreement between the apparent radii of gyration determined by this technique and the freely rotating polymer chain model is found with radii of gyration calculated from the intrinsic viscosity.

A new GROMOS force field for hexopyranose-based carbohydrates

A new parameter set (referred to as 45A4) is developed for the explicit-solvent simulation of hexopyranose-based carbohydrates. This set is compatible with the most recent version of the GROMOS force field for proteins, nucleic acids, and lipids, and the SPC water model. The parametrization procedure relies on: (1) reassigning the atomic partial charges based on a fit to the quantum-mechanical electrostatic potential around a trisaccharide; (2) refining the torsional potential parameters associated with the rotations of the hydroxymethyl, hydroxyl, and anomeric alkoxy groups by fitting to corresponding quantum-mechanical profiles for hexopyranosides; (3) adapting the torsional potential parameters determining the ring conformation so as to stabilize the (experimentally predominant) (4)C(1) chair conformation. The other (van der Waals and nontorsional covalent) parameters and the rules for third and excluded neighbors are taken directly from the most recent version of the GROMOS force field (except for one additional exclusion). The new set is general enough to define parameters for any (unbranched) hexopyranose-based mono-, di-, oligo- or polysaccharide. In the present article, this force field is validated for a limited set of monosaccharides (alpha- and beta-D-glucose, alpha- and beta-D-galactose) and disaccharides (trehalose, maltose, and cellobiose) in solution, by comparing the results of simulations to available experimental data. More extensive validation will be the scope of a forthcoming article. (c) 2005 Wiley Periodicals, Inc. J Comput Chem 26: 1400-1412, 2005.

Phase behavior of cationic hydroxyethyl cellulose-sodium dodecyl sulfate mixtures: effects of molecular weight and ethylene oxide side chain length of polymers

Novel cationic hydroxyethyl cellulose (HEC) polymers with different molecular weights (1.1 x 10(5) to 1.7 x 10(6) g/mol) and ethylene oxide (EO) side chain lengths (1.5-2.9 EO units) were mixed with sodium dodecyl sulfate (SDS) in aqueous solutions. The phase diagrams of cationic HEC-SDS complexes were determined in the dilute polymer concentration regime (< 0.5 wt %) with gradual addition of SDS molecules. The viscosity and structures of the complexes during the phase evolution were studied using rheometry and dynamic light scattering. The gradual addition of SDS first induced interchain associations with the bound SDS aggregates serving as cross-linkers to form an open network structure, producing a very broad size distribution and high viscosities of the complex solutions, and then condensed the network and induced a structure reorganization, resulting in globular aggregates with narrow size distributions. The growth of these globular aggregates in size eventually led to macroscopic sedimentation near charge neutralization. Further addition of SDS randomly broke the sedimentary aggregates into small particles and SDS micelles with low solution viscosities. The effects of molecular weight and EO side chain length of polymers on the phase boundary, viscosity, and structure of cationic HEC-SDS complexes were discussed.

Molecular dynamics simulations of alpha2 --> 8-linked disialoside: conformational analysis and implications for binding to proteins

Computational methods have played a key role in elucidating the various three-dimensional structures of oligosaccharides. Such structural information, together with other experimental data, leads to a better understanding of the role of oligosaccharide in various biological processes. The disialoside Neu5Ac-alpha2-->8-Neu5Ac appears as the terminal glycan in glycoproteins and glycolipids, and is known to play an important role in various events of cellular communication. Neurotoxins such as botulinum and tetanus require Neu5Ac-alpha2 --> 8-Neu5Ac for infecting the host. Glycoconjugates containing this disialoside and the enzymes catalyzing their biosynthesis are also regulated during cell growth, development, and differentiation. Unlike other biologically relevant disaccharides that have only two linkage bonds, the alpha2-->8-linked disialoside has four: C2-O, O-C8', C8'-C7', and C7'-C6'. The present report describes the results from nine 1 ns MD simulations of alpha2-->8-linked disialoside (Neu5Ac-alpha2-->8-Neu5Ac); simulations were run using GROMOS96 by explicitly considering the solvent molecules. Conformations around the O-C8' bond are restricted to the +sc/+ap regions due to stereochemical reasons. In contrast, conformations around the C2-O and C8'-C7' bonds were found to be largely unrestricted and all the three staggered regions are accessible. The conformations around the C7'-C6' bond were found to be in either the -sc or the anti region. These results are in excellent agreement with the available NMR and potential energy calculation studies. Overall, the disaccharide is flexible and adopts mainly two ensembles of conformations differing in the conformation around the C7'-C6' bond. The flexibility associated with this disaccharide allows for better optimization of intermolecular contacts while binding to proteins and this may partially compensate for the loss of conformational entropy that may be incurred due to disaccharide's flexibility.