Conformational selection of the AGA*IA(M) heparin pentasaccharide when bound to the fibroblast growth factor receptor

Characterization of Heparin's Conformational Ensemble by Molecular Dynamics Simulations and Nuclear Magnetic Resonance Spectroscopy

Heparin is a highly charged, polysulfated polysaccharide and serves as an anticoagulant. Heparin binds to multiple proteins throughout the body, suggesting a large range of potential therapeutic applications. Although its function has been characterized in multiple physiological contexts, heparin's solution conformational dynamics and structure-function relationships are not fully understood. Molecular dynamics (MD) simulations facilitate the analysis of a molecule's underlying conformational ensemble, which then provides important information necessary for understanding structure-function relationships. However, for MD simulations to afford meaningful results, they must both provide adequate sampling and accurately represent the energy properties of a molecule. The aim of this study is to compare heparin's conformational ensemble using two well-developed force fields for carbohydrates, known as GLYCAM06 and CHARMM36, using replica exchange molecular dynamics (REMD) simulations, and to validate these results with NMR experiments. The anticoagulant sequence, an ultra-low-molecular-weight heparin, known as Arixtra (fondaparinux, sodium), was simulated with both parameter sets. The results suggest that GLYCAM06 matches experimental nuclear magnetic resonance three-bond J-coupling values measured for Arixtra better than CHARMM36. In addition, NOESY and ROESY experiments suggest that Arixtra is very flexible in the sub-millisecond time scale and does not adopt a unique structure at 25 C. Moreover, GLYCAM06 affords a much more dynamic conformational ensemble for Arixtra than CHARMM36.

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.

Synthesis of type 1 Lewis b hexasaccharide antigen structures featuring flexible incorporation of l-[U-13C6]-fucose for NMR binding studies

While 13C-labelled proteins are common tools in NMR studies, lack of access to 13C-labelled carbohydrate structures has restricted their use. l-Fucose is involved in a wide range of physiological and pathophysiological processes in mammalian organisms. Here, l-[U-13C6]-Fuc labelled type I Lewis b (Leb) structures have been synthesised for use in NMR binding studies with the Blood-group Antigen Binding Adhesin (BabA), a membrane-bound protein from the bacterium Helicobacter pylori. As part of this work, an efficient synthesis of a benzylated l-[U-13C6]-Fuc thioglycoside donor from l-[U-13C6]-Gal has been developed. The design and synthesis of an orthogonally protected tetrasaccharide precursor enabled controlled introduction of one or two 13C-labelled or non-labelled fucosyl residues prior to global deprotection. NMR analysis showed that it is straightforward to assign the anomeric centres as well as the H-5 positions to the individual fucosyl residues which are relevant for NMR binding studies.

Synthetic Oligosaccharide Libraries and Microarray Technology: A Powerful Combination for the Success of Current Glycosaminoglycan Interactomics

Glycosaminoglycans (GAGs) are extracellular matrix and/or cell-surface sulfated glycans crucial to the regulation of various signaling proteins, the functions of which are essential in many pathophysiological systems. Because structural heterogeneity is high in GAG chains and purification is difficult, the use of structurally defined GAG oligosaccharides from natural sources as molecular models in both biophysical and pharmacological assays is limited. To overcome this obstacle, GAG-like oligosaccharides of well-defined structures are currently being synthesized by chemical and/or enzymatic means in many research groups around the world. These synthetic GAG oligosaccharides serve as useful molecular tools in studies of GAG-protein interactions. In this review, besides discussing the commonest routes used for the synthesis of GAG oligosaccharides, we also survey some libraries of these synthetic models currently available for research and discuss their activities in interaction studies with functional proteins, especially through the microarray approach.

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.

Exploiting Uniformly 13C-Labeled Carbohydrates for Probing Carbohydrate-Protein Interactions by NMR Spectroscopy

