NMR solution conformation of heparin-derived tetrasaccharide

Differential Solvent DEEP-STD NMR and MD Simulations Enable the Determinants of the Molecular Recognition of Heparin Oligosaccharides by Antithrombin to Be Disentangled

The interaction of heparin with antithrombin (AT) involves a specific sequence corresponding to the pentasaccharide GlcNAc/NS6S-GlcA-GlcNS3S6S-IdoA2S-GlcNS6S (AGA*IA). Recent studies have revealed that two AGA*IA-containing hexasaccharides, which differ in the sulfation degree of the iduronic acid unit, exhibit similar binding to AT, albeit with different affinities. However, the lack of experimental data concerning the molecular contacts between these ligands and the amino acids within the protein-binding site prevents a detailed description of the complexes. Differential epitope mapping (DEEP)-STD NMR, in combination with MD simulations, enables the experimental observation and comparison of two heparin pentasaccharides interacting with AT, revealing slightly different bound orientations and distinct affinities of both glycans for AT. We demonstrate the effectiveness of the differential solvent DEEP-STD NMR approach in determining the presence of polar residues in the recognition sites of glycosaminoglycan-binding proteins.

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.

Designing "high-affinity, high-specificity" glycosaminoglycan sequences through computerized modeling

The prediction of high-affinity and/or high-specificity protein-glycosaminoglycan (GAG) interactions is an inherently difficult task, due to several factors including the shallow nature of the typical GAG-binding site and the inherent size, flexibility, diversity, and polydisperse nature of the GAG molecules. Here, we present a generally applicable methodology termed Combinatorial Library Virtual Screening (CVLS) that can identify potential high-affinity, high-specificity protein-GAG interactions from very large GAG combinatorial libraries and a suitable GAG-binding protein. We describe the CVLS approach along with the rationale behind it and provide validation for the method using the well-known antithrombin-thrombin-heparin system.

Structure and mechanism of human UDP-xylose synthase: evidence for a promoting role of sugar ring distortion in a three-step catalytic conversion of UDP-glucuronic acid

UDP-xylose synthase (UXS) catalyzes decarboxylation of UDP-d-glucuronic acid to UDP-xylose. In mammals, UDP-xylose serves to initiate glycosaminoglycan synthesis on the protein core of extracellular matrix proteoglycans. Lack of UXS activity leads to a defective extracellular matrix, resulting in strong interference with cell signaling pathways. We present comprehensive structural and mechanistic characterization of the human form of UXS. The 1.26-Å crystal structure of the enzyme bound with NAD⁺ and UDP reveals a homodimeric short-chain dehydrogenase/reductase (SDR), belonging to the NDP-sugar epimerases/dehydratases subclass. We show that enzymatic reaction proceeds in three chemical steps via UDP-4-keto-d-glucuronic acid and UDP-4-keto-pentose intermediates. Molecular dynamics simulations reveal that the d-glucuronyl ring accommodated by UXS features a marked ⁴C₁chair to B_O,3boat distortion that facilitates catalysis in two different ways. It promotes oxidation at C₄ (step 1) by aligning the enzymatic base Tyr¹⁴⁷ with the reactive substrate hydroxyl and it brings the carboxylate group at C₅ into an almost fully axial position, ideal for decarboxylation of UDP-4-keto-d-glucuronic acid in the second chemical step. The protonated side chain of Tyr¹⁴⁷ stabilizes the enolate of decarboxylated C₄ keto species (²H₁half-chair) that is then protonated from the Si face at C₅, involving water coordinated by Glu¹²⁰. Arg²⁷⁷, which is positioned by a salt-link interaction with Glu¹²⁰, closes up the catalytic site and prevents release of the UDP-4-keto-pentose and NADH intermediates. Hydrogenation of the C₄ keto group by NADH, assisted by Tyr¹⁴⁷ as catalytic proton donor, yields UDP-xylose adopting the relaxed ⁴C₁chair conformation (step 3).

Heparin characterization: challenges and solutions

Although heparin is an important and widely prescribed pharmaceutical anticoagulant, its high degree of sequence microheterogeneity and size polydispersity make molecular-level characterization challenging. Unlike nucleic acids and proteins that are biosynthesized through template-driven assembly processes, heparin and the related glycosaminoglycan heparan sulfate are actively remodeled during biosynthesis through a series of enzymatic reactions that lead to variable levels of O- and N-sulfonation and uronic acid epimers. As summarized in this review, heparin sequence information is determined through a bottom-up approach that relies on depolymerization reactions, size- and charge-based separations, and sensitive mass spectrometric and nuclear magnetic resonance experiments to determine the structural identity of component oligosaccharides. The structure-elucidation process, along with its challenges and opportunities for future analytical improvements, is reviewed and illustrated for a heparin-derived hexasaccharide.

