Molecular modeling of the interaction between heparan sulfate and cellular growth factors: bringing pieces together

Heparan sulfate is a polysaccharide belonging to the glycaminoglycan family. It interacts with numerous proteins of the extracellular matrix, in particular cellular growth factors. The number of experimental protein-heparin sulfate complexes obtained by crystallography or nuclear magnetic resonance is limited. Alternatively, computational approaches can be employed. Generally, they restrain the conformation of the glycosidic rings and linkages in order to reduce the complexity of the problem. Modeling the interaction between protein and heparan sulfate is indeed challenging because of the large size of the fragment needed for a strong binding, the flexibility brought by the glycosidic rings and linkages and the high density of negative charges. We propose a two-step method based on molecular docking and molecular dynamics simulation. Molecular docking allows exploring the positioning of a rigid heparin sulfate fragment on the protein surface. Molecular dynamics refine selected docking models by explicitly representing solvent molecules and not restraining the polysaccharide backbone. The interaction of a hexamer of heparin sulfate was studied in interaction with fibroblast growth factor 2 and stromal cell-derived factor 1α. This approach shed light on the plasticity of the growth factors interacting with heparan sulfate. This approach can be extended to the study of other protein/glycosaminoglycan complexes.

Enhancement of plant and bacterial lectin binding affinities by three-dimensional organized cluster glycosides constructed on helical poly(phenylacetylene) backbones

A series of poly(phenylacetylene)s bearing diverse saccharide pendants--N-acetyl-D-glucosamine, D-lactose, and N-acetyl-D-neuraminic acid--were synthesized by rhodium-mediated polymerizations of the corresponding acetyl-protected glycosylated phenylacetylenes followed by deprotection. The circular dichroism spectra of these glycosylated poly(phenylacetylene)s each displayed split-type Cotton effects in the long absorption region of the conjugated polymer backbone (260-500 nm), thus indicating predominantly one-handed helical conformations in their backbones. The binding affinities of these glycosylated poly(phenylacetylene)s, and those of previously reported phenylacetylenes bearing D-galactose, towards plant and bacterial lectins were investigated by hemagglutination inhibition assay and isothermal titration calorimetry (ITC). The stoichiometries of binding vary strongly, depending on the lectin binding sites and the accessibilities of the carbohydrate residues in the helices. The measured affinities also vary, with the maximum value observed for the interaction between poly-PA-α-Gal and lectin I from Pseudomonas aeruginosa, with a K(d) value of 4 μM per monosaccharide representing a 200-fold increase relative to the corresponding monomer.

Biochemical characterization of the histidine triad protein PhtD as a cell surface zinc-binding protein of pneumococcus

Zinc homeostasis is critical for pathogen host colonization. Indeed, during invasion, Streptococcus pneumoniae has to finely regulate zinc transport to cope with a wide range of Zn(2+) concentrations within the various host niches. AdcAII was identified as a pneumococcal Zn(2+)-binding protein; its gene is present in an operon together with the phtD gene. PhtD belongs to the histidine triad protein family, but to date, its function has not been clarified. Using several complementary biochemical methods, we provide evidence that like AdcAII, PhtD is a metal-binding protein specific for zinc. When Zn(2+) binds (K(d) = 131 ± 10 nM), the protein displays substantial thermal stabilization. We also present the first direct evidence of a joint function of AdcAII and PhtD by demonstrating that their expression is corepressed by Zn(2+), that they interact directly in vitro, and that they are colocalized at the bacterial surface. These results suggest the common involvement of the AdcAII-PhtD system in pneumococcal zinc homeostasis.

NMR and molecular modeling reveal key structural features of synthetic nodulation factors

