  Eukaryotic cell motility involves a complex network of interactions between
biochemical components and mechanical processes. The cell employs this network
to polarize and induce shape changes that give rise to membrane protrusions and
retractions, ultimately leading to locomotion of the entire cell body. The
combination of a nonlinear reaction-diffusion model of cell polarization, noisy
bistable kinetics, and a dynamic phase field for the cell shape permits us to
capture the key features of this complex system to investigate several motility
scenarios, including amoeboid and fan-shaped forms as well as intermediate
states with distinct displacement mechanisms. We compare the numerical
simulations of our model to live cell imaging experiments of motile {\it
Dictyostelium discoideum} cells under different developmental conditions. The
dominant parameters of the mathematical model that determine the different
motility regimes are identified and discussed.
