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-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.     

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.     

Flexibility of the pyranose ring in α- and β-D-glucoses

Conformational energies of α- and β-D -glucopyranoses were computed by varying all the ring bond angles and torsional angles using semiempirical potential functions. Solvent accessibility calculations were also performed to obtain a measure of solvent interaction. The results indicate that the ⁴C₁ (D ) chair is the most favored conformation, both by potential energy and solvent accessibility criteria. The ⁴C₁ (D ) chair conformation is also found to be somewhat flexible, being able to accommodate variations up to 10° in the ring torsional angles without appreciable change in energy. Observed solid-state conformations of these sugars and their derivatives lie in the minimum-energy region, suggesting that the substituents and crystal field forces play a minor role in influencing the pyranose ring conformation. Theory also predicts the variations in the ring torsional angles, i.e., CCCC < CCCO < CCOC, in agreement with the experimental results. The boat and twist-boat conformations are found to be at least 5 kcal mol^?1 higher in energy compared to the ⁴C₁ (D ) chair, suggesting that these forms are unlikely to be present in a polysaccharide chain. The ¹C₄ (D ) chair has energy intermediate between that of the ⁴C₁ (D ) chair and that of the twist-boat conformation. The calculated energy barrier between ⁴C₁ (D ) and ¹C₄ (D ) conformations is high—about 11 kcal mol^?1.     

Solid and solution state conformations of (+/-)-3-O-acetyl-1,2:4,5-di-O-isopropylidene-allo-inositol and (+/-)-3-O-acetyl-1,2:4,5-di-O-isopropylidene-6-O-methyl-allo-inositol

The synthesis and conformational studies of (+/-)-3-O-acetyl-1,2:4,5-di-O-isopropylidene-allo-inositol and (+/-)-3-O-acetyl-1,2:4,5-di-O-isopropylidene-6-O-methyl-allo-inositol are described. Solid state conformations of the title compounds have been studied by solving their X-ray crystal structures. The inositol ring in both the compounds deviate considerably from the ideal chair conformation to flattened chair conformation in the solid state. Their conformations in solution were studied by the use of 1H NMR spectroscopy. These conformational analyses revealed that the title compounds adopt similar conformations in solid and solution states irrespective of the solvent polarity.     

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.     

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.     

Solid and solution state conformation of 1L-1-O-acetyl-2,3:5,6-di-O-isopropylidene-chiro-inositol

The X-ray crystal structure of 1L-1-O-acetyl-2,3:5,6-di-O-isopropylidene-chiro-inositol is described. The inositol ring deviates considerably from the ideal chair conformation to a flattened chair. A comparison of its conformation in solution with that in solid was made by the use of 1H NMR. This conformational analysis revealed that the title compound adopts similar conformations in solid state and in solution states irrespective of solvent polarity.     

4-Membered metallodithiophosphinate rings - flat or puckered? A comparison of two crystal structures with computational and literature data

Two crystalline forms of (dithiodiphenylphosphinate)(phenyl)(triphenylphosphine)-palladium(II) (C₃₆H₃₀P₂PdS₂), one without solvent, the other containing THF (C₄H₈O), are obtained after reaction of sodium diphenyldithiophosphinate with (phenyl) (bis-triphenylphosphine) palladium(II) chloride and crystallisation from two different solvent mixtures. The molecular structures, as determined by single crystal X-ray diffraction, differ in the planarity of the 4-membered palladium dithiophosphinate rings. The experimental conformations have been compared to the conformations of four-membered metal-S₂P rings reported in the Cambridge Structural Database. A flat conformation is more common than a puckered one. DFT calculations at the B3LYP level of theory indicate that the flat conformation of a model metallodithiophosphinate ring is very slightly lower in energy (1.2 kcal/mol) than the puckered conformation.     

Comparison of kifunensine and 1-deoxymannojirimycin binding to class I and II alpha-mannosidases demonstrates different saccharide distortions in inverting and retaining catalytic mechanisms

Mannosidases are key enzymes in the eukaryotic N-glycosylation pathway. These enzymes fall into two broad classes (I and II) and are characteristically different in catalytic mechanism, sequence, and structure. Kifunensine is an alkaloid that is a strong inhibitor against class I alpha-mannosidases but is only a weak inhibitor against class II alpha-mannosidases. In this paper, the 1.80 A resolution crystal structure of kifunensine bound to Drosophila melanogaster Golgi alpha-mannosidase II (dGMII) is presented. Kifunensine adopts a (1,4)B boat conformation in the class II dGMII, which contrasts the (1)C(4) chair conformation seen in class I human endoplasmic reticulum alpha1,2 mannosidase (hERMI, PDB ). The observed conformations are higher in conformational energy than the global minimum (4)C(1) conformation, although the conformation in hERMI is closer to the minimum, as supported by an energy calculation. Differing conformations of 1-deoxymannojirimycin were also observed: a (4)C(1) and (1)C(4) conformation in dGMII and hERMI, respectively. Thus, these two alpha-mannosidase classes distort these inhibitors in distinct manners. This is likely indicative of the binding characteristics of the two different catalytic mechanisms of these enzymes.     

