Replacement of leucine-93 by alanine or threonine slows down the decay of the N and O intermediates in the photocycle of bacteriorhodopsin: implications for proton uptake and 13-cis-retinal----all-trans-retinal reisomerization

Facade-Based Bicelles as a New Tool for Production of Active Membrane Proteins in a Cell-Free System

Integral membrane proteins are important components of a cell. Their structural and functional studies require production of milligram amounts of proteins, which nowadays is not a routine process. Cell-free protein synthesis is a prospective approach to resolve this task. However, there are few known membrane mimetics that can be used to synthesize active membrane proteins in high amounts. Here, we present the application of commercially available "Facade" detergents for the production of active rhodopsin. We show that the yield of active protein in lipid bicelles containing Facade-EM, Facade-TEM, and Facade-EPC is several times higher than in the case of conventional bicelles with CHAPS and DHPC and is comparable to the yield in the presence of lipid-protein nanodiscs. Moreover, the effects of the lipid-to-detergent ratio, concentration of detergent in the feeding mixture, and lipid composition of the bicelles on the total, soluble, and active protein yields are discussed. We show that Facade-based bicelles represent a prospective membrane mimetic, available for the production of membrane proteins in a cell-free system.

My remembrances of H.G. Khorana: exploring the mechanism of bacteriorhodopsin with site-directed mutagenesis and FTIR difference spectroscopy

H.G. Khorana's seminal contributions to molecular biology are well-known. He also had a lesser known but still major influence on current application of advanced vibrational spectroscopic techniques such as FTIR difference spectroscopy to explore the mechanism of bacteriorhodopsin and other integral membrane proteins. In this review, I provide a personal perspective of my collaborative research and interactions with Gobind, from 1982 to 1995 when our groups published over 25 papers together which resulted in an early picture of key features of the bacteriorhodopsin proton pump mechanism. Much of this early work served as a blueprint for subsequent advances based on combining protein bioengineering and vibrational spectroscopic techniques to study integral membrane proteins.

Structural characterization of proton-pumping rhodopsin lacking a cytoplasmic proton donor residue by X-ray crystallography

DTG/DTS rhodopsin, which was named based on a three-residue motif (DTG or DTS) that is important for its function, is a light-driven proton-pumping microbial rhodopsin using a retinal chromophore. In contrast to other light-driven ion-pumping rhodopsins, DTG/DTS rhodopsin does not have a cytoplasmic proton donor residue, such as Asp, Glu, or Lys. Because of the lack of cytoplasmic proton donor residue, proton directly binds to the retinal chromophore from the cytoplasmic solvent. However, mutational experiments that showed the complicated effects of mutations were not able to clarify the roles played by each residue, and the detail of proton uptake pathway is unclear because of the lack of structural information. To understand the proton transport mechanism of DTG/DTS rhodopsin, here we report the three-dimensional structure of one of the DTG/DTS rhodopsins, PspR from Pseudomonas putida, by X-ray crystallography. We show that the structure of the cytoplasmic side of the protein is significantly different from that of bacteriorhodopsin, the best-characterized proton-pumping rhodopsin, and large cytoplasmic cavities were observed. We propose that these hydrophilic cytoplasmic cavities enable direct proton uptake from the cytoplasmic solvent without the need for a specialized cytoplasmic donor residue. The introduction of carboxylic residues homologous to the cytoplasmic donors in other proton-pumping rhodopsins resulted in higher pumping activity with less pH dependence, suggesting that DTG/DTS rhodopsins are advantageous for producing energy and avoiding intracellular alkalization in soil and plant-associated bacteria.

The Photoreaction of the Proton-Pumping Rhodopsin 1 From the Maize Pathogenic Basidiomycete Ustilago maydis

Microbial rhodopsins have recently been discovered in pathogenic fungi and have been postulated to be involved in signaling during the course of an infection. Here, we report on the spectroscopic characterization of a light-driven proton pump rhodopsin (UmRh1) from the smut pathogen Ustilago maydis, the causative agent of tumors in maize plants. Electrophysiology, time-resolved UV/Vis and vibrational spectroscopy indicate a pH-dependent photocycle. We also characterized the impact of the auxin hormone indole-3-acetic acid that was shown to influence the pump activity of UmRh1 on individual photocycle intermediates. A facile pumping activity test was established of UmRh1 expressed in Pichia pastoris cells, for probing proton pumping out of the living yeast cells during illumination. We show similarities and distinct differences to the well-known bacteriorhodopsin from archaea and discuss the putative role of UmRh1 in pathogenesis.

Crystal structure of schizorhodopsin reveals mechanism of inward proton pumping

Schizorhodopsins (SzRs), a new rhodopsin family identified in Asgard archaea, are phylogenetically located at an intermediate position between type-1 microbial rhodopsins and heliorhodopsins. SzRs work as light-driven inward H⁺ pumps as xenorhodopsins in bacteria. Although E81 plays an essential role in inward H⁺ release, the H⁺ is not metastably trapped in such a putative H⁺ acceptor, unlike the other H⁺ pumps. It remains elusive why SzR exhibits different kinetic behaviors in H⁺ release. Here, we report the crystal structure of SzR AM_5_00977 at 2.1 Å resolution. The SzR structure superimposes well on that of bacteriorhodopsin rather than heliorhodopsin, suggesting that SzRs are classified with type-1 rhodopsins. The structure-based mutagenesis study demonstrated that the residues N100 and V103 around the β-ionone ring are essential for color tuning in SzRs. The cytoplasmic parts of transmembrane helices 2, 6, and 7 are shorter than those in the other microbial rhodopsins, and thus E81 is located near the cytosol and easily exposed to the solvent by light-induced structural change. We propose a model of untrapped inward H⁺ release; H⁺ is released through the water-mediated transport network from the retinal Schiff base to the cytosol by the side of E81. Moreover, most residues on the H⁺ transport pathway are not conserved between SzRs and xenorhodopsins, suggesting that they have entirely different inward H⁺ release mechanisms.

Mechanism of absorption wavelength shifts in anion channelrhodopsin-1 mutants

