Photoisomerization Of Rhodopsin Research Articles

Investigating the mechanism of photoisomerization in jellyfish rhodopsin with the counterion at an atypical position

The position of the counterion in animal rhodopsins plays a crucial role in maintaining visible light sensitivity and facilitating the photoisomerization of their retinal chromophore. The counterion displacement is thought to be closely related to the evolution of rhodopsins, with different positions found in invertebrates and vertebrates. Interestingly, box jellyfish rhodopsin (JelRh) acquired the counterion in transmembrane 2 (TM2) independently. This is a unique feature, as in most animal rhodopsins, the counterion is found in a different location. In this study, we used Fourier Transform Infrared spectroscopy to examine the structural changes that occur in the early photointermediate state of JelRh. We aimed to determine whether the photochemistry of JelRh is similar to that of other animal rhodopsins by comparing its spectra to those of vertebrate bovine rhodopsin (BovRh) and invertebrate squid rhodopsin (SquRh). We observed that the N-D stretching band of the retinal Schiff base was similar to that of BovRh, indicating the interaction between the Schiff base and the counterion is similar in both rhodopsins, despite their different counterion positions. Furthermore, we found that the chemical structure of the retinal in JelRh is similar to that in BovRh, including the changes in the hydrogen-out-of-plane band that indicates a retinal distortion. Overall, the protein conformational changes induced by the photoisomerization of JelRh yielded spectra that resemble an intermediate between BovRh and SquRh, suggesting a unique spectral property of JelRh, and making it the only animal rhodopsin with a counterion in TM2 and an ability to activate Gs protein.

Effects of cell size and bicarbonate on single photon response variability in retinal rods.

Accurate photon counting requires that rods generate highly amplified, reproducible single photon responses (SPRs). The SPR is generated within the rod outer segment (ROS), a multilayered structure built from membranous disks that house rhodopsin. Photoisomerization of rhodopsin at the disk rim causes a local depletion of cGMP that closes ion channels in the plasmalemma located nearby with relative rapidity. In contrast, a photoisomerization at the disk center, distant from the plasmalemma, has a delayed impact on the ion channels due to the time required for cGMP redistribution. Radial differences should be greatest in large diameter rods. By affecting membrane guanylate cyclase activity, bicarbonate could impact spatial inhomogeneity in cGMP content. It was previously known that in the absence of bicarbonate, SPRs are larger and faster at the base of a toad ROS (where the ROS attaches to the rest of the cell) than at the distal tip. Given that bicarbonate enters the ROS at the base and diffuses to the tip and that it expedites flash response recovery, there should be an axial concentration gradient for bicarbonate that would accentuate the base-to-tip SPR differences. Seeking to understand how ROS geometry and bicarbonate affect SPR variability, we used mathematical modeling and made electrophysiological recordings of single rods. Modeling predicted and our experiments confirmed minor radial SPR variability in large diameter, salamander rods that was essentially unchanged by bicarbonate. SPRs elicited at the base and tip of salamander rods were similar in the absence of bicarbonate, but when treated with 30 mM bicarbonate, SPRs at the base became slightly faster than those at the tip, verifying the existence of an axial gradient for bicarbonate. The differences were small and unlikely to undermine visual signaling. However, in toad rods with longer ROSs, bicarbonate somehow suppressed the substantial, axial SPR variability that is naturally present in the absence of bicarbonate. Modeling suggested that the axial gradient of bicarbonate might dampen the primary phototransduction cascade at the base of the ROS. This novel effect of bicarbonate solves a mystery as to how toad vision is able to function effectively in extremely dim light.

Problems of Quantum-Classical Modeling Of the Primary Photoreaction in Rhodopsin

A modified model of the primary photoreaction in rhodopsin, cis-trans photoisomerization of the chromophore (retinal), is studied. The quantum subsystem of the model includes three vibronic states: the ground state, the excited state, and the ground state of the primary photoproduct. These states correspond to three point masses in the classical subsystem. The modification consists in the exponential dependence of the electronic-vibrational coupling constant on the displacement of point masses. The properties of the optimal loci of the multiparameter space, which characterized by the best agreement with the experimental data, are studied. A rather small “multidimensional volume” of these loci shown in all ranges of the used values of the model parameters. Several ways to optimize the quantum-classical model of rhodopsin photoisomerization have been proposed.

Applications of Molecular Dynamics Simulation in Protein Study.

