A coarse-grained molecular model for glycosaminoglycans: application to chondroitin, chondroitin sulfate, and hyaluronic acid

Coarse-Grained Model of Glycosaminoglycans for Biomolecular Simulations

Proteoglycans contain glycosaminoglycans (GAGs) which are negatively charged linear polymers made of repeating disaccharide units of uronic acid and hexosamine units. They play vital roles in numerous physiological and pathological processes, particularly in governing cellular communication and attachment. Depending on their sulfonation state, acetylation, and glycosidic linkages, GAGs belong to different families. The high molecular weight, heterogeneity, and flexibility of GAGs hamper their characterization at atomic resolution, but this may be circumvented via coarse-grained (CG) approaches. In this work, we report a CG model for a library of common GAG types in their isolated or proteoglycan-linked states compatible with version 2.2 (v2.2) of the widely popular CG Martini force field. The model reproduces conformational and thermodynamic properties for a wide variety of GAGs, as well as matching structural and binding data for selected proteoglycan test systems. The parameters developed here may thus be employed to study a range of GAG-containing biomolecular systems, thereby benefiting from the efficiency and broad applicability of the Martini framework.

Chondroitin sulfate, dermatan sulfate, and hyaluronic acid differentially modify the biophysical properties of collagen-based hydrogels

Fibrillar collagens and glycosaminoglycans (GAGs) are structural biomolecules that are natively abundant to the extracellular matrix (ECM). Prior studies have quantified the effects of GAGs on the bulk mechanical properties of the ECM. However, there remains a lack of experimental studies on how GAGs alter other biophysical properties of the ECM, including ones that operate at the length scales of individual cells such as mass transport efficiency and matrix microstructure. This study focuses on the GAG molecules chondroitin sulfate (CS), dermatan sulfate (DS), and hyaluronic acid (HA). CS and DS are stereoisomers while HA is the only non-sulfated GAG. We characterized and decoupled the effects of these GAG molecules on the stiffness, transport, and matrix microarchitecture properties of type I collagen hydrogels using mechanical indentation testing, microfluidics, and confocal reflectance imaging, respectively. We complement these biophysical measurements with turbidity assays to profile collagen aggregate formation. Surprisingly, only HA enhanced the ECM indentation modulus, while all three GAGs had no effect on hydraulic permeability. Strikingly, we show that CS, DS, and HA differentially regulate the matrix microarchitecture of hydrogels due to their alterations to the kinetics of collagen self-assembly. In addition to providing information on how GAGs define key physical properties of the ECM, this work shows new ways in which stiffness measurements, microfluidics, microscopy, and turbidity kinetics can be used complementarily to reveal details of collagen self-assembly and structure. STATEMENT OF SIGNIFICANCE: Collagen and glycosaminoglycans (GAGs) are integral to the structure, function, and bioactivity of the extracellular matrix (ECM). Despite widespread interest in collagen-GAG composite hydrogels, there is a lack of quantitative understanding of how different GAGs alter the biophysical properties of the ECM across tissue, cellular, and subcellular length scales. Here we show using mechanical, microfluidic, microscopy, and analytical methods and measurements that the GAG molecules chondroitin sulfate, dermatan sulfate, and hyaluronic acid differentially regulate the mechanical, transport, and microstructural properties of hydrogels due to their alterations to the kinetics of collagen self-assembly. As such, these results will inform improved design and utilization of collagen-based scaffolds of tailored composition, mechanical properties, molecular availability due to mass transport, and microarchitecture.

Interpenetrating Networks of Collagen and Hyaluronic Acid That Serve as In Vitro Tissue Models for Assessing Macromolecular Transport

High-fidelity preclinical in vitro tissue models can reduce the failure rate of drugs entering clinical trials. Collagen and hyaluronic acid (HA) are major components of the extracellular matrix of many native tissues and affect therapeutic macromolecule diffusion and recovery through tissues. Although collagen and HA are commonly used in tissue engineering, the physical and mechanical properties of these materials are variable and depend highly on processing conditions. In this study, HA was chemically modified and crosslinked via hydrazone bonds to form interpenetrating networks of crosslinked HA (HAX) with collagen (Col). These networks enabled a wide range of mechanical properties, including stiffness and swellability, and microstructures, such as pore morphology and size, that can better recapitulate diverse tissues. We utilized these interpenetrating ColHAX hydrogels as in vitro tissue models to examine macromolecular transport and recovery for early-stage drug screening. Hydrogel formulations with varying collagen and HAX concentrations imparted different gel properties based on the ratio of collagen to HAX. These gels were stable and swelled up to 170% of their original mass, and the storage moduli of the ColHAX gels increased over an order of magnitude by increasing collagen and HA concentration. Interestingly, when HAX concentration was constant and collagen concentration increased, both the pore size and spatial colocalization of collagen and HA increased. HA in the system dominated the ζ-potentials of the gels. The hydrogel and macromolecule properties impacted the mass transport and recovery of lysozyme, β-lactoglobulin, and bovine serum albumin (BSA) from the ColHAX gels─large molecules were largely impacted by mesh size, whereas small molecules were influenced primarily by electrostatic forces. Overall, the tunable properties demonstrated by the ColHAX hydrogels can be used to mimic different tissues for early-stage assays to understand drug transport and its relationship to matrix properties.

Extension of the SUGRES-1P Coarse-Grained Model of Polysaccharides to Heparin

