In atomic force spectroscopic studies of the elastomeric protein ubiquitin, the β-strands 1-5 serve as the force clamp. Simulations show how the rupture force in the force-induced unfolding depends on the kinetics of water molecule insertion into positions where they can eventually form hydrogen bonding bridges with the backbone hydrogen bonds in the force-clamp region. The intrusion of water into this region is slowed down by the hydrophobic shielding effect of carbonaceous groups on the surface residues of β-strands 1-5, which thereby regulates water insertion prior to hydrogen bond breakage. The experiments show that the unfolding of the mechanically stressed protein is nonexponential due to static disorder. Our simulations show that different numbers and/ or locations of bridging water molecules give rise to a long-lived distribution of transition states and static disorder. We find that slowing down the translational (not rotational) motions of the water molecules by increasing the mass of their oxygen atoms, which leaves the force field and thereby the equilibrium structure of the solvent unchanged, increases the average rupture force; however, the early stages of the force versus time behavior are very similar for our "normal" and fictitious "heavy" water models. Finally, we construct six mutant systems to regulate the hydrophobic shielding effect of the surface residues in the force-clamp region. The mutations in the two termini of β-sheets 1-5 are found to determine a preference for different unfolding pathways and change mutant's average rupture force. T here is a growing interest in the kinetic mechanisms by which elastomeric proteins unfold under force, as these proteins are involved in a wide variety of biological processes. Most elastomeric proteins (e.g., ubiquitin, I27, and fnIII) (1-3) share the common structural feature of having β-rich structures, populated by parallel β-strands. In atomic force spectroscopic studies of the elastomeric protein ubiquitin, the β-strands 1-5 serve as the force clamp under shear loading. Steered molecular dynamics (SMD) simulations (2, 4) showed that the shear loading force needs to break all four backbone hydrogen bonds (bbHBs) between the parallel β-strands 1-5 of ubiquitin in order to trigger protein unfolding (1,(5)(6)(7)(8). In SMD simulations of I27, Lu and Schulten observed that in the early stages of the pulling process individual bbHBs occasionally break and quickly reform, due to thermal fluctuations. During these fluctuations water molecules interact with the groups of broken bbHBs, but after several picoseconds these water molecules leave the region and the bbHBs reform (9). The bbHBs in the clamp region were thought to break simultaneously right before the sheet ruptures (2, 9, 10). In our SMD studies we observe that in ubiquitin under the pulling forces the four bbHBs in the β-strands 1-5 break sequentially from either the N or C termini, accompanied by the insertion of water molecules into the broken bbHBs to form stable bridging hydrogen bonds with...