Thermal Half-Lives of Azobenzene Derivatives: Virtual Screening Based on Intersystem Crossing Using a Machine Learning Potential

Accessing a Diverse Set of Functional Red-Light Photoswitches by Selective Copper-Catalyzed Indigo N-Arylation

The ability to correlate the structure of a molecule with its properties is the key to the rational and accelerated design of new functional compounds and materials. Taking photoswitches as an example, the thermal stability of the metastable state is a crucial property that dictates their application in molecular systems. Indigos have recently emerged as an attractive motif for designing photoswitchable molecules due to their red-light addressability, which can be advantageous in biomedical and material applications. The lack of synthetic techniques to derivatize the abundant parent dye and a thorough understanding of the impact of structural factors on the photochemical and thermal properties hinder broad applications of this emerging photoswitch class. Herein, we report an efficient copper-catalyzed indigo N-arylation that enables the installation of a wide variety of aryl moieties carrying useful functional groups. The exclusive selectivity for monoarylation likely originates from a bimetallic cooperative mechanism through a binuclear copper-indigo intermediate. Functional N-aryl-N'-alkylindigos were prepared and shown to photoisomerize efficiently under red light. Moreover, this design allows for the modulation of thermal half-lives through N-aryl substituents, while the N'-alkyl groups enable the independent attachment of functional moieties without affecting the photochromic properties. A strong correlation between the structure of the N-aryl moiety and the thermal stability of the photogenerated Z-isomers was achieved by multivariate linear regression models obtained through a data-science workflow. This work thus builds an avenue leading to versatile red-light photoswitches and a general method for structure-property correlation that is expected to be broadly applicable to the design of photoresponsive molecules.

Detour to success: photoswitching via indirect excitation

Photoswitchable molecules that undergo nanoscopic changes upon photoisomerisation can be harnessed to control macroscopic properties such as colour, solubility, shape, and motion of the systems they are incorporated into. These molecules find applications in various fields of chemistry, physics, biology, and materials science. Until recently, research efforts have focused on the design of efficient photoswitches responsive to low-energy (red or near-infrared) irradiation, which however may compromise other molecular properties such as thermal stability and robustness. Indirect isomerisation methods enable photoisomerisation with low-energy photons without altering the photoswitch core, and also open up new avenues in controlling the thermal switching mechanism. In this perspective, we present the state of the art of five indirect excitation methods: two-photon excitation, triplet sensitisation, photon upconversion, photoinduced electron transfer, and indirect thermal methods. Each impacts our understanding of the fundamental physicochemical properties of photochemical switches, and offers unique application prospects in biomedical technologies and beyond.

Photoswitchable luminescent lanthanide complexes controlled and interrogated by four orthogonal wavelengths of light

Optical information storage requires careful control of excitation and emission wavelengths in a reversible and orthogonal manner to enable efficient reading, writing, and erasing of information. Photochromic systems, in which a photoswitch is typcially coupled to an emissive organic fluorophore, have much promise in this regard. However, these suffer from considerable spectral overlap between the switch and fluorophore, such that their emissive and photoswitchable properties are not orthogonal. Here, we overcome this limitation by coupling visible/NIR emissive lanthanide complexes with molecular photoswitches, enabling reversible and orthogonal photoswitching with visible light. Crucially, photoswitching does not lead to sensitised emission from the lanthanide, while excitation of the lanthanide does not induce photoswitching, enabling the state of the system to be probed without perturbation of the switch. This opens up the possibility of developing multi-colour read-write methods for information storage using emissive photoswitches.

Tuning properties of functionalized azobenzene photoswitches through variation of substitution and their position

This perspective article features the recent work published in this journal on functionalized azobenzene photoswitches, studied by Venkataramani and coworkers [Grewal, S. et al (2024) Photochem. Photobiol., 1-16]. Azobenzene-based photochromic molecules undergo reversible trans-to-cis isomerization under exposure to light. The reverse cis-to-trans isomerization also takes place under thermal conditions in the dark. The study explores how different functional groups at various positions affect the responses of these compounds to light, their thermal stabilities, the isomeric ratios in the photostationary states, and their supramolecular behavior. Results show varied thermal half-lives, influenced by the terminal groups and their positions. Aggregation studies reveal the formation of supramolecular architectures, from microcrystals to gels with these systems. Despite complexities with the large volume of work with 30 different compounds, the research provides insights into designing efficient photochromic supramolecular motifs.