Heparin is an unbranched periodic polysaccharide composed of negatively charged monomers and involved in key biological processes, including anticoagulation, angiogenesis, and inflammation. Its structure and dynamics have been studied extensively using experimental as well as theoretical approaches. The conventional approach of computational chemistry applied to the analysis of biomolecules is all-atom molecular dynamics, which captures the interactions of individual atoms by solving Newton's equation of motion. An alternative is molecular dynamics simulations using coarse-grained models of biomacromolecules, which offer a reduction of the representation and consequently enable us to extend the time and size scale of simulations by orders of magnitude. In this work, we extend the UNIfied COarse-gRaiNed (UNICORN) model of biological macromolecules developed in our laboratory to heparin. We carried out extensive tests to estimate the optimal weights of energy terms of the effective energy function as well as the optimal Debye-Hückel screening factor for electrostatic interactions. We applied the model to study unbound heparin molecules of polymerization degree ranging from 6 to 68 residues. We compare the obtained coarse-grained heparin conformations with models obtained from X-ray diffraction studies of heparin. The SUGRES-1P force field was able to accurately predict the general shape and global characteristics of heparin molecules.

A One-Bead-Per-Saccharide (1BPS) Model for Glycosaminoglycans

Glycosaminoglycans (GAGs) are polysaccharide compounds that play key roles in various biological processes. GAGs are important structural components of cartilage and the extracellular matrix of the brain. Due to the large size of these polysaccharides, coarse-grained approaches are indispensable for modeling these biopolymers. We develop a one-bead-per-saccharide model of chondroitin sulfates and hyaluronic acid based on an existing three-bead-per-saccharide coarse-grained model. Our coarse graining is carried out by using iterative Boltzmann inversion (IBI), including an additional coupling potential to incorporate the correlation between dihedral angles. The predictions of the model are verified against those of the existing three-bead-per-saccharin model and the experimental radius of gyration for hyaluronic acid.

Chondroitin sulfate, dermatan sulfate, and hyaluronic acid differentially modify the biophysical properties of collagen-based hydrogels

Fibrillar collagens and glycosaminoglycans (GAGs) are structural biomolecules that are natively abundant to the extracellular matrix (ECM). Prior studies have quantified the effects of GAGs on the bulk mechanical properties of the ECM. However, there remains a lack of experimental studies on how GAGs alter other biophysical properties of the ECM, including ones that operate at the length scales of individual cells such as mass transport efficiency and matrix microstructure. Here we characterized and decoupled the effects of the GAG molecules chondroitin sulfate (CS) dermatan sulfate (DS) and hyaluronic acid (HA) on the stiffness (indentation modulus), transport (hydraulic permeability), and matrix microarchitecture (pore size and fiber radius) properties of collagen-based hydrogels. We complement these biophysical measurements of collagen hydrogels with turbidity assays to profile collagen aggregate formation. Here we show that CS, DS, and HA differentially regulate the biophysical properties of hydrogels due to their alterations to the kinetics of collagen self-assembly. In addition to providing information on how GAGs play significant roles in defining key physical properties of the ECM, this work shows new ways in which stiffness measurements, microscopy, microfluidics, and turbidity kinetics can be used complementary to reveal details of collagen self-assembly and structure.

Single-molecule AFM study of hyaluronic acid softening in electrolyte solutions

Investigation of hyaluronic acid (HA) morphology and mechanical properties at a single-molecule level is important for the development of HA based biomaterials. We have developed the atomic force microscopy (AFM) based approach for quantitative characterization of conformation of HA molecules. HA molecules adsorbed on a modified graphitic surface form oriented linear segments. Conformation of HA molecules can be considered as two-dimensional quasi-projection of a three-dimensional conformation locally straightened by a substrate. The persistence length and Young's modulus of biomolecules estimated using wormlike chain model decrease from 15.7 to 9.9 nm, and from ∼21 to ∼13 GPa, respectively, when KCl concentration increases from 0 to 100 mM. The dependence of the persistence length on ionic strength supports the Odijk-Skolnick-Fixman model of polyelectrolyte stiffening in electrolyte solution. The obtained results represent a new insight into the conformation and mechanical characteristics of HA molecules and complement the characterization of this biopolymer by bulk methods.

Perineuronal nets restrict transport near the neuron surface: A coarse-grained molecular dynamics study

Perineuronal nets (PNNs) are mesh-like extracellular matrix structures that wrap around certain neurons in the central nervous system. They are hypothesized to stabilize memories in the brain and act as a barrier between cell and extracellular space. As a means to study the impact of PNNs on diffusion, the nets were approximated by negatively charged polymer brushes and simulated by coarse-grained molecular dynamics. Diffusion constants of single neutral and single charged particles were obtained in directions parallel and perpendicular to the brush substrate. The results for the neutral particle were compared to different theories of diffusion in a heuristic manner. Diffusion was found to be considerably reduced for brush spacings smaller than 10 nm, with a pronounced anisotropy for dense brushes. The exact dynamics of the chains was found to have a negligible impact on particle diffusion. The resistance of the brush proved small compared to typical values of the membrane resistance of a neuron, indicating that PNNs likely contribute little to the total resistance of an enwrapped neuron.

Atomic-Resolution Experimental Structural Biology and Molecular Dynamics Simulations of Hyaluronan and Its Complexes

This review summarizes the atomic-resolution structural biology of hyaluronan and its complexes available in the Protein Data Bank, as well as published studies of atomic-resolution explicit-solvent molecular dynamics simulations on these and other hyaluronan and hyaluronan-containing systems. Advances in accurate molecular mechanics force fields, simulation methods and software, and computer hardware have supported a recent flourish in such simulations, such that the simulation publications now outnumber the structural biology publications by an order of magnitude. In addition to supplementing the experimental structural biology with computed dynamic and thermodynamic information, the molecular dynamics studies provide a wealth of atomic-resolution information on hyaluronan-containing systems for which there is no atomic-resolution structural biology either available or possible. Examples of these summarized in this review include hyaluronan pairing with other hyaluronan molecules and glycosaminoglycans, with ions, with proteins and peptides, with lipids, and with drugs and drug-like molecules. Despite limitations imposed by present-day computing resources on system size and simulation timescale, atomic-resolution explicit-solvent molecular dynamics simulations have been able to contribute significant insight into hyaluronan's flexibility and capacity for intra- and intermolecular non-covalent interactions.

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.