Experimental and quantum chemical calculational studies on 2-[(4-Fluorophenylimino)methyl]-3,5-dimethoxyphenol

The Schiff base compound, 2-[(4-Fluorophenylimino)methyl]-3,5-dimethoxyphenol, has been synthesized and characterized by IR, electronic spectroscopy, and X-ray single-crystal determination. Molecular geometry from X-ray experiment of the title compound in the ground state have been compared using the Hartree-Fock (HF) and density functional method (B3LYP) with 6–31G(d) basis set. Calculated results show that density functional theory (DFT) and HF can well reproduce the structure of the title compound. The energetic behavior of the title compound in solvent media has been examined using B3LYP method with the 6–31G(d) basis set by applying the polarizable continuum model (PCM). The total energy of the title compound decrease with the increasing polarity of the solvent. By using TD-DFT and TD-HF methods, electronic absorption spectra of the title compound have been predicted and a good agreement with the TD-DFT method and the experimental ones is determined. In addition, DFT calculations of the title compound, molecular electrostatic potential (MEP), natural bond orbital (NBO), and thermodynamic properties were performed at B3LYP/6–31G(d) level of theory.     

An investigation of the pyranose ring interconversion path of alpha-L-idose calculated using density functional theory

The interconversion pathways of the pyranose ring conformation of alpha-L-idose from a (4)C1 chair to other conformations were investigated using density functional calculations. From these calculations, four different ring interconversion paths and their transition state structures from the (4)C1 chair to other conformations, such as B(3,O), and (1)S3, were obtained. These four transition-state conformations cover four possible combinations of the network patterns of the hydroxyl group hydrogen bonds (clockwise and counterclockwise) and the conformations of the primary alcohol group (tg and gg). The optimized conformations, transition states, and their intrinsic reaction coordinates (IRC) were all calculated at the B3LYP/6-31G** level. The energy differences among the structures obtained were evaluated at the B3LYP/6-311++G** level. The optimized conformations indicate that the conformers of (4)C1, (2)S(O), and B(3,O) have similar energies, while (1)S3 has a higher energy than the others. The comparison of the four transition states and their ring interconversion paths, which were confirmed using the IRC calculation, suggests that the most plausible ring interconversion of the alpha-L-idopyranose ring occurs between (4)C1 and B(3,O) through the E3 envelope, which involves a 5.21 kcal/mol energy barrier.     

The conformation of dolichol

An understanding of the natural conformation of dolichol is important for the elucidation of the mechanism of protein glycosylation and dolichol's other as yet undisclosed biological functions. Since the molecular mechanics method has been shown to be well suited for the prediction of alcohol and alkene conformations, we have employed it to study the conformations of apparent least energy of dolichol-19 and smaller polymers of isoprene, namely, squalene, trans,trans-farnesol, and cis,cis-farnesol. Additionally, the small-angle X-ray scattering (SAXS) method was employed to determine the validity of the apparent least energy conformer of dolichol-19 derived by the molecular mechanics method. The results indicate that the solution conformation of dolichol-19 is comprised of a central coiled region flanked by two arms. The central coiled region has two and a half turns of dimensions 9.84 x 16.55 x 51.66 A3. The arms of dimensions 3.99 x 5.89 x 17.47 A3 and 4.49 x 9.23 x 11.14 A3 are approximately diametrically opposed. Measurement of the intrinsic viscosity of dolichol in both isopentyl alcohol and oleyl alcohol showed that the natural conformation of dolichol is capable of increasing solution fluidity (i.e., lowering solution viscosity). Thus, while examination of the conformation of dolichol in a membrane-mimetic solvent by SAXS is not possible, the quantitative measure of the effect of dolichol on solution viscosity (and thus solution fluidity) is possible. The results are consistent with dolichol acting as a membrane-fluidizing agent and provide the first quantitative measure of the effect of dolichol on solution fluidity of a membrane-mimetic solvent.     

3-Amino-2-piperidinone-6-carboxylic acid. Molecular structures of cis and trans benzoyl derivatives