NMR of a uniformly ¹³C-labeled carbohydrate was used to elucidate the atomic details of a sugar-protein complex. The structure of the ¹³C-labeled Manα(1-2)Manα(1-2)ManαOMe trisaccharide ligand, when bound to cyanovirin-N (CV-N), was characterized and revealed that in the complex the glycosidic linkage torsion angles between the two reducing-end mannoses are different from the free trisaccharide. Distances within the carbohydrate were employed for conformational analysis, and NOE-based distance mapping between sugar and protein revealed that Manα(1-2)Manα(1-2)ManαOMe is bound more intimately with its two reducing-end mannoses into the domain A binding site of CV-N than with the nonreducing end unit. Taking advantage of the ¹³C spectral dispersion of ¹³C-labeled carbohydrates in isotope-filtered experiments is a versatile means for a simultaneous mapping of the binding interactions on both, the carbohydrate and the protein.

Synthesis of 3-O-Sulfated Oligosaccharides to Understand the Relationship between Structures and Functions of Heparan Sulfate

The sulfation at the 3-OH position of glucosamine is an important modification in forming structural domains for heparan sulfate to enable its biological functions. Seven 3-O-sulfotransferase isoforms in the human genome are involved in the biosynthesis of 3-O-sulfated heparan sulfate. As a rare modification present in heparan sulfate, the availability of 3-O-sulfated oligosaccharides is very limited. Here, we report the use of a chemoenzymatic synthetic approach to synthesize six 3-O-sulfated oligosaccharides, including three hexasaccharides and three octasaccharides. The synthesis was achieved by rearranging the enzymatic modification sequence to accommodate the substrate specificity of 3-O-sulfotransferase 3. We studied the impact of 3-O-sulfation on the conformation of the pyranose ring of 2-O-sulfated iduronic acid using NMR, and on the correlation between ring conformation and anticoagulant activity. We identified a novel octasaccharide that interacts with antithrombin and displays anti factor Xa activity. Interestingly, the octasaccharide displays a faster clearance rate than fondaparinux, an FDA-approved pentasaccharide drug, in a rat model, making this octasaccharide a potential short-acting anticoagulant drug candidate that could reduce bleeding risk. Having access to a set of critically important 3-O-sulfated oligosaccharides offers the potential to develop new heparan sulfate-based therapeutics.

Structure of the Complex between a Heparan Sulfate Octasaccharide and Mycobacterial Heparin-Binding Hemagglutinin

Heparin-binding hemagglutinin (HBHA) is a 199 amino acid virulence factor at the envelope of Mycobacterium tuberculosis that contributes to latent tuberculosis. The binding of HBHA to respiratory epithelial cells, which leads to extrapulmonary dissemination of the pathogen, is mediated by cell-surface heparan sulfate (HS). We report the structural characterization of the HBHA/HS complex by NMR spectroscopy. To develop a model for the molecular recognition, the first chemically synthesized uniformly ¹³ C- and ¹⁵ N-labeled HS octasaccharide and a uniformly ¹³ C- and ¹⁵ N-labeled form of HBHA were prepared. Residues 180-195 at the C-terminal region of HBHA show large chemical shift perturbation upon association with the octasaccharide. Molecular dynamics simulations conforming to the multidimensional NMR data revealed key electrostatic and even hydrophobic interactions between the binding partners that may aid in the development of agents targeting the binding event.

"Rules of Engagement" of Protein-Glycoconjugate Interactions: A Molecular View Achievable by using NMR Spectroscopy and Molecular Modeling

Understanding the dynamics of protein-ligand interactions, which lie at the heart of host-pathogen recognition, represents a crucial step to clarify the molecular determinants implicated in binding events, as well as to optimize the design of new molecules with therapeutic aims. Over the last decade, advances in complementary biophysical and spectroscopic methods permitted us to deeply dissect the fine structural details of biologically relevant molecular recognition processes with high resolution. This Review focuses on the development and use of modern nuclear magnetic resonance (NMR) techniques to dissect binding events. These spectroscopic methods, complementing X-ray crystallography and molecular modeling methodologies, will be taken into account as indispensable tools to provide a complete picture of protein-glycoconjugate binding mechanisms related to biomedicine applications against infectious diseases.

Current structural biology of the heparin interactome

