Retinal counterion switch mechanism in vision evaluated by molecular simulations

Photoisomerization of the retinylidene chromophore of rhodopsin is the starting point in the vision cascade. A counterion switch mechanism that stabilizes the retinal protonated Schiff base (PSB) has been proposed to be an essential step in rhodopsin activation. On the basis of vibrational and UV-visible spectroscopy, two counterion switch models have emerged. In the first model, the PSB is stabilized by Glu181 in the meta I state, while in the most recent proposal, it is stabilized by Glu113 as well as Glu181. We assess these models by conducting a pair of microsecond scale, all-atom molecular dynamics simulations of rhodopsin embedded in a 99-lipid bilayer of SDPC, SDPE, and cholesterol (2:2:1 ratio) varying the starting protonation state of Glu181. Theoretical simulations gave different orientations of retinal for the two counterion switch mechanisms, which were used to simulate experimental 2H NMR spectra for the C5, C9, and C13 methyl groups. Comparison of the simulated 2H NMR spectra with experimental data supports the complex-counterion mechanism. Hence, our results indicate that Glu113 and Glu181 stabilize the retinal PSB in the meta I state prior to activation of rhodopsin.     

Protein-induced bonding perturbation of the rhodopsin chromophore detected by double-quantum solid-state NMR

We have obtained carbon-carbon bond length data for the functional retinylidene chromophore of rhodopsin, with a spatial resolution of 3 pm. The very high resolution was obtained by performing double-quantum solid-state NMR on a set of noncrystalline isotopically labelled bovine rhodopsin samples. We detected localized perturbations of the carbon-carbon bond lengths of the retinylidene chromophore. The observations are consistent with a model in which the positive charge of the protonated Schiff base penetrates into the polyene chain and partially concentrates around the C13 position. This coincides with the proximity of a water molecule located between the glutamate-181 and serine-186 residues of the second extracellular loop, which is folded back into the transmembrane region. These measurements support the hypothesis that the polar residues of the second extracellular loop and the associated water molecule assist the rapid selective photoisomerization of the retinylidene chromophore by stabilizing a partial positive charge in the center of the polyene chain.     

牛视紫红质蛋白质中视黄醛的活性位点

视紫红质蛋白是一个跨膜蛋白, 视黄醛(RET)在该蛋白中的活性结合位点涉及到视觉过程机理, 与一些眼科疾病病理有关. 基于牛视紫红质蛋白1U19的蛋白质晶体结构数据, 采用密度泛函理论的B3LYP方法计算RET-Lys296残基与视黄醛分子周围半径为0.6 nm的空间范围30个氨基酸残基相互作用和结合能. 数值显示1U19蛋白中的残基Glu113、Glu181和Glu122是质子化的RET-Lys296残基的活性结合位点, 结合能分别为-333.38、-205.67和-194.56 kJ·mol-1. 这些氨基酸残基带有一个负电荷, 与质子化的RET-Lys296残基发生强烈的离子静电相互作用. 另外几个残基Ala292、Cys187、Phe293、Pro291以及Trp265等与质子化RET-Lys296残基也有相互吸引作用. 当RET-Lys296残基非质子化, 上述相互作用消失, 促使视黄醛分子与视蛋白分离. 研究发现残基Glu113和Glu181周围各自有一个结晶水分子通过双氢键形式起着稳定作用.     

Comparative Analysis of GPCR Crystal Structures†

