Peptoid conformational free energy landscapes from implicit-solvent molecular simulations in AMBER

To test the accuracy of existing AMBER force field models in predicting peptoid conformation and dynamics, we simulated a set of model peptoid molecules recently examined by Butterfoss et al. (JACS 2009, 131, 16798-16807) using QM methods as well as three peptoid sequences with experimentally determined structures. We found that AMBER force fields, when used with a Generalized Born/Surface Area (GBSA) implicit solvation model, could accurately reproduce the peptoid torsional landscape as well as the major conformers of known peptoid structures. Enhanced sampling by replica exchange molecular dynamics (REMD) using temperatures from 300 to 800 K was used to sample over cis-trans isomerization barriers. Compared to (Nrch)5 and cyclo-octasarcosyl, the free energy of N-(2-nitro-3-hydroxyl phenyl)glycine-N-(phenyl)glycine has the most "foldable" free energy landscape, due to deep trans-amide minima dictated by N-aryl sidechains. For peptoids with (S)-N (1-phenylethyl) (Nspe) side chains, we observe a discrepancy in backbone dihedral propensities between molecular simulations and QM calculations, which may be due to force field effects or the inability to capture n --> n* interactions. For these residues, an empirical phi-angle biasing potential can "rescue" the backbone propensities seen in QM. This approach can serve as a general strategy for addressing force fields without resorting to a complete reparameterization. Overall, this study demonstrates the utility of implicit-solvent REMD simulations for efficient sampling to predict peptoid conformational landscapes, providing a potential tool for first-principles design of sequences with specific folding properties.