Origin of Spectral Tuning in Rhodopsin—It Is Not the Binding Pocket

Color tuning in short wavelength-sensitive human and mouse visual pigments: ab initio quantum mechanics/molecular mechanics studies

We have investigated the protonation state and photoabsorption spectrum of Schiff-base (SB) nitrogen bound 11-cis-retinal in human blue and mouse UV cone visual pigments as well as in bovine rhodopsin by hybrid quantum mechanical/molecular mechanical (QM/MM) calculations. We have employed both multireference (MRCISD+Q, MR-SORCI+Q, and MR-DDCI2+Q) and single reference (TD-B3LYP and RI-CC2) QM methods. The calculated ground-state and vertical excitation energies show that UV-sensitive pigments have deprotonated SB nitrogen, while violet-sensitive pigments have protonated SB nitrogen, in agreement with some indirect experimental evidence. A significant blue shift of the absorption maxima of violet-sensitive pigments relative to rhodopsins arises from the increase in bond length alternation of the polyene chain of 11-cis-retinal induced by polarizing fields of these pigments. The main counterion is Glu113 in both violet-sensitive vertebrate pigments and bovine rhodopsin. Neither Glu113 nor the remaining pigment has a significant influence on the first excitation energy of 11-cis-retinal in the UV-sensitive pigments that have deprotonated SB nitrogen. There is no charge transfer between the SB and beta-ionone terminals of 11-cis-retinal in the ground and first excited states.

Electrostatic control of the photoisomerization efficiency and optical properties in visual pigments: on the role of counterion quenching

