As a direct simulation of a multistep proton transfer reaction involving protein residues, the proton relay shuttle between A and I forms of green fluorescent protein (GFP) is simulated in atomic detail by using a special molecular dynamics simulation technique. Electronic excitation of neutral chromophore in wild-type GFP is generally followed by excited-state proton transfer to a nearby glutamic acid residue via a water molecule and a serine residue. Here we show that the second and third transfer steps occur ultrafast on time scales of several tens of femtoseconds. Proton back-shuttle in the ground state is slower and occurs in a different sequence of events. The simulations provide atomic models of various intermediates and yield realistic rate constants for proton transfer events. In particular, we argue that the I form observed spectroscopically under equilibrium conditions may differ from the I form observed as a fast intermediate by an anti to syn rotation of the carboxyl proton of neutral Glu-222.B ecause of its unique photophysical properties-strong green fluorescence without need for an additional cofactorgreen fluorescent protein (GFP) has become a very powerful marker for gene expression, cellular localization, and dynamic intracellular events over the past 5 years (1). Its rich photophysical behavior was characterized quite well by a great variety of spectroscopic techniques (1-4), and its three-dimensional structure was determined by x-ray crystallography (5, 6). Wild-type GFP exhibits two absorption maxima ''A'' and ''B'' at 395 nanometers (nm) and at 475 nm (1) that are present roughly in a 6:1 ratio (1), and that correspond to forms of the protein with either neutral (A) or anionic (B) chromophores. Time-resolved fluorescence spectroscopy identified an additional form ''I*'' that is populated as a fast intermediate after excitation of the A form, A* 3 I* (2, 3). The mechanism of interconversion is likely caused by excited-state proton transfer from the chromophore to a nearby protein residue (2, 3) because the new form emits at a similar wave length as B* and because the transfer rate slows down upon deuteration (2); see Fig. 1 for an overview.Although ''I form'' originally denoted only this fast intermediate, it has become common to use ''I form'' as a term for GFP species that absorb at wavelengths red-shifted with respect to the B form. Surprisingly, recent hole-burning experiments demonstrated that an I form is already present in wild-type GFP samples at room temperature under equilibrium conditions (4). It is currently unclear whether this is the same I form observed in the time-resolved fluorescence experiments. Although the I form has so far not been described in atomic detail, a structural model of I was suggested on the basis of comparisons between the crystal structures of wild-type (A form) and mutant GFP proteins containing the B form (7, 8); see Fig. 2. In this model, Glu-222 is assumed to be protonated in the I form (7), thus requiring a rise of the pK a of Glu-222 by more than 4 ...