Heparin is the best known therapeutically active carbohydrate. It can bind and regulate multiple functional proteins such as coagulation cofactors, chemokines, and growth factors. This versatility has led to the recently developed concept of the heparin interactome--a group of proteins that, as the name implies, interact with heparin. The heparin interactome is structurally and functionally diverse. Though natural ligands of this class of proteins may be any of the glycosaminoglycans however, their structural biology is generally studied using heparin as a model compound. NMR spectroscopy contributes significantly to structural investigations of the resultant complexes in solution. This review aims therefore at discussing the current status in structural biology of the molecular complexes formed between heparin and its protein partners through the current concept of the heparin interactome.

Carbohydrate–Protein Interactions: A 3D View by NMR

NMR spectroscopy is a key tool for carbohydrate research. In studies with complex oligosaccharides there are limits to the amount of relevant structural information provided by these observables due to problems of signal overlapping, strong coupling and/or the scarcity of the key NOE information. Thus, there is an increasing need for additional parameters with structural information, such as residual dipolar couplings (RDCs), paramagnetic relaxation enhancements (PREs) or pseudo contact shifts (PCSs). Carbohydrates are rather flexible molecules. Therefore, NMR observables do not always correlate with a single conformer but with an ensemble of low free-energy conformers that can be accessed by thermal fluctuations. Depending on the system under study, different NMR approaches can be followed to characterize protein–carbohydrate interactions: the standard methodologies can usually be classified as “ligand-based” or “receptor-based”. The selection of the proper methodology is usually determined by the size of the receptor, the dissociation constant of the complex (KD), the availability of the labelled protein (15N, 13C) and the access to soluble receptors at enough concentration for NMR measurements.

NMR and molecular recognition. The application of ligand-based NMR methods to monitor molecular interactions

NMR allows the monitoring of molecular recognition processes in solution. Nowadays, a plethora of NMR methods are available to deduce the key features of the interaction from both the ligand or the receptor points of view.

Protein-carbohydrate interactions studied by NMR: from molecular recognition to drug design

Diseases that result from infection are, in general, a consequence of specific interactions between a pathogenic organism and the cells. The study of host-pathogen interactions has provided insights for the design of drugs with therapeutic properties. One area that has proved to be promising for such studies is the constituted by carbohydrates which participate in biological processes of paramount importance. On the one hand, carbohydrates have shown to be information carriers with similar, if not higher, importance than traditionally considered carriers as amino acids and nucleic acids. On the other hand, the knowledge on molecular recognition of sugars by lectins and other carbohydrate-binding proteins has been employed for the development of new biomedical strategies. Biophysical techniques such as X-Ray crystallography and NMR spectroscopy lead currently the investigation on this field. In this review, a description of traditional and novel NMR methodologies employed in the study of sugar-protein interactions is briefly presented in combination with a palette of NMR-based studies related to biologically and/or pharmaceutically relevant applications.

New targets for glycosaminoglycans and glycosaminoglycans as novel targets

Biological functions of a variety of proteins are mediated via their interaction with glycosaminoglycans (GAGs). The structural diversity within the wide GAG landscape provides individual interaction sites for a multitude of proteins involved in several pathophysiological processes. This 'GAG angle' of such proteins as well as their specific GAG ligands give rise to novel therapeutic concepts for drug development. Current glycomic technologies to elucidate the glycan structure-function relationships, methods to investigate the selectivity and specificity of glycan-protein interactions and existing therapeutic approaches to interfere with GAG-protein interactions are discussed.

Carbohydrate-aromatic interactions

