A time-nonlocal quantum master equation coupled with a perturbative scheme to evaluate the third-order polarization in the phase-matching direction k(s) = -k(1) + k(2) + k(3) is used to efficiently simulate three-pulse photon-echo signals. The present method is capable of describing photon-echo peak shifts including pulse overlap and bath memory effects. In addition, the method treats the non-Markovian evolution of the density matrix and the third-order polarization in a consistent manner, thus is expected to be useful in systems with rapid and complex dynamics. We apply the theoretical method to describe one-and two-color three-pulse photon-echo peak shift experiments performed on a bacterial photosynthetic reaction center and demonstrate that, by properly incorporating the pulse overlap effects, the method can be used to describe simultaneously all peak shift experiments and determine the electronic coupling between the localized Q(y) excitations on the bacteriopheophytin (BPhy) and accessory bateriochlorophyll (BChl) in the reaction center. A value of J 250 cm(-1) is found for the coupling between BPhy and BChl.