We report here a comparative study of the molecular conformational energy landscape generated using the mutually orthogonal Latin squares (MOLS) method, molecular dynamics (MD), and Monte Carlo (MC) simulation. The MOLS method, as described earlier from our laboratory, uses an experimental design technique to rapidly and exhaustively sample the low energy conformations of a molecule. MD and MC simulations have been used to perform similar tasks. In the comparison reported here, the three methods were applied to a pair of neuropeptides, namely Met- and Leu-enkephalin. A set of 1500 conformations of these enkephalins were generated using these methods with CHARMM22 force field, and the resulting samples were analyzed to determine the extent and nature of coverage of the conformational space. The results indicate that the MOLS method samples a larger number of possible conformations and identifies conformations closer to the experimental structures than the MD and MC simulations.