# Iminothioindoxyl as a molecular photoswitch with 100 nm band separation in the visible range

Mark Hoorens, Miroslav Medved', Adele D. Laurent, Mariangela Di Donato, Samuele Fanetti, Laura Slappendel, Michiel Hilbers, Ben Feringa, Wybren Jan Buma, Wiktor Szymanski

► To cite this version:

Mark Hoorens, Miroslav Medved', Adele D. Laurent, Mariangela Di Donato, Samuele Fanetti, et al. Iminothioindoxyl as a molecular photoswitch with 100 nm band separation in the visible range. Nature Communications, 2019, 10 (1), pp.17392-17402. 10.1038/s41467-019-10251-8. hal-02307621

HAL Id: hal-02307621

https://hal.science/hal-02307621

Submitted on 4 Jun 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

ARTICLE

https://doi.org/10.1038/s41467-019-10251-8

OPEN

# Iminothioindoxyl as a molecular photoswitch with 100 nm band separation in the visible range

Mark W.H. Hoorens<fnref n="1" />,², Miroslav Medved³,⁴, Adèle D. Laurent⁵, Mariangela Di Donato ID <fnref n="6" />,⁷, Samuele Fanetti ID ⁷,⁸, Laura Slappendel², Michiel Hilbers ID <fnref n="9" />, Ben L Feringa ID ², Wybren Jan Buma ID <fnref n="9" />,¹⁰ & Wiktor Szymanski ID <fnref n="1" />,²

Light is an exceptional external stimulus for establishing precise control over the properties and functions of chemical and biological systems, which is enabled through the use of molecular photoswitches. Ideal photoswitches are operated with visible light only, show large separation of absorption bands and are functional in various solvents including water, posing an unmet challenge. Here we show a class of fully-visible-light-operated molecular photo-switches, Iminothioindoxyls (ITIs) that meet these requirements. ITIs show a band separation of over 100 nm, isomerize on picosecond time scale and thermally relax on millisecond time scale. Using a combination of advanced spectroscopic and computational techniques, we provide the rationale for the switching behavior of ITIs and the influence of structural modifications and environment, including aqueous solution, on their photochemical properties. This research paves the way for the development of improved photo-controlled systems for a wide variety of applications that require fast responsive functions.

<footnote n="1">¹Department of Radiology, Medical Imaging Center, University Medical Center Groningen, University of Groningen, Hanzeplein 1, 9713 GZ Groningen, The Netherlands. ²Faculty of Science and Engineering, Centre for Systems Chemistry, Stratingh Institute for Chemistry, University of Groningen, Nijenborgh 7, 9747 AG Groningen, The Netherlands. ³Faculty of Science, Regional Centre of Advanced Technologies and Materials, Palacký University in Olomouc, Šlechtitelů 27, CZ-771 46 Olomouc, Czech Republic. ⁴Faculty of Natural Sciences, Department of Chemistry, Matej Bel University, Tajovského 40, SK-97400 Banská Bystrica, Slovak Republic. ⁵University of Nantes, CEISAM UMR CNRS 6230, BP 92208 2 Rue de la Houssiniere, 44322, Cedex 3 Nantes, France.</footnote>

<footnote n="6">⁶European Laboratory for Non Linear Spectroscopy (LENS) via N. Carrara 1, 50019 Sesto Fiorentino, Italy. ⁷INO, Istituto Nazionale di Ottica, Largo Fermi 6, 50125 Firenze, Italy. ⁸Van't Hoff Institute for Molecular Sciences, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands.</footnote>

<footnote n="9">⁹Institute for Molecules and Materials, FELIX Laboratory, Radboud University, Toernooiveld 7c, 6525 ED Nijmegen, The Netherlands. ¹⁰Department of Chemistry 'Ugo Shiff', University of Florence, via della Lastruccia 3-13, 50019 Sesto Fiorentino (FI), Italy. Correspondence and requests for materials should be addressed to B.L.F. (email: b.l.feringa@rug.nl) or to W.J.B. (email: W.J.Buma@uva.nl) or to W.S. (email: w.szymanski@umcg.nl)</footnote>
There is currently a growing interest in the development of responsive functional systems that can be controlled with light, which is a powerful, non-invasive external stimulus. Photochemical control is exerted at the molecular level through light-responsive chemical structures, i.e. photoswitches, which usually have two isomers that can be reversibly interconverted upon irradiation at different wavelengths¹,². Often, one of those isomers is less stable and thermally converts back over time to the stable isomer. The two photo-isomers of the switch differ in structure and chemical properties, which enables photochemical control of the systems in which they are embedded¹⁻⁴, including drugs and their protein targets⁵,⁶, drug delivery systems⁷,⁸, the function of hydrogels in regenerative medicine⁹, the conformation of peptides¹⁰ and nucleotides¹¹. Fascinating applications in bio-imaging¹²,¹³ and vision restoration¹⁴ are also emerging. However, for these applications, only a limited number of photoswitches is available, each with its own scope and limitations.

The selectivity in addressing the photoswitchable component in a complex functional system is crucial for its application. Because many molecular components of such systems absorb light in the UV range, a major challenge is to achieve selective switching through the design of photoswitches that can be operated in both directions using visible light. For example, in the emerging area of photopharmacology⁵,⁶,¹⁵⁻¹⁷, visible light switching is crucial to enable deep tissue penetration, especially in the 650–900 nm range³. However, most of the commonly used switches, such as diarylethenes, spiropyrans, Donor-Acceptor Stenhouse Adducts (DASAs) and fulgides, do not show absorption bands of both photo-isomers in the visible light region²,¹⁸,¹⁹. For switches that can be operated in both directions in the visible range, such as substituted azobenzenes¹ and indigoids such as indigo²⁰ and hemithioindigos²¹,²², the band separation becomes a challenge, limiting their selective bidirectional photoisomerization. Only recently, this problem has been addressed for azobenzenes by the groups of Woolley and Hecht, who developed fully-visible-light-responsive azobenzenes¹,³,²³, which - despite lower water solubility and challenging synthesis - have been successfully used for biological applications²⁴⁻²⁶. Yet, the band separation to achieve selectivity remains an unmet challenge.

In our continuous efforts to expand the limited repertoire of molecular photoswitches, we further focused on several characteristics that they should possess, besides the visible light operation with large band separation. Firstly, the photoswitch should be a small structural motif, in order to introduce it into the structure of a compound or material while affecting its original design only minimally. Secondly, it should be synthetically readily accessible. Thirdly, the parameters that control the rate of the thermal back isomerization reaction should be understood. Finally, for biological applications, the photoswitch should be able to operate under aqueous conditions. So far, realizing all these requirements in one molecular photoswitch has not been achieved.

Here we present the design, synthesis and evaluation of a class of photoswitches, which combine the photochromic dyes thioindigo and azobenzene into a photoswitch called Iminothioindoxyl (ITI). We demonstrate fully-visible (blue/orange) light switching of ITI in either direction and a large band separation between both isomers of over 100 nm. We furthermore investigate, through a comprehensive combination of synthesis, spectroscopy and theoretical calculations, the influence of the environment and chemical substitution on the switching process and re-isomerization speed of ITI. Also, we demonstrate that these spectacular photochemical properties are retained for aqueous solutions, which opens opportunities for applying ITI for reversibly controlling biological systems.

## Results

**Design and synthesis of ITI.** The design of iminothioindoxyl (ITI) is inspired by the structure of the visible-light-responsive molecular photoswitch hemithioindigo (HTI)²¹,²², which consists of half a thioindigo and half a stilbene moiety, featuring a photo-isomérizable C=C double bond. Yet, photo-isomerization is not limited to C=C double bonds. In particular, C=N photoisomerization has recently attracted attention in designing molecular photoswitches²⁷⁻³¹. Based on that, we envisioned that a molecular architecture combining azobenzene and indigoid photochromic unit could also show switching properties.

Already in the early 1900s, the chemical structures of ITI and similar compounds have been reported as dyes³². Back in 1910, Rudolf Pummerer reported the one-step synthesis of ITI by the condensation of thioindoxyl with nitrosobenzene³³. Nearly 100 years later, Soeta reported the synthesis of the same chemical structure using a Passerini-type [4 + 1] cycloaddition³⁴, also confirming through X-ray crystallography that the Z-form is the thermodynamically stable one. However, to the best of our knowledge, the behavior of these structures as molecular photoswitches has not been explored so far.

Here, we report the synthesis of six ITIs 1a-f by the condensation of thioindoxyl with substituted nitrosobenzene derivatives (Supplementary Fig. 1). Besides unsubstituted ITI 1a, two electron donating substituents (1b, 1c) and three electron withdrawing substituents (1d-1f) were placed at the R-position (Fig. 1a) to determine the influence of different substitution patterns on the photochemical properties of ITI, including absorption maxima and switching properties. Full experimental procedures and characterization is reported in Supplementary Methods and Supplementary Fig. 1-21.