Using a quantum mechanical/molecular mechanical approach, we show the mechanisms of how the protein environment of Guillardia theta anion channelrhodopsin-1 (GtACR1) can shift the absorption wavelength. The calculated absorption wavelengths for GtACR1 mutants, M105A, C133A, and C237A are in agreement with experimentally measured wavelengths. Among 192 mutant structures investigated, mutations at Thr101, Cys133, Pro208, and Cys237 are likely to increase the absorption wavelength. In particular, T101A GtACR1 was expressed in HEK293T cells. The measured absorption wavelength is 10 nm higher than that of wild type, consistent with the calculated wavelength. (i) Removal of a polar residue from the Schiff base moiety, (ii) addition of a polar or acidic residue to the β-ionone ring moiety, and (iii) addition of a bulky residue to increase the planarity of the β-ionone and Schiff base moieties are the basis of increasing the absorption wavelength.

Comparative Studies of the Fluorescence Properties of Microbial Rhodopsins: Spontaneous Emission Versus Photointermediate Fluorescence

Rhodopsins are seven-transmembrane photoreceptor proteins that bind to the retinal chromophore and have been utilized as a genetically encoded voltage indicator (GEVI). So far, archaerhodopsin-3 (AR3) has been successfully used as a GEVI, despite its low fluorescence intensity. We performed comparative and quantitative fluorescence analyses of 15 microbial rhodopsins to explore these highly fluorescent molecules and to clarify their fluorescence mechanism. These rhodopsins showed a wide range of fluorescence intensities in mouse hippocampal neurons. Some of them, GR, HwBR, IaNaR, MR, and NpHR, showed fluorescence intensities comparable with or higher than that of AR3, suggesting their potential for GEVIs. The fluorescence intensity in neurons correlated with that of the bright fluorescent photointermediate such as a Q-intermediate (R = 0.75), suggesting that the fluorescence in neurons originates from the fluorescence of the photointermediate. Our findings provide a crucial step for producing next-generation rhodopsin-based GEVIs.

Probing Structure and Reaction Dynamics of Proteins Using Time-Resolved Resonance Raman Spectroscopy

The mechanistic understanding of protein functions requires insight into the structural and reaction dynamics. To elucidate these processes, a variety of experimental approaches are employed. Among them, time-resolved (TR) resonance Raman (RR) is a particularly versatile tool to probe processes of proteins harboring cofactors with electronic transitions in the visible range, such as retinal or heme proteins. TR RR spectroscopy offers the advantage of simultaneously providing molecular structure and kinetic information. The various TR RR spectroscopic methods can cover a wide dynamic range down to the femtosecond time regime and have been employed in monitoring photoinduced reaction cascades, ligand binding and dissociation, electron transfer, enzymatic reactions, and protein un- and refolding. In this account, we review the achievements of TR RR spectroscopy of nearly 50 years of research in this field, which also illustrates how the role of TR RR spectroscopy in molecular life science has changed from the beginning until now. We outline the various methodological approaches and developments and point out current limitations and potential perspectives.

Analog Retinal Redshifts Visible Absorption of QuasAr Transmembrane Voltage Sensors into Near-infrared

Opsin‐based transmembrane voltage sensors (OTVSs) are increasingly important tools for neuroscience enabling neural function in complex brain circuits to be explored in live, behaving animals. However, the visible wavelengths required for fluorescence excitation of the current generation of OTVSs limit optogenetic imaging in the brain to depths of only a few mm due to the strong absorption and scattering of visible light by biological tissues. We report that substitution of the native A1 retinal chromophore of the widely used QuasAr1/2 OTVSs with the retinal analog MMAR containing a methylamino‐modified dimethylphenyl ring results in over a 100‐nm redshift of the maxima of the absorption and fluorescence emission bands to near 700 and 840 nm, respectively. FT‐Raman spectroscopy reveals that at pH 7 QuasAr1 with both the A1 and MMAR chromophores possess predominantly an all‐trans protonated Schiff base configuration with the MMAR chromophore exhibiting increased torsion of the polyene single‐/double‐bond system similar to the O‐intermediate of the BR photocycle. In contrast, the A1 and the MMAR chromophores of QuasAr2 exist partially in a 13‐cis PSB configuration. These results demonstrate that QuasArs containing the MMAR chromophore are attractive candidates for use as NIR‐OTVSs, especially for applications such as deep brain imaging.

Photocycle of Sensory Rhodopsin II from Halobacterium salinarum (HsSRII): Mutation of D103 Accelerates M Decay and Changes the Decay Pathway of a 13-cis O-like Species

Aspartic acid 103 (D103) of sensory rhodopsin II from Halobacterium salinarum (HsSRII, or also called phoborhodopsin) corresponds to D115 of bacteriorhodopsin (BR). This amino acid residue is functionally important in BR. This work reveals that a substitution of D103 with asparagine (D103N) or glutamic acid (D103E) can cause large changes in HsSRII photocycle. These changes include (1) shortened lifetime of the M intermediate in the following order: the wild-type > D103N > D103E; (2) altered decay pathway of a 13-cis O-like species. The 13-cis O-like species, tentatively named Px, was detected in HsSRII photocycle. Px appeared to undergo branched reactions at 0°C, leading to a recovery of the unphotolyzed state and formation of a metastable intermediate, named P370, that slowly decayed to the unphotolyzed state at room temperature. In wild-type HsSRII at 0°C, Px mainly decayed to the unphotolyzed state, and the decay reaction toward P370 was negligible. In mutant D103E at 0°C, Px decayed to P370, while the recovery of the unphotolyzed state became unobservable. In mutant D103N, the two reactions proceeded at comparable rates. Thus, D103 of HsSRII may play an important role in regulation of the photocycle of HsSRII.

Molecular details of the unique mechanism of chloride transport by a cyanobacterial rhodopsin

Microbial rhodopsins are well known as versatile and ubiquitous light-driven ion transporters and photosensors. While the proton transport mechanism has been studied in great detail, much less is known about various modes of anion transport. Until recently, only two main groups of light-driven anion pumps were known, archaeal halorhodopsins (HRs) and bacterial chloride pumps (known as ClRs or NTQs). Last year, another group of cyanobacterial anion pumps with a very distinct primary structure was reported. Here, we studied the chloride-transporting photocycle of a representative of this new group, Mastigocladopsis repens rhodopsin (MastR), using time-resolved spectroscopy in the infrared and visible ranges and site-directed mutagenesis. We found that, in accordance with its unique amino acid sequence containing many polar residues in the transmembrane region of the protein, its photocycle features a number of unusual molecular events not known for other anion-pumping rhodopsins. It appears that light-driven chloride ion transfers by MastR are coupled with translocation of protons and water molecules as well as perturbation of several polar sidechains. Of particular interest is transient deprotonation of Asp-85, homologous to the cytoplasmic proton donor of light-driven proton pumps (such as Asp-96 of bacteriorhodopsin), which may serve as a regulatory mechanism.

