Electrostatic Control of Photoisomerization in Channelrhodopsin 2

Theoretical Insight into the Fluorescence Spectral Tuning Mechanism: A Case Study of Flavin-Dependent Bacterial Luciferase

Bioluminescence of bacteria is widely applied in biological imaging, environmental toxicant detection, and many other situations. Understanding the spectral tuning mechanism not only helps explain the diversity of colors observed in nature but also provides principles for bioengineering new color variants for practical applications. In this study, time-dependent density functional theory (TD-DFT) and quantum mechanics and molecular mechanics (QM/MM) calculations have been employed to understand the fluorescence spectral tuning mechanism of bacterial luciferase with a focus on the electrostatic effect. The spectrum can be tuned by both a homogeneous dielectric environment and oriented external electric fields (OEEFs). Increasing the solvent polarity leads to a redshift of the fluorescence emission maximum, λF, accompanied by a substantial increase in density. In contrast, applying an OEEF along the long axis of the isoalloxazine ring (X-axis) leads to a significant red- or blue-shift in λF, depending on the direction of the OEEF, yet with much smaller changes in intensity. The effect of polar solvents is directionless, and the red-shifts can be attributed to the larger dipole moment of the S1 state compared with that of the S0 state. However, the effect of OEEFs directly correlates with the difference dipole moment between the S1 and S0 states, which is directional and is determined by the charge redistribution upon deexcitation. Moreover, the electrostatic effect of bacterial luciferase is in line with the presence of an internal electric field (IEF) pointing in the negative X direction. Finally, the key residues that contribute to this IEF and strategies for modulating the spectrum through site-directed point mutations are discussed.

Calcium-Sensitive Microbial Rhodopsin VirChR1: A Femtosecond to Second Photocycle Study

Viral rhodopsins are light-gated cation channels representing a novel class of microbial rhodopsins. For viral rhodopsin 1 subfamily members VirChR1 and OLPVR1, channel activity is abolished above a certain calcium concentration. Here we present a calcium-dependent spectroscopic analysis of VirChR1 on the femtosecond to second time scale. Unlike channelrhodopsin-2, VirChR1 possesses two intermediate states P1 and P2 on the ultrafast time scale, similar to J and K in ion-pumping rhodopsins. Subsequently, we observe multifaceted photocycle kinetics with up to seven intermediate states. Calcium predominantly affects the last photocycle steps, including the appearance of additional intermediates P6Ca and P7 representing the blocked channel. Furthermore, the photocycle of the counterion variant D80N is drastically altered, yielding intermediates with different spectra and kinetics compared to those of the wt. These findings demonstrate the central role of the counterion within the defined reaction sequence of microbial rhodopsins that ultimately defines the protein function.

Postoperatively Noninvasive Optogenetic Stimulation via Upconversion Nanoparticles Enhancing Sciatic Nerve Repair

The efficacy of electrical stimulation facilitating peripheral nerve regeneration is evidenced extensively, while the associated secondary damage resulting from repeated electrode invasion and indiscriminate stimulation is inevitable. Here, we present an optogenetics strategy that utilizes upconversion nanoparticles (UCNPs) to convert deeply penetrating near-infrared excitation into blue emission, which activates an adeno-associated virus-encoding ChR2 photoresponsive ion channel on cell membranes. The induced Ca2+ flux, similar to the ion flux in the electrical stimulation approach, efficiently regulates viability and proliferation, secretion of nerve growth factor, and neural function of RSC96 cells. Furthermore, deep near-infrared excitation is harnessed to stimulate autologous Schwann cells in situ via a UCNP-composited scaffold, which enhances nerve sprouting and myelination, consequently promoting functional recovery, electrophysiological restoration, and reinnervation of damaged nerves. This developed postoperatively noninvasive optogenetics strategy presents a novel, minimally traumatic, and enduring therapeutic stimulus to effectively promote peripheral nerve repair.

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.

Retinal photoisomerization versus counterion protonation in light and dark-adapted bacteriorhodopsin and its primary photoproduct