The recognition of saccharides by proteins has far reaching implications in biology, technology, and drug design. Within the past two decades, researchers have directed considerable effort toward a detailed understanding of these processes. Early crystallographic studies revealed, not surprisingly, that hydrogen-bonding interactions are usually involved in carbohydrate recognition. But less expectedly, researchers observed that despite the highly hydrophilic character of most sugars, aromatic rings of the receptor often play an important role in carbohydrate recognition. With further research, scientists now accept that noncovalent interactions mediated by aromatic rings are pivotal to sugar binding. For example, aromatic residues often stack against the faces of sugar pyranose rings in complexes between proteins and carbohydrates. Such contacts typically involve two or three CH groups of the pyranoses and the π electron density of the aromatic ring (called CH/π bonds), and these interactions can exhibit a variety of geometries, with either parallel or nonparallel arrangements of the aromatic and sugar units. In this Account, we provide an overview of the structural and thermodynamic features of protein-carbohydrate interactions, theoretical and experimental efforts to understand stacking in these complexes, and the implications of this understanding for chemical biology. The interaction energy between different aromatic rings and simple monosaccharides based on quantum mechanical calculations in the gas phase ranges from 3 to 6 kcal/mol range. Experimental values measured in water are somewhat smaller, approximately 1.5 kcal/mol for each interaction between a monosaccharide and an aromatic ring. This difference illustrates the dependence of these intermolecular interactions on their context and shows that this stacking can be modulated by entropic and solvent effects. Despite their relatively modest influence on the stability of carbohydrate/protein complexes, the aromatic platforms play a major role in determining the specificity of the molecular recognition process. The recognition of carbohydrate/aromatic interactions has prompted further analysis of the properties that influence them. Using a variety of experimental and theoretical methods, researchers have worked to quantify carbohydrate/aromatic stacking and identify the features that stabilize these complexes. Researchers have used site-directed mutagenesis, organic synthesis, or both to incorporate modifications in the receptor or ligand and then quantitatively analyzed the structural and thermodynamic features of these interactions. Researchers have also synthesized and characterized artificial receptors and simple model systems, employing a reductionistic chemistry-based strategy. Finally, using quantum mechanics calculations, researchers have examined the magnitude of each property's contribution to the interaction energy.

A perspective on the primary and three-dimensional structures of carbohydrates

Carbohydrates, in more biologically oriented areas referred to as glycans, constitute one of the four groups of biomolecules. The glycans, often present as glycoproteins or glycolipids, form highly complex structures. In mammals ten monosaccharides are utilized in building glycoconjugates in the form of oligo- (up to about a dozen monomers) and polysaccharides. Subsequent modifications and additions create a large number of different compounds. In bacteria, more than a hundred monosaccharides have been reported to be constituents of lipopolysaccharides, capsular polysaccharides, and exopolysaccharides. Thus, the number of polysaccharide structures possible to create is huge. NMR spectroscopy plays an essential part in elucidating the primary structure, that is, monosaccharide identity and ring size, anomeric configuration, linkage position, and sequence, of the sugar residues. The structural studies may also employ computational approaches for NMR chemical shift predictions (CASPER program). Once the components and sequence of sugar residues have been unraveled, the three-dimensional arrangement of the sugar residues relative to each other (conformation), their flexibility (transitions between and populations of conformational states), together with the dynamics (timescales) should be addressed. To shed light on these aspects we have utilized a combination of experimental liquid state NMR techniques together with molecular dynamics simulations. For the latter a molecular mechanics force field such as our CHARMM-based PARM22/SU01 has been used. The experimental NMR parameters acquired are typically (1)H,(1)H cross-relaxation rates (related to NOEs), (3)JCH and (3)JCCtrans-glycosidic coupling constants and (1)H,(13)C- and (1)H,(1)H-residual dipolar couplings. At a glycosidic linkage two torsion angles ϕ and ψ are defined and for 6-substituted residues also the ω torsion angle is required. Major conformers can be identified for which highly populated states are present. Thus, in many cases a well-defined albeit not rigid structure can be identified. However, on longer timescales, oligosaccharides must be considered as highly flexible molecules since also anti-conformations have been shown to exist with H-C-O-C torsion angles of ∼180°, compared to syn-conformations in which the protons at the carbon atoms forming the glycosidic linkage are in close proximity. The accessible conformational space governs possible interactions with proteins and both minor changes and significant alterations occur for the oligosaccharides in these interaction processes. Transferred NOE NMR experiments give information on the conformation of the glycan ligand when bound to the proteins whereas saturation transfer difference NMR experiments report on the carbohydrate part in contact with the protein. It is anticipated that the subtle differences in conformational preferences for glycan structures facilitate a means to regulate biochemical processes in different environments. Further developments in the analysis of glycan structure and in particular its role in interactions with other molecules, will lead to clarifications of the importance of structure in biochemical regulation processes essential to health and disease.