Nod factors are lipochitoligosaccharides originally produced by the soil bacteria Rhizobia that are involved in the symbiotic process with leguminous plants. Some synthetic analogs of the Nod factors present a strong biological activity, and the conformational behavior of these molecules is of interest for structure/function studies. Nod factor analogs containing an insertion of a phenyl group in the acyl chain at the oligosaccharidic non-reducing end were previously synthesized (Grenouillat N, Vauzeilles B, Bono J-J, Samain E, Beau J-M. 2004. Simple synthesis of nodulation-factor analogues exhibiting high affinity towards a specific binding protein. Angew Chem Int Ed Engl. 43:4644). Conformational studies of natural compounds and synthetic analogs have been performed combining molecular dynamics simulations in explicit water and NMR. Data revealed that the glycosidic head group can adopt only restricted conformations, whereas chemical modifications of the lipid chains, highly flexible in a water environment, influence the global shape of the molecules. Collected structural data could be used in the future to rationalize and understand their biological activity and affinity toward a putative receptor.

Synthesis of dodecavalent fullerene-based glycoclusters and evaluation of their binding properties towards a bacterial lectin

Multivalency is playing a major role in biological processes and particularly in lectin-carbohydrate interactions. The design of high-affinity ligands of lectins should provide molecules capable of interfering with these biological processes and potentially inhibit bacterial or viral infections. Azide-alkyne "click" chemistry was applied to the synthesis of dodecavalent fullerene-based glycoclusters. The conjugation could be efficiently performed from alkyne or azide functions on either partners (i.e. hexakis-fullerene adduct or glycoside). PA-IL is a bacterial lectin from the opportunistic pathogen Pseudomonas aeruginosa and is involved in the recognition of glycoconjugates on human tissues. The glycoclusters obtained were evaluated as ligands of PA-IL and for their potential for competing with its binding to glycosylated surfaces. The affinities measured by hemagglutination inhibition assay (HIA), enzyme-linked lectin assay (ELLA), and surface plasmon resonance (SPR) displayed a significant "glycoside cluster effect" with up to a 12,000-fold increase in binding when comparing a monovalent carbohydrate reference probe with a dodecavalent fullerene-based glycocluster, albeit with some differences depending on the analytical technique.

Nanoelectronic detection of lectin-carbohydrate interactions using carbon nanotubes

We have used single-walled carbon nanotube field-effect transistor (NTFET) devices to probe the interactions between carbohydrates and their recognition proteins called lectins. These interactions are involved in a wide range of biological processes, such as cell-cell recognition, cell-matrix interaction as well as viral and bacterial infections. In our experiments, NTFETs were functionalized noncovalently with porphyrin-based glycoconjugates synthesized using "click" azide-alkyne chemistry, and change in electrical conductance was measured upon specific binding of two bacterial lectins that present different carbohydrate preference, namely PA-IL, PA-IIL from Pseudomonas aeruginosa and a plant lectin Concanavalin A. However, no significant change in the device characteristics was observed when the devices were exposed to other lectins with different specificity. Detection of PA-IL binding to galactosylated NTFETs was highly sensitive (2 nM) with a measured dissociation constant (K(d) = 6.8 μM) corresponding to literature data. Fluorescence microscopy, atomic force microscopy, UV-vis-NIR spectroscopy, and several control measurements confirmed the NTFET response to selective interactions between carbohydrates and lectins.

Selectivity among two lectins: probing the effect of topology, multivalency and flexibility of "clicked" multivalent glycoclusters

The design of multivalent glycoconjugates has been developed over the past decades to obtain high-affinity ligands for lectin receptors. While multivalency frequently increases the affinity of a ligand for its lectin through the so-called "glycoside cluster effect", the binding profiles towards different lectins have been much less investigated. We have designed a series of multivalent galactosylated glycoconjugates and studied their binding properties towards two lectins, from plant and bacterial origins, to determine their potential selectivity. The synthesis was achieved through copper(I)-catalysed azide-alkyne cycloaddition (CuAAC) under microwave activation between propargylated multivalent scaffolds and an azido-functionalised carbohydrate derivative. The interactions of two galactose-binding lectins from Pseudomonas aeruginosa (PA-IL) and Erythrina cristagalli (ECA) with the synthesized glycoclusters were studied by hemagglutination inhibition assays (HIA), surface plasmon resonance (SPR) and isothermal titration microcalorimetry (ITC). The results obtained illustrate the influence of the scaffold's geometry on the affinity towards the lectin and also on the relative potency in comparison with a monovalent galactoside reference probe.