NMR and amber analysis of the neamine pharmacophore for the design of novel aminoglycoside antibiotics

The dependence of the solution structure of neamine on pH was determined by NMR and AMBER molecular dynamics methods at pD 3.3, pD 6.5, and pD 7.4 in D(2)O at 25°C. Unlike neamine structures at pD 3.3 and 6.5, which essentially showed only one conformer, slowly exchanging primary, P-state, and secondary, S-state, neamine conformers populated on the NMR time scale at ~80% and ~20%, respectively, were detected at pD 7.4 with kinetic constants k(on(P→S))=1.9771s(-1) and k(off(S→P))=1.1319s(-1). A tertiary, T-state, neamine species populated at ~3% was also detected by NMR at pD 7.4. The pKa values determined by NMR titration experiments are pKa1 6.44±0.13 for N3 of ring-II, pKa2 7.23±0.09 for N2' of ring-I, pKa3 7.77±0.19 for N1 of ring-II, and pKa4 8.08±0.15 for N6' of ring-I. Ring-I and ring-II of the P-state neamine and ring-I of the S and T-states of neamine possess the (4)C(1) chair conformation between pD 3.3 and pD=7.4. In contrast, ring-II of the S and T-states of neamine most likely adopt the (6)rH(1) half-chair conformation. The P and S-states of neamine exhibit a negative syn-ψ glycosidic geometry. The exocyclic aminomethyl group of ring-I adopts the gt exocyclic rotamer conformation around physiological pHs while the gg exocyclic rotamer conformation predominates in acidic solutions near and below pH 4.5. Neamine exists in the P-state as a mixture of tetra-/tri-/di-protonated species between pD 4.5 and pD 7.4, while the S-state neamine exist only in a di-protonated species around physiological pDs. The existence of the S-state neamine may facilitate binding of neamine-like aminoglycosides by favorable entropy of binding to the A-site of 16S ribosomal RNA, suggesting that novel aminoglycoside compounds carrying a S-state neamine pharmacophore can be developed.

NMR methods to monitor the enzymatic depolymerization of heparin

Heparin and the related glycosaminoglycan, heparan sulfate, are polydisperse linear polysaccharides that mediate numerous biological processes due to their interaction with proteins. Because of the structural complexity and heterogeneity of heparin and heparan sulfate, digestion to produce smaller oligosaccharides is commonly performed prior to separation and analysis. Current techniques used to monitor the extent of heparin depolymerization include UV absorption to follow product formation and size exclusion or strong anion exchange chromatography to monitor the size distribution of the components in the digest solution. In this study, we used ¹H nuclear magnetic resonance (NMR) survey spectra and NMR diffusion experiments in conjunction with UV absorption measurements to monitor heparin depolymerization using the enzyme heparinase I. Diffusion NMR does not require the physical separation of the components in the reaction mixture and instead can be used to monitor the reaction solution directly in the NMR tube. Using diffusion NMR, the enzymatic reaction can be stopped at the desired time point, maximizing the abundance of larger oligosaccharides for protein-binding studies or completion of the reaction if the goal of the study is exhaustive digestion for characterization of the disaccharide composition. In this study, porcine intestinal mucosa heparin was depolymerized using the enzyme heparinase I. The unsaturated bond formed by enzymatic cleavage serves as a UV chromophore that can be used to monitor the progress of the depolymerization and for the detection and quantification of oligosaccharides in subsequent separations. The double bond also introduces a unique multiplet with peaks at 5.973, 5.981, 5.990, and 5.998 ppm in the ¹H-NMR spectrum downfield of the anomeric region. This multiplet is produced by the proton of the C-4 double bond of the non-reducing end uronic acid at the cleavage site. Changes in this resonance were used to monitor the progression of the enzymatic digestion and compared to the profile obtained from UV absorbance measurements. In addition, in situ NMR diffusion measurements were explored for their ability to profile the different-sized components generated over the course of the digestion.

Figure

Conformational degeneracy restricts the effective information content of heparan sulfate

The linear, sulfated polysaccharide heparan sulfate occupies a pivotal position in intercellular signalling events, interacting with numerous proteins on the cell surface and in the extracellular matrix. Its complex sequences suggest high potential information content but, despite extensive efforts, a clear relationship between its substitution pattern and biological activity remains elusive. This results from technical limitations, compounded by attempts to correlate substitution pattern directly with activity without considering other conformational factors. For a series of systematically modified analogues of heparan sulfate, the relationship between substitution pattern and experimental (13)C NMR chemical shifts, which act as reporters of the presence of conformational change, particularly around the glycosidic linkages, was explored through chemometric analysis. From analysis of the experimental data it was evident that wide linkage variation arose from O-sulfation in iduronate and N-sulfation in glucosamine residues but, their effects were distinct, while 6-O-sulfation had much less impact. Models of saccharide sequences showed that the maximum spread of variation in glycosidic linkages occurred before maximum sequence diversity and revealed a highly degenerate system: a fraction of possible sequences is sufficient to provide diverse backbone conformations to satisfy particular protein binding requirements. The unique information content potentially available in HS sequences, defined ultimately by conformation, is vastly inferior to the potential sequence diversity.

