CARBOHYDRATES - UNITED ATOM AMBER-ASTERISK PARAMETERIZATION OF PYRANOSES AND SIMULATIONS YIELDING ANOMERIC FREE-ENERGIES

The success of molecular modeling using classical potential functions (i.e force field calculations) rests heavily on the availability of sp ecific, high-quality parameters that accurately describe the gas phase potential surface of the molecular system under study, on solvent mod els that reliably reproduce the effect of the medium, and on simulatio n methods that sample all significantly populated conformations of the entire system with the correct statistical weights. In this paper we present a set of molecular mechanics parameters that were developed us ing nb initio molecular orbital calculations to model pyranoses in the context of the AMBER force field in the molecular modeling package M acroModel 5.0. These parameters were tailored to reproduce the quantum mechanical conformational energies of certain small molecules that we re taken as models for common substructures in monosaccharides. Solven t was included as the GB/SA continuum model for water. The sampling pr oblem was solved for these systems using the recently described MC(JBW )/SD simulation method that facilitates interconversion of predetermin ed conformational minima by alternating between smart Monte Carlo and stochastic dynamics steps. A series of such MC(JBW)/SD simulations usi ng the new carbohydrate parameters was used to calculate anomeric alph a,beta ratios (and thus anomeric free energy differences) for tetrahyd ropyran derivatives and the pyranose monosaccharides glucose, methyl g lucoside, mannose, methyl mannoside, galactose, 2-deoxyglucose, and N- acetylglucosamine. In all cases, the simulations converged within 1 ns to yield anomeric free energies that are within 0.4 kcal/mol of the e xperimentally determined anomeric free energies in water.