Abstract Cytoplasmic dynein is a two-headed molecular motor that moves to the minus end of microtubule (MT) using ATP hydrolysis free energy. By employing its two heads (motor domains), cytoplasmic dynein shows various bipedal stepping motions; the inchworm and hand-over-hand motions, as well as non-alternate steps of one head. However, the molecular basis to achieve such diverse stepping manners remains obscure. Here, we propose a kinetic model for bipedal motions of cytoplasmic dynein and performed Gillespie Monte Carlo simulations that reproduces most experimental data obtained to date. The model represents status of each motor domain as five states according to conformations, nucleotide- and MT-binding conditions of the domain. Also, the relative positions of the two domains were approximated by three discrete states. Accompanied by ATP hydrolysis cycles, the model dynein stochastically and processively moved forward in multiple steps via diverse pathways, including inchworm and hand-over-hand motions, same as experimental data. The model reproduced key experimental motility-related parameters including velocity and run-length as functions of ATP concentration and external force. Our model reveals that, in a typical inchworm motion, the leading domain moves via the ATP-dependent power-stroke of the linker coupled with a small change in the stalk angle, whereas the lagging domain moves via diffusion dragged by the leading domain. Moreover, the hand-over-hand motion in the model dynein clearly differs from that of kinesin by the usage of the power-stroke. Author Summary Cytoplasmic dynein is a two-headed molecular motor, which moves linearly and transports intra-cellar organelles along microtubules driven by ATP hydrolysis free energy. In contrast to other better-known molecular motors, such as kinesin, dynein is known to take various stepping motions including motions akin to human walking and inchworm-like motions. However, molecular mechanisms underpinning the diverse stepping motions are unclear. Here, based on recent high-resolution structure information and single-molecule motility assay data, we designed a kinetic model that explicitly include two heads, each of which makes ATP hydrolysis cycles and moves along the microtubules. Using the model, we performed Monte Carlo simulations. The simulation reproduced most of currently available experimental results. More importantly, the simulation suggested molecular mechanisms of various stepping motions. While stepping motions apparently resemble to those proposed before, once looking into details, we found the resulting mechanisms distinct from previously proposed ones in the usage of ATP and protein conformation changes coupled with stepping motions.