Relative stability of alternative chair forms and hydroxymethyl conformations of β-d-glucopyranose

Abstract The relative energies of two hydroxymethyl conformers for each of the two chair forms (4C1 and 1C4) of β- d -glucose were calculated at much more complete levels of quantum mechanical (QM) electronic structure theory than previously, and relative free energies in solution were calculated by adding vibrational, rotational, and solvent effects. The gas-phase results are based on very large basis sets (up to 624 contracted basis functions), and the coupled cluster method for electron correlation. Solvation Model 4 was used to calculate the effects of hydration or nonpolar solvation. Molecular mechanics (MM) and QM electronic structure theory have been applied to analyze the factors contributing to the relative energies of these conformers. Relative energies varied widely (up to 35 kcal/mol) depending on theoretical level, and several levels of theory predict the experimentally unobserved 1C4 ring conformation to be the lower in energy. The highest level calculations predict the 4C1 chair to be lower in free energy by about 8 kcal/mol, and we also find that the gauche+ (gt) conformer of 4C1 is lower than the trans (tg) conformer. Low-energy structures optimized by either quantum mechanical or molecular mechanical methods were commonly characterized by multiple intramolecular hydrogen bonds. Superior hydrogen bonding geometries are available in the 13C4 chair, but are counteracted by increased steric repulsions between axial substituents; MM calculations also indicate increased torsional strain in the 1C4 chair. Manifestations of greater steric strain in the calculated 1C4 structures compared to the 4C1 structures include longer ring bonds, a larger bond angle at the ring oxygen atom, and smaller puckering amplitudes. The MM and QM 4C1 structures compare well with each other and with available X-ray diffraction data. The largest discrepancies between the two kinds of models occur for geometric parameters associated with the anomeric center - the QM structure agrees better with experiment. Greater differences between QM and MM structures are observed for 1C4 structures, especially in the relative orientations of hydroxyl groups serving as hydrogen bond acceptors. In water, the 4C1 chairs are better solvated than the 4C4 chairs by about 5 to 9 kcal/mol because of both larger polarization free energies and improved hydrogen bonding interactions with the first solvation shell. In (a hypothetical) n-hexadecane solution, the 4C1 chairs are better solvated by about 2 to 4 kcal/mol both because of larger polarization free energies and because the larger solvent accessible surface areas of the 4C1, conformers allow increased favorable dispersion interactions. The differential polarization free energies are associated primarily with the hydroxyl groups; the greater steric congestion in the 1C4 chairs reduces opportunities for favorable dielectric screening.

[1]  B. B. Neto,et al.  Comparative study of atomic charges derived from infrared intensities and from electrostatic potentials using AM1 and MNDO wavefunctions , 1992 .

[2]  G. Chałasiński,et al.  Perturbation analysis of the supermolecule interaction energy and the basis set superposition error , 1992 .

[3]  Gábor I. Csonka,et al.  Analysis of the core‐repulsion functions used in AM1 and PM3 semiempirical calculations: Conformational analysis of ring systems , 1993, J. Comput. Chem..

[4]  M. A. Fabian,et al.  The origin of the anomeric effect: conformational analysis of 2-methoxy-1,3-dimethylhexahydropyrimidine , 1994 .

[5]  N. P. Franks,et al.  Where do general anaesthetics act? , 1978, Nature.

[6]  J. T. Ham,et al.  The crystal and molecular structure of methyl β-cellobioside–methanol , 1970 .

[7]  K. Schulten,et al.  Erratum: Generalized moment expansion for observables of stochastic processes in dimensions d > 1: Application to Mössbauer spectra of proteins [J. Chem. Phys. 84, 4015 (1986)] , 1987 .

[8]  D. Truhlar,et al.  Direct dynamics calculations with NDDO (neglect of diatomic differential overlap) molecular orbital theory with specific reaction parameters , 1991 .

