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.

Nonadiabatic Dynamics Simulation of the Wavelength-Dependent Photochemistry of Azobenzene Excited to the nπ* and ππ* Excited States.

Azobenzene is one of the most ubiquitous photoswitches in photochemistry and a prototypical model for photoisomerizing systems. Despite this, its wavelength-dependent photochemistry has puzzled researchers for decades. Upon excitation to the higher energy ππ* excited state instead of the dipole-forbidden nπ* state, the quantum yield of isomerization from trans- to cis-azobenzene is halved. The difficulties associated with unambiguously resolving this effect both experimentally and theoretically have contributed to lasting controversies regarding the photochemistry of azobenzene. Here, we systematically characterize the dynamic photoreaction pathways of azobenzene by performing first-principles simulations of the nonadiabatic dynamics following excitation to both the ππ* and the nπ* states. We demonstrate that ground-state recovery is mediated by two distinct S1 decay pathways: a reactive twisting pathway and an unreactive planar pathway. Increased preference for the unreactive pathway upon ππ* excitation largely accounts for the wavelength-dependent behavior observed in azobenzene.

Hole-hole Tamm-Dancoff-approximated density functional theory: A highly efficient electronic structure method incorporating dynamic and static correlation.

The study of photochemical reaction dynamics requires accurate as well as computationally efficient electronic structure methods for the ground and excited states. While time-dependent density functional theory (TDDFT) is not able to capture static correlation, complete active space self-consistent field methods neglect much of the dynamic correlation. Hence, inexpensive methods that encompass both static and dynamic electron correlation effects are of high interest. Here, we revisit hole-hole Tamm-Dancoff approximated (hh-TDA) density functional theory for this purpose. The hh-TDA method is the hole-hole counterpart to the more established particle-particle TDA (pp-TDA) method, both of which are derived from the particle-particle random phase approximation (pp-RPA). In hh-TDA, the N-electron electronic states are obtained through double annihilations starting from a doubly anionic (N+2 electron) reference state. In this way, hh-TDA treats ground and excited states on equal footing, thus allowing for conical intersections to be correctly described. The treatment of dynamic correlation is introduced through the use of commonly employed density functional approximations to the exchange-correlation potential. We show that hh-TDA is a promising candidate to efficiently treat the photochemistry of organic and biochemical systems that involve several low-lying excited states-particularly those with both low-lying ππ* and nπ* states where inclusion of dynamic correlation is essential to describe the relative energetics. In contrast to the existing literature on pp-TDA and pp-RPA, we employ a functional-dependent choice for the response kernel in pp- and hh-TDA, which closely resembles the response kernels occurring in linear response and collinear spin-flip TDDFT.

First-Principles Characterization of the Elusive I Fluorescent State and the Structural Evolution of Retinal Protonated Schiff Base in Bacteriorhodopsin.

The conversion of light energy into work is essential to life on earth. Bacteriorhodopsin (bR), a light-activated proton pump in Archae, has served for many years as a model system for the study of this process in photoactive proteins. Upon absorption of a photon, its chromophore, the retinal protonated Schiff base (RPSB), isomerizes from its native all-trans form to a 13-cis form and pumps a proton out of the cell in a process that is coupled to eventual ATP synthesis. Despite numerous time-resolved spectroscopic studies over the years, the details of the photodynamics of bR on the excited state, particularly the characterization of the I fluorescent state, the time-resolved reaction mechanism, and the role of the counterion cluster of RPSB, remain uncertain. Here, we use ab initio multiple spawning (AIMS) with spin-restricted ensemble Kohn-Sham (REKS) theory to simulate the nonadiabatic dynamics of the ultrafast photoreaction in bR. The excited state dynamics can be partitioned into three distinct phases: (1) relaxation away from the Franck-Condon region dominated by changes in retinal bond length alternation, (2) dwell time on the excited state in the I fluorescent state featuring an untwisted, bond length inverted RPSB, and (3) rapid torsional evolution to the conical intersection after overcoming a small excited state barrier. We fully characterize the I fluorescent state and the excited state barrier that hinders direct evolution to the conical intersection following photoexcitation. We also find that photoisomerization is accompanied by weakening of the interaction between RPSB and its counterion cluster. However, in contradiction with a recent time-resolved X-ray experiment, hydrogen bond cleavage is not necessary to reproduce the observed photoisomerization dynamics.