Energetics and dynamics of a light-driven sodium-pumping rhodopsin

The conversion of light energy into ion gradients across biological membranes is one of the most fundamental reactions in primary biological energy transduction. Recently, the structure of the first light-activated Na⁺ pump, Krokinobacter eikastus rhodopsin 2 (KR2), was resolved at atomic resolution [Kato HE, et al. (2015) Nature 521:48-53]. To elucidate its molecular mechanism for Na⁺ pumping, we perform here extensive classical and quantum molecular dynamics (MD) simulations of transient photocycle states. Our simulations show how the dynamics of key residues regulate water and ion access between the bulk and the buried light-triggered retinal site. We identify putative Na⁺ binding sites and show how protonation and conformational changes gate the ion through these sites toward the extracellular side. We further show by correlated ab initio quantum chemical calculations that the obtained putative photocycle intermediates are in close agreement with experimental transient optical spectroscopic data. The combined results of the ion translocation and gating mechanisms in KR2 may provide a basis for the rational design of novel light-driven ion pumps with optogenetic applications.

Structural basis of the green-blue color switching in proteorhodopsin as determined by NMR spectroscopy

Proteorhodopsins (PRs) found in marine microbes are the most abundant retinal-based photoreceptors on this planet. PR variants show high levels of environmental adaptation, as their colors are tuned to the optimal wavelength of available light. The two major green and blue subfamilies can be interconverted through a L/Q point mutation at position 105. Here we reveal the structural basis behind this intriguing color-tuning effect. High-field solid-state NMR spectroscopy was used to visualize structural changes within green PR directly within the lipid bilayer upon introduction of the green-blue L105Q mutation. The observed effects are localized within the binding pocket and close to retinal carbons C14 and C15. Subsequently, magic-angle spinning (MAS) NMR spectroscopy with sensitivity enhancement by dynamic nuclear polarization (DNP) was applied to determine precisely the retinal structure around C14-C15. Upon mutation, a significantly stretched C14-C15 bond, deshielding of C15, and a slight alteration of the retinal chain's out-of-plane twist was observed. The L105Q blue switch therefore acts locally on the retinal itself and induces a conjugation defect between the isomerization region and the imine linkage. Consequently, the S0-S1 energy gap increases, resulting in the observed blue shift. The distortion of the chromophore structure also offers an explanation for the elongated primary reaction detected by pump-probe spectroscopy, while chemical shift perturbations within the protein can be linked to the elongation of late-photocycle intermediates studied by flash photolysis. Besides resolving a long-standing problem, this study also demonstrates that the combination of data obtained from high-field and DNP-enhanced MAS NMR spectroscopy together with time-resolved optical spectroscopy enables powerful synergies for in-depth functional studies of membrane proteins.

A color-determining amino acid residue of proteorhodopsin

Proteorhodopsin (PR) is a light-driven proton pump found in marine bacteria. More than 1000 PRs are classified as blue-absorbing (λmax ∼ 490 nm) and green-absorbing (λmax ∼ 525 nm) PRs. The color determinant is known to be at position 105, where blue-absorbing and green-absorbing PRs possess Gln and Leu, respectively. This suggests hydrophobicity at position 105 plays a key role in color tuning. Here we successfully introduced 19 amino acid residues into position 105 of green-absorbing PR in the membrane environment and investigated the absorption properties. High-performance liquid chromatography analysis shows that the isomeric composition of the all-trans form is >70% for all mutants, indicating little influence of different isomers on color tuning. Absorption spectra of the wild-type and 19 mutant proteins were well-characterized by the pH-dependent equilibria of the protonated and deprotonated counterion (Asp97) of the Schiff base, whereas the λmax values of these two states and the pKa value differed significantly among mutants. Although Gln and Leu are hydrophilic and hydrophobic residues, respectively, the λmax values of the two states and the pKa value did not correlate with the hydropathy index of residues. In contrast, the λmax and pKa were correlated with the volume of residues, though Gln and Leu possess similar volumes. This observation concludes that the λmax and pKa of Asp97 are determined by local and specific interactions in the Schiff base moiety, in which the volume of the residue at position 105 is more influential than its hydrophobicity. We suggest that the hydrogen-bonding network in the Schiff base moiety plays a key role in the λmax and pKa of Asp97, and the hydrogen-bonding network is significantly perturbed by large amino acid residues but may be preserved by additional water molecule(s) for small amino acid residues at position 105.

Directed evolution of Gloeobacter violaceus rhodopsin spectral properties

Proton-pumping rhodopsins (PPRs) are photoactive retinal-binding proteins that transport ions across biological membranes in response to light. These proteins are interesting for light-harvesting applications in bioenergy production, in optogenetics applications in neuroscience, and as fluorescent sensors of membrane potential. Little is known, however, about how the protein sequence determines the considerable variation in spectral properties of PPRs from different biological niches or how to engineer these properties in a given PPR. Here we report a comprehensive study of amino acid substitutions in the retinal-binding pocket of Gloeobacter violaceus rhodopsin (GR) that tune its spectral properties. Directed evolution generated 70 GR variants with absorption maxima shifted by up to ±80nm, extending the protein's light absorption significantly beyond the range of known natural PPRs. While proton-pumping activity was disrupted in many of the spectrally shifted variants, we identified single tuning mutations that incurred blue and red shifts of 42nm and 22nm, respectively, that did not disrupt proton pumping. Blue-shifting mutations were distributed evenly along the retinal molecule while red-shifting mutations were clustered near the residue K257, which forms a covalent bond with retinal through a Schiff base linkage. Thirty eight of the identified tuning mutations are not found in known microbial rhodopsins. We discovered a subset of red-shifted GRs that exhibit high levels of fluorescence relative to the WT (wild-type) protein.

Schiff base switch II precedes the retinal thermal isomerization in the photocycle of bacteriorhodopsin