The phototransduction cascade is perhaps the best understood model system for G protein‐coupled receptor (GPCR) signaling. Phototransduction links the absorption of a single photon of light to a decrease in cytosolic cGMP. Depletion of the cGMP pool induces closure of cGMP‐gated cation channels resulting in the hyperpolarization of photoreceptor cells and consequently a neuronal response. Many biochemical and both low‐ and high‐resolution structural approaches have been utilized to increase our understanding of rhodopsin, the key molecule of this signaling cascade. Rhodopsin, a member of the GPCR or seven‐transmembrane spanning receptor superfamily, is composed of a chromophore, 11‐cis‐retinal that is covalently bound by a protonated Schiff base linkage to the apo‐protein opsin at Lys²⁹⁶ (in bovine opsin). Upon absorption of a photon, isomerization of the chromophore to an all‐trans‐retinylidene conformation induces changes in the rhodopsin structure, ultimately converting it from an inactive to an activated state. This state allows it to activate the heterotrimeric G protein, transducin, by triggering nucleotide exchange. To fully understand the structural and functional aspects of rhodopsin it is necessary to critically examine crystal structures of its different photointermediates. In this review we summarize recent progress on the structure and activation of rhodopsin in the context of other GPCR structures.     

The effect of the protein environment on the structure and charge distribution of the retinal Schiff base in bacteriorhodopsin

The effects of the amino acid side chains of the binding pocket of bacteriorhodopsin (bR) and of a water molecule on the structure of the retinal Schiff base have been studied using Becke3LYP/6-31G* level of density functional theory. A model protonated Schiff base structure including six conjugated double bonds and methyl substituents was optimized in the presence of several amino acid side chains and of a water molecule, separately. The Schiff base structure was also calculated in the form of a neutral species. At each optimized complex geometry the atomic charges of the model Schiff base were calculated using Mulliken population analysis. In agreement with previously proposed counterion(s) of the protonated retinal Schiff base in bR, the results show that Asp₈₅ and Asp₂₁₂, which are present in the form of negatively charged groups, have significantly large effects on the structure and electronic configuration of both unprotonated and protonated model Schiff bases. The presence of a water molecule in the vicinity of the Schiff base demonstrates significant effects which are comparable to those of aspartate groups. Other side chains studied did not show any significant effect in this direction. Apart from the aspartate groups and the water molecule, in none of the other complexes studied are the atomic charges and the bond alternation of the model Schiff base significantly influenced by the presence of the neighboring amino acids. Received: 24 March 1998 / Accepted: 3 September 1998 / Published online: 10 December 1998     

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.     

A molecular spring for vision

Light absorption by the visual pigment rhodopsin leads to vision via a complex signal transduction pathway that is initiated by the ultrafast and highly efficient photoreaction of its chromophore, the retinal protonated Schiff base (RPSB). Here, we investigate this reaction in real time by means of unrestrained molecular dynamics simulations of the protein in a membrane mimetic environment, treating the chromophore at the density functional theory level. We demonstrate that a highly strained all-trans RPSB is formed starting from the 11-cis configuration (dark state) within approximately 100 fs by a minor rearrangement of the nuclei under preservation of the saltbridge with Glu113 and virtually no deformation of the binding pocket. Hence, the initial step of vision can be understood as the compression of a molecular spring by a minor change of the nuclear coordinates. This spring can then release its strain by altering the protein environment.     

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.     

ON THE MECHANISM OF WAVELENGTH REGULATION IN VISUAL PIGMENTS

Abstract— The contributions of different factors that might be responsible for the 500 nm absorption maximum of bovine rhodopsin are evaluated in detail. These include: (1) electrostatic interactions between the chromophore and a charged amino acid on the apoprotein; (2) exciton interactions with aromatic amino acids; (3) twisting about single bonds which have considerable double bond character; (4) weak interactions between the Schiff base and a putative counter-ion. Analysis of these mechanisms in terms of theoretical and experimental results suggests that(2–4) are all capable of contributing to the protein induced spectral shifts. However, the "external point charge" model proposed previously, i.e. mechanism (1), appears to provide the crucial interaction. In this paper, the theoretical basis for this model is discussed in detail. The model is briefly evaluated in light of the amino-acid sequence of bovine rhodopsin and possible implications for other visual pigments are considered.     

The photoreaction of vacuum-dried rhodopsin at low temperature: evidence for charge stabilization by water

