Fluorescent proteins exhibit a very diverse range of photochemical behaviour, from efficient fluorescence through photochromism to photochemical reactivity. Remarkably this very diverse behaviour arises from chromophores which have very similar structures. Here we describe measurements and modelling of the excited state dynamics in the chromophores of GFP (HBDI) and the kindling fluorescent protein, KFP (FHBMI). The methods are ultrafast fluorescence spectroscopy with sub 50 fs time resolution and the modelling is based on the Smoluchowski equation. The excited state decays of both chromophores are very fast, longer for their anions than for the neutral form and independent of wavelength. Detailed studies show the mean fluorescence wavelength to be independent of time. The excited state decay times are also observed to be a very weak function of solvent polarity and viscosity. These results are modelled utilising recently calculated potential energy surfaces for the ground and excited states as a function of the twist coordinates about the two bridging bonds of the chromophore. For FHBMI and the scarce data on the neutral HBDI the calculations are not successful suggesting the need for refinement of these potential energy surfaces. For HBDI in methanol the simulation is successful provided a strong dependence of the radiationless decay rate on the coordinate is assumed. Such dependence should be included in future calculations of excited state dynamics. When the simulations are extended to more viscous solvents they fail to reproduce the observed weak viscosity dependence. The implications of these results for the nature of the coordinate leading to radiationless decay in the chromophore and for the photodynamics of fluorescent proteins are discussed.