B3LYP/6-311++G** study of alpha- and beta-D-glucopyranose and 1,5-anhydro-D-glucitol: 4C1 and 1C4 chairs, (3,O)B and B(3,O) boats, and skew-boat conformations

Geometry optimization, at the B3LYP/6-311++G** level of theory, was carried out on 4C1 and 1C4 chairs, (3,O)B and B(3,O) boats, and skew-boat conformations of alpha- and beta-D-glucopyranose. Similar calculations on 1,5-anhydro-D-glucitol allowed examination of the effect of removal of the 1-hydroxy group on the energy preference of the hydroxymethyl rotamers. Stable minimum energy boat conformers of glucose were found, as were stable skew boats, all having energies ranging from approximately 4-15 kcal/mol above the global energy 4C1 chair conformation. The 1C4 chair electronic energies were approximately 5-10 kcal/mol higher than the 4C1 chair, with the 1C4 alpha-anomers being lower in energy than the beta-anomers. Zero-point energy, enthalpy, entropy, and relative Gibbs free energies are reported at the harmonic level of theory. The alpha-anomer 4C1 chair conformations were found to be approximately 1 kcal/mol lower in electronic energy than the beta-anomers. The hydroxymethyl gt conformation was of lowest electronic energy for both the alpha- and beta-anomers. The glucose alpha/beta anomer ratio calculated from the relative free energies is 63/37%. From a numerical Hessian calculation, the tg conformations were found to be approximately 0.4-0.7 kcal/mol higher in relative free energy than the gg or gt conformers. Transition-state barriers to rotation about the C-5-C-6 bond were calculated for each glucose anomer with resulting barriers to rotation of approximately 3.7-5.8 kcal/mol. No energy barrier was found for the path between the alpha-gt and alpha-gg B(3,O) boat forms and the equivalent 4C1 chair conformations. The alpha-tg conformation has an energy minimum in the 1S3 twist form. Other boat and skew-boat forms are described. The beta-anomer boats retained their starting conformations, with the exception of the beta-tg-(3,O)B boat that moved to a skew form upon optimization.     

DFT study of alpha- and beta-D-galactopyranose at the B3LYP/6-311++G** level of theory

Forty-one conformations of alpha- and beta-d-galactopyranose were geometry optimized using the B3LYP density functional and 6-311++G** basis set. Full geometry optimization was performed on different ring geometries and different hydroxymethyl rotamers (gg/gt/tg). Analytically derived Hessians were used to calculate zero point energy, enthalpy, and entropy. The lowest energy and free-energy conformation found is the alpha-gg-(4)C(1)-c chair conformation, which is of lower electronic and free energy than the lowest energy alpha-d-glucopyranose conformer because of favorable hydrogen-bonding interactions. The in vacuo calculations showed considerable ( approximately 2.2kcal/mol) energetic preference for the alpha over the beta anomer for galactopyranose in both the (4)C(1) and (1)C(4) chair conformations. Results are compared to glucopyranose and mannopyranose calculations in vacuo. Boat and skew-boat forms were found that remained stable upon gradient optimization, although many starting conformations moved to other boat forms upon optimization. As with glucopyranose and mannopyranose, the orientation and interaction of the hydroxyl groups make the most significant contributions to the conformation-energy relationship in vacuo.     

DFT study of alpha- and beta-D-mannopyranose at the B3LYP/6-311++G** level

Thirty-five conformations of alpha- and beta-d-mannopyranose, the C-2 substituted epimer of glucopyranose, were geometry optimized using the density functional (B3LYP), and basis set (6-311++G**). Full geometry optimization was performed on the hydroxymethyl rotamers (gg/gt/tg) and an analytical hessian program was used to calculate the harmonic vibrational frequencies, zero point energy, enthalpy, and entropy. The lowest energy conformation investigated is the beta-tg in the (4)C(1) chair conformation. The in vacuo calculations showed little energetic preference for either the alpha or beta anomer for mannopyranose in the (4)C(1) chair conformation. Results are compared to similar glucopyranose calculations in vacuo where the alpha anomer is approximately 1kcal/mol lower in electronic energy than the beta anomer. In the case of the generally higher energy (1)C(4) chair conformations, one low-energy, low-entropy beta-gg-(1)C(4) chair conformation was identified that is within approximately 1.4kcal/mol of the lowest energy (4)C(1) conformation of mannopyranose. Other (1)C(4) chair conformations in our investigation are approximately 2.9-7.9kcal/mol higher in overall energy. Many of the (3,O)B, B(3,O), (1,4)B, and B(1,4) boat forms passed through transitions without barriers to (1)S(3), (5)S(1), (1)S(5) skew forms with energies between approximately 3.6 and 8.9kcal/mol higher in energy than the lowest energy conformation of mannopyranose. Boat forms were found that remained stable upon gradient optimization. As with glucopyranose, the orientation and interaction of the hydroxy groups make a significant contribution to the conformation/energy relationship in vacuo.     

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.     