Residual dipolar coupling investigation of a heparin tetrasaccharide confirms the limited effect of flexibility of the iduronic acid on the molecular shape of heparin

The solution conformation of a fully sulfated heparin-derived tetrasaccharide, I, was studied in the presence of a 4-fold excess of Ca²⁺. Proton–proton and proton–carbon residual dipolar couplings (RDCs) were measured in a neutral aligning medium. The order parameters of two rigid hexosamine rings of I were determined separately using singular value decomposition and ab initio structures of disaccharide fragments of I. The order parameters were very similar implying that a common order tensor can be used to analyze the structure of I. Using one order tensor, RDCs of both hexosamine rings were used as restraints in molecular dynamics simulations. RDCs of the inner iduronic acid were calculated for every point of the molecular dynamics trajectory. The fitting of the calculated RDCs of the two forms of the iduronic acid to the experimental values yielded a population of ¹C₄ and ²S_o conformers of iduronic acid that agreed well with the analysis based on proton–proton scalar coupling constants. The glycosidic linkage torsion angles in RDC-restrained molecular dynamics (MD) structures of I are consistent with the interglycosidic three-bond proton–carbon coupling constants. These structures also show that the shape of heparin is not affected dramatically by the conformational flexibility of the iduronic acid ring. This is in line with conclusions of previous studies based on MD simulations and the analysis of ¹H-¹H NOEs. Our work therefore demonstrates the effectiveness of RDCs in the conformational analysis of glycosaminoglycans.

Free energy calculations of glycosaminoglycan-protein interactions

Glycosaminoglycans (GAGs) are complex highly charged linear polysaccharides that have a variety of roles in biological processes. We report the first use of molecular dynamics (MD) free energy calculations using the MM/PBSA method to investigate the binding of GAGs to protein molecules, namely the platelet endothelial cell adhesion molecule 1 (PECAM-1) and annexin A2. Calculations of the free energy of the binding of heparin fragments of different sizes reveal the existence of a region of low GAG-binding affinity in domains 5-6 of PECAM-1 and a region of high affinity in domains 2-3, consistent with experimental data and ligand-protein docking studies. A conformational hinge movement between domains 2 and 3 was observed, which allows the binding of heparin fragments of increasing size (pentasaccharides to octasaccharides) with an increasingly higher binding affinity. Similar simulations of the binding of a heparin fragment to annexin A2 reveal the optimization of electrostatic and hydrogen bonding interactions with the protein and protein-bound calcium ions. In general, these free energy calculations reveal that the binding of heparin to protein surfaces is dominated by strong electrostatic interactions for longer fragments, with equally important contributions from van der Waals interactions and vibrational entropy changes, against a large unfavorable desolvation penalty due to the high charge density of these molecules.

Conformational analysis of a dermatan sulfate-derived tetrasaccharide by NMR, molecular modeling, and residual dipolar couplings

The solution conformation behavior of a dermatan-derived tetrasaccharide--Delta HexA-(1-->3)-GalNAc4S-beta-(1-->4)-IdoA-alpha-(1-->3)-red-GalNAc4S (S is a sulfate group)--has been explored by means of NMR spectroscopy, especially by NOE-based conformational analysis. The tetrasaccharide was present as four species, two of which are chemically different in the anomeric orientation of the reducing 2-deoxy-2-acetamido-galactose (red-GalNAc) residue, while the other two are the result of different conformations of the iduronic acid (IdoA) unit. The two alpha-beta-interconverting anomers were present in a 0.6:1 ratio. Ring conformations have been defined by analysis of (3)J(H,H) coupling constants and interresidual NOE contacts. Both 2-deoxy-2-acetamido-galactose (GalNAc) residues were found in the (4)C(1) chair conformation, the unsaturated uronic acid (Delta-Hex A) adopts a strongly predominant half-chair (1)H(2) conformation, while the IdoA residue exists either in the (1)C(4) chair or in the (2)S(0) skewed boat geometries, in a 4:1 ratio. There is a moderate flexibility of Phi and Psi torsions as suggested by nuclear Overhauser effects (NOEs), molecular modeling (MM), and molecular dynamics (MD) studies. This was further investigated by residual dipolar couplings (RDCs). One-bond C--H RDCs ((1)D(C,H)) and long-range H-H ((3)D(H,H)) RDCs were measured for the tetrasaccharide in a phage solution and interpreted in combination with restrained MD simulation. The RDC-derived data substantially confirmed the validity of the conformer distribution resulting from the NOE-derived simulations, but allowed an improved definition of the conformational behavior of the oligosaccharides in solution. In summary, the data show a moderate flexibility of the four tetrasaccharide species at the central glycosidic linkage. Differences in the shapes of species with the IdoA in skew and in chair conformations and in the distribution of the sulfate groups have also been highlighted.

