Molecular dynamics (MD) simulation using the reference potential spatial warping algorithm (REPSWA) is a promising method for computing the conformational equilibrium of a system with a rugged energy surface. Its effectiveness has previously been demonstrated using only a simplified model system of a hydrocarbon chain omitting nonbonded interaction terms from the potential energy function. To evaluate the applicability of REPSWA MD simulation to more realistic problems, we applied it to small peptide systems in an aqueous environment. Difficulties were encountered, however, forcing us to devise several modifications. We investigated their effectiveness in comparison to conventional constant-temperature and multicanonical MD simulations. We found that the sampling efficiency of the modified REPSWA MD after a careful optimization of its parameters was better than that of the constant-temperature MD and comparable to that of the multicanonical MD in several cases.