The photoreaction of vacuum-dried rhodopsin was monitored by UV-visible absorption spectroscopy. The results indicate that in dry rhodopsin, isorhodopsin and lumirhodopsin a protonation equilibrium exists between the protonated and the non-protonated Schiff base. On hydration the water stabilizes the protonated forms. In metarhodopsin-I the protein itself is able to stabilize the protonated Schiff base. The direct involvement of water in the retinal binding site was demonstrated by measuring the rhodopsin-bathorhodopsin FTIR difference spectra of rhodopsin hydrated with H2O and H2(18)O. The results are discussed with respect to the problem of charge stabilization and energy storage.     

Modulating rhodopsin receptor activation by altering the pKa of the retinal Schiff base

The visual pigment rhodopsin is a seven-transmembrane (7-TM) G protein-coupled receptor (GPCR). Activation of rhodopsin involves two pH-dependent steps: proton uptake at a conserved cytoplasmic motif between TM helices 3 and 6, and disruption of a salt bridge between a protonated Schiff base (PSB) and its carboxylate counterion in the transmembrane core of the receptor. Formation of an artificial pigment with a retinal chromophore fluorinated at C14 decreases the intrinsic pKa of the PSB and thereby destabilizes this salt bridge. Using Fourier transform infrared difference and UV-visible spectroscopy, we characterized the pH-dependent equilibrium between the active photoproduct Meta II and its inactive precursor, Meta I, in the 14-fluoro (14-F) analogue pigment. The 14-F chromophore decreases the enthalpy change of the Meta I-to-Meta II transition and shifts the Meta I/Meta II equilibrium toward Meta II. Combining C14 fluorination with deletion of the retinal beta-ionone ring to form a 14-F acyclic artificial pigment uncouples disruption of the Schiff base salt bridge from transition to Meta II and in particular from the cytoplasmic proton uptake reaction, as confirmed by combining the 14-F acyclic chromophore with the E134Q mutant. The 14-F acyclic analogue formed a stable Meta I state with a deprotonated Schiff base and an at least partially protonated protein counterion. The combination of retinal modification and site-directed mutagenesis reveals that disruption of the protonated Schiff base salt bridge is the most important step thermodynamically in the transition from Meta I to Meta II. This finding is particularly important since deprotonation of the retinal PSB is known to precede the transition to the active state in rhodopsin activation and is consistent with models of agonist-dependent activation of other GPCRs.     

Effect of charge distribution on electrostatic chromophore—protein interactions in Bacteriorhodopsin

Charge distributions of a protonated and unprotonated Schiff base model compound are determined using different quantum chemical methods. After fitting the model molecule onto the protonated retinal Schiff base in Bacteriorhodopsin, electrostatic interaction energies between the model molecule and protein are calculated. Interaction energies as well as the calculated pK_1/2 values of the model molecule are shown to depend considerably on the chosen charge distribution. Electrostatic potential derived partial charges determined at different ab initio levels reveal interaction energies between the model molecule and nearby residues such as ARG-82, ASP-85, and ASP-212, which are relatively method independent. Consequently, such charge distributions also result in pK_1/2 values for the model molecule that are very similar. Larger deviations in the electrostatic interaction energies, however, are found in the case of charge distributions derived according to the Mulliken population analysis. Nevertheless, some sets of Mulliken derived partial charges predicted pK_1/2 values for the model molecule that are close to those determined with electrostatic potential derived partial charges. This agreement, however, is only achieved because the individual errors of the contributing terms are approximately compensated. The use of the extended atom model is shown to be problematic. Although potential derived charges can correctly describe electrostatic interaction energies, they fail to predict pK_1/2 values. On the basis of the present investigation a new set of partial charges for the protonated and unprotonated retinal Schiff base is proposed to be used in molecular dynamics simulations and electrostatics calculations. © 1997 by John Wiley & Sons, Inc.     

