Osmotic pressure of aqueous chondroitin sulfate solution: a molecular modeling investigation

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.

Viscoelastic Chondroitin Sulfate and Hyaluronic Acid Double-Network Hydrogels with Reversible Cross-Links

Viscoelastic hydrogels are gaining interest as they possess necessary requirements for bioprinting and injectability. By means of reversible, dynamic covalent bonds, it is possible to achieve features that recapitulate the dynamic character of the extracellular matrix. Dually cross-linked and double-network (DN) hydrogels seem to be ideal for the design of novel biomaterials and bioinks, as a wide range of properties required for mimicking advanced and complex tissues can be achieved. In this study, we investigated the fabrication of chondroitin sulfate/hyaluronic acid (CS/HA)-based DN hydrogels, in which two networks are interpenetrated and cross-linked with the dynamic covalent bonds of very different lifetimes. Namely, Diels-Alder adducts (between methylfuran and maleimide) and hydrazone bonds (between aldehyde and hydrazide) were chosen as cross-links, leading to viscoelastic hydrogels. Furthermore, we show that viscoelasticity and the dynamic character of the resulting hydrogels could be tuned by changing the composition, that is, the ratio between the two types of cross-links. Also, due to a very dynamic nature and short lifetime of hydrazone cross-links (∼800 s), the DN hydrogel is easily processable (e.g., injectable) in the first stages of gelation, allowing the material to be used in extrusion-based 3D printing. The more long-lasting and robust Diels-Alder cross-links are responsible for giving the network enhanced mechanical strength and structural stability. Being highly charged and hydrophilic, the cross-linked CS and HA enable a high swelling capacity (maximum swelling ratio ranging from 6 to 12), which upon confinement results in osmotically stiffened constructs, able to mimic the mechanical properties of cartilage tissue, with the equilibrium moduli ranging from 0.3 to 0.5 MPa. Moreover, the mesenchymal stromal cells were viable in the presence of the hydrogels, and the effect of the degradation products on the macrophages suggests their safe use for further translational applications. The DN hydrogels with dynamic covalent cross-links hold great potential for the development of novel smart and tunable viscoelastic materials to be used as biomaterial inks or bioinks in bioprinting and regenerative medicine.

Ionic effects in semi-dilute biopolymer solutions: A small angle scattering study

Systematic investigations using neutron and X-ray small angle scattering in near-physiological salt solutions were made to reveal the effect of polymer concentration, pH, and calcium ion concentration on the structure of semi-dilute solutions of four model biopolymers [polyaspartic acid, DNA, chondroitin sulfate, and hyaluronic acid (HA)] representing typical backbone structures. In the low q range (<0.01 Å-1), the scattering response I(q) is dominated by scattering from large clusters. In the intermediate q range, I(q) varies approximately as q -1, exposing the linear nature of the scatterers. In these polyelectrolyte solutions, the correlation length L displays a power law dependence on the polymer concentration c that resembles that of neutral polymer solutions. L increases with increasing calcium chloride concentration and with decreasing pH. The effect of the different divalent cations, Ba, Mg, Ca, Sr, and Mn, on the structure of DNA solutions is practically identical. However, in mixed salt conditions at the same ionic strength, the combined effect of mono- and divalent counter-ions on the structure of the polymer solutions deviates significantly from additivity. Anomalous small angle X-ray scattering observations on both DNA and HA solutions reveal that the divalent strontium counter-ions form a tight sheath around the polymer chain. The shape of the divalent ion cloud is similar in these two systems.

Osmolality predictive models of different polymers as tools in parenteral and ophthalmic formulation development

