Comparative quantum-classical dynamics of natural and synthetic molecular rotors show how vibrational synchronization modulates the photoisomerization quantum efficiency.

We use quantum-classical trajectories to investigate the origin of the different photoisomerization quantum efficiency observed in the dim-light visual pigment Rhodopsin and in the light-driven biomimetic molecular rotor para-methoxy N-methyl indanylidene-pyrrolinium (MeO-NAIP) in methanol. Our results reveal that effective light-energy conversion requires, in general, an auxiliary molecular vibration (called promoter) that does not correspond to the rotary motion but synchronizes with it at specific times. They also reveal that Nature has designed Rhodopsin to exploit two mechanisms working in a vibrationally coherent regime. The first uses a wag promoter to ensure that ca. 75% of the absorbed photons lead to unidirectional rotations. The second mechanism ensures that the same process is fast enough to avoid directional randomization. It is found that MeO-NAIP in methanol is incapable of exploiting the above mechanisms resulting into a 50% quantum efficiency loss. However, when the solvent is removed, MeO-NAIP rotation is predicted to synchronize with a ring-inversion promoter leading to a 30% increase in quantum efficiency and, therefore, biomimetic behavior.

Picosecond quantum-classical dynamics reveals that the coexistence of light-induced microbial and animal chromophore rotary motion modulates the isomerization quantum yield of heliorhodopsin.

Rhodopsins are light-responsive proteins forming two vast and evolutionary distinct superfamilies whose functions are invariably triggered by the photoisomerization of a single retinal chromophore. In 2018 a third widespread superfamily of rhodopsins called heliorhodopsins was discovered using functional metagenomics. Heliorhodopsins, with their markedly different structural features with respect to the animal and microbial superfamilies, offer an opportunity to study how evolution has manipulated the chromophore photoisomerization to achieve adaptation. One question is related to the mechanism of such a reaction and how it differs from that of animal and microbial rhodopsins. To address this question, we use hundreds of quantum-classical trajectories to simulate the spectroscopically documented picosecond light-induced dynamics of a heliorhodopsin from the archaea thermoplasmatales archaeon (TaHeR). We show that, consistently with the observations, the trajectories reveal two excited state decay channels. However, inconsistently with previous hypotheses, only one channel is associated with the -C13C14- rotation of microbial rhodopsins while the second channel is characterized by the -C11C12- rotation typical of animal rhodopsins. The fact that such -C11C12- rotation is aborted upon decay and ground state relaxation, explains why illumination of TaHeR only produces the 13-cis isomer with a low quantum efficiency. We argue that the documented lack of regioselectivity in double-bond excited state twisting motion is the result of an "adaptation" that could be completely lost via specific residue substitutions modulating the steric hindrance experienced along the isomerization motion.

Resonant multiphoton processes and excitation limits to structural dynamics.

Understanding the chemical reactions that give rise to functional biological systems is at the core of structural biology. As techniques are developed to study the chemical reactions that drive biological processes, it must be ensured that the reaction occurring is indeed a biologically relevant pathway. There is mounting evidence indicating that there has been a propagation of systematic error in the study of photoactive biological processes; the optical methods used to probe the structural dynamics of light activated protein functions have failed to ensure that the photoexcitation prepares a well-defined initial state relevant to the biological process of interest. Photoexcitation in nature occurs in the linear (one-photon per chromophore) regime; however, the extreme excitation conditions used experimentally give rise to biologically irrelevant multiphoton absorption. To evaluate and ensure the biological relevance of past and future experiments, a theoretical framework has been developed to determine the excitation conditions, which lead to resonant multiphoton absorption (RMPA) and thus define the excitation limit in general for the study of structural dynamics within the 1-photon excitation regime. Here, we apply the theoretical model to bacteriorhodopsin (bR) and show that RMPA occurs when excitation conditions exceed the linear saturation threshold, well below typical excitation conditions used in this class of experiments. This work provides the guidelines to ensure excitation in the linear 1-photon regime is relevant to biological and chemical processes.

