A material point method (MPM) for the analysis of two-step heating of metals subjected to ultrafast laser irradiation is developed based on the weak formulation of the two-temperature model (TTM). The electron and lattice subsystems in the metallic targets are both represented by a finite number of Lagrangian material points to discretize the TTM equations, while a background Eulerian grid mesh is used to solve the spatially discrete TTM equations for the electron and lattice temperatures carried by the material points. The verification, convergence, and robustness of the proposed MPM method are demonstrated using representative examples for femto/picosecond laser heating of a gold thin film. In addition to avoiding the numerical difficulties due to the mesh distortion and entanglement in the mesh-based method, the presented MPM algorithm automatically handles the adiabatic boundary condition without requiring additional boundary treatments and could achieve higher computational efficiency than the mesh-based approaches at a comparable level of accuracy, showing its advantage over the mesh-based methods for modeling energy transfer in metals irradiated by ultrafast lasers.