Molecular dynamics of a tetrasaccharide subunit of chondroitin 4-sulfate in water

A review on revolutionizing ophthalmic therapy: Unveiling the potential of chitosan, hyaluronic acid, cellulose, cyclodextrin, and poloxamer in eye disease treatments

Ocular disorders, encompassing both common ailments like dry eye syndrome and more severe situations for instance age-related macular degeneration, present significant challenges to effective treatment due to the intricate architecture and physiological barriers of the eye. Polysaccharides are emerging as potential solutions for drug delivery to the eyes due to their compatibility with living organisms, natural biodegradability, and adhesive properties. In this review, we explore not only the recent advancements in polysaccharide-based technologies and their transformative potential in treating ocular illnesses, offering renewed optimism for both patients and professionals but also anatomy of the eye and the significant obstacles hindering drug transportation, followed by an investigation into various drug administration methods and their ability to overcome ocular-specific challenges. Our focus lies on biological adhesive polymers, including chitosan, hyaluronic acid, cellulose, cyclodextrin, and poloxamer, known for their adhesive characteristics enhancing drug retention on ocular surfaces and increasing bioavailability. A detailed analysis of material designs used in ophthalmic formulations, such as gels, lenses, eye drops, nanofibers, microneedles, microspheres, and nanoparticles, their advantages and limitations, the potential of formulations in improving therapeutic outcomes for various eye conditions. Moreover, we underscore the discovery of novel polysaccharides and their potential uses in ocular drug delivery.

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.

Formation of Polyion Complex Aggregate Formed from a Cationic Block Copolymer and Anionic Polysaccharide

Block copolymers (PmMn; P20M101 and P100M98) comprising poly(2-(methacryloyloxy)ethylphosphorylcholine) (PMPC, P) containing biocompatible phosphorylcholin pendants and cationic poly((3-acryloylaminopropyl) trimethylammonium chloride) (PMAPTAC, M) were synthesized via a controlled radical polymerization method. The degrees of polymerization of the PMPC and PMAPTAC segments are denoted by subscripts (PmMn). The mixture of cationic PmMn and anionic sodium chondroitin sulfate C (CS) with the pendant anionic carboxylate and sulfonate groups formed polyion complex (PIC) aggregates in phosphate-buffered saline. A charge-neutralized mixture of P20M101 with CS formed P20M101/CS PIC vesicles with a hydrodynamic radius (Rh) of 97.2 nm, zeta potential of ca. 0 mV, and aggregation number (Nagg) of 23,044. PMPC shells covered the surface of the PIC vesicles. The mixture of P100M98 and CS formed PIC spherical micelles with the PIC core and hydrophilic PMPC shells. The Rh, zeta potential, and Nagg of the PIC micelles were 26.4 nm, ca. 0 mV, and 404, respectively. At pH < 4, the carboxylate anions in CS were protonated. Thus, the charge balance in the PIC micelles shifted to decrease the core density owing to the electrostatic repulsions of the excess cations in the core. The PIC micelles dissociated at a NaCl concentration ≥0.6 M owing to the charge screening effect. The positively charged PIC micelles with excess P100M98 can encapsulate anionic dyes owing to electrostatic interaction.

Sulfation Effects of Chondroitin Sulfate to Bind a Chemokine in Aqueous Medium: Conformational Heterogeneity and Dynamics from Molecular Simulation

The sulfation patterns and degree of sulfation of chondroitin sulfate (CS), an important class of glycosaminoglycans (GAG), and their interactions with chemokines are accountable for various diseases. To realize the underlying mechanism of such complex biological phenomena at a molecular level and their application in rational drug design, a study on conformations and dynamics of CSs is necessary. To explore this, in this study, we performed a series of atomistic molecular dynamics (MD) simulations with different sulfated variants of octadecasaccharide CS, like CS-C, CS-E, and CS-T, in their free forms and when bound to the protein chemokine CXCL8 dimer in an aqueous medium. The calculated binding free energy of CSs with the CXCL8 dimer is favorable, and the degree of sulfation favors the complexation process further with prominent hydrophobic and hydrogen-bonded interactions. We find that the recognition is associated with the configurational entropy loss of the CS molecules as calculated from the Gaussian mixture approach, which supports that the degree of sulfation regulates the process. Cluster analysis through the k-means algorithm and end-to-end distance measurement revealed that although the free CS molecules adopted linear conformations, the nonlinear conformations during binding with protein were noted. Adaptation of nonlinear forms in the bound forms is noteworthy for the less-sulfated CS-C and CS-E. Apart from favorable 4C1 conformations, the occasional appearance of skew-boat forms from the free-energy map of ring pucker for the GlcUA unit was observed, which remains unaffected by the sulfation. We find that during recognition, the average relaxation time of intra-CS and inter-CS-CXCL8 hydrogen bonds (HBs) is about a magnitude lesser than that of CS-water HBs, most prominent on the involvement of higher sulfated CS-T analogues. The translational motion of surrounded water molecules in CSs exhibited sublinear diffusion, and the degree of sublinearity increases around the heavily sulfated molecules due to the hindrance created by them as well as the presence of the chemokine and exhibited markedly slow heterogeneous diffusion.