Novel perspectives of sodium handling in type 2 diabetes mellitus

INTRODUCTION

As a key regulator of body water, sodium homeostasis forms an essential component of human physiology. Type 2 Diabetes Mellitus (T2D)-associated sodium overload stems from chronic renal retention of sodium, contributing toward the development of adverse cardiovascular sequelae.

AREAS COVERED

Our traditional model of sodium regulation invokes two compartments: extracellular fluid (ECF [plasma and interstitial fluid]) and intracellular fluid (ICF). Data from the Mars program reveal inconsistencies with this two-space model, including mismatches between net body sodium and water. Recent data utilizing 23Na magnetic resonance imaging (MRI) show a preponderance of bound sodium within human dermis, consistent with a third space repository and providing compelling evidence to support a three-space model in which dermal sodium binding facilitates sodium homeostasis within the ECF and ICF. This buffer is impaired in T2D, with diminishment of dermal bound sodium that may promote deleterious sequelae of sodium overload within the ECF and ICF.

EXPERT OPINION

Future studies should focus on novel therapeutic opportunities for sodium regulation in T2D and other conditions of sodium dysregulation. The ratio of free:bound dermal sodium (reflecting sodium storage capacity) could be utilized as a clinical biomarker for salt and water balance, to improve diagnostic accuracy and facilitate clinical decision-making.

Optimizing the sensitivity and resolution of hyaluronan analysis with solid-state nanopores

Hyaluronan (HA) is an essential carbohydrate in vertebrates that is a potentially robust bioindicator due to its critical roles in diverse physiological functions in health and disease. The intricate size-dependent function that exists for HA and its low abundance in most biological fluids have highlighted the need for sensitive technologies to provide accurate and quantitative assessments of polysaccharide molecular weight and concentration. We have demonstrated that solid state (SS-) nanopore technology can be exploited for this purpose, given its molecular sensitivity and analytical capacity, but there remains a need to further understand the impacts of experimental variables on the SS-nanopore signal for optimal interpretation of results. Here, we use model quasi-monodisperse HA polymers to determine the dependence of HA signal characteristics on a range of SS-nanopore measurement conditions, including applied voltage, pore diameter, and ionic buffer asymmetry. Our results identify important factors for improving the signal-to-noise ratio, resolution, and sensitivity of HA analysis with SS-nanopores.

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.

Poly-phosphocholination of liposomes leads to highly-extended retention time in mice joints

Surface-attached layers of phosphatidylcholine (PC) lipid vesicles (liposomes) may reduce the friction coefficient μ (= force-to-slide/load) between the sliding surfaces down to μ ≈ 10⁻³–10⁻⁴ up to tens of atm contact pressures, as high as those in the major joints (hips or knees). Such friction reduction is attributed to hydration lubrication by the highly-hydrated phosphocholine head-groups exposed at the outer vesicle surfaces. It has been suggested therefore that intra-articular (IA) administration of liposomes as potential boundary lubricants may alleviate degenerative, friction-associated joint conditions such as osteoarthritis (OA), which is associated with insufficient lubrication at the articular cartilage surface. To overcome the problem, common to all nanoparticles, of rapid removal by the mononuclear phagocyte system, as well as to ensure long-term colloidal stability during storage, functionalizing liposomes with poly(ethylene glycol) moieties, PEGylation, is often used. Here we describe a different liposome functionalization approach, using poly(2-methacryloyloxyethyl phosphorylcholine), PMPC, moieties (strictly, lipid–PMPC conjugates), and compare the retention time in mice joints of such PMPCylated liposomes with otherwise-identical but PEGylated vesicles following IA administration. We find, using fluorescence labeling and in vivo optical imaging, that when PMPC-stabilized liposomes are injected into mice knee joints, there is a massive increase of the vesicles’ retention half-life in the joints of about (4–5)-fold (ca. 300–400% increase in retention time) compared with the PEGylated liposomes (and some 100-fold longer than the retention time of intra-articularly injected hyaluronan or HA). Such PMPCylated liposomes are therefore promising candidates as potential long-lived boundary lubricants at the articular cartilage surface, with implication for friction-associated pathologies. Moreover, as lipid vesicles are well known to be efficient drug carriers, such long retention in the joints may enable analgesic or anti-inflammatory agents for joint pathologies to be more efficiently delivered via IA administration using PMPCylated liposomal vehicles relative to PEGylated ones.

PMPCylated liposomes injected into mice joints show a massive increase in retention half-life compared with PEGylated liposomes (or hyaluronan, HA), making them promising candidates as boundary lubricants at articular cartilage, or as drug carriers.

Theory and Practice of Coarse-Grained Molecular Dynamics of Biologically Important Systems

Molecular dynamics with coarse-grained models is nowadays extensively used to simulate biomolecular systems at large time and size scales, compared to those accessible to all-atom molecular dynamics. In this review article, we describe the physical basis of coarse-grained molecular dynamics, the coarse-grained force fields, the equations of motion and the respective numerical integration algorithms, and selected practical applications of coarse-grained molecular dynamics. We demonstrate that the motion of coarse-grained sites is governed by the potential of mean force and the friction and stochastic forces, resulting from integrating out the secondary degrees of freedom. Consequently, Langevin dynamics is a natural means of describing the motion of a system at the coarse-grained level and the potential of mean force is the physical basis of the coarse-grained force fields. Moreover, the choice of coarse-grained variables and the fact that coarse-grained sites often do not have spherical symmetry implies a non-diagonal inertia tensor. We describe selected coarse-grained models used in molecular dynamics simulations, including the most popular MARTINI model developed by Marrink's group and the UNICORN model of biological macromolecules developed in our laboratory. We conclude by discussing examples of the application of coarse-grained molecular dynamics to study biologically important processes.

Sulfation and Calcium Favor Compact Conformations of Chondroitin in Aqueous Solutions

