Split green fluorescent protein (GFP) is a powerful tool for imaging of protein protein interactions in living cells, but molecular mechanisms of the folding and the assembly of split GFPs are poorly understood. Here, using a simple Go model that is based on the energy landscape theory, we performed comprehensive folding simulations of six split GFPs with different split points. Of the six, the fluorescence recovery was reported in four but not in the other two. In the simulations, we found that when the complete folding and were observed, the N-terminal fragment always folded earlier than the C-terminal fragment. The in silico folding rates of the split GFPs were larger for the four split GFPs that the fluorescence recovery was reported in literature. The stability of standalone N-terminal fragments were well-correlated with the folding rates of split GFPs. These suggest that the efficient folding and assembly of split GFPs are realized when the N-terminal fragment folds spontaneously with the central alpha-helix as a nucleation core and that the C-terminal fragment folding is coupled to the assembly to the preformed N-terminal fragment.