**Solvent effects of ITI photo-isomerization.** To determine the influence of the medium on the photochemical properties of unsubstituted ITI 1a, absorption spectra were recorded in five solvents with different polarity (Fig. 1b, Table 1). In all solvents, the Z-isomer of ITI has an absorption band in the 400–500 nm region, with only limited solvatochromism. No clear correlation between solvent polarity and $\lambda_{max,Z}$ was observed within the group of polar solvents examined (Supplementary Fig. 61), similarly to the hemithioindigo switch³⁵. Time-dependent density functional theory (TD-DFT) calculations at the TD-M06-2X/6-311++G(2df,2p) level³⁶,³⁷, in combination with the universal solvation model based on density (SMD)³⁸ (see Supplementary Information) predicted that the band corresponds to the $S_0 \rightarrow S_2$ transition with prevailing $\pi \rightarrow \pi^*$ (HOMO $\rightarrow$ LUMO) character, while the first excited state $S_1$ is a mixed state with a significant $n \rightarrow \pi^*$ (HOMO-4 $\rightarrow$ LUMO) contribution (Supplementary Note 2, Supplementary Tables 1-3, Supplementary Figs. 23-25). In fact, due to twisting of the phenyl group out of the molecular plane (see $\theta_2$ in Fig. 2a), both excited states are partially mixed.

The photo-isomerization of 1a was followed by transient absorption spectroscopy (TA) in the millisecond time range, which revealed changes in the absorption spectrum upon irradiation at a short timescale. The transient spectra show a red-shifted absorption band, assigned to the thermally unstable E-isomer of the unsubstituted ITI 1a (Fig. 1c, Supplementary Figs. 39-48) in the 500 to 600 nm region, where Z-ITI 1a does not absorb. In all solvents, the spectrum of the E-isomer has two maxima (506–517 and at 549–554 nm), of which the most intense has been highlighted in bold (Table 1). ITI thus shows a large $\Delta\lambda_{max}$ between the two photo-isomers of over 100 nm. In comparison, HTIs usually show $\Delta\lambda_{max}$ of only 10 to 50 nm²²,³⁹.

The experimentally observed large $\Delta\lambda_{max}$ values are reproduced by the TD-DFT calculations, which further support the
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Design and absorption of ITI. **a** The structure of Iminothioindoxyl (ITI) is a hybrid of thioindigo (purple) and azobenzene (orange). The R group indicates various substituents to study the electronic effects on the photochemical properties. **b** Absorption spectra of 40 µM ITI **1a** in cyclohexane, toluene, chloroform, MeOH and DMSO. **c** Millisecond transient absorption of 400 µM ITI **1a** in MeOH at room temperature. The sample was irradiated with a 430 nm light pulse, upon which the spectrum was recorded with 1 ms delay steps. The color bar represents increased delay of transient absorption spectroscopy and the purple line represents the spectrum of 40 µM of Z-ITI **1a** in MeOH after thermal equilibration</figcaption></figure>

assignment of the absorption bands. Based on the Molecular Orbital (MO) analysis, the absorption band of the *E*-isomer corresponds to the $S_0 \rightarrow S_1$ transition with a predominant $\pi \rightarrow \pi^*$ character and a small $n \rightarrow \pi^*$ contribution (Table 1 and Supplementary Table 2). The huge bathochromic shift observed upon photoisomerization can be explained by the twist around the central double bond (C2 = N4) in the *E* isomer (see $\theta_1$ and $\theta_2$ in Fig. 2a). In the more twisted structure (*E*), the $\pi$ orbital (HOMO) is destabilized (due to less efficient overlap of 2p orbitals of C2 and N4 atoms, see Fig. 2b) leading to a smaller energy gap in the *E* isomer.

The half-life for the *E* isomer of ITI **1a** in the thermal re-isomerization process was determined at room temperature to be in the millisecond time range, which is much shorter than found for HTI²². This finding can be ascribed to the presence of a nitrogen atom in ITI that can undergo inversion (Supplementary Fig. 30), a thermal relaxation mechanism also observed for azobenzenes⁴⁰ and imine photoswitches⁴¹. The rate of nitrogen inversion is medium-dependent, with polar solvents increasing the reaction barrier⁴², which is consistent with our experimental data (Table 1, Supplementary Figs. 62, 63).

Theoretical observations of the thermal half-life are in line with the experimental ones, taking into account the limitations of continuum models to accurately describe the protic nature of MeOH. The calculations reveal that in all solvents the phenyl group is perpendicular to the molecular plane in the transition state for back isomerization from *E* to *Z*, although a concurrent (less stable) transition state with planar structure was identified in less polar solvents as well (Supplementary Note 4, Supplementary Table 5). The preference for the twisted structure is apparently related to the higher polarity of this conformation compared to the planar one (Supplementary Table 6, Supplementary Figs. 27–30) favoring its interactions with solvent molecules.

The isomerization was further studied with low-temperature NMR experiments at –60 °C. NMR spectra (Fig. 3a) showed that, upon irradiation with 455 nm light, the signals of the *Z*-isomer decreased with a concomitant rise of new signals that can be assigned to the *E*-isomer, reaching a photostationary state (PSS) of 65%. The upfield shift of proton signals upon photo-isomerization of **1a** is also predicted by calculations (see Supplementary Note 9 and Supplementary Table 12), further supporting our structure assignment. Thermal relaxation at –60 °C resulted again in the formation of the *Z*-isomer with a half-life of 6.8 ± 0.5 min without any observable degradation. An Eyring analysis, based on the determination of the back-isomerization rate at different temperatures by NMR, allowed for the calculation of the thermodynamic properties of the *E-Z* re-isomerization step (Supplementary Figs. 78–82), showing $\Delta H^\ddagger = 61.8 \pm 5.2 \text{ kJmol}^{-1}$ and $\Delta S^\ddagger = 81.6 \pm 23.4 \text{ JK}^{-1}\text{mol}^{-1}$, which results in a $\Delta G^\ddagger = 86.1 \pm 8.7 \text{ kJmol}^{-1}$ (at 298 K).

An important feature of a photoswitch is the ability to be operated photochemically in both directions exclusively with visible light. To test whether the reverse *E-Z* isomerization can be achieved photochemically, ITI **1a** in CD₃OD at –60 °C was switched to the *E*-isomer by irradiation with 455 nm (blue) light, and the rate of back-isomerization was then determined either without or with $\lambda = 595 \text{ nm}$ (orange) light irradiation. An approximately two-fold increase in the back-isomerization rate was observed under irradiation (Fig. 3b)⁴³, showing that **1a** is indeed both a T- and P-type photoswitch, while the heating effect of irradiation could be excluded (Supplementary Fig. 77). Yet it must be noted that the observation of photochemical *E* to *Z* isomerization is not of additional value at room temperature, because of the fast thermal re-isomerization.

