sunlight absorption by biological molecules plays a crucial role in life development and sustains the earths energy cycle ; as a counterpart , photoprotective mechanisms are necessary to minimize the potential damage of the uv portion of the solar spectrum that reaches organisms.1 the interest in natural uv photoprotection research has increased lately , given substantial losses in the stratospheric ozone layer and the consequent increase in uv radiation levels in the troposphere.2 ,  3 one of the major strategies used by organisms to minimize uv - induced damage is screening by intra- or extracellular compounds that are able to block uv radiation.1 ,  4 ,  5 mycosporine - like amino acids ( maas ) and gadusol ( 3,5,6-trihydroxy-5-hydroxymethyl-2-methoxycyclohex-2-en-1-one ) belong to a series of natural uv - absorbing compounds that have been related to photoprotective and antioxidant functions in aquatic organisms.68 the photophysics and photochemistry of some of these secondary metabolites have been examined in aqueous solution.914 the main findings are : negligible fluorescence yields , high photostability , and efficient dissipation of the absorbed radiation energy as heat .
in fact , various reports indicate that maas are biosynthesized and accumulated in response to uv radiation , thus preventing damaging photons from reaching potential cellular targets and avoiding photosensitization processes by the favorable competition of thermal de - excitation pathways.57 ,  15 ,  16 however , the precise connection between these properties and the chemical structures of these compounds remains to be fully addressed . in that sense
, it is interesting to examine some trends in the correlation of the available photophysicochemical data with the structural features of this family of natural molecules .
taurine , and mycosporine  glutaminol  glucoside , all of which share a cyclohexenone core , are efficient antioxidants,6 ,  1719 whereas the imino - mycosporines shinorine and porphyra-334 are rather resistant to oxidation.20 the photodecomposition quantum yield of gadusol in aqueous solution is low , but is 100-fold greater than that of the conjugate anion gadusolate , the prevailing species under physiological ph.14 it is known that cyclic ,-unsaturated ketones produce relatively stabilized triplets if substituted at the -carbon atom;21 however , no triplet state from gadusolate could be observed by direct laser photolysis .
moreover , interaction between gadusolate and some triplet sensitizers leads to reductive quenching of the latter.14 on the contrary , sensitization experiments have given evidence for only the energy transfer to yield the excited triplet states of the imino - maas shinorine , porphyra-334 , and palythine.9 ,  10 ,  12 in this context , calculation of the deactivation pathways can provide valuable information for rationalization of the experimental results and consequently potentiate the future design of new synthetic photoprotectors .
recently , a computational study using the caspt2//casscf strategy was carried out on palythine , the simplest structure among imino - maas.22 the fast relaxation on the excited - state energy surface via a conical intersection ( ci ) point is evident from this study .
the main deformation of the structure corresponds to an out - of - plane movement of the imine moiety , while the adjacent carbon atoms approach each other and the alkene moiety remains almost planar .
this result rules out the previously proposed deactivation pathway by simple c = n bond isomerization.12 it is well established now that conical intersections mediate a variety of chemical events and have enormous relevance in biology , as they are responsible for the efficient internal conversion channels that recover dna bases after sunlight activation in an ultrafast way.23 ,  24 this suggests that dna bases and maa structures were naturally selected for their ability to dissipate energy.25 gadusol represents the basic molecular structure of the oxo - compounds mentioned above . at the same time
, it allows exploration of the influence of the substituents in the conjugated system by simply shifting the enol 
enolate equilibrium.26 a time - dependent density functional theory ( td - dft ) calculation of the transition energies for the enol form in water was carried out by arbeloa et  al.14 the results support the  * assignment of the excited singlet state accordingly with the generalized solvatochromic shifts of the absorption spectra .
however , no chemical calculation has been performed on gadusolate , the enolate form , which dominates the acid  base equilibrium under biological ph conditions .
herein we present the results of a theoretical study carried out with caspt2//casscf methodology in order to describe the underlying features responsible for the photoprotective role of gadusol . in turn
, this information could be useful to understand the photoprotective mechanisms selected by natural evolution and for the design and synthesis of new and more efficient artificial sunscreens .
the influence of the medium in the uv absorption of gadusol was recently reported.14 a bathochromic shift was observed as the polarity of the solvent increases .
for instance , the maximum in the absorption band changes from 268  nm in water to 263  nm in acetonitrile ; ph dependence was also observed in aqueous solution . for acidic solutions ( ph  2.5 )
maximal absorption was found at 268  nm , whereas at higher ph values ( ph7 ) an increase in band intensity was measured together with a bathochromic shift to 296  nm . these observations were explained on the basis of a displacement in the acid  base equilibrium between the neutral form ( gadusol ) present in acidic solutions and the enolate species ( gadusolate ) that dominates at ph  7 , although the exact structure of the gadusolate remained unclear .
therefore , as an initial step we computed the uv spectrum of gadusol in the gas phase , and the results are summarized in table  1 .
experimental and franck  condon vertical caspt2 excitation energies , orbital transitions , and oscillator strengths ( f ) for gadusol and gadusolate in the gas phase .
it follows from these results that the computed caspt2 uv spectrum for gadusol includes two n * transitions to s1 and s3 with low oscillator strength values , as is usual for this type of transition .
the relevant band corresponds to a  * transition located at 237  nm with an oscillator strength of 0.35 .
thus , the bright state in gadusol is s2 , and this state will be mainly populated after radiation absorption . the experimental band maximum was reported to be 268  nm in water .
therefore , to determine the relevance of the solvent in our calculations , we computed the uv spectrum using the polarizable continuum model ( pcm ) with water and acetonitrile ( table  2 ) . in both solvents ,
the state order was the same as that found in the gas phase , and the bright state was again s2 with oscillator strength values of 0.25 ( water ) and 0.26 ( acetonitrile ) , both similar to the value of 0.35 obtained in the gas phase .
the computed maximum in water was located at 242  nm , whereas in acetonitrile the maximum appeared at 239  nm .
the effect of the basis set was also explored by recomputing the absorption spectra in solvent with the ano - l - vdz basis set . in this case , the absorption in water was found at 245  nm , whereas in acetonitrile it appeared at 241  nm . as can be seen from the results in tables  1 and 2 , the caspt2 values with 6 - 31 g * in gas phase are qualitatively good enough to reproduce the experimental data .
next , we explored the excited - state deactivation of gadusol after light absorption and population of the bright - state s2 .
the minimum energy path ( mep ) computed for the relaxation along the s2 potential energy surface ( pes ) is reported in figure  1 .
condon region , very fast deactivation of s2 is expected , as no minima or energy barriers were found before crossing between s2 and s1 at a ci point .
once the ci s2/s1 point is reached , decay to s1 can be easily achieved .
the main molecular distortion after relaxation in s2 implies an out - of - plane movement of the oxygen atoms in the c = o and oh moieties of the chromophore .
comparison between the structures of gadusol in the min s0 and ci s2/s1 structures shows that the molecule remains substantially planar in the ground state to minimize the energy through conjugation , whereas at the ci point the molecule is not planar due to partial breakage of the  system .
the c = o bond has an initial length in the ground state of 1.23   , while the length in the ci s2/s1 point is 1.34   .
also , the c = c bond length changes from 1.35 to 1.49   .
in addition , the c = o and oh moieties at positions  1 and 3 of the chromophore move out of the plane of the molecule in the ci s2/s1 point ( 64 and 51 , respectively ) in agreement with the  * nature of the excited state . similarly , relaxation along the s1 pes leads directly to a different ci point connecting with the ground state , as shown in figure  2 .
thus , ultrafast deactivation takes place in s1 as well , as was predicted for s2 .
again , the main molecular distortions are related to breakage of the  system together with the out - of - plane movement of the oxygen atoms in the carbonyl and hydroxy groups that are part of the chromophore .
this implies that the same molecular movements are involved in the deactivation in both s2 and s1 . in turn
, this allows an ultrafast deactivation from the photon absorption to the ground - state recovery .
the initial c = o bond subsequently lengthens along the s1 pes and at the ci s1/s0 point has a value of 1.58   .
correspondingly , the initial c = c bond reaches a value of 1.51   .
the substituents at positions  1 and 3 are displaced from the plane of the molecule by 34 ( c = o moiety ) and 65 ( oh at position  3 ) .
the nuclear motion that lifts the degeneracy between the two states is represented by the gradient difference ( gd ) and derivative coupling ( dc ) vectors shown in figure  3 .
a )  derivative coupling ( dc ) and b )  gradient difference ( gd ) vectors for gadusol ci s1/s0 .
the nuclear motion pointed by the direction of the vectors suggests gadusol recovery once the ground state is reached . to check this , we optimized four different geometries in the ground state obtained from a slight distortion of the ci s1/s0 geometry in the directions marked by the dc and gd vectors . in all cases ,
the same geometry of the min s0 was found , probing the high photostability of gadusol .
this is consistent with the low photodecomposition quantum yield experimentally determined in solution.14 this , together with the occurrence of two energetically accessible conical intersection points ( ci ) , provides gadusol with excellent features as a photoprotective compound .
these conical intersection points connect states s2/s1 and s1/s0 that can be reached by small geometrical changes , thus enabling the ultrafast decay from the franck 
energies are in kcal  mol , relative to the ground - state minimum . as stated above , the photophysical properties of gadusol were found to be ph dependent.14 to explain this behavior , we aimed at studying the enolate forms that could be present in neutral solutions . from the various alternatives that could be found in a neutral medium
these two structures correspond to deprotonation at two sites , both of which are considered to be capable of forming relatively stable anions .
each of the structures in figure  5 was computed in the gas phase and in water using the pcm method . in both cases structure
a was found to be more stable than b. this is not surprising , as the enolate a can delocalize the negative charge along the  system , whereas the charge in enolate b is localized at the oxygen atom .
this causes structure a to have similar co ( 1.23   in both cases ) and cc ( 1.42 and 1.40   ) distances .
however , structure b features clearly distinct co ( 1.21 vs. 1.34   ) and cc ( 1.47 and 1.33   ) bond lengths .
in addition , the relative energies of these two species show that no mixture of compounds should be present in neutral solutions of gadusol , as the energy difference is , in any case , very high .
results from cas(16,10)/6 - 31 g * with pcm indicate that enolate a is 42.5  kcal  mol more stable than b. thus , in the study of gadusol in aqueous solution , we performed our calculations exclusively on enolate a. the gas - phase absorption spectrum was computed , and the results are summarized in table  1 . a strong absorption ( f=0.97 ) centered at 248  nm and characterized by a ( , * ) transition was found for s1 , while s2 shows ( n, * ) character with a very small oscillator strength .
this is consistent with the experimental band maximum found for gadusol at ph  7.14 the inclusion of water in the calculations by using pcm does not affect the qualitative picture , and the bright state is also s1 , with an oscillator strength value of 0.91 and a band centered at 263  nm .
the use of the ano - l - vdz basis set yields a band at 273  nm . comparison between gadusol and gadusolate absorption spectra denotes some differences .
first , in the gas phase the band maximum for gadusolate appears at 248  nm , whereas for gadusol the maximum is located at 237  nm . also , gadusolate has a much more intense absorption ( f=0.97 vs. 0.38 ) .
thus , gadusolate , the main species present in neutral aqueous solutions under pseudophysiological conditions , shows both a higher and red - shifted absorption relative to gadusol .
these two characteristics suggest an increased efficiency of gadusolate in terms of the photoprotective role in a natural environment .
also , as s1 is the bright state , gadusolate could have a far simpler deactivation pathway to dissipate the light energy , thus implying an improved performance .
so far this observation is in agreement with the experimentally observed lower photodecomposition quantum yield for gadusolate than for gadusol.14 to determine the photoprotective capacity of gadusolate , we computed the reaction pathway starting from the franck 
energies are relative to ground - state minimum . after light absorption , relaxation of gadusolate along
interestingly , the path connecting the franck  condon region with the ci point implies an ultrafast deactivation of the excited gadusolate , as no energy barrier is located along the reaction path .
the main geometry deformation is the out - of - plane motion of the negatively charged oxygen atom of the chromophore ( 71 ) together with the c = c bond elongation ( from 1.40   in min s0 to 1.45   at the ci point ) .
in contrast , the c = o bond length remains unaltered ( 1.23   in both geometries ) with a small out - of - plane displacement of 45. this deformation is similar to the geometry found for the ci points in related compounds , such as gadusol or palythine.22 the branching plane for this ci is characterized by the dc and gd vectors shown in figure  7 . as expected , the nuclear motion that lifts the degeneracy between the states corresponds to the arrangement of the nuclei that make part of the chromophore , keeping the rest of the molecule almost unaltered .
thus , it seems that both the excitation and geometry deformation are mainly located in just a part of the molecule . from a photophysical point of view
, the photostability of gadusol is mainly due to the substituted cyclohexenone moiety , and the rest of the molecule is not relevant for the relaxation pathways of the excited state .
a )  derivative coupling ( dc ) and b )  gradient difference ( gd ) vectors for gadusolate ci s1/s0 .
deformation of the ci point geometry in the directions indicated by the two vectors leads to direct recovery of the starting material .
thus , no side - products are expected after irradiation of gadusolate which may account for the photodecomposition of this compound .
as shown by the theoretical calculations , although the deprotonated form of gadusol is responsible for the main absorption at ph  7 and thus it is the relevant species under physiological conditions , both gadusol and gadusolate share some photochemical features that confer these compounds a remarkable uv - protective capacity .
also , the absence of any intermediates along the reaction paths is in agreement with the lack of fluorescence observed for these compounds .
this fact , together with the recovery of the starting material after irradiation supported by the low decomposition quantum yields , imparts high photostability to these molecules , which in turn allows for efficient energy wastage through subsequent photocycles .
the photoprotective mechanism of gadusol and gadusolate is very similar to that previously described for maas.9 ,  10 ,  12 ,  14 ,  22 this also suggests the relation between gadusol and maas in terms of metabolic routes or even as different steps in the evolution of natural photoprotective compounds .
the structural and mechanistic information provided by the molecular - level study of the photoprotective capacities of these and related compounds could encourage the design of new and more efficient products with interesting properties as sunscreens for commercial applications .
computational details : the most stable ground - state conformation was selected as the starting structure for the complete photochemical study , performed at multi - state complete active space perturbation to second order with a complete active space self - consistent field reference wavefunction ( i.e. , ms - caspt2//sa - casscf methodology).27 ,  28 this methodology has already shown success in describing the photochemistry of various molecular systems.29 ,  30 for this study , the state averaged - complete active space self - consistent field ( sa - casscf ) method27 was used for calculation of the electronic transitions and the computed minimum energy paths ( meps ) , including four states ( s0 , s1 , s2 , and s3 ) with the same weight in the averaged wavefunction .
the active space chosen comprises all the  and  * orbitals together with the n orbitals of the oxygen atoms that influence the chromophore ( 16 electrons in 10 orbitals ) .
both casscf and caspt2 calculations were performed with the standard 6 - 31 g * basis set .
the uv spectra were also computed by using the ano - l - vdz basis set .
meps were computed at the casscf level with the methodology present in molcas-6.4.31 meps representing steepest descendent minimum - energy reaction paths were built through a series of geometry optimizations , each requiring minimization of the potential energy on a hyperspherical cross - section of the potential energy surface ( pes ) centered on a given reference geometry and characterized by a predefined radius .
the paths were computed by taking a value of 0.1  amu for the steps .
starting from the franck  condon ( fc ) structure , the path was followed to a conical intersection point . once each new structure was obtained ,
this was taken as the new hypersphere center , and the procedure was iterated until the bottom of the energy surface was reached .
bulk solvent effects on the uv spectra were included by using the polarizable continuum model ( pcm)32 as implemented in molcas-6.4 .
the molecule is considered as included in a cavity surrounded by an infinite medium with the dielectric constant corresponding to the specific solvent .
the standard values of 78.4 for water and 38.8 for acetonitrile were considered in these calculations .
the uv spectra were computed under non - equilibrium conditions , that is , only solvent electronic polarization is in equilibrium with excited - state electron density .
this kind of calculation is more adequate to compute vertical excitation energies , as those needed for the uv spectra .
calculations of the gradient difference ( gd ) and derivative coupling ( dc ) vectors were performed with the gaussian 03 software package,33 and molcas 6.4 was used for calculation of the meps and caspt2 single - point energy corrections .
as a service to our authors and readers , this journal provides supporting information supplied by the authors .
such materials are peer reviewed and may be re - organized for online delivery , but are not copy - edited or typeset . technical support issues arising from supporting information ( other than missing files ) should be addressed to the authors .
 gadusol shows one of the simplest structures among a series of natural uv - absorbing compounds that have been related to the photoprotective and antioxidant functions in aquatic organisms . 
 caspt2//casscf methodology was used to carry out a theoretical study on this basic structure in order to describe the underlying features responsible for the photoprotective capacity of the molecule . 
 the influence of the enol  enolate equilibrium on the photophysical properties was explored . 
 the results confirm that both forms undergo a rapid deactivation , which very efficiently dissipates light energy as heat . 
 this work highlights the potential of molecular - level studies to provide an understanding of natural photoprotective mechanisms and gives support to the future design of structurally related new synthetic sunscreens . 
Introduction
Results and Discussion
Conclusions
Experimental Section
Supporting Information