The material point method (MPM), as proposed by Sulsky et al. (1994), has been developed to simulate large deformation problems. A continuum body is divided into a finite number of subregions represented by Lagrangian material points, while the governing equations are formulated and solved with the Eulerian computational grid. Since this grid can be chosen arbitrarily, mesh tangling does not appear in the MPM. As evidenced in the literature, however, the original MPM has not been used to simulate fluid-saturated porous media response under diverse loading conditions. In the present study, the MPM is coupled with the finite difference method (FDM) for simulating fluid-saturated porous solid response under quasi-static loading. The governing differential equations are formulated based on the equations of motion for the solid-fluid mixture and the continuity equation for the pore fluid, in which the equations of motion and the continuity equation are discretized by the MPM and FDM, respectively. To demonstrate the proposed numerical scheme, deformation analyses of water-saturated porous elastic solid are performed, and the numerical results are discussed via a comparison with those obtained by the finite element method and the analytical approach.