During the development of parenteral dosage forms, different physicochemical studies are required to ensure stable, effective and safe formulations. The osmolality of this kind of dosage forms should bear a close similarity to the body fluids to prevent local irritation, pain or even more significant side effects like endothelial damage. The osmotic studies performed in Polyethylene glycol 400 (PEG 400), Polyethylene glycol 4000 (PEG 4000), Poloxamer 407 (P407), Sodium Hyaluronate (SH), Chondroitin Sulphate Sodium (CS), Cremophor RH 40 (CRE40) and Polyvinyl alcohol (PVA) aqueous solutions, showed that the theoretical determination of the osmolality based on their molecular weight as the only determinant factor did not agree with the values obtained by the measurement of colligative properties such as the freezing point depression. The data obtained from this study and its analysis, provided predictive equations that can be used as tools in the primary development to estimate formulation's osmolality at different concentrations; and its evolution over a period at the hypothetical worst-case scenario of storage temperature.

Effect of crosslinking in cartilage-like collagen microstructures

The mechanical performance of biological tissues is underpinned by a complex and finely balanced structure. Central to this is collagen, the most abundant protein in our bodies, which plays a dominant role in the functioning of tissues, and also in disease. Based on the collagen meshwork of articular cartilage, we have developed a bottom-up spring-node model of collagen and examined the effect of fibril connectivity, implemented by crosslinking, on mechanical behaviour. Although changing individual crosslink stiffness within an order of magnitude had no significant effect on modelling predictions, the density of crosslinks in a meshwork had a substantial impact on its behaviour. Highly crosslinked meshworks maintained a 'normal' configuration under loading, with stronger resistance to deformation and improved recovery relative to sparsely crosslinked meshwork. Stress on individual fibrils, however, was higher in highly crosslinked meshworks. Meshworks with low numbers of crosslinks reconfigured to disease-like states upon deformation and recovery. The importance of collagen interconnectivity may provide insight into the role of ultrastructure and its mechanics in the initiation, and early stages, of diseases such as osteoarthritis.

Predicting Knee Osteoarthritis

Treatment options for osteoarthritis (OA) beyond pain relief or total knee replacement are very limited. Because of this, attention has shifted to identifying which factors increase the risk of OA in vulnerable populations in order to be able to give recommendations to delay disease onset or to slow disease progression. The gold standard is then to use principles of risk management, first to provide subject-specific estimates of risk and then to find ways of reducing that risk. Population studies of OA risk based on statistical associations do not provide such individually tailored information. Here we argue that mechanistic models of cartilage tissue maintenance and damage coupled to statistical models incorporating model uncertainty, united within the framework of structural reliability analysis, provide an avenue for bridging the disciplines of epidemiology, cell biology, genetics and biomechanics. Such models promise subject-specific OA risk assessment and personalized strategies for mitigating or even avoiding OA. We illustrate the proposed approach with a simple model of cartilage extracellular matrix synthesis and loss regulated by daily physical activity.

Altered mechano-chemical environment in hip articular cartilage: effect of obesity

The production of extracellular matrix (ECM) components of articular cartilage is regulated, among other factors, by an intercellular signaling mechanism mediated by the interaction of cell surface receptors (CSR) with insulin-like growth factor-1 (IGF-1). In ECM, the presence of binding proteins (IGFBP) hinders IGF-1 delivery to CSR. It has been reported that levels of IGF-1 and IGFBP in obese population are, respectively, lower and higher than those found in normal population. In this study, an experimental-numerical approach was adopted to quantify the effect of this metabolic alteration found in obese population on the homeostasis of femoral hip cartilage. A new computational model, based on the mechano-electrochemical mixture theory, was developed to describe competitive binding kinetics of IGF-1 with IGFBP and CSR, and associated glycosaminoglycan (GAG) biosynthesis. Moreover, a gait analysis was carried out on obese and normal subjects to experimentally characterize mechanical loads on hip cartilage during walking. This information was deployed into the model to account for effects of physiologically relevant tissue deformation on GAG production in ECM. Numerical simulations were performed to compare GAG biosynthesis in femoral hip cartilage of normal and obese subjects. Results indicated that the lower ratio of IGF-1 to IGFBP found in obese population reduces cartilage GAG concentration up to 18 % when compared to normal population. Moreover, moderate physical activity, such as walking, has a modest beneficial effect on GAG production. The findings of this study suggest that IGF-1/IGFBP metabolic unbalance should be accounted for when considering the association of obesity with hip osteoarthritis.

