From a Newtonian potential–density pair we construct axisymmetric and spherical relativistic sources for a Majumdar–Papapetrou type conformastatic spacetime. As a first example, we considerer a family of Plummer–Hernquist type relativistic spherical sources, and as a second application we construct relativistic galaxy models from the Miyamoto–Nagai potential–density pair. The models also include dark matter. We study the equatorial circular motion of test particles around such configurations. Also, the stability of the orbits is analyzed for radial perturbation using an extension of the Rayleigh criterion. The models considered satisfy all the energy conditions.