Nonadiabatic Photodynamics of Retinal Protonated Schiff Base in Channelrhodopsin 2

Formulation of transition dipole gradients for non-adiabatic dynamics with polaritonic states

A general formulation of the strong coupling between photons confined in a cavity and molecular electronic states is developed for the state-interaction state-average spin-restricted ensemble-referenced Kohn-Sham method. The light-matter interaction is included in the Jaynes-Cummings model, which requires the derivation and implementation of the analytical derivatives of the transition dipole moments between the molecular electronic states. The developed formalism is tested in the simulations of the nonadiabatic dynamics in the polaritonic states resulting from the strong coupling between the cavity photon mode and the ground and excited states of the penta-2,4-dieniminium cation, also known as PSB3. Comparison with the field-free simulations of the excited-state decay dynamics in PSB3 reveals that the light-matter coupling can considerably alter the decay dynamics by increasing the excited state lifetime and hindering photochemically induced torsion about the C=C double bonds of PSB3. The necessity of obtaining analytical transition dipole gradients for the accurate propagation of the dynamics is underlined.

Deciphering Photoreceptors Through Atomistic Modeling from Light Absorption to Conformational Response

In this review, we discuss the successes and challenges of the atomistic modeling of photoreceptors. Throughout our presentation, we integrate explanations of the primary methodological approaches, ranging from quantum mechanical descriptions to classical enhanced sampling methods, all while providing illustrative examples of their practical application to specific systems. To enhance the effectiveness of our analysis, our primary focus has been directed towards the examination of applications across three distinct photoreceptors. These include an example of Blue Light-Using Flavin (BLUF) domains, a bacteriophytochrome, and the orange carotenoid protein (OCP) employed by cyanobacteria for photoprotection. Particular emphasis will be placed on the pivotal role played by the protein matrix in fine-tuning the initial photochemical event within the embedded chromophore. Furthermore, we will investigate how this localized perturbation initiates a cascade of events propagating from the binding pocket throughout the entire protein structure, thanks to the intricate network of interactions between the chromophore and the protein.

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.

Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

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.

Steric and Electronic Origins of Fluorescence in GFP and GFP-like Proteins

Fluorescent proteins have become routine tools for biological imaging. However, their nanosecond lifetimes on the excited state present computational hurdles to a full understanding of these photoactive proteins. In this work, we simulate approximately 0.5 nanoseconds of ab initio molecular dynamics to elucidate steric and electronic features responsible for fluorescent protein behavior. Using green fluorescent protein (GFP) and Dronpa2─widely used fluorescent proteins with contrasting functionality─as case studies, we leverage previous findings in the gas phase and solution to explore the deactivation mechanisms available to these proteins. Starting with ground-state analyses, we identify steric (the distribution of empty pockets near the chromophore) and electronic (electric fields exerted on chromophore moieties) factors that offer potential avenues for rational design. Picosecond timescale simulations on the excited state reveal that the chromophore can access twisted structures in Dronpa2, while the chromophore is largely confined to planarity in GFP. We couple ab initio multiple spawning (AIMS) and enhanced sampling simulations to discover and characterize conical intersection seams that facilitate internal conversion, which is a rare event in both systems. Our AIMS simulations correctly capture the relative fluorescence profiles of GFP and Dronpa2 within the first few picoseconds, and we attribute the diminished fluorescence intensity of Dronpa2, relative to GFP, to flexible chromophore intermediates on the excited state. Furthermore, we predict that twisted chromophore intermediates produce red-shifted intensities in the Dronpa2 fluorescence spectrum. If confirmed experimentally, this spectroscopic signature would provide valuable insights when screening and developing novel fluorescent proteins.

Electronic Structure Methods for the Description of Nonadiabatic Effects and Conical Intersections

Nonadiabatic effects are ubiquitous in photophysics and photochemistry, and therefore, many theoretical developments have been made to properly describe them. Conical intersections are central in nonadiabatic processes, as they promote efficient and ultrafast nonadiabatic transitions between electronic states. A proper theoretical description requires developments in electronic structure and specifically in methods that describe conical intersections between states and nonadiabatic coupling terms. This review focuses on the electronic structure aspects of nonadiabatic processes. We discuss the requirements of electronic structure methods to describe conical intersections and nonadiabatic couplings, how the most common excited state methods perform in describing these effects, and what the recent developments are in expanding the methodology and implementing nonadiabatic couplings.

Structural Factors Determining the Absorption Spectrum of Channelrhodopsins: A Case Study of the Chimera C1C2

