Combined QM/MM study of the opsin shift in bacteriorhodopsin

Combined quantum mechanical and molecular mechanical (QM/MM) calculations and molecular dynamics simulations of bacteriorhodopsin (bR) in the membrane matrix have been carried out to determine the factors that make significant contributions to the opsin shift. We found that both solvation and interactions with the protein significantly shifts the absorption maximum of the retinal protonated Schiff base, but the effects are much more pronounced in polar solvents such as methanol, acetonitrile, and water than in the protein environment. The differential solvatochromic shifts of PSB in methanol and in bR leads to a bathochromic shift of about 1800 cm(-1). Because the combined QM/MM configuration interaction calculation is essentially a point charge model, this contribution is attributed to the extended point-charge model of Honig and Nakanishi. The incorporation of retinal in bR is accompanied by a change in retinal conformation from the 6-s-cis form in solution to the 6-s-trans configuration in bR. The extension of the pi-conjugated system further increases the red-shift by 2400 cm(-1). The remaining factors are due to the change in dispersion interactions. Using an estimate of about 1000 cm(-1) in the dispersion contribution by Houjou et al., we obtained a theoretical opsin shift of 5200 cm(-1) in bR, which is in excellent agreement with the experimental value of 5100 cm(-1). Structural analysis of the PSB binding site revealed the specific interactions that make contributions to the observed opsin shift. The combined QM/MM method used in the present study provides an opportunity to accurately model the photoisomerization and proton transfer reactions in bR.     

Color vision: "OH-site" rule for seeing red and green

Eyes gather information, and color forms an extremely important component of the information, more so in the case of animals to forage and navigate within their immediate environment. By using the ONIOM (QM/MM) (ONIOM = our own N-layer integrated molecular orbital plus molecular mechanics) method, we report a comprehensive theoretical analysis of the structure and molecular mechanism of spectral tuning of monkey red- and green-sensitive visual pigments. We show that interaction of retinal with three hydroxyl-bearing amino acids near the β-ionone ring part of the retinal in opsin, A164S, F261Y, and A269T, increases the electron delocalization, decreases the bond length alternation, and leads to variation in the wavelength of maximal absorbance of the retinal in the red- and green-sensitive visual pigments. On the basis of the analysis, we propose the "OH-site" rule for seeing red and green. This rule is also shown to account for the spectral shifts obtained from hydroxyl-bearing amino acids near the Schiff base in different visual pigments: at site 292 (A292S, A292Y, and A292T) in bovine and at site 111 (Y111) in squid opsins. Therefore, the OH-site rule is shown to be site-specific and not pigment-specific and thus can be used for tracking spectral shifts in any visual pigment.     

Rhodopsin regeneration is accelerated via noncovalent 11-cis retinal-opsin complex--a role of retinal binding pocket of opsin

The regeneration of bovine rhodopsin from its apoprotein opsin and the prosthetic group 11-cis retinal involves the formation of a retinylidene Schiff base with the epsilon-amino group of the active lysine residue of opsin. The pH dependence of a Schiff base formation in solution follows a typical bell-shaped profile because of the pH dependence of the formation and the following dehydration of a 1-aminoethanol intermediate. Unexpectedly, however, we find that the formation of rhodopsin from 11-cis retinal and opsin does not depend on pH over a wide pH range. These results are interpreted by the Matsumoto and Yoshizawa (Nature 258 [1975] 523) model of rhodopsin regeneration in which the 11-cis retinal chromophore binds first to opsin through the beta-ionone ring, followed by the slow formation of the retinylidene Schiff base in a restricted space. We find the second-order rate constant of the rhodopsin formation is 6100+/-300 mol(-1) s(-1) at 25 degrees C over the pH range 5-10. The second-order rate constant is much greater than that of a model Schiff base in solution by a factor of more than 10(7). A previous report by Pajares and Rando (J Biol Chem 264 [1989] 6804) suggests that the lysyl epsilon-NH(2) group of opsin is protonated when the beta-ionone ring binding site is unoccupied. The acceleration of the Schiff base formation in rhodopsin is explained by stabilization of the deprotonated form of the lysyl epsilon-NH(2) group which might be induced when the beta-ionone ring binding site is occupied through the noncovalent binding of 11-cis retinal to opsin at the initial stage of rhodopsin regeneration, followed by the proximity and orientation effect rendered by the formation of noncovalent 11-cis retinal-opsin complex.     

