Saturation transfer difference (STD) 1H-NMR experiments and in silico docking experiments to probe the binding of N-acetylneuraminic acid and derivatives to Vibrio cholerae sialidase

NMR investigations of glycan conformation, dynamics, and interactions

Glycans are ubiquitous in nature, decorating our cells and serving as the initial points of contact with any visiting entities. These glycan interactions are fundamental to host-pathogen recognition and are related to various diseases, including inflammation and cancer. Therefore, understanding the conformations and dynamics of glycans, as well as the key features that regulate their interactions with proteins, is crucial for designing new therapeutics. Due to the intrinsic flexibility of glycans, NMR is an essential tool for unravelling these properties. In this review, we describe the key NMR parameters that can be extracted from the different experiments, and which allow us to deduce the necessary geometry and molecular motion information, with a special emphasis on assessing the internal motions of the glycosidic linkages. We specifically address the NMR peculiarities of various natural glycans, from histo-blood group antigens to glycosaminoglycans, and also consider the special characteristics of their synthetic analogues (glycomimetics). Finally, we discuss the application of NMR protocols to study glycan-related molecular recognition events, both from the carbohydrate and receptor perspectives, including the use of stable isotopes and paramagnetic NMR methods to overcome the inherent degeneracy of glycan chemical shifts.

Human disease glycomics: technology advances enabling protein glycosylation analysis - part 1

INTRODUCTION

Protein glycosylation is recognized as an important post-translational modification, with specific substructures having significant effects on protein folding, conformation, distribution, stability and activity. However, due to the structural complexity of glycans, elucidating glycan structure-function relationships is demanding. The fine detail of glycan structures attached to proteins (including sequence, branching, linkage and anomericity) is still best analysed after the glycans are released from the purified or mixture of glycoproteins (glycomics). The technologies currently available for glycomics are becoming streamlined and standardized and many features of protein glycosylation can now be determined using instruments available in most protein analytical laboratories. Areas covered: This review focuses on the current glycomics technologies being commonly used for the analysis of the microheterogeneity of monosaccharide composition, sequence, branching and linkage of released N- and O-linked glycans that enable the determination of precise glycan structural determinants presented on secreted proteins and on the surface of all cells. Expert commentary: Several emerging advances in these technologies enabling glycomics analysis are discussed. The technological and bioinformatics requirements to be able to accurately assign these precise glycan features at biological levels in a disease context are assessed.

The use of glycoinformatics in glycochemistry

Glycoinformatics is a small but growing branch of bioinformatics and chemoinformatics. Various resources are now available that can be of use to glycobiologists, but also to chemists who work on the synthesis or analysis of carbohydrates. This article gives an overview of existing glyco-specific databases and tools, with a focus on their application to glycochemistry: Databases can provide information on candidate glycan structures for synthesis, or on glyco-enzymes that can be used to synthesize carbohydrates. Statistical analyses of glycan databases help to plan glycan synthesis experiments. 3D-Structural data of protein-carbohydrate complexes are used in targeted drug design, and tools to support glycan structure analysis aid with quality control. Specific problems of glycoinformatics compared to bioinformatics for genomics or proteomics, especially concerning integration and long-term maintenance of the existing glycan databases, are also discussed.

NMR characterization and molecular modeling of fucoidan showing the importance of oligosaccharide branching in its anticomplementary activity

Fucoidan is a potent inhibitor of the human complement system whose activity is mediated through interactions with certain proteins belonging to the classical pathway, particularly the protein C4. Branched fucoidan oligosaccharides displayed a higher anticomplementary activity as compared to linear structures. Nuclear magnetic resonance (NMR) characterization of the branched oligosaccharides and saturation transfer difference-NMR experiment of the interaction with the protein C4 allowed the identification of the glycan residues in close contact with the target protein. Transferred nuclear Overhauser effect spectroscopy experiment and molecular modeling of fucoidan oligosaccharides indicated that the presence of side chains reduces the flexibility of the oligosaccharide backbone, which thus adopts a conformation which is very close to the one recognized by the protein C4. Together, these results suggest that branching of fucoidan oligosaccharides, determining their conformational state, has a major impact on their anticomplementary activity.

