Molecular motors are responsible for active transport and organization in the cell, underlying an enormous number of crucial biological processes. Dynein is more complicated in its structure and function than other motors. Recent experiments have found that, unlike other motors, dynein can take different size steps along microtubules depending on load and ATP concentration. We use Monte Carlo simulations to model the molecular motor function of cytoplasmic dynein at the single-molecule level. The theory relates dynein's enzymatic properties to its mechanical force production. Our simulations reproduce the main features of recent single-molecule experiments that found a discrete distribution of dynein step sizes, depending on load and ATP concentration. The model reproduces the large steps found experimentally under high ATP and no load by assuming that the ATP binding affinities at the secondary sites decrease as the number of ATP bound to these sites increases. Additionally, to capture the essential features of the step-size distribution at very low ATP concentration and no load, the ATP hydrolysis of the primary site must be dramatically reduced when none of the secondary sites have ATP bound to them. We make testable predictions that should guide future experiments related to dynein function. molecular motors ͉ theory ͉ simulations M olecular motors are responsible for active transport and organization in the cell, underlying an enormous number of crucial biological processes (1). There are three classes of molecular motors, kinesin, myosin, and dynein. Myosins move along actin filaments, kinesin moves toward the plus end of a microtubule (MT), and dynein moves toward the MT minus end. To understand these motors better, experimental studies of the function of the motors at the single-molecule level have been conducted (2-9), and these quantitative measurements have then been modeled theoretically by using coupled differential rate equations (10-15).Dynein is more complicated than kinesin or myosin. This complexity can be seen experimentally in the step-size distribution. The distribution of step sizes for kinesin and myosin are centered about a single value. For kinesin, the average step size is 8 nm (16), and for myosin-V it is 37 nm (17). Dynein, in contrast, displays four different step sizes (8,16,24, and 32 nm), depending on the load and ATP concentration (18). By ''load'' we mean an external backward force applied to a polystyrene bead by using an optical trap, when a dynein molecule is attached to the bead and hauls it along a microtubule. If there is no load, at low ATP the dynein moves with a mixture of 24-and 32-nm steps (18). When load is present, if sufficient ATP is available, the dynein can decrease its step size to 8 nm and produce force up to 1.1 pN. If the load is large enough so that the motor is no longer able to move the bead, we refer to the load as the stalling force. The stalling force increases linearly with ATP concentration before saturating at 1.1 pN (18). In this work, we use Monte...