The effects of sulfation and calcium cations (Ca²⁺) on the atomic-resolution conformational properties of chondroitin carbohydrate polymers in aqueous solutions are not well studied owing to experimental challenges. Here, we compare all-atom explicit-solvent molecular dynamics simulations results for pairs of O-type (nonsulfated) and A-type (GlcNAc 4-O-sulfated) chondroitin 20-mers in 140 mM NaCl with and without Ca²⁺ and find that both sulfation and Ca²⁺ favor more compact polymer conformations. We also show that subtle differences in force-field parametrization can have dramatic effects on Ca²⁺ binding to chondroitin carboxylate and sulfate functional groups and thereby determine Ca²⁺-mediated intra- and interstrand association. In addition to providing an atomic-resolution picture of the interaction of Ca²⁺ with sulfated and nonsulfated chondroitin polymers, the molecular dynamics data emphasize the importance of careful force-field parametrization to balance ion-water and ion-chondroitin interactions and suggest additional parametrization efforts to tune interactions involving sulfate.

Investigation of the structure of regulatory proteins interacting with glycosaminoglycans by combining NMR spectroscopy and molecular modeling - the beginning of a wonderful friendship

The interaction of regulatory proteins with extracellular matrix or cell surface-anchored glycosaminoglycans (GAGs) plays important roles in molecular recognition, wound healing, growth, inflammation and many other processes. In spite of their high biological relevance, protein-GAG complexes are significantly underrepresented in structural databases because standard tools for structure determination experience difficulties in studying these complexes. Co-crystallization with subsequent X-ray analysis is hampered by the high flexibility of GAGs. NMR spectroscopy experiences difficulties related to the periodic nature of the GAGs and the sparse proton network between protein and GAG with distances that typically exceed the detection limit of nuclear Overhauser enhancement spectroscopy. In contrast, computer modeling tools have advanced over the last years delivering specific protein-GAG docking approaches successfully complemented with molecular dynamics (MD)-based analysis. Especially the combination of NMR spectroscopy in solution providing sparse structural constraints with molecular docking and MD simulations represents a useful synergy of forces to describe the structure of protein-GAG complexes. Here we review recent methodological progress in this field and bring up examples where the combination of new NMR methods along with cutting-edge modeling has yielded detailed structural information on complexes of highly relevant cytokines with GAGs.

Molecular architecture and domain arrangement of the placental malaria protein VAR2CSA suggests a model for carbohydrate binding

VAR2CSA is the placental-malaria-specific member of the antigenically variant Plasmodium falciparum erythrocyte membrane protein 1 (PfEMP1) family. It is expressed on the surface of Plasmodium falciparum-infected host red blood cells and binds to specific chondroitin-4-sulfate chains of the placental proteoglycan receptor. The functional ∼310 kDa ectodomain of VAR2CSA is a multidomain protein that requires a minimum 12-mer chondroitin-4-sulfate molecule for specific, high affinity receptor binding. However, it is not known how the individual domains are organized and interact to create the receptor-binding surface, limiting efforts to exploit its potential as an effective vaccine or drug target. Using small angle X-ray scattering and single particle reconstruction from negative-stained electron micrographs of the ectodomain and multidomain constructs, we have determined the structural architecture of VAR2CSA. The relative locations of the domains creates two distinct pores that can each accommodate the 12-mer of chondroitin-4-sulfate, suggesting a model for receptor binding. This model has important implications for understanding cytoadherence of infected red blood cells and potentially provides a starting point for developing novel strategies to prevent and/or treat placental malaria.

A multiscale coarse-grained model to predict the molecular architecture and drug transport properties of modified chitosan hydrogels

Hydrogels constructed with functionalized polysaccharides are of interest in a multitude of applications, chiefly the design of therapeutic and regenerative formulations. Tailoring the chemical modification of polysaccharide-based hydrogels to achieve specific drug release properties involves the optimization of many tunable parameters, including (i) the type, degree (χ), and pattern of the functional groups, (ii) the water-polymer ratio, and (iii) the drug payload. To guide the design of modified polysaccharide hydrogels for drug release, we have developed a computational toolbox that predicts the structure and physicochemical properties of acylated chitosan chains, and their impact on the transport of drug molecules. Herein, we present a multiscale coarse-grained model to investigate the structure of networks of chitosan chains modified with acetyl, butanoyl, or heptanoyl moieties, as well as the diffusion of drugs doxorubicin (Dox) and gemcitabine (Gem) through the resulting networks. The model predicts the formation of different network structures, in particular the hydrophobically-driven transition from a uniform to a cluster/channel morphology and the formation of fibers of chitin chains. The model also describes the impact of structural and physicochemical properties on drug transport, which was confirmed experimentally by measuring Dox and Gem diffusion through an ensemble of modified chitosan hydrogels.

The need to freeze-Dehydration during specimen preparation for electron microscopy collapses the endothelial glycocalyx regardless of fixation method

OBJECTIVE

The endothelial glycocalyx covers the luminal surface of the endothelium and plays key roles in vascular function. Despite its biological importance, ideal visualization techniques are lacking. The current study aimed to improve the preservation and subsequent imaging quality of the endothelial glycocalyx.

METHODS

In mice, the endothelial glycocalyx was contrasted with a mixture of lanthanum and dysprosium (LaDy). Standard chemical fixation was compared with high-pressure frozen specimens processed with freeze substitution. Also, isolated brain microvessels and cultured endothelial cells were high-pressure frozen and by transmission soft x-rays, imaged under cryogenic conditions.

RESULTS

The endothelial glycocalyx was in some tissues significantly more voluminous from chemically fixed specimens compared with high-pressure frozen specimens. LaDy labeling introduced excessive absorption contrast, which impeded glycocalyx measurements in isolated brain microvessels when using transmission soft x-rays. In non-contrasted vessels, the glycocalyx was not resolved. LaDy-contrasted, cultured brain endothelial cells allowed to assess glycocalyx volume in vitro.

CONCLUSIONS

