Modeling of deoxy- and dideoxyaldohexopyranosyl ring puckering with MM3(92)

Probing deoxysugar conformational preference: A comprehensive computational study investigating the effects of deoxygenation

Deoxysugars are intrinsic components in a number of antibiotics, antimicrobials, and therapeutic agents that often dictate receptor binding, improve efficacy, and provide a diverse toolbox in modifying glycoconjugate function due to an extensive number of unique isomers and inherent conformational flexibility. Hence, this work provides a comprehensive examination of the conformational effects associated with deoxygenation of the pyranose ring. Both the location and degree of deoxygenation were evaluated by interrogating the energetic landscape for a number of mono- and dideoxyhexopyranose derivatives using DFT methods (M05-2X/cc-pVTZ(-f)). Both anomeric forms and in some cases, the alternate chair form, have been investigated in the gas phase. As was documented in a preceding study, variation of the C-6 oxidation state has been shown to affect the anomeric preference of select glucose stereoisomers. Similar results were also observed for several deoxysugar isomers in this work, wherein the alternate anomer was favored upon reduction to the 6-deoxyhexose derivative or oxidation to the hexonic acid. Additionally, comparison of relative Gibbs free energies revealed C-3 deoxygenation imparts greater instability compared to C-2 or C-4 deoxygenation, as indicated by an increase in free energy for 3-deoxysugars. A polarizable continuum solvation model was also applied to empirically validate theoretical results for several deoxysugars, wherein good agreement with both carbon (σ = 1.6 ppm) and proton (σ = 0.20 ppm) NMR shifts was observed for the majority of isomers. Solvated and gas phase anomeric ratios were also calculated and compared favorably to reported literature values, although some discrepancies are noted.

Structure of l-rhamnose isomerase in complex with l-rhamnopyranose demonstrates the sugar-ring opening mechanism and the role of a substrate sub-binding site

l-Rhamnose isomerase (l-RhI) catalyzes the reversible isomerization of l-rhamnose to l-rhamnulose. Previously determined X-ray structures of l-RhI showed a hydride-shift mechanism for the isomerization of substrates in a linear form, but the mechanism for opening of the sugar-ring is still unclear. To elucidate this mechanism, we determined X-ray structures of a mutant l-RhI in complex with l-rhamnopyranose and d-allopyranose. Results suggest that a catalytic water molecule, which acts as an acid/base catalyst in the isomerization reaction, is likely to be involved in pyranose-ring opening, and that a newly found substrate sub-binding site in the vicinity of the catalytic site may recognize different anomers of substrates.

Analysis of the reaction coordinate of alpha-L-fucosidases: a combined structural and quantum mechanical approach

The enzymatic hydrolysis of alpha-L-fucosides is of importance in cancer, bacterial infections, and fucosidosis, a neurodegenerative lysosomal storage disorder. Here we show a series of snapshots along the reaction coordinate of a glycoside hydrolase family GH29 alpha-L-fucosidase unveiling a Michaelis (ES) complex in a (1)C(4) (chair) conformation and a covalent glycosyl-enzyme intermediate in (3)S(1) (skew-boat). First principles metadynamics simulations on isolated alpha-L-fucose strongly support a (1)C(4)<-->(3)H(4)<-->(3)S(1) conformational itinerary for the glycosylation step of the reaction mechanism and indicate a strong "preactivation" of the (1)C(4) complex to nucleophilic attack as reflected by free energy, C1-O1/O5-C1 bond length elongation/reduction, C1-O1 bond orientation, and positive charge development around the anomeric carbon. Analysis of an imino sugar inhibitor is consistent with tight binding of a chair-conformed charged species.

Synthesis of and molecular dynamics simulations on a tetrasaccharide corresponding to the repeating unit of the capsular polysaccharide from Salmonella enteritidis

Syntheses of two oligosaccharides as methyl glycosides related to the repeating unit of S. enteritidis capsular polysaccharide (CPS) are presented. The trisaccharide corresponds to the backbone of the CPS whereas the tetrasaccharide is a model for the repeating unit which has a branched structure. Molecular dynamics simulations investigating their flexibility and dynamics revealed that the oligosaccharides populate several conformational states and indicate that conformational averaging should be used in describing the accessible conformational space.

Differentiation of the fucoidan sulfated l-fucose isomers constituents by CE-ESIMS and molecular modeling

Alpha-L-fucose, the monosaccharide component of fucoidan, is found in the polysaccharide mainly as its sulfated form where sulfate groups are in position 2 and/or 4 and/or 3. The correlation between biological activities and structure of fucoidan requires the determination of the sulfation pattern of the fucose residues. Therefore, it is of importance to discriminate between the isobaric sulfated fucose isomers. For this purpose, the three isomers 2-O-, 3-O-, and 4-O-sulfated fucose have been analyzed using electrospray ion trap mass spectrometry and capillary electrophoresis. The results reported herein show that it is possible to differentiate between these three positional isomers of sulfated fucose based on their fragmentation pattern upon MS/MS experiments. 3-O-Sulfated fucose was characterized by the loss of the hydrogenosulfate anion HSO4- as the main fragmentation product, while the two other isomers 2-O-, and 4-O-sulfated fucose exhibited cross-ring fragmentation yielding to distinctive (0,2)X and (0,2)A daughter ion, respectively. A computational study of the conformation of the sulfated fucose isomers was carried out providing an understanding of the fragmentation pattern with respect to the position of the sulfate group.

Spectral signatures and structural motifs in isolated and hydrated monosaccharides: phenyl α- and β-l-fucopyranoside

