Insights on the conformational properties of hyaluronic acid by using NMR residual dipolar couplings and MD simulations

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.

NMR Characterization of the Interactions Between Glycosaminoglycans and Proteins

Glycosaminoglycans (GAGs) constitute a considerable fraction of the glycoconjugates found on cellular membranes and in the extracellular matrix of virtually all mammalian tissues. The essential role of GAG-protein interactions in the regulation of physiological processes has been recognized for decades. However, the underlying molecular basis of these interactions has only emerged since 1990s. The binding specificity of GAGs is encoded in their primary structures, but ultimately depends on how their functional groups are presented to a protein in the three-dimensional space. This review focuses on the application of NMR spectroscopy on the characterization of the GAG-protein interactions. Examples of interpretation of the complex mechanism and characterization of structural motifs involved in the GAG-protein interactions are given. Selected families of GAG-binding proteins investigated using NMR are also described.

Computational and experimental approaches for chitosan-based nano PECs design: Insights on a deeper comprehension of nanostructure formation

Nanostructured polyelectrolyte complexes (nano PECs) based on biopolymers are an important technological strategy to target drugs to the action and/or absorption site in a more effective way. In this work, computational studies were performed to predict the ionization, spatial arrangement and interaction energies of chitosan (CS), hyaluronic acid (HA), and hypromellose phthalate (HP), for the design of nano PEC carriers for methotrexate (MTX). The optimal pH range (5.0-5.5) for preparing nano PECs was selected by experimental and computational methodologies, favoring the polymers interactions. CS, HA, HP and MTX addition order was also rationalized, maximizing their interactions and MTX entrapment. Spherical nano-sized particles (256-575 nm, by dynamic light scattering measurement) with positive surface charge (+25.5 to +29.2 mV) were successfully prepared. The MTX association efficiency ranged from 20 to 32 %. XRD analyses evidenced the formation of a new material with an organized structure, in relation to raw polymers.

Constructing 3-Dimensional Atomic-Resolution Models of Nonsulfated Glycosaminoglycans with Arbitrary Lengths Using Conformations from Molecular Dynamics

Glycosaminoglycans (GAGs) are the linear carbohydrate components of proteoglycans (PGs) and are key mediators in the bioactivity of PGs in animal tissue. GAGs are heterogeneous, conformationally complex, and polydisperse, containing up to 200 monosaccharide units. These complexities make studying GAG conformation a challenge for existing experimental and computational methods. We previously described an algorithm we developed that applies conformational parameters (i.e., all bond lengths, bond angles, and dihedral angles) from molecular dynamics (MD) simulations of nonsulfated chondroitin GAG 20-mers to construct 3-D atomic-resolution models of nonsulfated chondroitin GAGs of arbitrary length. In the current study, we applied our algorithm to other GAGs, including hyaluronan and nonsulfated forms of dermatan, keratan, and heparan and expanded our database of MD-generated GAG conformations. Here, we show that individual glycosidic linkages and monosaccharide rings in 10- and 20-mers of hyaluronan and nonsulfated dermatan, keratan, and heparan behave randomly and independently in MD simulation and, therefore, using a database of MD-generated 20-mer conformations, that our algorithm can construct conformational ensembles of 10- and 20-mers of various GAG types that accurately represent the backbone flexibility seen in MD simulations. Furthermore, our algorithm efficiently constructs conformational ensembles of GAG 200-mers that we would reasonably expect from MD simulations.

Effect of aggrecan degradation on the nanomechanics of hyaluronan in extra-fibrillar matrix of annulus fibrosus: A molecular dynamics investigation