In bacteriorhodopsin, the order of molecular events that control the cytoplasmic or extracellular accessibility of the Schiff bases (SB) are not well understood. We use molecular dynamics simulations to study a process involved in the second accessibility switch of SB that occurs after its reprotonation in the N intermediate of the photocycle. We find that once protonated, the SB C15 = NZ bond switches from a cytoplasmic facing (13-cis, 15-anti) configuration to an extracellular facing (13-cis, 15-syn) configuration on the pico to nanosecond timescale. Significantly, rotation about the retinal’s C13 = C14 double bond is not observed. The dynamics of the isomeric state transitions of the protonated SB are strongly influenced by the surrounding charges and dielectric effects of other buried ions, particularly D96 and D212. Our simulations indicate that the thermal isomerization of retinal from 13-cis back to all-trans likely occurs independently from and after the SB C15 = NZ rotation in the N-to-O transition.

Probing a polar cluster in the retinal binding pocket of bacteriorhodopsin by a chemical design approach

Bacteriorhodopsin has a polar cluster of amino acids surrounding the retinal molecule, which is responsible for light harvesting to fuel proton pumping. From our previous studies, we have shown that threonine 90 is the pivotal amino acid in this polar cluster, both functionally and structurally. In an attempt to perform a phenotype rescue, we have chemically designed a retinal analogue molecule to compensate the drastic effects of the T90A mutation in bacteriorhodopsin. This analogue substitutes the methyl group at position C₁₃ of the retinal hydrocarbon chain by and ethyl group (20-methyl retinal). We have analyzed the effect of reconstituting the wild-type and the T90A mutant apoproteins with all-trans-retinal and its 20-methyl derivative (hereafter, 13-ethyl retinal). Biophysical characterization indicates that recovering the steric interaction between the residue 90 and retinal, eases the accommodation of the chromophore, however it is not enough for a complete phenotype rescue. The characterization of these chemically engineered chromoproteins provides further insight into the role of the hydrogen bond network and the steric interactions involving the retinal binding pocket in bacteriorhodopsin and other microbial sensory rhodopsins.

Ultra light-sensitive and fast neuronal activation with the Ca²+-permeable channelrhodopsin CatCh

The light-gated cation channel channelrhodopsin-2 (ChR2) has rapidly become an important tool in neuroscience, and its use is being considered in therapeutic interventions. Although wild-type and known variant ChR2s are able to drive light-activated spike trains, their use in potential clinical applications is limited by either low light sensitivity or slow channel kinetics. We present a new variant, calcium translocating channelrhodopsin (CatCh), which mediates an accelerated response time and a voltage response that is ~70-fold more light sensitive than that of wild-type ChR2. CatCh's superior properties stem from its enhanced Ca²(+) permeability. An increase in [Ca²(+)](i) elevates the internal surface potential, facilitating activation of voltage-gated Na(+) channels and indirectly increasing light sensitivity. Repolarization following light-stimulation is markedly accelerated by Ca²(+)-dependent BK channel activation. Our results demonstrate a previously unknown principle: shifting permeability from monovalent to divalent cations to increase sensitivity without compromising fast kinetics of neuronal activation. This paves the way for clinical use of light-gated channels.

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).

Long-distance proton transfer with a break in the bacteriorhodopsin active site

Bacteriorhodopsin is a proton-pumping membrane protein found in the plasma membrane of the archaeon Halobacterium salinarium. Light-induced isomerization of the retinal chromophore from all-trans to 13-cis leads to a sequence of five conformation-coupled proton transfer steps and the net transport of one proton from the cytoplasmic to the extracellular side of the membrane. The mechanism of the long-distance proton transfer from the primary acceptor Asp85 to the extracellular proton release group during the O --> bR is poorly understood. Experiments suggest that this long-distance transfer could involve a transient state [O] in which the proton resides on the intermediate carrier Asp212. To assess whether the transient protonation of Asp212 participates in the deprotonation of Asp85, we performed hybrid Quantum Mechanics/Molecular Mechanics proton transfer calculations using different protein structures and with different retinal geometries and active site water molecules. The structural models were assessed by computing UV-vis excitation energies and C=O vibrational frequencies. The results indicate that a transient [O] conformer with protonated Asp212 could indeed be sampled during the long-distance proton transfer to the proton release group. Our calculations suggest that, in the starting proton transfer state O, the retinal is strongly twisted and at least three water molecules are present in the active site.

Raman spectroscopy reveals direct chromophore interactions in the Leu/Gln105 spectral tuning switch of proteorhodopsins

Proteorhodopsins are an extensive family of photoactive membrane proteins found in proteobacteria distributed throughout the world's oceans which are often classified as green- or blue-absorbing (GPR and BPR, respectively) on the basis of their visible absorption maxima. GPR and BPR have significantly different properties including photocycle lifetimes and wavelength dependence on pH. Previous studies revealed that these different properties are correlated with a single residue, Leu105 in GPR and Gln105 in BPR, although the molecular basis for the different properties of GPR and BPR has not yet been elucidated. We have studied the unexcited states of GPR and BPR using resonance Raman spectroscopy which enhances almost exclusively chromophore vibrations. We find that both spectra are remarkably similar, indicating that the retinylidene structure of GPR and BPR are almost identical. However, the frequency of a band assigned to the retinal C13-methyl-rock vibration is shifted from 1006 cm (-1) in GPR to 1012 cm (-1) in BPR. A similar shift is observed in the GPR mutant L105Q indicating Leu and Gln residues interact differently with the retinal C13-methyl group. The environment of the Schiff base of GPR and BPR differ as indicated by differences in the H/D induced down-shift of the Schiff base vibration. Residues located in transmembrane helices (D-G) do not contribute to the observed differences in the protein-chromophore interaction between BPR and GPR based on the Raman spectra of chimeras. These results support a model whereby the substitution of the hydrophilic Gln105 in BPR with the smaller hydrophobic Leu105 in GPR directly alters the environment of both the retinal C13 group and the Schiff base.

Tryptophan 171 in Pharaonis phoborhodopsin (sensory rhodopsin II) interacts with the chromophore retinal and its substitution with alanine or threonine slowed down the decay of M- and O-intermediate