Both chemical and cryogenic fixation followed by dehydration lead to substantial collapse of the glycocalyx. Cryogenic fixation without freeze substitution could be a way forward although transmission soft x-ray tomography based solely on amplitude contrast seems unsuitable.

Genetic Control of Radical Cross-linking in a Semisynthetic Hydrogel

Enhancing materials with the qualities of living systems, including sensing, computation, and adaptation, is an important challenge in designing next-generation technologies. Living materials address this challenge by incorporating live cells as actuating components that control material function. For abiotic materials, this requires new methods that couple genetic and metabolic processes to material properties. Toward this goal, we demonstrate that extracellular electron transfer (EET) from Shewanella oneidensis can be leveraged to control radical cross-linking of a methacrylate-functionalized hyaluronic acid hydrogel. Cross-linking rates and hydrogel mechanics, specifically storage modulus, were dependent on various chemical and biological factors, including S. oneidensis genotype. Bacteria remained viable and metabolically active in the networks for a least 1 week, while cell tracking revealed that EET genes also encode control over hydrogel microstructure. Moreover, construction of an inducible gene circuit allowed transcriptional control of storage modulus and cross-linking rate via the tailored expression of a key electron transfer protein, MtrC. Finally, we quantitatively modeled hydrogel stiffness as a function of steady-state mtrC expression and generalized this result by demonstrating the strong relationship between relative gene expression and material properties. This general mechanism for radical cross-linking provides a foundation for programming the form and function of synthetic materials through genetic control over extracellular electron transfer.

Chondroitin-Sulfate-A-Coated Magnetite Nanoparticles: Synthesis, Characterization and Testing to Predict Their Colloidal Behavior in Biological Milieu

Biopolymer coated magnetite nanoparticles (MNPs) are suitable to fabricate biocompatible magnetic fluid (MF). Their comprehensive characterization, however, is a necessary step to assess whether bioapplications are feasible before expensive in vitro and in vivo tests. The MNPs were prepared by co-precipitation, and after careful purification, they were coated by chondroitin-sulfate-A (CSA). CSA exhibits high affinity adsorption to MNPs (H-type isotherm). We could only make stable MF of CSA coated MNPs (CSA@MNPs) under accurate conditions. The CSA@MNP was characterized by TEM (size ~10 nm) and VSM (saturation magnetization ~57 emu/g). Inner-sphere metal–carboxylate complex formation between CSA and MNP was proved by FTIR-ATR and XPS. Electrophoresis and DLS measurements show that the CSA@MNPs at CSA-loading > 0.2 mmol/g were stable at pH > 4. The salt tolerance of the product improved up to ~0.5 M NaCl at pH~6.3. Under favorable redox conditions, no iron leaching from the magnetic core was detected by ICP measurements. Thus, the characterization predicts both chemical and colloidal stability of CSA@MNPs in biological milieu regarding its pH and salt concentration. MTT assays showed no significant impact of CSA@MNP on the proliferation of A431 cells. According to these facts, the CSA@MNPs have a great potential in biocompatible MF preparation for medical applications.

Perspective on computational simulations of glycosaminoglycans

Glycosaminoglycans (GAGs) represent a formidable frontier for chemists, biochemists, biologists, medicinal chemists and drug delivery specialists because of massive structural complexity. GAGs are arguably the most complex, natural linear biopolymers with theoretical diversity orders of magnitude higher than proteins and nucleic acids. Yet, this diversity remains generally untapped. Computational approaches offer major routes to understand GAG structure and dynamics so as to enable novel applications of these biopolymers. In fact, computational algorithms, softwares, online tools and techniques have reached a level of sophistication that help understand atomistic details of conformational variation and protein recognition of individual GAG sequences. This review describes current approaches and challenges in computational study of GAGs. It presents a history of major findings since the earliest mention of GAGs (the 1960s), the development of parameters and force fields specific for GAGs, and the application of these tools in understanding GAG structure-function relationship. This review also presents a section on how to perform simulation of GAGs, which is directed toward researchers interested in entering this promising field with potential to impact therapy.

Flowable Polyethylene Glycol Hydrogels Support the in Vitro Survival and Proliferation of Dermal Progenitor Cells in a Mechanically Dependent Manner

Cell-based therapies have garnered considerable interest largely because of their potential utility for tissue regeneration in a variety of organs, including skin. Designing vehicles that enable optimal delivery and purposeful integration of donor cells within tissues will be critical for their success. Here, we investigate the utility of an injectable, self-polymerizing, fully synthetic hydrogel in supporting the survival, proliferation, and function of cultured adult dermal progenitor cells (DPCs) which may serve as a source of renewable cells to repair severe skin injuries or restore hair growth. We show that modifying the stiffness of these transglutaminase cross-linked poly(ethylene glycol) (TG-PEG) hydrogels significantly alters DPC behavior and phenotype; increasing stiffness promotes their differentiation and migration whereas softer gels maintained them in a proliferative state. We found that 2-3% TG-PEG was optimal to promote cell expansion and survival. Unexpectedly, DPCs grown in all conditions maintained their inductive function and thus generated de novo hair follicles. Our data suggests that TG-PEG hydrogels may be a versatile platform for stem and progenitor cell transplantation and fate specification while maintaining functional competence.

Molecular weight impact on the mechanical forces between hyaluronan and its receptor

Hyaluronan (HA) possesses manifold mechanical and signaling properties in the body. Most of these activities are largely regulated by its molecular weight, which often triggers opposing functions. However the molecular basis for such function distinction between HA size categories remains unclear. Using a combination of biophysical techniques, we measured the physical forces between HA ligand and its specific receptor CD44 in both normal and lateral directions, at different HA molecular weights and bound states. It was found that the impact of HA multivalency is more than just the sum of separate monovalent bindings. The HA-CD44 specific interaction enhances with HA molecular weight and the maximum binding occurs at ∼1000 kD, possibly due to the balance between multivalent HA zipping effect and conformational entropy. High friction patches, probably from CD44 protein clustering, was observed in friction force microscopy (FFM) upon HA shearing, which is also dependent on HA molecular weight. These results could help to understand the biophysical mechanism of HA in regulating CD44-induced physiological activities and thus facilitate the new design of HA-based material in fine tuning the receptor responses.