B3LYP/6-311++G** study of structure and spin-spin coupling constant in methyl 2-O-sulfo-alpha-L-iduronate

Structures of three most stable conformers ((1)C4, (4)C1, (2)S0) of methyl 2-O-sulfo-alpha-L-iduronate monosodium salt have been analyzed by DFT using the B3LYP/6-311++G** method. The optimized geometries confirmed the influence of both 2-O-sulfate and carboxylate groups upon the pyranose ring geometry. The computed energies showed that the chair (1)C4 form is the most stable one. Time-averaged DFT-calculated proton-proton and proton-carbon spin-spin coupling constants agree with the experimental data and indicate that only two chair forms contribute to the conformational equilibrium of methyl 2-O-sulfo-alpha-L-iduronate monosodium salt. The influence of the charged groups upon the magnitudes of spin-spin coupling constants is also discussed.     

Conformation and absolute configuration of 11,12-cyclopropyl retinal analogs – an RHF/DFT/CIS study

Ab initio RHF and DFT/B3LYP calculations at the 6-31G** level have been performed to study possible conformations of the cyclopropyl retinal Schiff base analog 3 of known absolute configuration. In both the free base and the protonated form, the geometries are determined on the diene side by optimum conjugative interaction with the three-membered ring, on the triene side by repulsive interaction with the 9-methyl group. There are three low energy conformations, in which the seven-membered ring is either in a chair or in a twist-chair conformation. To decide between these alternatives, chiroptical parameters were calculated employing the GAUSSIAN/CIS routines and compared with the CD spectrum obtained by Nakanishi et al. Of the energy-minimized geometries only two fit the experimental data. In both, the dihedral angle C12-C13, which is indicative of the relative orientation of the two chromophores, is positive.     

B3LYP/6-311++G* * study of structure and spin-spin coupling constant in heparin disaccharide

Structures of heparin disaccharide have been analyzed by DFT using the B3LYP/6-311++G( * *) method. The optimized geometries of two forms of this disaccharide, differing in the conformation ((1)C(4) and (2)S(0)) of the IdoA2S residue, confirmed considerable influences of the sulfate and the carboxylate groups upon the pyranose ring geometries. The computed energies showed that disaccharide having the (1)C(4) form of the IdoA2S residue is more stable than that with the (2)S(0) form. Interatomic distances, bond and torsion angles showed that interconversion of the IdoA2S residue results in geometry changes in the GlcN,6S residue as well. Three-bond proton-proton and proton-carbon spin-spin coupling constants computed for both forms agree with the experimental data and indicate that only two chair forms contribute to the conformational equilibrium in disaccharide. Influences of the charged groups upon the magnitudes of spin-spin coupling constants are also discussed.     

A theoretical investigation on the proton transfer tautomerization mechanisms of 2-thioxanthine within microsolvent and long range solvent