Computational model for the analysis of cartilage and cartilage tissue constructs

We propose a new non-linear poroelastic model that is suited to the analysis of soft tissues. In this paper the model is tailored to the analysis of cartilage and the engineering design of cartilage constructs. The proposed continuum formulation of the governing equations enables the strain of the individual material components within the extracellular matrix (ECM) to be followed over time, as the individual material components are synthesized, assembled and incorporated within the ECM or lost through passive transport or degradation. The material component analysis developed here naturally captures the effect of time-dependent changes of ECM composition on the deformation and internal stress states of the ECM. For example, it is shown that increased synthesis of aggrecan by chondrocytes embedded within a decellularized cartilage matrix initially devoid of aggrecan results in osmotic expansion of the newly synthesized proteoglycan matrix and tension within the structural collagen network. Specifically, we predict that the collagen network experiences a tensile strain, with a maximum of ~2% at the fixed base of the cartilage. The analysis of an example problem demonstrates the temporal and spatial evolution of the stresses and strains in each component of a self-equilibrating composite tissue construct, and the role played by the flux of water through the tissue.

Contribution of proteoglycan osmotic swelling pressure to the compressive properties of articular cartilage

The negatively charged proteoglycans (PG) provide compressive resistance to articular cartilage by means of their fixed charge density (FCD) and high osmotic pressure (π(PG)), and the collagen network (CN) provides the restraining forces to counterbalance π(PG). Our objectives in this work were to: 1), account for collagen intrafibrillar water when transforming biochemical measurements into a FCD-π(PG) relationship; 2), compute π(PG) and CN contributions to the compressive behavior of full-thickness cartilage during bovine growth (fetal, calf, and adult) and human adult aging (young and old); and 3), predict the effect of depth from the articular surface on π(PG) in human aging. Extrafibrillar FCD (FCD(EF)) and π(PG) increased with bovine growth due to an increase in CN concentration, whereas PG concentration was steady. This maturation-related increase was amplified by compression. With normal human aging, FCD(EF) and π(PG) decreased. The π(PG)-values were close to equilibrium stress (σ(EQ)) in all bovine and young human cartilage, but were only approximately half of σ(EQ) in old human cartilage. Depth-related variations in the strain, FCD(EF), π(PG), and CN stress profiles in human cartilage suggested a functional deterioration of the superficial layer with aging. These results suggest the utility of the FCD-π(PG) relationship for elucidating the contribution of matrix macromolecules to the biomechanical properties of cartilage.

Nanomechanics of the Cartilage Extracellular Matrix

Cartilage is a hydrated biomacromolecular fiber composite located at the ends of long bones that enables proper joint lubrication, articulation, loading, and energy dissipation. Degradation of extracellular matrix molecular components and changes in their nanoscale structure greatly influence the macroscale behavior of the tissue and result in dysfunction with age, injury, and diseases such as osteoarthritis. Here, the application of the field of nanomechanics to cartilage is reviewed. Nanomechanics involves the measurement and prediction of nanoscale forces and displacements, intra- and intermolecular interactions, spatially varying mechanical properties, and other mechanical phenomena existing at small length scales. Experimental nanomechanics and theoretical nanomechanics have been applied to cartilage at varying levels of material complexity, e.g., nanoscale properties of intact tissue, the matrix associated with single cells, biomimetic molecular assemblies, and individual extracellular matrix biomolecules (such as aggrecan, collagen, and hyaluronan). These studies have contributed to establishing a fundamental mechanism-based understanding of native and engineered cartilage tissue function, quality, and pathology.

Calcium diffusion enhanced after cleavage of negatively charged components of brain extracellular matrix by chondroitinase ABC