Structures of heparin-derived tetrasaccharide bound to cobra cardiotoxins: heparin binding at a single protein site with diverse side chain interactions

Cobra cardiotoxins (CTXs) are three-fingered polypeptides with positively charged domains that have been shown to bind to anionic ligands of snake venom citrate, glycosaminoglycans, sulfoglycosphingolipid, and nucleotide triphosphate with various biochemical effects including toxin dimerization, cell surface retention, membrane pore formation, cell internalization and blocking of enzymatic activities of kinase and ATPase. The reported anionic binding sites, however, are found to be different among different CTX homologues for potentially different CTX activities. Herein, by NMR studies of the binding of inorganic phosphate, dATP (stable form of ATP), and heparin-derived tetrasaccharide to Naja atra CTX A1, a novel CTX molecule exhibiting in vivo necrotic activity on skeletal muscle, we demonstrate that diverse ligands binding to CTXs could also occur at a single protein site with flexible side chain interactions. The flexibility of such an interaction is also illustrated by the available heparin-CTX A3 complex structures with different heparin chain lengths binding at the same site. Our results provide a likely structural explanation on how the interaction between heparan sufate and proteins depends more on the overall charge cluster organization rather than on their fine structures. We also suggest that the ligand binding site of CTX homologues can be fine-tuned by nonconserved residues near the binding pocket because of their flexible side chain interaction and dimerization ability, even for the rigid CTX molecules tightened by four disulfide bonds.

DMT-MM mediated functionalisation of the non-reducing end of glycosaminoglycans

Efficient functionalisation of the non-reducing end of uronic acid derivatives and glycosaminoglycan-derived disaccharides using peptide coupling has been achieved, mediated by the water-soluble agent DMT-MM.

Conformational Behavior of Basic Monomeric Building Units of Glycosaminoglycans: Isolated Systems and Solvent Effect

Our work reports in detail the results of systematic large-scale theoretical investigations of the glycosaminoglycan building units 1-OMe DeltaIdoA-2SNa2 (1; 2H1 and 1H2 forms), 1-OMe GlcN-S6SNa2 (2), 1,4-DiOMe GlcNa (3), 1,4-DiOMe GlcN-S3S6SNa3 (4), 1,4-DiOMe IdoA-2SNa2 (5; 4C1, 1C4, and 2So conformations), and 1,4-DiOMe GlcN-S6SNa2 (6) at the BP86/TZ2P level of the density functional theory. The optimized geometries indicate that the most stable structure of these monomeric units in the neutral state is stabilized via "multifurcated" sodium bonds. The equilibrium structure of the biologically active anionic forms of the glycosaminoglycans studied changed considerably in a water solution. The computed interaction energies, DeltaE, of sodium coordinated systems 1-6 are negative and span a rather broad energy interval (from -130 to -590 kcal mol-1). Computations that include the effect of solvation indicated that in water the relative stability of Na+...ligand ionic bonds is considerably diminished. The computed interaction energy in water is small (from -20 to -53 kcal mol-1) and negative, that is, stabilizing.

An investigation of the pyranose ring interconversion path of α-l-idose calculated using density functional theory

The interconversion pathways of the pyranose ring conformation of alpha-L-idose from a (4)C1 chair to other conformations were investigated using density functional calculations. From these calculations, four different ring interconversion paths and their transition state structures from the (4)C1 chair to other conformations, such as B(3,O), and (1)S3, were obtained. These four transition-state conformations cover four possible combinations of the network patterns of the hydroxyl group hydrogen bonds (clockwise and counterclockwise) and the conformations of the primary alcohol group (tg and gg). The optimized conformations, transition states, and their intrinsic reaction coordinates (IRC) were all calculated at the B3LYP/6-31G** level. The energy differences among the structures obtained were evaluated at the B3LYP/6-311++G** level. The optimized conformations indicate that the conformers of (4)C1, (2)S(O), and B(3,O) have similar energies, while (1)S3 has a higher energy than the others. The comparison of the four transition states and their ring interconversion paths, which were confirmed using the IRC calculation, suggests that the most plausible ring interconversion of the alpha-L-idopyranose ring occurs between (4)C1 and B(3,O) through the E3 envelope, which involves a 5.21 kcal/mol energy barrier.