A relative complete study on the mechanisms of the proton transfer reactions of 2-thioxanthine was carried out with density functional theory. The models were designed with monohydrated and dihydrated microsolvent catalyses either with or without the presence of water solvent considered with the polarized continuum model (PCM). A total number of 114 complexes and 67 transition states were found with the B3LYP/6-311+G** calculations. The energies were refined with both B3LYP/aug-cc-pVTZ and PCM-B3LYP/aug-cc-pVTZ methods. The activation energies were reported with respect to the Gibbs free energies obtained in conjunction with the standard statistical thermodynamics. Possible reaction pathways were confirmed with the intrinsic reaction coordinates. Pathways via C8 atom on the imidazole ring, via the bridged C4 and C5 atoms between pyrimidine and imidazole rings and via N, O and S atom on the pyrimidine ring were examined. The results show that the most feasible pathway is the proton transfers within the long range solvent surrounding via the N, O and S atoms in the pyrimidine ring with di-hydrated catalysis: N(7)H?+?2H₂O?→?IM89?→?IM90?→?P13?+?2H₂O?→?IM91?→?IM92?→?P6?+?2H₂O?→?IM71?→?IM72?→?P7?+?2H₂O?→?IM107?→?IM108?→?P18?+?2H₂O?→?IM111?→?IM112?→?P19?+?2H₂O?→?IM113?→?IM114?→?P17?+?2H₂O?→?IM105?→?IM106?→?N(9)H?+?2H₂O that has the highest energy barrier of 44.0 kJ mol^?1 in the transition of IM89 to IM90 via TS54. The small energy barrier is in good agreement with the experimental observation that 2-TX tautomerizes at room temperature in water. In the aqueous phase, the most stable intermediate is found to be IM21 [N(7)H?+?2H₂O] and the possible co-existing species are the monohydrated IM1, IM9, IM39 and IM46, and the di-hydrated IM5, IM8, IM13, IM16, IM81, IM89, IM90, IM91 and IM106 complexes that have a relative concentration larger than 10^?6 (1 ppm) with respect to IM21.

Study of the hydrogen bond in different orientations of adenine-thymine base pairs: An <Emphasis Type="Italic">ab initio</Emphasis> study

In order to gain deeper insight into structure, charge distribution, and energies of A-T base pairs, we have performed quantum chemical ab initio and density functional calculations at the HF (Hartree-Fock) and B3LYP levels with 3-21G*, 6-31G*, 6-31G**, and 6-31++G** basis sets. The calculated donor-acceptor atom distances in the Watson-Crick A-T base pair are in good agreement with the experimental mean values obtained from an analysis of 21 high resolution DNA structures. In addition, for further correction of interaction energies between adenine and thymine, the basis set superposition error (BSSE) associated with the hydrogen bond energy has been computed via the counterpoise method using the individual bases as fragments. In the Watson-Crick A-T base pair there is a good agreement between theory and experimental results. The distances for (N2...H23-N19), (N8-H13...O24), and (C1...O18) are 2.84, 2.94, and 3.63 A, respectively, at B3LYP/6-31G** level, which is in good agreement with experimental results (2.82, 2.98, and 3.52 A). Interaction energy of the Watson-Crick A-T base pair is -13.90 and -10.24 kcal/mol at B3LYP/6-31G** and HF/6-31G** levels, respectively. The interaction energy of model (9) A-T base pair is larger than others, -18.28 and -17.26 kcal/mol, and for model (2) is the smallest value, -13.53 and -13.03 kcal/mol, at B3LYP/6-31G** and B3LYP/6-31++G** levels, respectively. The computed B3LYP/6-31G** bond enthalpies for Watson-Crick A-T pairs of -14.4 kcal/mol agree well with the experimental results of -12.1 kcal/mol deviating by as little as -2.3 kcal/mol. The BSSE of some cases is large (9.85 kcal/mol) and some is quite small (0.6 kcal/mol).     

Several transition states from C1 to skew conformations of β-d-glucopyranose

The interconversion pathways in the ring distortion of β-d-glucopyranose were investigated using density functional calculations. We examined the energies of several conformers of β-d-glucopyranose and tried to obtain the transition-state conformation and determine the pathway between a ⁴C₁ chair and some distorted ring conformers. The results showed that two E₃/²H₃ conformations and one E₃/⁴H₃ conformation were transition states in such ring puckering. The transition state with the lowest energy conformation is the E₃/²H₃ ring conformation with the side-chain conformation of r-ggG⁺. Intrinsic reaction coordinate calculations indicated that the E₃/²H₃ conformation with the lowest conformational energy is a transition state of the ring interconversion path between the conformations of ⁴C₁ and ²S_O/B_3,O. The energy barrier of this interconversion was 6.13 kcal/mol. As far as we know, this is the first example of finding pathways for an interconversion of glucopyranose ring puckering at the level of a quantum chemical calculation.     

Stereochemical studies on nucleic acid analogues. I. Conformations of alpha-nucleosides and alpha-nucleotides: interconversion of sugar puckers via O4'-exo.