Pharaonis phoborhodopsin (ppR), also called pharaonis sensory rhodopsin II, NpSRII, is a photoreceptor for the photophobic response of Natronomonas pharaonis. Tryptophan 182 (W182) of bacteriorhodopsin (bR) is near the chromophore retinal and has been suggested to interact with retinal during the photoreaction and also to be involved in the hydrogen-bonding network around the retinal. W182 of bR is conserved in ppR as tryptophan 171 (W171). To elucidate whether W171 of ppR interacts with retinal during the photoreaction and/or is involved in the hydrogen-bonding network as in bR, we formed W171-substituted mutants of ppR, W171A and W171T. Our low-temperature spectroscopic study has revealed that the substitution of W171 to Ala or Thr resulted in the stabilization of M- and O-intermediates. The stability of M and absorption spectral changes during the M-decay were different depending on the substituted residue. These findings suggest that W171 in ppR interacts with retinal and the degree of the interaction depends on the substituted residues, which might be rate determining in the M-decay. In addition, the involvement of W171 in the hydrogen-bonding network is suggested by the O-decay. We also found that glycerol slowed the decay of M and not of O.

Crystal Structure of the L Intermediate of Bacteriorhodopsin: Evidence for Vertical Translocation of a Water Molecule during the Proton Pumping Cycle

For structural investigation of the L intermediate of bacteriorhodopsin, a 3D crystal belonging to the space group P622 was illuminated with green light at 160 K and subsequently with red light at 100 K. This yielded a approximately 1:4 mixture of the L intermediate and the ground-state. Diffraction data from such crystals were collected using a low flux of X-rays ( approximately 2 x 10(15) photons/mm2 per crystal), and their merged data were compared with those from unphotolyzed crystals. These structural data, together with our previous data, indicate that the retinal chromophore, which is largely twisted in the K-intermediate, takes a more planar 13-cis, 15-anti configuration in the L intermediate. This configurational change, which is accompanied by re-orientation of the Schiff base N-H bond towards the intracellular side, is coupled with a large rotation of the side-chain of an amino acid residue (Leu93) making contact with the C13 methyl group of retinal. Following these motions, a water molecule, at first hydrogen-bonded to the Schiff base and Asp85, is dragged to a space that is originally occupied by Leu93. Diffraction data from a crystal containing the M intermediate showed that this water molecule moves further towards the intracellular side in the L-to-M transition. It is very likely that detachment of this water molecule from the protonated Schiff base causes a significant decrease in the pKa of the Schiff base, thereby facilitating the proton transfer to Asp85. On the basis of these observations, we argue that the vertical movement of a water molecule in the K-to-L transition is a key event determining the directionality of proton translocation in the protein.

Proteorhodopsin genes are distributed among divergent marine bacterial taxa

Proteorhodopsin (PR) is a retinal-binding bacterial integral membrane protein that functions as a light-driven proton pump. The gene encoding this photoprotein was originally discovered on a large genome fragment derived from an uncultured marine gamma-proteobacterium of the SAR86 group. Subsequently, many variants of the PR gene have been detected in marine plankton, via PCR-based gene surveys. It has not been clear, however, whether these different PR genes are widely distributed among different bacterial groups, or whether they have a restricted taxonomic distribution. We report here comparative analyses of PR-bearing genomic fragments recovered directly from planktonic bacteria inhabiting the California coast, the central Pacific Ocean, and waters offshore the Antarctica Peninsula. Sequence analysis of an Antarctic genome fragment harboring PR (ANT32C12) revealed moderate conservation in gene order and identity, compared with a previously reported PR-containing genome fragment from a Monterey Bay gamma-proteobacterium (EBAC31A08). Outside the limited region of synteny shared between these clones, however, no significant DNA or protein identity was evident. Analysis of a third PR-containing genome fragment (HOT2C01) from the North Pacific subtropical gyre showed even more divergence from the gamma-proteobacterial PR-flanking region. Subsequent phylogenetic and comparative genomic analyses revealed that the Central North Pacific PR-containing genome fragment (HOT2C01) originated from a planktonic alpha-proteobacterium. These data indicate that PR genes are distributed among a variety of divergent marine bacterial taxa, including both alpha- and gamma-proteobacteria. Our analyses also demonstrate the utility of cultivation-independent comparative genomic approaches for assessing gene content and distribution in naturally occurring microbes.

Structural insights into the mechanism of proton pumping by bacteriorhodopsin

For over three decades, bacteriorhodopsin has served as a paradigm for the study of the mechanisms underlying ion pumping across biological membranes. It is perhaps among the simplest known ion pumps, which functions by converting light energy into an electrochemical gradient by pumping protons out of the cytoplasm. The combination of spectroscopic, biochemical and crystallographic studies on bacteriorhodopsin provides a unique opportunity to dissect the principal elements underlying the mechanism of transmembrane proton transport. Here, we provide a brief review of recent developments related to the determination of the structural changes during proton transport using crystallographic approaches. Taken together with previous spectroscopic and biochemical investigations, these studies allow the description of a detailed molecular mechanism of the main steps in vectorial proton transport by bacteriorhodopsin.

Late events in the photocycle of bacteriorhodopsin mutant L93A

In the photocycle of bacteriorhodopsin (bR) from Halobacterium salinarum mutant L93A, the O-intermediate accumulates and the cycling time is increased approximately 200 times. Nevertheless, under continuous illumination, the protein pumps protons at near wild-type rates. We excited the mutant L93A in purple membrane with single or triple laser flashes and quasicontinuous illumination, (i.e., light for a few seconds) and recorded proton release and uptake, electric signals, and absorbance changes. We found long-living, correlated, kinetic components in all three measurements, which-with exception of the absorbance changes-had not been seen in earlier investigations. At room temperature, the O-intermediate decays to bR in two transitions with rate constants of 350 and 1800 ms. Proton uptake from the cytoplasmic surface continues with similar kinetics until the bR state is reestablished. An analysis of the data from quasicontinuous illumination and multiple flash excitation led to the conclusion that acceleration of the photocycle in continuous light is due to excitation of the N-component in the fast N<-->O equilibrium, which is established at the beginning of the severe cycle slowdown. This conclusion was confirmed by an action spectrum.

Photoexcitation of the O-intermediate in bacteriorhodopsin mutant L93A