Gas-Phase and Solution Conformations of the α-L-Iduronic Acid Structural Unit of Heparin

The IdoA2S structural unit of heparin (subunit G) may oscillate among the three conformations (4C1, 1C4, and 2So). Only the twisted boat conformation allowed the biologically active pentasaccharide unit of heparin (DEFGH) to bind to antithrombin. Our work reports, in detail, the results of systematic large-scale theoretical investigations of the three basic conformations (4C1, 1C4, and 2So) of the IdoA2S structural unit of heparin, its anionic forms, and its sodium salt using the B3LYP/6-311++G(d, p) and B3LYP/6-31+G(d) model chemistries. According to our calculations, the most stable structure of these molecules corresponds to the 2So skew-boat conformation. This form is also the most stable in a water solution. The 2So conformation of neutral molecules is not maintained in the anionic species. With anions, both 1C4 and 4C1 conformations are present. The relative stability of individual species of the substituted iduronic acid affects extra stabilization by means of intramolecular hydrogen bonds. The calculated macroscopic pKa of 1,4-DiOMe IdoA2S are as follows: pKa = 0.25 for the terminal C(2)-OSO3H group, pKa = 3.67 for the terminal C(5)-CO2H group, and pKa = 14.00 for the C(3)-OH hydroxyl group. The computed Gibbs interaction energies, DeltaGdegrees , for the reaction 1,4-DiOMe IdoA2S(2-) + 2Na+ <==> 1,4-DiOMe IdoA2SNa2 (4C1, 1C4, and 2So conformations) are negative and span a rather small energy interval (from -1244 to -1290 kJ mol(-1)).

Molecular dynamics and atomic charge calculations in the study of heparin conformation in aqueous solution

HF/6-31G** and molecular dynamics (MD) simulations were used to evaluate the performance of different atomic charge basis sets (i.e., Mulliken, Lowdin, and Electrostatic Potential Derived Charges--ESP) in heparin simulations. HF/3-21 G calculations were also used to study the NMR conformation of the IdoA residue. The results thus obtained indicated that ESP and Lowdin charges gave the better results in heparin simulations, followed by Mulliken charges, and that the minimum-energy conformation of IdoA can be different from that observed by NMR spectroscopy by less than 1 Angstrom. However, it was found that this small conformational modification is capable of inducing a change of almost 200 kJ/mol in the interactions of heparin with the surrounding environment, which is a meaningful amount of energy in the context of ligand-receptor interactions. This information can be potentially of great relevance in the design of heparin-derived antithrombotic compounds.

Dynamic properties of biologically active synthetic heparin-like hexasaccharides

A complete study of the dynamics of two synthetic heparin-like hexasaccharides, D-GlcNHSO3-6-SO4-alpha-(1-->4)-L-IdoA-2-SO4-alpha-(1-->4)-D-GlcNHSO3-6-SO4-alpha-(1-->4)-L-IdoA-2-SO4-alpha-(1-->4)-D-GlcNHSO3-6-SO4-alpha-(1-->4)-L-IdoA-2-SO4-alpha-1-->iPr (1) and -->4)-L-IdoA-2-SO4-alpha-(1-->4)-D-GlcNHAc-6-SO4-alpha-(1-->4)-L-IdoA-alpha-(1-->4)-D-GlcNHSO3-alpha-(1-->4)-L-IdoA-2-SO4-alpha-1-->iPr (2), has been performed using 13C-nuclear magnetic resonance (NMR) relaxation parameters, T1, T2, and heteronuclear nuclear Overhauser effect (NOEs). Compound 1 is constituted from sequences corresponding to the major polysaccharide heparin region, while compound 2 contains a sequence never found in natural heparin. They differ from each other only in sulphation patterns, and are capable of stimulating fibroblast growth factors (FGFs)-1 induced mitogenesis. Both oligosaccharides exhibit a remarkable anisotropic overall motion in solution as revealed by their anisotropic ratios (tau /tau||), 4.0 and 3.0 respectively. This is a characteristic behaviour of natural glycosaminoglycans (GAG) which has also been observed for the antithrombin (AT) binding pentasaccharide D-GlcNHSO3-6-SO4-alpha-(1-->4)-D-GlcA-beta-(1-->4)-D-GlcNHSO3-(3,6-SO4)-alpha-(1-->4)-L-IdoA-2-SO4-alpha-(1-->4)-D-GlcNHSO3-6-SO4-alpha-1-->Me (3) (Hricovíni, M., Guerrini, M., Torri, G., Piani, S., and Ungarelli, F. (1995) Conformational analysis of heparin epoxide in aqueous solution. An NMR relaxation study. Carbohydr. Res., 277, 11-23). The motional properties observed for 1 and 2 provide additional support to the suitability of these compounds as heparin models in agreement with previous structural (de Paz, J.L., Angulo, J., Lassaletta, J.M., Nieto, P.M., Redondo-Horcajo, M., Lozano, R.M., Jiménez-Gallego, G., and Martín-Lomas, M. (2001) The activation of fibroblast growth factors by heparin: synthesis, structure and biological activity of heparin-like oligosaccharides. Chembiochem, 2, 673-685; Ojeda, R., Angulo, J., Nieto, P.M., and Martin-Lomas. M. (2002) The activation of fibroblast growth factors by heparin: synthesis and structural study of rationally modified heparin-like oligosaccharides. Can. J. Chem,. 80, 917-936; Lucas, R., Angulo, J., Nieto, P.M., and Martin-Lomas, M. (2003) Synthesis and structural studies of two new heparin-like hexasaccharides. Org. Biomol. Chem., 1, 2253-2266) and biological data (Angulo, J., Ojeda, R., de Paz, J.L., Lucas, R., Nieto, P.M., Lozano, R.M., Redondo-Horcajo, M., Giménez-Gallego, G., and Martín-Lomas, M. (2004) The activation of fibroblast growth factors (FGFs) by glycosaminoglycans: influence of the sulphation pattern on the biological activity of FGF-1. Chembiochem, 5, 55-61). Fast internal motions observed for the less sulphated compound 2, as compared with 1, may be related to their different behavior in stimulating FGF1-induced mitogenic activity.

