Structure of archaerhodopsin-2 at 1.8 Å resolution

Existence of two substates in the O intermediate of the bacteriorhodopsin photocycle

The proton pumping cycle of bacteriorhodopsin (bR) is initiated when the retinal chromophore with the 13-trans configuration is photo-isomerized into the 13-cis configuration. To understand the recovery processes of the initial retinal configuration that occur in the late stage of the photocycle, we have performed a comprehensive analysis of absorption kinetics data collected at various pH levels and at different salt concentrations. The result of analysis revealed the following features of the late stages of the trans photocycle. i) Two substates occur in the O intermediate. ii) The visible absorption band of the first substate (O1) appears at a much shorter wavelength than that of the late substate (O2). iii) O1 is in rapid equilibrium with the preceding state (N), but O1 becomes less stable than N when an ionizable residue (X₁) with a pKa value of 6.5 (in 2 M KCl) is deprotonated. iv) At a low pH and at a low salt concentration, the decay time constant of O2 is longer than those of the preceding states, but the relationship between these time constants is altered when the medium pH or the salt concentration is increased. On the basis of the present observations and previous studies on the structure of the chromophore in O, we suspect that the retinal chromophore in O1 takes on a distorted 13-cis configuration and the O1-to-O2 transition is accompanied by cis-to-trans isomerization about C13C14 bond.

QuasAr Odyssey: the origin of fluorescence and its voltage sensitivity in microbial rhodopsins

Rhodopsins had long been considered non-fluorescent until a peculiar voltage-sensitive fluorescence was reported for archaerhodopsin-3 (Arch3) derivatives. These proteins named QuasArs have been used for imaging membrane voltage changes in cell cultures and small animals. However due to the low fluorescence intensity, these constructs require use of much higher light intensity than other optogenetic tools. To develop the next generation of sensors, it is indispensable to first understand the molecular basis of the fluorescence and its modulation by the membrane voltage. Based on spectroscopic studies of fluorescent Arch3 derivatives, we propose a unique photo-reaction scheme with extended excited-state lifetimes and inefficient photoisomerization. Molecular dynamics simulations of Arch3, of the Arch3 fluorescent derivative Archon1, and of several its mutants have revealed different voltage-dependent changes of the hydrogen-bonding networks including the protonated retinal Schiff-base and adjacent residues. Experimental observations suggest that under negative voltage, these changes modulate retinal Schiff base deprotonation and promote a decrease in the populations of fluorescent species. Finally, we identified molecular constraints that further improve fluorescence quantum yield and voltage sensitivity.

The authors present an in-depth investigation of excited state dynamics and molecular mechanism of the voltage sensing in microbial rhodopsins. Using a combination of spectroscopic investigations and molecular dynamics simulations, the study proposes the voltage-modulated deprotonation of the chromophore as the key event in the voltage sensing. Thus, molecular constraints that may further improve the fluorescence quantum yield and the voltage sensitivity are presented.

Two substates in the O intermediate of the light-driven proton pump archaerhodopsin-2

The proton pumping cycle of archaerhodopsin-2 (aR2) was investigated over a wide pH range and at different salt concentrations. We have found that two substates, which are spectroscopically and kinetically distinguishable, occur in the O intermediate. The first O-intermediate (O1) absorbs maximumly at ~580 nm, whereas the late O-intermediate (O2) absorbs maximumly at 605 nm. At neutral pH, O1 is in rapid equilibrium with the N intermediate. When the medium pH is increased, O1 becomes less stable than N and, in proportion to the amount of O1 in the dynamic equilibrium between N and O1, the formation rate of O2 decreases. By contrast, the decay rate of O2 increases ~100 folds when the pH of a low-salt membrane suspension is increased from 5.5 to 7.5 or when the salt concentration is increased to 2 M KCl. Together with our recent study on two substates in the O intermediate of bacteriorhodopsin (bR), the present study suggests that the thermally activated re-isomerization of the retinylidene chromophore into the initial all-trans configuration takes place in the O1-to-O2 transition; that is, O1 contains a distorted 13-cis chromophore. It is also found that the pKa value of the key ionizable residue (Asp101^aR2, Asp96^bR) in the proton uptake channel is elevated in the O1 state of aR2 as compared to the O1 state of bR. This implies that the structural property of O1 in the aR2 photocycle can be investigated over a wide pH range.

Conserved hydrogen-bond motifs of membrane transporters and receptors

