Can the stereochemical outcome of glycosylation reactions be controlled by the conformational preferences of the glycosyl donor?

Previous static and dynamical density functional theory studies of the 2,6-di-O-acetyl-3,4-O-isopropylidene-D-galactopyranosyl cations and their methanol adducts has led to an hypothesis that these cations exist in two families of conformers characterized as (2)S(O) and B(2,5), respectively. These families differ by ring inversion, each with its own reactivity. New calculations on the 2,6-di-O-acetyl-3,4-di-O-methyl-D-galactopyranosyl cation confirmed these trends. Removing the isopropylidene group allows more flexibility, but two families of conformers can be discerned with the monocyclic oxocarbenium ions in the E(3) conformation and the bicyclic dioxolenium ions in the (4)H(5) conformation. Attack on the beta-face of these monocyclic cations is favored by hydrogen bonding and the anomeric effect. The experimentally observed high beta-stereoselectivity of mannopyranosyl donors and high alpha-stereoselectivity of glucopyranosyl donors with the 4,6-O-benzylidene protecting groups can be rationalized assuming that the trans-fused 1,3-dioxane ring allows population of only one family of conformers. The combination of hydrogen bonding and conformational changes of the pyranose ring in response to the C-5[bond]O-5[bond]C-1[bond]C-2 torsion angle changes are identified as key factors in stereoselectivity. Based on these observations a strategy to design face discriminated glycosyl donors that exist predominantly in only one family of conformers is proposed.     

Metadynamics modelling of the solvent effect on primary hydroxyl rotamer equilibria in hexopyranosides

Accurate modelling of rotamer equilibria for the primary hydroxyl groups of monosaccharides continues to be a great challenge of computational glycochemistry. The metadynamics technique was applied to study the conformational free energy surfaces of methyl α-d-glucopyranoside and methyl α-d-galactopyranoside, employing the glycam06 force field. For both molecules, seven to eight conformational free-energy minima, differing in the ω (O-5–C-5–C-6–O-6) and χ (C-3–C-4–O-4–HO-4) dihedral angles, were identified in vacuum or in a water environment. The calculated rotamer equilibrium of the primary hydroxyl group is significantly different in vacuum than in water. The major effect of a water environment is the destabilisation of a hydrogen bond between O-4–HO-4 and O-6–HO-6 groups. It was possible to calculate the free-energy differences of individual rotamers with an accuracy of better than 2 kJ/mol. The calculated gg, gt and tg rotamer populations in water are in close agreement with experimental measurements, and therefore support the theoretical background of metadynamics.     

Conformational features of crystal-surface cellulose from higher plants

Native cellulose in higher plants forms crystalline fibrils a few nm across, with a substantial fraction of their glucan chains at the surface. The accepted crystal structures feature a flat-ribbon 21 helical chain conformation with every glucose residue locked to the next by hydrogen bonds from O-3' to O-5 and from O-2 to O-6'. Using solid-state NMR spectroscopy we show that the surface chains have a different C-6 conformation so that O-6 is not in the correct position for the hydrogen bond from O-2. We also present evidence consistent with a model in which alternate glucosyl residues are transiently or permanently twisted away from the flat-ribbon conformation of the chain, weakening the O-3' - 0-5 hydrogen bond. Previous molecular modelling and the modelling studies reported here indicate that this 'translational' chain conformation is energetically feasible and does not preclude binding of the surface chains to the interior chains, because the surface chains share the axial repeat distance of the 21 helix. Reduced intramolecular hydrogen bonding allows the surface chains to form more hydrogen bonds to external molecules in textiles, wood, paper and the living plant.     

Crystal and molecular structure of methyl 4-O-methyl-beta-D-glucopyranosyl-(1-->4)-beta-D-glucopyranoside

