in the past 15 years , green fluorescent protein ( gfp ) has changed from a nearly unknown protein to a commonly used molecular imaging tool in biology , chemistry , genetics , and medicine . in 2006 , more than 10  000 papers
gfps and gfp - like proteins ( i.e. , chromoproteins and fluorescent proteins ) are particularly useful due to their stability and the fact that the chromophore ( see figure 1 ) is formed in an autocatalytic cyclization of the 65syg67 sequence that does not require a cofactor .
this means that unlike most other bioluminescent reporters , gfp fluoresces in the absence of any other proteins , substrates , or cofactors .
furthermore , it appears that fusion of gfp to a protein does not alter the function or location of the protein . by changing residue 66 and/or the amino acid residues around the chromophore
(n1c1c2c3 ) and (c1c2c3c4 ) dihedral angles of the gfp chromophore . in the protein ,
r1 is gly67 and r2 is ser65 , and in hbdi , an often - used model compound , r1 = r2 = ch3 . in  one - bond flips ( -obf ) the dihedral rotation occurs around the  torsional angle , in a -obf it is around the  dihedral angle , and in a positively correlated hula - twist ( + ht ) the  and  dihedral angles concertedly rotate in the same direction ( as shown above ) , while in a negatively correlated hula twist ( ht ) they concertedly rotate in opposite directions . a plot of the  and  dihedrals for a perfectly correlated negative ht will have a slope of 1 .
if the chromophore cavity is complementary with a planar chromophore , then the  and  best fit line will pass through the origin and all the  and  angles will be centered around the origin .
a nonzero intercept along the  or  axis ( see , for example , figure 4 ) or  and  dihedrals centered in quadrant ii (  > 0 ;
 < 0 ) or quadrant iv (  < 0 ;  > 0 ; see , for example , figure 8) are indications of a cavity that is not complementary with a planar chromophore . the fluorescent emission of the chromophore within gfp occurs with high efficiency ( quantum yield fl = 0.8 ) and a respectable fluorescence lifetime (   3 ns).(19 ) when the protein is denatured the fluorescence yield decreases by at least 3 orders of magnitude.(20 ) model compounds of the chromophore do not fluoresce in solution ( quantum yield fl < 10 ) , unless the rotation of the arylalkene bond is restrained.(21 ) fluorescence can , however , be obtained by lowering the temperature to 77 k ; this freezes the solution .
