The unfolding behavior of ubiquitin under the influence of a stretching force recently was investigated experimentally by single-molecule constant-force methods. Many observed unfolding traces had a simple two-state character, whereas others showed clear evidence of intermediate states. Here, we use Monte Carlo simulations to investigate the force-induced unfolding of ubiquitin at the atomic level. In agreement with experimental data, we find that the unfolding process can occur either in a single step or through intermediate states. In addition to this randomness, we find that many quantities, such as the frequency of occurrence of intermediates, show a clear systematic dependence on the strength of the applied force. Despite this diversity, one common feature can be identified in the simulated unfolding events, which is the order in which the secondary-structure elements break. This order is the same in two-and three-state events and at the different forces studied. The observed order remains to be verified experimentally but appears physically reasonable. all-atom model ͉ force-induced unfolding ͉ Monte Carlo simulation T he 76-residue protein ubiquitin fulfills many important regulatory functions in eukaryotic cells through its covalent attachment to other proteins (1, 2). In many cases, the ubiquitin tag consists of a chain of ubiquitin domains (polyubiquitin), which is formed by linkages between an exposed lysine side chain of the last ubiquitin of a growing chain and the C terminus of a new ubiquitin. The fate of a polyubiquitin-tagged protein depends on the linkage. For example, Lys-48-C-linked polyubiquitin marks the protein substrate for proteasomal degradation (3).Recently, Fernandez and coworkers (4-7) and Chyan et al. (8) investigated the mechanical properties of polyubiquitin by single-molecule force spectroscopy. It was shown that Lys-48-Clinked as well as end-to-end (N-C)-linked polyubiquitin can withstand a stretching force; the average unfolding force was 85 pN for Lys-48-C linkage and Ϸ200 pN for N-C linkage (4). In these experiments, the polyubiquitin chains were pulled with a constant velocity. In another experiment on N-C-linked polyubiquitin, the stretching force was kept constant (6). At constant force, the fraction of unfolded ubiquitin domains was found to show an approximately single-exponential time dependence, as expected if the unfolding of individual domains is a simple Markovian two-state process. Nevertheless, the unfolding of individual domains sometimes occurred through intermediate states. The precise nature of these different unfolding pathways, including the structure of the intermediate states, remains to be determined. We stress that these intermediates are states along forced unfolding trajectories. To what extent there are significant folding intermediates for small proteins is a debated (9, 10), but different, issue. Here, we use Monte Carlo (MC) simulations to examine the unfolding of ubiquitin under a constant stretching force in atomic detail. Our calculations are b...