The less stable *E*-isomer was also further characterized by measuring *E-Z* difference FTIR spectra obtained upon irradiating the sample at $\lambda = 405 \text{ nm}$ at 184 K (Fig. 3c). Importantly, these
<table><caption>Table 1 Computational studies on solvents on ITI photo-isomerization</caption><thead><tr><th rowspan="3">Solvent (ε<sub>r</sub>)</th><th colspan="2">Z-isomer</th><th colspan="3">Transition</th></tr><tr><th colspan="2">λ<sub>max,z</sub> (nm)</th><th rowspan="2">ππ*/nπ*</th><th rowspan="2">θ<sub>1</sub>(C1-C2-N4-C5)/θ<sub>2</sub>(C2-N4-C5-C6)</th><th rowspan="2">Δμ<sub>ES-GS,z</sub>(D)</th></tr><tr><th>exp.</th><th>calc.</th></tr></thead><tbody><tr><td>Cyclohexane (2.02)</td><td>416</td><td>373</td><td>S<sub>0</sub>→S<sub>2</sub><br>0.61/0.18</td><td>179.8/49.8</td><td>1.77</td></tr><tr><td>Toluene (2.37)</td><td>430</td><td>374</td><td>S<sub>0</sub>→S<sub>2</sub><br>0.61/0.20</td><td>179.8/50.2</td><td>1.86</td></tr><tr><td>CHCl<sub>3</sub> (4.71)</td><td>435</td><td>378</td><td>S<sub>0</sub>→S<sub>2</sub><br>0.61/0.20</td><td>180.0/51.4</td><td>2.26</td></tr><tr><td>MeOH (32.61)</td><td>429</td><td>398</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.58/0.34</td><td>-179.8/54.1</td><td>1.55</td></tr><tr><td>DMSO (46.83)</td><td>432</td><td>379</td><td>S<sub>0</sub>→S<sub>2</sub><br>0.61/0.24</td><td>180.0/54.1</td><td>2.82</td></tr></tbody></table><table><thead><tr><th rowspan="3">Solvent (ε<sub>r</sub>)</th><th colspan="2">E-isomer</th><th colspan="3">Transition</th></tr><tr><th colspan="2">λ<sub>max,E</sub> (nm)</th><th rowspan="2">ππ*/nπ*</th><th rowspan="2">θ<sub>1</sub>(C1-C2-N4-C5)/θ<sub>2</sub>(C2-N4-C5-C6)</th><th rowspan="2">Δμ<sub>ES-GS,E</sub>(D)</th></tr><tr><th>exp.</th><th>calc.</th></tr></thead><tbody><tr><td>Cyclohexane (2.02)</td><td><b>517</b>, 554</td><td>520</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.55/0.31</td><td>9.5/60.6</td><td>-2.85</td></tr><tr><td>Toluene (2.37)</td><td><b>510</b>, 551</td><td>519</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.55/0.31</td><td>9.5/61.7</td><td>-2.96</td></tr><tr><td>CHCl<sub>3</sub> (4.71)</td><td><b>506</b>, 549</td><td>513</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.55/0.31</td><td>9.3/62.9</td><td>-3.48</td></tr><tr><td>MeOH (32.61)</td><td><b>515</b>, 552</td><td>505</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.54/0.30</td><td>9.3/66.0</td><td>-4.17</td></tr><tr><td>DMSO (46.83)</td><td><b>514</b>, 553</td><td>503</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.55/0.31</td><td>9.0/66.9</td><td>-4.05</td></tr></tbody></table><table><thead><tr><th rowspan="2">Solvent (ε<sub>r</sub>)</th><th colspan="2">Δλ<sub>max</sub> (nm)</th><th rowspan="2">t<sub>1/2</sub> (ms)</th><th colspan="2">ΔG<sup>‡</sup><sub>Z-E</sub> (kcal/mol)</th><th colspan="2">Transition state</th></tr><tr><th>exp.</th><th>calc.</th><th>exp.</th><th>calc.</th><th>θ<sub>1</sub>(C1-C2-N4-C5)/θ<sub>2</sub>(C2-N4-C5-C6)</th><th>μ<sub>GS,TS</sub>(D)</th></tr></thead><tbody><tr><td>Gas phase (1.00)</td><td>-</td><td>-</td><td>-</td><td>-</td><td>NA (13.2)</td><td>NA (0.0/0.0)</td><td>NA (1.20)</td></tr><tr><td>Cyclohexane (2.02)</td><td><b>101</b>, 138</td><td>147</td><td>9.5 ± 0.4</td><td>14.1</td><td>12.8 (12.8)</td><td>0.0/90.4 (-0.1/0.1)</td><td>3.85 (1.30)</td></tr><tr><td>Toluene (2.37)</td><td>80, <b>121</b></td><td>145</td><td>12.4 ± 0.9</td><td>14.2</td><td>12.7 (12.9)</td><td>0.0/90.4 (-0.1/0.1)</td><td>3.97 (1.33)</td></tr><tr><td>CHCl<sub>3</sub> (4.71)</td><td>71, <b>114</b></td><td>135</td><td>16.9 ± 1.2</td><td>14.4</td><td>13.3 (13.7)</td><td>0.0/90.4 (-0.1/0.1)</td><td>4.44 (1.52)</td></tr><tr><td>MeOH (32.61)</td><td><b>86</b>, 123</td><td>107</td><td>18.5 ± 1.4</td><td>14.4</td><td>14.4 (NA)</td><td>0.0/87.5 (NA)</td><td>5.25 (NA)</td></tr><tr><td>DMSO (46.83)</td><td><b>82</b>, 121</td><td>124</td><td>23.3 ± 2.0</td><td>14.6</td><td>13.5 (NA)</td><td>0.0/90.2 (NA)</td><td>4.92 (NA)</td></tr></tbody></table>

Solvatochromic shifts of $\lambda_{max}$ for the Z (Top) and E (Middle) isomers of ITI 1a. Experimental $\lambda_{max,E}$ values are obtained from TA that show two absorption maxima which are both reported and the maximum, of which the one with the highest absorption is highlighted in bold. Theoretical $\lambda_{max}$ values and the difference of GS and ES dipole moments ($\Delta\mu_{ES-GS}$) were obtained at the SMD-TD-M06-2X/6-311++G(2df,2p) level using the SMD-M06-2X/6-31+G(d) geometries, from which also twisting angles $\theta_1$ and $\theta_2$ were derived (see Fig. 2a). Bottom: Thermal relaxation of ITI 1a. Experimental half-lives were calculated from ms TA. The GS dipole moments for the transition state ($\mu_{GS,TS}$) were obtained at the SMD-M06-2X/6-31+G(d) level, at which also the twisting angles $\theta_1$ and $\theta_2$ as well as the activation barriers for thermal relaxation were derived (see Fig. 2a). The data in parentheses refer to a planar TS structure

spectra were acquired with the sample in a KBr pellet, demonstrating that isomerization also occurs at the solid state. The main spectral features related to structural differences between the two isomers are fairly well reproduced by the DFT calculations (see band assignment in Supplementary Table 11, Supplementary Figs. 35 and 36 and Supplementary Note 7.).

**Z-E isomerization of ITI is a fast process.** Transient absorption measurements with sub-picosecond time resolution were performed to determine the timescale of forward Z to E isomerization of ITI, which is expected to be very fast, based on structural analogies with HTIs and azobenzenes$^{22,43}$. For unsubstituted ITI 1a, the spectra recorded immediately after excitation with $\lambda = 400$ nm light are dominated by a very broad excited state absorption band with an intensity that rapidly decays, leaving a constant weak differential signal as shown in the time-resolved spectra reported in Fig. 4a and the kinetic traces in Fig. 4b. Importantly, the long-living signal matches the one measured on the millisecond timescale (Fig. 1c), and can thus assigned unambiguously be as the Z-E difference spectrum. The very fast decay of the excited state absorption band indicates that isomerization itself is a very fast process, since the system has to reach the conical intersection (CI) leading to the formation of the Z and E isomers in their respective ground states before the deactivation of the excited states. In order to get additional kinetic information on the process, we measured the pump-probe anisotropy by recording the transient spectra with parallel and perpendicular polarization of the pump beam with respect to the probe. Interestingly, the resulting

anisotropy signal, reported in Fig. 4c, shows a fast rise component, on a timescale of a few hundred femtoseconds, and a slower decay, occurring within 12–16 ps. The timescale of the anisotropy decay is in line with what has been observed for azobenzene in solution$^{44}$. The rise of the anisotropy within the initial 500 fs indicates that a significant charge redistribution rapidly occurs once the molecule starts to move on the excited state potential energy surface towards the conical intersection region, in line with the computed large difference in transition dipole moments for the Z and E forms (Table 1). It is worth noticing that a similar rise in the anisotropy in a few hundred fs has been previously observed for rhodopsin, which is known to isomerize on an ultrafast timescale and interpreted in terms of rapid and substantial change in the charge distribution of the molecule due to the activation of the vibrational modes leading to isomerization$^{45}$.

Our calculations indicate that the bright state of ITI is the $S_2$ state. Taking into account the observed fast excited state decay, we therefore envisioned the excited state relaxation pathway to be similar to that of azobenzene. To extract the time constants describing the photodynamics of the system, we fitted the transient isotropic data with the kinetic scheme shown in Fig. 4e, retrieving the lifetimes reported therein and the Species-Associated Difference Spectra (SADS) of the transient intermediates (Fig. 4d). Upon excitation to $S_2$, the system rapidly undergoes internal conversion towards $S_1$, with a time constant below the time resolution of our measurements. This results in an unreasonable spectral shape for this state, which is not shown in Fig. 4d. The remaining SADS are assigned to the $S_1$ state (black line), to the hot Z isomer (red curve) and the E isomer (blue
<figure><img src="image_4.png" /><figcaption>**Fig. 2** Computational studies on solvents on ITI photo-isomerization. **a** Angles $\theta_1$ (top) and $\theta_2$ (bottom). **b** Structures of the Z and E forms of ITI **1a** in MeOH with the numbering of atoms in the central part of a molecule, molecular orbitals involved in the observed electronic transition (energies in Hartrees) and electron density difference (EDD) plots showing the decrease (blue) and increase (red) of the electron density upon excitation obtained at the SMD-TD-M06-2X/6-311++G(2df,2p)//SMD-M06-2X/6-31+G(d) level of theory</figcaption></figure>

curve). The very short $S_2$ lifetime is again similar to what is known for azobenzene, for which a value of 50 fs has been recently determined$^{44,46}$. The decay of the broad $S_1$ excited state band within 320 fs and the rise of anisotropy on the same timescale indicate that ITI reaches the conical intersection region on a time scale competing with vibrational relaxation in $S_1$. From there, the molecule relaxes to the ground state of either the Z and E isomers, where vibrational cooling takes place on a time scale of 10 ps.

Support to our hypothesis that isomerization starts from a hot $S_1$ state comes from the computation of the forces acting on the individual atoms of ITI in $S_2$ and $S_1$ after vertical excitation, showing that the molecule undergoes more pronounced structural changes in the $S_1$ state (for more details see Supplementary Fig. 26, Supplementary Note 3 and Supplementary Table 4). The presence of a nitrogen atom in the isomerizing double bond opens the possibility for isomerization to occur through either an inversion or rotation mechanism. The negligible change in the excited state relaxation time scale observed in solvents with different viscosity (see Supplementary Fig. 38) in first instance favors an inversion mechanism, although most probably the simple vision of motion along a single reaction coordinate is not realistic, as recently pointed out for azobenzene$^{44}$.

