Two-photon spectroscopy of fluorescent proteins is a powerful bio-imaging tool characterized by deep tissue penetration and little damage. However, two-photon spectroscopy has lower sensitivity than one- photon microscopy alternatives and hence a protein with a large two-photon absorption cross-section is needed. We use time-dependent density functional theory (TD-DFT) at the B3LYP/6-31+G(d,p) level of theory to screen twenty-two possible chromophores that can be formed upon replacing the amino-acid Tyr66 that forms the green fluorescent protein (GFP) chromophore with a non-canonical amino acid. A proposed chromophore with a nitro substituent was found to have a large two-photon absorption cross-section (29 GM) compared to other fluorescent protein chromophores as determined at the same level of theory. Classical molecular dynamics are then performed on a nitro-modified fluorescent protein to test its stability and study the effect of the conformational flexibility of the chromophore on its two-photon absorption cross-section. The theoretical results show that the large cross-section is primarily due to the difference between the permanent dipole moments of the excited and ground states of the nitro-modified chromophore. This large difference is maintained through the various conformations assumed by the chromophore in the protein cavity. The nitro-derived protein appears to be very promising as a two-photon absorption probe.