Heparosan-coated liposomes for drug delivery

Liposomal encapsulation is a useful drug delivery strategy for small molecules, especially chemotherapeutic agents such as doxorubicin. Doxil® is a doxorubicin-containing liposome ("dox-liposome") that passively targets drug to tumors while reducing side effects caused by free drug permeating and poisoning healthy tissues. Polyethylene glycol (PEG) is the hydrophilic coating of Doxil® that protects the formulation from triggering the mononuclear phagocyte system (MPS). Evading the MPS prolongs dox-liposome circulation time thus increasing drug deposition at the tumor site. However, multiple doses of Doxil® sometimes activate an anti-PEG immune response that enhances liposome clearance from circulation and causes hypersensitivity, further limiting its effectiveness against disease. These side effects constrain the utility of PEG-coated liposomes in certain populations, justifying the need for investigation into alternative coatings that could improve drug delivery for better patient quality of life and outcome. We hypothesized that heparosan (HEP; [-4-GlcA-β1-4-GlcNAc-α1-]n) may serve as a PEG alternative for coating liposomes. HEP is a natural precursor to heparin biosynthesis in mammals. Also, bacteria expressing an HEP extracellular capsule during infection escape detection and are recognized as "self," not a foreign threat. By analogy, coating drug-carrying liposomes with HEP should camouflage the delivery vehicle from the MPS, extending circulation time and potentially avoiding immune-mediated clearance. In this study, we characterize the postmodification insertion of HEP-lipids into liposomes by dynamic light scattering and coarse-grain computer modeling, test HEP-lipid immunogenicity in rats, and compare the efficacy of drug delivered by HEP-coated liposomes to PEG-coated liposomes in a human breast cancer xenograft mouse model.

Glycosaminoglycans in extracellular matrix organisation: are concepts from soft matter physics key to understanding the formation of perineuronal nets?

Conventional wisdom has it that proteins fold and assemble into definite structures, and that this defines their function. Glycosaminoglycans (GAGs) are different. In most cases the structures they form have a low degree of order, even when interacting with proteins. Here, we discuss how physical features common to all GAGs-hydrophilicity, charge, linearity and semi-flexibility-underpin the overall properties of GAG-rich matrices. By integrating soft matter physics concepts (e.g. polymer brushes and phase separation) with our molecular understanding of GAG-protein interactions, we can better comprehend how GAG-rich matrices assemble, what their properties are, and how they function. Taking perineuronal nets (PNNs)-a GAG-rich matrix enveloping neurons-as a relevant example, we propose that microphase separation determines the holey PNN anatomy that is pivotal to PNN functions.

Multiscale modeling of glycosaminoglycan structure and dynamics: current methods and challenges

Glycosaminoglycans are long unbranched and complex polysaccharides that are an essential component of mammalian extracellular matrices. Characterization of their molecular structure, dynamics and interactions are essential to understand important biological phenomena in health and disease, and will lead to novel therapeutics and medical devices. However, this has proven to be a challenge experimentally and theoretical techniques are needed to develop new hypotheses, and interpret experiments. This review aims to examine the current theoretical (rather than experimental) methods used by researchers to investigate glycosaminoglycan structure, dynamics and interactions, from the monosaccharide to the macromolecular scale. It will consider techniques such as quantum mechanics, molecular mechanics, molecular dynamics, coarse graining and docking.

An accurate coarse-grained model for chitosan polysaccharides in aqueous solution

Computational models can provide detailed information about molecular conformations and interactions in solution, which is currently inaccessible by other means in many cases. Here we describe an efficient and precise coarse-grained model for long polysaccharides in aqueous solution at different physico-chemical conditions such as pH and ionic strength. The Model is carefully constructed based on all-atom simulations of small saccharides and metadynamics sampling of the dihedral angles in the glycosidic links, which represent the most flexible degrees of freedom of the polysaccharides. The model is validated against experimental data for Chitosan molecules in solution with various degree of deacetylation, and is shown to closely reproduce the available experimental data. For long polymers, subtle differences of the free energy maps of the glycosidic links are found to significantly affect the measurable polymer properties. Therefore, for titratable monomers the free energy maps of the corresponding links are updated according to the current charge of the monomers. We then characterize the microscopic and mesoscopic structural properties of large chitosan polysaccharides in solution for a wide range of solvent pH and ionic strength, and investigate the effect of polymer length and degree and pattern of deacetylation on the polymer properties.

Brain tissue stiffness is a sensitive marker for acidosis

•

We introduce tissue stiffness as a sensitive indicator for pathophysiological changes to CNS tissue.

•

We applied atomic force microscopy to investigate tissue stiffness.

•

We found that CO₂ overexposure-induced acidosis changes brain properties.

Background

Carbon dioxide overdose is frequently used to cull rodents for tissue harvesting. However, this treatment may lead to respiratory acidosis, which potentially could change the properties of the investigated tissue.

New method

Mechanical tissue properties often change in pathological conditions and may thus offer a sensitive generic readout for changes in biological tissues with clinical relevance. In this study, we performed force-indentation measurements with an atomic force microscope on acute cerebellar slices from adult rats to test if brain tissue undergoes changes following overexposure to CO₂ compared to other methods of euthanasia.

Results

The pH significantly decreased in brain tissue of animals exposed to CO₂. Concomitant with the drop in pH, cerebellar grey matter significantly stiffened. Tissue stiffening was reproduced by incubation of acute cerebellar slices in acidic medium.

Comparison with existing methods

Tissue stiffness provides an early, generic indicator for pathophysiological changes in the CNS. Atomic force microscopy offers unprecedented high spatial resolution to detect such changes.