Hybrid QM(CASPT2//CASSCF/6-31G*)/MM(Amber) computations have been used to map the photoisomerization path of the retinal chromophore in Rhodopsin and explore the reasons behind the photoactivity efficiency and spectral control in the visual pigments. It is shown that while the electrostatic environment plays a central role in properly tuning the optical properties of the chromophore, it is also critical in biasing the ultrafast photochemical event: it controls the slope of the photoisomerization channel as well as the accessibility of the S(1)/S(0) crossing space triggering the ultrafast decay. The roles of the E113 counterion, the E181 residue, and the other amino acids of the protein pocket are explicitly analyzed: it appears that counterion quenching by the protein environment plays a key role in setting up the chromophore's optical properties and its photochemical efficiency. A unified scenario is presented that discloses the relationship between spectroscopic and mechanistic properties in rhodopsins and allows us to draw a solid mechanism for spectral tuning in color vision pigments: a tunable counterion shielding appears as the elective mechanism for L<-->M spectral modulation, while a retinal conformational control must dictate S absorption. Finally, it is suggested that this model may contribute to shed new light into mutations-related vision deficiencies that opens innovative perspectives for experimental biomolecular investigations in this field.

Color-changing mutation in the E-F loop of proteorhodopsin

It is usually assumed that only amino acids located near the retinal chromophore are responsible for color tuning of rhodopsins. However, we recently found that replacement of Ala178 with Arg in the E-F loop of proteorhodopsin (PR), an archaeal-type rhodopsin in marine bacteria, shifts the lambda(max) from 525 to 545 nm at neutral pH [Yoshitsugu, M., Shibata, M., Ikeda, D., Furutani, Y., and Kandori, H. (2008) Angew. Chem., Int. Ed. 47, 3923-3926]. Since the location of Ala178 is distant from the retinal chromophore (approximately 25 A), the molecular mechanism of the unusual mutation effect on color tuning is intriguing. Here we studied this mechanism by using additional mutations and some analytical methods. Introduction of Arg into the corresponding amino acid in bacteriorhodopsin (BR, M163R mutant) does not change the absorption spectra, indicating that the effect is specific to PR. Introduction of Arg into the A-B or C-D loop yields little (3 nm) or no color change, respectively. T177R and P180R mutants exhibited absorption spectra identical to that of the wild type, while N176R and S179R mutants exhibit lambda(max) values of 528 and 535 nm, respectively. Therefore, the observed color change is position-specific, being fully effective at position 178 and half-effective at position 179. Salt affects the absorption spectra of wild-type and A178R PR similarly. FTIR spectroscopy at 77 K indicated similar chromophore structures for wild-type and A178R PR, and A178R PR pumps protons normally. We infer that the E-F loop has a unique structure in PR and the mutation of Ala178 disrupts the structure that includes the transmembrane region, leading to the observed changes in color and pK(a).

Water-mediated spectral shifts in rhodopsin and bathorhodopsin

The role of water molecules in spectral tuning of proteins has been left largely unexplored. This topic is important because changing hydrogen bond patterns during the activation process may lead to spectral shifts which can be of diagnostic value for the underlying structures. Arguments put forward in this article are based on spectral shift calculations of the rhodopsin and bathorhodopsin chromophore due to wat2a and 2b in the presence and absence of the counterion and of the amino acids lining the rhodopsin binding pocket. They show, among others, that a single water molecule can shift the absorbance by up to 0.1 eV or 34 nm depending on the environment of the chromophore.

Spectral tuning in visual pigments: an ONIOM(QM:MM) study on bovine rhodopsin and its mutants

We have investigated geometries and excitation energies of bovine rhodopsin and some of its mutants by hybrid quantum mechanical/molecular mechanical (QM/MM) calculations in ONIOM scheme, employing B3LYP and BLYP density functionals as well as DFTB method for the QM part and AMBER force field for the MM part. QM/MM geometries of the protonated Schiff-base 11- cis-retinal with B3LYP and DFTB are very similar to each other. TD-B3LYP/MM excitation energy calculations reproduce the experimental absorption maximum of 500 nm in the presence of native rhodopsin environment and predict spectral shifts due to mutations within 10 nm, whereas TD-BLYP/MM excitation energies have red-shift error of at least 50 nm. In the wild-type rhodopsin, Glu113 shifts the first excitation energy to blue and accounts for most of the shift found. Other amino acids individually contribute to the first excitation energy but their net effect is small. The electronic polarization effect is essential for reproducing experimental bond length alternation along the polyene chain in protonated Schiff-base retinal, which correlates with the computed first excitation energy. It also corrects the excitation energies and spectral shifts in mutants, more effectively for deprotonated Schiff-base retinal than for the protonated form. The protonation state and conformation of mutated residues affect electronic spectrum significantly. The present QM/MM calculations estimate not only the experimental excitation energies but also the source of spectral shifts in mutants.

Quantum mechanical/molecular mechanical studies on spectral tuning mechanisms of visual pigments and other photoactive proteins

The protein environments surrounding the retinal tune electronic absorption maximum from 350 to 630 nm. Hybrid quantum mechanical/molecular mechanical (QM/MM) methods can be used in calculating excitation energies of retinal in its native protein environments and in studying the molecular basis of spectral tuning. We hereby review recent QM/MM results on the phototransduction of bovine rhodopsin, bacteriorhodopsin, sensory rhodopsin II, nonretinal photoactive yellow protein and their mutants.

An opsin shift in rhodopsin: retinal S0-S1 excitation in protein, in solution, and in the gas phase

We considered a series of model systems for treating the photoabsorption of the 11-cis retinal chromophore in the protonated Schiff-base form in vacuum, solutions, and the protein environment. A high computational level, including the quantum mechanical-molecular mechanical (QM/MM) approach for solution and protein was utilized in simulations. The S0-S1 excitation energies in quantum subsystems were evaluated by means of an augmented version of the multiconfigurational quasidegenerate perturbation theory (aug-MCQDPT2) with the ground-state geometry parameters optimized in the density functional theory PBE0/cc-pVDZ approximation. The computed positions of absorption bands lambdamax, 599(g), 448(s), and 515(p) nm for the gas phase, solution, and protein, respectively, are in excellent agreement with the corresponding experimental data, 610(g), 445(s), and 500(p) nm. Such consistency provides a support for the formulated qualitative conclusions on the role of the chromophore geometry, environmental electrostatic field, and the counterion in different media. An essentially nonplanar geometry conformation of the chromophore group in the region of the C14-C15 bond was obtained for the protein, in particular, owing to the presence of the neighboring charged amino acid residue Glu181. Nonplanarity of the C14-C15 bond region along with the influence of the negatively charged counterions Glu181 and Glu113 are found to be important to reproduce the spectroscopic features of retinal chromophore inside the Rh cavity. Furthermore, the protein field is responsible for the largest bond-order decrease at the C11-C12 double bond upon excitation, which may be the reason for the 11-cis photoisomerization specificity.

The Tautomeric Half-reaction of BphD, a C-C Bond Hydrolase

BphD of Burkholderia xenovorans LB400 catalyzes an unusual C-C bond hydrolysis of 2-hydroxy-6-oxo-6-phenylhexa-2,4-dienoic acid (HOPDA) to afford benzoic acid and 2-hydroxy-2,4-pentadienoic acid (HPD). An enol-keto tautomerization has been proposed to precede hydrolysis via a gem-diol intermediate. The role of the canonical catalytic triad (Ser-112, His-265, Asp-237) in mediating these two half-reactions remains unclear. We previously reported that the BphD-catalyzed hydrolysis of HOPDA (lambda(max) is 434 nm for the free enolate) proceeds via an unidentified intermediate with a red-shifted absorption spectrum (lambda(max) is 492 nm) (Horsman, G. P., Ke, J., Dai, S., Seah, S. Y. K., Bolin, J. T., and Eltis, L. D. (2006) Biochemistry 45, 11071-11086). Here we demonstrate that the S112A variant generates and traps a similar intermediate (lambda(max) is 506 nm) with a similar rate, 1/tau approximately 500 s(-1). The crystal structure of the S112A:HOPDA complex at 1.8-A resolution identified this intermediate as the keto tautomer, (E)-2,6-dioxo-6-phenyl-hex-3-enoate. This keto tautomer did not accumulate in either the H265A or the S112A/H265A double variants, indicating that His-265 catalyzes tautomerization. Consistent with this role, the wild type and S112A enzymes catalyzed tautomerization of the product HPD, whereas H265A variants did not. This study thus identifies a keto intermediate, and demonstrates that the catalytic triad histidine catalyzes the tautomerization half-reaction, expanding the role of this residue from its purely hydrolytic function in other serine hydrolases. Finally, the S112A:HOPDA crystal structure is more consistent with hydrolysis occurring via an acyl-enzyme intermediate than a gem-diol intermediate as solvent molecules have poor access to C6, and the closest ordered water is 7 A away.