Molecular Dynamics (MD) Simulations is increasingly used as a powerful tool to study protein structure-related questions. Starting from the early simulation study on the photoisomerization in rhodopsin in 1976, MD Simulations has been used to study protein function, protein stability, protein–protein interaction, enzymatic reactions and drug–protein interactions, and membrane proteins. In this review, we provide a brief review for the history of MD Simulations application and the current status of MD Simulations applications in protein studies.

A state-averaged orbital-optimized hybrid quantum–classical algorithm for a democratic description of ground and excited states

In the noisy intermediate-scale quantum (NISQ) era, solving the electronic structure problem from chemistry is considered as the ‘killer application’ for near-term quantum devices. In spite of the success of variational hybrid quantum/classical algorithms in providing accurate energy profiles for small molecules, careful considerations are still required for the description of complicated features of potential energy surfaces. Because the current quantum resources are very limited, it is common to focus on a restricted part of the Hilbert space (determined by the set of active orbitals). While physically motivated, this approximation can severely impact the description of these complicated features. A perfect example is that of conical intersections (i.e. a singular point of degeneracy between electronic states), which are of primary importance to understand many prominent reactions. Designing active spaces so that the improved accuracy from a quantum computer is not rendered useless is key to finding useful applications of these promising devices within the field of chemistry. To answer this issue, we introduce a NISQ-friendly method called ‘state-averaged orbital-optimized variational quantum eigensolver’ which combines two algorithms: (1) a state-averaged orbital-optimizer, and (2) a state-averaged VQE. To demonstrate the success of the method, we classically simulate it on a minimal Schiff base model (namely the formaldimine molecule CH2NH) relevant also for the photoisomerization in rhodopsin—a crucial step in the process of vision mediated by the presence of a conical intersection. We show that merging both algorithms fulfil the necessary condition to describe the molecule’s conical intersection, i.e. the ability to treat degenerate (or quasi-degenerate) states on the same footing.

Position of rhodopsin photoisomerization on the disk surface confers variability to the rising phase of the single photon response in vertebrate rod photoreceptors.

Retinal rods function as accurate photon counters to provide for vision under very dim light. To do so, rods must generate highly amplified, reproducible responses to single photons, yet outer segment architecture and randomness in the location of rhodopsin photoisomerization on the surface of an internal disk introduce variability to the rising phase of the photon response. Soon after a photoisomerization at a disk rim, depletion of cGMP near the plasma membrane closes ion channels and hyperpolarizes the rod. But with a photoisomerization in the center of a disk, local depletion of cGMP is distant from the channels in the plasma membrane. Thus, channel closure is delayed by the time required for the reduction of cGMP concentration to reach the plasma membrane. Moreover, the local fall in cGMP dissipates over a larger volume before affecting the channels, so response amplitude is reduced. This source of variability increases with disk radius. Using a fully space-resolved biophysical model of rod phototransduction, we quantified the variability attributable to randomness in the location of photoisomerization as a function of disk structure. In mouse rods that have small disks bearing a single incisure, this variability was negligible in the absence of the incisure. Variability was increased slightly by the incisure, but randomness in the shutoff of rhodopsin emerged as the main source of single photon response variability at all but the earliest times. Variability arising from randomness in the transverse location of photoisomerization increased in magnitude and persisted over a longer period in the photon response of large salamander rods. A symmetric arrangement of multiple incisures in the disks of salamander rods greatly reduced this variability during the rising phase, but the incisures had the opposite effect on variability arising from randomness in rhodopsin shutoff at later times.

Simulating conical intersection dynamics in the condensed phase with hybrid quantum master equations

We present a framework for simulating relaxation dynamics through a conical intersection of an open quantum system that combines methods to approximate the motion of degrees of freedom with disparate time and energy scales. In the vicinity of a conical intersection, a few degrees of freedom render the nuclear dynamics nonadiabatic with respect to the electronic degrees of freedom. We treat these strongly coupled modes by evolving their wavepacket dynamics in the absence of additional coupling exactly. The remaining weakly coupled nuclear degrees of freedom are partitioned into modes that are fast relative to the nonadiabatic coupling and those that are slow. The fast degrees of freedom can be traced out and treated with second-order perturbation theory in the form of the time-convolutionless master equation. The slow degrees of freedom are assumed to be frozen over the ultrafast relaxation and treated as sources of static disorder. In this way, we adopt the recently developed frozen-mode extension to second-order quantum master equations. We benchmark this approach to numerically exact results in models of pyrazine internal conversion and rhodopsin photoisomerization. We use this framework to study the dependence of the quantum yield on the reorganization energy and the characteristic time scale of the bath in a two-mode model of photoisomerization. We find that the yield is monotonically increasing with reorganization energy for a Markovian bath but monotonically decreasing with reorganization energy for a non-Markovian bath. This reflects the subtle interplay between dissipation and decoherence in conical intersection dynamics in the condensed phase.