Conclusions

Our results indicate that the stiffness particularly of grey matter strongly correlates with changes of the pH in the cerebellum. Furthermore, the method of tissue harvesting and preparation may not only change tissue stiffness but very likely also other physiologically relevant parameters, highlighting the importance of appropriate sample preparation.

Computational analysis of interactions in structurally available protein-glycosaminoglycan complexes

Glycosaminoglycans represent a class of linear anionic periodic polysaccharides, which play a key role in a variety of biological processes in the extracellular matrix via interactions with their protein targets. Computationally, glycosaminoglycans are very challenging due to their high flexibility, periodicity and electrostatics-driven nature of the interactions with their protein counterparts. In this work, we carry out a detailed computational characterization of the interactions in protein-glycosaminoglycan complexes from the Protein Data Bank (PDB), which are split into two subsets accounting for their intrinsic nature: non-enzymatic-protein-glycosaminoglycan and enzyme-glycosaminoglycan complexes. We apply molecular dynamics to analyze the differences in these two subsets in terms of flexibility, retainment of the native interactions in the simulations, free energy components of binding and contributions of protein residue types to glycosaminoglycan binding. Furthermore, we systematically demonstrate that protein electrostatic potential calculations, previously found to be successful for glycosaminoglycan binding sites prediction for individual systems, are in general very useful for proposing protein surface regions as putative glycosaminoglycan binding sites, which can be further used for local docking calculations with these particular polysaccharides. Finally, the performance of six different docking programs (Autodock 3, Autodock Vina, MOE, eHiTS, FlexX and Glide), some of which proved to perform well for particular protein-glycosaminoglycan complexes in previous work, is evaluated on the complete protein-glycosaminoglycan data set from the PDB. This work contributes to widen our knowledge of protein-glycosaminoglycan molecular recognition and could be useful to steer a choice of the strategies to be applied in theoretical studies of these systems.

Chondrocytes, Mesenchymal Stem Cells, and Their Combination in Articular Cartilage Regenerative Medicine

Articular cartilage (AC) is a highly organized connective tissue lining, covering the ends of bones within articulating joints. Its highly ordered structure is essential for stable motion and provides a frictionless surface easing load transfer. AC is vulnerable to lesions and, because it is aneural and avascular, it has limited self-repair potential which often leads to osteoarthritis. To date, no fully successful treatment for osteoarthritis has been reported. Thus, the development of innovative therapeutic approaches is desperately needed. Autologous chondrocyte implantation, the only cell-based surgical intervention approved in the United States for treating cartilage defects, has limitations because of de-differentiation of articular chondrocytes (AChs) upon in vitro expansion. De-differentiation can be abated if initial populations of AChs are co-cultured with mesenchymal stem cells (MSCs), which not only undergo chondrogenesis themselves but also support chondrocyte vitality. In this review we summarize studies utilizing AChs, non-AChs, and MSCs and compare associated outcomes. Moreover, a comprehensive set of recent human studies using chondrocytes to direct MSC differentiation, MSCs to support chondrocyte re-differentiation and proliferation in co-culture environments, and exploratory animal intra- and inter-species studies are systematically reviewed and discussed in an innovative manner allowing side-by-side comparisons of protocols and outcomes. Finally, a comprehensive set of recommendations are made for future studies.

Tension Amplification in Tethered Layers of Bottle-Brush Polymers

Molecular dynamics simulations of a coarse-grained bead-spring model have been used to study the effects of molecular crowding on the accumulation of tension in the backbone of bottle-brush polymers tethered to a flat substrate. The number of bottle-brushes per unit surface area, Σ, as well as the lengths of the bottle-brush backbones Nbb (50 ≤ Nbb ≤ 200) and side chains Nsc (50 ≤ Nsc ≤ 200) were varied to determine how the dimensions and degree of crowding of bottle-brushes give rise to bond tension amplification along the backbone, especially near the substrate. From these simulations, we have identified three separate regimes of tension. For low Σ, the tension is due solely to intramolecular interactions and is dominated by the side chain repulsion that governs the lateral brush dimensions. With increasing Σ, the interactions between bottle-brush polymers induce compression of the side chains, transmitting increasing tension to the backbone. For large Σ, intermolecular side chain repulsion increases, forcing side chain extension and reorientation in the direction normal to the surface and transmitting considerable tension to the backbone.

Enzyme responsive GAG-based natural-synthetic hybrid hydrogel for tunable growth factor delivery and stem cell differentiation

We describe an enzymatically formed chondroitin sulfate (CS) and poly(ethylene glycol) (PEG) based hybrid hydrogel system, which by tuning the architecture and composition of modular building blocks, allows the application-specific tailoring of growth factor delivery and cellular responses. CS, a negatively charged sulfate-rich glycosaminoglycan of the extracellular matrix (ECM), known for its growth factor binding and stem cell regulatory functions, is used as a starting material for the engineering of this biomimetic materials platform. The functionalization of CS with transglutaminase factor XIII specific substrate sequences is utilized to allow cross-linking of CS with previously described fibrin-mimetic TG-PEG hydrogel precursors. We show that the hydrogel network properties can be tuned by varying the degree of functionalization of CS as well as the ratio and concentrations of PEG and CS precursors. Taking advantage of TG-PEG hydrogel, compatible tagged bio-functional building blocks, including RGD peptides or matrix metalloproteinase sensitive domains, can be incorporated on demand allowing the three-dimensional culture and expansion of human bone marrow mesenchymal stem cells (BM-MSCs). The binding of bone morphogenetic protein-2 (BMP-2) in a CS concentration dependent manner and the BMP-2 release mediated osteogenic differentiation of BM-MSCs indicate the potential of CS-PEG hybrid hydrogels to promote regeneration of bone tissue. Their modular design allows facile incorporation of additional signaling elements, rendering CS-PEG hydrogels a highly flexible platform with potential for multiple biomedical applications.