Comparative Study of Uracil Excited-State Photophysics in Water and Acetonitrile via RMS-CASPT2-Driven Quantum-Classical Trajectories.

We present a nonadiabatic molecular dynamics study of the ultrafast processes occurring in uracil upon UV light absorption, leading to electronic excitation and subsequent nonradiative decay. Previous studies have indicated that the mechanistic details of this process are drastically different depending on whether the process takes place in the gas phase, acetonitrile, or water. However, such results have been produced using quantum chemical methods that did not incorporate both static and dynamic electron correlation. In order to assess the previously proposed mechanisms, we simulate the photodynamics of uracil in the three environments mentioned above using quantum-classical trajectories and, for solvated uracil, hybrid quantum mechanics/molecular mechanics (QM/MM) models driven by the rotated multistate complete active space second-order perturbation (RMS-CASPT2) method. To do so, we exploit the gradient recently made available in OpenMolcas and compare the results to those obtained using the complete active space self-consistent field (CASSCF) method only accounting for static electron correlation. We show that RMS-CASPT2 produces, in general, a mechanistic picture different from the one obtained at the CASSCF level but confirms the hypothesis advanced on the basis of previous ROKS and TDDFT studies thus highlighting the importance of incorporating dynamic electron correlation in the investigation of ultrafast electronic deactivation processes.

Assessment of the Electron Correlation Treatment on the Quantum-Classical Dynamics of Retinal Protonated Schiff Base Models: XMS-CASPT2, RMS-CASPT2, and REKS Methods.