Membrane transporters and receptors often rely on conserved hydrogen bonds to assemble transient paths for ion transfer or long-distance conformational couplings. For transporters and receptors that use proton binding and proton transfer for function, inter-helical hydrogen bonds of titratable protein sidechains that could change protonation are of central interest to formulate hypotheses about reaction mechanisms. Knowledge of hydrogen bonds common at sites of potential interest for proton binding could thus inform and guide studies on functional mechanisms of protonation-coupled membrane proteins. Here we apply graph-theory approaches to identify hydrogen-bond motifs of carboxylate and histidine sidechains in a large data set of static membrane protein structures. We find that carboxylate-hydroxyl hydrogen bonds are present in numerous structures of the dataset, and can be part of more extended H-bond clusters that could be relevant to conformational coupling. Carboxylate-carboxyamide and imidazole-imidazole hydrogen bonds are represented in comparably fewer protein structures of the dataset. Atomistic simulations on two membrane transporters in lipid membranes suggest that many of the hydrogen bond motifs present in static protein structures tend to be robust, and can be part of larger hydrogen-bond clusters that recruit additional hydrogen bonds.

Crystallization of Microbial Rhodopsins

Microbial rhodopsins are light-sensitive transmembrane proteins, evolutionary adapted by various organisms like archaea, bacteria, simple eukaryote, and viruses to utilize solar energy for their survival. A complete understanding of functional mechanisms of these proteins is not possible without the knowledge of their high-resolution structures, which can be primarily obtained by X-ray crystallography. This technique, however, requires high-quality crystals, growing of which is a great challenge especially in case of membrane proteins. In this chapter, we summarize methods applied for crystallization of microbial rhodopsins with the emphasis on crystallization in lipidic mesophases, also known as in meso approach. In particular, we describe in detail the methods of crystallization using lipidic cubic phase to grow both large crystals optimized for traditional crystallographic data collection and microcrystals for serial crystallography.

Structures of the archaerhodopsin-3 transporter reveal that disordering of internal water networks underpins receptor sensitization

Many transmembrane receptors have a desensitized state, in which they are unable to respond to external stimuli. The family of microbial rhodopsin proteins includes one such group of receptors, whose inactive or dark-adapted (DA) state is established in the prolonged absence of light. Here, we present high-resolution crystal structures of the ground (light-adapted) and DA states of Archaerhodopsin-3 (AR3), solved to 1.1 Å and 1.3 Å resolution respectively. We observe significant differences between the two states in the dynamics of water molecules that are coupled via H-bonds to the retinal Schiff Base. Supporting QM/MM calculations reveal how the DA state permits a thermodynamic equilibrium between retinal isomers to be established, and how this same change is prevented in the ground state in the absence of light. We suggest that the different arrangement of internal water networks in AR3 is responsible for the faster photocycle kinetics compared to homologs.

NeoR, a near-infrared absorbing rhodopsin

The Rhizoclosmatium globosum genome encodes three rhodopsin-guanylyl cyclases (RGCs), which are predicted to facilitate visual orientation of the fungal zoospores. Here, we show that RGC1 and RGC2 function as light-activated cyclases only upon heterodimerization with RGC3 (NeoR). RGC1/2 utilize conventional green or blue-light-sensitive rhodopsins (λmax = 550 and 480 nm, respectively), with short-lived signaling states, responsible for light-activation of the enzyme. The bistable NeoR is photoswitchable between a near-infrared-sensitive (NIR, λmax = 690 nm) highly fluorescent state (QF = 0.2) and a UV-sensitive non-fluorescent state, thereby modulating the activity by NIR pre-illumination. No other rhodopsin has been reported so far to be functional as a heterooligomer, or as having such a long wavelength absorption or high fluorescence yield. Site-specific mutagenesis and hybrid quantum mechanics/molecular mechanics simulations support the idea that the unusual photochemical properties result from the rigidity of the retinal chromophore and a unique counterion triad composed of two glutamic and one aspartic acids. These findings substantially expand our understanding of the natural potential and limitations of spectral tuning in rhodopsin photoreceptors.

a-ARM: Automatic Rhodopsin Modeling with Chromophore Cavity Generation, Ionization State Selection, and External Counterion Placement