Nuclear Magnetic Resonance Insight into the Multiple Glycosaminoglycan Binding Modes of the Link Module from Human TSG-6

Tumor necrosis factor-stimulated gene-6 (TSG-6) is a hyaluronan (HA)-binding protein that is essential for stabilizing and remodeling the extracellular matrix (ECM) during ovulation and inflammatory disease processes such as arthritis. The Link module, one of the domains of TSG-6, is responsible for binding hyaluronan and other glycosaminoglycans found in the ECM. In this study, we used a well-defined chondroitin sulfate (CS) hexasaccharide (ΔC444S) to determine the structure of the Link module, in solution, in its chondroitin sulfate-bound state. A variety of nuclear magnetic resonance techniques were employed, including chemical shift perturbation, residual dipolar couplings (RDCs), nuclear Overhauser effects, spin relaxation measurements, and paramagnetic relaxation enhancements from a spin-labeled analogue of ΔC444S. The binding site for ΔC444S on the Link module overlapped with that of HA. Surprisingly, ΔC444S binding induced dimerization of the Link module (as confirmed by analytical ultracentrifugation), and a second weak binding site that partially overlapped with a previously identified heparin site was detected. A dimer model was generated using chemical shift perturbations and RDCs as restraints in the docking program HADDOCK. We postulate that the molecular cross-linking enhanced by the multiple binding modes of the Link module might be critical for remodeling the ECM during inflammation/ovulation and might contribute to other functions of TSG-6.

Polyplex Micelles with Double-Protective Compartments of Hydrophilic Shell and Thermoswitchable Palisade of Poly(oxazoline)-Based Block Copolymers for Promoted Gene Transfection

Improving the stability of polyplex micelles under physiological conditions is a critical issue for promoting gene transfection efficiencies. To this end, hydrophobic palisade was installed between the inner core of packaged plasmid DNA (pDNA) and the hydrophilic shell of polyplex micelles using a triblock copolymer consisting of hydrophilic poly(2-ethyl-2-oxazoline), thermoswitchable amphiphilic poly(2-n-propyl-2-oxazoline) (PnPrOx) and cationic poly(L-lysine). The two-step preparation procedure, mixing the triblock copolymer with pDNA below the lower critical solution temperature (LCST) of PnPrOx, followed by incubation above the LCST to form a hydrophobic palisade of the collapsed PnPrOx segment, induced the formation of spatially aligned hydrophilic-hydrophobic double-protected polyplex micelles. The prepared polyplex micelles exhibited significant tolerance against attacks from nuclease and polyanions compared to those without hydrophobic palisades, thereby promoting gene transfection. These results corroborated the utility of amphiphilic poly(oxazoline) as a molecular thermal switch to improve the stability of polyplex gene carriers relevant for physiological applications.

Aggrecan nanoscale solid-fluid interactions are a primary determinant of cartilage dynamic mechanical properties

Poroelastic interactions between interstitial fluid and the extracellular matrix of connective tissues are critical to biological and pathophysiological functions involving solute transport, energy dissipation, self-stiffening and lubrication. However, the molecular origins of poroelasticity at the nanoscale are largely unknown. Here, the broad-spectrum dynamic nanomechanical behavior of cartilage aggrecan monolayer is revealed for the first time, including the equilibrium and instantaneous moduli and the peak in the phase angle of the complex modulus. By performing a length scale study and comparing the experimental results to theoretical predictions, we confirm that the mechanism underlying the observed dynamic nanomechanics is due to solid-fluid interactions (poroelasticity) at the molecular scale. Utilizing finite element modeling, the molecular-scale hydraulic permeability of the aggrecan assembly was quantified (kaggrecan = (4.8 ± 2.8) × 10(-15) m(4)/N·s) and found to be similar to the nanoscale hydraulic permeability of intact normal cartilage tissue but much lower than that of early diseased tissue. The mechanisms underlying aggrecan poroelasticity were further investigated by altering electrostatic interactions between the molecule's constituent glycosaminoglycan chains: electrostatic interactions dominated steric interactions in governing molecular behavior. While the hydraulic permeability of aggrecan layers does not change across species and age, aggrecan from adult human cartilage is stiffer than the aggrecan from newborn human tissue.

Paracrine signaling mediated at cell-cell contacts

Recent findings in several organ systems show that cytoneme-mediated signaling transports signaling proteins along cellular extensions and targets cell-to-cell exchanges to synaptic contacts. This mechanism of paracrine signaling may be a general one that is used by many (or all) cell types in many (or all) organs. We briefly review these findings in this perspective. We also describe the properties of several signaling systems that have previously been interpreted to support a passive diffusion mechanism of signaling protein dispersion, but can now be understood in the context of the cytoneme mechanism. Also watch the Video Abstract.

Formulation of in situ chemically cross-linked hydrogel depots for protein release: from the blob model perspective

The fast release rate and the undesirable covalent binding are two major problems often encountered in formulating in situ chemically cross-linked hydrogel as protein release depot, particularly when prolonged release over months is desirable. In this study, we applied the De Gennes' blob theory to analyze and tackle these two problems using a vinylsulfone-thiol (VS-SH) reaction based in situ hydrogel system. We showed that the simple scaling relation ξb ≈ Rg(c/c*)(-v/(3v-1)) is applicable to the in situ hydrogel and the mesh size estimated from the precursor polymer parameters is a reasonable match to experimental results. On the other hand, as predicted by the theory and confirmed by experiments, the drug diffusion within hydrogel depends mainly on polymer concentration but not the degree of modification (DM). The covalent binding was found to be caused by the mismatch of location between the reactive groups and the entanglement points. The mismatch and, thus, the protein binding were minimized by increasing the DM and concentration of the SH polymer relative to the VS polymer, as predicted by theory. Using these principles, an in situ hydrogel system for the controlled release of an antiangiogenic antibody therapeutics bevacizumab for 3 months was developed.