The effect of substituent position and solvent on thermal Z‒E isomerization of dihydroquinolylazotetrazole dyes: kinetic, thermodynamic, and spectral approaches

Kinetic and thermodynamic parameters have been investigated for the thermal Z‒E isomerization of dihydroquinolylazotetrazole dyes with alkyl substituents (Me, t-Bu, and Adm) at positions 1 (dyes 2) and 2 (dyes 3) of the tetrazole moiety in two solvents of different polarity, acetonitrile (MeCN) and toluene. The experimental results show crucial dependence of these parameters on a substituent position in the tetrazole moiety and on a solvent. For dyes 2, Eact and ΔH‡ are lower in MeCN than in toluene that results in a high increase in the lifetimes of the Z isomers: from milliseconds in MeCN to minutes in toluene. For dyes 3, the difference in Eact and ΔH‡ in the two solvents is opposite: Eact and ΔH‡ are by more than 20 kJ mol-1 higher in MeCN, nevertheless, the rate constants for 3 in toluene are comparable with those in MeCN at the ambient temperature and the difference in the behavior is determined by the value of negative entropy of activation. Quantum-chemical calculations of the thermal Z‒E isomerization show the possibility of the process to occur via crossing from the S0 to the thermally induced T1 state. The contribution of this path is highest for 3 in toluene. The analysis of the absorption spectra demonstrates that for the E isomers, the n‒π* and π‒π* transitions are within the long-wavelength absorption band and their positions relative each other are opposite in the solvents: the n‒π* transition is blue-shifted relative to the π‒π* transition in MeCN and is red-shifted in toluene.

Generative organic electronic molecular design informed by quantum chemistry

Generative molecular design strategies have emerged as promising alternatives to trial-and-error approaches for exploring and optimizing within large chemical spaces. To date, generative models with reinforcement learning approaches have frequently used low-cost methods to evaluate the quality of the generated molecules, enabling many loops through the generative model. However, for functional molecular materials tasks, such low-cost methods are either not available or would require the generation of large amounts of training data to train surrogate machine learning models. In this work, we develop a framework that connects the REINVENT reinforcement learning framework with excited state quantum chemistry calculations to discover molecules with specified molecular excited state energy levels, specifically molecules with excited state landscapes that would serve as promising singlet fission or triplet-triplet annihilation materials. We employ a two-step curriculum strategy to first find a set of diverse promising molecules, then demonstrate the framework's ability to exploit a more focused chemical space with anthracene derivatives. Under this protocol, we show that the framework can find desired molecules and improve Pareto fronts for targeted properties versus synthesizability. Moreover, we are able to find several different design principles used by chemists for the design of singlet fission and triplet-triplet annihilation molecules.

Singlet and Triplet Pathways Determine the Thermal Z/E Isomerization of an Arylazopyrazole-Based Photoswitch

Understanding the thermal isomerization mechanism of azobenzene derivatives is essential to designing photoswitches with tunable half-lives. Herein, we employ quantum chemical calculations, nonadiabatic transition state theory, and photosensitized experiments to unravel the thermal Z/E isomerization of a heteroaromatic azoswitch, the phenylazo-1,3,5-trimethylpyrazole. In contrast to the parent azobenzene, we predict two pathways to be operative at room temperature. One is a conventional ground-state reaction occurring via inversion of the aryl group, and the other is a nonadiabatic process involving intersystem crossing to the lowest-lying triplet state and back to the ground state, accompanied by a torsional motion around the azo bond. Our results illustrate that the fastest reaction rate is not controlled by the mechanism involving the lowest activation energy, but the size of the spin-orbit couplings at the crossing between the singlet and the triplet potential energy surfaces is also determinant. It is therefore mandatory to consider all of the multiple reaction pathways in azoswitches in order to predict experimental half-lives.

On the Simulation of Thermal Isomerization of Molecular Photoswitches in Biological Systems