Structure-function similarities between a plant receptor-like kinase and the human interleukin-1 receptor-associated kinase-4

Phylogenetic analysis has previously shown that plant receptor-like kinases (RLKs) are monophyletic with respect to the kinase domain and share an evolutionary origin with the animal interleukin-1 receptor-associated kinase/Pelle-soluble kinases. The lysin motif domain-containing receptor-like kinase-3 (LYK3) of the legume Medicago truncatula shows 33% amino acid sequence identity with human IRAK-4 over the kinase domain. Using the structure of this animal kinase as a template, homology modeling revealed that the plant RLK contains structural features particular to this group of kinases, including the tyrosine gatekeeper and the N-terminal extension α-helix B. Functional analysis revealed the importance of these conserved features for kinase activity and suggests that kinase activity is essential for the biological role of LYK3 in the establishment of the root nodule nitrogen-fixing symbiosis with rhizobia bacteria. The kinase domain of LYK3 has dual serine/threonine and tyrosine specificity, and mass spectrometry analysis identified seven serine, eight threonine, and one tyrosine residue as autophosphorylation sites in vitro. Three activation loop serine/threonine residues are required for biological activity, and molecular dynamics simulations suggest that Thr-475 is the prototypical phosphorylated residue that interacts with the conserved arginine in the catalytic loop, whereas Ser-471 and Thr-472 may be secondary sites. A threonine in the juxtamembrane region and two threonines in the C-terminal lobe of the kinase domain are important for biological but not kinase activity. We present evidence that the structure-function similarities that we have identified between LYK3 and IRAK-4 may be more widely applicable to plant RLKs in general.

Current trends in the structure-activity relationships of sialyltransferases