Discovered over 50 years ago, bacteriorhodopsin is the first recognized and most widely studied microbial retinal protein. Serving as a light-activated proton pump, it represents the archetypal ion-pumping system. Here we compare the photochemical dynamics of bacteriorhodopsin light and dark-adapted forms with that of the first metastable photocycle intermediate known as "K". We observe that following thermal double isomerization of retinal in the dark from bio-active all-trans 15-anti to 13-cis, 15-syn, photochemistry proceeds even faster than the ~0.5 ps decay of the former, exhibiting ballistic wave packet curve crossing to the ground state. In contrast, photoexcitation of K containing a 13-cis, 15-anti chromophore leads to markedly multi-exponential excited state decay including much slower stages. QM/MM calculations, aimed to interpret these results, highlight the crucial role of protonation, showing that the classic quadrupole counterion model poorly reproduces spectral data and dynamics. Single protonation of ASP212 rectifies discrepancies and predicts triple ground state structural heterogeneity aligning with experimental observations. These findings prompt a reevaluation of counter ion protonation in bacteriorhodopsin and contribute to the broader understanding of its photochemical dynamics.

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.

Optogenetics: Illuminating the Future of Hearing Restoration and Understanding Auditory Perception

Hearing loss is a prevalent sensory impairment significantly affecting communication and quality of life. Traditional approaches for hearing restoration, such as cochlear implants, have limitations in frequency resolution and spatial selectivity. Optogenetics, an emerging field utilizing light-sensitive proteins, offers a promising avenue for addressing these limitations and revolutionizing hearing rehabilitation. This review explores the methods of introducing Channelrhodopsin- 2 (ChR2), a key light-sensitive protein, into cochlear cells to enable optogenetic stimulation. Viral- mediated gene delivery is a widely employed technique in optogenetics. Selecting a suitable viral vector, such as adeno-associated viruses (AAV), is crucial in efficient gene delivery to cochlear cells. The ChR2 gene is inserted into the viral vector through molecular cloning techniques, and the resulting viral vector is introduced into cochlear cells via direct injection or round window membrane delivery. This allows for the expression of ChR2 and subsequent light sensitivity in targeted cells. Alternatively, direct cell transfection offers a non-viral approach for ChR2 delivery. The ChR2 gene is cloned into a plasmid vector, which is then combined with transfection agents like liposomes or nanoparticles. This mixture is applied to cochlear cells, facilitating the entry of the plasmid DNA into the target cells and enabling ChR2 expression. Optogenetic stimulation using ChR2 allows for precise and selective activation of specific neurons in response to light, potentially overcoming the limitations of current auditory prostheses. Moreover, optogenetics has broader implications in understanding the neural circuits involved in auditory processing and behavior. The combination of optogenetics and gene delivery techniques provides a promising avenue for improving hearing restoration strategies, offering the potential for enhanced frequency resolution, spatial selectivity, and improved auditory perception.

Computational Characterization of the Reactivity of Compound I in Unspecific Peroxygenases

Unspecific peroxygenases (UPOs) are emerging as promising biocatalysts for selective oxyfunctionalization of unactivated C-H bonds. However, their potential in large-scale synthesis is currently constrained by suboptimal chemical selectivity. Improving the selectivity of UPOs requires a deep understanding of the molecular basis of their catalysis. Recent molecular simulations have sought to unravel UPO's selectivity and inform their design principles. However, most of these studies focused on substrate-binding poses. Few researchers have investigated how the reactivity of CpdI, the principal oxidizing intermediate in the catalytic cycle, influences selectivity in a realistic protein environment. Moreover, the influence of protein electrostatics on the reaction kinetics of CpdI has also been largely overlooked. To bridge this gap, we used multiscale simulations to interpret the regio- and enantioselective hydroxylation of the n-heptane substrate catalyzed by Agrocybe aegerita UPO (AaeUPO). We comprehensively characterized the energetics and kinetics of the hydrogen atom-transfer (HAT) step, initiated by CpdI, and the subsequent oxygen rebound step forming the product. Notably, our approach involved both free energy and potential energy evaluations in a quantum mechanics/molecular mechanics (QM/MM) setting, mitigating the dependence of results on the choice of initial conditions. These calculations illuminate the thermodynamics and kinetics of the HAT and oxygen rebound steps. Our findings highlight that both the conformational selection and the distinct chemical reactivity of different substrate hydrogen atoms together dictate the regio- and enantio-selectivity. Building on our previous study of CpdI's formation in AaeUPO, our results indicate that the HAT step is the rate-limiting step in the overall catalytic cycle. The subsequent oxygen rebound step is swift and retains the selectivity determined by the HAT step. We also pinpointed several polar and charged amino acid residues whose electrostatic potentials considerably influence the reaction barrier of the HAT step. Notably, the Glu196 residue is pivotal for both the CpdI's formation and participation in the HAT step. Our research offers in-depth insights into the catalytic cycle of AaeUPO, which will be instrumental in the rational design of UPOs with enhanced properties.