Primary events in dim light vision: a chemical and spectroscopic approach toward understanding protein/chromophore interactions in rhodopsin

The visual pigment rhodopsin (bovine) is a 40 kDa protein consisting of 348 amino acids, and is a prototypical member of the subfamily A of G protein-coupled receptors (GPCRs). This remarkably efficient light-activated protein (quantum yield = 0.67) binds the chromophore 11-cis-retinal covalently by attachment to Lys296 through a protonated Schiff base. The 11-cis geometry of the retinylidene chromophore keeps the partially active opsin protein locked in its inactive state (inverse agonist). Several retinal analogs with defined configurations and stereochemistry have been incorporated into the apoprotein to give rhodopsin analogs. These incorporation results along with the spectroscopic properties of the rhodopsin analogs clarify the mode of entry of the chromophore into the apoprotein and the biologically relevant conformation of the chromophore in the rhodopsin binding site. In addition, difference UV, CD, and photoaffinity labeling studies with a 3-diazo-4-oxo analog of 11-cis-retinal have been used to chart the movement of the retinylidene chromophore through the various intermediate stages of visual transduction.     

STRUCTURE OF HYPSORHODOPSIN: ANALYSIS BY FOURIER TRANSFORM INFRARED SPECTROSCOPY AT 10 K

Abstract— Vibrational bands of hypsorhodopsin in the difference Fourier transform infrared spectra were identified as the bands which arose after formation of isorhodopsin by successive irradiations of bovine rhodopsin at 10 K with >500 nm light, and also as the bands disappeared upon conversion to bathorhodopsin by warming. The chromophore bands were assigned by the bands which shifted upon deuterium substitution of the polyene chain of the retinylidene chromophore. The presence of chromophore bands which shift by D₂O exchange clearly shows that the Schiff base chromophore of hypsorhodopsin is protonated. The amide I bands and several other protein bands of hypsorhodopsin appeared at the same frequencies as those of bathorhodopsin, but they are different from those of rhodopsin and isorhodopsin. Furthermore, like bathorhodopsin, hypsorhodopsin displays the C_l—H out-of-plane bending mode which is weakly coupled with C₁₂-_-–H out-of-plane mode. These facts show that hypsorhodopsin has a chromophore conformation and chromophore-opsin interaction more similar to bathorhodopsin than to rhodopsin and isorhodopsin.     

Rhodopsin deactivation is affected by mutations of Tyrl91

The 9-methyl group of 11-cis retinal is important in the efficient formation of the active conformation of rhodopsin, Meta II. Here, Tyrl91 rhodopsin mutants were generated because of its proximity to that methyl group in the dark structure. If photoactivation results in movement of the 9-methyl group toward Tyrl91, the steric interactions involved with activation and/or deactivation might not be as tightly coupled in mutant proteins with smaller amino acids at that position. Tyrl91 mutations have no effect on the dark pigment. However, after photobleaching, the lifetime of Meta II is shorter; Meta II decays quickly into two inactive species: (1) a Meta III or Meta III-like species and (2) opsin and free retinal. The Meta III-like fraction maintains the covalent Schiff base linkage of the chromophore much longer than the wild type. On the other hand, the fast chromophore release is similar to cone pigments. Taken together, the data suggest that the role of the 9-methyl group after photo-isomerization is not only to form Meta II efficiently, but also to maintain its active conformation and allow for the timely hydrolysis of the Schiff base. Perturbation of this interaction effects changes in the hydrolysis of the Schiff base and for the case of the Y191A mutation the folded structure of the protein after photobleaching.     

Engineered Functional Recovery of Microbial Rhodopsin Without Retinal‐Binding Lysine

