Computational analysis of two similar neuropeptides yields distinct conformational ensembles

Conformational states and thermodynamic properties for two similar neuropeptides, GDPFLRF-NH(2) and GYPFLRF-NH(2), have been computed by Monte Carlo simulated annealing (MCSA) conformational searches and Metropolis Monte Carlo (MMC) calculations. These peptides were recently shown to have dramatically different conformations in solution by NMR [Edison et al., J Neuroscience 1999;19:6318-6326]. Final conformations of multiple independent MCSA runs were the starting points for MMC calculations, and conformations saved at intervals during MMC runs were characterized in terms of total energy, configuration entropy, side-chain fraction population, and ensemble average inter-nuclear distances. Without the use of any NMR data-generated pseudo-potentials, the present calculations were in excellent qualitative agreement with all previous NMR experimental data and provided a foundation by which to more quantitatively interpret the experimental NMR results. Proteins 2000;40:367-377.