Conformation of glycosaminoglycans by ion mobility mass spectrometry and molecular modelling

We have performed conformational analyses of heparin-derived oligosaccharide ions in the gas phase using a combination of ion-mobility mass spectrometry and molecular modelling. Negative mode electrospray ionisation was used to generate singly (disaccharide, [C12H15NO19S3Na3]-) and doubly charged (tetrasaccharides, [C24H30N2O38S6Na6]2- and [C24H31N2O35S5Na5]2-) ions containing three and six Na+ ions, respectively. Good agreement was obtained between the experimental and theoretical cross sections. The latter were obtained using modelled structures generated by the AMBER-based force field. Analysis of the conformations of the oligosaccharide ions shows that sodium cations play a major role in stabilizing these ions in the gas phase. This was seen in the formation of oligomers of the disaccharide ion and "compact" structures of tetrasaccharide ions. Interestingly, the gas phase conformations of the three tetrasaccharide ions with different primary structures were significantly different.

Molecular dynamics simulation of a decasaccharide fragment of heparin in aqueous solution

Molecular dynamics (MD) simulations on heparin-water-sodium systems were carried out in order to establish a simulation protocol able to represent heparin solution conformation under physiological conditions. Atomic charges suitable for heparin oligosaccharides were obtained from ab initio quantum-mechanical computations, at the 6-31G(**) level. The GROMACS forcefield, the SPC, and SPC/E water models were employed. Also heparin was simulated with IdoA residues in 1C(4) or 2S(0) conformational states. The results of the performed MD simulations are in agreement with the available experimental data, suggesting that this approach can be applied for the study of heparin interactions with its target proteins and thus play a role in the development of new antithrombotic agents.

Heparin-binding domains in vascular biology

Heparin is a major anticoagulant with activity mediated primarily through its interaction with antithrombin (AT). Heparan sulfate (HS), structurally related to heparin, binds a wide range of proteins of different functionality, taking part in various physiological and pathological processes. The heparin-AT complex, the most well understood facet of anticoagulation, serves as a prototypical example of the important role of heparin/HS in vascular biology. Extensive studies have identified common structural features in heparin/HS-binding sites of proteins. These include the elucidation of consensus sequences in proteins, patterns of clusters of basic and nonbasic residues, and common spatial arrangements of basic amino acids in the heparin-binding sites. Although these studies have provided valuable information, heparin/HS-binding proteins differ widely in structure. The prediction of heparin/HS-binding proteins from sequence information is not currently possible, and elucidation of protein-binding sites requires the individual study of each glycosaminoglycan-protein complex. Thus, x-ray crystallography and site-directed mutagenesis experiments are among the most powerful tools, providing accurate structural information, facilitating the characterization of heparin-protein complexes. Heparin and structurally related heparan sulfate bind a large number of proteins, taking part in a wide range of biological processes, particularly ones involved in vascular biology. Heparin-binding domains share certain common structural features, but there is no absolute dependency on specific sequences or protein folds.

Synthesis and structural study of two new heparin-like hexasaccharides

Two new heparin-like hexasaccharides, 5 and 6, have been synthesised using a convergent block strategy and their solution conformations have been determined by NMR spectroscopy and molecular modelling. Both hexasaccharides contain the basic structural motif of the regular region of heparin but with negative charge distributions which have been designed to get insight into the mechanism of fibroblast growth factors (FGFs) activation.