it has been suggested that twisting between the phenolate and imadazolidinone groups of the chromophore is the mechanism for the ultrafast fluorescence quenching internal conversion process .
neutral form ( a ) of the chromophore , with the phenolic oxygen protonated , can convert to the anionic species ( b ) by going through an intermediate state ( i ) .
the change from forms a to i is solely a protonation change , while the change from i to b is a conformational change with most changes occurring at thr203 .
upon excitation of the a state , an excited - state proton transfer ( espt ) occurs in which the proton is transferred from the chromophore to glu222 in a time scale of the order of picoseconds . following radiative relaxation from the excited - state intermediate ( i * ) ,
the systems returns to the ground state a through the ground - state intermediates i1 and i2.(29 ) excitation of the anionic b state results in direct emission from the excited state ( b * ) at 482 nm . recently
, a nonfluorescent dark state , state c , has been observed that is distinct from states a and b and absorbs at higher energies.(30 ) the c state , perhaps the neutral trans form of the chromophore , may be populated by nonradiative decay from a * and it may be depopulated by excitation to the excited c * state with transcis isomerization to repopulate state a. fluorescent blinking has been ascribed to nonadiabatic crossing and conversions between the neutral and anionic states,(31 ) or a possible dark ( z ) zwitterionic state . while the ground - state minimum of the gfp chromophore is close to planar , this is not necessarily so for the excited state ; in fact , in some cases the excited state has an energy minimum with a twisted chromophore in which both rings are 90 to each other.(32 ) according to quantum mechanical calculations , the ground and excited states for the  one - bond flip ( obf ) and hula twist ( ht ) in the neutral form ( a ) and the  obf in the zwitterionic form come very close to each other .
it has been proposed that in the absence of the protein matrix , which surrounds the chromophore and prevents twisting , this process can lead to fluorescence quenching internal crossing;(32 ) see figure 3 .
recent calculations on the gfp chromophore model compound hbdi suggest that the anionic form of hbdi may also undergo a -obf that leads to a favored radiationless decay channel , which is particularly efficient in solvent.(34 ) ( a ) model compounds of the gfp chromophore in the ground state ( s0 ) can be excited to the first singlet state ( s1 ) in which a ht or obf can freely occur . upon reaching the perpendicularly twisted conformation , fluorescence quenching nac ( nonadiabatic crossing ) occurs .
( b ) in the ground state ( s0 ) , the residues surrounding the gfp chromophore exert a twisting force on the chromophore (  ) . upon excitation , the conjugation across the ethylenic bridge of the chromophore is reduced and it will twist ; however , the protein matrix prevents the chromophore from reaching the perpendicularly twisted conformation (  ) and fluorescence quenching internal crossing is prevented .
recently , the results from a number of interesting experiments that provide evidence for chromophore twisting and/or fluorescence quenching internal crossing ( ic ) have been published .
for example , a molecular dynamics simulation of the chromophore of cyan fluorescent protein found that the average  dihedral angle was about 0 ; however , when the surrounding protein was considered in the simulation the average shifted by approximately 5 , showing that the protein matrix of cfp twists the chromophore . on the basis of their calculations , the authors concluded that  the driving force for this twist comes from the strong short - range repulsion by four residues ( ile167 , val150 , phe165 , and thr203 ) surrounding the average position of the chromophore.(46 ) photoswitching fluorescent proteins can be switched back and forth between the naturally occurring green state and a dark state by 405 nm irradiation ( e.g. , dronpa , mtfp0.7 , kfp1 ) .
a cistrans isomerization of the chromophore has been proposed as the structural basis for the photoswitching observed in dronpa.(47 ) this is supported by the fact that mutating either val157 or met159 with smaller residues accelerates photoswitching , presumably by decreasing steric hindrance to cistrans isomerization.(48 )
the m159 t and v157 g mutations also decrease the quantum yield of dronpa from 0.85 to quantum yields of 0.23 and 0.77 , respectively.(48 ) recently , it has been suggested that adoption of a trans configuration can not solely be responsible for the nonfluorescent form.(49 ) on the basis of nmr analyses , miyawaki et al .
propose that  the fluorescence of the protein is regulated by the degree of flexibility of the chromophore but is not necessarily accompanied by cistrans isomerization.(49 ) interestingly , gfp is not unique  in photoactive yellow protein ( pyp ) the protein matrix also prevents the chromophore from adopting a completely planar structure . in pyp
the asymmetric proteinchromophore interaction probably serves as the initial accelerant for the light induced photocycle,(50 ) which ultimately leads to a cistrans isomerization.(51 ) the chromophore in wild - type gfp is planar due its extended  system ( see figure 1 ) ; however , the energy barrier to deformation is low and the protein matrix in wild - type gfp exerts some strain away from planarity .
when the chromophore is computationally permitted to freely rotate , it will adopt a conformation that complements the protein matrix .
recently , we have used computational methods to show that wild - type gfp is not an anomaly and that all gfp and gfp - like proteins in the protein databank have a protein matrix that is not complementary with a planar chromophore.(54 ) in most cases the freely rotating chromophore undergoes  rotations of at least 20. in some cases these rotations are accompanied by an equal but opposite rotation of the  dihedral angle ( a negatively correlated ht ) .
none of the proteins examined have a cavity that causes a rotation solely around the  dihedral angle.(54 ) these calculations were done by minimizing , with freely rotating  and  dihedral angles , the crystal structure of 38 gfp analogues and mutants found in the pdb .
they found the energy minimum conformation of a freely rotating chromophore in the protein matrix of the gfp mutant or gfp - like protein examined ; however , they did not provide any information about the range of low - energy conformations available to a freely rotating chromophore . to get this information
, we ran molecular dynamics simulations of some of the interesting gfp mutants in the protein databank.(56 ) by running molecular dynamics simulations , with freely rotating  and  dihedrals , we have been able to determine the range of conformations available to chromophores with complete rotational freedom .
the coordinates of three crystal structures ( 1gfl,(57 ) 1myw,(58 ) 2emd(59 ) ) were obtained from the protein data bank ( pdb),(56 ) and hydrogen atoms were added to protein and solvent atoms as required . the opls_2005 force field of macromodel v9.0016(60 ) was used .
starting structures for mutants for which no crystal structure has been determined were calculated by graphically mutating a known structure and undertaking a conformational search .
the flexible dihedral angles of all the side chains of residues that are within 8.00  of the chromophore were randomly rotated by between 0 and 180 , and all solvent molecules in that sphere were randomly rotated and translated by between 0 and 1.00  in each monte carlo ( mc ) step.(63 ) 15  000 mc steps were taken in each search .
structures within 50 kj / mol of the lowest energy minimum were kept and a usage - directed method(62 ) was used to select structures for subsequent mc steps .
structures found in the conformational search were considered unique if the least - squared superimposition of equivalent non - hydrogen atoms found one or more pairs separated by 0.25  or more . the lowest energy structure obtained in the monte carlo torsional and molecular position variation search
the final structures obtained from the mc search ( monte carlo torsional and molecular position variation and large - scale low mode ) or fully minimized pdb structures were used to initiate md simulations with freely rotating  and  dihedral angles ( v1 = v2 = v3 = 0.000 ) .
the predynamics simulation was set for 100 ps and the full md simulation for 5000 ps .
the cambridge structure database ( csd ) v5.29 was released in january 2008 ; it comprises 436  384 small molecule crystal structures .
the area of the convex hull of the  and  graphs was calculated as was the smallest convex set containing the region .
molecular mechanics and dynamics calculations have been used to examine the steric environment of the chromophore in gfp in its ground state .
they are techniques that are based on classical physics and were designed to model structural and not electronic properties ; therefore , molecular mechanics and dynamics simulations of gfp can not examine the excited state of the chromophore .
quantum calculations are an excellent technique to determine the energy profiles of the ground and excited states of the chromophore ; however , they are cpu - intensive and it would be cost prohibitive to do conformational searches on the excited state of a series of gfp and gfp - like proteins .
therefore , in an attempt to supplement the quantum calculations , we have examined the conformational space available to the chromophore within gfp using a freely rotating chromophore that is an approximation ( based on qm calculations ) of the conformational space available to the chromophore in its excited state . as shown in figure 2 ,
gfp adopts two states , the neutral state ( a ) of the chromophore with the phenolic oxygen protonated , and the anionic species ( b )
the change from forms a to i is solely a protonation change , while the change from i to b is a conformational change with most changes occurring at thr203.(67 ) in order to see whether the gfp chromophore has different dihedral freedom in forms a , b , and i , we conducted molecular dynamical simulations with a freely rotating chromophore , as described in the .
not surprisingly , table 1 and figure 4 show that forms a , b , and i of gfp have similar dihedral freedom and that all undergo similar negatively correlated ( same slopes of best fit lines in figure 4 ) hula twists .
plot of the  vs  dihedral angles ( see figure 1 for nomenclature ) for the 2000 gfp - a ( pink ) , gfp - b ( green ) , and gfp - i ( light - blue ) structures obtained in the freely rotating molecular dynamics simulation . however , there are subtle differences ; forms a and i are closer to the planar ground - state conformation , while the b form adopts conformations further from planarity when given rotational freedom around the  and  dihedrals ( intercept further from  = 0 ) ; and the anionic b form has the most freedom ( largest area , and  and  range ) , while the neutral a form has the least freedom .
the a form has robust hydrogen bonds between both arg96 and gln94 and the imidazolidinone carbonyl group that remain intact throughout the simulations . in the intermediate ( i )
form these hydrogen bonds are supplemented by an additional hydrogen bond between the phenolic oxygen of the chromophore and his148 ( see table 2 ) .
the main difference between the anionic b form and the intermediate form is that thr203 has to rotate in order to hydrogen bond to the phenolic oxygen in form i , otherwise they have the same hydrogen - bonding interactions throughout the simulation .
the y66h mutant of gfp exhibits blue fluorescence and has therefore been termed blue fluorescent protein .
the crystal structure of bfp ( not to be confused with blue fluorescent protein from aequorin bfp - aq(68 ) ) has been solved at ph 4.5(27 ) and ph 8.5;(59 ) the overall fold of the protein is identical to wild - type gfp and consists of a chromophore surrounded by an 11-stranded -barrel .
while gfp has absorption maxima at 395 and 475 nm and emits at 508 nm , bfp absorbs at 382 nm and emits at 448 nm.(9 ) unfolding of bfp results in an absorption red shift of 15 nm,(27 ) and quantum mechanical calculations suggest that the 15 nm shift might be due to the protein - induced nonplanarity of the chromophore.(69 ) blue fluorescent protein ( bfp ) has a much lower fluorescence quantum yield than gfp ( fl = 0.20 vs 0.80 ) .
it has been suggested that this is due to the fact that his66 ( bfp - chromophore ) forms fewer hydrogen bonds with the surrounding protein than tyr66 ( gfp - chromophore ) does , and that the smaller imadazole ring ( his66 ) in bfp may have more conformational freedom than the larger phenol ( tyr66 ) , which leads to more intersystem crossing.(27 ) in a very elegant series of experiments , boxer et al .
have shown that the fluorescence quantum yield of blue fluorescent protein increases from fl = 0.20 to 0.35 when the pressure is increased from atmospheric pressure to 570 mpa.(70 )  analysis of the fluorescence lifetimes in the picosecond and nanosecond regimes reveals that the enhancement of the fluorescence quantum yield is due to the inhibition of fast quenching processes .
temperature - dependent fluorescence measurements reveal two barriers ( 19 and 3 kj / mol , respectively ) for the transition into nonfluorescing states .
these steps are probably linked with dissociation of the hydrogen bond between the chromophore and his148 or an intervening water molecule and to the barrier for chromophore twisting in the excited state , respectively.(70 ) in order to establish the consequences of the y66h mutation on the flexibility of the chromophore , molecular dynamics simulations of the ph 8.5 bfp crystal structure ( pdb code 2emd with neutral imidazole rings for his66 and his148 ) with a freely rotating chromophore were run .
[ we thank one of the referees for suggesting the ph 8.5 ( 2emd ) over the ph 4.5 ( 1bfp ) structure .
] figure 5 and table 1 show that the chromophore in bfp has significantly more rotational freedom than that available to the chromophore in gfp , and that the bfp cavity is less complementary to a planar chromophore than the gfp cavity is ; that is , the bfp sampled structures are further from a planar chromophore (  =  = 0 ) and the intercept for the best fit line through the bfp structures is further from the origin than that of gfp .
plot of the  vs  dihedral angles for the 2000 bfp ( blue ) and gfp - a ( green ) structures obtained in the freely rotating molecular dynamics simulation .
analysis of the 2000 bfp structures revealed a hydrogen bond between the nh of the his66 imidazole ring and glu222 , which is retained through - out the entire simulation ( see figure 7 ) .
the smaller imadazole ring of the bfp chromophore results in more dihedral freedom for the chromophore in bfp than phenol in gfp ; however , there seems to be no large difference in the number and stability of hydrogen bonds formed by the chromophore in the two fps .
recently , two new blue fluorescent proteins ( azurite and a5 ) with enhanced brightness and photostability were created .
the methodology applied to find the brighter bfp mutants was based on the concept that replacing the residues surrounding the chromophore with bulkier amino acids would constrain the chromophore s motion and thereby increase the proteins brightness .
the crystal structure of a5 has yet to be solved . in order to find a starting conformation for the md simulation , the 2emd crystal structure
was graphically mutated to a5 and a thorough conformational search was undertaken ( see ) .
given the high structural similarity between all the solid state structures of gfp mutants in the protein databank , a thorough conformational search should find the lowest energy conformations of a5 by making the relevant mutations to the 2emd structure .
figure 6 and table 1 show that the conformational space available to the chromophore in bfp is larger than that in a5 and that the bfp conformations tend to be further from planarity ; i.e. , for a5 median = 2.92 , median = 2.22 and intercept of best - fit line is 4.62 vs for bfp median = 27.59 , median = 10.18 and intercept of best - fit line 14.27 .
similar results were found for md simulations carried out from all of the five lowest energy conformational families sampled in the conformational search of a5 and in a md simulation initiated from an a5 structure obtained by taking the 2emd structure , graphically mutating it to a5 , and thoroughly minimizing the structure without undertaking a conformational search .
plot of the  vs  dihedral angles for the 2000 bfp ( dark blue ) and a5 ( pink ) structures obtained in the freely rotating molecular dynamics simulation .
the n of the his66 imidazole ring hydrogen bonds with water305 ( 2emd numbering ) in all the sampled structures .
this water connects his66 to the surface of the protein via a robust hydrogen - bonding network to his148 and is also hydrogen bonded to glu222 . in half the structures sampled , the n of the his66 imidazole ring hydrogen bonds to water302 which is in an extensive hydrogen - bonding network with val68 , arg224 , and gln69 .
in contrast to 2emd there are no hydrogen - bonding interactions with the n hydrogen ; see figure 7 .
hydrogen - bonding interactions ( --- ) for 2emd ( top ) and a5 ( bottom ) chromophore .
these results seem to indicate that the enhanced brightness of a5 is at least in part due to the fact that its chromophore movement is restricted relative to bfp .
this restriction can be due to both steric factors and the increased hydrogen - bonding networks in a5 . in the first rationally designed mutant based on the crystal structure of gfp - s65 t ,
it was decided to mutate t203 into a tyrosine so that it could  stack with the phenolic group in the chromophore .
( 73 ) the resultant mutant , yfp , is red - shifted by 16 nm relative to gfp - s65 t and does indeed have a  stacking interaction between the chromophore and tyr203.(74 ) on the basis of the crystal structure of yfp , remington et al .
proposed that the red shift was due to the additional polarizability of the -stacked tyr203 , the hydrogen - bond pattern around the chromophore , or an out - of - plane distortion of the chromophore ( in analogy with out - of - plane distortions in porphyrin systems).(74 ) comparison with the t203v mutant revealed that the t203y substitution leads to a significant i - form population and only about 10 nm of the shift can be ascribed to the  interaction.(75 ) a new varient of yfp , venus , with improved brightness and maturation properties as well as a reduced environmental sensitivity was developed and crystallized.(58 ) since some researchers have argued that the presence of tyr203 leads to a decrease in chromophore flexibility ( van der waals volume for tyr = 141  vs 93  for thr ) , and others argue that this is not necessarily so,(32 ) we have examined the conformational space available to a freely rotating chromophore in yfp ( pdb code 1myw ) with the phenol of tyr66 in both the protonated and unprotonated forms . both table 1 and figure 8
show that the chromophore ( in the phenolic form ) in gfp and yfp has differing rotational freedom . in yfp
the chromophore has approximately 1.2 times more dihedral freedom than in gfp and although the dihedrals undergo negatively correlated hula - twisting in both proteins , in yfp the average conformation sampled is further from planarity than that found in gfp ; i.e. , the majority of the sampled conformations are located in quadrant iv .
plot of the  vs  dihedral angles ( see figure 1 for nomenclature ) for the 2000 structures obtained in the freely rotating molecular dynamics simulation of gfp ( dark blue ) and yfp ( pink ) in their phenolic forms .
the majority of the yfp conformations sampled are nonplanar ( i.e. , in quadrant iv  < 0 and  > 0 ) . during the entire simulation the hydrogen bonds between tyr66 and ser205 , arg96 and the imidazolone carbonyl , as well as between tyr203 and water354 remain stable , while those between tyr203 and gln69 were less stable .
we also examined the distance between the centroid of tyr66 ( the choromophore ) and tyr203 .
the distance is fairly short and does not change much during the simulation ( minimum = 3.43  , maximum = 4.18  , average = 3.75  , sd = 0.103  ) while the angle between the planes of the phenols varies between 0.0 and 22.5 ( average = 6.9 , sd = 4.11 ) .
the short distance between the phenol rings is indicative of  stacking , which may be responsible for the red shift observed in yfp ; however , the variability in the angle between the planes of the phenols shows that tyr203 does not significantly restrict the rotational freedom of the excited chromophore .
the cambridge structural database ( csd)(77 ) v5.29 comprises 436  384 small molecules crystal structures .
it has 5551 structures with two phenol rings in separate molecules that have centroids within 5.00  of each other .
the closest ones are 3.16  from each other , and only 832 structures are within 3.75  ; i.e. , the average distance between the phenol centroids in the freely rotating yfp .
a plot of the distance vs the angle between the phenol planes in our freely rotating simulation and the csd structures shows that the protein matrix of yfp restricts the distance between tyr66 and tyr203 from being larger than 4.18  , but that the angles between the two phenols have as much freedom in yfp as in the small molecule structures with intermolecular distances of less than 4.18  ( see figure 9 ) .
it does not seem as if the increased fluorescence lifetime of yfp is due to a decrease in the dihedral freedom of the chromophore .
it is more likely due to electronic effects such as those proposed by jung et al.(78 ) they suggested that the anionic form of the chromophore can be described by two mesomeric forms , the benzoidal and the quinoidal structure .
the benzoidal form with the majority of its negative charge on the phenol is stabilized by hydrogen bonding to thr203 , which is not possible in yfp and therefore the t203y mutation favors the quinoidal form , which may have a longer fluorescence lifetime.(78 ) distance vs angle between the phenol planes in yfp (  ) and all structures in the csd that have two unconnected phenols separated by less than 5.00  ( + ) .
some have suggested that the crystal structures of strongly fluorescent gfp and gfp - like proteins have their chromophores in a cis configuration , while those of the nonfluorescent proteins are in a trans configuration .
however , the trans - planar form of eqfp611 is an important exception : it is fluorescent.(82 ) remington et al .
have suggested that coplanarity might be more important than isomer configuration in determining the efficiency of fluorescence .
most current theories ( see ) suggest that a restriction in the rotational freedom of the fluorescent protein chromophore will lead to an increase in fluorescence brightness and/or quantum yield .
our calculations show that this is the case for the systems examined . for bfp , a5 , yfp , and gfp
, there is an inverse correlation between the dihedral freedom of the chromophore and the quantum yield ; see table 1 . in all simulations ,
the protein matrix is not complementary with a planar chromophore , and in all cases the freely rotating chromophore undergoes a negatively correlated hula twist ( also known as a bottom hula twist mechanism ) .
 green fluorescent protein ( gfp ) and gfp - like fluorescent proteins owe their photophysical properties to an autocatalytically formed intrinsic chromophore . according to quantum mechanical calculations , 
 the excited state of chromophore model systems has significant dihedral freedom , which may lead to fluorescence quenching intersystem crossing . 
 molecular dynamics simulations with freely rotating chromophoric dihedrals were performed on green , yellow , and blue fluorescent proteins in order to model the dihedral freedom available to the chromophore in the excited state . 
 most current theories suggest that a restriction in the rotational freedom of the fluorescent protein chromophore will lead to an increase in fluorescence brightness and/or quantum yield . according to our calculations , the dihedral freedom of the systems studied ( bfp > a5 > yfp > gfp ) increases in the inverse order to the quantum yield . in all simulations , 
 the chromophore undergoes a negatively correlated hula twist ( also known as a bottom hula twist mechanism ) . 
Introduction
Experimental Section
Results and Discussion
Conclusion