The preferred conformations of ribo and deoxyribo alpha-nucleosides and alpha-nucleotides, the stereoisomers of the naturally occurring beta-isomers, are worked out by minimizing the conformational energy as a function of all the major parameters including the sugar ring conformations along the pseudorotation path. The results of the studies bring out certain distinct conformational features that are at variance with their beta counterparts. The range of glycosyl conformations are found to be not only quite restricted here but favor predominantly the anti conformation. The syn glycosyl conformation for the entire region of P values are found to be energetically less favorable, with the barrier to anti in equilibrium with syn interconversion being higher especially in alpha-ribonucleosides. The energetically preferred range of pseudorotation phase angles (P) is also considerably restricted and P values corresponding to the C1'-exo range of sugars are highly unfavorable for alpha-nucleosides, in sharp contrast to the broad range of sugar ring conformations favored by beta-isomers. While both trans congruent to 180 degrees and skew congruent to 270 degrees conformations around the C3'-O3' (phi') bond are favored for alpha-3'-nucleotides with deoxyribose sugars, ribose sugars seem to favor only the skew values of phi'. Most interestingly and in sharp contrast to beta-stereoisomers, an energy barrier is encountered at P values corresponding to O4'-endo sugars. This suggests that the possible sugar pucker interconversion between C2'-endo/C3'-exo and C3'-endo/C2'-exo in alpha-anomers could take place only through the O4'-exo region. Likewise the possible path of anti in equilibrium with syn interconversion in alpha-nucleosides is not via high anti, in sharp contrast to beta-nucleosides. These observations should be borne in mind while assigning the sugar ring conformers in alpha-nucleosides and those containing them from nmr investigations. Comparison of the results with beta-anomers seem to suggest on the whole a lack of conformational variability or the restricted nature of alpha-stereoisomers. This could be one of the reasons for its nonselection in the naturally occurring nucleic acids.     

Electronic structure simulations of 2,6-dimethyl-2,5-heptadien-4-one by FTIR,FT-Raman,NMR, UV–vis,NBO and density functional theory

Density functional theory (DFT) (B3LYP and B3PW91) calculations have been carried out for 2,6-dimethyl-2,5-heptadien-4-one (DMHD4O) using 6–311++ G** basis set. Complete vibrational assignment and analysis of the fundamental modes of the compound were carried out from the FTIR and FT-Raman spectral data. The theoretical electronic absorption has been calculated by using time-dependent DFT (TD-DFT) methods and compared with the experimental spectra. The theoretically computed Frontier energy gaps and TD-DFT calculations are in good agreement with the experimental UV–vis spectral absorption. The chemical hardness measured from the Frontier molecular orbital energies of DMHD4O is 0.0693 eV. Electronic stability of the compound arising from hyperconjugative interactions and charge delocalisation were also investigated based on the natural bond orbital (NBO) analysis. Effective stabilisation energy E ⁽²⁾ associated with the interactions of the π and the lone pair of electrons was determined by the NBO analysis. ¹³C and ¹H NMR chemical shifts of the compound have been calculated by means of Gauge-Invariant Atomic Orbital using B3LYP/6–311++ G** method. The partial ionic character of the carbonyl group due to resonance render a partially positive charge to the carbonyl carbon, and thus C4 chemical shift lie in the very downfield 191.6 ppm. Comparison between the experimental and the theoretical results indicates that B3LYP method is able to provide satisfactory results for predicting vibrational, electronic and NMR properties.     

An epoxidation mechanism of carbamazepine by CYP3A4

Human CYP3A4 catalyzes the 10,11-epoxidation of carbamazepine (CBZ). However, the epoxide is less stable in terms of potential energy than hydroxides of the six-membered aromatic ring. To clarify the reason why CYP3A4 produces such an energetically unfavorable compound, the mechanism of epoxidation of CBZ by CYP3A4 was investigated by theoretical calculations. The reaction consisted of two elementary processes in which two C–O bonds were generated stepwise. The rate-determining step was the first one and the activation energy was 21.3 kcal/mol at the DFT (B3LYP/6-31G^**) level. The activation energy level of the first step of the 10,11-epoxidation was lower than that of the hydroxylation of the aromatic ring. For this reason, 10,11-epoxidation is more probable than hydroxylation of the aromatic ring, and only 10,11-epoxide is formed.     

The influence of the peptide bond on the conformation of amino acids: a theoretical and FT-IR matrix-isolation study of N-acetylproline