Insights into the glycosaminoglycan-mediated cytotoxic mechanism of eosinophil cationic protein revealed by NMR

Protein-glycosaminoglycan interactions are essential in many biological processes and human diseases, yet how their recognition occurs is poorly understood. Eosinophil cationic protein (ECP) is a cytotoxic ribonuclease that interacts with glycosaminoglycans at the cell surface; this promotes the destabilization of the cellular membrane and triggers ECP's toxic activity. To understand this membrane destabilization event and the differences in the toxicity of ECP and its homologues, the high resolution solution structure of the complex between full length folded ECP and a heparin-derived trisaccharide (O-iPr-α-D-GlcNS6S-α(1-4)-L-IdoA2S-α(1-4)-D-GlcNS6S) has been solved by NMR methods and molecular dynamics simulations. The bound protein retains the tertiary structure of the free protein. The (2)S(0) conformation of the IdoA ring is preferably recognized by the protein. We have identified the precise location of the heparin binding site, dissected the specific interactions responsible for molecular recognition, and defined the structural requirements for this interaction. The structure reveals the contribution of Arg7, Gln14, and His15 in helix α1, Gln40 in strand β1, His64 in loop 4, and His128 in strand β6 in the recognition event and corroborates the previously reported participation of residues Arg34-Asn39. The participation of the catalytic triad (His15, Lys38, His128) in recognizing the heparin mimetic reveals, at atomic resolution, the mechanism of heparin's inhibition of ECP's ribonucleolytic activity. We have integrated all the available data to propose a molecular model for the membrane interaction process. The solved NMR complex provides the structural model necessary to design inhibitors to block ECP's toxicity implicated in eosinophil pathologies.

Toward an accurate conformational modeling of iduronic acid

Iduronic acid (IdoA), unlike most other monosaccharides, can adopt different ring conformations, depending on the context of the molecular structure. Accurate modeling of this building block is essential for understanding the role of glycosaminoglycans and other glycoconjugates. Here, we use metadynamics to predict equilibria of (1)C(4), (4)C(1) and (2)S(O) conformations of α-L-IdoA-OMe and α-L-IdoA2S-OMe. Different schemes of scaling of atoms separated by three bonds (1-4 interaction) were tested. It was found that scaling (reduction) of 1-4 electrostatic interactions significantly changes conformational preferences toward the (4)C(1) conformation. More interestingly, scaling of 1-4 van der Waals interaction favors skew-boat conformations. This shows that a minor modification of noncovalent 1-4 interactions parameters can provide a good agreement between populations of conformers of iduronic acid in water from simulations and experiments.

NMR spectroscopic analysis reveals extensive binding interactions of complex xyloglucan oligosaccharides with the Cellvibrio japonicus glycoside hydrolase family 31 α-xylosidase

The study of the interaction of glycoside hydrolases with their substrates is fundamental to diverse applications in medicine, food and feed production, and biomass-resource utilization. Recent molecular modeling of the α-xylosidase CjXyl31A from the soil saprophyte Cellvibrio japonicus, together with protein crystallography and enzyme-kinetic analysis, has suggested that an appended PA14 protein domain, unique among glycoside hydrolase family 31 members, may confer specificity for large oligosaccharide fragments of the ubiquitous plant polysaccharide xyloglucan (J. Larsbrink, A. Izumi, F.M. Ibatullin, A. Nakhai, H.J. Gilbert, G.J. Davies, H. Brumer, Biochem. J. 2011, 436, 567-580). In the present study, a combination of NMR spectroscopic techniques, including saturation transfer difference (STD) and transfer NOE (TR-NOE) spectroscopy, was used to reveal extensive interactions between CjXyl31A active-site variants and xyloglucan hexa- and heptasaccharides. The data specifically indicate that the enzyme recognizes the entire cello-tetraosyl backbone of the substrate and product in positive enzyme subsites and makes further significant interactions with internal pendant α-(1→6)-linked xylosyl units. As such, the present analysis provides an important rationalization of previous kinetic data on CjXyl31A and unique insight into the role of the PA14 domain, which was not otherwise obtainable by protein crystallography.