Journal of Chemical Physics, Vol.119, No.17, 9274-9279, 2003
Exploring repulsive interactions in a model helical peptide: A parallel tempering Monte Carlo study
By implementing the parallel-tempering algorithm to the canonical ensemble, the conformational changes of an isolated Ac-W(RAAAR)(5)A-NH2 model peptide were determined. The interparticle interactions were modeled using a minimalist potential, i.e., a beadlike model that uses harmonic oscillators to describe covalent interactions and modified Lennard-Jones potentials to model nonbonding interactions. In particular, the interactions between arginines are modeled by repulsive interactions, causing a stabilization of the alpha-helix structure at low temperatures. The conformational changes were identified by anomalies in the constant volume heat capacity as a function of temperature. Namely, the temperature at which the constant volume heat capacity reached a maximum in the transition region was associated with the temperature at which a conformational change occurred. The transitions were also characterized by computing the radius of gyration of the peptide and the most probable isomeric structure obtained at a given temperature. Three changes were observed at low temperatures and one at high temperature. The low temperature transitions were analogous to the peptide folding, whereas the high temperature transition was related to the peptide unfolding. The results obtained were compared with experimental data generated from isotope edited Fourier transform infrared spectroscopy and two-dimensional correlation analysis for a similar peptide containing salt bridge interactions. (C) 2003 American Institute of Physics.