Impact of protein-chromophore interaction on the retinal excited state and photocycle of Gloeobacter rhodopsin: role of conserved tryptophan residues

The function of microbial as well as mammalian retinal proteins (aka rhodopsins) is associated with a photocycle initiated by light excitation of the retinal chromophore of the protein, covalently bound through a protonated Schiff base linkage. Although electrostatics controls chemical reactions of many organic molecules, attempt to understand its role in controlling excited state reactivity of rhodopsins and, thereby, their photocycle is scarce. Here, we investigate the effect of highly conserved tryptophan residues, between which the all-trans retinal chromophore of the protein is sandwiched in microbial rhodopsins, on the charge distribution along the retinal excited state, quantum yield and nature of the light-induced photocycle and absorption properties of Gloeobacter rhodopsin (GR). Replacement of these tryptophan residues by non-aromatic leucine (W222L and W122L) or phenylalanine (W222F) does not significantly affect the absorption maximum of the protein, while all the mutants showed higher sensitivity to photobleaching, compared to wild-type GR. Flash photolysis studies revealed lower quantum yield of trans-cis photoisomerization in W222L as well as W222F mutants relative to wild-type. The photocycle kinetics are also controlled by these tryptophan residues, resulting in altered accumulation and lifetime of the intermediates in the W222L and W222F mutants. We propose that protein-retinal interactions facilitated by conserved tryptophan residues are crucial for achieving high quantum yield of the light-induced retinal isomerization, and affect the thermal retinal re-isomerization to the resting state.

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

This review discusses how ultrafast organic photochemical reactions are controlled by conical intersections, highlighting that decay to the ground-state at multiple points of the intersection space results in their multi-mode character.

Clustered tree regression to learn protein energy change with mutated amino acid

Accurate and effective prediction of mutation-induced protein energy change remains a great challenge and of great interest in computational biology. However, high resource consumption and insufficient structural information of proteins severely limit the experimental techniques and structure-based prediction methods. Here, we design a structure-independent protocol to accurately and effectively predict the mutation-induced protein folding free energy change with only sequence, physicochemical and evolutionary features. The proposed clustered tree regression protocol is capable of effectively exploiting the inherent data patterns by integrating unsupervised feature clustering by K-means and supervised tree regression using XGBoost, and thus enabling fast and accurate protein predictions with different mutations, with an average Pearson correlation coefficient of 0.83 and an average root-mean-square error of 0.94kcal/mol. The proposed sequence-based method not only eliminates the dependence on protein structures, but also has potential applications in protein predictions with rare structural information.

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.

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.

Time-resolved photoacoustics of channelrhodopsins: early energetics and light-driven volume changes

In biological photoreceptors, the energy stored in early transient species is a key feature to drive the photocycle or a chain of reactions. Time-resolved photoacoustics (PA) can explore the energy landscape of transient species formed within few ns after photoexcitation, as well as volumetric changes (ΔV) of these intermediates with respect to the parental state. In this work, PA identified these important parameters for several channelrhodopsins, namely CaChR1 from Chlamydomonas augustae and CrChR2 from Chlamydomonas reinhardtii and various variants. PA has access to the sub-ns formation of the early photoproduct P1 and to its relaxation, provided that this latter process occurs within a few μs. We found that ΔVP1 for CaChR1 is ca. 12 mL/mol, while it is much smaller for CrChR2 (4.7 mL/mol) and for H. salinarum bacteriorhodopsin (HsBR, ΔVK = 2.8 mL/mol). PA experiments on variants strongly indicate that part of this large ΔVP1 value for CaChR1 is caused by the protonation dynamics of the Schiff base counterion complex involving E169 and D299. PA data further show that the energy level of P1 is higher in CrChR2 (ca. 96 kJ/mol) than in CaChr1 (ca. 46 kJ/mol), comparable to the energy level of the K state of HsBR (60 kJ/mol). Instrumental to gain these molecular values from the raw PA data was the estimation of the quantum yield (Φ) for P1 formation via transient spectroscopy; for both channelrhodopsins, ΦP2 was evaluated as ca. 0.4.