The conformation and structure of phenyl-alpha-l-fucopyranoside (alpha-PhFuc), phenyl-beta-L-fucopyranoside (beta-PhFuc) and their singly hydrated complexes (alpha,beta-PhFuc.H(2)O) isolated in a molecular beam, have been investigated by means of resonant two photon ionization (R2PI) spectroscopy and ultraviolet and infrared ion-dip spectroscopy. Conformational and structural assignments have been based on comparisons between their experimental and computed near IR spectra, calculated using density functional theory (DFT) and their relative energies, determined from ab initio (MP2) calculations. The near IR spectra of "free" and hydrated alpha- and beta-PhFuc, and many other mono- and di-saccharides, provide extremely sensitive probes of hydrogen-bonded interactions which can be finely tuned by small (or large) changes in the molecular conformation. They provide characteristic "signatures" which reflect anomeric, or axial vs. equatorial differences, both revealed through comparisons between alpha/beta-PhFuc and alpha/beta-PhXyl; or similarities, revealed through comparisons between fucose (6-deoxy galactose) and galactose; or binding motifs, for example, "insertion" vs. "addition" structures in hydrated complexes. At the monosaccharide level (the first step in the carbohydrate hierarchy), these trends appear to be general. In contrast to the monohydrates of galactose (beta-PhGal) and glucose (beta-PhGlc), the conformations of alpha- and beta-PhFuc are unaffected by the binding of a single water molecule though changes in the R2PI spectra of multiply hydrated alpha-PhFucW(n) however, may reflect a conformational transformation when n> or = 3.

Modeling ring puckering in strained systems: application to 3,6-anhydroglycosides

Different conformations of methyl 3,6-anhydroglycosides with the beta-D-galacto, alpha-D-galacto, and beta-D-gluco configurations were studied by molecular mechanics (using the program mm3) and by quantum mechanical (QM) methods at the HF/- and B3LYP/6-31+G** levels, with and without solvent emulation. Using molecular mechanics, the energies were plotted against the phi, theta puckering coordinates of Cremer and Pople. In such strained systems, only two extreme conformations of the six-membered ring are likely: (1)C(4) and B(1,4), or any one close to either of them. Results show the preponderance of a distorted chair conformation over that of the distorted boat, though the energy difference is lower and the distortions are larger for the compound with the beta-D-galacto configuration. For derivatives of this compound, experimental data in solution indicate both chair and boat forms, depending on the compound and the solvent, whereas for the remaining compounds, experimental data always show the preponderance of the chair conformation. The more accurate DFT calculations lead to the lower energy differences, suggesting that HF and MM3 underestimate the stability of the boat-like conformations. Similar studies on model compounds depict the importance of the anomeric effect in the conformational preferences.

Comparative performance of MM3(92) and two TINKER? MM3 versions for the modeling of carbohydrates

The 1992 version of MM3 was largely used for modeling mono-, di-, and trisaccharides. In later versions of MM3 improvements were made in some parameters that may be important for carbohydrates. This corrected MM3 force field is part of the Tinker package, freely available (as its 4.1 version), and included in the Chem 3D Ultra 8.0 package (as the 3.7 version). The latter version lacks the corrections to the standard bond lengths produced by electronegativity and anomeric effects, whereas the Tinker 4.1 version only lacks the latter correction. The present work compares the performance of the three MM3 versions (and in some cases, DFT and/or HF/ab initio procedures) on several carbohydrate model problems as the chair and rotamer equilibria in 2-hydroxy- and 2-methoxytetrahydropyran, hydrogen bonding in cis-2,3-dihydroxytetrahydropyran, and the potential energy surfaces around the glycosidic bonds of two sulfated disaccharides and two trisaccharides. Tinker MM3 can be used accurately to estimate carbohydrate energies and geometries, and-with the help of some programming-to pursue studies on the potential energy surfaces of di- and trisaccharides. In most cases results obtained using the three MM3 versions are similar, although large energy differences are obtained when comparing a rotameric distribution around a O-C-O-H dihedral, which is almost forced to the exo-anomeric position by the Tinker versions. In other systems smaller energy differences are found, but they can nevertheless lead to a different global minimum when comparing conformers of similar energy. MM3(92) establishes better the differences between the bond lengths in both anomers, as an expected expression of the anomeric correction.

Guanosine diphosphate-4-keto-6-deoxy-d-mannose reductase in the pathway for the synthesis of GDP-6-deoxy-d-talose in Actinobacillus actinomycetemcomitans

The serotype a-specific polysaccharide antigen of Actinobacillus actinomycetemcomitans is an unusual sugar, 6-deoxy-d-talose. Guanosine diphosphate (GDP)-6-deoxy-d-talose is the activated sugar nucleotide form of 6-deoxy-d-talose, which has been identified as a constituent of only a few microbial polysaccharides. In this paper, we identify two genes encoding GDP-6-deoxy-d-talose synthetic enzymes, GDP-alpha-d-mannose 4,6-dehydratase and GDP-4-keto-6-deoxy-d-mannose reductase, in the gene cluster required for the biosynthesis of serotype a-specific polysaccharide antigen from A. actinomycetemcomitans SUNYaB 75. Both gene products were produced and purified from Escherichia coli transformed with plasmids containing these genes. Their enzymatic reactants were analysed by reversed-phase HPLC (RP-HPLC). The sugar nucleotide produced from GDP-alpha-d-mannose by these enzymes was purified by RP-HPLC and identified by electrospray ionization-MS, 1H nuclear magnetic resonance, and GC/MS. The results indicated that GDP-6-deoxy-d-talose is produced from GDP-alpha-d-mannose. This paper is the first report on the GDP-6-deoxy-d-talose biosynthetic pathway and the role of GDP-4-keto-6-deoxy-d-mannose reductase in the synthesis of GDP-6-deoxy-d-talose.