The solution conformation of methyl alpha-lactoside has been studied through molecular mechanics calculations using the AMBER/Homans, CVFF and CFF91 force fields, and compared to NMR nuclear Overhauser data. Steady-state and transient nuclear Overhauser effects (NOEs) have been interpreted in terms of the ensemble average distribution of conformers. The NOEs have been analysed using the complete relaxation matrix approach for a rigid and isotropic motion model. The molecular mechanics calculations have been performed at two dielectric constants (i.e. epsilon = 1 and 80 debyes, or epsilon = r and 80 debyes) in an exhaustive way, and, in some cases, have been complemented by specific calculations at intermediate epsilon values. Relaxed energy maps and adiabatic surfaces have been generated for the different dielectric constants. The probability distribution of conformers has been estimated from these steric energy maps. Molecular dynamics simulations in vacuo have also been performed. Our results indicate that the beta-(1-->4) glycosidic linkage shows some fluctuations between three low-energy regions, although it spends about 90% of its time in the region close to the global minimum. The observed conformation of methyl alpha-lactoside seems to be closer to that predicted by CVFF, although the AMBER/Homans results are also in qualitative agreement with the experimental data.