During the extended lifetime of the O-state in bacteriorhodopsin (bR) mutant L93A, two substates have been distinguished. The first O-intermediate (OI) is in rapid equilibrium with N and apparently still has a 13-cis chromophore. OI undergoes a photoreaction with a small absorbance change, positive charge transport in the pumping direction, and proton release and uptake. None of these effects was detected after photoexcitation of the late O (OII). The most likely interpretation of the effects seen is an accelerated return of the molecule from the OI- to the bR-state. However, with a lifetime approximately 140 ms, the reaction cannot account for the observed high pumping efficiency of the mutant under continuous illumination. We suggest that OII corresponds to the O-intermediate with a twisted all-trans chromophore seen in the photocycle of wild-type bR, where the 13-cis OI-intermediate under the usual conditions does not accumulate in easily detectable amounts and, therefore, has generally been overlooked. Both the OI- and OII-decays are apparently strongly inhibited in the mutant.

Diversification and spectral tuning in marine proteorhodopsins

Proteorhodopsins, ubiquitous retinylidene photoactive proton pumps, were recently discovered in the cosmopolitan uncultured SAR86 bacterial group in oceanic surface waters. Two related proteorhodopsin families were found that absorb light with different absorption maxima, 525 nm (green) and 490 nm (blue), and their distribution was shown to be stratified with depth. Using structural modeling comparisons and mutagenesis, we report here on a single amino acid residue at position 105 that functions as a spectral tuning switch and accounts for most of the spectral difference between the two pigment families. Furthermore, looking at natural environments, we found novel proteorhodopsin gene clusters spanning the range of 540-505 nm and containing changes in the same identified key switch residue leading to changes in their absorption maxima. The results suggest a simultaneous diversification of green proteorhodopsin and the new key switch variant pigments. Our observations demonstrate that this single-residue switch mechanism is the major determinant of proteorhodopsin wavelength regulation in natural marine environments.

Time-resolved FTIR studies of sensory rhodopsin II (NpSRII) from Natronobacterium pharaonis: implications for proton transport and receptor activation

The photocycle of the photophobic receptor from Natronobacterium pharaonis, NpSRII, is studied by static and time-resolved step-scan Fourier transform infrared spectroscopy. Both low-temperature static and time-resolved spectra resolve a K-like intermediate, and the corresponding spectra show little difference within the noise of the time-resolved data. As compared to intermediate K of bacteriorhodopsin, relatively large amide I bands indicate correspondingly larger distortions of the protein backbone. The time-resolved spectra identify an intermediate L-like state with surprisingly small additional molecular alterations. With the formation of intermediate M, the Schiff-base proton is transferred to the counterion Asp-75. This state is characterized by larger amide bands indicating larger distortions of the protein. We can identify a second M state that differs only in small-protein bands. Reisomerization of the chromophore to all-trans occurs with the formation of intermediate O. The accelerated decay of intermediate M caused by azide results in another red-shifted intermediate with a protonated Schiff base. The chromophore in this state, however, still has 13-cis geometry. Nevertheless, the reisomerization is still as slow as under the conditions without azide. The results are discussed with respect to mechanisms of the observed proton pumping and the possible roles of the intermediates in receptor activation.

Bacteriorhodopsin: a high-resolution structural view of vectorial proton transport

Recent 3-D structures of several intermediates in the photocycle of bacteriorhodopsin (bR) provide a detailed structural picture of this molecular proton pump in action. In this review, we describe the sequence of conformational changes of bR following the photoisomerization of its all-trans retinal chromophore, which is covalently bound via a protonated Schiff base to Lys216 in helix G, to a 13-cis configuration. The initial changes are localized near the protein's active site and a key water molecule is disordered. This water molecule serves as a keystone for the ground state of bR since, within the framework of the complex counter ion, it is important both for stabilizing the structure of the extracellular half of the protein, and for maintaining the high pK(a) of the Schiff base (the primary proton donor) and the low pK(a) of Asp85 (the primary proton acceptor). Subsequent structural rearrangements propagate out from the active site towards the extracellular half of the protein, with a local flex of helix C exaggerating an early movement of Asp85 towards the Schiff base, thereby facilitating proton transfer between these two groups. Other coupled rearrangements indicate the mechanism of proton release to the extracellular medium. On the cytoplasmic half of the protein, a local unwinding of helix G near the backbone of Lys216 provides sites for water molecules to order and define a pathway for the reprotonation of the Schiff base from Asp96 later in the photocycle. A steric clash of the photoisomerized retinal with Trp182 in helix F drives an outward tilt of the cytoplasmic half of this helix, opening the proton transport channel and enabling a proton to be taken up from the cytoplasm. Although bR is the first integral membrane protein to have its catalytic mechanism structurally characterized in detail, several key results were anticipated in advance of the structural model and the general framework for vectorial proton transport has, by and large, been preserved.

Environment around the chromophore in pharaonis phoborhodopsin: mutation analysis of the retinal binding site

Phoborhodopsin (pR or sensory rhodopsin II, sRII) and pharaonis phoborhodopsin (ppR or pharaonis sRII, psRII) have a unique absorption maximum (lambda(max)) compared with three other archaeal rhodopsins: lambda(max) of pR and ppR is approx. 500 nm and of others (e.g. bacteriorhodopsin, bR) is 560-590 nm. To determine the residue contributing to the opsin shift from ppR to bR, we constructed various ppR mutants, in which a single residue was substituted for a residue corresponding to that of bR. The residues mutated were those which differ from that of bR and locate within 5 A from the conjugated polyene chain of the chromophore or any methyl group of the polyene chain. The shifts of lambda(max) of all mutants were small, however. We constructed a mutant in which all residues which differ from those of bR in the retinal binding site were simultaneously substituted for those of bR, but the shift was only from 499 to 509 nm. Next, we constructed a mutant in which 10 residues located within 5 A from the polyene as described above were simultaneously substituted. Only 44% of the opsin shift (lambda(max) of 524 nm) from ppR to bR was obtained even when all amino acids around the chromophore were replaced by the same residues as bR. We therefore conclude that the structural factor is more important in accounting for the difference of lambda(max) between ppR and bR rather than amino acid substitutions. The possible structural factors are discussed.

Electrophysiological characterization of specific interactions between bacterial sensory rhodopsins and their transducers