Molecular photoswitches offer precise, reversible photocontrol over biomolecular functions and are promising light-regulated drug candidates with minimal side effects. Quantifying thermal isomerization rates of photoswitches in their target biomolecules is essential for fine-tuning their light-controlled drug activity. However, the effects of protein binding on isomerization kinetics remain poorly understood, and simulations are crucial for filling this gap. Challenges in the simulation include describing multireference electronic structures near transition states, disentangling competing reaction pathways, and sampling protein-ligand interactions. To overcome these challenges, we used multiscale simulations to characterize the thermal isomerization of photostatins (PSTs), which are light-regulated microtubule inhibitors for potential cancer phototherapy. We employed a new ab initio multireference electronic structure method in a quantum mechanics/molecular mechanics setting and combined it with enhanced sampling techniques to characterize the cis to trans free-energy profiles of three PSTs in a vacuum, aqueous solution, and tubulin dimer. The significant advantage of our novel approach is the efficient treatment of the multireference character in PSTs' electronic wavefunction throughout the conformational sampling of protein-ligand interactions along their isomerization pathways. We also benchmarked our calculations using high-level ab initio multireference electronic structure methods and explored the competing isomerization pathways. Notably, calculations in a vacuum and implicit solvent models cannot predict the order of the PSTs' thermal half-lives in the aqueous solution observed in the experiment. Only by explicitly treating the solvent molecules can the correct order of isomerization kinetics be reproduced. Protein binding perturbs free-energy barriers due to hydrogen bonding between PSTs and nearby polar residues. Our work generates comprehensive, high-quality benchmark data and offers guidance for selecting computational methods to study the thermal isomerization of photoswitches. Ab initio multireference free-energy calculations in explicit molecular environments are crucial for predicting the effects of substituents on the thermal half-lives of photoswitches in biological systems.

Mapping the Space of Photoswitchable Ligands and Photodruggable Proteins with Computational Modeling

Light-activated drugs are a promising way to localize biological activity and minimize side effects. However, their development is complicated by the numerous photophysical and biological properties that must be simultaneously optimized. To accelerate the design of photoactive drugs, we describe a procedure that combines ligand-protein docking with chemical property prediction based on machine learning (ML). We apply this procedure to 58 proteins and 9000 photo-drug candidates based on azobenzene cis-trans isomerism. We find that most proteins display a preference for trans isomers over cis and that the binding affinities of nominally active/inactive pairs are in fact highly correlated. These findings have significant value for photopharmacology research, and reinforce the need for virtual screening to identify compounds with rare desirable properties. Further, we combine our procedure with quantum chemical validation to identify promising candidates for the photoactive inhibition of PARP1, an enzyme that is over-expressed in cancer cells. The top compounds are predicted to have long-lived active forms, differential bioactivity, and absorption in the near-infrared therapeutic window.

Harnessing Data Augmentation and Normalization Preprocessing to Improve the Performance of Chemical Reaction Predictions of Data-Driven Model

As a template-free, data-driven methodology, the molecular transformer model provides an alternative by which to predict the outcome of chemical reactions and design the route of the retrosynthetic plane in the field of organic synthesis and polymer chemistry. However, in consideration of the small datasets of chemical reactions, the data-driven model suffers from the difficulty of low accuracy in the prediction tasks of chemical reactions. In this contribution, we integrate the molecular transformer model with the strategies of data augmentation and normalization preprocessing to accomplish the three tasks of chemical reactions, including the forward predictions of chemical reactions, and single-step retrosynthetic predictions with and without the reaction classes. It is clearly demonstrated that the prediction accuracy of the molecular transformer model can be significantly raised by the use of proposed strategies for the three tasks of chemical reactions. Notably, after the introduction of the 40-level data augmentation and normalization preprocessing, the top-1 accuracy of the forward prediction increases markedly from 71.6% to 84.2% and the top-1 accuracy of the single-step retrosynthetic prediction with additional reaction class increases from 53.2% to 63.4%. Furthermore, it is found that the superior performance of the data-driven model originates from the correction of the grammatical errors of the SMILES strings, especially for the case of the reaction classes with small datasets.

Solving the Azobenzene Entropy Puzzle: Direct Evidence for Multi-State Reactivity

A solution to the azobenzene "entropy puzzle" [ J. Phys.: Condens. Matter, 2017, 29, 314002] is provided. Previous computational studies of the thermal Z → E (back-)isomerization of azobenzene could not describe the experimentally observed large negative activation entropies. Here it is shown that the experimental results are only compatible with a more complicated multistate rotation mechanism that involves a triplet excited state. Using nonadiabatic transition state theory, close to perfect agreement is achieved between all calculated and experimental Eyring parameters. We also provide new experiments that indicate the presence of a noticeable external heavy-atom effect, which is a direct result of spin-orbit coupling effects being important in the proposed mechanism. These results suggest a reexamination of the mechanisms of related thermal double bond isomerizations in other systems in cases when an excited state of triplet (or other) multiplicity becomes thermally accessible during a rotation process.