The cis (2a) and trans (2b) isomers of methyl 3-benzamido-2-piperidinone-6-carboxylate (Apca) were prepared and separated by fractional recrystallizations. Proton n.m.r. studies in dimethylsulfoxide solution indicate that the six-membered lactam ring adopts a distorted chair conformation with an equatorially oriented benzamido substituent in both 2a and 2b. The carboxyl function also is equatorially oriented in the trans isomer 2b, but is disposed axially in the cis isomer 2a. In the crystal structure, the six-membered lactam ring of 2a is clearly in a boat conformation with the benzamido and carboxyl functions attached to the two apex carbon atoms equatorially. The trans isomer, 2b, exists as two crystallographically independent, conformationally distinct molecules in one unit cell. The lactam ring in both molecules adopts a distorted chair conformation, as is the case in solution, with both the benazamido and carboxyl functions attached equatorially. The rotameric orientation for the endocyclic lactam differs between the two molecules. Both structures show evidence of C-H...O hydrogen bond formation intermolecularly in the solid state. This ability, along with the distinctive conformational features of Apca, may be exploitable in the design of unique features of polypeptides.     

Molecular modeling of poly(ethylene oxide) model cofactors; 1,3,6-tri-O-galloyl-β-<Emphasis Type="SmallCaps">d</Emphasis>-glucose and corilagin

The most stable structures of two poly(ethylene oxide) (PEO) model cofactors, beta-1-O-galloyl-3,6-( R)-hexahydroxydiphenoyl- d-glucose (corilagin) and 1,3,6-tri-O-galloyl-beta- d-glucose (TGG), are calculated using molecular modeling and PM3 semiempirical molecular orbital theories. The theoretical PM3 structures agree with interpreted structures from experimental NMR; the glucopyranose ring of corilagin has a boat and TGG a chair conformation, for which the heats of formation, torsion angles, distances, van der Waals surface, and the infrared spectra are calculated.     

Can current force fields reproduce ring puckering in 2-O-sulfo-α-l-iduronic acid? A molecular dynamics simulation study

The monosaccharide 2-O-sulfo-α-l-iduronic acid (IdoA2S) is one of the major components of glycosaminoglycans. The ability of molecular mechanics force fields to reproduce ring-puckering conformational equilibrium is important for the successful prediction of the free energies of interaction of these carbohydrates with proteins. Here we report unconstrained molecular dynamics simulations of IdoA2S monosaccharide that were carried out to investigate the ability of commonly used force fields to reproduce its ring conformational flexibility in aqueous solution. In particular, the distribution of ring conformer populations of IdoA2S was determined. The GROMOS96 force field with the SPC/E water potential can predict successfully the dominant skew-boat to chair conformational transition of the IdoA2S monosaccharide in aqueous solution. On the other hand, the GLYCAM06 force field with the TIP3P water potential sampled transitional conformations between the boat and chair forms. Simulations using the GROMOS96 force field showed no pseudorotational equilibrium fluctuations and hence no inter-conversion between the boat and twist boat ring conformers. Calculations of theoretical proton NMR coupling constants showed that the GROMOS96 force field can predict the skew-boat to chair conformational ratio in good agreement with the experiment, whereas GLYCAM06 shows worse agreement. The omega rotamer distribution about the C5-C6 bond was predicted by both force fields to have torsions around 10°, 190°, and 360°.     

Molecular dynamics simulations of synthetic peptide folding

Because the time scale of protein folding is much greater than that of the widely used simulations of native structures, a detailed report of molecular dynamics simulations of folding has not been available. In this study, we Included the average solvent effect in the potential functions to simplify the calculation of the solvent effect and carried out long molecular dynamics simulations of the alanine-based synthetic peptides at 274 K. From either an extended or a randomly generated conformation, the simulations approached a helix-coil equilibrium in about 3 ns. The multiple minima problem did not prevent helix folding. The calculated helical ratio of Ac-AAQAAAAQAAAAQAAY-NH₂ was 47%, in good agreement with the circular dichroism measurement (about 50%). A helical segment with frayed ends was the most stable conformation, but the hydrophobic interaction favored the compact, distorted helix-turn-helix conformations. The transition between the two types of conformations occurred in a much larger time scale than helix propagation. The transient hydrogen bonds between the glutamine side chain and the backbone carbonyl group could reduce the free energy barrier of helix folding and unfolding. The substitution of a single alanine residue in the middle of the peptide with valine or glycine decreased the average helical ratio significantly, in agreement with experimental observations. © 1996 Wiley-Liss, Inc.     

Inclusion of solvation free energy with molecular mechanics energy: alanyl dipeptide as a test case.

A combined force field of molecular mechanics and solvation free energy is tested by carrying out energy minimization and molecular dynamics on several conformations of the alanyl dipeptide. Our results are qualitatively consistent with previous experimental and computational studies, in that the addition of solvation energy stabilizes the C5 conformation of the alanyl dipeptide relative to the C7.