5,6-Dihydroxyindole Oxidation in Phosphate Buffer/Polyvinyl Alcohol: A New Model System for Studies of Visible Chromophore Development in Synthetic Eumelanin Polymers

The determinants of the broadband absorption spectrum of eumelanins are still largely unknown. Herein we report a novel approach to investigate eumelanin chromophore which is based on the biomimetic oxidation of the key monomer precursor, 5,6-dihydroxyindole (DHI, 1 ), with peroxidase/hydrogen peroxide in phosphate buffer, pH 7, containing 1–5% polyvinylalcohol (PVA, 27 000 Da). This approach relies on the discovery that as low as 1% PVA can prevent precipitation of the growing melanin polymer thus allowing investigation of the chromophoric phases accompanying oxidation of DHI without confounding scattering effects. Spectrophotometric monitoring showed the initial development of a band around 530 nm persisting for about 1 h before gradually changing into the typical broadband spectrum of eumelanin. Reductive treatment caused a significant absorbance decrease in the visible region without affecting an absorption band around 320 nm. Initial product analysis indicated an altered formation ratio of 2,4′-biindolyl ( 2 ) and 2,7′-biindolyl ( 3 ) relative to control experiments. Overall, these results demonstrate for the first time that the development in solution of visible chromophores since the early oligomer stages is independent of strong aggregation/precipitation phenomena.     

Intermolecular π‐Electron Perturbations Generate Extrinsic Visible Contributions to Eumelanin Black Chromophore in Model Polymers with Interrupted Interring Conjugation

The key structural factors underlying the unique black chromophore of eumelanin biopolymers have so far defied elucidation. Capitalizing on the ability of 1% polyvinylalcohol (PVA) to prevent pigment precipitation during melanogenesis in vitro, we have investigated the visible chromophore properties of soluble eumelanin‐like polymers produced by biomimetic oxidation of 5,6‐dihydroxyindole (DHI) and 5,6‐dihydroxyindole‐2‐carboxylic (DHICA) in 1% PVA‐containing buffer at pH 7. Upon dilution DHI‐eumelanin solutions exhibited almost linear visible absorbance changes, whereas DHICA‐eumelanin displayed a remarkable deviation from linearity in simple buffer, but not in PVA‐containing buffer. It is suggested that in DHICA polymers, exhibiting repeated interruptions of interring conjugation due to lack of planar conformations, the black chromophore is not due to an overlap of static entities defined intrinsically by the conjugation length across the carbon frame, but results largely from aggregation‐related intermolecular perturbations of the π‐electron systems which are extrinsic in character.     

Effect of stacking and redox state on optical absorption spectra of melanins -- comparison of theoretical and experimental results

In this work the effect of aggregation and oxidation on the optical absorption of eumelanin oligomeric sheets is investigated by applying quantum mechanics and atomistic simulation studies to a simplified eumelanin structural model that includes 1-3 sheets of hexameric oligomer sheets. The oligomeric hypothesis is supported by AFM characterizations of synthetic eumelanins, formed by auto-oxidation or electrochemical oxidation of dihydroxyindole (DHI). Comparison of calculated absorption spectra to experimental spectra demonstrates a red shift in absorption with oxidation and stacking of the eumelanin and validates the theoretical results.     

Theoretical study on the excited-state intramolecular proton transfer in the aromatic schiff base salicylidene methylamine: an electronic structure and quantum dynamical approach

The proton-transfer dynamics in the aromatic Schiff base salicylidene methylamine has been theoretically analyzed in the ground and first singlet (pi,pi) excited electronic states by density functional theory calculations and quantum wave-packet dynamics. The potential energies obtained through electronic calculations that use the time-dependent density functional theory formalism, which predict a barrierless excited-state intramolecular proton transfer, are fitted to a reduced three-dimensional potential energy surface. The time evolution in this surface is solved by means of the multiconfiguration time-dependent Hartree algorithm applied to solve the time-dependent Schr?dinger equation. It is shown that the excited-state proton transfer occurs within 11 fs for hydrogen and 25 fs for deuterium, so that a large kinetic isotope effect is predicted. These results are compared to those of the only previous theoretical work published on this system [Zgierski, M. Z.; Grabowska, A. J. Chem. Phys. 2000, 113, 7845], reporting a configuration interaction singles barrier of 1.6 kcal mol(-1) and time reactions of 30 and 115 fs for the hydrogen and deuterium transfers, respectively, evaluated with the semiclassical instanton approach.     

The ring-opening reaction of chromenes: a photochemical mode-dependent transformation

The efficiency of the photochemical ring-opening of chromenes (or benzopyrans) depends on the vibronic transition selected by the chosen excitation wavelength. In the present work, ab initio CASPT2//CASSCF calculations are used to determine the excited-state ring-opening reaction coordinate for 2H-chromene (C) and 2,2-diethyl-2H-chromene (DEC) and provide an explanation for such an unusual mode-dependent behavior. It is shown that excited-state relaxation and decay occur via a multimodal and barrierless (or nearly barrierless) reaction coordinate. In particular, the relaxation out of the Franck-Condon involves a combination of in-plane skeletal stretching and out-of-plane modes, while the second part of the reaction coordinate is dominated exclusively by a different out-of-plane mode. Population of this last mode is shown to be preparatory with respect to both C-O bond breaking and decay via an S(1)/S(0) conical intersection. The observed mode-dependent ring-opening efficiency is explained by showing that the vibrational mode corresponding to the most efficient vibronic transition has the largest projection onto the out-of-plane mode of the reaction coordinate. To support the computationally derived mechanism, we provide experimental evidence that the photochemical ring-opening reaction of 2,2-dimethyl-7,8-benzo(2H)chromene, that similarly to DEC exhibits a mode-dependent photoreaction, has a low ( approximately 1 kcal mol(-1)) activation energy barrier.     

Exploring the frontiers of synthetic eumelanin polymers by high‐resolution matrix‐assisted laser/desorption ionization mass spectrometry