The activation of fibroblast growth factors by heparin: Synthesis and structural study of rationally modified heparin-like oligosaccharides

Heparin-like hexasaccharide 3 and octasaccharide 4 have been synthesized using a convergent block strategy and their solution conformations have been determined by NMR spectroscopy. Both oligosaccharides contain the basic structural motif of the regular region of heparin but have been constructed as to display negatively charged sulfate groups only on one side of their solution helical structures. This charge distribution along the saccharide chain has been designed to get insight into the proposed mechanism for fibroblast growth factors (FGFs) activation that involves heparin-induced FGF dimerization.Key words: heparin oligosaccharides, synthesis design, conformational analysis, FGF activation.

Heparin-protein interactions

Heparin, a sulfated polysaccharide belonging to the family of glycosaminoglycans, has numerous important biological activities, associated with its interaction with diverse proteins. Heparin is widely used as an anticoagulant drug based on its ability to accelerate the rate at which antithrombin inhibits serine proteases in the blood coagulation cascade. Heparin and the structurally related heparan sulfate are complex linear polymers comprised of a mixture of chains of different length, having variable sequences. Heparan sulfate is ubiquitously distributed on the surfaces of animal cells and in the extracellular matrix. It also mediates various physiologic and pathophysiologic processes. Difficulties in evaluating the role of heparin and heparan sulfate in vivo may be partly ascribed to ignorance of the detailed structure and sequence of these polysaccharides. In addition, the understanding of carbohydrate-protein interactions has lagged behind that of the more thoroughly studied protein-protein and protein-nucleic acid interactions. The recent extensive studies on the structural, kinetic, and thermodynamic aspects of the protein binding of heparin and heparan sulfate have led to an improved understanding of heparin-protein interactions. A high degree of specificity could be identified in many of these interactions. An understanding of these interactions at the molecular level is of fundamental importance in the design of new highly specific therapeutic agents. This review focuses on aspects of heparin structure and conformation, which are important for its interactions with proteins. It also describes the interaction of heparin and heparan sulfate with selected families of heparin-binding proteins.

Synthesis of conformationally locked L-iduronic acid derivatives: direct evidence for a critical role of the skew-boat 2S0 conformer in the activation of antithrombin by heparin

We have used organic synthesis to understand the role of L-iduronic acid conformational flexibility in the activation of antithrombin by heparin. Among known synthetic analogues of the genuine pentasaccharidic sequence representing the antithrombin binding site of heparin, we have selected as a reference compound the methylated anti-factor Xa pentasaccharide 1. As in the genuine original fragment, the single L-iduronic acid moiety of this molecule exists in water solution as an equilibrium between three conformers 1C4, 4C1 and 2S0. We have thus synthesized three analogues of 1, in which the L-iduronic acid unit is locked in one of these three fixed conformations. A covalent two atom bridge between carbon atoms two and five of L-iduronic acid was first introduced to lock the pseudorotational itinerary of the pyranoid ring around the 2S0 form. A key compound to achieve this connection was the D-glucose derivative 5 in which the H-5 hydrogen atom has been replaced by a vinyl group, which is a progenitor of the carboxylic acid. Selective manipulations of this molecule resulted in the 2S0-type pentasaccharide 23. Starting from the D-glucose derivative 28, a covalent two atom bridge was now built up between carbon atoms three and five to lock the L-iduronic acid moiety around the 1C4 chair form conformation, and the 1C4-type pentasaccharide 43 was synthesized. Finally the L-iduronic acid containing disaccharide 58 which, due to the presence of the methoxymethyl substituent at position five adopts a 4C1 conformation, was directly used to synthesize the 4C1-type pentasaccharide 61. The locked pentasaccharide 23 showed about the same activity as the reference compound 1 in an antithrombin-mediated anti-Xa assay, whereas the two pentasaccharides 43 and 61 displayed very low activity. These results clearly establish the critical importance of the 2S0 conformation of L-iduronic acid in the activation of antithrombin by heparin.

Analytical method of heparan sulfates using high-performance liquid chromatography turbo-ionspray ionization tandem mass spectrometry

We established a highly sensitive quantitative analytical method of heparan sulfates (HS) by LC-MS-MS. It became possible to determine the unsaturated disaccharides produced by the enzyme digestion of HS, and to perform the whole analyses on one sample within 3 min by use of a short column of CAPCELL PAK NH2 UG80 (35 mm x 2 mm I.D.). The assay method was validated and showed the satisfactory sensitivity, precision and accuracy, which enabled the quantitation up to picomol level. By employing this method, we performed the analyses of HS in mouse brain and liver, and tumor tissues of tumor-bearing mouse transplanted subcutaneously with Meth A fibrosarcoma cells. The compositions of the unsaturated disaccharide units derived from HS were found to be somewhat different among those tissues. It is assumed that the site of sulfation in HS may be controlled by certain regulatory mechanisms. The quantitative method developed in this study is believed to be a very useful method for the determination of compositional profiles of constitutive disaccharide units of tissue HS.

