An ensemble of exciton Hamiltonians for the amide-I band of the folded and unfolded states of a helical beta -heptapeptide is generated using a molecular dynamics (RID) simulation. The correlated fluctuations of its parameters and their signatures in two-dimensional (2D) vibrational echo spectroscopy are computed. This technique uses infrared pulse sequences to provide ultrafast snapshots of molecular structural fluctuations, in analogy with multidimensional NMR. The present study demonstrates that, by combining a method of calculating the vibrational Hamiltonian from MD snapshots and the nonlinear exciton equations (NEE), it may be possible to simulate realistic multidimensional IR spectra of chemically and biologically interesting systems.