Nanomanipulation of biomolecules by using single-molecule methods and computer simulations has made it possible to visualize the energy landscape of biomolecules and the structures that are sampled during the folding process. We use simulations and single-molecule force spectroscopy to map the complex energy landscape of GFP that is used as a marker in cell biology and biotechnology. By engineering internal disulfide bonds at selected positions in the GFP structure, mechanical unfolding routes are precisely controlled, thus allowing us to infer features of the energy landscape of the wild-type GFP. To elucidate the structures of the unfolding pathways and reveal the multiple unfolding routes, the experimental results are complemented with simulations of a self-organized polymer (SOP) model of GFP. The SOP representation of proteins, which is a coarse-grained description of biomolecules, allows us to perform forced-induced simulations at loading rates and time scales that closely match those used in atomic force microscopy experiments. By using the combined approach, we show that forced unfolding of GFP involves a bifurcation in the pathways to the stretched state. After detachment of an N-terminal ␣-helix, unfolding proceeds along two distinct pathways. In the dominant pathway, unfolding starts from the detachment of the primary N-terminal -strand, while in the minor pathway rupture of the last, C-terminal -strand initiates the unfolding process. The combined approach has allowed us to map the features of the complex energy landscape of GFP including a characterization of the structures, albeit at a coarse-grained level, of the three metastable intermediates.AFM experiments ͉ coarse-grained simulations ͉ cross-link mutants ͉ pathway bifurcation ͉ plasticity of energy landscape P rotein structures, which are astounding examples of selforganization in living systems, reach their folded states by navigating through a rugged energy landscape. Considerable progress has been made in understanding the folding mechanisms of small, single-domain proteins by using a combination of theory and experiments (1-5). Folding of many of these proteins can be approximately described as being two-state-like, that is, their energy landscape does not exhibit pronounced local minima corresponding to partially folded or misfolded structures. However, the folding energy landscapes of larger proteins, with complex topology, can be difficult to characterize because of the presence of multiple metastable intermediate structures (6,7). A detailed characterization of the structures of the intermediate states and the associated energetics is a challenge for the experimentalist. In part, the difficulty arises because complex proteins generally fold slowly and tend to aggregate in bulk experiments (8), a problem that can be avoided in mechanical unfolding of single proteins. Indeed, single-molecule mechanical methods have recently provided new possibilities for directly probing the energy landscapes of proteins and RNA because they can monito...