The Automatic Rhodopsin Modeling (ARM) protocol has recently been proposed as a tool for the fast and parallel generation of basic hybrid quantum mechanics/molecular mechanics (QM/MM) models of wild type and mutant rhodopsins. However, in its present version, input preparation requires a few hours long user's manipulation of the template protein structure, which also impairs the reproducibility of the generated models. This limitation, which makes model building semiautomatic rather than fully automatic, comprises four tasks: definition of the retinal chromophore cavity, assignment of protonation states of the ionizable residues, neutralization of the protein with external counterions, and finally congruous generation of single or multiple mutations. In this work, we show that the automation of the original ARM protocol can be extended to a level suitable for performing the above tasks without user's manipulation and with an input preparation time of minutes. The new protocol, called a-ARM, delivers fully reproducible (i.e., user independent) rhodopsin QM/MM models as well as an improved model quality. More specifically, we show that the trend in vertical excitation energies observed for a set of 25 wild type and 14 mutant rhodopsins is predicted by the new protocol better than when using the original. Such an agreement is reflected by an estimated (relative to the probed set) trend deviation of 0.7 ± 0.5 kcal mol-1 (0.03 ± 0.02 eV) and mean absolute error of 1.0 kcal mol-1 (0.04 eV).

Redshifted and Near-infrared Active Analog Pigments Based upon Archaerhodopsin-3

Archaerhodopsin‐3 (AR3) is a member of the microbial rhodopsin family of hepta‐helical transmembrane proteins, containing a covalently bound molecule of all‐trans retinal as a chromophore. It displays an absorbance band in the visible region of the solar spectrum (λmax 556 nm) and functions as a light‐driven proton pump in the archaeon Halorubrum sodomense. AR3 and its mutants are widely used in neuroscience as optogenetic neural silencers and in particular as fluorescent indicators of transmembrane potential. In this study, we investigated the effect of analogs of the native ligand all‐trans retinal A1 on the spectral properties and proton‐pumping activity of AR3 and its single mutant AR3 (F229S). While, surprisingly, the 3‐methoxyretinal A2 analog did not redshift the absorbance maximum of AR3, the analogs retinal A2 and 3‐methylamino‐16‐nor‐1,2,3,4‐didehydroretinal (MMAR) did generate active redshifted AR3 pigments. The MMAR analog pigments could even be activated by near‐infrared light. Furthermore, the MMAR pigments showed strongly enhanced fluorescence with an emission band in the near‐infrared peaking around 815 nm. We anticipate that the AR3 pigments generated in this study have widespread potential for near‐infrared exploitation as fluorescent voltage‐gated sensors in optogenetics and artificial leafs and as proton pumps in bioenergy‐based applications.

Multiscale Simulations of Biological Membranes: The Challenge To Understand Biological Phenomena in a Living Substance

Biological membranes are tricky to investigate. They are complex in terms of molecular composition and structure, functional over a wide range of time scales, and characterized by nonequilibrium conditions. Because of all of these features, simulations are a great technique to study biomembrane behavior. A significant part of the functional processes in biological membranes takes place at the molecular level; thus computer simulations are the method of choice to explore how their properties emerge from specific molecular features and how the interplay among the numerous molecules gives rise to function over spatial and time scales larger than the molecular ones. In this review, we focus on this broad theme. We discuss the current state-of-the-art of biomembrane simulations that, until now, have largely focused on a rather narrow picture of the complexity of the membranes. Given this, we also discuss the challenges that we should unravel in the foreseeable future. Numerous features such as the actin-cytoskeleton network, the glycocalyx network, and nonequilibrium transport under ATP-driven conditions have so far received very little attention; however, the potential of simulations to solve them would be exceptionally high. A major milestone for this research would be that one day we could say that computer simulations genuinely research biological membranes, not just lipid bilayers.

X-ray structure analysis of bacteriorhodopsin at 1.3 Å resolution

Bacteriorhodopsin (bR) of Halobacterium salinarum is a membrane protein that acts as a light-driven proton pump. bR and its homologues have recently been utilized in optogenetics and other applications. Although the structures of those have been reported so far, the resolutions are not sufficient for elucidation of the intrinsic structural features critical to the color tuning and ion pumping properties. Here we report the accurate crystallographic analysis of bR in the ground state. The influence of X-rays was suppressed by collecting the data under a low irradiation dose at 15 K. Consequently, individual atoms could be separately observed in the electron density map at better than 1.3 Å resolution. Residues from Thr5 to Ala233 were continuously constructed in the model. The twist of the retinal polyene was determined to be different from those in the previous models. Two conformations were observed for the proton release region. We discuss the meaning of these fine structural features.

Functional roles of tyrosine 185 during the bacteriorhodopsin photocycle as revealed by in situ spectroscopic studies

