Direct QM/MM simulation of photoexcitation dynamics in bacteriorhodopsin and halorhodopsin

Particle-Breaking Hartree-Fock Theory for Open Molecular Systems

In this work we present the particle-breaking Hartree-Fock (PBHF) model which is a mean-field approach to open molecular systems. The interaction of a system with the environment is parametrized through a particle-breaking term in the molecular Hamiltonian. The PBHF wave function is constructed through an exponential unitary transformation of a Slater determinant with a given number of electrons. We consider only the closed-shell formalism. The parametrization results in a linear combination of Slater determinants with different numbers of electrons, i.e., the PBHF wave function is not an eigenfunction of the number operator. As a result, the density matrix may have fractional occupations which are, because of the unitary parametrization, always between 0.0 and 2.0. The occupations are optimized simultaneously with the orbitals, using the trust-region optimization procedure. In the limit of a particle-conserving Hamiltonian, the PBHF optimization will converge to a standard Hartree-Fock wave function. We show that the average number of electrons may be decreased or increased depending on whether the particle-breaking term affects occupied or virtual orbitals.

A Unified View on Varied Ultrafast Dynamics of the Primary Process in Microbial Rhodopsins

All-trans to 13-cis photoisomerization of the protonated retinal Schiff base (PRSB) chromophore is the primary step that triggers various biological functions of microbial rhodopsins. While this ultrafast primary process has been extensively studied, it has been recognized that the relevant excited-state relaxation dynamics differ significantly from one rhodopsin to another. To elucidate the origin of the complicated ultrafast dynamics of the primary process in microbial rhodopsins, we studied the excited-state dynamics of proteorhodopsin, its D97N mutant, and bacteriorhodopsin by femtosecond time-resolved absorption (TA) spectroscopy in a wide pH range. The TA data showed that their excited-state relaxation dynamics drastically change when pH approaches the pK_a of the counterion residue of the PRSB chromophore in the ground state. This result reveals that the varied excited-state relaxation dynamics in different rhodopsins mainly originate from the difference of the ground-state heterogeneity (i.e., protonation/deprotonation of the PRSB counterion).

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.

Early-stage dynamics of chloride ion-pumping rhodopsin revealed by a femtosecond X-ray laser

Light-driven rhodopsin proteins pump ions across cell membranes. They have applications in optogenetics and can potentially be used to develop solar energy–harvesting devices. A detailed understanding of rhodopsin dynamics and functions may therefore assist research in medicine, health, and clean energy. This time-resolved crystallography study carried out with X-ray free-electron lasers reveals detailed dynamics of chloride ion–pumping rhodopsin (ClR) within 100 ps of light activation. It shows the dissociation of Cl⁻ from the Schiff base binding site upon light-triggered retinal isomerization. This Cl⁻ dissociation is followed by diffusion toward the intracellular direction. The results hint at a common ion-pumping mechanism across rhodopsin families.

Chloride ion–pumping rhodopsin (ClR) in some marine bacteria utilizes light energy to actively transport Cl⁻ into cells. How the ClR initiates the transport is elusive. Here, we show the dynamics of ion transport observed with time-resolved serial femtosecond (fs) crystallography using the Linac Coherent Light Source. X-ray pulses captured structural changes in ClR upon flash illumination with a 550 nm fs-pumping laser. High-resolution structures for five time points (dark to 100 ps after flashing) reveal complex and coordinated dynamics comprising retinal isomerization, water molecule rearrangement, and conformational changes of various residues. Combining data from time-resolved spectroscopy experiments and molecular dynamics simulations, this study reveals that the chloride ion close to the Schiff base undergoes a dissociation–diffusion process upon light-triggered retinal isomerization.

Quantum and Quantum-Classical Studies of the Photoisomerization of a Retinal Chromophore Model

We report an in-depth analysis of the photo-induced isomerization of the 2-cis-penta-2,4-dieniminium cation: a minimal model of the 11-cis retinal protonated Schiff base chromophore of the dim-light photoreceptor rhodopsin. Based on recently developed three-dimensional potentials parametrized on ab initio multi-state multi-configurational second-order perturbation theory data, we perform quantum-dynamical studies. In addition, simulations based on various quantum-classical methods, among which Tully surface hopping and the coupled-trajectory approach derived from the exact factorization, allow us to validate their performance against vibronic wavepacket propagation and, therefore, a purely quantum treatment. Quantum-dynamics results uncover qualitative differences with respect to the two-dimensional Hahn-Stock potentials, widely used as model potentials for the isomerization of the same chromophore, due to the increased dimensionality and three-mode correlation. Quantum-classical simulations show, instead, that three-dimensional model potentials are capable of capturing a number of features revealed by atomistic simulations and experimental observations. In particular, a recently reported vibrational phase relationship between double-bond torsion and hydrogen-out-of-plane modes critical for rhodopsin isomerization efficiency is correctly reproduced.

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.

Fluorescence quenching by photoinduced electron transfer between 7-methoxycoumarin and guanine base facilitated by hydrogen bonds: an in silico study

In this study, the effects of hydrogen bond (H-bond) formation on fluorescence quenching of 7-methoxycoumarin (7MC) via photo-induced electron transfer from a guanine base (Gua) are investigated using a combined quantum mechanics/molecular mechanics simulation. The electronic structure is calculated by the floating occupation molecular orbital complete active space configuration interaction modification on a semiempirical method. Then the full multiple spawning method is employed for the dynamics simulations on multiple electronic states. The methods employed here are validated by simulating direct dynamics of 7MC (without Gua) and compared with available experimental results. Our computational results are in good agreement with the previously reported experimental results in terms of spectroscopic properties of 7MC. In the case of a H-bonded 7MC-Gua complex, the results from constrained dynamics simulations and single-point calculations suggest that the electron transfer occurs on the second excited state and it depends not only on the H-bond length but also on the intermolecular planarity between 7MC and Gua. Moreover, a proton coupled electron transfer can occur at ≈1 Å of H-bond length, where a proton from Gua is also transferred together with the electron to 7MC. The obtained simulations are expected to be greatly beneficial for designing effective fluorescently labeled nucleotide probes as well as providing information for precise fluorescence signal interpretation.

Excited States and Photochemistry of Chromophores in the Photoactive Proteins Explored by the Combined Quantum Mechanical and Molecular Mechanical Calculations

A photoactive protein usually contains a unique chromophore that is responsible for the initial photoresponse and functions of the photoactive protein are determined by the interaction between the chromophore and its protein surroundings. The combined quantum mechanical and molecular mechanical (QM/MM) approach is demonstrated to be a very useful tool for exploring structures and functions of a photoactive protein with the chromophore and its protein surroundings treated by the QM and MM methods, respectively. In this review, we summarize the basic formulas of the QM/MM approach and emphasize its applications to excited states and photoreactions of chromophores in rhodopsin protein, photoactive yellow protein, and green fluorescent protein.