Annexin V--heparin oligosaccharide complex suggests heparan sulfate--mediated assembly on cell surfaces

BACKGROUND

Annexin V, an abundant anticoagulant protein, has been proposed to exert its effects by self-assembling into highly ordered arrays on phospholipid membranes to form a protective anti-thrombotic shield at the cell surface. The protein exhibits very high-affinity calcium-dependent interactions with acidic phospholipid membranes, as well as specific binding to glycosaminoglycans (GAGs) such as heparin and heparan sulfate, a major component of cell surface proteoglycans. At present, there is no structural information to elucidate this interaction or the role it may play in annexin V function at the cell surface.

RESULTS

We report the 1.9 A crystal structure of annexin V in complex with heparin-derived tetrasaccharides. This structure represents the first of a heparin oligosaccharide binding to a protein where calcium ions are essential for the interaction. Two distinct GAG binding sites are situated on opposite protein surfaces. Basic residues at each site were identified from the structure and site-directed mutants were prepared. The heparin binding properties of these mutants were measured by surface plasmon resonance. The results confirm the roles of these mutated residues in heparin binding, and the kinetic and thermodynamic data define the functionally distinct character of each distal binding surface.

CONCLUSION

The annexin V molecule, as it self-assembles into an organized array on the membrane surface, can bind the heparan sulfate components of cell surface proteoglycans. A novel model is presented in which proteoglycan heparan sulfate could assist in the localization of annexin V to the cell surface membrane and/or stabilization of the entire molecular assembly to promote anticoagulation.

An NMR and molecular modeling study of the site-specific binding of histamine by heparin, chemically modified heparin, and heparin-derived oligosaccharides

The diprotonated form of histamine binds site-specifically to heparin, a highly sulfated 1-->4 linked repeating copolymer comprised predominantly of 2-O-sulfo-alpha-L-iduronic acid (the I ring) and 2-deoxy-2-sulfamido-6-O-sulfo-alpha-D-glucopyranosyl (the A ring). The binding is mediated by electrostatic interactions. The structural features of histamine and heparin, which are required for the site-specific binding, have been identified from the results of (1)H NMR studies of the binding of histamine by six heparin-derived oligosaccharides and four chemically modified heparins and molecular modeling studies. The results indicate that the imidazolium ring of diprotonated histamine is critical for directing site-specific binding, while the ammonium group increases the binding affinity. The imidazolium ring binds within a cleft, with the A ring of an IAI triad at the top of the cleft, and the I rings forming the two sides. The H3 proton of the A ring is in the shielding cone of the imidazolium ring. The carboxylate group of the I-ring at the reducing end of the IAI triad and possibly the sulfamido group of the A-ring are essential for site-specific binding, whereas the 2-O-sulfate group of the I ring and the 6-O-sulfate group of the A ring are not. The results indicate that histamine binds to the IAI triad with the I rings in the (1)C(4) conformation. Also, the configuration of the carboxylate group is critical, as indicated by the absence of site-specific binding of histamine by the related IAG sequence, where G is alpha-D-glucuronic acid. The molecular modeling results indicate that the N1H and N3H protons of the imidazolium ring of site-specifically bound histamine are hydrogen bonded to the carboxylates of the I rings at the nonreducing and reducing ends of the IAI trisaccharide sequence.

NMR solution conformation of heparin-derived hexasaccharide

The solution conformation of homogeneous, heparin-derived hexasaccharide (residues A, B, C, D, E, F) has been investigated by using 1H-NMR spectroscopy. Intra-ring conformations have been defined by J-coupling constants and inter-proton nuclear Overhauser effects (NOEs), and the orientation of one ring with respect to the other has been defined by inter-ring NOEs. NOE-based conformational modelling has been done by using the iterative relaxation matrix approach (IRMA), restrained energy minimization to refine structures and to distinguish between minor structural differences and equilibria between various intra-ring forms. All glucosamine residues B, D and F are in the 4C1 chair conformation. The uronate (A) residue is mostly represented by the 1H2 form, whereas internal iduronates (C and E) exist in equilibrium between the chair and skewed boat forms. Deviations in some NOEs indicate a minor contribution of the 2H1 form to the A ring. Glycosidic dihedral angles, which define the overall oligosaccharide conformation, were further refined by combining in vacuo energy map calculations and restrained energy minimization in explicit solvent water. Conformational stability was further assessed by subjecting NOE and IRMA-derived structures to 600 ps of unrestrained molecular dynamics in explicit solvent.