QM/MM study of dehydro and dihydro β-ionone retinal analogues in squid and bovine rhodopsins: implications for vision in salamander rhodopsin

Visual pigment rhodopsin provides a decisive crossing point for interaction between organisms and environment. Naturally occurring visual pigments contain only PSB11 and 3,4-dehydro-PSB11 as chromophores. Therefore, the ability of visual opsin to discriminate between the retinal geometries is investigated by means of QM/MM incorporation of PSB11, 6-s-cis and 6-s-trans forms of 3,4-dehydro-PSB11, and 3,4-dehydro-5,6-dihydro-PSB11 and 5,6-dihydro-PSB11 analogues into squid and bovine rhodopsin environments. The analogue-protein interaction reveals the binding site of squid rhodopsin to be malleable and ductile, while that of bovine rhodopsin is rigid and stiff. On the basis of these studies, a tentative model of the salamander rhodopsin binding site is also proposed.     

THE CONFORMATIONAL ANALYSIS AND PHOTOISOMERIZATION OF RETINOCHROME ANALOGS WITH POLYENALS

Abstract— 3, 7-Dimethyl-2, 4, 6, 8, 10-dodecapentaenal was synthesized for reconstitution of the retinochrome analog. Its opsin shift was 1000 cm ¹ smaller than that of native retinochrome, whose chromophore contains the same number of double bonds. The conformational change from 6-s-trans to 6-s-cis , as figured in a retinal molecule, plays an important role in the formation of the retinochrome analog, based on the estimation of opsin shifts for retinal analogs locked in the 6-s conformation. Thus the conformation of the 6–7 single bond in the native retinochrome was suggested to be 6 -cis . Analysis of the circular dichroic spectra of retinochrome analogs revealed that the 6-s conformation is independent of the appearance of the β-band. The stereoselectivity in the photoisomerization of the retinal analogs by a retinochrome template depends on the hydrophobic binding in the region of the β-ionone ring.     

THE 'OPSIN SHIFT' IN BACTERIORHODOPSIN: STUDIES WITH ARTIFICIAL BACTERIORHODOPSINS

Abstract— The difference (in cm⁻¹) in absorption maxima between the protonated Schiff base of retinals and the pigment derived therefrom has been defined as the opsin shift. It represents the influence of the opsin binding site on the chromophore. The analysis of the opsin shifts of a series of dihydrobacteriorhodopsins has led to the external point-charge model, which in addition to a counter anion near the Schiff base ammonium, carries another negative charge in the vicinity of the β-ionone ring. This is in striking contrast to the external point-charge model proposed earlier for the bovine visual pigment. The absorption maxima of rhodopsins formed from bromo- and phenyl retinals support the two models. A retinal carrying a photoaffinity label has yielded a nonbleachable bacteriorhodopsin.     

Photochemistry of a retinal protonated schiff-base analogue mimicking the opsin shift of bacteriorhodopsin