Intervertebral Disc (IVD) Degeneration is one of the primary causes of low back pain among the adult population - the most significant cause being the degradation of aggrecan present in the extra-fibrillar matrix (EFM). Aggrecan degradation is closely associated with loss of water content leading to an alteration in the mechanical behaviour of the IVD. The loss in water content has a significant impact on the chemo-mechanical interplay of IVD biochemical constituents at the fundamental level. This work presents a mechanistic understanding of the effect of hydration, closely associated with aggrecan degradation, on the nanoscale mechanical behaviour of the hyaluronan present in the EFM of the Annulus Fibrosus. For this purpose, explicit three-dimensional molecular dynamics analyses of tensile and compressive tests are performed on a representative atomistic model of the hyaluronan present in the EFM. To account for the degradation of aggrecan, hydration levels are varied from 0 to 75% by weight of water. Analyses show that an increase in the hydration levels decreases the elastic modulus of hyaluronan in tension from ~4.6 GPa to ~2.1 GPa. On the other hand, the increase in hydration level increases the elastic moduli in axial compression from ~1.6 GPa in un-hydrated condition to ~6 GPa in 50% hydrated condition. But as the hydration levels increase to 75%, the elastic modulus reduces to ~3.5 GPa signifying a shift in load-bearing characteristic, from the solid hyaluronan component to the fluid component. Furthermore, analyses show a reduction in the intermolecular energy between hyaluronan and water, under axial tensile loading, indicating a nanoscale intermolecular debonding between hyaluronan and water molecules. This is attributed to the ability of hyaluronan to form stabilizing intra-molecular hydrogen bonds between adjacent residues. Compressive loading, on the other hand, causes intensive coiling of hyaluronan molecule, which traps more water through hydrogen bonding and aids in bearing compressive loads. Overall, study shows that hydration level has a strong influence on the atomistic level interactions between hyaluronan molecules and hyaluronan and water molecules in the EFM which influences the nanoscale mechanics of the Annulus Fibrosus.

The role of hydrophobic modification on hyaluronic acid dynamics and self-assembly

The advent of nanomedicine has rejuvenated the need for increased understanding of the fundamental physicochemical properties of polymeric amphiphiles. Hyaluronic acid (HA) is a hydrophilic polysaccharide that is frequently conjugated to hydrophobic moieties and then used to entrap dyes and therapeutics. Here, we develop computational models to examine the effects of the hydrophobic modification on supramolecular behavior among three systematically designed HA derivatives substituted with alkyl chains of increasing length. Our simulations coalesce with experimentally obtained results to demonstrate the dependence of supramolecular behavior on intramolecular forces. We show that the formation of clearly defined hydrophobic domains in samples of octadecylamine-modified HA compared to HA conjugates with shorter alkyl chains is a result of more favorable hydrophobic interactions. Trends in hydrodynamic radius and polydispersity are observed in experimental results that coalesce with theoretical calculations, suggesting that supramolecular properties are dependent on the physicochemical characteristics of individual polymer strands.

Computational Chemistry Tools in Glycobiology: Modelling of Carbohydrate–Protein Interactions

Molecular modelling provides a major impact in the field of glycosciences, helping in the characterisation of the molecular basis of the recognition between lectins from pathogens and human glycoconjugates, and in the design of glycocompounds with anti-infectious properties. The conformational properties of oligosaccharides are complex, and therefore, the simulation of these properties is a challenging task. Indeed, the development of suitable force fields is required for the proper simulation of important problems in glycobiology, such as the interatomic interactions responsible for oligosaccharide and glycoprotein dynamics, including O-linkages in oligo- and polysaccharides, and N- and O-linkages in glycoproteins. The computational description of representative examples is discussed, herein, related to biologically active oligosaccharides and their interaction with lectins and other proteins, and the new routes open for the design of glycocompounds with promising biological activities.

Q.E.COSY: determining sign and size of small deuterium residual quadrupolar couplings using an extended E.COSY principle

Residual quadrupolar couplings contain important structural information comparable with residual dipolar couplings. However, the measurement of sign and size of especially small residual quadrupolar couplings is difficult. Here, we present an extension of the E.COSY principle to spin systems consisting of a Spin 1 coupled to a spin ½ nucleus, which allows the determination of the sign of the quadrupolar coupling of the Spin 1 nucleus relative to the heteronuclear coupling between the spins. The so-called Q.E.COSY approach is demonstrated with its sign-sensitivity using variable angle NMR, stretched gels and liquid crystalline phases applied to various CD and CD3 groups. Especially the sign-sensitive measurement of residual quadrupolar couplings that remain unresolved in conventional deuterium 1D spectra is shown.

Specific excitonic interactions in the aggregates of hyaluronic acid and cyanine dyes with different lengths of methine group