Malignancies in Immunoglobulin G4-related ophthalmic disease.

PURPOSE: Clinical phenotypes in Immunoglobulin G4-related disease (IgG4-RD) according to the patterns of affecting organs have different risks of malignancies. We attempt to determine the association of malignancies with IgG4-related ophthalmic disease (IgG4-ROD). DESIGN: Retrospective cohort study. METHODS: Review of medical records, orbital images and histopathology reports in a territory-wide cohort of biopsy proven IgG4-ROD patients from 2005-2019. FINDINGS: Among 122 patients who had biopsies taken from adnexal lesions including lacrimal glands (n = 108), orbital mass (n = 30), infiltrated orbital fat (n = 10), conjunctiva (n = 2) or extraocular muscles (n = 3), 13% (16/122) developed malignancies over 73 ± 48months' follow-up. There were 9 cases of ocular adnexal lymphoma (OAL) and 7 extra-orbital malignancies. Compared with the general population, the incidence of OAL was significantly higher (standardized incidence ratios, SIRs = 10.0, 95%CI = 4.5-17.6) while that of extra-orbital malignancies was similar. The SIRs was highest within the first year (SIR = 46.7, 95%CI = 18.5-87.6) when 7 OAL were concomitantly diagnosed. Patients who developed OAL or extra-orbital malignancies were older than other patients at IgG4-ROD diagnosis (64.9 ± 7.1, 68.3 ± 8.5 versus 55.2 ± 15.0 years, P < 0.05). Asymmetric lacrimal gland enlargement (78% versus 13%), lack of frontal (0% versus 12%) or infraorbital nerve enlargement (0% versus 36%) were associated with OAL (all P < 0.05). Pre-treatment serum IgG4 level or extra-orbital IgG4-RD involvement was similar among patients with or without malignancies. CONCLUSION: In this biopsy-proven IgG4-ROD cohort, 7% developed OAL which was 10 times higher than the general population. Patients with asymmetric lacrimal gland enlargement or without trigeminal nerves involvement radiologically were associated with OAL.

Predictions of Pre-edge Features in Time-Resolved Near-Edge X-ray Absorption Fine Structure Spectroscopy from Hole-Hole Tamm-Dancoff-Approximated Density Functional Theory.

Time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS) spectroscopy is a powerful technique for studying photochemical reaction dynamics with femtosecond time resolution. In order to avoid ambiguity in TR-NEXAFS spectra from nonadiabatic dynamics simulations, core- and valence-excited states must be evaluated on equal footing and those valence states must also define the potential energy surfaces used in the nonadiabatic dynamics simulation. In this work, we demonstrate that hole-hole Tamm-Dancoff-approximated density functional theory (hh-TDA) is capable of directly simulating TR-NEXAFS spectroscopies. We apply hh-TDA to the excited-state dynamics of acrolein. We identify two pre-edge features in the oxygen K-edge TR-NEXAFS spectrum associated with the S2 (ππ*) and S1 (nπ*) excited states. We show that these features can be used to follow the internal conversion dynamics between the lowest three electronic states of acrolein. Due to the low, O(N2) apparent computational complexity of hh-TDA and our GPU-accelerated implementation, this method is promising for the simulation of pre-edge features in TR-NEXAFS spectra of large molecules and molecules in the condensed phase.

Ab Initio Nonadiabatic Molecular Dynamics with Hole-Hole Tamm-Dancoff Approximated Density Functional Theory.

