A procedure is demonstrated for summing the Moller-Plesset many-body perturbation expansion based on the ability of quadratic summation approximants to locate branch point singularities in the complex plane of the perturbation parameter. Accuracy comparable to that from CCSDT coupled-cluster calculations is obtained using fourth-order perturbation theory.