Sialyltransferases (STs) represent an important group of enzymes that transfer N-acetylneuraminic acid (Neu5Ac) from cytidine monophosphate-Neu5Ac to various acceptor substrates. In higher animals, sialylated oligosaccharide structures play crucial roles in many biological processes but also in diseases, notably in microbial infection and cancer. Cell surface sialic acids have also been found in a few microorganisms, mainly pathogenic bacteria, and their presence is often associated with virulence. STs are distributed into five different families in the CAZy database (http://www.cazy.org/). On the basis of crystallographic data available for three ST families and fold recognition analysis for the two other families, STs can be grouped into two structural superfamilies that represent variations of the canonical glycosyltransferase (GT-A and GT-B) folds. These two superfamilies differ in the nature of their active site residues, notably the catalytic base (a histidine or an aspartate residue). The observed structural and functional differences strongly suggest that these two structural superfamilies have evolved independently.

Polyester nanoparticles presenting mannose residues: toward the development of new vaccine delivery systems combining biodegradability and targeting properties

We report the synthesis of fully biodegradable polymeric nanoparticles presenting mannose residues at their surface and their interaction with lectins. A simple and versatile method was used to reach the surface functionalization of poly(D,L-lactic acid) (PLA) nanoparticles by mannose moieties: It consists in using an amphiphilic mannosylated poly(ethylene oxide)-b-poly(E-caprolactone) (PEO-b-PCL) diblock copolymer as a bioresorbable surface modifier in a simple nanoprecipitation-evaporation procedure. The size and zeta potential of the nanoparticles were found to depend on the molar copolymer/PLA ratio, demonstrating the influence of the copolymer on the formation of the nanoparticles. The bioavailability of the mannose residues as specific recognition sites on the nanoparticle surface could be demonstrated by a modified enzyme-linked lectin assay (ELLA) using biotin-labeled lectins which interact specifically with alpha-D-mannopyrannoside derivatives. Besides specific interaction by lectin-mannose complex formation, nonspecific adsorption of the proteins on the nanoparticle surface was observed. These results were fully supported by isothermal titration calorimetry experiments which suggested that the balance between specific and nonspecific interactions can be controlled by the amount of glycosylated polymer used for the preparation of the nanoparticles. Such nanoparticles are expected to be specifically recognized by mannose receptors, which are highly expressed in cells of the immune system. The targeting properties of these carrier systems combined with their potential adjuvant effects due to their size in the range of 200-300 nm make them attractive candidates as vaccine delivery systems.

Role of water molecules in structure and energetics of Pseudomonas aeruginosa lectin I interacting with disaccharides

Calcium-dependent lectin I from Pseudomonas aeruginosa (PA-IL) binds specifically to oligosaccharides presenting an alpha-galactose residue at their nonreducing end, such as the disaccharides alphaGal1-2betaGalOMe, alphaGal1-3betaGalOMe, and alphaGal1-4betaGalOMe. This provides a unique model for studying the effect of the glycosidic linkage of the ligands on structure and thermodynamics of the complexes by means of experimental and theoretical tools. The structural features of PA-IL in complex with the three disaccharides were established by docking and molecular dynamics simulations and compared with those observed in available crystal structures, including PA-IL.alphaGal1-2betaGalOMe complex, which was solved at 2.4 A resolution and reported herein. The role of a structural bridge water molecule in the binding site of PA-IL was also elucidated through molecular dynamics simulations and free energy calculations. This water molecule establishes three very stable hydrogen bonds with O6 of nonreducing galactose, oxygen from Pro-51 main chain, and nitrogen from Gln-53 main chain of the lectin binding site. Binding free energies for PA-IL in complex with the three disaccharides were investigated, and the results were compared with the experimental data determined by titration microcalorimetry. When the bridge water molecule was included in the free energy calculations, the simulations predicted the correct binding affinity trends with the 1-2-linked disaccharide presenting three times stronger affinity ligand than the other two. These results highlight the role of the water molecule in the binding site of PA-IL and indicate that it should be taken into account when designing glycoderivatives active against P. aeruginosa adhesion.

Discoidin I from Dictyostelium discoideum and Interactions with oligosaccharides: specificity, affinity, crystal structures, and comparison with discoidin II

Discoidin I (DiscI) and discoidin II (DiscII) are N-acetylgalactosamine (GalNAc)-binding proteins from Dictyostelium discoideum. They consist of two domains: an N-terminal discoidin domain and a C-terminal H-type lectin domain. They were cloned and expressed in high yield in recombinant form in Escherichia coli. Although both lectins bind galactose (Gal) and GalNAc, glycan array experiments performed on the recombinant proteins displayed strong differences in their specificity for oligosaccharides. DiscI and DiscII bind preferentially to Gal/GalNAcbeta1-3Gal/GalNAc-containing and Gal/GalNAcbeta1-4GlcNAcbeta1-6Gal/GalNAc-containing glycans, respectively. The affinity of the interaction of DiscI with monosaccharides and disaccharides was evaluated using isothermal titration calorimetry experiments. The three-dimensional structures of native DiscI and its complexes with GalNAc, GalNAcbeta1-3Gal, and Galbeta1-3GalNAc were solved by X-ray crystallography. DiscI forms trimers with involvement of calcium at the monomer interface. The N-terminal discoidin domain presents a structural similarity to F-type lectins such as the eel agglutinin, where an amphiphilic binding pocket suggests possible carbohydrate-binding activity. In the C-terminal H-type lectin domain, the GalNAc residue establishes specific hydrogen bonds that explain the observed affinity (K(d)=3x10(-4) M). The different specificities of DiscI and DiscII for oligosaccharides were rationalized from the different structures obtained by either X-ray crystallography or molecular modeling.

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

The conformational features of hyaluronic acid, a key polysaccharide with important biological properties, have been determined through the combined used of nuclear magnetic resonance (NMR) spectroscopy and molecular modeling techniques. A decasaccharide fragment of sodium hyaluronate (HA) was submitted to 3.5 ns of molecular dynamics in explicit water environment form. The same decasaccharide was prepared by hyaluronidase digestion for the experimental study. The approach consisted in the measurements of NMR residual dipolar coupling (RDC) which were used to filter the molecular dynamics data by retaining those structures which were in agreement with the experimental observations. Further analysis of the new conformer ensemble (HA(RDC)) and clustering the molecules with respect to their overall length led to seven representative structures, which were described in terms of their secondary motifs, namely the best fitting helix geometry. As a result, this protocol permitted the assessment that hyaluronic acid can adopt two different arrangements, which can be described by a three- or four-folded left-handed helix, with a higher occurrence of the first one.

Structural studies of langerin and Birbeck granule: a macromolecular organization model

Dendritic cells, a sentinel immunity cell lineage, include different cell subsets that express various C-type lectins. For example, epidermal Langerhans cells express langerin, and some dermal dendritic cells express DC-SIGN. Langerin is a crucial component of Birbeck granules, the Langerhans cell hallmark organelle, and may have a preventive role toward HIV, by its internalization into Birbeck granules. Since langerin carbohydrate recognition domain (CRD) is crucial for HIV interaction and Birbeck granule formation, we produced the CRD of human langerin and solved its structure at 1.5 A resolution. On this basis gp120 high-mannose oligosaccharide binding has been evaluated by molecular modeling. Hydrodynamic studies reveal a very elongated shape of recombinant langerin extracellular domain (ECD). A molecular model of the langerin ECD, integrating the CRD structure, has been generated and validated by comparison with hydrodynamic parameters. In parallel, Langerhans cells were isolated from human skin. From their analysis by electron microscopy and the langerin ECD model, an ultrastructural organization is proposed for Birbeck granules. To delineate the role of the different langerin domains in Birbeck granule formation, we generated truncated and mutated langerin constructs. After transfection into a fibroblastic cell line, we highlighted, in accordance with our model, the role of the CRD in the membrane zipping occurring in BG formation as well as some contribution of the cytoplasmic domain. Finally, we have shown that langerin ECD triggering with a specific mAb promotes global rearrangements of LC morphology. Our results open the way to the definition of a new membrane deformation mechanism.

Molecular basis of arabinobio-hydrolase activity in phytopathogenic fungi: crystal structure and catalytic mechanism of Fusarium graminearum GH93 exo-alpha-L-arabinanase

The phytopathogenic fungus Fusarium graminearum secretes a very diverse pool of glycoside hydrolases (GHs) aimed at degrading plant cell walls. alpha-l-Arabinanases are essential GHs participating in the complete hydrolysis of hemicellulose, a natural resource for various industrial processes, such as bioethanol or pharmaceuticals production. Arb93A, the exo-1,5-alpha-l-arabinanase of F. graminearum encoded by the gene fg03054.1, belongs to the GH93 family, for which no structural data exists. The enzyme is highly active (1065 units/mg) and displays a strict substrate specificity for linear alpha-1,5-l-arabinan. Biochemical assays and NMR experiments demonstrated that the enzyme releases alpha-1,5-l-arabinobiose from the nonreducing end of the polysaccharide. We determined the crystal structure of the native enzyme and its complex with alpha-1,5-l-arabinobiose, a degradation product of alpha-Me-1,5-l-arabinotetraose, at 1.85 and 2.05A resolution, respectively. Arb93A is a monomeric enzyme, which presents the six-bladed beta-propeller fold characteristic of sialidases of clan GHE. The configuration of the bound arabinobiose is consistent with the retaining mechanism proposed for the GH93 family. Catalytic residues were proposed from the structural analysis, and site-directed mutagenesis was used to validate their role. They are significantly different from those observed for GHE sialidases.

Combination of several bioinformatics approaches for the identification of new putative glycosyltransferases in Arabidopsis

Approximately 450 glycosyltransferase (GT) sequences have been already identified in the Arabidopsis genome that organize into 40 sequence-based families, but a vast majority of these gene products remain biochemically uncharacterized open reading frames. Given the complexity of the cell wall carbohydrate network, it can be inferred that some of the biosynthetic genes have not yet been identified by classical bioinformatics approaches. With the objective to identify new plant GT genes, we designed a bioinformatic strategy that is based on the use of several remote homology detection methods that act at the 1D, 2D, and 3D level. Together, these methods led to the identification of more than 150 candidate protein sequences. Among them, 20 are considered as putative glycosyltransferases that should further be investigated since known GT signatures were clearly identified.

Examination of the effect of structural variation on the N-glycosidic torsion (PhiN) among N-(beta-D-glycopyranosyl)acetamido and propionamido derivatives of monosaccharides based on crystallography and quantum chemical calculations

GlcNAcbetaAsn linkage is conserved in the N-glycoproteins of all eukaryotes. l-Glutamine (Gln), which is a one carbon higher homolog of Asn, is never glycosylated. X-ray crystallographic study of several beta-1-N-acetamido- and propionamido derivatives of monosaccharides has earlier shown that the N-glycosidic torsion, Phi(N), is influenced to a larger extent by the structural variation of the sugar part than that of the aglycon moiety. In order to examine the influence of the carbohydrate pendent groups on the conformational preference of the N-glycosidic linkage with respect to Phi(N,) several models and analogs with gluco and manno configuration have been studied in the present work by computational chemistry. The crystal structure of XylbetaNHPr is reported here and its molecular packing compared with related analogs. The conjunction of combining Crystallographic and computational studies allows to demonstrate the strong influence that the group at C2, and environmental factors particularly inter- and intramolecular interactions involving regular hydrogen bonds and the weak C-H...O contacts, have on the energy preference of the Phi(N) torsion angle.

Catalytic key amino acids and UDP-sugar donor specificity of a plant glucuronosyltransferase, UGT94B1: molecular modeling substantiated by site-specific mutagenesis and biochemical analyses

The plant UDP-dependent glucosyltransferase (UGT) BpUGT94B1 catalyzes the synthesis of a glucuronosylated cyanidin-derived flavonoid in red daisy (Bellis perennis). The functional properties of BpUGT94B1 were investigated using protein modeling, site-directed mutagenesis, and analysis of the substrate specificity of isolated wild-type and mutated forms of BpUGT94B1. A single unique arginine residue (R25) positioned outside the conserved plant secondary product glycosyltransferase region was identified as crucial for the activity with UDP-glucuronic acid. The mutants R25S, R25G, and R25K all exhibited only 0.5% to 2.5% of wild-type activity with UDP-glucuronic acid, but showed a 3-fold increase in activity with UDP-glucose. The model of BpUGT94B1 also enabled identification of key residues in the acceptor pocket. The mutations N123A and D152A decreased the activity with cyanidin 3-O-glucoside to less than 15% of wild type. The wild-type enzyme activity toward delphinidin-3-O-glucoside was only 5% to 10% of the activity with cyanidin 3-O-glucoside. Independent point mutations of three residues positioned near the acceptor B ring were introduced to increase the activity toward delphinidin-3-O-glucoside. In all three mutant enzymes, the enzymatic activity toward both acceptors was reduced to less than 15% of wild type. The model of BpUGT94B1 allowed for correct identification of catalytically important residues, within as well as outside the plant secondary product glycosyltransferase motif, determining sugar donor and acceptor specificity.