Variation of packing density in particle deforming from spheres to cubes is studied. A new model is presented to describe particle deformation between different particle shapes. Deformation is simulated by relative motion of component spheres in the sphere assembly model of a particle. Random close packings of particles in deformation form spheres to cubes are simulated with an improved relaxation algorithm. Packings in both 2D and 3D cases are simulated. With the simulations, we find that the packing density increases while the particle sphericity decreases in the deformation. Spheres and cubes give the minimum (0.6404) and maximum (0.7755) of packing density in the deformation respectively. In each deforming step, packings starting from a random configuration and from the final packing of last deforming step are both simulated. The packing density in the latter case is larger than the former in two dimensions, but is smaller in three dimensions. The deformation model can be applied to other particle shapes as well.