New trends in material science and nanotechnologies have spurred growing interest in eumelanins black insoluble biopolymers derived by tyrosinase‐catalysed oxidation of tyrosine via 5,6‐dihydroxyindole (DHI) and its 2‐carboxylic acid (DHICA). Efficient antioxidant and photoprotective actions, associated with peculiar optoelectronic properties, are recognised as prominent functions of eumelanin macromolecules within the human and mammalian pigmentary system, making them unique candidates for the realisation of innovative bio‐inspired functional soft materials, with structure‐based physical–chemical properties. An unprecedented breakthrough into the mechanism of synthetic eumelanin buildup has derived from a detailed investigation of the oxidative polymerization of DHI and its N‐methyl derivative (NMDHI) by linear and reflectron matrix‐assisted laser/desorption ionization mass spectrometry. Regular collections of oligomers of increasing masses, spanning the entire m/z ranges up to 5000 Da (>30‐mer) and 8000 Da (> 50‐mer) for the two building blocks, respectively, were disclosed. It is the first time that the in vitro polymerisation of dihydroxyindoles to form synthetic eumelanins is explored up to its high mass limits, giving at the same time information on the polymerisation mode, whether it follows a stepwise pattern (being this the conclusion in our case) or a staking sequencing of small‐sized entities. It also highlighted the influence of the N‐methyl substituent on the polymerization process; this opens the way to the production of N‐functionalized, synthetic eumelanin‐inspired soft materials, for possible future technological applications. Copyright © 2012 John Wiley & Sons, Ltd.     

Exciton interactions in helical crystals of a hydrogen-bonded eumelanin monomer

Eumelanin, a naturally occurring group of heterogeneous polymers/aggregates providing photoprotection to living organisms, consist of 5,6-dihydroxyindole (DHI) and 5,6-dihydroxyindole-2-carboxylic acid (DHICA) building blocks. Despite their prevalence in the animal world, the structure and therefore the mechanism behind the photoprotective broadband absorption and non-radiative decay of eumelanin remain largely unknown. As a small step towards solving the incessant mystery, DHI is crystallized in a non-protic solvent environment to obtain DHI crystals having a helical packing motif. The present approach reflects the solitary directional effect of hydrogen bonds between the DHI chromophores for generating the crystalline assembly and filters out any involvement of the surrounding solvent environment. The DHI single crystals having an atypical chiral packing motif (P2₁2₁2₁ Sohncke space group) incorporate enantiomeric zig-zag helical stacks arranged in a herringbone fashion with respect to each other. Each of the zig-zag helical stacks originates from a bifurcated hydrogen bonding interaction between the hydroxyl substituents in adjacent DHI chromophores which act as the backbone structure for the helical assembly. Fragment-based excited state analysis performed on the DHI crystalline assembly demonstrates exciton delocalization along the DHI units that connect each enantiomeric helical stack while, within each stack, the excitons remain localized. Fascinatingly, over the time evolution for generation of single-crystals of the DHI-monomer, mesoscopic double-helical crystals are formed, possibly attributed to the presence of covalently connected DHI trimers in chloroform solution. The oligomeric DHI (in line with the chemical disorder model) along with the characteristic crystalline packing observed for DHI provides insights into the broadband absorption feature exhibited by the chromophore.

Single crystals of DHI monomer, a eumelanin precursor, adopt an atypical chiral packing arrangement incorporating enantiomeric zig-zag helical stacks while its covalently connected DHI trimer forms double-helical crystals in the mesoscopic scale.