[9]  Klaus Gundertofte,et al.  A comparison of conformational energies calculated by molecular mechanics (MM2(85), Sybyl 5.1, Sybyl 5.21, and ChemX) and semiempirical (AM1 and PM3) methods , 1991 .

[10]  C. Cramer,et al.  Molecular orbital theory calculations of aqueous solvation effects on chemical equilibria [Erratum to document cited in CA115(19):207215n] , 1991 .

[11]  P. C. Hariharan,et al.  The influence of polarization functions on molecular orbital hydrogenation energies , 1973 .

[12]  R. Franck.  A revision of the value for the anomeric effect , 1983 .

[13]  D. Rohrer,et al.  The structure of gentiobiose , 1980 .

[14]  R M Venable,et al.  Molecular dynamics simulations of a lipid bilayer and of hexadecane: an investigation of membrane fluidity. , 1993, Science.

[15]  G. A. Jeffrey,et al.  The refinement of the crystal structures of -D-glucose and cellobiose , 1968 .

[16]  Carlos Aleman,et al.  Suitability of the PM3‐derived molecular electrostatic potentials , 1993, J. Comput. Chem..

[17]  J. Tomasi,et al.  The influence of the solvent on the conformational energy differences due to the anomeric effect , 1991 .

[18]  George C. Shields,et al.  Ability of the PM3 quantum‐mechanical method to model intermolecular hydrogen bonding between neutral molecules , 1993, J. Comput. Chem..

[19]  Ernest L. Eliel,et al.  Stereochemistry of Organic Compounds , 1962 .

[20]  A. Cerezo,et al.  Use of a general purpose force-field (MM2) for the conformational analysis of the disaccharide α-D-galactopyranosyl-(1→3)-β-D-galactopyranose , 1994 .

[21]  J. Kroon,et al.  Hydrogen-bond geometry around sugar molecules: comparison of crystal statistics with simulated aqueous solutions☆ , 1990 .

[22]  J. Praly,et al.  Influence of solvent on the magnitude of the anomeric effect , 1987 .

[23]  James J. P. Stewart,et al.  Bond indices and valency , 1973 .

[24]  Saul Wolfe,et al.  On the magnitudes and origins of the “anomeric effects”, “exo-anomeric effects”, “reverse anomeric effects”, and CX and CY bond -lengths in XCH2YH molecules☆ , 1979 .

[25]  R. Dwek,et al.  Analysis of glycoprotein-associated oligosaccharides. , 1993, Annual review of biochemistry.

[26]  J. Pople,et al.  Self‐Consistent Molecular‐Orbital Methods. I. Use of Gaussian Expansions of Slater‐Type Atomic Orbitals , 1969 .

[27]  A. Vasella,et al.  Glycosylidene Carbenes. Part 14. Glycosidation of partially protected galactopyranose-, glucopyranose-, and mannopyranose-derived vicinal diols , 1994 .

[28]  Vincenzo Mollica,et al.  Group contributions to the thermodynamic properties of non-ionic organic solutes in dilute aqueous solution , 1981 .

[29]  J. Stewart.  Optimization of parameters for semiempirical methods II. Applications , 1989 .

[30]  Michael H. Abraham,et al.  Thermodynamics of solute transfer from water to hexadecane , 1990 .

[31]  D. Weaver,et al.  An examination of intermolecular and intramolecular hydrogen bonding in biomolecules by AM1 and MNDO/M semiempirical molecular orbital studies , 1991 .

[32]  A. Finkelstein,et al.  Water and nonelectrolyte permeability of lipid bilayer membranes , 1976, The Journal of general physiology.

[33]  J. Pople,et al.  Self‐Consistent Molecular‐Orbital Methods. IX. An Extended Gaussian‐Type Basis for Molecular‐Orbital Studies of Organic Molecules , 1971 .

[34]  Leo Radom,et al.  The application of ab initio molecular orbital theory to the anomeric effect. A comparison of theoretical predictions and experimental data on conformations and bond lengths in some pyranoses and methyl pyranosides , 1972 .