The halobacterial phototaxis receptors sensory rhodopsin I and II (SRI, SRII) enable the bacteria to seek optimal light conditions for ion pumping by bacteriorhodopsin and/or halorhodopsin. The incoming signal is transferred across the plasma membrane by means of receptor-specific transducer proteins that bind tightly to their corresponding photoreceptors. To investigate the receptor/transducer interaction, advantage is taken of the observation that both SRI and SRII can function as proton pumps. SRI from Halobacterium salinarum, which triggers the positive phototaxis, the photophobic receptor SRII from Natronobacterium pharaonis (pSRII), as well as the mutant pSRII-F86D were expressed in Xenopus oocytes. Voltage-clamp studies confirm that SRI and pSRII function as light-driven, outwardly directed proton pumps with a much stronger voltage dependence than the ion pumps bacteriorhodopsin and halorhodopsin. Coexpression of SRI and pSRII-F86D with their corresponding transducers suppresses the proton transport, revealing a tight binding and specific interaction of the two proteins. These latter results may be exploited to further analyze the binding interaction of the photoreceptors with their downstream effectors.

Electrophysiological characterization of specific interactions between bacterial sensory rhodopsins and their transducers

The halobacterial phototaxis receptors sensory rhodopsin I and II (SRI, SRII) enable the bacteria to seek optimal light conditions for ion pumping by bacteriorhodopsin and/or halorhodopsin. The incoming signal is transferred across the plasma membrane by means of receptor-specific transducer proteins that bind tightly to their corresponding photoreceptors. To investigate the receptor/transducer interaction, advantage is taken of the observation that both SRI and SRII can function as proton pumps. SRI from Halobacterium salinarum, which triggers the positive phototaxis, the photophobic receptor SRII from Natronobacterium pharaonis (pSRII), as well as the mutant pSRII-F86D were expressed in Xenopus oocytes. Voltage-clamp studies confirm that SRI and pSRII function as light-driven, outwardly directed proton pumps with a much stronger voltage dependence than the ion pumps bacteriorhodopsin and halorhodopsin. Coexpression of SRI and pSRII-F86D with their corresponding transducers suppresses the proton transport, revealing a tight binding and specific interaction of the two proteins. These latter results may be exploited to further analyze the binding interaction of the photoreceptors with their downstream effectors.

Atomic resolution structures of bacteriorhodopsin photocycle intermediates: the role of discrete water molecules in the function of this light-driven ion pump

High-resolution X-ray crystallographic studies of bacteriorhodopsin have tremendously advanced our understanding of this light-driven ion pump during the last 2 years, and emphasized the crucial role of discrete internal water molecules in the pump cycle. In the extracellular region an extensive three-dimensional hydrogen-bonded network of protein residues and seven water molecules leads from the buried retinal Schiff base via water 402 and the initial proton acceptor Asp85 to the membrane surface. Near Lys216 where the retinal binds, transmembrane helix G contains a pi-bulge that causes a non-proline kink. The bulge is stabilized by hydrogen bonding of the main chain carbonyl groups of Ala215 and Lys216 with two buried water molecules located in the otherwise very hydrophobic region between the Schiff base and the proton donor Asp96 in the cytoplasmic region. The M intermediate trapped in the D96N mutant corresponds to a late M state in the transport cycle, after protonation of Asp85 and release of a proton to the extracellular membrane surface, but before reprotonation of the deprotonated retinal Schiff base. The M intermediate from the E204Q mutant corresponds to an earlier M, as in this mutant the Schiff base deprotonates without proton release. The structures of these two M states reveal progressive displacements of the retinal, main chain and side chains induced by photoisomerization of the retinal to 13-cis,15-anti, and an extensive rearrangement of the three-dimensional network of hydrogen-bonded residues and bound water that accounts for the changed pK(a)s of the Schiff base, Asp85, the proton release group and Asp96. The structure for the M state from E204Q suggests, moreover, that relaxation of the steric conflicts of the distorted 13-cis,15-anti retinal plays a critical role in the reprotonation of the Schiff base by Asp96. Two additional waters now connect Asp96 to the carbonyl of residue 216, in what appears to be the beginning of a hydrogen-bonded chain that would later extend to the retinal Schiff base. Based on the ground state and M intermediate structures, models of the molecular events in the early part of the photocycle are presented, including a novel model which proposes that bacteriorhodopsin pumps hydroxide (OH(-)) ions from the extracellular to the cytoplasmic side.

Crystallographic studies of the conformational changes that drive directional transmembrane ion movement in bacteriorhodopsin

Recent advances in the determination of the X-ray crystallographic structures of bacteriorhodopsin, and some of its photointermediates, reveal the nature of the linkage between the relaxation of electrostatic and steric conflicts at the retinal and events elsewhere in the protein. The transport cycle can be now understood in terms of specific and well-described displacements of hydrogen-bonded water, and main-chain and side-chain atoms, that lower the pK(a)s of the proton release group in the extracellular region and Asp-96 in the cytoplasmic region. Thus, local electrostatic conflict of the photoisomerized retinal with Asp-85 and Asp-212 causes deprotonation of the Schiff base, and results in a cascade of events culminating in proton release to the extracellular surface. Local steric conflict of the 13-methyl group with Trp-182 causes, in turn, a cascade of movements in the cytoplasmic region, and results in reprotonation of the Schiff base. Although numerous questions concerning the mechanism of each of these proton (or perhaps hydroxyl ion) transfers remain, the structural results provide a detailed molecular explanation for how the directionality of the ion transfers is determined by the configurational relaxation of the retinal.

Coupling photoisomerization of retinal to directional transport in bacteriorhodopsin

In order to understand how isomerization of the retinal drives unidirectional transmembrane ion transport in bacteriorhodopsin, we determined the atomic structures of the BR state and M photointermediate of the E204Q mutant, to 1.7 and 1.8 A resolution, respectively. Comparison of this M, in which proton release to the extracellular surface is blocked, with the previously determined M in the D96N mutant indicates that the changes in the extracellular region are initiated by changes in the electrostatic interactions of the retinal Schiff base with Asp85 and Asp212, but those on the cytoplasmic side originate from steric conflict of the 13-methyl retinal group with Trp182 and distortion of the pi-bulge of helix G. The structural changes suggest that protonation of Asp85 initiates a cascade of atomic displacements in the extracellular region that cause release of a proton to the surface. The progressive relaxation of the strained 13-cis retinal chain with deprotonated Schiff base, in turn, initiates atomic displacements in the cytoplasmic region that cause the intercalation of a hydrogen-bonded water molecule between Thr46 and Asp96. This accounts for the lowering of the pK(a) of Asp96, which then reprotonates the Schiff base via a newly formed chain of water molecules that is extending toward the Schiff base.