Eumelanin, which represents a broad class of natural pigments found in the animal kingdom, acts as a biological shield for protecting the skin cells against harsh UV radiation.¹ Eumelanin, a black coloured pigment, obtained from the oxidative polymerization of 5,6-dihydroxyindoles (DHIs) and 5,6-dihydroxyindole carboxylic acid (DHICA) is one of the extensively explored archetypes of the melanin family (Fig. 1a).² Eumelanin has synergistic merits of possessing broadband UV absorption and proficient dissipation of the excessive electronic energy via non-radiative deactivation of the excited states, thereby resulting in the photoprotective nature of the pigment.³ Apart from the photoprotective behaviour, melanin possesses exceptional antioxidant activity via its free radical scavenging traits.^4–7 However, much less has been understood about the fundamental photophysics and structural features of eumelanin due to the enormous heterogeneity in the molecular framework^8,9 coupled with poor solubility in common solvents.¹⁰ Recent years have witnessed a growing interest towards unravelling the excited state processes occurring in the eumelanin pigment upon interaction with light.^11,12 A better correlation between the structure–property relationship and photoexcited state processes in eumelanin can guide the development of inspired functional materials for potential application in biomedical and dermo-cosmetic fields.^1,13–15Open in a separate windowFig. 1(a) Chemical diagrams of 5,6-dihydroxyindole (DHI) and 5,6-dihydroxyindole carboxylic acid (DHICA). (b) Various unconventional hydrogen bonding interactions identified in DHI crystals.In the natural world, the chromophoric architecture of eumelanin can symbolize an organized and efficient organic system for photoprotection that nature developed through evolution. The outcome of the research done so far indicates the presence of continuous π-stacks of oligomers in eumelanin which induce different levels of aggregation to construct the eumelanin framework.¹⁶ Furthermore, eumelanin has been reported to exhibit weak fluorescence which indicates the presence of competing non-radiative channels that provide efficient de-excitation pathways for repopulating the ground states.^17–20 The broadband absorption of eumelanin has theoretically and experimentally been evaluated, in part, to be a result of π-stacking interactions between the oligomers of DHI/DHICA in multiple oxidation states.^21–24 While dealing with biomacromolecules, non-covalent interactions such as hydrogen bonding and π–π stacking take the centre stage in controlling the supramolecular architecture, especially in the 3-dimensional structures of proteins and DNA. Modulating the balance between each of these noncovalent interactions over another will produce significant changes not only in the structure but also in the functional properties.Hydrogen bonding is simultaneously both ubiquitous and diverse and therefore its significance in biochemical systems comes as no surprise particularly due to the surrounding aqueous environment. Apart from the classical hydrogen bonding interactions, an array of hydrogen bond-like weak interactions which include a delocalized π-system acting as the acceptor group to the X–H hydrogen donor (X = O, N, C) is identified to provide additional contributions for stabilizing the biomolecular structure and controlling intrinsic functions (Fig. 1b).^25,26 Investigations aimed at identification (using X-ray crystallography) followed by energetic quantification of the major stabilizing interactions such as those with the aromatic π-rings in biological complexes are of paramount importance for developments in diverse areas including drug design. There has been extensive research conducted on eumelanin building blocks showcasing their ability to form hydrogen bonds through the –OH and –NH functional groups.^27,28 Most reports almost exclusively focus on the hydrogen bonding with the solvent environment surrounding the eumelanin monomer units.²⁹ Findings from these theoretical studies have demonstrated the role of several deactivation pathways in the presence of a protic solvent, namely –OH and –NH bond elongation and 5-/6-membered ring puckerings.¹⁸Chemical and spectral evidence from the eumelanin polymeric structures identified so far points to five main levels of chemical disorder leading to the supramolecular structure, which includes (i) disorder from the simultaneous presence of different building blocks; (ii) molecular size disorder; (iii) disorder from the position of coupling; (iv) electronic/redox disorder of the constituent units and (v) supramolecular disorder.³⁰ Given the complex structure of melanin, a bottom-up approach using the building blocks or basic constituent molecules of eumelanin is a pertinent strategy for the mechanistic study of the photoproperties of eumelanin. This can be followed by understanding the more intricate structures of melanin formed from the constituents with less complex approaches. Due to the abundant presence of water in the natural media, the corresponding solute–solvent interactions can have profound significance in driving the fast polymerization and the consequent heterogeneity of natural melanin. The tedious task of extracting melanin from natural sources and the lack of solubility of the polymeric melanin material in organic/aqueous solvents have called for basic model systems to understand the complex eumelanin architecture. In this regard, we have adopted a facile approach to decode the perpetual puzzle by single crystal X-ray crystallographic and spectroscopic analyses of DHI crystalline aggregates derived from a non-aqueous environment (chloroform). Due to the highly autooxidative³¹ nature of eumelanin precursors even in the slight presence of protic solvents, the simple model implemented here precludes the contribution of solute–solvent hydrogen bonding interaction towards the formation and resultant structure of DHI crystalline aggregates. In chloroform, each DHI molecule experiences weak interactions solely from the neighbouring DHI chromophores thereby leading to helical aggregation.Our efforts towards recognizing and monitoring the photogenerated excitons and charge-transfer dynamics in crystalline and contorted polyaromatic assemblies^32–36 prompted us to explore the structure–optical property relationship in the eumelanin precursor molecule DHI. Unlike the commonly understood π-stacking in eumelanin derivatives, the single crystals of DHI arrange in a helical zig-zag fashion with a completely edge-to-face aggregate structure driven by both conventional and unconventional hydrogen bonds (Fig. 1b). The structural heterogeneity imposed by the different hydrogen bonds has led to varied levels of exciton delocalization between the neighbouring chromophores in the crystalline DHI aggregates. Along with the diffracting single crystals of monomeric DHI, covalently connected trimeric units of DHI are also identified in chloroform solution, which form double-helical crystals in the mesoscopic scale. Such double-helical architectures are omnipresent in nature as exemplified by the DNA structure.DHI was synthesized by following a previously reported procedure having l-dopa as the starting material (Scheme S1, Appendix C1–C4 and C7, ESI^†).³⁷ Slow evaporation from dry chloroform solution of DHI produced colourless diamond shaped single crystals of DHI (Fig. S1a^†). Interestingly, the DHI molecule with no chiral centre atypically crystallized in the Sohncke space group, P2₁2₁2₁ (Table S1^†). Single-crystal X-ray structure analysis revealed the presence of conventional and unconventional hydrogen bonds directing the crystalline self-assembly of DHI chromophores about a zig-zag helical backbone. Fig. 2 presents the four different types of non-covalent dimers (D1–4) distinguished within the DHI crystal. The zig-zag helical stacks proceed along the crystallographic a-axis (Fig. S2^†) and are fabricated by bifurcated O–H⋯O hydrogen bonds (d_O1⋯H = 2.25 Å, d_O2⋯H = 2.34 Å, <O1–H–O2 = 67.89°) between the two hydroxyl substituents in the DHI chromophore (D3 in Fig. 2). Such bifurcated hydrogen bonded assemblies are prevalent in the secondary and tertiary structures of proteins.³⁸ Interestingly, in the DHI crystal, enantiomeric helical stacks (Fig. S2^†) that are arranged with respect to different screw axes are observed, wherein each stack aligns in a herringbone fashion to the other zig-zag helix (as represented by the dimers D1, D2 and D4). The stacks are interconnected majorly through the unconventional hydrogen bonds such as C–H⋯π (d_C–H⋯π = 2.66–2.99 Å), O–H⋯π (d_O–H⋯π = 2.59 Å), C–H⋯N (d_C–H⋯N = 2.72 Å), C–H⋯O (d_C–H⋯O = 2.66 Å) and the classical N–H⋯O (d_N–H⋯O = 2.72 Å) hydrogen bonds. The absence of π–π stacking interaction is validated by the Hirshfeld surface analysis wherein the formation of the DHI crystalline assembly is majorly stabilized by the C⋯H (40.5%), H⋯H (29.7%), N⋯H (4.3%) and O⋯H (25.4%) noncovalent interactions (Fig. S3 and Table S2^†).Open in a separate windowFig. 2Different orientations of DHI (D1–4) and the directing hydrogen bonds observed in the single crystal.Detailed examination of the interchromophoric interactions supporting the zig-zag helical stacks in DHI crystals using Bader''s quantum theory of atoms in molecules (QTAIM) analysis revealed the presence of supramolecular synthons in the crystal system (Fig. S4^†). This is evidenced by the (3,+1) ring critical points in each of the representative dimers. For a molecular self-assembly to occur efficiently, recognition between the intermolecular functionalities is important, which often culminates in the formation of smaller repetitive units or supramolecular synthons.³⁹ The recognition information which is then carried by these units forms the kernel of self-aggregation or crystallization processes. In the case of DHI crystals, all the dimer assemblies D1–D4 display synthon formation with the D1 and D3 synthons showing greater energetic stabilities. The dimer unit representing the helical backbone, D3, forms two supramolecular synthons orchestrated by the bifurcated O–H⋯O hydrogen bonds and a weak C–H⋯O interaction (Fig. S4^†). The five- and six-membered rings so formed fabricate the helical zig-zag backbone of the DHI crystal. Similarly, the dimer D1 also forms two synthons through a classical N–H⋯O interaction along with the weak C–H⋯π and O–H⋯π interactions. D2 and D4 dimers are stabilized by one synthon each, materialized by C–H⋯π and C–H⋯N interactions in D2 and N–H⋯O and C–H⋯O interactions in D4.The synthon formation and the concomitant electron delocalization involving the π-rings in dimers D1 and D2 have resulted in the aromatic stabilization of the π-rings upon comparison with the monomer DHI. The nucleus independent chemical shift (NICS(1)) values evaluated for D1–4 and monomer DHI in the ground state indicate the aromatic stabilization of the π-surface leading to the favourable alternate stacking of the enantiomeric zig-zag helices facilitated by the unconventional hydrogen bonds. The negative NICS(1) values for the six- and five-membered (6C, 5C) rings of molecule A (Fig. S5 and Table S3^†) increased to −27.22 ppm (6C) and −29.13 ppm (5C) when compared to the monomer DHI (6C: −26.14 ppm, 5C: −28.03 ppm). Similarly, in D2, the π-surface of molecule B that is involved in the weak interaction undergoes aromatic stabilization (6C: −29.43 ppm, 5C: −31.48 ppm). Truncated symmetry adapted perturbation theory (SAPT(0)) analysis⁴⁰ performed on the DHI dimers shows higher stabilization for D1 (E^SAPT_int = −9.70 kcal mol⁻¹) and D3 (E^SAPT_int = −6.57 kcal mol⁻¹) dimers, which could be attributed to the two supporting supramolecular synthons in both the dimers (Table S4^†). The total stabilization of D1 and D3 orientations is facilitated by the higher contributions of electrostatic energy (E⁽¹⁾_elc = −6.17 to −6.14 kcal mol⁻¹) and induction energy (E⁽¹⁾_elc = −1.60 to −1.38 kcal mol⁻¹) towards the total SAPT energy. The dominant role of the stronger classical hydrogen bonds in the fabrication of D1 and D3 synthons when compared to the other DHI orientations (having equal contribution from the weak unconventional interactions) explains the observed energy distribution in the SAPT(0) analysis.The crystalline architecture of the DHI precursor molecule identified herein could provide a sound model for understanding the inherent nature of the excited energy states leading to the characteristic photo-function of the eumelanin pigment. Several experimental and theoretical investigations on eumelanin aggregates revealed the occurrence of excitation energy transfer within the aggregates.^41,42 The extent of energy delocalization within the four dimer orientations in the DHI crystal structure is determined using the fragment-based excited state analysis developed by Plasser, executed in the TheoDORE package.^43–45 The expanse of excitation delocalization amongst the fragment units is described by the value of participation ratio (PR). The contribution towards the exciton delocalization from the fragments is defined by the mean position or the POS value, wherein the indicated number shows the involvement of one or more units. The charge-transfer or Frenkel character of the excited states is defined by the CT number which assumes values closer to one for pure CT states and closer to zero for pure Frenkel states. In dimers D1 and D2, the first singlet excited state of highest oscillator strength (S2) shows Frenkel exciton character with effective delocalization of the excitons between the individual constituting units during the excitation process (Tables S5 and S6^†). The low CT values of the S2 states in D1 (PR = 1.79, POS = 1.33, CT = 0.08) and D2 (PR = 1.98, POS = 1.45, CT = 0.02) along with a PR close to a value of two indicated the delocalization of the Frenkel excitons on the two monomers. Hence, there exists a possible excitation energy delocalization along the adjacent enantiomeric stacks in the DHI crystal. However, for D3 with the bifurcated hydrogen bonding interaction between the monomeric units, both S1 and S2 states have significant oscillator strength. The Frenkel excitons in S1 and S2 states remain localized on only one fragment of D3 (Table S7^†), while in D4, there exists a partial delocalization of the Frenkel excitons in the S1 state (Table S8^†). Thus, it is understood that within each enantiomeric stack the initial Frenkel excitons remain localized on one fragment. The natural transition orbitals (NTOs) of dimers D1–4 give an idea about the nature of the dominant orbital transitions for the allowed electronic excitations. The allowed electronic excitation in all the dimers could only be well represented by taking two distinct orbital transitions with significant coefficients into consideration. The absence of a single dominant orbital transition in the dimers proposed the need for fragment-based hole–electron analysis for the better understanding of the electronic energy delocalization. The hole–electron isosurface analysis provides a pictorial representation of the corresponding delocalized and localized nature of the initial Frenkel excitons in the four DHI orientations (Fig. 3). The delocalization of the hole and electron density on both the fragments puts forward the possibility of effective energy transfer in dimers D1, D2 and D4.⁴⁶ In dimer D3, a very weak delocalization of the electron–hole densities was observed suggesting a localized exciton formation.Open in a separate windowFig. 3Hole–electron isosurface plots of the DHI orientations in the crystal. (a) D1, (b) D2, (c) D3 and (d) D4.Several factors including slight exposure to air,⁴⁷ humidity and/or light have been observed to cause the autooxidative polymerization of DHI, wherein lowered temperatures decrease the kinetics of this solid-state polymerization. Although, a major fraction of the solid oligomer mixture remains as the DHI monomer (≥80%), oligomeric units up to DHI-hexamers have been identified with varying solubilities in different solvents. MALDI-MS spectra (Appendix C5, ESI^†) sequentially collected for the oligomer mixture in CHCl₃ and DMSO indicated the major presence of DHI trimers (in CHCl₃; m/z = 442.138) and DHI-hexamers (in DMSO; m/z = 882.382) along with the smaller counterparts. Interestingly, over the time course for nucleation and the subsequent growth mechanism of the single crystals of monomeric DHI in chloroform, oligomerization of DHI is found to occur concurrently to form the covalently connected DHI trimer (DHI-T). Along with the diffracting single crystals of DHI, small right-handed double-helical crystals (non-diffracting, Fig. S1b–f and S6^†) are observed for the first time, which could be attributed to the self-assembled morphology of DHI-T (Fig. 4a). Observed only in chloroform, we speculate that the prolonged exposure of the chlorinated solvent plays a significant role in the chemical transformation of the DHI monomer to the covalent trimer, possibly through a radical initiated reaction.⁴⁸ The ¹H-NMR spectrum obtained for the bulk crystalline sample dissolved in CDCl₃ along with the observed MALDI-MS data of the DHI trimer evidenced the formation of DHI-T, although in very low yields when compared to the monomer (Appendix C5 and C6^†).Open in a separate windowFig. 4(a) Optical microscopy images of the right-handed double-helical crystals of the DHI trimer. (b) Normalized absorbance and excitation spectra showing DHI-T formation in chloroform solution and the Kubelka–Munk transformed diffuse reflectance spectrum of the DHI bulk crystal. (c) Optimized structure of DHI-T at the CAM-B3LYP/6-311g+(d,p) level of theory.With the understanding that DHI readily aggregates/crystallizes in chloroform, the directionality of the hydrogen bonding interactions (to the π-ring) within the definite spatial arrangements of D1 and D2 orientations hints towards the mode/position of coupling for the associated generation of DHI-T (Appendix C9^†). The molecular structure of the DHI trimer that best fits the observed characterization data (Fig. S7a and Appendix C8, ESI^†) is in line with the chemical disorder model having semiquinone and catechol units connected covalently as an extension of the D1 and D2 noncovalent interactions. Geometry optimization of the predicted structure of the DHI-T performed using the CAM-B3LYP/6-311g+(d,p) level of theory in Gaussian 16 suite shows a twisted conformation having the possibility of forming intramolecular hydrogen bonds from the –OH and –NH functionalities (Fig. 4 and S7b^†). Separation of DHI-T from DHI is a real challenge since exposure of DHI to the adsorbent in column chromatography accelerates the oxidative polymerization of DHI, resulting in a black, insoluble material difficult to characterize. Also, the presence of higher oligomer units of DHI (hexamers etc.) was not identified in the multiple data sets collected for the DHI sample dissolved in CHCl₃.Solvent-dependent steady-state UV-vis absorption and fluorescence emission measurements of DHI were performed and the line shapes of the absorption spectra of monomeric DHI in different solvents match consistently (Fig. S8^†). Two major absorption bands at λ_a1 ∼ 270 nm and λ_a2 ∼ 300 nm form the characteristic absorption spectrum of the DHI monomer. The fluorescence emission of monomeric DHI exhibits a single broad spectral feature peaking at λ_em ∼ 330 nm in a majority of the solvents. The relative fluorescence quantum yields are found to be exceptionally low in chloroform, dichloromethane, THF and water (Table S9^†) indicating the presence of non-radiative decay channels for dissipating the excitation energy. In chloroform, the emergence of a red-shifted tail in the absorption spectrum of the DHI-monomer is observed over time possibly signifying the onset of DHI-T formation. The fluorescence emission in CHCl₃ also shows a new band peaking at λ_em ∼ 460 nm along with the emission band at λ_em ∼ 335 nm (Fig. S9^†). A broad red-shifted band arising at 370 nm in the excitation spectrum of DHI solution collected in chloroform at λ_em ∼ 460 nm (Fig. 4b) evidenced the presence of the DHI trimer. A similar decrease in fluorescence quantum yield in chloroform and the concomitant emergence of new bands in the fluorescence emission profiles have been noted previously in tryptophan and other indole species synthesized for eumelanin investigation.^49,50 In such cases, the photoionization of excited indole leads to the ejection of a solvated electron which attacks the chloroform molecule, releasing a chloride ion, and further undergoes reactions to yield photoproducts.Spectroscopic investigation of the crystalline DHI (containing both monomer single crystals and covalent trimer crystals) sample showed broad, red-shifted absorption bands spanning from 210 to 560 nm (Fig. 4b and S10^†). The solid-state absorption spectrum shows two prominent bands at λ_a1 ∼ 280 nm and λ_a2 ∼ 305 nm which could be attributed to the red-shifted absorption bands of the crystalline DHI monomer (compared to the monomer absorption bands in the solution state). The observed red-shift in the absorption band of the DHI crystal arises from the nonplanar packing motif and the ensuing intermolecular interactions in the solid state. The presence of a broad shoulder band centred at λ_a3 ∼ 375 nm in the absorption spectrum could be assigned to the double-helical crystal of the covalent DHI trimer. The crystalline state fluorescence emission spectrum of DHI spans from 390 to 490 nm (Fig. S10^†). The excitation energies and the allowed vertical transitions of the monomer and DHI-T have been computed at the CAM-B3LYP/6-311g+(d,p) level of theory. Unlike the precursor DHI monomer which undergoes higher energy electronic transitions (at 270 nm and 300 nm), the favorable transition in the covalent trimer is red-shifted with the S₀ → S₁ electronic excitation occurring at ∼434 nm. Hence, the spectroscopic and theoretical investigation of bulk crystalline DHI indicates that the broad absorption profile of the eumelanin precursor could be ascribed to the combined effects of the non-planar chromophore stacking and the presence of covalent DHI trimers that exist as double-helical aggregates. The photoprotective nature of eumelanin arises from the signature broadband absorption of eumelanin which spans throughout the UV and visible region tailing around 800 nm as explained by the chemical disorder model. In line with this understanding, the spectroscopic data of the DHI crystal also exhibit broadband absorption which expands up to 600 nm unlike the DHI monomer. Moreover, the low relative fluorescence quantum yields of DHI suggest the presence of non-radiative decay pathways within the DHI units. The oligomeric trimer which in itself shows structural heterogeneity aggregates as double helical structures and shows a red-shifted absorption band which is comparable to the computed TDDFT vertical excitation energies. Thus, our report on the characterization of DHI and the oligomeric trimer could possibly be beneficial in advancing melanin structure characterization and elucidating the photoprotective function of eumelanin.The solid-state CD spectrum of crystalline DHI (in KBr, Fig. S11a^†) showed the signatures for the presence of helical packing^51–55 (possibly from the zig-zag helical motif along with the double helical arrangement). However, the basis of the CD couplet of significant intensity spanning from ∼250 to 600 nm (including ranges outside of the absorption maxima) could not be exclusively assigned to chiral absorption from the chromophoric packing.⁵⁶ In the case of macromolecular systems having long-range organization, differential scattering of incident left and right circularly polarized light can provide significant contributions to the observed circular dichroism.^57–59 The occurrence of broad CD bands outside the absorption bands of the macromolecule can signify the possible role of differential scattering in the circular dichroism.⁶⁰ Although, for the DHI sample, the characteristic CD spectrum has been found to be concurrent for the different data sets collected using freshly prepared crystalline samples on different days (Fig. S11a^†), the ratios of the intensities of the positive and negative bands have been observed to vary. Such a heterogeneity in the ellipticity values of the positive and negative bands could be attributed to the possible presence of different enantiomeric assemblies that exhibit varying abilities to undergo chiral absorption and differential scattering.⁵⁷ The possibility of having linear dichroism (LD) artifacts in the CD data was evaluated for the DHI sample (Fig. S11b^†).⁶¹ The LD artifact fell within the error bar of the order of 10⁻³–10⁻⁴ mdeg and hence, the contribution of LD to the strong CD signal of DHI could be ignored.^62,63 Also, the idea of having a chiral nucleation centre, probably from any conformationally chiral DHI oligomer units, leading to the double-helical aggregation and the consequent mesoscopic chirality could not be ignored while assessing the origin of the observed CD spectrum. The existing uncertainties in solving the source of the double-helical aggregation of the DHI chromophores and identifying the intermolecular forces acting behind the same remain a challenge that requires detailed examination in future studies.     