Quantum dynamics of multi-dimensional rhodopsin photoisomerization models: Approximate versus accurate treatment of the secondary modes

We numerically investigate the quantum dynamics of the high-dimensional rhodopsin photoisomerization models of Hahn and Stock using the multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) method. The original two-state two-mode model is augmented with additional Raman-active modes observed experimentally and by a bath of low frequency modes that accounts in a generic fashion for the effect the protein and solvent environment. A previously used approximate method in which the two primary modes are treated using an accurate quantum dynamical method while the secondary modes are treated using the time-dependent Hartree (TDH) method is tested against ML-MCTDH calculations. It is shown that the former method does not capture the main effect of these modes on the electronic populations, on the coherent torsional dynamics and on the quantum yield of the photoisomerization reaction. Our ML-MCTDH calculations predict small photoisomerization quantum yields, in contradiction with the experiments, and call for a revision of the model.

Floquet Study of Quantum Control of the Cis-Trans Photoisomerization of Rhodopsin.

Understanding how to control reaction dynamics of polyatomic systems by using ultrafast laser technology is a fundamental challenge of great technological interest. Here, we report a Floquet theoretical study of the effect of light-induced potentials on the ultrafast cis-trans photoisomerization dynamics of rhodopsin. The Floquet Hamiltonian involves an empirical 3-state 25-mode model with frequencies and excited-state gradients parametrized to reproduce the rhodopsin electronic vertical excitation energy, the resonance Raman spectrum, and the photoisomerization time and efficiency as probed by ultrafast spectroscopy. We simulate the excited state relaxation dynamics using the time-dependent self-consistent field method, as described by a 3-state 2-mode nuclear wavepacket coupled to a Gaussian ansatz of 23 vibronic modes. We analyze the reaction time and product yield obtained with pulses of various widths and intensity profiles, defining 'dressed states' where the perturbational effect of the pulses is naturally decoupled along the different reaction channels. We find pulses that delay the excited-state photoisomerization for hundreds of femtoseconds, and we gain insights on the underlying control mechanisms. The reported findings provide understanding of quantum control, particularly valuable for the development of ultrafast optical switches based on visual pigments.

Quantum-classical modeling of rhodopsin photoisomerization

Quantum-classical modeling of rhodopsin photoisomerization

The phototransduction machinery in the rod outer segment has a strong efficacy gradient

Rod photoreceptors consist of an outer segment (OS) and an inner segment. Inside the OS a biochemical machinery transforms the rhodopsin photoisomerization into electrical signal. This machinery has been treated as and is thought to be homogenous with marginal inhomogeneities. To verify this assumption, we developed a methodology based on special tapered optical fibers (TOFs) to deliver highly localized light stimulations. By using these TOFs, specific regions of the rod OS could be stimulated with spots of light highly confined in space. As the TOF is moved from the OS base toward its tip, the amplitude of saturating and single photon responses decreases, demonstrating that the efficacy of the transduction machinery is not uniform and is 5-10 times higher at the base than at the tip. This gradient of efficacy of the transduction machinery is attributed to a progressive depletion of the phosphodiesterase along the rod OS. Moreover we demonstrate that, using restricted spots of light, the duration of the photoresponse along the OS does not increase linearly with the light intensity as with diffuse light.

Considering channel activation and inhibition

Considering channel activation and inhibition

Mechanistic Origin of the Vibrational Coherence Accompanying the Photoreaction of Biomimetic Molecular Switches

The coherent photoisomerization of a chromophore in condensed phase is a rare process in which light energy is funneled into specific molecular vibrations during electronic relaxation from the excited to the ground state. In this work, we employed ultrafast spectroscopy and computational methods to investigate the molecular origin of the coherent motion accompanying the photoisomerization of indanylidene-pyrroline (IP) molecular switches. UV/Vis femtosecond transient absorption gave evidence for an excited- and ground-state vibrational wave packet, which appears as a general feature of the IP compounds investigated. In close resemblance to the coherent photoisomerization of rhodopsin, the sudden onset of a far-red-detuned and rapidly blue-shifting photoproduct signature indicated that the population arriving on the electronic ground state after nonadiabatic decay through the conical intersection (CI) is still very focused in the form of a vibrational wave packet. Semiclassical trajectories were employed to investigate the reaction mechanism. Their analysis showed that coupled double-bond twisting and ring inversions, already populated during the excited-state reactive motion, induced periodic changes in π-conjugation that modulate the ground-state absorption after the non-adiabatic decay. This prediction further supports that the observed ground-state oscillation results from the reactive motion, which is in line with a biomimetic, coherent photoisomerization scenario. The IP compounds thus appear as a model system to investigate the mechanism of mode-selective photomechanical energy transduction. The presented mechanism opens new perspectives for energy transduction at the molecular level, with applications to the design of efficient molecular devices.