Molecular Dynamics-Based Comparative Analysis of Chondroitin and Dermatan Sulfates

Glycosaminoglycans (GAGs) are a class of linear anionic periodic polysaccharides containing disaccharide repetitive units. These molecules interact with a variety of proteins in the extracellular matrix and so participate in biochemically crucial processes such as cell signalling affecting tissue regeneration as well as the onset of cancer, Alzheimer's or Parkinson's diseases. Due to their flexibility, periodicity and chemical heterogeneity, often termed "sulfation code", GAGs are challenging molecules both for experiments and computation. One of the key questions in the GAG research is the specificity of their intermolecular interactions. In this study, we make a step forward to deciphering the "sulfation code" of chondroitin sulfates-4,6 (CS4, CS6, where the numbers correspond to the position of sulfation in NAcGal residue) and dermatan sulfate (DS), which is different from CSs by the presence of IdoA acid instead of GlcA. We rigorously investigate two sets of these GAGs in dimeric, tetrameric and hexameric forms with molecular dynamics-based descriptors. Our data clearly suggest that CS4, CS6 and DS are substantially different in terms of their structural, conformational and dynamic properties, which contributes to the understanding of how these molecules can be different when they bind proteins, which could have practical implications for the GAG-based drug design strategies in the regenerative medicine.

In-Depth Molecular Dynamics Study of All Possible Chondroitin Sulfate Disaccharides Reveals Key Insight into Structural Heterogeneity and Dynamism

GAGs exhibit a high level of conformational and configurational diversity, which remains untapped in terms of the recognition and modulation of proteins. Although GAGs are suggested to bind to more than 800 biologically important proteins, very few therapeutics have been designed or discovered so far. A key challenge is the inability to identify, understand and predict distinct topologies accessed by GAGs, which may help design novel protein-binding GAG sequences. Recent studies on chondroitin sulfate (CS), a key member of the GAG family, pinpointing its role in multiple biological functions led us to study the conformational dynamism of CS building blocks using molecular dynamics (MD). In the present study, we used the all-atom GLYCAM06 force field for the first time to explore the conformational space of all possible CS building blocks. Each of the 16 disaccharides was solvated in a TIP3P water box with an appropriate number of counter ions followed by equilibration and a production run. We analyzed the MD trajectories for torsional space, inter- and intra-molecular H-bonding, bridging water, conformational spread and energy landscapes. An in-house phi and psi probability density analysis showed that 1→3-linked sequences were more flexible than 1→4-linked sequences. More specifically, phi and psi regions for 1→4-linked sequences were held within a narrower range because of intra-molecular H-bonding between the GalNAc O5 atom and GlcA O3 atom, irrespective of sulfation pattern. In contrast, no such intra-molecular interaction arose for 1→3-linked sequences. Further, the stability of 1→4-linked sequences also arose from inter-molecular interactions involving bridged water molecules. The energy landscape for both classes of CS disaccharides demonstrated increased ruggedness as the level of sulfation increased. The results show that CS building blocks present distinct conformational dynamism that offers the high possibility of unique electrostatic surfaces for protein recognition. The fundamental results presented here will support the development of algorithms that help to design longer CS chains for protein recognition.

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.

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.

Individual hippocampal subfield assessment indicates that matrix macromolecules and gliosis are key elements for the increased T2 relaxation time seen in temporal lobe epilepsy

OBJECTIVE

Increased T2 relaxation time is often seen in temporal lobe epilepsy (TLE) with hippocampal sclerosis. Water content directly affects the effective T2 in a voxel. Our aim was to evaluate the relation between T2 values and two molecules associated with brain water homeostasis aquaporin 4 (AQP4) and chondroitin sulfate proteoglycan (CSPG), as well as cellular populations in the hippocampal region of patients with TLE.

METHODS

Hippocampal T2 imaging and diffusion tensor imaging (DTI) were obtained from 42 drug-resistant patients with TLE and 20 healthy volunteers (radiologic controls, RCs). A similar protocol (ex vivo) was applied to hippocampal sections from the same TLE cases and 14 autopsy control hippocampi (histologic and radiologic controls, HRCs), and each hippocampal subfield was evaluated. Hippocampal sections from TLE cases and HRC controls were submitted to immunohistochemistry for neurons (neuron nuclei [NeuN]), reactive astrocytes (glial fibrillary acidic protein [GFAP]), activated microglia (human leukocyte antigen-D-related [HLA-DR]), polarized AQP4, and CSPG.

RESULTS