Excited-State Proton Transfer and Decay in Hydrogen-Bonded Oxazole System:MS-CASPT2//CASSCF Study

Herein we have employed high-level multi-reference CASSCF and MS-CASPT2 electronic structure methods to systematically study the photochemical mechanism of intramolecularly hydrogen-bonded 2-(2'-hydroxyphenyl)-4-methyloxazole. At the CASSCF level, we have optimized minima, conical intersections, minimum-energy reaction paths relevant to the excited-state intramolecular proton transfer (ESIPT), rotation, photoisomerization, and the excited-state deactivation pathways. The energies of all structures and paths are refined by the MS-CASPT2 method. On the basis of the present results, we found that the ESIPT process in a conformer with the OH…N hydrogen bond is essentially barrierless process; whereas, the ESIPT process is inhibited in the other conformer with the OH…O hydrogen bond. The central single-bond rotation of the S₁ enol species is energetically unfavorable due to a large barrier. In addition, the excited-state deactivation of the S₁ keto species, as a result of the ultrafast ESIPT, is very efficient because of the existence of two easily-approached keto S₁/S₀ conical intersections. In stark contrast to the S₁ keto species, the decay of the S₁ enol species is almostly blocked. The present theoretical study contributes valuable knowledge to the understanding of photochemistry of similar intramolecularly hydrogen-bonded molecular and biological systems.     

Mechanistic Insights into the Photophysics of Ortho-hydroxyl GFP Core Chromophores

Herein we have employed the MS-CASPT2//CASSCF method to study the S₁ excited-state intramolecular proton transfers (ESIPTs) of recently synthesized ortho-hydroxyl GFP core chromophores, i.e. OHIM, CHBDI, and MHBID, and their excited-state relaxation pathways. We have found that in OHIM and CHBDI, the ESIPT process is associated with small barriers of 3.4 and 4.2 kcal/mol; while, in MHBDI, it becomes essentially barrierless. Moreover, we have found two main S₁ excited-state radiationless channels. In the first one, the enol S₁ species decays to the S₀ state via the enol S₁/S₀ conical intersection after overcoming considerable barriers of 7.0 and 7.7 kcal/mol in OHIM and CHBDI (however, in MHBDI, it is nearly barrierless). In the second one, the keto S₁ species is first generated through the ESIPT event; then, it is de-excited into the S₀ state in the vicinity of the keto S₁/S₀ conical intersection. These energetically allowed excited-state decay channels rationalize experimentally observed ultralow fluorescence quantum yields. The insights gained from the present work may help to guide the design of new ortho-hydroxyl GFP core chromophores with improved fluorescence emission and brightness.     

Excited-state hydrogen relay along a blended-alcohol chain as a model system of a proton wire: deuterium effect on the reaction dynamics

The excited-state deuteron transfer (ESDT) of deuterated 7-hydroxyquinoline (7DQ) along a heterogeneous hydrogen (H)-bonded chain composed of two deuterated alcohol (ROD) molecules having different acidities, as a model system of a proton wire consisting of diverse amino acids, has been investigated. To understand dynamic differences between deuteron transfer and proton transfer, solvent-inventory experiments have been performed with variation of the combination as well as the composition of alcohols in a H-bonded mixed-alcohol chain. Deuteron transfer from the adjacent ROD molecule to the basic imino group of 7DQ via tunneling, which is the rate-determining step, initiates ESDT, and subsequent barrierless deuteron relay from the acidic enolic group of 7DQ to the alkoxide moiety along the H-bonded chain completes ESDT. Whereas the acceleration of the reaction has been observed in excited-state proton transfer because of the accumulated proton-donating abilities of two alcohol molecules in a H-bonded chain by a push-ahead effect, such acceleration is not observed in ESDT. Because the energy barrier of deuteron relay is much higher than that of proton relay due to the low zero-point energy of 7DQ·(ROD)(2) and a deuteron is twice as heavy as a proton, it is hard for a deuteron to pass through the barrier via tunneling. Moreover, both the H-bonding ability and the acidity of ROD molecules are so weak that their deuteron-donating abilities cannot be accumulated at the rate-determining step of ESDT. Consequently, the rate constant of ESDT is determined mostly by the acidity of the ROD molecule H-bonded directly to the imino group of 7DQ.     

Solvochromic effects in model eumelanin compounds

We have created an indolic compound which is ideally suited to the study of the relationship between structure and function in eumelanin formation. N-methyl-5-hydroxy-6-methoxyindole (MHMI) is stable in solid and liquid states, highly soluble in a variety of solvents and forms a dimer only through the 4-4' positions. The limited binding possibilities are due to functional groups strategically placed to inhibit chemical interactions through the 2 and 7 positions. It forms a crystal structure with a remarkable packing arrangement, with four monomers grouped in parallel pairs spaced 3.5 A apart within each unit cell. Optical spectra reveal a multi-peaked absorbance profile similar to 5,6-dihydroxyindole (DHI) and N-acetyl-tryptophanamide (NATA), and strong fluorescence emission with radiative quantum yields of 29% and 33% in benzene and acetonitrile, respectively. The quantum yield is similar to that of 5,6-dihydroxyindole-2-carboxylic acid (DHICA) and shows that solvent aromaticity by itself does not affect the yield. Solution in chloroform results in an almost complete quenching of the fluorescence but an increase in emission is observed with photoactivation. Crystallographic results shown here suggest new structural possibilities for eumelanin and the controlled binding possibilities make this an excellent model for monitoring changes in function with increasing oligomer size in eumelanin formation.     

Ab initio studies on the photophysics of guanine tautomers: out-of-plane deformation and NH dissociation pathways to conical intersections

The radiationless decay mechanisms of the S1 excited states of the 7H-keto-amino, 7H-enol-amino, and 7H-keto-imino tautomers of guanine have been investigated with the CASPT2//CASSCF method. Out-of-plane deformation of the six-membered ring or the imino group as well as dissociation of NH bonds have been considered as photochemical pathways leading to conical intersections with the electronic ground state. It has been found that all three tautomers can reach S0-S1 conical intersections by out-of-plane deformation. However, only in the 7H-keto-amino tautomer the reaction path leading to the conical intersection is barrierless. This tautomer also has the lowest energy barrier for hydrogen detachment via the (1)pi sigma* state, whose potential energy surface intersects that of the (1)pi pi* state as well as that of the ground state. The other tautomers of guanine exhibit substantial energy barriers on their S1 potential energy surfaces with respect to both reaction mechanisms. These findings suggest that the 7H-keto-amino tautomer exhibits the shortest excited-state lifetime of the three tautomers due to particularly fast nonradiative deactivation processes through S0-S1 conical intersections. The computational results explain the remarkable observation that the energetically most stable 7H-keto-amino tautomer is missing in the resonant two-photon ionization spectrum of guanine in a supersonic jet. The results also explain that the energetically less stable 7H-enol-amino and 7H-keto-imino tautomers have longer excited-state lifetimes and are thus detectable by resonant two-photon ionization.     

Excited-state potential energy surface for the photophysics of adenine

The decay paths on the singlet excited-state surface of 9H-adenine and the associated energy barriers have been calculated at the CAS-PT2//CASSCF level. There are three fundamental paths for the photophysics: two paths for the (1)L(b) state which are virtually barrierless at the present level of theory and correspond to formation of the (n,pi) intermediate and direct decay to the ground state and a third path for ground-state decay of the (n,pi) state with an activation barrier of approximately 0.1 eV. The (1)L(a) state, which has the largest oscillator strength, either decays directly to the ground state or contributes indirectly to the excited-state lifetime by populating the two other states. The results are used to interpret the photophysics in terms of an excited-state plateau for the (1)L(b) state that corresponds to the short-lived excited-state component (approximately 0.1 ps) and a well (i.e., a proper minimum) for the (n,pi) state that gives rise to the long component (1 ps or more). The direct decay to the ground state of the (1)L(b) state is probably the decay channel invoked to explain the experimental wavelength dependence of the relative amplitudes of the two components. In addition to that, the excited-state component in the nanosecond range detected in the time-resolved photoelectron spectrum is proposed to be a triplet (pi,pi) state formed after intersystem crossing from the singlet (n,pi) state.     

Ultrafast Radical Photogeneration Pathways in Eumelanin†

Eumelanin is a ubiquitous biological pigment that rapidly and efficiently deactivates excited states created by UV or visible radiation. Paradoxically, photoirradiation of eumelanin also generates radicals and harmful reactive oxygen species, but the relationship between these pathways and excited-state deactivation is uncertain. Here, greatly expanding the excitation tuning range (225–620 nm) and probing window (400–1500 nm) in femtosecond transient absorption spectroscopy of the synthetic eumelanin, DOPA melanin, enables the detection of photogenerated radials with ultrafast time resolution for the first time. Despite its heterogeneous nature, the transient absorption signals can be modeled by two spectral components assigned to solvated electrons and photogenerated radicals. Radical absorbance measured several nanoseconds after excitation increases exponentially with increasing photon energy, matching the trend in radical yields measured in electron paramagnetic resonance spectroscopy experiments. Spectral modeling of the transient signals reveals two radical generation mechanisms: (1) photoionization by UV light; and (2) photoinduced charge transfer among eumelanin chromophores by UVA and visible wavelengths capable of reaching the pigment in skin. Concurrent ultrafast relaxation and radical generation underlie the ability of eumelanin to be both photoprotective and photodamaging, and the branching between these pathways likely depends on the wavelength of the absorbed light.     

Watching ultrafast barrierless excited-state isomerization of pseudocyanine in real time

The photoinduced excited-state processes in 1,1'-diethyl-2,2'-cyanine iodine are investigated using femtosecond time-resolved pump-probe spectroscopy. Using a broad range of probe wavelengths, the relaxation of the initially prepared excited-state wavepacket can be followed down to the sink region. The data directly visualize the directed downhill motion along the torsional reaction coordinate and suggest a barrierless excited-state isomerization in the short chain cyanine dye. Additionally, ultrafast ground-state hole and excited-state hole replica broadening is observed. While the narrow excited-state wavepacket broadens during pump-probe overlap, the ground-state hole burning dynamics takes place on a significantly longer time-scale. The experiment reported can be considered as a direct monitoring of the shape and the position of the photoprepared wavepacket on the excited-state potential energy surface.     

Photophysics of organic photostabilizers. Ab initio study of the excited-state deactivation mechanisms of 2-(2'-hydroxyphenyl)benzotriazole

Excited-state reaction paths and the corresponding energy profiles of 2-(2'-hydroxyphenyl)benzotriazole (TIN-H) have been determined with the CC2 (simplified singles-and-doubles coupled-cluster) ab initio method. Hydrogen transfer along the intramolecular hydrogen bond, torsion of the aromatic rings and pyramidization of the central nitrogen atom are identified as the most relevant photochemical reaction coordinates. The keto-type planar S(1) state reached by barrierless intramolecular hydrogen transfer is found to be unstable with respect to torsion. The latter mode, together with a moderate pyramidization of the central nitrogen atom, provides barrierless access to a S(1)-S(0) conical intersection. Only the pi-type orbitals of the aromatic rings are involved in the open-shell structures. The S(1)-S(0) conical intersection, which occurs for perpendicular geometry of the aromatic rings, is a pure biradical. From the conical intersection, a barrierless reaction path steers the system back to the enol-type minimum of the S(0) potential-energy surface, thus closing the photocycle. This photophysical pathway accounts for the remarkable photostability of the molecule.     

Photophysics of intramolecularly hydrogen-bonded aromatic systems: ab initio exploration of the excited-state deactivation mechanisms of salicylic acid

Excited state reaction paths and the corresponding energy profiles of salicylic acid have been determined with the CC2 method, which is a simplified version of singles-and-doubles coupled cluster theory. At crucial points of the potential energy hypersurfaces, single-point energy calculations have been performed with the CASPT2 method (second-order perturbation theory based on the complete active space self-consistent field reference). Hydrogen transfer along the intramolecular hydrogen bond as well as torsion and pyramidization of the carboxy group have been identified as the most relevant photochemical reaction coordinates. The keto-type planar S(1) state reached by barrierless intramolecular hydrogen transfer represents a local minimum of the S(1) energy surface, which is separated by a very low barrier from a reaction path leading to a low-lying S(1)-S(0) conical intersection via torsion and pyramidization of the carboxy group. The S(1)-S(0) conical intersection, which occurs for perpendicular geometry of the carboxy group, is a pure biradical. From the conical intersection, a barrierless reaction path steers the system back to the two known minima of the S(0) potential energy surface (rotamer I, rotamer II). A novel structure, 7-oxa-bicyclo[4.2.0]octa-1(6),2,4-triene-8,8-diol, has been identified as a possible transient intermediate in the photophysics of salicylic acid.     

Nonradiative decay mechanisms of the biologically relevant tautomer of guanine

We present the excited-state potential energy profiles of the biologically relevant 9H-keto-amino tautomer of guanine with respect to the radiationless decay via the out-of-plane deformation of the six-membered ring as well as the dissociation of NH bonds. The CASPT2//CASSCF method is employed for the reaction-path calculations. The reaction path for the out-of-plane deformation in the (1)pi pi* state leads in a barrierless way to a conical intersection with the electronic ground state. For the NH dissociation via the (1)pi sigma* state, the 9H-keto-amino tautomer is shown to have lower energy barriers than the 7H tautomers which we have studied recently. These two radiationless decay mechanisms explain the unexpected missing of the biologically relevant form in the resonant two-photon ionization spectrum of guanine in a supersonic jet. It is suggested that these ultrafast deactivation processes provide the biologically relevant tautomer of guanine with a high degree of photostability.     

Nonadiabatic dynamics of a truncated indigo model

Indigo (1) is stable when exposed to ultraviolet light. We employ electronic structure calculations and nonadiabatic trajectory surface-hopping dynamics simulations to study the photoinduced processes and the photoprotection mechanism of an indigo model, bispyrroleindigo (2). Consistent with recent static ab initio calculations on 1 and 2 (Phys. Chem. Chem. Phys., 2011, 13, 1618), we find an efficient deactivation process that proceeds as follows. After vertical photoexcitation, the S(1)(ππ*) state undergoes an essentially barrierless intramolecular single proton transfer and relaxes to the minimum of an S(1) tautomer, which is structurally and energetically close to a nearby conical intersection that acts as a funnel to the S(0) state; after this internal conversion, a reverse single hydrogen transfer leads back to the equilibrium structure of the most stable S(0) tautomer. This deactivation process is completely dominant in our semiempirical OM2/MRCI nonadiabatic dynamics simulations. The other two mechanisms considered previously, namely excited-state intramolecular double proton transfer and trans-cis double bond isomerization, are not seen in any of the 325 trajectories of the present surface-hopping simulations. On the basis of the computed time-dependent populations of the S(1) state, we estimate an S(1) lifetime of about 700 fs for 2 in the gas phase.