The concentration of extracellular calcium plays a critical role in synaptic transmission and neuronal excitability as well as other physiological processes. The time course and extent of local fluctuations in the concentration of this ion largely depend on its effective diffusion coefficient (D*) and it has been speculated that fixed negative charges on chondroitin sulphate proteoglycans (CSPGs) and other components of the extracellular matrix may influence calcium diffusion because it is a divalent cation. In this study we used ion-selective microelectrodes combined with pressure ejection or iontophoresis of ions from a micropipette to quantify diffusion characteristics of neocortex and hippocampus in rat brain slices. We show that D* for calcium is less than the value predicted from the behaviour of the monovalent cation tetramethylammonium (TMA), a commonly used diffusion probe, but D* for calcium increases in both brain regions after the slices are treated with chondroitinase ABC, an enzyme that predominantly cleaves chondroitin sulphate glycans. These results suggest that CSPGs do play a role in determining the local diffusion properties of calcium in brain tissue, most likely through electrostatic interactions mediating rapid equilibrium binding. In contrast, chondroitinase ABC does not affect either the TMA diffusion or the extracellular volume fraction, indicating that the enzyme does not alter the structure of the extracellular space and that the diffusion of small monovalent cations is not affected by CSPGs in the normal brain ionic milieu. Both calcium and CSPGs are known to have many distinct roles in brain physiology, including brain repair, and our study suggests they may be functionally coupled through calcium diffusion properties.

Structure and interactions of aggrecans: statistical thermodynamic approach

Weak polyelectrolytes tethered to cylindrical surfaces are investigated using a molecular theory. These polymers form a model system to describe the properties of aggrecan molecules, which is one of the main components of cartilage. We have studied the structural and thermodynamical properties of two interacting aggrecans with a molecular density functional theory that incorporates the acid-base equilibrium as well as the molecular properties: including conformations, size, shape, and charge distribution of all molecular species. The effect of acidity and salt concentration on the behavior is explored in detail. The repulsive interactions between two cylindrical-shaped aggrecans are strongly influenced by both the salt concentration and the pH. With increasing acidity, the polyelectrolytes of the aggrecan acquire charge and with decreasing salt concentration those charges become less screened. Consequently the interactions increase in size and range with increasing acidity and decreasing salt concentration. The size and range of the forces offers a possible explanation to the aggregation behavior of aggrecans and for their ability to resist compressive forces in cartilage. Likewise, the interdigitation of two aggrecan molecules is strongly affected by the salt concentration as well as the pH. With increasing pH, the number of charges increases, causing the repulsions between the polymers to increase, leading to a lower interdigitation of the two cylindrical polymer layers of the aggrecan molecules. The low interdigitation in charged polyelectrolytes layers provides an explanation for the good lubrication properties of polyelectrolyte layers in general and cartilage in particular.

Cartilage aggrecan can undergo self-adhesion

Here it is reported that aggrecan, the highly negatively charged macromolecule in the cartilage extracellular matrix, undergoes Ca(2+)-mediated self-adhesion after static compression even in the presence of strong electrostatic repulsion in physiological-like solution conditions. Aggrecan was chemically end-attached onto gold-coated planar silicon substrates and gold-coated microspherical atomic force microscope probe tips (end radius R approximately 2.5 mum) at a density ( approximately 40 mg/mL) that simulates physiological conditions in the tissue ( approximately 20-80 mg/mL). Colloidal force spectroscopy was employed to measure the adhesion between opposing aggrecan monolayers in NaCl (0.001-1.0 M) and NaCl + CaCl(2) ([Cl(-)] = 0.15 M, [Ca(2+)] = 0 - 75 mM) aqueous electrolyte solutions. Aggrecan self-adhesion was found to increase with increasing surface equilibration time upon compression (0-30 s). Hydrogen bonding and physical entanglements between the chondroitin sulfate-glycosaminoglycan side chains are proposed as important factors contributing to aggrecan self-adhesion. Self-adhesion was found to significantly increase with decreasing bath ionic strength (and hence, electrostatic double-layer repulsion), as well as increasing Ca(2+) concentration due to the additional ion-bridging effects. It is hypothesized that aggrecan self-adhesion, and the macromolecular energy dissipation that results from this self-adhesion, could be important factors contributing to the self-assembled architecture and integrity of the cartilage extracellular matrix in vivo.