Graphical Abstract

Variational Density Functional Calculations of Excited States: Conical Intersection and Avoided Crossing in Ethylene Bond Twisting

Theoretical studies of photochemical processes require a description of the energy surfaces of excited electronic states, especially near degeneracies, where transitions between states are most likely. Systems relevant to photochemical applications are typically too large for high-level multireference methods, and while time-dependent density functional theory (TDDFT) is efficient, it can fail to provide the required accuracy. A variational, time-independent density functional approach is applied to the twisting of the double bond and pyramidal distortion in ethylene, the quintessential model for photochemical studies. By allowing for symmetry breaking, the calculated energy surfaces exhibit the correct topology around the twisted-pyramidalized conical intersection even when using a semilocal functional approximation, and by including explicit self-interaction correction, the torsional energy curves are in close agreement with published multireference results. The findings of the present work point to the possibility of using a single determinant time-independent density functional approach to simulate nonadiabatic dynamics, even for large systems where multireference methods are impractical and TDDFT is often not accurate enough.

Effects of Enzyme-Ligand Interactions on the Photoisomerization of a Light-Regulated Chemotherapeutic Drug

Molecular photoswitches permit using light to control protein activity with high spatiotemporal resolutions, thereby alleviating the side effects of conventional chemotherapy. However, due to the challenges in probing ultrafast photoisomerization reactions in biological environments, it remains elusive how the protein influences the photochemistry of the photoswitches, which hampers the rational design of light-regulated therapeutics. To overcome this challenge, we employed first-principles nonadiabatic dynamics simulations to characterize the photodynamics of the phototrexate (PTX), a recently developed photoswitchable anticancer chemotherapeutic that reversibly inhibits its target enzyme dihydrofolate reductase (DHFR). Our simulations show that the protein environment impedes the trans to cis photoisomerization of the PTX. The confinement in the ligand-binding cavity slows down the isomerization kinetics and quantum yield of the photoswitch by reshaping its conical intersection, increasing its excited-state free-energy barrier and quenching its local density fluctuations. Also, the protein environment results in a suboptimal binding mode of the photoproduct that needs to undergo large structural rearrangement to effectively inhibit the enzyme. Therefore, we predict that the PTX's trans → cis photoisomerization in solution precedes its binding with the protein, despite the favorable binding energy of the trans isomer. Our findings highlight the importance of the protein environment on the photochemical reactions of the molecular photoswitches. As such, our work represents an important step toward the rational design of light-regulated drugs in photopharmacology.

Electrostatic profiling of photosynthetic pigments: implications for directed spectral tuning

Photosynthetic pigment-protein complexes harvest solar energy with a high quantum efficiency. Protein scaffolds are known to tune the spectral properties of embedded pigments principally through structured electrostatic environments. Although the physical nature of electrostatic tuning is straightforward, the precise spatial principles of electrostatic preorganization remain poorly explored for different protein matrices and incompletely characterized with respect to the intrinsic properties of different photosynthetic pigments. In this work, we study the electronic structure features associated with the lowest excited state of a series of eight naturally occurring (bacterio)chlorophylls and pheophytins to describe the precise topological differences in electrostatic potentials and hence determine intrinsic differences in the expected mode and impact of electrostatic tuning. The difference electrostatic potentials between the ground and first excited states are used as fingerprints. Both the spatial profile and the propensity for spectral tuning are found to be unique for each pigment, indicating spatially and directionally distinct modes of electrostatic tuning. The results define a specific partitioning of the protein matrix around each pigment as an aid to identify regions with a maximal impact on spectral tuning and have direct implications for dimensionality reduction in protein design and engineering. Thus, a quantum mechanical basis is provided for understanding, predicting, and ultimately designing sequence-modified or pigment-exchanged biological systems, as suggested for selected examples of pigment-reconstituted proteins.