A retinal Schiff base analogue which artificially mimics the protein-induced red shifting of absorption in bacteriorhodopsin (BR) has been investigated with femtosecond multichannel pump probe spectroscopy. The objective is to determine if the catalysis of retinal internal conversion in the native protein BR, which absorbs at 570 nm, is directly correlated with the protein-induced Stokes shifting of this absorption band otherwise known as the "opsin shift". Results demonstrate that the red shift afforded in the model system does not hasten internal conversion relative to that taking place in a free retinal-protonated Schiff base (RPSB) in methanol solution, and stimulated emission takes place with biexponential kinetics and characteristic timescales of approximately 2 and 10.5 ps. This shows that interactions between the prosthetic group and the protein that lead to the opsin shift in BR are not directly involved in reducing the excited-state lifetime by nearly an order of magnitude. A sub-picosecond phase of spectral evolution, analogues of which are detected in photoexcited retinal proteins and RPSBs in solution, is observed after excitation anywhere within the intense visible absorption band. It consists of a large and discontinuous spectral shift in excited-state absorption and is assigned to electronic relaxation between excited states, a scenario which might also be relevant to those systems as well. Finally, a transient excess bleach component that tunes with the excitation wavelength is detected in the data and tentatively assigned to inhomogeneous broadening in the ground state absorption band. Possible sources of such inhomogeneity and its relevance to native RPSB photochemistry are discussed.     

Exploring the molecular mechanism for color distinction in humans

We examine here the role of the red, green, and blue human opsin structures in modulating the absorption properties of 11-cis-retinal bonded to the protein via a protonated Schiff base (PSB). We built the three-dimensional structures of the human red, green, and blue opsins using homology modeling techniques with the crystal structure of bovine rhodopsin as the template. We then used quantum mechanics (QM) combined with molecular mechanics (MM) (denoted as QM/MM) techniques in conjunction with molecular dynamics to determine how the room temperature molecular structures of the three human color opsin proteins modulate the absorption frequency of the same bound 11-cis-retinal chromophore to account for the differences in the observed absorption spectra. We find that the conformational twisting of the 11-cis-retinal PSB plays an important role in the green to blue opsin shift, whereas the dipolar side chains in the binding pocket play a surprising role of red-shifting the blue opsin with respect to the green opsin, as a fine adjustment to the opsin shift. The dipolar side chains play a large role in the opsin shift from red to green.     

THE ORIGIN OF THE RED-SHIFTED ABSORPTION MAXIMUM OF THE M412 INTERMEDIATE IN THE BACTERIORHODOPSIN PHOTOCYCLE

The factors that red shift the absorption maximum of the retinal Schiff base chromophore in the M₄₁₂ intermediate of bacteriorhodopsin photocycle relative to absorption in solution were investigated using a series of artificial pigments and studies of model compounds in solution. The artificial pigments derived from retinal analogs that perturb chromophore-protein interactions in the vicinity of the ring moiety indicate that a considerable part of the red shift may originate from interactions in the vicinity of the Schiff base linkage. Studies with model compounds revealed that hydrogen bonding to the Schiff base moiety can significantly red shift the absorption maximum. Furthermore, it was demonstrated that although s-trans ring-chain planarity prevails in the M₄₁₂ intermediate it does not contribute significantly (only ca 750 cm⁻¹) to the opsin shift observed in M₄₁₂. It is suggested that in M₄₁₂, the Schiff base linkage is hydrogen bonded to bound water and/or protein residues inducing a considerable red shift in the absorption maximum of the retinal chromophore.     

Absorption studies of neutral retinal Schiff base chromophores

The neutral retinal Schiff base is connected to opsin in UV sensing pigments and in the blue-shifted meta-II signaling state of the rhodopsin photocycle. We have designed and synthesized two model systems for this neutral chromophore and have measured their gas-phase absorption spectra in the electrostatic storage ring ELISA with a photofragmentation technique. By comparison to the absorption spectrum of the protonated retinal Schiff base in vacuo, we found that the blue shift caused by deprotonation of the Schiff base is more than 200 nm. The absorption properties of the UV absorbing proteins are thus largely determined by the intrinsic properties of the chromophore. The effect of approaching a positive charge to the Schiff base was also studied, as well as the susceptibility of the protonated and unprotonated chromophores to experience spectral shifts in different solvents.     

CONFORMATIONAL ANALYSIS OF THE RHODOPSIN CHROMOPHORE USING BICYCLIC RETINAL ANALOGUES