Channelrhodopsins are photosensitive proteins that trigger flagella motion in single-cell algae and have been successfully utilized in optogenetic applications. In optogenetics, light is used to activate neural cells in living organisms, which can be achieved by exploiting the ion channel signaling of channelrhodopsins. Tailoring channelrhodopsins for such applications includes the tuning of the absorption maximum. In order to establish rational design and to obtain a desired spectral shift, a basic understanding of the absorption spectrum is required. We have studied the chimera C1C2 as a representative of this protein family and the first member with an available crystal structure. For this purpose, we sampled the conformations of C1C2 using quantum mechanical/molecular mechanical molecular dynamics and subjected the resulting snapshots of the trajectory to excitation energy calculations using ADC(2) and simplified time-dependent density functional theory. In contrast to previous reports, we found that different hydrogen-bonding networks-involving the retinal protonated Schiff base, the putative counterions E162 and D292, and water molecules-had only a small impact on the absorption spectrum. However, in the case of deprotonated E162, increasing the distance to the Schiff base hydrogen-bonding partner led to a systematic blue shift. The β-ionone ring rotation was identified as another important contributor. Yet the most important factors were found to be the bond length alternation and bond order alternation that were linearly correlated to the absorption maximum by up to 62 and 82%, respectively. We ascribe this novel insight into the structural basis of the absorption spectrum to our enhanced protein setup that includes membrane embedding as well as long and extensive sampling.

First-Principles Nonadiabatic Dynamics Simulation of Azobenzene Photodynamics in Solutions

The photoisomerization of azobenzene is a prototypical reaction of various light-activated processes in material and biomedical sciences. However, its reaction mechanism has been under debate for decades, partly due to the challenges in computational simulations to accurately describe the molecule's photodynamics. A recent study (J. Am. Chem. Soc. 2020, 142 (49), 20,680-20,690) addressed the challenges by combining the hole-hole Tamm-Dancoff Approximated (hh-TDA) density functional theory (DFT) method with the ab initio multiple spawning (AIMS) algorithm. The hh-TDA-DFT/AIMS method was applied to first-principles nonadiabatic dynamics simulation of azobenzene's photodynamics in the vacuum. However, it remains necessary to benchmark this new method in realistic molecular environments against experimental data. In the current work, the hh-TDA-DFT/AIMS method was employed in a quantum mechanics/molecular mechanics setting to characterize the trans azobenzene's photodynamics in explicit methanol and n-hexane solvents, following both the S₁ (nπ*) and S₂ (ππ*) excitations. The simulated absorption and fluorescence spectra following the S₂ excitation quantitatively agree with the experiments. However, the hh-TDA-DFT method overestimates the torsional barrier on the S₁ state, leading to an overestimation of the S₁ state lifetime. The excited-state population decays to the ground state through two competing channels. The reactive channel partially yields the cis azobenzene photoproduct, and the unreactive channel exclusively leads to the reactant. The S₂ excitation increases the decay through the unreactive channel and thus decreases the isomerization quantum yield compared to the S₁ excitation. The solvent slows down the azobenzene's torsional dynamics on the S₁ state, but its polarity minimally affects the reaction kinetics and quantum yields. Interestingly, the dynamics of the central torsion and angles of azobenzene play a critical role in determining the final isomer of the azobenzene. This benchmark study validates the hh-TDA-DFT/AIMS method's accuracy for simulating the azobenzene's photodynamics in realistic molecular environments.

Electrostatic Control of Photoisomerization in Channelrhodopsin 2

Channelrhodopsin 2 (ChR2) is the most commonly used tool in optogenetics. Because of its faster photocycle compared to wild-type (WT) ChR2, the E123T mutant of ChR2 is a useful optogenetic tool when fast neuronal stimulation is needed. Interestingly, in spite of its faster photocycle, the initial step of the photocycle in E123T (photoisomerization of retinal protonated Schiff base or RPSB) was found experimentally to be much slower than that of WT ChR2. The E123T mutant replaces the negatively charged E123 residue with a neutral T123 residue, perturbing the electric field around the RPSB. Understanding the RPSB photoisomerization mechanism in ChR2 mutants will provide molecular-level insights into how ChR2 photochemical reactivity can be controlled, which will lay the foundation for improving the design of optogenetic tools. In this work, we combine ab initio nonadiabatic dynamics simulation, excited state free energy calculation, and reaction path search to comprehensively characterize the RPSB photoisomerization mechanism in the E123T mutant of ChR2. Our simulation agrees with previous experiments in predicting a red-shifted absorption spectrum and significant slowdown of photoisomerization in the E123T mutant. Interestingly, our simulations predict similar photoisomerization quantum yields for the mutant and WT despite the differences in excited-state lifetime and absorption maximum. Upon mutation, the neutralization of the negative charge on the E123 residue increases the isomerization barrier, alters the reaction pathway, and changes the relative stability of two fluorescent states. Our findings provide new insight into the intricate role of the electrostatic environment on the RPSB photoisomerization mechanism in microbial rhodopsins.