Conformational free energy surface of alpha-N-acetylneuraminic acid: an interplay between hydrogen bonding and solvation

The conformational free energy surface of alpha-N-acetylneuraminic acid (Neu5Ac, sialic acid) in the space of ring-puckering coordinates was calculated using the metadynamics method. Free energy surfaces in vacuum and with an explicit solvent were calculated in GLYCAM 06 force field. In vacuum three structures are almost equivalently populated, namely, the (2)C(5) chair and the B(3,6)/(2)S(6) and (O)S(3) boat/skew-boat conformations. The B(3,6)/(2)S(6) structure is stabilized by an ionic hydrogen bond between the amide N-H bond and the carboxylic group. However, this structure is unfavorable in a water environment in which the experimentally observed (2)C(5) chair conformation is predicted to be more stable than the other structures. These results indicate that environment significantly influences conformation of Neu5Ac and that Neu5Ac-processing enzymes might modify a conformation of their substrates solely by a changing polarity of the environment. The structure of Neu5Ac bound in influenza neuraminidase ((4)S(2)/B(2,5)) belongs to conformations preferred in a water environment. The free energy penalty of this conformational change was calculated (relative to (2)C(5)) as 10.2 +/- 2.0 and 17.3 +/- 2.0 kJ/mol for (4,O)B/(O)S(3) and (4)S(2), respectively. This result indicates that mimicking of the enzyme-bound conformation is likely to be a viable strategy for the design of neuraminidase inhibitors.

Saturation transfer difference NMR spectroscopy as a technique to investigate protein-carbohydrate interactions in solution

Saturation transfer difference (STD) Nuclear Magnetic Resonance (NMR) spectroscopy is a powerful method for studying protein-ligand interactions in solution. The STD NMR method is capable of identifying the binding epitope of a ligand when bound to its receptor protein. Ligand protons that are in close contact with the receptor protein receive a higher degree of saturation, and as a result stronger STD NMR signals can be observed. Protons that are either less or not involved in the binding process reveal no STD NMR signals. Therefore, the STD NMR method is an excellent tool to investigate how a binding ligand interacts with its receptor molecule. The STD NMR experiment is easy to implement and only small amounts of native protein are required. This chapter comprises a detailed experimental protocol to acquire STD NMR spectra and determine the binding epitope of a ligand bound to its target protein. As representative examples the ligands uridyl-triphosphate (UTP) and uridyl-glucose-diphosphate (UDP-glucose) when bound to the Leishmania major UDP-glucose-pyrophosphorylase (UGP) as target protein are examined.

Structural glycobiology: a game of snakes and ladders

Oligo- and polysaccharides are infamous for being extremely flexible molecules, populating a series of well-defined rotational isomeric states under physiological conditions. Characterization of this heterogeneous conformational ensemble has been a major obstacle impeding high-resolution structure determination of carbohydrates and acting as a bottleneck in the effort to understand the relationship between the carbohydrate structure and function. This challenge has compelled the field to develop and apply theoretical and experimental methods that can explore conformational ensembles by both capturing and deconvoluting the structural and dynamic properties of carbohydrates. This review focuses on computational approaches that have been successfully used in combination with experiment to detail the three-dimensional structure of carbohydrates in a solution and in a complex with proteins. In addition, emerging experimental techniques for three-dimensional structural characterization of carbohydrate-protein complexes and future challenges in the field of structural glycobiology are discussed. The review is divided into five sections: (1) The complexity and plasticity of carbohydrates, (2) Predicting carbohydrate-protein interactions, (3) Calculating relative and absolute binding free energies for carbohydrate-protein complexes, (4) Emerging and evolving techniques for experimental characterization of carbohydrate-protein structures, and (5) Current challenges in structural glycoscience.