Abstract— For investigation of the chromophore conformation around the trimethyl cyclohexene ring and of the origin of the induced β-circular dichroism band in rhodopsin, two C₆-C₇ single bond-fixed retinal analogues, 6s-cb- and 6s-trans-locked bicyclic retinals (6 and 7, respectively) were synthesized and incorporated into bovine opsin in CHAPS-PC mixture. 6s-cb- and 6s-tram-Locked rhodopsin analogues (8 and 9 ) with A max at 539 and 545 nm, respectively, were formed. Interestingly, both 8 and 9 displayed α- and β-circular dichroism bands. The ellipticity of α-bands are similar in each other, while the β-band of 8 was about three times stronger than that of 9. Irradiation of 6s-trans-locked rhodopsin, 9, in the presence of hyroxylamine, resulted in the formation of only one of the enantiomers of 6s-rrans-locked retinal oxime showing a positive circular dichroism signal at around 390 nm. This fact strongly suggests that the retinal binding site of rhodopsin shows a chiral discrimination. From these experimental results, the interactions between the trimethyl cyclohexene ring portion in the chromophore and the neighbouring protein moiety in the rhodopsin molecule are discussed.     

The Effective Conjugation Length Is Responsible for the Red/Green Spectral Tuning in the Cyanobacteriochrome Slr1393g3

The origin of the spectral shift from a red‐ to a green‐absorbing form in a cyanobacteriochrome, Slr1393g3, was identified by combined quantum mechanics/molecular mechanics simulations. This protein, related to classical phytochromes, carries the open‐chain tetrapyrrole chromophore phycocyanobilin. Our calculations reveal that the effective conjugation length in the chromophore becomes shorter upon conversion from the red to the green form. This is related to the planarity of the entire chromophore. A large distortion was found for the terminal pyrrole rings A and D; however, the D ring contributes more strongly to the photoproduct tuning, despite a larger change in the twist of the A ring. Our findings implicate that the D ring twist can be exploited to regulate the absorption of the photoproduct. Hence, mutations that affect the D ring twist can lead to rational tuning of the photoproduct absorption, allowing the tailoring of cyanobacteriochromes for biotechnological applications such as optogenetics and bioimaging.     

Counterion controlled photoisomerization of retinal chromophore models: a computational investigation

CASPT2//CASSCF photoisomerization path computations have been used to unveil the effects of an acetate counterion on the photochemistry of two retinal protonated Schiff base (PSB) models: the 2-cis-penta-2,4-dieniminium and the all-trans-epta-2,4,6-trieniminium cations. Different positions/orientations of the counterion have been investigated and related to (i) the spectral tuning and relative stability of the S0, S1, and S2 singlet states; (ii) the selection of the photochemically relevant excited state; (iii) the control of the radiationless decay and photoisomerization rates; and, finally, (iv) the control of the photoisomerization stereospecificity. A rationale for the results is given on the basis of a simple (electrostatic) qualitative model. We show that the model readily explains the computational results providing a qualitative explanation for different aspects of the experimentally observed "environment" dependent PSB photochemistry. Electrostatic effects likely involved in controlling retinal photoisomerization stereoselectivity in the protein are also discussed under the light of these results, and clues for a stereocontrolled electrostatically driven photochemical process are presented. These computations provide a rational basis for the formulation of a mechanistic model for photoisomerization electrostatic catalysis.     

Color tuning in rhodopsins: the mechanism for the spectral shift between bacteriorhodopsin and sensory rhodopsin II

The mechanism of color tuning in the rhodopsin family of proteins has been studied by comparing the optical properties of the light-driven proton pump bacteriorhodopsin (bR) and the light detector sensory rhodopsin II (sRII). Despite a high structural similarity, the maximal absorption is blue-shifted from 568 nm in bR to 497 nm in sRII. The molecular mechanism of this shift is still a matter of debate, and its clarification sheds light onto the general mechanisms of color tuning in retinal proteins. The calculations employ a combined quantum mechanical/molecular mechanical (QM/MM) technique, using a DFT-based method for ground state properties and the semiempirical OM2/MRCI method and ab initio SORCI method for excited state calculations. The high efficiency of the methodology has allowed us to study a wide variety of aspects including dynamical effects. The absorption shift as well as various mutation experiments and vibrational properties have been successfully reproduced. Our results indicate that several sources contribute to the spectral shift between bR and sRII. The main factors are the counterion region at the extracellular side of retinal and the amino acid composition of the binding pocket. Our analysis allows a distinction and identification of the different effects in detail and leads to a clear picture of the mechanism of color tuning, which is in good agreement with available experimental data.     