The study of photoinduced dynamics in chemical systems necessitates accurate and computationally efficient electronic structure methods, especially as the systems of interest grow larger. The linear response hole-hole Tamm-Dancoff approximated (hh-TDA) density functional theory method was recently proposed to satisfy such demands. The N-electron electronic states are obtained by means of double annihilations on a doubly anionic (N + 2)-electron reference state, allowing for the ground and excited states to be formed on the same footing and thus enabling the correct description of conical intersections. Dynamic electron correlation effects are incorporated by means of the exchange-correlation functional. The accuracy afforded by the simultaneous treatment of static and dynamic correlation in addition to the relatively low computational cost, comparable to that of time-dependent density functional theory (TDDFT), makes it a promising ab initio electronic structure method for on-the-fly generation of potential energy surfaces in nonadiabatic dynamics simulations of photochemical systems, particularly those for which the nπ* and ππ* electronic excitations are most relevant. Here, we apply the hh-TDA method to nonadiabatic dynamics simulations of prototypical photochemical processes. First, we demonstrate the ability of hh-TDA to adequately describe conical intersection geometries. We next examine its ability to describe the ultrafast excited state dynamics of photoexcited ethylene through an ab initio multiple spawning (AIMS) dynamics simulation. Finally, we present an alternative variant of the hh-TDA method, which uses orbitals from a fractional occupation number Kohn-Sham (FON-KS) calculation applied to an ensemble with N-electrons. The resulting method is termed floating occupation molecular orbital hh-TDA (FOMO-hh-TDA). This scheme allows us to combine hh-TDA with global hybrid functionals and allows us to avoid unbound valence orbitals that may result from an (N + 2)-electron self-consistent field (SCF) procedure. FOMO-hh-TDA-BHLYP faithfully reproduces the nonadiabatic dynamics of trans-azobenzene (TAB) and is used to characterize the excited state decay pathways from the first (nπ*) excited state.

Simultaneous observation of nuclear and electronic dynamics by ultrafast electron diffraction.

Simultaneous observation of nuclear and electronic motion is crucial for a complete understanding of molecular dynamics in excited electronic states. It is challenging for a single experiment to independently follow both electronic and nuclear dynamics at the same time. Here we show that ultrafast electron diffraction can be used to simultaneously record both electronic and nuclear dynamics in isolated pyridine molecules, naturally disentangling the two components. Electronic state changes (S1→S0 internal conversion) were reflected by a strong transient signal in small-angle inelastic scattering, and nuclear structural changes (ring puckering) were monitored by large-angle elastic diffraction. Supported by ab initio nonadiabatic molecular dynamics and diffraction simulations, our experiment provides a clear view of the interplay between electronic and nuclear dynamics of the photoexcited pyridine molecule.

Circadian rhythms of Per2::Luc in individual primary mouse hepatocytes and cultures.

BACKGROUND: Hepatocytes, the parenchymal cells of the liver, express core clock genes, such as Period2 and Cryptochrome2, which are involved in the transcriptional/translational feedback loop of the circadian clock. Whether or not the liver is capable of sustaining rhythms independent of a central pacemaker is controversial. Whether and how circadian information may be shared among cells in the liver in order to sustain oscillations is currently unknown. RESULTS: In this study we isolated primary hepatocytes from transgenic Per2(Luc) mice and used bioluminescence as a read-out of the state of the circadian clock. Hepatocytes cultured in a collagen gel sandwich configuration exhibited persistent circadian rhythms for several weeks. The amplitude of the rhythms damped, but medium changes consistently reset the phase and amplitude of the cultures. Cry2(-/-) Per2(Luc) cells oscillated robustly and expressed a longer period. Co-culturing with wildtype cells did not significantly shorten the period, indicating that coupling among hepatocytes is insufficient to synchronize cells with significantly differing periods. However, spatial patterns revealed by cellular imaging of wildtype cultures provided evidence of weak local coupling among the hepatocytes. CONCLUSIONS: Our results with primary hepatocyte cultures demonstrate that cultured hepatocytes are weakly coupled. While this coupling is not sufficient to sustain global synchrony, it does increase local synchrony, which may stabilize the circadian rhythms of peripheral oscillators, such as the liver, against noise in the entraining signals.