**Substituent effects on ITI photo-isomerization.** The influence of the substituents on photoswitching of ITI was studied using a small library of ITIs with either an electron donating (**1b,c**) or an electron withdrawing group (**1d-f**). As shown in Fig. 5, electron donating groups (EDG) result in a slight red-shift of $\lambda_{max,Z}$ and increased absorption, while electron withdrawing groups (EWG) result in a slight blue-shift of $\lambda_{max,Z}$ and decreased absorption (Supplementary Fig. 64). Theoretical calculations reproduce this trend and show that the auxochromic effects are mainly due to the twist around the =N-C- central single bond ($\theta_2$, Fig. 2a). Indeed, $\theta_2$ is smaller for **1b,c**, leading to a more planar structure and favoring the electron delocalization (Supplementary Fig. 31 and Supplementary Table 8) upon excitation and increasing $\lambda_{max,Z}$. In the ground state, EDGs increase the electron density on the phenyl ring which tends to “planarize” to increase conjugation with the thioindoxyl moiety in accordance with similar auxochromic effects have been observed in HTIs$^{47}$.

Isomerization of the differently substituted ITIs was measured in MeOH upon irradiation with $\lambda = 430$ nm light (Fig. 5b, Supplementary Figs. 49–58). A new absorption band was found for all the substituted ITIs and for electron donating ITIs **1b** and **c** an impressively large $\Delta\lambda_{max}$ of over 100 nm was observed. ITI **1b** was dissolved in MeOH and irradiated with 400 nm while cooled to –60 °C (Fig. 5d). Compared to the thermally adapted state, isomerization resulted in a clear change in color. Switching for several cycles of **1b** in MeOH did not result in observable degradation (Fig. 5e). For all ITIs, the quantum yield for forward switching was estimated to be between 4 and 6%, which is relatively low compared to many other photoswitches$^{21}$. No clear correlation between Hammett parameter R and the quantum yield (Supplementary Note 13, Supplementary Table 14) for the single studied position was found, meaning that both electron withdrawing and electron donating groups are tolerated.

Our calculations show that the auxochromic effects on $\Delta\lambda_{max}$ can be explained by a combination of geometrical and electronic effects (Supplementary Note 5). While $\theta_2$ is governing the auxochromic effects for the Z and E isomers in the same way ($\theta_2$ is larger for E than for Z but the extent to which E and Z are influenced by a substituent is similar), a twist around the C = N central double bond ($\theta_1$) is only observed for the E isomer. The $\theta_1$ twist, being more pronounced for EDG substituents (**1b,c**), leads to a stronger destabilization of the $\pi$ orbital (HOMO) of the E isomer for these substituents compared to the Z isomer. Such geometrical feature partly contributes to the decrease of the $\Delta\lambda_{max}$ when going from **1b,c** to **1a,d,e,f**. In addition, the change of the dipole moment upon excitation for the E form decreases from 2.37 D (**1b**) to –5.85 D (**1f**) in methanol following the nature of the substituents (Table 2). We have found that the more negative $\Delta\mu$, the larger destabilization of the ES with respect to GS. This electronic effect also contributes to a smaller $\Delta\lambda_{max}$ for EWG substituents (Supplementary Table 7, 8 and Supplementary Fig. 31).

Apart from changes in the absorption spectra of Z and E, substituents also influence the rate of thermal relaxation of the E isomer (Table 2). No clear correlation between the Hammett parameter and the half-lives of the E isomer was observed, albeit the data suggested a trend in EWG groups results in faster
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 3** NMR and IR spectroscopy. **a** NMR spectra of ITI 1a in CD₃OD at -60 °C for the thermally adapted, irradiated and again thermally adapted sample **b**: E-Z isomerization of ITI 1a at -60 °C in CD₃OD, recorded without (thermal) and with $\lambda$ = 595 nm irradiation (Supplementary Fig. 77, Supplementary Note 12) **c**: E-Z FTIR difference spectrum recorded upon irradiation at 405 nm in KBr at 184 K for ITI 1a. Comparison of experimental and theoretical IR difference spectra of 1a. Experimental FTIR difference spectrum of the compound 1a was obtained from the spectra in the dark and under 405 nm light measured at 184 K in a KBr pellet (Supplementary Figs. 84–86). Simulated difference spectrum was obtained from scaled harmonic GS IR spectra (scaling factor $f$ = 0.98) of the E- and Z-isomers of 1a in acetonitrile calculated with at the SMD-B3LYP/6-31++G(d,p) level. The experimental FTIR spectra are also reported in Supplementary Fig. 84 for better visualization. Further IR characterization can be found in Supplementary Notes 7,8, Supplementary Figs. 33–37 and Supplementary Table 11</figcaption></figure>

re-isomerization (Supplementary Figs. 32, 65, 66, Supplementary Note 6, and Supplementary Table 9). The same correlation between Hammett parameter R and the half-lives of the *E* isomer was observed at -60 °C upon 455 nm irradiation in the NMR experiment (Supplementary Notes 9 and 11, Supplementary table 13, Supplementary Figs. 67–76). DFT results were in

line with these observations, revealing that the weak correlation of activation energy with the Hammett constants could be caused by qualitatively different relaxation paths for the EDG- and EWG-substituted (and neutral) ITIs. Whereas the *E-Z* relaxation proceeded through a planar TS structure in the case of **1b-c**, **1a,d-f** adopted a twisted conformation in the TS (Supplementary Fig. 31). The different behavior is a result of interplay between the stabilization of the TS due to π-electron delocalization (favoring the planar conformation) and the stabilization due to polarity of the TS (favoring the more polar twisted structure). By decreasing the electron density on the phenyl ring, EWG substituents enhance the interaction of the 2p orbital on nitrogen with π-orbitals of the phenyl ring favoring the twisted structure (Supplementary Figs. 27, 30 and 32).

**Isomerization of ITI in aqueous solutions.** In the field of photopharmacology, photo-control over the stereochemistry of a double bond is used to establish a difference in biological activity between both photo-isomers, as has been demonstrated for azobenzene and hemithioindigo photoswitches⁶,⁴⁸. For such biological applications of photoswitches, solubility at medicinally relevant conditions and photo-isomerization under aqueous and physiological conditions are crucial, yet are rarely observed for fully-visible-light switches. For example, photo-isomerization of HTI at physiological conditions has not been reported. To evaluate the performance of ITI in aqueous solutions, unsubstituted ITI 1a was dissolved in phosphate buffered saline (PBS, pH 7.4, 1.7% DMSO) at ~30 µM. Irradiation with 400 nm light did not results in observable degradation (Supplementary Fig. 88). We also demonstrated that ITI has resistance against glutathione (GSH), which is found in concentrations up to 10 mM in cells and is the key factor for degradation of other molecular photoswitches⁴⁹.

Isomerization of ITIs in aqueous PBS (pH 7.4, 6.7% DMSO) was studied using the most red-light shifted p-MeO-ITI **1b** (Fig. 5C, Supplementary Figs. 59, 60) with ms transient absorption spectroscopy. The Z isomer of **1b** has an absorption maximum at 459 nm. Upon irradiation with blue light, the *E* isomer was observed with an absorption maximum at 560 nm, demonstrating that a spectacular difference of absorption maxima is also maintained in aqueous solutions (Fig. 5c). From the same experiment, the half-life of the *E* isomer was found to be 10.0 ± 0.8 ms at room temperature.

**Discussion**

For application in biological systems, new and improved switches are needed. This is underlined e.g. by a recent report by the group of Thorn-Seshold⁴⁸, in which the first HTI-based photo-controlled pharmacophore was reported. This study demonstrates both the potential of indigoid-based photoswitches as well as the need for improved band separation of photo-isomers and improved water solubility.

Here we described the discovery of Iminothioindoxyls, a class of small, synthetically accessible visible-light photoswitches with excellent photochemical properties, showing very fast switching and an absorption band separation of photo-isomers of over 100 nm. Importantly, ITIs switch in solid state and in solvents ranging in polarity from cyclohexane to water, being therefore suitable for a very wide range of applications, varying from responsive materials to photopharmacology.

ITIs show unique properties when compared to other fully-visible-light-responsive photoswitches. A promising feature of ITIs is the millisecond half-life, making them useful for applications requiring fast responses. Indeed, many biological processes, such as signal transduction and neuronal
<figure><img src="image_7.png" /><figcaption>**Fig. 4** Ultra-fast Transient Spectroscopy of ITI. **a** Transient absorption spectra of unsubstituted ITI 1a recorded in methanol with excitation at 400 nm. **b** Representative kinetic traces (open symbols) and fits obtained from target analysis (continuous line), **c** Time-resolved anisotropy, the initial 3 ps are shown in the inset, **d** Species-Associated Decay Spectra (SADS), obtained by analyzing the kinetic traces with the kinetic model depicted on the right-bottom side of the figure. The black curve represents the S₁ state, the red curve hot Z isomer and the blue curve the E isomer. **e** Proposed model for photo-isomerization</figcaption></figure>