Raman spectroelectrochemical study of Meldola blue, adsorbed and electropolymerized at a gold electrode

Surface-enhanced Raman spectroscopy (SERS) was used to probe the structure of adsorbed and electropolymerized Meldola blue (MB) films on roughened gold surfaces in solutions with pH 1.0 and 7.0 by using 785 nm excitation wavelength. Spectral bands were assigned based on density functional theory (DFT) calculations at B3LYP/6-311+G(2d, p) level. The most characteristic band of the oxidized MB form was found to be the breathing vibration of the central ring containing heteroatoms at 596 cm(-1). Based on a red shift of bands assigned to vibrations of double C=N(C(2)H(6)) bonds and adjacent ring C=C bonds in surface spectra as compared with solution 1 it was suggested that polymerization and interaction with an electrode surface proceed through these moieties. The presence of out-of-plane bands in SERS spectra was attributed to "flat" or slightly "tilted" orientation of aromatic rings at the interface. Potential-dependent spectral changes were followed by SERS spectroscopy. Raman spectra of the reduced MB form were obtained in both pH 1.0 and pH 7.0 solutions by analysis of the potential-difference SERS spectra. Reduced MB form can be recognized by characteristic bands near 1620, 1574, 1374, and 1234 cm(-1). By comparing the intensities of 1637 cm(-1) (oxidized MB form) and 1374 cm(-1) (reduced MB form) bands in experimental spectra of polymerized MB in pH 1.0 solution, a reduction-induced decrease by factor of 7 was estimated. A similar tendency in intensity changes showed calculations indicating that this effect is associated with reduction-induced changes in the molecular structure of the dye.     

The role of the beta-ionone ring in the photochemical reaction of rhodopsin

Time-dependent density functional theory (TDDFT) calculations on the photoabsorption process of the 11-cis retinal protonated Schiff base (PSB) chromophore show that the Franck-Condon relaxation of the first excited state of the chromophore involves a torsional twist motion of the beta-ionone ring relative to the conjugated retinyl chain. For the ground state, the beta-ionone ring and the retinyl chain of the free retinal PSB chromophore form a -40 degrees dihedral angle as compared to -94 degrees for the first excited state. The double bonds of the retinal are shorter for the fully optimized structure of the excited state than for the ground state suggesting a higher cis-trans isomerization barrier for the excited state than for the ground state. According to the present TDDFT calculations, the excitation of the retinal PSB chromophore does not primarily lead to a reaction along the cis-trans torsional coordinate at the C11-C12 bond. The activation of the isomerization center seems to occur at a later stage of the photo reaction. The results obtained at the TDDFT level are supported by second-order M?ller-Plesset (MP2) and approximate singles and doubles-coupled cluster (CC2) calculations on retinal chromophore models; the MP2 and CC2 calculations yield for them qualitatively the same ground state and excited-state structures as obtained in the density functional theory and TDDFT calculations.     

Structure, spectroscopy, and spectral tuning of the gas-phase retinal chromophore: the beta-ionone "handle" and alkyl group effect