A combined experimental matrix-isolation FT-IR and theoretical study has been performed to investigate the conformational behavior of N-acetylproline. The conformational landscape of N-acetylproline was explored using successively higher computational methods, i.e. HF, DFT(B3LYP) and finally MP2. The exploration resulted in 10 conformations with a relative energy difference smaller than 22 kJ.mol(-1) at the HF/3-21G level of theory. These conformations led to six different conformations after DFT(B3LYP) optimizations. Further optimization at the MP2/6-31++G** level of theory resulted in the same six conformations, all of them with an energy difference smaller than 11.5kJ.mol(-1). One conformation with an intramolecular H-bond was found which was energetically the most favorable conformation. The vibrational and thermodynamical features were calculated using the DFT and MP2 methodologies. In the experimental matrix-isolation FT-IR spectrum, the most stable conformation was dominant and at least two non-H-bonded conformations could be identified. An experimental rotamerization constant between the H-bonded and the other non-H-bonded conformations was estimated and appeared to agree reasonably well with the theoretical MP2 predictions. Some new spectral features of N-acetylproline compared to proline were discovered which might be used to discriminate between the acetylated and non-acetylated form.     

DFT/MM modeling of the five-membered ring in 3,6-anhydrogalactose derivatives and its influence on disaccharide adiabatic maps

Different conformations of methyl 3,6-anhydro-4-O-methyl-alpha-d-galactoside (1) and 3,6-anhydro-4-O-methylgalactitol (2) were studied by molecular mechanics (using the program mm3) and by quantum mechanical (QM) methods at the B3LYP/6-31+G( * *) and MP2/6-311++G( * *) levels, with and without solvent emulation. In 2, where the five-membered ring is free to move, two main stable conformations of this ring were found, identified as North (N) and South (S). The latter appears to be more stable, by either calculation, though the energy difference is reduced when emulating solution behavior. In order to find out the possible influence of a glycosidic bond over its shape, and to explain the marked NMR chemical shift displacements observed by opening of the ring, the adiabatic maps of two disaccharides carrying an analog of beta-galactoside linked to O-4 of 1 and 2 were generated. It was shown that the characteristics of the 3,6-AnGal terminal influence the characteristics of the map, especially at lower dielectric constants. On the other hand, different glycosidic angles also promote distinct stable conformations of the five-membered ring, changing from N to S, or even variants. Comparison with experimental results leads to the idea of highly flexible disaccharides, with variable values for both the five-membered ring and the glycosidic angles.     

Quantum mechanical and NMR spectroscopy studies on the conformations of the hydroxymethyl and methoxymethyl groups in aldohexosides

The potential energy surfaces of the hydroxymethyl and methoxymethyl groups in methyl hexopyranosides have been extensively studied, employing quantum mechanical calculations and high resolution NMR data. The structure and energy of the C-5-C-6 rotamers were calculated at the B3LYP level of the density functional theory (DFT). For all, geometry optimizations were carried out for 264 conformers of 16 methyl D-gluco- and methyl D-galactopyranoside derivatives 1-16 at the B3LYP/6-31G** level. For all calculated minima, single-point calculations were performed at the B3LYP/6-311++G** level. Solvent effects were considered using a self-consistent reaction field method. Values of the vicinal coupling constants 3J(H-5-H-6R), 3J(H-5-H-6S), 3J(C-4-H-6R), and 3J(C-4-H-6S) for methyl D-glucopyranosides, methyl D-galactopyranosides and their 6-O-methyl derivatives 9-16 were measured in two solvents, methanol and water. The calculated gg, gt, and tg rotamer populations of the hydroxymethyl and methoxymethyl groups in 9-16 agreed well with experimental data. The results clearly showed that the population of gg, gt, and tg rotamers is sensitive to solvent effects. It was concluded that the preference of rotamers in 1-16 is due to the hydrogen bonding and solvent effects.     

DFT studies of the conversion of four mesylate esters during reaction with ammonia

The energetics of the Menshutkin-like reaction between four mesylate derivatives and ammonia have been computed using B3LYP functional with the 6-31+G** basis set. Additionally, MPW1K/6-31+G** level calculations were carried out to estimate activation barrier heights in the gas phase. Solvent effect corrections were computed using PCM/B3LYP/6-31+G** level. The conversion of the reactant complexes into ion pairs is accompanied by a strong energy decrease in the gas phase and in all solvents. The ion pairs are stabilized with two strong hydrogen bonds in the gas phase. The bifurcation at C2 causes a significant activation barrier increase. Also, bifurcation at C5 leads to noticeable barrier height differentiation. Both B3LYP/6-31+G** and MPW1K/6-31+G** activation barriers suggest the reaction 2 (2a?+?NH₃) to be the fastest in the gas phase. The reaction 4 is the slowest one in all environments.