Photomechanical Responses in Drosophila Photoreceptors

Phototransduction in Drosophila microvillar photoreceptor cells is mediated by a G protein-activated phospholipase C (PLC). PLC hydrolyzes the minor membrane lipid phosphatidylinositol 4,5-bisphosphate (PIP(2)), leading by an unknown mechanism to activation of the prototypical transient receptor potential (TRP) and TRP-like (TRPL) channels. We found that light exposure evoked rapid PLC-mediated contractions of the photoreceptor cells and modulated the activity of mechanosensitive channels introduced into photoreceptor cells. Furthermore, photoreceptor light responses were facilitated by membrane stretch and were inhibited by amphipaths, which alter lipid bilayer properties. These results indicate that, by cleaving PIP(2), PLC generates rapid physical changes in the lipid bilayer that lead to contractions of the microvilli, and suggest that the resultant mechanical forces contribute to gating the light-sensitive channels.

QM/MM Trajectory Surface Hopping Approach to Photoisomerization of Rhodopsin and Isorhodopsin: The Origin of Faster and More Efficient Isomerization for Rhodopsin

The photoinduced cis-trans isomerization dynamics of rhodopsin and isorhodopsin are studied using a newly developed hybrid QM/MM trajectory surface hopping MD scheme based on the Zhu-Nakamura theory for nonadiabatic transitions. Rhodopsin and isorhodopsin have 11-cis and 9-cis forms of retinal as chromophore and the two proteins are isomerized to bathorhodopsin enclosing the all-trans form. The simulation reproduced faster and more efficient isomerization in rhodopsin than in isorhodopsin. In the excited state, rhodopsin shows a straightforward dynamics, whereas isorhodopsin dynamics is rather complicated and in a back-and-forth manner. The latter complicated dynamics would be mainly due to a narrow space near the active dihedral angle ═C8-C9═C10-C11═ (ϕ9) created by Thr 118 and Tyr 268 in opsin. Rhodopsin gives bathorhodopsin only while isorhodopsin yields a byproduct. The rigorous selectivity in rhodopsin would be another reason why rhodopsin is selected biologically. Comparison with our previous opsin-free investigations reveals that opsin tends to confine the twist of the active dihedral to only one direction and funnels transitions into the vicinity of minimum energy conical intersections (MECI). The twist-confinement totally blocks simultaneous twisting of ϕ9 and ϕ11 (═C10-C11═C12-C13═) and enhances the quantum yields. The opposite rotation of ϕ9 and ϕ11 ("wring-a-wet-towel" motion) takes place upon photoexcitation, which also does without opsin. The wring-a-wet-towel motion is dynamically enhanced in comparison with the one expected from locations of the MECI. The present simulation reveals that the Weiss-Warshel model for cis-trans photoisomerization is not applicable for rhodopsin because the branching ratio after transition is crucial.

Arrestin Translocation Is Stoichiometric to Rhodopsin Isomerization and Accelerated by Phototransduction in Drosophila Photoreceptors

Upon illumination, visual arrestin translocates from photoreceptor cell bodies to rhodopsin and membrane-rich photosensory compartments, vertebrate outer segments or invertebrate rhabdomeres, where it quenches activated rhodopsin. Both the mechanism and function of arrestin translocation are unresolved and controversial. In dark-adapted photoreceptors of the fruitfly Drosophila, confocal immunocytochemistry shows arrestin (Arr2) associated with distributed photoreceptor endomembranes. Immunocytochemistry and live imaging of GFP-tagged Arr2 demonstrate rapid reversible translocation to stimulated rhabdomeres in stoichiometric proportion to rhodopsin photoisomerization. Translocation is very rapid in normal photoreceptors (time constant <10 s) and can also be resolved in the time course of electroretinogram recordings. Genetic elimination of key phototransduction proteins, including phospholipase C (PLC), Gq, and the light-sensitive Ca2+-permeable TRP channels, slows translocation by 10- to 100-fold. Our results indicate that Arr2 translocation in Drosophila photoreceptors is driven by diffusion, but profoundly accelerated by phototransduction and Ca2+ influx.