communication, operate at the millisecond scale and their
photomodulation has been achieved with quickly re-
isomerizing switches⁵⁰,⁵¹. Furthermore, ITIs forward switch-
ing is faster and shows better band separation than hemi-
thioindigo, while also operating on a completely different
mechanism for thermal relaxation. Finally, photo-isomerization
of HTI in aqueous solutions at physiological pH has so far not
been realized, while for ITI it could be readily observed. Also if
compared to red-shifted azobenzenes, ITIs present favorable
properties: they are slightly smaller in structure and syntheti-
cally more accessible, showing faster switching and a lar-
ger absorption band separation between the two isomers, high
stability under irradiation and under heavily reducing condi-
tions such as those encountered in living cells.

Currently, the fast re-isomerization of ITIs prevents the use of
their bi-directional photochemical isomerization at room tem-
perature. To fully exploit the various properties of this class of
photoswitches, an increased build-up and a longer lifetime of the
E isomer is needed. This could be achieved through judiciously
substitution patterns that improve the quantum yield and
increasing the thermal barrier of re-isomerization. Similar situa-
tions have occurred in the past when other types of switches have
been developed. In view of the successful studies that have

followed to optimize these switches, we are confident that also for
ITIs this will be a realistic target. We therefore consider the
discovery of ITIs a break-through in the field of photocontrol,
providing the starting point for developing improved photo-
switches, resulting in major opportunities towards responsive
systems well beyond those offered by the current very limited
repertoire of all-visible light switches.

**Methods**

**Organic synthesis.** All reported starting materials, chemical reagents and organic solvents in this study were bought from Sigma-Aldrich, Acros, Fluka, Fischer, TCI and were used as received. Dry DCM was purified by passage through an MБraun SPS-800 solvent purification column. All aqueous solutions were prepared using deionized water. Kieselgel 60, F₂₅₄ silica gel plates (Merck, TLC silica gel 60 F₂₅₄) were used for TLC (Thin Layer Chromatography) analysis and UV light of 254 nm and potassium permanganate solution (KMnO₄) were used for the detection of compounds. Drying of solutions was performed using dry MgSO₄ and solvents and other volatiles were removed using a rotary evaporator.

**Analytical procedures.** Nuclear Magnetic Resonance (NMR) spectra were recorded using an Agilent Technologies 400-MR (400/54 Premium Shielded) spectrometer (400 MHz), at room temperature (22–24 °C), unless indicated otherwise. The multiplicities of the signals are reported as follows: *s* (singlet), *d* (doublet), *t* (triplet), *q* (quartet) or *m* (multiplet). All ¹³C-NMR spectra are ¹H-broadband
<figure><img src="image_8.png" /><figcaption>**Fig. 5** Spectroscopy studies on the substituent effects on ITI photo-isomerization. **a** Absorption spectra of 40 µM ITIs **1a-f** in MeOH. **b** Transient absorption spectra of ITIs **1a-f** in MeOH after irradiation at 430 nm after 3 ms delay. **c** Transient absorption spectroscopy of 120 µM p-MeO-ITI **1b** in aqueous PBS buffer (6.7% DMSO), irradiated with a 10 ns 430 nm light pulse and spectra recorded with 1 ms delay steps. The purple line indicates the absorption spectrum of 120 µM Z-ITI **1b** in aqueous PBS buffer (6.7% DMSO). The color bar represents increased delay in transient absorption spectroscopy. **d** Cuvettes 1 and 2 contain 200 µM ITI **1b** in MeOH. Left: both thermally adapted. Middle: cuvette 2 irradiated with 400 nm light while cooled at -60 °C in acetone bath. Right: reheating of cuvette 2 to room temperature. **e** Three cycles of photo-isomerisation of 100 µM **1b** in MeOH, thermally adapted and switched with 400 nm light, while cooled at -60 °C in acetone bath (Supplementary Fig. 89)</figcaption></figure>

decoupled. Melting points (Mp) were measured using a Stuart analogue capillary melting point SMP11 apparatus. High-resolution mass spectrometric (HRMS) measurements were performed using a Thermo scientific LTQ OrbitrapXL spectrometer, which is equipped with ESI ionization. In the experimental procedures, the mass of the molecule-ion [M + H]⁺ are reported in m/z-units. Absorption spectra were measured using an Agilent 8453 UV/Vis diode array. All solutions for absorption spectra were prepared in Uvasol® grade solvents and were measured in quartz cuvettes with a 1 cm path-length. Purity was determined using LCMS, for which the following setup was used: Column: ACQUITY UPLC® HSS T3 1.8 µm, 2.1 × 150 mm; Detection: Total Ion Count (TIC), λ₁ = 254 nm, λ₂ = 430 nm; Flow: 0.3 mL/min; Eluent A: 0.1% formic acid in HPLC grade demineralized H₂O; Eluent B: 0.1% formic acid in acetonitrile; Gradient Program: (0-1 min) 5% eluent B; (1-8 min) linear gradient to 90% eluent B; (8-11 min) 90% eluent B; (11-12 min) linear gradient to 5% eluent B; (12-17 min) 5% eluent B.

**Computational studies.** The ground state (GS) structures of the Z/E-isomers and the GS transition state (TS) of the backward reaction (E→Z; thermal relaxation process) for ITIs (**1a-f**) were optimized at the M06-2X level³⁶ using the 6-31+G(d) atomic basis set³⁷, since this exchange-correlation functional is known to perform well not only for the GS thermochemistry, but also in describing excited states⁵². In addition, ITIs are not subject to the known TD-DFT limitations such as charge-transfer (see Supplementary Table 6), double excitations (see t1 and t2 amplitudes), singlet-triplet transition, etc. All minima were checked against the presence of imaginary frequencies. The TS structures were obtained by geometry optimization starting from a structure with the angle C2-N4-C5 set to very close to 180°. This choice was based on the potential energy scan for the out-of-plane distortion from the in-plane-TS structure showing that the distortion is energetically unfavourable (Supplementary Fig. 28). The optimized TS structures (first-order saddle points) were checked against the presence of a single imaginary frequency. The optimized GS structures are presented in Supplementary Note 1 and Supplementary Fig. 22. The solvent effects were considered employing the solvation Model based on Density (SMD)³⁸. Cyclohexane (CHX), toluene (TOL), chloroform (CHL), methanol and dimethylsulfoxide (DMSO) are used consistently with experimental data. The IR spectra were simulated at the B3LYP/6-31++G(d,p) level⁵³,⁵⁴ which was found to provide a reasonable agreement with the experimental FTIR spectra for the Z-isomers (Supplementary Note 7, Supplementary Figs. 33-36). The IR band assignment was based on the potential energy distribution (PED) analysis⁵⁵ by using the VEDA 4 program⁵⁶. Vertical excitation energies (VEE) were obtained with a larger basis set, namely 6-31++G(2df,2p). SMD was combined with the