The low-lying singlet states (i.e. S0, S1, and S2) of the chromophore of rhodopsin, the protonated Schiff base of 11-cis-retinal (PSB11), and of its all-trans photoproduct have been studied in isolated conditions by using ab initio multiconfigurational second-order perturbation theory. The computed spectroscopic features include the vertical excitation, the band origin, and the fluorescence maximum of both isomers. On the basis of the S0-->S1 vertical excitation, the gas-phase absorption maximum of PSB11 is predicted to be 545 nm (2.28 eV). Thus, the predicted absorption maximum appears to be closer to that of the rhodopsin pigment (2.48 eV) and considerably red-shifted with respect to that measured in solution (2.82 eV in methanol). In addition, the absorption maxima associated with the blue, green, and red cone visual pigments are tentatively rationalized in terms of the spectral changes computed for PSB11 structures featuring differently twisted beta-ionone rings. More specifically, a blue-shifted absorption maximum is explained in terms of a large twisting of the beta-ionone ring (with respect to the main conjugated chain) in the visual S-cone (blue) pigment chromophore. In contrast, the chromophore of the visual L-cone (red) pigment is expected to have a nearly coplanar beta-ionone ring yielding a six double bond fully conjugated framework. Finally, the M-cone (green) chromophore is expected to feature a twisting angle between 10 and 60 degrees. The spectroscopic effects of the alkyl substituents on the PSB11 spectroscopic properties have also been investigated. It is found that they have a not negligible stabilizing effect on the S1-S0 energy gap (and, thus, cause a red shift of the absorption maximum) only when the double bond of the beta-ionone ring conjugates significantly with the rest of the conjugated chain.     

Butyl conformational reorganization as a possible explanation for the longitudinal flexibility of the binding site of bacteriorhodopsin. The azulene and C-22 retinoid analogs.

The UV-VIS absorption data of four bacteriorhodopsin (BR) analogs formed from azulene-retinals of varying polyene chain length show that the one-bond-shortened to one-bond-lengthened analogs possess comparable opsin shift values to that of BR. A two-bond-shortened analog exhibited a much smaller opsin shift. These data, combined with those reported for the C-22 retinal analog (Tokunaga et al., 1977, Biophys. J. 19, 191-198) were analyzed by molecular modelling and computer graphics in terms of a model where conformational flexibility of the appended butyl is the controlling factor in determining ease of pigment formation and protein/substrate interaction.     

Promoting vibrations in human purine nucleoside phosphorylase. A molecular dynamics and hybrid quantum mechanical/molecular mechanical study

Crystallographic studies of human purine nucleoside phosphorylase (hPNP) with several transition-state (TS) analogues in the immucillin family showed an unusual geometric arrangement of the atoms O-5', O-4', and O(P), the nucleophilic phosphate oxygen, lying in a close three-oxygen stack. These observations were corroborated by extensive experimental kinetic isotope effect analysis. We propose that protein-facilitated dynamic modes in hPNP cause this stack, centered on the ribosyl O-4' oxygen, to squeeze together and push electrons toward the purine ring, stabilizing the oxacarbenium character of the TS. As the N-ribosidic bond is cleaved during the reaction, the pK(a) values of N-7 and O-6 increase by the electron density expelled by the oxygen-stack compression toward the purine ring. Increased electron density in the purine ring improves electrostatic interactions with nearby residues and facilitates the abstraction of a proton from a solvent proton or an unidentified general acid, making the purine a better leaving group, and accelerating catalysis. Classical and mixed quantum/classical molecular dynamics (MD) simulations of the Michaelis complex of hPNP with the substrates guanosine and phosphate were performed to assess the existence of protein-promoting vibrations (PPVs). Analogous simulations were performed for the substrates in aqueous solution. In the catalytic site, the O-5', O-4', and O(P) oxygens vibrate at frequencies of ca. 125 and 465 cm(-1), as opposed to 285 cm(-1) in the absence of hPNP. The hybrid quantum mechanical/molecular mechanical method was used to assess whether this enzymatic vibration pushing the oxygens together is coupled to the reaction coordinate, and thus has a direct positive impact on catalysis. The potential energy surface for the phosphorolysis reaction for several snapshots taken from the classical MD simulation showed substantial differences in oxygen compression. Our calculations showed the existence of PPVs coupled to the reaction coordinate, which effect electronic alterations in the active site by pushing the three oxygen centers together in proximity, and accelerate substrate turnover in the phosphorolysis reaction catalyzed by hPNP.