Light Penetration and Photoisomerization in Rhodopsin studied by Numerical Simulations and Double-Quantum Solid-State NMR Spectroscopy

The penetration of light into optically thick samples containing the G-protein-coupled receptor rhodopsin is studied by numerical finite-element simulations and double-quantum solid-state NMR experiments. Illumination with white light leads to the generation of the active bathorhodopsin photostate in the outer layer of the sample but generates a large amount of the side product, isorhodopsin, in the sample interior. The overall yield of bathorhodopsin is improved by using monochromatic 420 nm illumination and by mixing the sample with transparent glass beads. The implications of these findings on the interpretation of previously published rhodopsin NMR data are discussed.

Coherent processes in formation of primary products of rhodopsin photolysis

194 The absorption of the light quantum by rhodopsin (R 498 ), a visual pigment, results in the cis → trans isomerization of its chromophore, 11cis -retinal, over approximately 200 fs [1, 2] with a quantum yield of 0.65 [3]. The physical meaning of this fast and efficient photoreaction is in the maximal utilization of the absorbed light quantum energy for isomerization of 11cis -retinal rather than its dissipation into heat or loss as fluorescence. The isomerization of 11cis -retinal leads to conformational rearrangement in the protein moiety of rhodopsin molecule; as a result, rhodopsin as a G-protein binding receptor acquires the ability to interact with G-protein (transducin), thereby initiating phototransduction in the photoreceptor cell. The dynamics of the 11cis -retinal photoisomerization in rhodopsin has been studied by laser kinetic spectroscopy with a femtosecond resolution [1, 2, 4–6]. It has been demonstrated that the first reaction product, photorhodopsin (photo 570 ), is formed over 200 fs to transform into the next product, bathorhodopsin (batho 535 ), in 2–3 ps [1, 4, 5]. The oscillations of the absorption signals of photo 570 , batho 535 , and initial state R 498 , observed during several picoseconds after excitation, have been also discovered. A high rate and quantum yield of 11cis -retinal photoisomerization in rhodopsin as well as the oscillations of absorption signals of the reaction products suggested that this is a coherent reaction. This means that the photoreaction proceeds via passing of the system through nonstationary oscillation states. In terms of the model of two states [1, 7], when two potential energy surfaces (PESs) are involved in photochemical reaction of rhodopsin (Fig. 1), the corresponding processes can be represented as follows. The action of a femtosecond excitation pulse results in a coherent population of certain oscillatory states in the electron-excited molecule, the set of which is named a coherent oscillatory wave packet (hereinafter, wave packet). The coherent photoreaction can be represented as the movement of wave packet over PES S 1 of the excited state along the reaction coordinate. The reaction products are formed Coherent Processes in Formation of Primary Products of Rhodopsin Photolysis

Protein Assistance in the Photoisomerization of Rhodopsin and 9-cis-RhodopsinInsights from Experiment and Theory

The artificial isomeric 9-cis form of retinal has been investigated within the chromophore binding pocket of rhodopsin by x-ray crystallography and quantum-mechanics with particular emphasis on the geometry imposed on the chromophore by the protein. Comparison with the native 11-cis-rhodopsin structure reveals similarities and differences providing valuable insights into the photophysical properties of the visual pigment.

Purification of visual arrestin from squid photoreceptors and characterization of arrestin interaction with rhodopsin and rhodopsin kinase

Invertebrate visual signal transduction involves photoisomerization of rhodopsin, activating a guanine nucleotide binding protein (G protein) of the G(q) class, iG(q), which stimulates a phospholipase C, increasing intracellular Ca2+. Arrestin binding to photoactivated rhodopsin is a key mechanism of desensitization. We have previously reported the cloning of a retina-specific arrestin cDNA from Loligo pealei displaying 56-64% sequence similarity to other reported arrestin sequences. Here, we report the purification of the 55-kDa squid visual arrestin. Purified squid visual arrestin is able to inhibit light-activated GTPase activity dose-dependently in arrestin-depleted rhabdomeric membranes and associate with the membrane in a light-dependent manner. Membrane association can be partially inhibited by inositol 1,2,3,4,5,6-hexakisphosphate (IP6), a soluble analog of the membrane lipid phosphatidylinositol 3,4,5-triphosphate. In reconstitution assays, we demonstrate arrestin phosphorylation by squid rhodopsin kinase, a novel function among the G protein-coupled receptor kinase family. Phosphorylation of purified arrestin requires squid rhodopsin kinase, membranes, light-activation, and the presence of Ca2+. This is the first large-scale purification of an invertebrate arrestin and biochemical demonstration of arrestin function in the invertebrate visual system.