Patients with TLE had higher in vivo and ex vivo hippocampal T2 relaxation time. Hippocampi from epilepsy cases had lower neuron density, higher gliosis, decreased AQP4 polarization, and increased CSPG immunoreactive area. In vivo relaxation correlated with astrogliosis in the subiculum and extracellular CSPG in the hilus. Ex vivo T2 relaxation time correlated with astrogliosis in the hilus, CA4, and subiculum, and with microgliosis in CA1. The difference between in vivo and ex vivo relaxation ratio correlated with mean diffusivity and with the immunopositive area for CSPG in the hilus.

SIGNIFICANCE

Our data indicate that astrogliosis, microgliosis, and CSPG expression correlate with the increased T2 relaxation time seen in the hippocampi of patients with TLE.

Sulfation and cation effects on the conformational properties of the glycan backbone of chondroitin sulfate disaccharides

Chondroitin sulfate (CS) is one of several glycosaminoglycans that are major components of proteoglycans. A linear polymer consisting of repeats of the disaccharide -4GlcAβ1-3GalNAcβ1-, CS undergoes differential sulfation resulting in five unique sulfation patterns. Because of the dimer repeat, the CS glycosidic "backbone" has two distinct sets of conformational degrees of freedom defined by pairs of dihedral angles: (ϕ1, ψ1) about the β1-3 glycosidic linkage and (ϕ2, ψ2) about the β1-4 glycosidic linkage. Differential sulfation and the possibility of cation binding, combined with the conformational flexibility and biological diversity of CS, complicate experimental efforts to understand CS three-dimensional structures at atomic resolution. Therefore, all-atom explicit-solvent molecular dynamics simulations with Adaptive Biasing Force sampling of the CS backbone were applied to obtain high-resolution, high-precision free energies of CS disaccharides as a function of all possible backbone geometries. All 10 disaccharides (β1-3 vs β1-4 linkage × five different sulfation patterns) were studied; additionally, ion effects were investigated by considering each disaccharide in the presence of either neutralizing sodium or calcium cations. GlcAβ1-3GalNAc disaccharides have a single, broad, thermodynamically important free-energy minimum, whereas GalNAcβ1-4GlcA disaccharides have two such minima. Calcium cations but not sodium cations bind to the disaccharides, and binding is primarily to the GlcA -COO(-) moiety as opposed to sulfate groups. This binding alters the glycan backbone thermodynamics in instances where a calcium cation bound to -COO(-) can act to bridge and stabilize an interaction with an adjacent sulfate group, whereas, in the absence of this cation, the proximity of a sulfate group to -COO(-) results in two like charges being both desolvated and placed adjacent to each other and is found to be destabilizing. In addition to providing information on sulfation and cation effects, the present results can be applied to building models of CS polymers and as a point of comparison in studies of CS polymer backbone dynamics and thermodynamics.

The SPASIBA force field for chondroitin sulfate: vibrational analysis of D-glucuronic and N-acetyl-D-galactosamine 4-sulfate sodium salts

Normal-mode analyses were carried out on the two components of the chondroitin 4-sulfate linear glycosaminoglycan, a copolymer implying alternate D-glucuronate beta-(1-->3) and N-acetyl-D-galactosamine 4-sulfate beta-(1-->4) (hereafter named D-galactosamine 4-sulfate) residues. Scaled quantum mechanical calculations (SQM) using the density functional theory approach at different levels of theory (B3LYP/6-31G** and B3LYP/6-31++G**) were performed to obtain correct vibrational assignments. The SPASIBA empirical force field parameters were then obtained from both theoretical predictions and observed IR and Raman data. It is shown that calculations including diffuse functions at the B3LP/6-31++G** level and the introduction of the Na+ counterion are necessary to give correct assignments of the CO2- symmetric (nu(s)) and antisymmetric (nu(a)) stretching modes for the glucuronic carboxylate residue.

Conformational behavior of chondroitin and chondroitin sulfate in relation to their physical properties as inferred by molecular modeling

Chondroitin and chondroitin sulfates belong to the family of glycosaminoglycans. They are most widely distributed in animal tissues, where they are involved in structural functions and in cell-cell communication. Their basic structures consist of a disaccharidic repeating unit of beta-D-glucuronic acid (GlcA) and 2-acetamido-2-deoxy-beta-D-galactose (GalNAc), this latter being sulfated at different positions. Molecular mechanics has been applied to calculate the adiabatic energy maps for each of the constituting disaccharides of chondroitin, chondroitin 4-sulfate, and chondroitin 6-sulfate using the MM3 force field. Based on these maps, higher levels of structural organization have been simulated. On one hand, the disordered state is studied through a Metropolis-based algorithm; the resulting chains present a behavior of semirigid polymers, with an order of stiffness: chondroitin 4-sulfate > chondroitin > chondroitin 6-sulfate. On the other hand, the exploration of the stable ordered forms leads to numerous helical conformations of comparable energies. Several of these conformations correspond to the experimentally observed ones. The ability of coordination with cations has also been explored, resulting in a preferential stereospecificity for calcium ions when compared to sodium ions.