Boat conformations synthesis, NMR spectroscopy, and molecular dynamics of methyl 4,6-O-benzylidene-3-deoxy-3-phthalimido-alpha-D-altropyranoside derivatives

Addition of the elements of phthalimide to methyl 2,3-anhydro-4,6-O-benzylidene-alpha-D-mannopyranoside (1) under fusion conditions has yielded methyl 4,6-O-benzylidene-3-deoxy-3-phthalimido-alpha-D-altropyranoside (2). The conformation of the pyranose ring of 2 has been shown to be non-chair by 1H NMR spectroscopy, in contrast to the conformations of related derivatives having smaller substituents at C-3. Molecular dynamics simulations of 2 in explicit chloroform-d solvent have indicated four principal conformational possibilities. Of these, the 7C5/1S5 chair/skew boat form 2d has the lowest potential energy, and is largely consistent with the observed vicinal 1H-1H NMR coupling constants.     

Conformational study of trinucleoside tetraphosphate d(pCpGpCp): Transition of right-handed form to left-handed form

Using the semiempirical potential functions, conformational energies of the model compounds DMP^?, d(pCp), d(pGp), and d(pCpGpCp) are calculated, and the B → Z transition is discussed along the pseudorotational path of the sugar ring. For dimethylmonophosphate anion, DMP^?, the energy contour map is presented and the stabilities of the phosphodiester conformations discussed. For the sugar ring without the base attached, the minimum energies for each sugar-puckering form are calculated along the pseudorotational path. The energy barrier of the interconversion between the C(3′)-endo form and the C(2′)-endo form is calculated to be about 2.0 kcal/mol. From the conformational energy calculations of the interconversions of mononucleoside diphosphates, d(pCp) and d(pGp), between the C(2′)-endo conformer and the C(3′)-endo conformer, the purine sugar segment is known to be more convertible than the pyrimidine sugar segment. The results also support the finding that the pseudorotational transition occurred with the O(1′)-endo form more easily than with the O(1′)-exo form. Based on the results of conformational studies of DMP^?, d(pCp), and d(pGp), a topological transition of the handedness of the model compound, d(pCpGpCp), is studied. The left-handed Z-form is found to be less stable by about 8.5 kcal/mol than is the right-handed B-form. The energy barrier of the Z → B transition is calculated to be about 17.4 kcal/mol. The contributions of the electrostatic and nonbonded energies to the energy barrier are discussed in connection with the conformation changes of the model compound, d(pCpGpCp).     

Docking studies on glycoside hydrolase Family 47 endoplasmic reticulum alpha-(1-->2)-mannosidase I to elucidate the pathway to the substrate transition state

Alpha-(1-->2)-mannosidase I from the endoplasmic reticulum (ERManI), a Family 47 glycoside hydrolase, is a key enzyme in the N-glycan synthesis pathway. Catalytic-domain crystal structures of yeast and human ERMan1s have been determined, the former with a hydrolytic product and the latter without ligands, with the inhibitors 1-deoxymannojirimycin and kifunensine, and with a thiodisaccharide substrate analog. Both inhibitors were bound at the base of the funnel-shaped active site as the unusual 1C4 conformer, while the substrate analog glycon is a 3S1 conformer. In the current study, AutoDock was used to dock alpha-D-mannopyranosyl-(1-->2)-alpha-D-mannopyranose with its glycon in chair (1C4,4C1), half-chair (3H2,3H4,4H3), skew-boat (OS2,3S1,5S1), boat (2,5B,3,OB,B1,4,B2,5), and envelope (3E,4E,E3,E4) conformations into the yeast ERManI active site. Both docked energies and forces on docked ligand atoms were calculated to determine how the ligand distorts to the transition state. From these, we can conclude that (1) both 1C4 and OS2 can be the starting conformers; (2) the most likely binding pathway is 1C4-->3H2-->OS2-->3,OB-->3S1-->3E; (3) the transition state is likely to be close to a 3E conformation.