[35]  G. D. Hawkins,et al.  Erratum: Comment by Prof. Cramer, Mr. Hawkins and Prof. Truhlar in Faraday symposium 29 on potential-energy surfaces and organic reaction paths (Journal of the Chemical Society, Faraday Transactions (1994) 90 (1802-1804) , 1994 .

[36]  Bruce Tidor,et al.  Solvent effect on the anomeric equilibrium in D-glucose: a free energy simulation analysis , 1991 .

[37]  R. Hooft,et al.  Molecular dynamics study of conformational and anomeric equilibria in aqueous D-glucose , 1993 .

[38]  J. Phillip Bowen,et al.  An ab initio study of dimethoxymethane protonation and its relevance to glycoside hydrolysis , 1989 .

[39]  Eamonn F. Healy,et al.  Development and use of quantum mechanical molecular models. 76. AM1: a new general purpose quantum mechanical molecular model , 1985 .

[40]  C. Cramer.  Anomeric and reverse anomeric effects in the gas phase and aqueous solution , 1992 .

[41]  A. Abe,et al.  Conformation and conformational energies of dimethoxymethane and 1,1-dimethoxyethane , 1990 .

[42]  E. Davidson,et al.  A possible definition of basis set superposition error , 1994 .

[43]  D. Davies,et al.  Intramolecular hydrogen bonding in 1'-sucrose derivatives determined by SIMPLE proton NMR spectroscopy. , 1986, Journal of the American Chemical Society.

[44]  John W. Brady,et al.  Molecular dynamics simulations of .alpha.-D-glucose in aqueous solution , 1989 .

[45]  Norman L. Allinger,et al.  Molecular mechanics parameters , 1994 .

[46]  K. Merz,et al.  Conformational preferences for hydroxyl groups in substituted tetrahydropyrans , 1992 .

[47]  Miguel A. Ríos,et al.  Semiempirical study of compounds with O‐HO intramolecular hydrogen bond , 1992 .

[48]  Alfred D. French,et al.  Analysis of the ring‐form tautomers of psicose with MM3 (92) , 1994, J. Comput. Chem..

[49]  C. Cramer.  Where is the unpaired electron in the phosphoranyl radicals H3PS− and H3PSH? , 1993 .

[50]  Timothy Clark,et al.  Multicenter point charge model for high‐quality molecular electrostatic potentials from AM1 calculations , 1993, J. Comput. Chem..

[51]  Gabriel Cuevas,et al.  Recent studies of the anomeric effect , 1992 .

[52]  Peter J. Reilly,et al.  Modeling of aldopyranosyl ring puckering with MM3 (92) , 1994 .

[53]  Hiroshi Ohrui,et al.  1H-NMR studies of (6r)- and (6s)-deuterated d-hexoses: assignment of the preferred rotamers about C5C6 bond of D-glucose and D-galactose derivatives in solutions , 1984 .

[54]  Norman L. Allinger,et al.  Molecular mechanics. The MM3 force field for hydrocarbons. 2. Vibrational frequencies and thermodynamics , 1989 .

[55]  Kenneth M. Merz,et al.  A force field for monosaccharides and (1 → 4) linked polysaccharides , 1994, J. Comput. Chem..

[56]  A. French,et al.  Overlapping anomeric effects in a sucrose analogue. , 1993, Carbohydrate research.

[57]  Walter Thiel,et al.  Semiempirical methods: current status and perspectives , 1988 .

[58]  Donald G. Truhlar,et al.  Polarization of the nucleic acid bases in aqueous solution , 1992 .

[59]  W. L. Jorgensen,et al.  Monte-Carlo Results for the Effect of Solvation on the Anomeric Equilibrium for 2-Methoxytetrahydropyran , 1994 .

[60]  George P. Ford,et al.  Atomic charges derived from a fast and accurate method for electrostatic potentials based on modified AM1 calculations , 1994, J. Comput. Chem..

[61]  S. F. Boys,et al.  The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors , 1970 .

[62]  G. A. Jeffrey,et al.  Stereographic representation of the cremer-pople ring-puckering parameters for pyranoid rings☆ , 1979 .

[63]  I. Tvaroška,et al.  Ab initio molecular orbital calculation of carbohydrate model compounds. 2. Conformational analysis of axial and equatorial 2-methoxytetrahydropyrans , 1994 .

[64]  N. L. Allinger,et al.  Molecular Mechanics (MM3). Calculations of Furan, Vinyl Ethers, and Related Compounds , 1993 .

[65]  Norman L. Allinger,et al.  A molecular mechanics force field (MM3) for alcohols and ethers , 1990 .

[66]  Structure and biosynthesis of prokaryotic glycoproteins. , 1988, Biochimie.

[67]  J. V. Hunt,et al.  Oxidative alterations in the experimental glycation model of diabetes mellitus are due to protein-glucose adduct oxidation. Some fundamental differences in proposed mechanisms of glucose oxidation and oxidant production. , 1993, The Biochemical journal.

[68]  A. Aabloo,et al.  Miniature crystal models of cellulose polymorphs and other carbohydrates. , 1993, International journal of biological macromolecules.

[69]  O. Runquist,et al.  Conformational analysis of some 2-alkoxytetrahydropyrans , 1968 .

[70]  C. S. Ewig,et al.  Ab Initio computed molecular structures and energies of the conformers of glucose , 1992 .

[71]  Influence of polar medium on the anomeric effect.: Quantum chemical study using the ab initio continuum model , 1991 .

[72]  L. Poppe,et al.  The rigidity of sucrose : just an illusion ? , 1992 .

[73]  F M Richards,et al.  Areas, volumes, packing and protein structure. , 1977, Annual review of biophysics and bioengineering.

[74]  J. Vliegenthart,et al.  A 1H-NMR and MD study of intramolecular hydrogen bonds in methyl β-cellobioside , 1992 .

[75]  Donald G. Truhlar,et al.  Quantum Chemical Conformational Analysis of 1,2-Ethanediol: Correlation and Solvation Effects on the Tendency To Form Internal Hydrogen Bonds in the Gas Phase and in Aqueous Solution , 1994 .

[76]  Norman L. Allinger,et al.  Molecular mechanics. The MM3 force field for hydrocarbons. 1 , 1989 .

[77]  R. Bartlett,et al.  A full coupled‐cluster singles and doubles model: The inclusion of disconnected triples , 1982 .

[78]  J. Phillip Bowen,et al.  An ab initio study (6-31G*) of transition states in glycoside hydrolysis based on axial and equatorial 2-methoxytetrahydropyrans , 1991 .

[79]  Jan Kroon,et al.  Solvent effect on the conformation of the hydroxymethyl group established by molecular dynamics simulations of methyl‐β‐D‐glucoside in water , 1990 .

[80]  N. Yasuoka,et al.  The crystal and molecular structure of a 3:2 mixture of laminarabiose and O-α-d-glucopyranosyl-(1→3)-β-d-glucopyranose , 1977 .

[81]  Warren J. Hehre,et al.  AB INITIO Molecular Orbital Theory , 1986 .

[82]  David J. Giesen,et al.  A SEMIEMPIRICAL QUANTUM MECHANICAL SOLVATION MODEL FOR SOLVATION FREE ENERGIES IN ALL ALKANE SOLVENTS , 1995 .

[83]  Cornelis Altona,et al.  The anomeric effect: Ab‐initio studies on molecules of the type XCH2OCH3 , 1990 .

[84]  S. Angyal.  Conformational analysis in carbohydrate chemistry. I. Conformational free energies. The conformations and α : β ratios of aldopyranoses in aqueous solution , 1968 .

[85]  A. French,et al.  AB INITIO-MIA AND MOLECULAR MECHANICS STUDIES OF THE DISTORTED SUCROSE LINKAGE OF RAFFINOSE , 1994 .

[86]  J. Pople,et al.  Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules , 1972 .

[87]  M. Márquez,et al.  The energy components of the anomeric effect for 2-methoxytetrahydropyran. An experimental comparison of the gas phase and solutions , 1994 .

[88]  G. A. Jeffrey,et al.  The application of ab initio molecular orbital theory to structural moieties of carbohydrates , 1974 .

[89]  S. J. Angyal,et al.  The Composition and Conformation of Sugars in Solution , 1969 .

[90]  Norman L. Allinger,et al.  Directional hydrogen bonding in the MM3 force field. I , 1994 .

[91]  P. Kollman,et al.  Comparison of ab initio, semiempirical, and molecular mechanics calculations for the conformational analysis of ring systems , 1992 .

[92]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[93]  B. Lee,et al.  The interpretation of protein structures: estimation of static accessibility. , 1971, Journal of molecular biology.

[94]  G. Thatcher.  The Anomeric effect and associated stereoelectronic effects , 1993 .

[95]  I. Tvaroška,et al.  Ab Initio Molecular Orbital Calculation on Carbohydrate Model Compounds. 1. The Anomeric Effect in Fluoro and Chloro Derivatives of Tetrahydropyran , 1994 .

[96]  U. Salzner,et al.  Generalized anomeric effects and hyperconjugation in CH2(OH)2, CH2(SH)2, CH2(SeH)2, and CH2(TeH)2 , 1993 .

[97]  A. J. Kirby,et al.  The Anomeric Effect and Related Stereoelectronic Effects at Oxygen , 1982 .

[98]  George P. Ford,et al.  New approach to the rapid semiempirical calculation of molecular electrostatic potentials based on the am1 wave function: Comparison with ab initio hf/6‐31g* results , 1993, J. Comput. Chem..

[99]  Donald G. Truhlar,et al.  Systematic study of basis set superposition errors in the calculated interaction energy of two HF molecules , 1985 .

[100]  L. Lerner,et al.  Observation of hydroxyl protons of sucrose in aqueous solution: no evidence for persistent intramolecular hydrogen bonds , 1992 .

[101]  K. Hirotsu,et al.  The Crystal and Molecular Structure of -Lactose , 1974 .

[102]  David J. Giesen,et al.  Class IV charge models: A new semiempirical approach in quantum chemistry , 1995, J. Comput. Aided Mol. Des..

[103]  C. Cramer,et al.  General parameterized SCF model for free energies of solvation in aqueous solution , 1991 .

[104]  J. Brady,et al.  A revised potential-energy surface for molecular mechanics studies of carbohydrates. , 1988, Carbohydrate research.

[105]  David J. Giesen,et al.  General Semiempirical Quantum Mechanical Solvation Model for Nonpolar Solvation Free Energies. n-Hexadecane , 1995 .

[106]  Kenneth M. Merz,et al.  Study of hydrogen bonding interactions relevant to biomolecular structure and function , 1992 .

[107]  G. A. Jeffrey,et al.  Crystal Structure of Glucosamine Hydrobromide , 1939, Nature.

[108]  Robert B. Hermann,et al.  Theory of hydrophobic bonding. II. Correlation of hydrocarbon solubility in water with solvent cavity surface area , 1972 .

[109]  R. Zhbankov.  Vibrational spectra and structure of mono- and polysaccharides , 1992 .

[110]  R. E. Reeves.  The Shape of Pyranoside Rings , 1950 .

[111]  T. H. Dunning.  Gaussian basis sets for use in correlated molecular calculations. I. The atoms boron through neon and hydrogen , 1989 .

[112]  D. Cremer,et al.  General definition of ring puckering coordinates , 1975 .

[113]  Donald G. Truhlar,et al.  Quantum Chemical Conformational Analysis of Glucose in Aqueous Solution , 1993 .