Putting Photomechanical Switches to Work: An Ab Initio Multiple Spawning Study of Donor-Acceptor Stenhouse Adducts

Photomechanical switches are light sensitive molecules capable of transducing the energy of a photon into mechanical work via photodynamics. In this Letter, we present the first atomistic investigation of the photodynamics of a novel class of photochromes called donor-acceptor Stenhouse adducts (DASA) using state-of-the-art ab initio multiple spawning interfaced with state-averaged complete active-space self-consistent field theory. Understanding the Z/E photoisomerization mechanism in DASAs at the molecular level is crucial in designing new derivatives with improved photoswitching capabilities. Our dynamics simulations show that the actinic step consists of competing nonradiative relaxation pathways that collectively contribute to DASAs' low (21% in toluene) photoisomerization quantum yield. Furthermore, we highlight the important role the intramolecular hydrogen bond plays in the selectivity of photoisomerization in DASAs, identifying it as a possible structural element to tune DASA properties. Our fully ab initio simulations reveal the key degrees of freedom involved in the actinic step, paving the way for the rational design of new generations of DASAs with improved quantum yield and efficiency.

TeraChem: Accelerating electronic structure and ab initio molecular dynamics with graphical processing units

Developed over the past decade, TeraChem is an electronic structure and ab initio molecular dynamics software package designed from the ground up to leverage graphics processing units (GPUs) to perform large-scale ground and excited state quantum chemistry calculations in the gas and the condensed phase. TeraChem's speed stems from the reformulation of conventional electronic structure theories in terms of a set of individually optimized high-performance electronic structure operations (e.g., Coulomb and exchange matrix builds, one- and two-particle density matrix builds) and rank-reduction techniques (e.g., tensor hypercontraction). Recent efforts have encapsulated these core operations and provided language-agnostic interfaces. This greatly increases the accessibility and flexibility of TeraChem as a platform to develop new electronic structure methods on GPUs and provides clear optimization targets for emerging parallel computing architectures.

Computation of Molecular Ionization Energies Using an Ensemble Density Functional Theory Method

Computation of the ionization energies and of the respective Dyson orbitals based on the use of the extended Koopmans theorem (EKT) is implemented in connection with an ensemble density functional theory (eDFT) method, the state-interaction state-averaged spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS or SSR) method. The new methodology enables fast computation of the ionization energies and evaluation of the respective Dyson orbitals, the square norms of which are related with the ionization probabilities, in the ground and excited electronic states of molecules. As the application of EKT recycles the intermediate quantities from the SSR analytical energy gradient, evaluation of the ionization energies and probabilities can be carried out on-the-fly during the nonadiabatic molecular dynamics simulations. This opens up a perspective for fast theoretical simulation of the time-resolved photoelectron spectroscopy observations. In the present work, the new methodology is tested in the computation of the ionization energies and Dyson orbitals of several molecules in the ground and excited electronic states, including strongly correlated species, such as the ozone molecule, dissociating chemical bonds, and conical intersections.

Combining Graphics Processing Units, Simplified Time-Dependent Density Functional Theory, and Finite-Difference Couplings to Accelerate Nonadiabatic Molecular Dynamics

Starting from our recently published implementation of nonadiabatic molecular dynamics (NAMD) on graphics processing units (GPUs), we explore further approaches to accelerate ab initio NAMD calculations at the time-dependent density functional theory (TDDFT) level of theory. We employ (1) the simplified TDDFT schemes of Grimme et al. and (2) the Hammes-Schiffer-Tully approach to obtain nonadiabatic couplings from finite-difference calculations. The resulting scheme delivers an accurate physical picture while virtually eliminating the two computationally most demanding steps of the algorithm. Combined with our GPU-based integral routines for SCF, TDDFT, and TDDFT derivative calculations, NAMD simulations of systems of a few hundreds of atoms at a reasonable time scale become accessible on a single compute node. To demonstrate this and to present a first, illustrative example, we perform TDDFT/MM-NAMD simulations of the rhodopsin protein.

Observe while it happens: catching photoactive proteins in the act with non-adiabatic molecular dynamics simulations

Organisms use photo-receptors to react to light. The first step is usually the absorption of a photon by a prosthetic group embedded inside the photo-receptor, often a conjugated chromophore. The electronic changes in the chromophore induced by photo-absorption can trigger a cascade of structural or chemical transformations that culminate into a response to light. Understanding how these proteins have evolved to mediate their activation process has remained challenging because the required time and spacial resolutions are notoriously difficult to achieve experimentally. Therefore, mechanistic insights into photoreceptor activation have been predominantly obtained with computer simulations. Here we briefly outline the challenges associated with such computations and review the progress made in this field.