Definition of rhodopsin is the retinal‐binding membrane protein with the Schiff base linkage at a lysine on the 7^th transmembrane helix. However, ~ 600 microbial rhodopsins lack retinal‐binding lysine at the corresponding position (Rh‐noK) among ~ 5500 known microbial rhodopsins, suggesting that Rh‐noK has each functional role without chromophore. Here, we report successful functional recovery of Rh‐noK. Two Rh‐noKs from bacteria were heterologously expressed in Escherichia coli, which exhibited no color. When retinal‐binding lysine was introduced, one of them gained visible color. Additional mutation of the Schiff base counterion further gained proton‐pumping activity. Successful engineered functional recovery such as visible color and proton‐pump activity suggests that the Rh‐noK protein forms a characteristic structure of microbial rhodopsins.     

RESONANCE RAMAN STUDIES OF BACTERIORHODOPSIN ANALOGUES

Abstract— We present the results of resonance Raman measurements on a series of bacteriorhodopsin (bR) analogues formed from synthetic retinals which have replaced the native chromophore in the active site. Specifically, 5,6-dihydro-bR, 13-desmethyl-bR, 10-methyl-bR, 14-methyl-bR, and 10.14-dimethyl-bR have been studied. All five analogues bind and form Schiff base retinal-apoprotein linkages. While the Schiff base linkages of 5,6-dihydro-bR, 13-desmethyl-bR, and 10-methyl-bR are protonated, like the native chromophore, the 14-methyl-bR, and 10,14-dimethyl-bR Schiff bases are unprotonated. These results suggest that the binding site of bacteriorhodopsin near the Schiff base moiety is different from that of rhodopsin. The protonated Schiff base -C=NH- stretching frequency of 5.6-dihydro-bR lies at 1660 cm^-1 which is unusually high for a bacteriorhodopsin based pigment. The downward shift upon deuteration is 16 cm^-1, essentially identical to that measured for bacteriorhodopsin. This and the other analogue results strongly reinforce our previous arguments that the Schiff base stretching frequency is determined in large part by two factors, the C=N force constant and the stretch interaction with C=N-H bend. On the other hand, the deuterium isotope effect is determined primarily by the stretch-bend interaction.     

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.     

Structural forms of green fluorescent protein by quantum mechanics/molecular mechanics calculations

The molecular modeling of structural forms of the green fluorescent protein (GFP) with the Ser65Thr single-site mutation was performed by the quantum mechanics/molecular mechanics (QM/MM) method. Two model systems were constructed based on the crystallographic structure from the Protein Data Bank (PDB entry code 1EMA.) The model systems differ in the initial protonation state of the side chain of the amino acid residue Glu222 near the chromophore. The atomic coordinates of the protein macromolecule corresponding to the equilibrium geometric configurations were determined by total energy minimization using the QM/MM method within the density functional theory approximation PBE0/cc-pVDZ for the quantum subsystem that consists of the chromophore, a water molecule, and the side chains of Arg96, Glu222, and Ser205, and with the parameters of the AMBER force field for the molecular mechanics subsystem. In the analysis of the results, particular attention was given to the hydrogen bond redistribution in the chromophore-containing region of the protein caused by a change in the protonation state of the chromophore. The results obtained from the model containing the initially protonated side chain of Glu222 suggest a new interpretation of the photophysical processes in the green fluorescent protein.     

Electron density deformations provide new insights into the spectral shift of rhodopsins

Spectral shifts of rhodopsin, which are related to variations of the electron distribution in 11‐cis‐retinal, are investigated here using the method of deformed atoms in molecules. We found that systems carrying the M207R and S186W mutations display large perturbations of the π‐conjugated system with respect to wild‐type rhodopsins. These changes agree with the predicted behavior of the bond length alternation (BLA) and the blue shifts of vertical excitation energies of these systems. The effect of the planarity of the central and Schiff‐base regions of retinal chain on the electronic structure of the chromophore is also investigated. By establishing nonlinear polynomial relations between BLA, chain distortions, and vertical excitation energies, we are also able to provide a semiquantitative approach for the understanding of the mechanisms regulating spectral shifts in rhodopsin and its mutants. © 2013 Wiley Periodicals, Inc.