We compare the performance of three different multiconfigurational wave function-based electronic structure methods and two implementations of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method. The study is characterized by three features: (i) it uses a small set of quantum-classical trajectories rather than potential energy surface mapping, (ii) it focuses, exclusively, on the photoisomerization of retinal protonated Schiff base models, and (iii) it probes the effect of both methyl substitution and the increase in length of the conjugate π-system. For each tested method, the corresponding analytical gradients are used to drive the quantum-classical (Tully's FSSH method) trajectory propagation, including the recent multistate XMS-CASPT2 and RMS-CASPT2 gradients. It is shown that while CASSCF, XMS-CASPT2, and RMS-CASPT2 yield consistent photoisomerization dynamics descriptions, REKS produces, in some of these systems, qualitatively different behavior that is attributed to a flatter and topographically different excited state potential energy surface. The origin of this behavior can be traced back to the effect of the employed density functional approximation. The above studies are further expanded by benchmarking, at the CASSCF and REKS levels, the electronic structure methods using a QM/MM model of the visual pigment rhodopsin.

Automatic Rhodopsin Modeling with Multiple Protonation Microstates.

Automatic Rhodopsin Modeling (ARM) is a simulation protocol providing QM/MM models of rhodopsins capable of reproducing experimental electronic absorption and emission trends. Currently, ARM is restricted to a single protonation microstate for each rhodopsin model. Herein, we incorporate an extension of the minimal electrostatic model (MEM) into the ARM protocol to account for all relevant protonation microstates at a given pH. The new ARM+MEM protocol determines the most important microstates contributing to the description of the absorption spectrum. As a test case, we have applied this methodology to simulate the pH-dependent absorption spectrum of a toy model, showing that the single-microstate picture breaks down at certain pH values. Subsequently, we applied ARM+MEM toAnabaenasensory rhodopsin, confirming an improved description of its absorption spectrum when the titration of several key residues is considered.

Experimental Assessment of the Electronic and Geometrical Structure of a Near-Infrared Absorbing and Highly Fluorescent Microbial Rhodopsin.

The recently discovered Neorhodopsin (NeoR) exhibits absorption and emission maxima in the near-infrared spectral region, which together with the high fluorescence quantum yield makes it an attractive retinal protein for optogenetic applications. The unique optical properties can be rationalized by a theoretical model that predicts a high charge transfer character in the electronic ground state (S0) which is otherwise typical of the excited state S1 in canonical retinal proteins. The present study sets out to assess the electronic structure of the NeoR chromophore by resonance Raman (RR) spectroscopy since frequencies and relative intensities of RR bands are controlled by the ground and excited state's properties. The RR spectra of NeoR differ dramatically from those of canonical rhodopsins but can be reliably reproduced by the calculations carried out within two different structural models. The remarkable agreement between the experimental and calculated spectra confirms the consistency and robustness of the theoretical approach.

Correction: From a one-mode to a multi-mode understanding of conical intersection mediated ultrafast organic photochemical reactions.

Correction for 'From a one-mode to a multi-mode understanding of conical intersection mediated ultrafast organic photochemical reactions' by Yorrick Boeije et al., Chem. Soc. Rev., 2023, 52, 2643-2687, https://doi.org/10.1039/D2CS00719C.

From a one-mode to a multi-mode understanding of conical intersection mediated ultrafast organic photochemical reactions.

Over the last few decades, conical intersections (CoIns) have grown from theoretical curiosities into common mechanistic features of photochemical reactions, whose function is to funnel electronically excited molecules back to their ground state in regions where the potential energy surfaces (PESs) of two electronic states become degenerate. Analogous to transition states in thermal chemistry, CoIns appear as transient structures providing a kinetic bottleneck along a reaction coordinate. However, such a bottleneck is not associated with the probability of crossing an energy barrier but rather with an excited state decay probability along a full "line" of transient structures connected by non-reactive modes, the intersection space (IS). This article will review our understanding of the factors controlling CoIn mediated ultrafast photochemical reactions, taking a physical organic chemist approach by discussing a number of case studies for small organic molecules and photoactive proteins. Such discussion will be carried out by first introducing the "standard" one-mode model based on Landau-Zener (LZ) theory to describe a reactive excited state decay event intercepting, locally, a single CoIn along a single direction, and then by providing a modern perspective based on the effects of the phase matching of multiple modes on the same local event, thus redefining and expanding the description of the excited state reaction coordinate. The direct proportionality between the slope (or velocity) along one mode and decay probability at a single CoIn is a widely applied fundamental principle that follows from the LZ model, yet it fails to provide a complete understanding of photochemical reactions whose local reaction coordinate changes along the IS. We show that in these situations, in particular by focussing on rhodopsin double bond photoisomerization, it is mandatory to consider additional molecular modes and their phase relationship approaching the IS, hence providing a key mechanistic principle of ultrafast photochemistry based on the phase matching of those modes. We anticipate that this qualitative mechanistic principle should be considered in the rational design of any ultrafast excited state process, impacting various fields of research ranging from photobiology to light-driven molecular devices.

Automated QM/MM Screening of Rhodopsin Variants with Enhanced Fluorescence.

We present a computational protocol for the fast and automated screening of excited-state hybrid quantum mechanics/molecular mechanics (QM/MM) models of rhodopsins to be used as fluorescent probes based on the automatic rhodopsin modeling protocol (a-ARM). Such "a-ARM fluorescence screening protocol" is implemented through a general Python-based driver, PyARM, that is also proposed here. The implementation and performance of the protocol are benchmarked using different sets of rhodopsin variants whose absorption and, more relevantly, emission spectra have been experimentally measured. We show that, despite important limitations that make unsafe to use it as a black-box tool, the protocol reproduces the observed trends in fluorescence and it is capable of selecting novel potentially fluorescent rhodopsins. We also show that the protocol can be used in mechanistic investigations to discern fluorescence enhancement effects associated with a near degeneracy of the S1/S2 states or, alternatively, with a barrier generated via coupling of the S0/S1 wave functions.

Retinal chromophore charge delocalization and confinement explain the extreme photophysics of Neorhodopsin.

The understanding of how the rhodopsin sequence can be modified to exactly modulate the spectroscopic properties of its retinal chromophore, is a prerequisite for the rational design of more effective optogenetic tools. One key problem is that of establishing the rules to be satisfied for achieving highly fluorescent rhodopsins with a near infrared absorption. In the present paper we use multi-configurational quantum chemistry to construct a computer model of a recently discovered natural rhodopsin, Neorhodopsin, displaying exactly such properties. We show that the model, that successfully replicates the relevant experimental observables, unveils a geometrical and electronic structure of the chromophore featuring a highly diffuse charge distribution along its conjugated chain. The same model reveals that a charge confinement process occurring along the chromophore excited state isomerization coordinate, is the primary cause of the observed fluorescence enhancement.

On the fluorescence enhancement of arch neuronal optogenetic reporters.

The lack of a theory capable of connecting the amino acid sequence of a light-absorbing protein with its fluorescence brightness is hampering the development of tools for understanding neuronal communications. Here we demonstrate that a theory can be established by constructing quantum chemical models of a set of Archaerhodopsin reporters in their electronically excited state. We found that the experimentally observed increase in fluorescence quantum yield is proportional to the computed decrease in energy difference between the fluorescent state and a nearby photoisomerization channel leading to an exotic diradical of the protein chromophore. This finding will ultimately support the development of technologies for searching novel fluorescent rhodopsin variants and unveil electrostatic changes that make light emission brighter and brighter.

Towards the engineering of a photon-only two-stroke rotary molecular motor.

The rational engineering of photoresponsive materials, e.g., light-driven molecular motors, is a challenging task. Here, we use structure-related design rules to prepare a prototype molecular rotary motor capable of completing an entire revolution using, exclusively, the sequential absorption of two photons; i.e., a photon-only two-stroke motor. The mechanism of rotation is then characterised using a combination of non-adiabatic dynamics simulations and transient absorption spectroscopy measurements. The results show that the rotor moiety rotates axially relative to the stator and produces, within a few picoseconds at ambient T, an intermediate with the same helicity as the starting structure. We discuss how such properties, that include a 0.25 quantum efficiency, can help overcome the operational limitations of the classical overcrowded alkene designs.

Torsionally broken symmetry assists infrared excitation of biomimetic charge-coupled nuclear motions in the electronic ground state.

The concerted interplay between reactive nuclear and electronic motions in molecules actuates chemistry. Here, we demonstrate that out-of-plane torsional deformation and vibrational excitation of stretching motions in the electronic ground state modulate the charge-density distribution in a donor-bridge-acceptor molecule in solution. The vibrationally-induced change, visualised by transient absorption spectroscopy with a mid-infrared pump and a visible probe, is mechanistically resolved by ab initio molecular dynamics simulations. Mapping the potential energy landscape attributes the observed charge-coupled coherent nuclear motions to the population of the initial segment of a double-bond isomerization channel, also seen in biological molecules. Our results illustrate the pivotal role of pre-twisted molecular geometries in enhancing the transfer of vibrational energy to specific molecular modes, prior to thermal redistribution. This motivates the search for synthetic strategies towards achieving potentially new infrared-mediated chemistry.

Design, synthesis and biological evaluation of light-driven on-off multitarget AChE and MAO-B inhibitors.

Neurodegenerative diseases are multifactorial disorders characterized by protein misfolding, oxidative stress, and neuroinflammation, finally resulting in neuronal loss and cognitive dysfunctions. Nowadays, an attractive strategy to improve the classical treatments is the development of multitarget-directed molecules able to synergistically interact with different enzymes and/or receptors. In addition, an interesting tool to refine personalized therapies may arise from the use of bioactive species able to modify their activity as a result of light irradiation. To this aim, we designed and synthesized a small library of cinnamic acid-inspired isomeric compounds with light modulated activity able to inhibit acetylcholinesterase (AChE) and monoamine oxidase B (MAO-B), with remarkable selectivity over butyrylcholinesterase (BChE) and MAO-A, which have been investigated as the enzyme targets related to Alzheimer's disease (AD). The inhibitory activities were evaluated for the pure E-diastereomers and the E/Z-diastereomer mixtures, obtained upon UV irradiation. Molecular docking studies were carried out to rationalize the differences in the inhibition potency of the E and Z diastereomers of the best performing analogue 1c. Our preliminary findings may open-up the way for developing innovative multitarget photo-switch drugs against neurodegenerative diseases.

Design, Synthesis and Characterization of a Visible-Light-Sensitive Molecular Switch and Its PEGylation Towards a Self-Assembling Molecule.

HBDI-like chromophores represent a novel set of biomimetic switches mimicking the fluorophore of the green fluorescent protein that are currently studied with the hope to expand the molecular switch/motor toolbox. However, until now members capable of absorbing visible light in their neutral (i. e. non-anionic) form have not been reported. In this contribution we report the preparation of an HBDI-like chromophore based on a 3-phenylbenzofulvene scaffold capable of absorbing blue light and photoisomerizing on the picosecond timescale. More specifically, we show that double-bond photoisomerization occurs in both the E-to-Z and Z-to-E directions and that these can be controlled by irradiating with blue and UV light, respectively. Finally, as a preliminary applicative result, we report the incorporation of the chromophore in an amphiphilic molecule and demonstrate the formation of a visible-light-sensitive nanoaggregated state in water.

Quantum-classical simulations of rhodopsin reveal excited-state population splitting and its effects on quantum efficiency.

The activation of rhodopsin, the light-sensitive G-protein-coupled receptor responsible for dim-light vision in vertebrates, is driven by an ultrafast excited-state double-bond isomerization with a quantum efficiency of almost 70%. The origin of such light sensitivity is not understood and a key question is whether in-phase nuclear motion controls the quantum efficiency value. In this study we used hundreds of quantum-classical trajectories to show that, 15 fs after light absorption, a degeneracy between the reactive excited state and a neighbouring state causes the splitting of the rhodopsin population into subpopulations. These subpopulations propagate with different velocities and lead to distinct contributions to the quantum efficiency. We also show here that such splitting is modulated by protein electrostatics, thus linking amino acid sequence variations to quantum efficiency modulation. Finally, we discuss how such a linkage that in principle could be exploited to achieve higher quantum efficiencies would simultaneously increase the receptor thermal noise leading to a trade-off that may have played a role in rhodopsin evolution.

Evolution of the Automatic Rhodopsin Modeling (ARM) Protocol.

In recent years, photoactive proteins such as rhodopsins have become a common target for cutting-edge research in the field of optogenetics. Alongside wet-lab research, computational methods are also developing rapidly to provide the necessary tools to analyze and rationalize experimental results and, most of all, drive the design of novel systems. The Automatic Rhodopsin Modeling (ARM) protocol is focused on providing exactly the necessary computational tools to study rhodopsins, those being either natural or resulting from mutations. The code has evolved along the years to finally provide results that are reproducible by any user, accurate and reliable so as to replicate experimental trends. Furthermore, the code is efficient in terms of necessary computing resources and time, and scalable in terms of both number of concurrent calculations as well as features. In this review, we will show how the code underlying ARM achieved each of these properties.

Xanthopsin-Like Systems via Site-Specific Click-Functionalization of a Retinoic Acid Binding Protein.

The use of light-responsive proteins to control both living or synthetic cells, is at the core of the expanding fields of optogenetics and synthetic biology. It is thus apparent that a richer reaction toolbox for the preparation of such systems is of fundamental importance. Here, we provide a proof-of-principle demonstration that Morita-Baylis-Hillman adducts can be employed to perform a facile site-specific, irreversible and diastereoselective click-functionalization of a lysine residue buried into a lipophilic binding pocket and yielding an unnatural chromophore with an extended π-system. In doing so we effectively open the path to the in vitro preparation of a library of synthetic proteins structurally reminiscent of xanthopsin eubacterial photoreceptors. We argue that such a library, made of variable unnatural chromophores inserted in an easy-to-mutate and crystallize retinoic acid transporter, significantly expand the scope of the recently introduced rhodopsin mimics as both optogenetic and "lab-on-a-molecule" tools.