corrected linear response (cLR) approach⁵⁷ to model VEE within the non-equilibrium regime. (TD)-DFT calculations were performed using the Gaussian09 and Gaussian16 programs⁵⁸,⁵⁹. All Gaussian default thresholds and algorithms were used except for improving optimization. In the latter case, a threshold of 10⁻⁵ a.u. on average residual forces was imposed, a self-consistent field convergence criterion of 10⁻¹⁰ a.u., and the use of the ultrafine DFT integration grid. Gas phase CC2 and ADC(2) calculations of the excitation energies were performed using aug-cc-pVTZ basis set with the Turbomole program [TURBOMOLE V6.6 2014, a development of University of Karlsruhe and Forschungszentrum Karlsruhe GmbH, 198-2007, TURBOMOLE GmbH, since 2007; available from http://www.turbomole.com]. NMR shieldings for the protons in Z/E-isomers of **1a** were obtained with the Gauge-Independent Atomic Orbital (GIAO) method⁶⁰ with the B3LYP functional and the 6-31++G(d,p) basis set. Chloroform environment was treated using the SMD model. Proton shielding for TMS selected as reference was computed under the same conditions using the M06-2X/6-31+G(d) geometry.

**Ultrafast spectroscopy.** Ultrafast transient absorption spectra of unsubstituted Iminothiodoxyl **1a** were measured on a system consisting of a Ti:sapphire laser oscillator (Spectra Physics Tsunami) and regenerative amplifier system (BMI Alpha 1000) which produced pulses of 100 femtosecond at 800 nm with an average output power of 450 to 500 mW. Excitation pulses at a wavelength of 400 nm were obtained by second harmonic generation of the fundamental laser output in a 2 mm thick BBO crystal. For all measurements in methanol, the pump beam polarization was set either to perpendicular or parallel with respect to the parallel probe beam by rotating a λ/2 plate. For measurement in solvents other than methanol, polarization was set to magic angle so as to exclude rotational contributions to the transient signal. From the parallel and perpendicular intensities the anisotropy r(t) is calculated using Eq. (1).

$$r(t) = \frac{I_{\parallel} - I_{\perp}}{I_{\parallel} + 2I_{\perp}} \qquad (1)$$

where $I_{\parallel}$ and $I_{\perp}$ are the signal intensity respectively recorded with parallel and perpendicular pump polarization. The isotropic signal in methanol is obtained from the parallel and perpendicular signals using Eq. (2).

$$I_{Iso} = \frac{I_{\parallel} + 2I_{\perp}}{3} \qquad (2)$$

The excitation powers were on the order of 50 to 100 nJ. The probe pulses were generated upon focusing the 800 nm radiation beam partially on a 2 mm thick
<table><caption>Table 2 Computational studies on substituent effects on ITI photo-isomerization</caption><thead><tr><th rowspan="2"></th><th rowspan="2">R (Hammet constant σ)</th><th colspan="2">Z-isomer</th><th colspan="4">Transition</th></tr><tr><th>exp.</th><th>calc.</th><th>ππ*/nπ*</th><th>θ<sub>1</sub>(C1-C2-N4-C5)/θ<sub>2</sub>(C2-N4-C5-C6)</th><th>Δμ<sub>ES-GS,Z</sub> (D)</th><th>ε (mol<sup>-1</sup>cm<sup>-1</sup>)</th></tr></thead><tbody><tr><td>1a</td><td>H (0.00)</td><td>429</td><td>398</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.58/0.38</td><td>-179.8/54.1</td><td>1.55</td><td>4300</td></tr><tr><td>1b</td><td>MeO (-0.27)</td><td>448</td><td>413</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.62/0.29</td><td>-179.1/38.6</td><td>4.43</td><td>11000</td></tr><tr><td>1c</td><td>Me (-0.17)</td><td>434</td><td>406</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.61/0.32</td><td>-179.9/50.3</td><td>2.66</td><td>5700</td></tr><tr><td>1d</td><td>COOMe (0.45)</td><td>427</td><td>399</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.56/-0.15</td><td>-179.1/60.8</td><td>0.85</td><td>2300</td></tr><tr><td>1e</td><td>CF<sub>3</sub> (0.54)</td><td>424</td><td>391</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.54/0.37</td><td>-179.1/62.0</td><td>-0.62</td><td>2100</td></tr><tr><td>1f</td><td>NO<sub>2</sub> (0.78)</td><td>-</td><td>390</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.53/-0.33</td><td>-178.9/67.3</td><td>-0.26</td><td>2600</td></tr></tbody></table>

<table><thead><tr><th rowspan="2"></th><th rowspan="2">R (Hammet constant σ)</th><th colspan="2">E-isomer</th><th colspan="4">Transition</th></tr><tr><th>exp.</th><th>calc.</th><th>ππ*/nπ*</th><th>θ<sub>1</sub>(C1-C2-N4-C5)/θ<sub>2</sub>(C2-N4-C5-C6)</th><th>Δμ<sub>ES-GS,E</sub> (D)</th></tr></thead><tbody><tr><td>1a</td><td>H (0.00)</td><td>515, 552</td><td>505</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.54/0.30</td><td>9.3/66.0</td><td>-4.16</td></tr><tr><td>1b</td><td>MeO (-0.27)</td><td>553</td><td>533</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.58/-0.32</td><td>12.4/52.9</td><td>2.37</td></tr><tr><td>1c</td><td>Me (-0.17)</td><td>511, 548</td><td>519</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.58/-0.31</td><td>9.9/63.2</td><td>-1.01</td></tr><tr><td>1d</td><td>COOMe (0.45)</td><td>503</td><td>484</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.52/-0.31</td><td>0.26/93.0</td><td>-4.45</td></tr><tr><td>1e</td><td>CF<sub>3</sub> (0.54)</td><td>500</td><td>482</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.52/-0.31</td><td>2.23/86.1</td><td>-5.81</td></tr><tr><td>1f</td><td>NO<sub>2</sub> (0.78)</td><td>501</td><td>470</td><td>S<sub>0</sub>→S<sub>1</sub><br>0.54/0.30</td><td>0.00/-92.9</td><td>-5.85</td></tr></tbody></table>

<table><thead><tr><th rowspan="2"></th><th rowspan="2">R (Hammet constant σ)</th><th colspan="2">Δλ<sub>max</sub> (nm)</th><th colspan="2">t<sub>1/2</sub> (ms)</th><th colspan="2">E<sub>a,E-Z</sub> (kcal/mol)</th><th colspan="4">Transition state</th></tr><tr><th>exp.</th><th>calc.</th><th>exp.</th><th>calc.</th><th>exp.</th><th>calc.</th><th>α(C2-N4-C5)/θ<sub>2</sub>(C2-N4-C5-C6)</th><th>ΔμGS-TS,Z (D)</th><th>φ<sub>Z-E</sub> (%)</th></tr></thead><tbody><tr><td>1a</td><td>H (0.00)</td><td>86, 123</td><td>107</td><td>18.5 ± 1.4</td><td>14.4</td><td>14.4</td><td>177.2/87.5</td><td>-0.02</td><td>6.2</td></tr><tr><td>1b</td><td>MeO (-0.27)</td><td>105</td><td>120</td><td>12.7 ± 0.5</td><td>14.2</td><td>13.0</td><td>177.7/0.0</td><td>-2.45</td><td>4.5</td></tr><tr><td>1c</td><td>Me (-0.17)</td><td>77, 114</td><td>113</td><td>21.1 ± 1.2</td><td>14.5</td><td>14.0</td><td>177.7/0.0</td><td>-2.68</td><td>5.4</td></tr><tr><td>1d</td><td>COOMe (0.45)</td><td>76</td><td>85</td><td>4.0 ± 0.3</td><td>13.6</td><td>13.1</td><td>177.6/92.0</td><td>0.38</td><td>6.3</td></tr><tr><td>1e</td><td>CF<sub>3</sub> (0.54)</td><td>76</td><td>91</td><td>9.9 ± 1.0</td><td>14.1</td><td>13.6</td><td>177.6/90.1</td><td>0.94</td><td>4.9</td></tr><tr><td>1f</td><td>NO<sub>2</sub> (0.78)</td><td>-</td><td>80</td><td>2.8 ± 0.5</td><td>13.4</td><td>12.0</td><td>177.8/90.1</td><td>2.48</td><td>4.1</td></tr></tbody></table>

Shifts of λ<sub>max</sub> for the Z (Top) and E (Middle) isomers of ITIs 1a-f in MeOH. Experimental λ<sub>max,E</sub> values are obtained from TA that show two absorption maxima which are both reported with the maximum with the highest absorption highlighted in bold. Theoretical λ<sub>max</sub> values and the difference of GS and ES dipole moments (Δμ<sub>ES-GS</sub>) were obtained at the SMD-TD-M06-2X/6-311++G(2df,2p) level using the SMD-M06-2X/6-31+G(d) geometries, from which also twisting angles (θ<sub>1</sub> and θ<sub>2</sub>, Fig. 2a) were derived. Bottom: Thermal relaxation of ITIs 1a-f in MeOH. Experimental half-lives were calculated from ms TA spectroscopy. The differences of dipole moment of the transition state and that of the Z-form in their GS (Δμ<sub>GS-TS,Z</sub>) were obtained at the SMD-M06-2X/6-31+G(d) level, at which also the transition state twisting angles (θ<sub>1</sub> and θ<sub>2</sub>) as well as the activation barriers for thermal relaxation were derived

sapphire window, after which it was passed through the sample. Subsequently the white light probe was sent to a flat field monochromator which was coupled to a home-made CCD detector (http://lens.unifi.it/ew). Transient spectra were recorded in a time interval spanning up to 500 ps. All measurements were performed in a quartz cell (2 mm thick) mounted on a movable stage in order to refresh the solution and avoid undesired photochemical degradation of the sample.

Analysis of the Transient data was performed using Singular Value Decomposition (SVD)⁶¹ and global analysis⁶², which allows the simultaneous fit at all the measured wavelengths with a combination of exponential decay functions. The kinetic scheme employed for data analysis, involving fast internal conversion among two close-lying excited states and excited state decay associated to partial Z-E isomerization, is shown in Fig. 4 of the main text. Data analysis has been performed using the software GLOTARAN⁶³.

**Nanosecond transient absorption spectroscopy.** Nanosecond transient absorptions were recorded with an in-house assembled setup. For all ITIs and all solvents,

an excitation wavelength of 430 nm was used. The excitation wavelength of 430 nm was generated using a tunable Nd:YAG-laser system (NT342B, Ekspla) comprising the pump laser (NL300) with harmonics generators (SHG, THG) producing 355 nm to pump an optical parametric oscillator (OPO) with SHG connected in a single device. The laser system was operated at a repetition rate of 5 Hz. The probe light running at 10 Hz was generated by a high-stability short arc xenon flash lamp (FX-1160, Excelitas Technologies) using a modified PS302 controller (EG&G). Using a 50/50 beam splitter, the probe light was split equally into a signal beam and a reference beam with and focused on the entrance slit of a spectrograph (SpectraPro-150, Princeton Instruments). The probe beam (A = 1 mm²) was passed through the sample cell and orthogonally overlapped with the excitation beam on a 1 mm × 1 cm area. The excitation energy was recorded by measuring the excitation power at the back of an empty sample holder. In order to correct for fluctuations in the flash lamp spectral intensity, the reference was used to normalize the signal. Both beams were recorded simultaneously using a gated intensified CCD camera (PI-MAX3, Princeton Instruments) which has an adjustable gate of minimal 2.9 ns. A delay generator (DG535, Stanford Research Systems, Inc.) was used to time the
excitation pulse, the flash lamp, and the gate of the camera. The setup was controlled by an in-house written LabView program.

**In situ NMR irradiation experiments.** NMR spectra were recorded with an Agilent Technologies Inova 500 Spectrometer (500 MHz), and for in situ irradiation, a set-up based on LED and an option fiber were used, according to a reported system⁶⁴. The fiber-optic cable (M28L05; Ø400 µm, 0.39 NA, SMA-SMA Fiber Patch Cable, 5m) and the LEDs were purchased from Thorlabs: royal blue 455 nm Fiber-coupled LED(M455F1, 11.0 mW); amber 595 nm Fiber-coupled LED (M595F2, 11.0 mW). NMR tubes were purchased from Wilmad-LabGlass (SP Scienceware):WGS-5BL, Coaxial Insert for 5 mm NMR Sample Tube and 535-PP-7, 5 mm Thin Wall Precision NMR Sample Tube 7" L, 600 MHz.

**FTIR.** Low-temperature FTIR spectra were recorded on a FTIR Bruker IFS 120 HR spectrometer with maximum resolution 0.002 cm⁻¹. For current measurements spectra were registered with 1 cm⁻¹ spectral resolution. The instrument is equipped with a globar IR source and a MCT detector. The sample has been cooled using a liquid helium cold tip closed cycle cryostat (minimal nominal temperature 5 K), temperature has been monitored at the sample position using a K-type thermo-couple (reading error 0.1 K)⁶⁵.

The in situ irradiation source was a 80 mW laser diode, with a spot size of 6 × 4 mm, centered at 405 nm (FWHM ~10 nm). The sample was prepared as a KBr pellet, and contained in a home-made cell equipped with two calcium fluoride windows. Spectra without and under irradiation were measured at 184 K.

**Data availability**

The authors declare that the data supporting the findings of this study are available within the paper and its supplementary information files. Additional data on methods used are available from the corresponding author upon reasonable request.

Received: 15 January 2019 Accepted: 25 April 2019

Published online: 03 June 2019

**References**

1. Beharry, A. A. & Woolley, G. A. Azobenzene photoswitches for biomolecules. *Chem. Soc. Rev.* **40**, 4422–4437 (2011).

2. Szymanski, W., Beierle, J. M., Kistemaker, H. A. V., Velema, W. A. & Feringa, B. L. Reversible photocontrol of biological systems by the incorporation of molecular photoswitches. *Chem. Rev.* **113**, 6114–6178 (2013).

3. Bléger, D. & Hecht, S. Visible-light-activated molecular switches. *Angew. Chem. Int. Ed.* **54**, 11338–11349 (2015).

4. Harris, J. D., Moran, M. J. & Aprahamian, I. New molecular switch architectures. *Proc. Natl. Acad. Sci.* **115**, 9414–9422 (2018).

5. Hoorens, M. W. H. & Szymanski, W. Reversible, spatial and temporal control over protein activity using light. *Trends Biochem. Sci.* **43**, 567–575 (2018).

6. Hüll, K., Morstein, J. & Trauner, D. In vivo photopharmacology. *Chem. Rev.* **118**, 10710–10747 (2018).

7. Yu, J. J. et al. Photo-powered stretchable nano-containers based on well-defined vesicles formed by an overcrowded alkene switch. *Chem. Commun.* **52**, 12056–12059 (2016).

8. Senthilkumar, T. et al. Conjugated polymer nanoparticles appending photo-responsive units for controlled drug delivery, release and imaging. *Angew. Chemie Int. Ed.* **57**, 13114–13119. (2018).

9. Lee, I. N. et al. Photoresponsive hydrogels with photoswitchable mechanical properties allow time-resolved analysis of cellular responses to matrix stiffening. *ACS Appl. Mater. Interfaces* **10**, 7765–7776 (2018).

10. Mart, R. J. & Alemann, R. K. Azobenzene photocontrol of peptides and proteins. *Chem. Commun.* **52**, 12262–12277 (2016).

11. Lubbe, A. S. et al. Photoswitching of DNA hybridization using a molecular motor. *J. Am. Chem. Soc.* **140**, 5069–5076 (2018).

12. Zhang, X. et al. Highly photostable, reversibly photoswitchable fluorescent protein with high contrast ratio for live-cell superresolution microscopy. *Proc. Natl Acad. Sci. USA* **113**, 10364–10369 (2016).

13. Laptenok, S. P. et al. Infrared spectroscopy reveals multi-step multi-timescale photoactivation in the photoconvertible protein archetype dronpa. *Nat. Chem.* **10**, 845–852 (2018).

14. Tochitsky, I., Kienzler, M. A., Isacoff, E. & Kramer, R. H. Restoring vision to the blind with chemical photoswitches. *Chem. Rev.* **118**, 10748–10773 (2018).

15. Velema, W. A., Szymanski, W. & Feringa, B. L. Photopharmacology: beyond proof of principle. *J. Am. Chem. Soc.* **136**, 2178–2191 (2014).

16. Broichhagen, J., Frank, J. A. & Trauner, D. A roadmap to success in photopharmacology. *Acc. Chem. Res.* **48**, 1947–1960 (2015).

17. Lerch, M. M. et al. Emerging targets in photopharmacology. *Angew. Chem. Int. Ed.* **55**, 10978–10999 (2016).

18. Kneuttinger, A. C. et al. Artificial light regulation of an allosteric enzyyme complex by a photosensitive ligand. *ChemBioChem* **19**, 1750–1757 (2018).

19. Lerch, M. M., Szymański, W. & Feringa, B. L. The (photo)chemistry of Stenhouse photoswitches: guiding principles and system design. *Chem. Soc. Rev.* **47**, 1910–1937 (2018).

20. Huang, C. Y. et al. N,N'-disubstituted indigos as readily available red-light photoswitches with tunable thermal half-Lives. *J. Am. Chem. Soc.* **139**, 15205–15211 (2017).

21. Wiedbrauk, S. & Dube, H. Hemithioindigo—an emerging photoswitch. *Tetrahedron Lett.* **56**, 4266–4274 (2015).

22. Petermayer, C. & Dube, H. Indigoid photoswitches: visible light responsive molecular tools. *Acc. Chem. Res.* **51**, 1153–1163 (2018).

23. Sadovski, O., Beharry, A. A., Zhang, F. & Woolley, G. A. Spectral tuning of azobenzene photoswitches for biological applications. *Angew. Chem. Int. Ed.* **48**, 1484–1486 (2009).

24. Dong, M. et al. Near-infrared photoswitching of azobenzenes under physiological conditions. *J. Am. Chem. Soc.* **139**, 13483–13486 (2017).

25. Wegener, M., Hansen, M. J., Driessen, A. J. M., Szymanski, W. & Feringa, B. L. Photocontrol of antibacterial activity: shifting from UV to red light activation. *J. Am. Chem. Soc.* **139**, 17979–17986 (2017).

26. Passlick, S., Richers, M. & Ellis-Davies, G. C. R. Thermodynamically stable, photoresversible pharmacology in neurons with one- and two-photon excitation. *Angew. Chem. Int. Ed.* **57**, 12554–12557 (2018).

27. Greb, L. & Lehn, J. M. Light-driven molecular motors: Imines as four-step or two-step unidirectional rotors. *J. Am. Chem. Soc.* **136**, 13114–13117 (2014).

28. Greb, L., Eichhöfer, A. & Lehn, J. M. Synthetic molecular motors: thermal N inversion and directional photoinduced C=N bond rotation of camphorquinone imines. *Angew. Chem. - Int. Ed.* **54**, 14345–14348 (2015).

29. Van Dijken, D. J., Kovafićek, P., Ihrig, S. P. & Hecht, S. Acylhydrazones as widely tunable photoswitches. *J. Am. Chem. Soc.* **137**, 14982–14991 (2015).

30. Qian, H., Pramanik, S. & Aprahamian, I. Photochromic hydrazone switches with extremely long thermal half-lives. *J. Am. Chem. Soc.* **139**, 9140–9143 (2017).

31. Li, Q., Qian, H., Shao, B., Hughes, R. P. & Aprahamian, I. Building strain with large macrocycles and using it to tune the thermal half-lives of hydrazone photochromes. *J. Am. Chem. Soc.* **140**, 11829–11835 (2018).

32. Bezdrik, A., Friedländer, P. & Koeniger, P. Über einige derivate des thionaphthens. *Chem. Ber.* **41**, 227–242 (1908).

33. Pummerer, R. Über isatin-anile, derivate des thionaphtenchinons. *Chem. Ber.* **43**, 1370–1376 (1910).

34. Soeta, T., Shitaya, S., Okuno, T., Fujinami, S. & Ukaji, Y. Efficient synthesis of benzothiophenes by [4+1] cycloaddition of 2-mercaptobenzaldehyde derivatives with isocyanides. *Tetrahedron* **72**, 7901–7905 (2016).

35. Wiedbrauk, S. et al. Twisted hemithioindigo photoswitches: solvent polarity determines the type of light-induced rotations. *J. Am. Chem. Soc.* **138**, 12219–12227 (2016).

36. Zhao, Y. & Truhlar, D. G. The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other function. *Theor. Chem. Acc.* **120**, 215–241 (2008).

37. Ditchfield, R., Hehe, W. J. & Pople, J. A. Self-consistent molecular-orbital methods. ix. an extended gaussian-type basis for molecular-orbital studies of organic molecules. *J. Chem. Phys.* **54**, 724–728 (1971).

38. Marenich, A. V., Cramer, C. J. & Truhlar, D. G. Unvierval solvation mode based on solute electron density and a contiuum model of the solvent defind by the bulk dielectric constant and atomic surface tensions. *J. Phys. Chem. B.* **113**, 6378–6396 (2009).

39. Zweig, J. E. & Newhouse, T. R. Isomer-specific hydrogen bonding as a design principle for bidirectionally quantitative and redshifted hemithioindigo photoswitches. *J. Am. Chem. Soc.* **139**, 10956–10959 (2017).

40. Shaabani, A. & Zahedi, M. Semiempirical molecular orbital calculation of azobenzene: Stability study of isomers and mechanism of E/Z isomerization. *J. Mol. Struct.* **506**, 257–261 (2000).

41. Lehn, J. M. Conjecture: imines as unidirectional photodriven molecular motors-motional and constitutional dynamic devices. *Chem. A Eur. J.* **12**, 5910–5915 (2006).

42. Lehn, J. M. Nitrogen inversion. *Dyn. Stereochem. Fortschr. der Chem. Forsch.* **15**, 311–377 (1970).

43. Kitzig, S., Thilemann, M., Cordes, T. & Rück-braun, K. Light-switchable peptides with a hemithioindigo unit: peptide design, photochromism, and optical Spectrosc. Chem. Phys. **17**, 1252–1263 (2016).

44. Otolski, C. J., Raj, A. M., Ramamurthy, V. & Elles, G. G. Ultrafast dynamics of encapsulated molecules reveals new insights on the photoisomerization mechanism for azobenzenes. *J. Phys. Chem. Lett.* **10**, 121–127 (2019).

45. Johnson, P. J. M. et al. Local vibrational coherences drive the primary photochemistry of vision. *Nat. Chem.* **7**, 980–986 (2015).
46. Nenov, A. et al. UV-light-induced vibrational coherences; the key to understand kasha rulo violation in trans-azobenzene. *J. Phys. Chem. Let.* **9**, 1534–1541 (2018).

47. Maerz, B. et al. Making fast photoswitches faster—using hammett analysis to understand the limit of donor-acceptor approaches for faster hemithioindigo photoswitches. *Chem. Eur. J.* **20**, 13984–13992 (2014).

48. Sailer, A. et al., Hemithioindigos for cellular photopharmacology: desymmetrised molecular switch scaffolds enabling design control over the isomer-dependency of potent antimitotic bioactivity. *Chem. Bio. Chem.* https://doi.org/10.1002/cbic.201800752 (2019).

49. Levine, W. G. Metabolism of AZO dyes: implication for detoxication and activation. *Drug Metab. Rev.* **23**, 253–309 (1991).

50. Kienzler, M. A. et al. A red-shifted, fast-relaxing azobenzene photoswitch for visible light control of an ionotropic glutamate receptor. *J. Am. Chem. Soc.* **135**, 17683–17686 (2013).

51. Izquierdo-Serra, M. et al. Two-photon neuronal and astrocytic stimulation with azobenzene-based photoswitches. *J. Am. Chem. Soc.* **136**, 8693–8701 (2014).

52. Laurent, A. D. & Jacquemin, D. TD-DFT benchmarks: a review. *Int. J. Quantum Chem.* **113**, 2019–2039 (2013).

53. Becke, A. D. Density-functional exchange-energy approximation with correct asymptotic behavior. *Phys. Rev. A* **38**, 3098–3100 (1988).

54. Lee, C., Yang, W. & Parr, G. R. Development of the Colic-Salvetti correlation-energy into a functional of the electron density. *Am. Phys. Soc.* **37**, 785–789 (1988).

55. Jamróz, M. H., Dobrowolski, J. C. & Brzozowski, R. Vibrational modes of 2,6-, 2,7-, and 2,3-diisopropylnaphthalene. A DFT study. *J. Mol. Struct.* **787**, 172–183 (2006).

56. Jamróz, M. H. Vibrational energy distribution analysis (VEDA): Scopes and limitations. *Spectrochim. Acta. A Mol. Biomol. Spectrosc.* **114**, 220–230 (2013).

57. Caricato, M. et al. Formation and relaxation of excited states in solution: a new time dependent polarizable continuum model based on time dependent density functional theory. *J. Chem. Phys.* **124**, 124520-1-13 (2006).

58. Frisch, M. et al. Gaussian09.D01. (Gaussian Inc, Wallingford, 2009).

59. Frisch, M. et al. Gaussian16.A03. (Gaussian Inc, Wallingford, 2009).

60. Cheeseman, J. R., Trucks, G. W., Keith, T. A. & Frisch, M. J. A comparison of models for calculating nuclear magnetic resonance shielding tensors. *J. Chem. Phys.* **104**, 5497–5509 (1996).

61. Henry, E. R. The use of matrix methods in the modeling of spectroscopic data sets. *Biophys. J.* **72**, 652–673 (1997).

62. Van Stokkum, I. H. M., Larsen, D. S. & Van Grondelle, R. Global and target analysis of time-resolved spectra. *Biochim. Biophys. Acta. Bioenerg.* **1657**, 82–104 (2004).

63. Snellenburg, J. J., Laptenok, S. P., Seger, R., Mullen, K. M. & van Stokkum, I. H. M. Glotaran: a java-based graphical user interface for the R package TIMP. *J. Stat. Softw.* **49**, https://doi.org/10.18637/jss.v049.i03 (2012).

64. Feldmeier, C., Bartling, H., Riedle, E. & Gschwind, R. M. LED based NMR illumination device for mechanistic studies on photochemical reactions—versatile and simple, yet surprisingly powerful. *J. Magn. Reson.* **232**, 39–44 (2013).

65. Bini, R., Ballerini, R., Pratesi, G. & Jodl, H. J. Experimental setup for Fourier transform infrared spectroscopy studies in condensed matter at high pressure and low temperatures. *Rev. Sci. Instrum.* **68**, 3154–3160 (1997).

## Acknowledgements

The support of the Netherlands Organization for Scientific Research (NWO-CW VIDI grant 723.014.001 to W.S.) and the European Union Horizon 2020 Research and Innovation

Programme (grant agreement: "Laserlab-Europe", H2020 EC-GA 654148) is kindly acknowledged. M.M. acknowledges the ERDF/ESF project "Nanotechnologies for Future" (CZ.02.1.01/0.0/0.0/16_019/0000754), the Slovak Research and Development Agency (project no. APVV-15-0105) and CMST COST Action CM1405 MOLIM: MOleCules In Motion. This research used resources of (1) the GENCI-CINES/IDRIS (Grants A0020805 17), (2) CCIPL (Centre de Calcul Intensif des Pays de Loire), (3) the HPCC of the Matej Bel University in Banska Bystrica (ITMS 26230120002 and 26210120002 supported by the Research and Development Operational Programme funded by the ERDF). This work was supported financially by the European Research Council (ERC; advanced grant no. 694345 to B.L.F.) and the Ministry of Education, Culture and Science (Gravitation program no. 024.001.035). We thank Pieter van der Meulen for assistance with the in NMR irradiation experiments. M.M. and A.D.L. thank Denis Jacquemin for careful advice and fruitful discussions. We thank Mark Koenis for recording room-temperature IR spectra.

## Author contributions

M.W.H.H and W.S. conceived the project and designed the molecules. M.W.H.H. and L.S. performed the synthesis. Nanosecond TA spectroscopy was performed by M.W.H.H., M.H. and W.J.B., while M.D.D. performed the femtosecond TA spectroscopic experiments. UV-VIS experiments were performed by M.W.H.H. NMR experiments were done by M.W.H.H and W.S.; low-temperature FT-IR experiments were done by S.F. and M.D.D. All calculations were done by M.M. and A.D.L. The research was guided by B.L.F., W.J.B. and W.S. The manuscript was written by M.W.H.H, M.M, A.D.L, M.D.D., B.L.F, W.J.B. and W.S. All authors discussed the results and progress in all stages.

## Additional information

Supplementary Information accompanies this paper at https://doi.org/10.1038/s41467-019-10251-8.

**Competing interests:** The authors declare no competing interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**Journal peer review information:** *Nature Communications* thanks the anonymous reviewer(s) for their contribution to the peer review of this work.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019