Inhibition of bacterial cell division protein FtsZ by cinnamaldehyde

Cinnamaldehyde is a natural product from spices that inhibits cell separation in Bacillus cereus. Cell division is regulated by FtsZ, a prokaryotic homolog of tubulin. FtsZ assembles into the Z-ring at the site of cell division. Here, we report the effect of cinnamaldehyde on FtsZ and hence on the cell division apparatus. Cinnamaldehyde decreases the in vitro assembly reaction and bundling of FtsZ. It is found that cinnamaldehyde perturbs the Z-ring morphology in vivo and reduces the frequency of the Z ring per unit cell length of Escherichia coli. In addition, GTP dependent FtsZ polymerization is inhibited by cinnamaldehyde. Cinnamaldehyde inhibits the rate of GTP hydrolysis and binds FtsZ with an affinity constant of 1.0+/-0.2 microM(-1). Isothermal titration calorimetry reveals that binding of cinnamaldehyde to FtsZ is driven by favorable enthalpic interactions. Further, we map the cinnamaldehyde binding region of FtsZ, using the saturation transfer difference-nuclear magnetic resonance and an in silico docking model. Both predict the cinnamaldehyde binding pocket at the C terminal region involving the T7 loop of FtsZ. Our results show that cinnamaldehyde binds FtsZ, perturbs the cytokinetic Z-ring formation and inhibits its assembly dynamics. This suggests that cinnamaldehyde, a small molecule of plant origin, is a potential lead compound that can be developed as an anti-FtsZ agent towards drug design.

STD NMR spectroscopy and molecular modeling investigation of the binding of N-acetylneuraminic acid derivatives to rhesus rotavirus VP8* core

The VP8* subunit of rotavirus spike protein VP4 contains a sialic acid (Sia)-binding domain important for host cell attachment and infection. In this study, the binding epitope of the N-acetylneuraminic acid (Neu5Ac) derivatives has been characterized by saturation transfer difference (STD) nuclear magnetic resonance (NMR) spectroscopy. From this STD NMR data, it is proposed that the VP8* core recognizes an identical binding epitope in both methyl alpha-D-N-acetylneuraminide (Neu5Acalpha2Me) and the disaccharide methyl S-(alpha-D-N-acetylneuraminosyl)-(2-->6)-6-thio-beta-D-galactopyranoside (Neu5Ac-alpha(2,6)-S-Galbeta1Me). In the VP8*-disaccharide complex, the Neu5Ac moiety contributes to the majority of interaction with the protein, whereas the galactose moiety is solvent-exposed. Molecular dynamics calculations of the VP8*-disaccharide complex indicated that the galactose moiety is unable to adopt a conformation that is in close proximity to the protein surface. STD NMR experiments with methyl 9-O-acetyl-alpha-D-N-acetylneuraminide (Neu5,9Ac(2)alpha2Me) in complex with rhesus rotavirus (RRV) VP8* revealed that both the N-acetamide and 9-O-acetate moieties are in close proximity to the Sia-binding domain, with the N-acetamide's methyl group being saturated to a larger extent, indicating a closer association with the protein. RRV VP8* does not appear to significantly recognize the unsaturated Neu5Ac derivative [2-deoxy-2,3-didehydro-D-N-acetylneuraminic acid (Neu5Ac2en)]. Molecular modeling of the protein-Neu5Ac2en complex indicates that key interactions between the protein and the unsaturated Neu5Ac derivative when compared with Neu5Acalpha2Me would not be sustained. Neu5Acalpha2Me, Neu5Ac-alpha(2,6)-S-Galbeta1Me, Neu5,9Ac(2)alpha2Me, and Neu5Ac2en inhibited rotavirus infection of MA104 cells by 61%, 35%, 30%, and 0%, respectively, at 10 mM concentration. NMR spectroscopic, molecular modeling, and infectivity inhibition results are in excellent agreement and provide valuable information for the design of inhibitors of rotavirus infection.