The absorption and circular dichroism (CD) spectra of three types of cyanine dyes with different lengths of methine group (3,3'-diethylthiadicarbocyanine iodide, DTDC; 3,3'-diethylthiacarbocyanine iodide, DTC; and 3,3'-diethylthiacyanine iodide, DTTHC) in an aqueous solution were compared with and without hyaluronic acid (HA), which has a helical structure. DTDC forms chiral H- and J-aggregates, whereas DTC and DTTHC are unable to form any aggregates. DTDC also forms H- and J-aggregates in the presence of sodium polyacrylate (PA) with a random-coil structure; however, the PA-DTDC aggregates exhibit no chirality. These results suggest that the chirality of HA-DTDC aggregates is induced by the helical structure of HA. In 2.4 vol% and 10 vol% methanol, HA-DTDC aggregates displayed different patterns of temperature dependence, whereas no aggregation was observed in 30 vol% methanol. The solubility of DTDC in a mixed solvent of water and methanol is generally improved by the addition of methanol, which prevents the aggregation of DTDC.

Review: Use of residual dipolar couplings to determine the structure of carbohydrates

Solution nuclear magnetic resonance spectroscopy is especially useful in the carbohydrate field. The measurement of residual dipolar couplings provides long-range structural information, a valuable complement for the structural study of carbohydrates either in its free form or in the bound state to proteins. They permit to deduce the geometry and the flexibility of the glycosidic linkages, which have a major influence on the conformation of carbohydrates and their overall shape. This article reviews the current application of the residual dipolar couplings methodology to carbohydrates.

Critical parameters for the analysis of anionic oligosaccharides by desorption electrospray ionization mass spectrometry

Sulfated oligosaccharides derived from glycosaminoglycans (GAGs) are fragile compounds, highly polar and anionic. We report here on the rare but successful application of desorption electrospray ionization (DESI) - LTQ-Orbitrap mass spectrometry (MS) to the high-resolution analysis of anionic and sulfated oligosaccharides derived from the GAGs hyaluronic acid and heparin. For that purpose, key parameters affecting DESI performance, comprising the geometric parameters of the DESI source, the probed surface and the spraying conditions, applied spray voltage, flow rates and solvent composition were investigated. Under suitable conditions, the DESI technique allows the preservation of the structural integrity of such fragile compounds. DESI enabled the sensitive detection of anionic hyaluronic acid and heparin oligosaccharides with a limit of detection (LOD) down to 5 fmol (≈10 pg) for the hyaluronic acid decasaccharide. Detection of hyaluronic acid oligosaccharides in urine sample was also successfully achieved with LOD values inferior to the ng range. Multistage tandem mass spectrometry (MS(n) ) through the combination of the DESI source with a hybrid linear ion trap-orbitrap mass spectrometer allowed the discrimination of isomeric sulfated oligosaccharides and the sequence determination of a hyaluronic acid decasaccharide. These results open promising ways in glycomic and glycobiology fields where structure-activity relationships of bioactive carbohydrates are currently questioned.

Molecular conformations in the pentasaccharide LNF-1 derived from NMR spectroscopy and molecular dynamics simulations

The conformational dynamics of the human milk oligosaccharide lacto-N-fucopentaose (LNF-1), α-L-Fucp-(1 → 2)-β-D-Galp-(1 → 3)-β-D-GlcpNAc-(1 → 3)-β-D-Galp-(1 → 4)-D-Glcp, has been analyzed using NMR spectroscopy and molecular dynamics (MD) computer simulations. Employing the Hadamard (13)C-excitation technique and the J-HMBC experiment, (1)H,(13)C trans-glycosidic J coupling constants were obtained, and from one- and two-dimensional (1)H,(1)H T-ROESY experiments, proton-proton cross-relaxation rates were determined in isotropic D(2)O solution. In the lyotropic liquid-crystalline medium consisting of ditetradecylphosphatidylcholine, dihexylphosphatidylcholine, N-cetyl-N,N,N-trimethylammonium bromide, and D(2)O, (1)H, (1)H and one-bond (1)H, (13)C residual dipolar couplings (RDCs), as well as relative sign information on homonuclear RDCs, were determined for the pentasaccharide. Molecular dynamics simulations with explicit water were carried out from which the internal isomerization relaxation time constant, τ(N), was calculated for transitions at the ψ torsion angle of the β-(1 → 3) linkage to the lactosyl group in LNF-1. Compared to the global reorientation time, τ(M), of ∼0.6 ns determined experimentally in D(2)O solution, the time constant for the isomerization relaxation process, τ(N(scaled)), is about one-third as large. The NMR parameters derived from the isotropic solution show very good agreement with those calculated from the MD simulations. The only notable difference occurs at the reducing end, which should be more flexible than observed by the molecular simulation, a conclusion in complete agreement with previous (13)C NMR relaxation data. A hydrogen-bond analysis of the MD simulation revealed that inter-residue hydrogen bonds on the order of ∼30% were present across the glycosidic linkages to sugar ring oxygens. This finding highlights that intramolecular hydrogen bonds might be important in preserving well-defined structures in otherwise flexible molecules. An analysis including generalized order parameters obtained from nuclear spin relaxation experiments was performed and successfully shown to limit the conformational space accessible to the molecule when the number of experimental data are too scarce for a complete conformational analysis.