Two groups control light-induced Schiff base deprotonation and the proton affinity of Asp85 in the Arg82 his mutant of bacteriorhodopsin

Arg(82) is one of the four buried charged residues in the retinal binding pocket of bacteriorhodopsin (bR). Previous studies show that Arg(82) controls the pK(a)s of Asp(85) and the proton release group and is essential for fast light-induced proton release. To further investigate the role of Arg(82) in light-induced proton pumping, we replaced Arg(82) with histidine and studied the resulting pigment and its photochemical properties. The main pK(a) of the purple-to-blue transition (pK(a) of Asp(85)) is unusually low in R82H: 1.0 versus 2.6 in wild type (WT). At pH 3, the pigment is purple and shows light and dark adaptation, but almost no light-induced Schiff base deprotonation (formation of the M intermediate) is observed. As the pH is increased from 3 to 7 the M yield increases with pK(a) 4.5 to a value approximately 40% of that in the WT. A transition with a similar pK(a) is observed in the pH dependence of the rate constant of dark adaptation, k(da). These data can be explained, assuming that some group deprotonates with pK(a) 4.5, causing an increase in the pK(a) of Asp(85) and thus affecting k(da) and the yield of M. As the pH is increased from 7 to 10.5 there is a further 2.5-fold increase in the yield of M and a decrease in its rise time from 200 micros to 75 micros with pK(a) 9. 4. The chromophore absorption band undergoes a 4-nm red shift with a similar pK(a). We assume that at high pH, the proton release group deprotonates in the unphotolyzed pigment, causing a transformation of the pigment into a red-shifted "alkaline" form which has a faster rate of light-induced Schiff base deprotonation. The pH dependence of proton release shows that coupling between Asp(85) and the proton release group is weakened in R82H. The pK(a) of the proton release group in M is 7.2 (versus 5.8 in the WT). At pH < 7, most of the proton release occurs during O --> bR transition with tau approximately 45 ms. This transition is slowed in R82H, indicating that Arg(82) is important for the proton transfer from Asp(85) to the proton release group. A model describing the interaction of Asp(85) with two ionizable residues is proposed to describe the pH dependence of light-induced Schiff base deprotonation and proton release.

Structural changes in bacteriorhodopsin during ion transport at 2 angstrom resolution

Crystal structures of the Asp96 to Asn mutant of the light-driven proton pump bacteriorhodopsin and its M photointermediate produced by illumination at ambient temperature have been determined to 1.8 and 2.0 angstroms resolution, respectively. The trapped photoproduct corresponds to the late M state in the transport cycle-that is, after proton transfer to Asp85 and release of a proton to the extracellular membrane surface, but before reprotonation of the deprotonated retinal Schiff base. Its density map describes displacements of side chains near the retinal induced by its photoisomerization to 13-cis,15-anti and an extensive rearrangement of the three-dimensional network of hydrogen-bonded residues and bound water that accounts for the changed pKa values (where Ka is the acid constant) of the Schiff base and Asp85. The structural changes detected suggest the means for conserving energy at the active site and for ensuring the directionality of proton translocation.

Bacteriorhodopsin

Bacteriorhodopsin is a seven-transmembrane helical protein that contains all-trans retinal. In this light-driven pump, a reaction cycle initiated by photoisomerization to 13-cis causes translocation of a proton across the membrane. Local changes in the geometry of the protonated Schiff base and the proton acceptor Asp85, and the proton conductivities of the half channels that lead from this active site to the two membrane surfaces, interact so as to allow timely proton transfers that result in proton release on the extracellular side and proton uptake on the cytoplasmic one. The details of the steps in this photocycle, and the underlying principles that ensure unidirectionality of the movement of a proton across the protein, provide strong clues to how ion pumps function.

Understanding structure and function in the light-driven proton pump bacteriorhodopsin

The atomic structure of bacteriorhodopsin and the outlines of its proton transport mechanism are now available. Photoisomerization of the retinal in the chromophore creates a steric and electrostatic conflict at the retinal binding site. The free energy gain sets off a sequence of reactions in which directed proton transfers take place between the protonated retinal Schiff base, Asp-85, and Asp-96. These internal steps, and other proton transfers at and near the two aqueous interfaces, add up to the translocation of a proton from the cytoplasmic to the extracellular side of the membrane. Bound water plays a crucial role in proton conduction in both extracellular and cytoplasmic regions, but the means by which the protons move from site to site differ. Proton release to the extracellular surface is through interaction of a hydrogen-bonded chain of identified aspartic acid, arginine, water, and glutamic acid residues with Asp-85, while proton uptake from the cytoplasmic surface utilizes a single aspartic acid, Asp-96, whose protonation state appears to be regulated by the protein conformation dependent hydration of this region. The directionality of the translocation is ensured by the accessibility of the Schiff base to the extracellular and cytoplasmic directions after the retinal is photoisomerized, as well as the changing proton affinities of the acceptor Asp-85 and donor Asp-96.

Binding pathway of retinal to bacterio-opsin: a prediction by molecular dynamics simulations

Formation of bacteriorhodopsin (bR) from apoprotein and retinal has been studied experimentally, but the actual pathway, including the point of entry, is little understood. Molecular dynamics simulations provide a surprisingly clear prediction. A window between bR helices E and F in the transmembrane part of the protein can be identified as an entry point for retinal. Steered molecular dynamics, performed by applying a series of external forces in the range of 200-1000 pN over a period of 0.2 ns to retinal, allows one to extract this chromophore from bR once the Schiff base bond to Lys216 is cleaved. Extraction proceeds until the retinal tail forms a hydrogen bond network with Ala144, Met145, and Ser183 side groups lining the exit/entry window. The manipulation induces a distortion with a fitted root mean square deviation of coordinates (ignoring retinal, water, and hydrogen atoms) of less than 1.9 A by the time the retinal carbonyl reaches the protein surface. The forces needed to extract retinal are due to friction and do not indicate significant potential barriers. The simulations therefore suggest a pathway for the binding of retinal. Water molecules are found to play a crucial role in the binding process.