The cellulose model compound methyl 4-O-methyl-beta-D-glucopyranosyl-(1-->4)-beta-D-glucopyranoside (6) was synthesised in high overall yield from methyl beta-D-cellobioside. The compound was crystallised from methanol to give colourless prisms, and the crystal structure was determined. The monoclinic space group is P2(1) with Z=2 and unit cell parameters a=6.6060 (13), b=14.074 (3), c=9.3180 (19) A, beta=108.95(3) degrees. The structure was solved by direct methods and refined to R=0.0286 for 2528 reflections. Both glucopyranoses occur in the 4C(1) chair conformation with endocyclic bond angles in the range of standard values. The relative orientation of both units described by the interglycosidic torsional angles [phi (O-5' [bond] C-1' [bond] O-4 [bond] C-4) -89.1 degrees, Phi (C-1' [bond] O-4 [bond] C-4 [bond] C-5) -152.0 degrees] is responsible for the very flat shape of the molecule and is similar to those found in other cellodextrins. Different rotamers at the exocyclic hydroxymethyl group for both units are present. The hydroxymethyl group of the terminal glucose moiety displays a gauche-trans orientation, whereas the side chain of the reducing unit occurs in a gauche-gauche conformation. The solid state (13)C NMR spectrum of compound 6 exhibits all 14 carbon resonances. By using different cross polarisation times, the resonances of the two methyl groups and C-6 carbons can easily be distinguished. Distinct differences of the C-1 and C-4 chemical shifts in the solid and liquid states are found.     

The conformational properties of the glycosidic linkage

Stereochemical properties of the glycosidic linkage have been studied by the quantum-chemical PCILO method, using 2-methoxytetrahydropyran as a model. Calculations of the two-dimensional, conformational (Φ, Ψ) maps showed that the rotation around the C-1---O-1 bond is more hindered than that around the O-1---C-6 bond, and that there are differences in the shape of the energy curve for the axial and equatorial forms of 2-methoxytetrahydropyran. The observed population of the five stable conformers at equilibrium (GG:GT:TG¹:TG₂:TT = 70.8:6.0:19.9:2.0:1.3) is consistent with the prediction of the anomeric and exo-anomeric effects. The calculated abundance (76.8%) of the axial form of 2-methoxytetrahydropyran is comparable with experimental results (77–80%) obtained by n.m.r. measurements in non-polar solvents. The energies found for individual conformers made it possible to calculate the magnitude of the anomeric effect (3 kJ/mol) and to determine, for the first time, the values of the exo-anomeric effect for axial (6 kJ/mol) and equatorial 2-methoxytetrahydropyran (7 kJ/mol). The calculated variations of the geometry arising from rotation around the C-1---O-1 bond are consistent with results obtained by statistical analysis of experimental data for - and β-glycosides. The results obtained, indicating that the energy, geometry, and electronic structure of glycosides are largely affected by the conformation of the acetal segment, are discussed from the point of view of conformational analysis of oligo- and poly-saccharides.     

Conformational preferences in glycosylamines. implications for the exo-anomeric effect

The conformational preferences about the C-N bond in N-(4-methoxyphenyl)-2,3,4,6-tetra-O-acetyl-alpha (1) and beta-D-glucopyranosylamine (2), in the solid state and in solution, have been investigated. The crystal structure of the axially substituted alpha anomer (1) indicates a conformational preference about the C-1-N bond in which nN-->sigma*C-O exo-anomeric interactions may be expressed, although this conformational preference is not displayed in solution. The solution conformation relieves steric interactions that result from expression of the exo-anomeric effect in the solid-state conformation. The conformational preference in the equatorially substituted beta anomer (2) both in solution and in the solid state is similar and permits expression of nN-->sigma*C-O exo-anomeric interactions. The structural data for 1 and 2 indicate significant differences in O-5-C-1-N-1 bond angles but insignificant differences in each of the O-5-C-1 or C-1-N-1 bond lengths. The J(C-1-H-1 coupling constants in 1 and 2 indicate a greater coupling constant for the alpha anomer that is consistent with a dominant nO-->sigma*C-H orbital interaction in the beta anomer that weakens the C-1-H-1 bond.     

Conformational comparison of mu-selective endomorphin-2 with its C-terminal free acid in DMSO solution, by 1H NMR spectroscopy and molecular modeling calculation.

In order to make clear the structural role of the C-terminal amide group of endomorphin-2 (EM2, H-Tyr-Pro-Phe-Phe-NH2), an endogenous mu-receptor ligand, in the biological function, the solution conformations of endomorphin-2 and its C-terminal free acid (EM2OH, H-Tyr-Pro-Phe-Phe-OH), studied using two-dimensional 1H NMR measurements and molecular modeling calculations, were compared. Both peptides were in equilibrium between the cis and trans isomers around the Tyr-Pro omega bond in a population ratio of approximately/= 1:2. The lack of significant temperature and concentration dependence of NH protons suggested that the NMR spectra reflected the conformational features of the respective molecules themselves. Fifty possible 3D structures for the each isomer were generated by the dynamical simulated annealing method under the proton-proton distance constraints derived from the ROE cross-peaks. These energy-minimized conformers, which were all in the phi torsion angles estimated from J(NHCalphaH) coupling constants within +/- 30 degrees, were then classified in groups one or two according to the folding backbone structures. All trans and cis EM2 conformers adopt an open conformation in which their extended backbone structures are twisted at the Pro2-Phe3 moiety. In contrast, the trans and cis conformers of EM2OH show conformational variation between the 'bow'-shaped extended and folded backbone structures, although the cis conformers of its zwitterionic form are refined into the folded structure of the close disposition of C- and N-terminal groups. These results indicate clearly that the substitution of carboxyl group for C-terminal amide group makes the peptide flexible. The conformational requirement for mu-receptor activation has been discussed based on the active form proposed for endomorphin-1 and by comparing conformational features of EM2 and EM2OH.     

A comparative study of the O-3 reactivity of isomeric N-dimethylmaleoyl-protected D-glucosamine and D-allosamine acceptors

Four isomeric N-dimethylmaleoyl 4,6-O-benzylidene-protected d-hexosamine acceptors (2, 3, 4, and 5) with all possible configurations at C-1 and C-3 (e.g., derived from d-glucosamine and D-allosamine) were prepared, and the assessment of their O-3 relative reactivity through competition experiments using the known per-O-acetylated D-galactopyranosyl trichloroacetimidate donor (15) was then carried out. The reactivities are in the order 4?2>5>3. The analysis of the NMR spectra of 2-5 at different temperature and modeling experiments carried out on analogs of 2-5 (DFT) and on the acceptors themselves (MM) are coincident, and have helped to establish the stability of the different hydrogen bonds, and of the conformers which carry them. The whole results suggest that the electronic effects (hydrogen bonds) are required to explain the observed trend, in spite of the axial conformation of the most reactive hydroxyl group. The steric effects appear only when hydrogen bonds are weak.     

8-Chloroguanosine: solid-state and solution conformations and their biological implications

The three-dimensional structure of 8-chloroguanosine dihydrate was determined by X-ray crystallography. The crystals belong to the orthorhombic space group P2(1)2(1)2(1), and the cell dimensions are a = 4.871 (1) A, b = 12.040 (1) A, and c = 24.506 (1) A. The structure was determined by direct methods, and least-squares refinement, which included all hydrogen atoms, converged at R = 0.031 for 1599 observed reflections. The conformation about the glycosidic bond is syn with chi CN = -131.1 degrees. The ribose ring has a C(2')-endo/C-(1')-exo (2T1) pucker, and the gauche+ conformation of the -CH2OH side chain is stabilized by an intramolecular O-(5')-H...N(3) hydrogen bond. Conformational analysis by means of 1H NMR spectroscopy showed that, in dimethyl sulfoxide, the sugar ring exhibits a marked preference for the C(2')-endo conformation (approximately 70%) and a conformation about the glycosidic bond predominantly syn (approximately 90%), hence similar to that in the solid state. However, the conformation of the exocyclic 5'-CH2OH group exhibits only a moderate preference for the gauche+ rotamer (approximately 40%), presumably due to the inability to form the intramolecular hydrogen bond to N(3) in a polar medium. The conformational features are examined in relation to the behavior of 8-substituted purine nucleosides in several enzymatic systems, with due account taken of the steric bulk and electronegativities of the 8-substituents.     

β‐Hairpin folding and stability: molecular dynamics simulations of designed peptides in aqueous solution

The structural properties of a 10‐residue and a 15‐residue peptide in aqueous solution were investigated by molecular dynamics simulation. The two designed peptides, SYINSDGTWT and SESYINSDGTWTVTE, had been studied previously by NMR at 278 K and the resulting model structures were classified as 3:5 β‐hairpins with a type I + G1 β‐bulge turn. In simulations at 278 K, starting from the NMR model structure, the 3:5 β‐hairpin conformers proved to be stable over the time period evaluated (30 ns). Starting from an extended conformation, simulations of the decapeptide at 278 K, 323 K and 353 K were also performed to study folding. Over the relatively short time scales explored (30 ns at 278 K and 323 K, 56 ns at 353 K), folding to the 3:5 β‐hairpin could only be observed at 353 K. At this temperature, the collapse to β‐hairpin‐like conformations is very fast. The conformational space accessible to the peptide is entirely dominated by loop structures with different degrees of β‐hairpin character. The transitions between different types of ordered loops and β‐hairpins occur through two unstructured loop conformations stabilized by a single side‐chain interaction between Tyr2 and Trp9, which facilitates the changes of the hydrogen‐bond register. In agreement with previous experimental results, β‐hairpin formation is initially driven by the bending propensity of the turn segment. Nevertheless, the fine organization of the turn region appears to be a late event in the folding process. Copyright © 2004 European Peptide Society and John Wiley & Sons, Ltd.     

Determination of the DNA sugar pucker using 13C NMR spectroscopy

Solid-state 13C NMR spectroscopy of a series of crystalline nucleosides and nucleotides allows direct measurement of the effect of the deoxyribose ring conformation on the carbon chemical shift. It is found that 3'-endo conformers have 3' and 5' chemical shifts significantly (5-10 ppm) upfield of comparable 3'-exo and 2'-endo conformers. The latter two conformers may be distinguished by smaller but still significant differences in the carbon chemical shifts at the C-2' and C-4' positions. High-resolution solid-state NMR of three modifications of fibrous calf thymus DNA shows that these trends are maintained in high-molecular-weight DNA and confirms that the major ring pucker in A-DNA is 3'-endo, while both B-DNA and C-DNA are largely 2'-endo. The data show that 13C NMR spectroscopy is a straightforward and useful probe of DNA ring pucker in both solution and the solid state.     

The Free Energy Profile of Tubulin Straight-Bent Conformational Changes,with Implications for Microtubule Assembly and Drug Discovery

αβ-tubulin dimers need to convert between a ‘bent’ conformation observed for free dimers in solution and a ‘straight’ conformation required for incorporation into the microtubule lattice. Here, we investigate the free energy landscape of αβ-tubulin using molecular dynamics simulations, emphasizing implications for models of assembly, and modulation of the conformational landscape by colchicine, a tubulin-binding drug that inhibits microtubule polymerization. Specifically, we performed molecular dynamics, potential-of-mean force simulations to obtain the free energy profile for unpolymerized GDP-bound tubulin as a function of the ∼12° intradimer rotation differentiating the straight and bent conformers. Our results predict that the unassembled GDP-tubulin heterodimer exists in a continuum of conformations ranging between straight and bent, but, in agreement with existing structural data, suggests that an intermediate bent state has a lower free energy (by ∼1 kcal/mol) and thus dominates in solution. In agreement with predictions of the lattice model of microtubule assembly, lateral binding of two αβ-tubulins strongly shifts the conformational equilibrium towards the straight state, which is then ∼1 kcal/mol lower in free energy than the bent state. Finally, calculations of colchicine binding to a single αβ-tubulin dimer strongly shifts the equilibrium toward the bent states, and disfavors the straight state to the extent that it is no longer thermodynamically populated.     

Etude cristallographique des β-d-glucopyranosyl- et β-d-mannopyranosyl-benzènes acétylés,analogues de nucléosides pyrimidiques

The structures of the two title C-glycopyranosylarene nucleosides have been determined by X-ray diffraction. The aim of this work was to relate the conformation around the extracyclic C-1C-7 bond to steric hindrance between the pyranose and benzene rings. The torsion angles observed in the two compounds (O-5C-1C-7C-8: +61,7° for 1, ?13,4° for 2) signify of a C-2 configurational modification. Moreover, the interaction between O-5 and an o-phenyl hydrogen could explain the particular conformation of the aryl substituent in 2.     

Identification of dynamical hinge points of the L1 ligase molecular switch

The L1 ligase is an in vitro selected ribozyme that uses a noncanonically base-paired ligation site to catalyze regioselectively and regiospecifically the 5′ to 3′ phosphodiester bond ligation, a reaction relevant to origin of life hypotheses that invoke an RNA world scenario. The L1 ligase crystal structure revealed two different conformational states that were proposed to represent the active and inactive forms. It remains an open question as to what degree these two conformers persist as stable conformational intermediates in solution, and along what pathway are they able to interconvert. To explore these questions, we have performed a series of molecular dynamics simulations in explicit solvent of the inactive–active conformational switch in L1 ligase. Four simulations were performed departing from both conformers in both the reactant and product states, in addition to a simulation where local unfolding in the active state was induced. From these simulations, along with crystallographic data, a set of four virtual torsion angles that span two evolutionarily conserved and restricted regions were identified as dynamical hinge points in the conformational switch transition. The ligation site visits three distinct states characterized by hydrogen bond patterns that are correlated with the formation of specific contacts that may promote catalysis. The insights gained from these simulations contribute to a more detailed understanding of the coupled catalytic/conformational switch mechanism of L1 ligase that may facilitate the design and engineering of new catalytic riboswitches.     

Backbone modifications in cyclic peptides. Conformational analysis of a cyclic pseudopentapeptide containing a thiomethylene ether amide bond replacement

NMR and X-ray crystallographic studies have shown that cyclic pentapeptides of the general structure cyclo(D-Xxx-Pro-Gly-Pro-Gly) possess beta- and gamma-turn intramolecular hydrogen bonds. As part of our continuing series surveying the compatibility of various amide bond replacements on peptide structure, we have synthesized cyclo(D-Phe-Pro psi[CH2S]Gly-Pro-Gly). The pseudopeptide was prepared by solid phase methods and cleaved from the resin by a new procedure involving phase transfer catalysis using K2CO3 and tetrabutylammonium hydrogen sulfate. Cyclization was carried out with the use of DPPA, HOBt, and DMAP to afford the product in 69% yield. The conformational behavior of the pseudopeptide was analyzed by 1H and 13C (1D and 2D) NMR techniques. The backbone modification replaced the amide bond that is involved in a gamma-turn intramolecular hydrogen bond in the all-amide structure. In CDCl3, the pseudopeptide adopted the same all-trans conformation as its parent, although the remaining beta-turn hydrogen bond was weaker according to delta delta/delta TNH measurements. In DMSO-d6, the all-trans conformer and a second conformer were observed in a ratio of 55:45. These conformers, which slowly interconverted on the NMR time scale, could be separately assigned; peaks due to chemical exchange were readily distinguishable by the ROESY technique as reported earlier by others. 13C and ROESY experiments suggested the minor conformer contained one cis amide bond at the Gly1-Pro2 position. Thus, both the location and type of amide surrogate are important determinants affecting the compatibility of the replacement with a particular conformational feature.     

O-2 Substituted pyranosyl oxacarbenium ions are C-2-O-2 2-fold rotors with a strong syn preference

The substituent at O-2 of glycopyranosides is known to have a pronounced effect on both the formation and the cleavage of glycosides at C-1. This is primarily attributed to stereoelectronic effects on the formation and stability of the related glycopyranosyl oxacarbenium ions. Previous QM studies of 2-O-methyl substituted manno and gluco configured pyranosyl oxacarbenium ions found a preference for the methyl carbon to be syn to the CH-2 methine. This study examines the conformational preference of variously substituted O-2 tetrahydropyranosyl oxacarbenium ions and confirms this syn preference. Neutral analogues are shown to have the expected 3-fold rotation whereas the charged species exhibit 2-fold rotation about C-2-O-2. Natural bond order (NBO) calculations suggest that the dominant stabilizing interaction is a unimodal O-2 lone pair to C-1-O-5 pi-bond hyperconjugative interaction. This syn conformational preference has important implications for mimics of glycopyranosyl oxacarbenium ion transition states. It also suggests a conformational based mechanism that can be exploited to tune the reactivity of glycopyranosyl donors in the glycosylation reaction.     

A comparative study of the O-3 reactivity of isomeric N-dimethylmaleoyl-protected d-glucosamine and d-allosamine acceptors

Four isomeric N-dimethylmaleoyl 4,6-O-benzylidene-protected d-hexosamine acceptors (2, 3, 4, and 5) with all possible configurations at C-1 and C-3 (e.g., derived from d-glucosamine and d-allosamine) were prepared, and the assessment of their O-3 relative reactivity through competition experiments using the known per-O-acetylated d-galactopyranosyl trichloroacetimidate donor (15) was then carried out. The reactivities are in the order 4 ? 2 > 5 > 3. The analysis of the NMR spectra of 2–5 at different temperature and modeling experiments carried out on analogs of 2–5 (DFT) and on the acceptors themselves (MM) are coincident, and have helped to establish the stability of the different hydrogen bonds, and of the conformers which carry them. The whole results suggest that the electronic effects (hydrogen bonds) are required to explain the observed trend, in spite of the axial conformation of the most reactive hydroxyl group. The steric effects appear only when hydrogen bonds are weak.     

Structural analysis of two 2-deoxy analogues of alpha- and beta-KDO and the methyl alpha- and beta-glycosides of KDO, and determination of their metal-ion-binding properties

Ammonium 2,6-anhydro-3-deoxy-D-glycero-D-talo-octonate (1), a potent inhibitor of the enzyme CMP-KDO synthetase, its C-2 epimer 2, and the methyl beta- (3) and alpha-glycoside (4) of KDO were studied by 1H- and 13C-n.m.r. spectroscopy. Compound 1 was also analysed by X-ray crystallography. Each compound adopted a 5C2 chair conformation with the side chain equatorial. The preponderant side-chain conformation of 1 in solution was the same as that in the crystal and was stabilised by an intramolecular hydrogen bond from HO-8 to the carboxylate group. This hydrogen bond appeared to be present also in 3. However, the side-chain conformation of 2 and 4 was different from that in 1 and 3. The metal-ion-binding properties, determined on the basis of the line-broadening effects of Mn2+ on the 13C-n.m.r. signals, showed that the carboxylate group was involved in the binding with O-8 in 1 and 3 and with O-6 and O-8 in 2 and 4.     

NMR and molecular modeling characterization of RGD containing peptides.

The tripeptide sequence arginine-glycine-aspartic acid (RGD) has been shown to be the key recognition segment in numerous cell adhesion proteins. The solution conformation and dynamics in DMSO-d6 of the cyclic pentapeptides, [formula: see text], a potent fibrinogen receptor antagonist, and [formula: see text], a weak fibrinogen receptor antagonist, have been characterized by nuclear magnetic resonance (NMR) spectroscopy and molecular modeling. 1H-1H distance constraints derived from two-dimensional NOE spectroscopy and torsional angle constraints obtained from 3JNH-H alpha coupling constants, combined with computer-assisted modeling using conformational searching algorithms and energy minimization have allowed several low energy conformations of the peptides to be determined. Low temperature studies in combination with molecular dynamics simulations suggest that each peptide does not exist in a single, well-defined conformation, but as an equilibrating mixture of conformers in fast exchange on the NMR timescale. The experimental results can be fit by considering pairs of low energy conformers. Despite this inherent flexibility, distinct conformational preferences were found which may be related to the biological activity of the peptides.     

Conformational analysis and molecular dynamics simulations of maltose

Constrained conformational energy minimizations have been used to calculate an adiabatic (Φ, ψ) potential energy surface for the disaccharide β-maltose. The inclusion of molecular flexibility in the conformational energy analysis of the disaccharide was found to significantly lower the barriers to conformational transitions, as has been observed previously for other systems. Several low energy wells were identified on the adiabatic surface which differ in energy by small amounts and with low absolute barriers separating them, indicating the possibility of a non-negligible equilibrium population distribution in each well. If such a distribution of conformations existed in the physical system, the conformation observed by NMR NOE measurements would thus be a “virtual” conformation. Molecular dynamics simulations of the motions of this molecule in vacuum were also conducted and indicate that the rate of relaxation of the molecule to the adiabatic surface may be slower than the typical timescale of conformational fluctuations. This effect is apparently due to an unphysical persistence of hydrogen bond patterns in vacuum which does not occur in aqueous solution. Trajectories undergoing transitions between wells were calculated and the effects of such conformational transitions upon the ensemble mean structure, such as might be observed in an NMR experiment, were demonstrated.