Determination of the Bound Conformation of a Competitive Nanomolar Inhibitor ofMycobacterium tuberculosis Type II Dehydroquinase by NMR Spectroscopy

The synergy between tuberculosis and the AIDS epidemic, along with the surge of multidrug-resistant isolates of M. tuberculosis, has reaffirmed tuberculosis as a primary public health threat. It is therefore necessary to discover new, safe, and more efficient antibiotics against this disease. On the other hand, mapping the dynamic interactions of inhibitors of a target protein can provide information for the development of more potent inhibitors and consequently, more potent potential drugs. In this context, the conformational binding of our previously reported nanomolar inhibitor of M. tuberculosis type II dehydroquinase, the 3-nitrophenyl derivative 1, was studied using saturation transfer difference (STD) and transferred NOESY experiments. These studies have shown that in the bound state, one conformation of those present in solution of the competitive nanomolar inhibitor 3-nitrophenyl derivative 1 is selected. In the bound conformation, the aromatic ring is slightly shifted from coplanarity, with the double bond and the nitro group of 1 oriented towards the double bond side.

A 1H STD NMR spectroscopic investigation of sialylnucleoside mimetics as probes of CMP-Kdn synthetase

CMP-Kdn synthetase catalyses the reaction of sialic acids (Sia) and CTP to the corresponding activated sugar nucleotide CMP-Sia and pyrophosphate PP( i ). Saturation Transfer Difference (STD) NMR spectroscopy has been employed to investigate the sub-structural requirements of the enzyme's binding domain. Sialylnucleoside mimetics, where the sialic acid moiety has been replaced by a carboxyl group and a hydrophobic moiety, have been used in NMR experiments, to probe the tolerance of the CMP-Kdn synthetase to such replacements. From our data it would appear that unlike another sialylnucleotide-recognising protein, the CMP-Neu5Ac transport protein, either a phosphate group or other functional groups on the sialic acid framework may play important roles in recognition by the synthetase.

Probing a CMP-Kdn synthetase by 1H, 31P, and STD NMR spectroscopy

CMP-Kdn synthetase catalyses the reaction of sialic acids (Sia) and cytidine-5'-triphosphate (CTP) to the corresponding activated sugar nucleotide CMP-Sia and pyrophosphate PP(i). STD NMR experiments of a recombinant nucleotide cytidine-5'-monophosphate-3-deoxy-d-glycero-d-galacto-nonulosonic acid synthetase (CMP-Kdn synthetase) were performed to map the binding epitope of the substrate CTP and the product CMP-Neu5Ac. The STD NMR analysis clearly shows that the anomeric proton of the ribose moiety of both investigated compounds is in close proximity to the protein surface and is likely to play a key role in the binding process. The relative rates of the enzyme reaction, derived from (1)H NMR signal integrals, show that Kdn is activated at a rate 2.5 and 3.1 faster than Neu5Ac and Neu5Gc, respectively. Furthermore, proton-decoupled (31)P NMR spectroscopy was successfully used to follow the enzyme reaction and clearly confirmed the appearance of CMP-Sia and the inorganic pyrophosphate by-product.

An investigation of the activity of recombinant rat skeletal muscle cytosolic sialidase

Rat cytosolic sialidase is expressed at elevated levels in skeletal muscle and is believed to play a role in the myogenic differentiation of muscle cells. Here, we observed varying levels of enhancement of sialidase activity in the presence a range of divalent cations. In particular, a significant enhancement of activity was observed in the presence of Ca2+. Conversely, inhibition of the sialidase activity was found when the enzyme was incubated in the presence of Cu2+, EDTA, and a range of carbohydrate-based inhibitors. Finally, an investigation of the enzymatic hydrolysis of a synthetic substrate, 4-methylumbelliferyl N-acetyl-alpha-D-neuraminide, by 1H NMR spectroscopy revealed that the reaction catalysed by rat skeletal muscle cytosolic sialidase proceeds with overall retention of anomeric configuration. This result further supports the notion that all sialidases appear to be retaining enzymes.