Tyrosine 185 (Y185), one of the aromatic residues within the retinal (Ret) chromophore binding pocket in helix F of bacteriorhodopsin (bR), is highly conserved among the microbial rhodopsin family proteins. Many studies have investigated the functions of Y185, but its underlying mechanism during the bR photocycle remains unclear. To address this research gap, in situ two-dimensional (2D) magic-angle spinning (MAS) solid-state NMR (ssNMR) of specifically labelled bR, combined with light-induced transient absorption change measurements, dynamic light scattering (DLS) measurements, titration analysis and site-directed mutagenesis, was used to elucidate the functional roles of Y185 during the bR photocycle in the native membrane environment. Different interaction modes were identified between Y185 and the Ret chromophore in the dark-adapted (inactive) state and M (active) state, indicating that Y185 may serve as a rotamer switch maintaining the protein dynamics, and plays an important role in the efficient proton-pumping mechanism in the bR purple membrane.

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.

Directed Evolution of a Bright Near-Infrared Fluorescent Rhodopsin Using a Synthetic Chromophore

By engineering a microbial rhodopsin, Archaerhodopsin-3 (Arch), to bind a synthetic chromophore, merocyanine retinal, in place of the natural chromophore all-trans-retinal (ATR), we generated a protein with exceptionally bright and unprecedentedly red-shifted near-infrared (NIR) fluorescence. We show that chromophore substitution generates a fluorescent Arch complex with a 200-nm bathochromic excitation shift relative to ATR-bound wild-type Arch and an emission maximum at 772 nm. Directed evolution of this complex produced variants with pH-sensitive NIR fluorescence and molecular brightness 8.5-fold greater than the brightest ATR-bound Arch variant. The resulting proteins are well suited to bacterial imaging; expression and stability have not been optimized for mammalian cell imaging. By targeting both the protein and its chromophore, we overcome inherent challenges associated with engineering bright NIR fluorescence into Archaerhodopsin. This work demonstrates an efficient strategy for engineering non-natural, tailored properties into microbial opsins, properties relevant for imaging and interrogating biological systems.

Structural basis for the slow photocycle and late proton release in Acetabularia rhodopsin I from the marine plant Acetabularia acetabulum

Although many crystal structures of microbial rhodopsins have been solved, those with sufficient resolution to identify the functional water molecules are very limited. In this study, the Acetabularia rhodopsin I (ARI) protein derived from the marine alga A. acetabulum was synthesized on a large scale by the Escherichia coli cell-free membrane-protein production method, and crystal structures of ARI were determined at the second highest (1.52–1.80 Å) resolution for a microbial rhodopsin, following bacteriorhodopsin (BR). Examinations of the photochemical properties of ARI revealed that the photocycle of ARI is slower than that of BR and that its proton-transfer reactions are different from those of BR. In the present structures, a large cavity containing numerous water molecules exists on the extracellular side of ARI, explaining the relatively low pK a of Glu206ARI, which cannot function as an initial proton-releasing residue at any pH. An interhelical hydrogen bond exists between Leu97ARI and Tyr221ARI on the cytoplasmic side, which facilitates the slow photocycle and regulates the pK a of Asp100ARI, a potential proton donor to the Schiff base, in the dark state.

Disordered regions in transmembrane proteins

The functions of transmembrane proteins in living cells are widespread; they range from various transport processes to energy production, from cell-cell adhesion to communication. Structurally, they are highly ordered in their membrane-spanning regions, but may contain disordered regions in the cytosolic and extra-cytosolic parts. In this study, we have investigated the disordered regions in transmembrane proteins by a stringent definition of disordered residues on the currently available largest experimental dataset, and show a significant correlation between the spatial distributions of positively charged residues and disordered regions. This finding suggests a new role of disordered regions in transmembrane proteins by providing structural flexibility for stabilizing interactions with negatively charged head groups of the lipid molecules. We also find a preference of structural disorder in the terminal--as opposed to loop--regions in transmembrane proteins, and survey the respective functions involved in recruiting other proteins or mediating allosteric signaling effects. Finally, we critically compare disorder prediction methods on our transmembrane protein set. While there are no major differences between these methods using the usual statistics, such as per residue accuracies, Matthew's correlation coefficients, etc.; substantial differences can be found regarding the spatial distribution of the predicted disordered regions. We conclude that a predictor optimized for transmembrane proteins would be of high value to the field of structural disorder.