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

Bacteriorhodopsin proton-pumping mechanism: Successes and challenges in computational approaches

Bacteriorhodopsin (bR) is perhaps the best-studied proton pump. Over about four decades, research on this fascinating photocyclic light-driven protein inspired the development of key experimental and computational methodologies that are now widely used in membrane protein studies. We review here failures and successes in computational approaches that have been applied to study the bR proton-transfer steps. Conflict between experimental results pertaining to the proton transfer mechanisms in the early photocycle intermediates was resolved by detailed quantum mechanical/molecular mechanical computation, the results of which were confirmed more than a decade later. Key to this approach was the realization that, to understand how the pump works and achieves directional transfer of protons, the individual reaction steps-proton transfer and reorganization of the internal hydrogen-bond network-needed to be considered within the context of the energy landscape of the complete reaction cycle.

Hydrogen-Bonding and Hydrophobic Interaction Networks as Structural Determinants of Microbial Rhodopsin Function

Microbial pump rhodopsins are highly versatile light-driven membrane proteins that couple protein conformational dynamics with ion translocation across the cell membranes. Understanding how microbial pump rhodopsins use specific amino acid residues at key functional sites to control ion selectivity and ion pumping direction is of general interest for membrane transporters, and could guide site-directed mutagenesis for optogenetics applications. To enable direct comparisons between proteins with different sequences we implement, for the first time, a unique numbering scheme for the microbial pump rhodopsin residues, NS-mrho. We use NS-mrho to show that distinct microbial pump rhodopsins typically have hydrogen-bond networks that are less conserved than anticipated from the amino acid residue conservation, whereas their hydrophobic interaction networks are largely conserved. To illustrate the role of the hydrogen-bond networks as structural elements that determine the functionality of microbial pump rhodopsins, we performed experiments, atomic-level simulations, and hydrogen bond network analyses on GR, the outward proton pump from Gloeobacter violaceus, and KR2, the outward sodium pump from Krokinobacter eikastus. The experiments indicate that multiple mutations that recover KR2 amino acid residues in GR not only fail to convert it into a sodium pump, but completely inactivate GR by abolishing photoisomerization of the retinal chromophore. This observation could be attributed to the drastically altered hydrogen-bond interaction network identified with simulations and network analyses. Taken together, our findings suggest that functional specificity could be encoded in the collective hydrogen-bond network of microbial pump rhodopsins.

Appraising protein conformational changes by resampling time-resolved serial x-ray crystallography data

With the development of serial crystallography at both x-ray free electron laser and synchrotron radiation sources, time-resolved x-ray crystallography is increasingly being applied to study conformational changes in macromolecules. A successful time-resolved serial crystallography study requires the growth of microcrystals, a mechanism for synchronized and homogeneous excitation of the reaction of interest within microcrystals, and tools for structural interpretation. Here, we utilize time-resolved serial femtosecond crystallography data collected from microcrystals of bacteriorhodopsin to compare results from partial occupancy structural refinement and refinement against extrapolated data. We illustrate the domain wherein the amplitude of refined conformational changes is inversely proportional to the activated state occupancy. We illustrate how resampling strategies allow coordinate uncertainty to be estimated and demonstrate that these two approaches to structural refinement agree within coordinate errors. We illustrate how singular value decomposition of a set of difference Fourier electron density maps calculated from resampled data can minimize phase bias in these maps, and we quantify residual densities for transient water molecules by analyzing difference Fourier and Polder omit maps from resampled data. We suggest that these tools may assist others in judging the confidence with which observed electron density differences may be interpreted as functionally important conformational changes.

Water, Protons, and the Gating of Voltage-Gated Potassium Channels

Ion channels are ubiquitous throughout all forms of life. Potassium channels are even found in viruses. Every cell must communicate with its surroundings, so all cells have them, and excitable cells, in particular, especially nerve cells, depend on the behavior of these channels. Every channel must be open at the appropriate time, and only then, so that each channel opens in response to the stimulus that tells that channel to open. One set of channels, including those in nerve cells, responds to voltage. There is a standard model for the gating of these channels that has a section of the protein moving in response to the voltage. However, there is evidence that protons are moving, rather than protein. Water is critical as part of the gating process, although it is hard to see how this works in the standard model. Here, we review the extensive evidence of the importance of the role of water and protons in gating these channels. Our principal example, but by no means the only example, will be the Kv1.2 channel. Evidence comes from the effects of D2O, from mutations in the voltage sensing domain, as well as in the linker between that domain and the gate, and at the gate itself. There is additional evidence from computations, especially quantum calculations. Structural evidence comes from X-ray studies. The hydration of ions is critical in the transfer of ions in constricted spaces, such as the gate region and the pore of a channel; we will see how the structure of the hydrated ion fits with the structure of the channel. In addition, there is macroscopic evidence from osmotic experiments and streaming current measurements. The combined evidence is discussed in the context of a model that emphasizes the role of protons and water in gating these channels.

Protons taken hostage: Dynamic H-bond networks of the pH-sensing GPR68

Proton-sensing G Protein Coupled Receptors (GPCRs) sense changes in the extracellular pH to effect cell signaling for cellular homeostasis. They tend to be overexpressed in solid tumors associated with acidic extracellular pH, and are of direct interest as drug targets. How proton-sensing GPCRs sense extracellular acidification and activate upon protonation change is important to understand, because it may guide the design of therapeutics. Lack of publicly available experimental structures make it challenging to discriminate between conflicting mechanisms proposed for proton-binding, as main roles have been assigned to either an extracellular histidine cluster or to an internal carboxylic triad. Here we present a protocol to derive and evaluate structural models of the proton-sensing GPR68. This approach integrates state-of-the-art homology modeling with microsecond-timescale atomistic simulations, and with a detailed assessment of the compatibility of the structural models with known structural features of class A GPCRs. To decipher structural elements of potential interest for protonation-coupled conformational changes of GPR68, we used the best-compatible model as a starting point for independent atomistic simulations of GPR68 with different protonation states, and graph computations to characterize the response of GPR68 to changes in protonation. We found that GPR68 hosts an extended hydrogen-bond network that inter-connects the extracellular histidine cluster to the internal carboxylic triad, and which can even reach groups at the cytoplasmic G-protein binding site. Taken together, results suggest that GPR68 relies on dynamic, hydrogen-bond networks to inter-connect extracellular and internal proton-binding sites, and to elicit conformational changes at the cytoplasmic G-protein binding site.

Detailed analysis of distorted retinal and its interaction with surrounding residues in the K intermediate of bacteriorhodopsin

The K intermediate of proton pumping bacteriorhodopsin is the first intermediate generated after isomerization of retinal to the 13-cis form. Although various structures have been reported for the K intermediate until now, these differ from each other, especially in terms of the conformation of the retinal chromophore and its interaction with surrounding residues. We report here an accurate X-ray crystallographic analysis of the K structure. The polyene chain of 13-cis retinal is observed to be S-shaped. The side chain of Lys216, which is covalently bound to retinal via the Schiff-base linkage, interacts with residues, Asp85 and Thr89. In addition, the Nζ-H of the protonated Schiff-base linkage interacts with a residue, Asp212 and a water molecule, W402. Based on quantum chemical calculations for this K structure, we examine the stabilizing factors of distorted conformation of retinal and propose a relaxation manner to the next L intermediate.

Mechanism of Absorption Wavelength Shift of Bacteriorhodopsin During Photocycle

Bacteriorhodopsin, a light-driven proton pump, alters the absorption wavelengths in the range of 410-617 nm during the photocycle. Here, we report the absorption wavelengths, calculated using 12 bacteriorhodopsin crystal structures (including the BR, BR_13-cis, J, K₀, K_E, K_L, L, M, N, and O state structures) and a combined quantum mechanical/molecular mechanical/polarizable continuum model (QM/MM/PCM) approach. The QM/MM/PCM calculations reproduced the experimentally measured absorption wavelengths with a standard deviation of 4 nm. The shifts in the absorption wavelengths can be explained mainly by the following four factors: (i) retinal Schiff base deformation/twist induced by the protein environment, leading to a decrease in the electrostatic interaction between the protein environment and the retinal Schiff base; (ii) changes in the protonation state of the protein environment, directly altering the electrostatic interaction between the protein environment and the retinal Schiff base; (iii) changes in the protonation state; or (iv) isomerization of the retinal Schiff base, where the absorption wavelengths of the isomers originally differ.

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.

True-atomic-resolution insights into the structure and functional role of linear chains and low-barrier hydrogen bonds in proteins

Hydrogen bonds are fundamental to the structure and function of biological macromolecules and have been explored in detail. The chains of hydrogen bonds (CHBs) and low-barrier hydrogen bonds (LBHBs) were proposed to play essential roles in enzyme catalysis and proton transport. However, high-resolution structural data from CHBs and LBHBs is limited. The challenge is that their 'visualization' requires ultrahigh-resolution structures of the ground and functionally important intermediate states to identify proton translocation events and perform their structural assignment. Our true-atomic-resolution structures of the light-driven proton pump bacteriorhodopsin, a model in studies of proton transport, show that CHBs and LBHBs not only serve as proton pathways, but also are indispensable for long-range communications, signaling and proton storage in proteins. The complete picture of CHBs and LBHBs discloses their multifunctional roles in providing protein functions and presents a consistent picture of proton transport and storage resolving long-standing debates and controversies.

Microbial Rhodopsins

The first microbial rhodopsin, a light-driven proton pump bacteriorhodopsin from Halobacterium salinarum (HsBR), was discovered in 1971. Since then, this seven-α-helical protein, comprising a retinal molecule as a cofactor, became a major driver of groundbreaking developments in membrane protein research. However, until 1999 only a few archaeal rhodopsins, acting as light-driven proton and chloride pumps and also photosensors, were known. A new microbial rhodopsin era started in 2000 when the first bacterial rhodopsin, a proton pump, was discovered. Later it became clear that there are unexpectedly many rhodopsins, and they are present in all the domains of life and even in viruses. It turned out that they execute such a diversity of functions while being "nearly the same." The incredible evolution of the research area of rhodopsins and the scientific and technological potential of the proteins is described in the review with a focus on their function-structure relationships.

Quantum-Mechanical Molecular Dynamics Simulations on Secondary Proton Transfer in Bacteriorhodopsin Using Realistic Models

Bacteriorhodopsin (BR) transports a proton from intracellular to extracellular (EC) sites through five proton transfers. The second proton transfer is the release of an excess proton stored in BR into the EC medium, and an atomistic understanding of this whole process has remained unexplored due to its ubiquitous environment. Here, fully quantum mechanical (QM) molecular dynamics (MD) and metadynamics (MTD) simulations for this process were performed at the divide-and-conquer density-functional tight-binding level using realistic models (∼50000 and ∼20000 atoms) based on the time-resolved photointermediate structures from an X-ray free electron laser. Regarding the proton storage process, the QM-MD/MTD simulations confirmed the Glu-shared mechanism, in which an excess proton is stored between Glu194 and Glu204, and clarified that the activation occurs by localizing the proton at Glu204 in the photocycle. Furthermore, the QM-MD/MTD simulations elucidated a release pathway from Glu204 through Ser193 to the EC water molecules and clarified that the proton release starts at ∼250 μs. In the ubiquitous proton diffusion in the EC medium, the transient proton receptors predicted experimentally were assigned to carboxylates in Glu9 and Glu74. Large-scale QM-MD/MTD simulations beyond the conventional sizes, which provided the above findings and confirmations, were possible by adopting our Dcdftbmd program.

Proton-Binding Motifs of Membrane-Bound Proteins: From Bacteriorhodopsin to Spike Protein S

Membrane-bound proteins that change protonation during function use specific protein groups to bind and transfer protons. Knowledge of the identity of the proton-binding groups is of paramount importance to decipher the reaction mechanism of the protein, and protonation states of prominent are studied extensively using experimental and computational approaches. Analyses of model transporters and receptors from different organisms, and with widely different biological functions, indicate common structure-sequence motifs at internal proton-binding sites. Proton-binding dynamic hydrogen-bond networks that are exposed to the bulk might provide alternative proton-binding sites and proton-binding pathways. In this perspective article I discuss protonation coupling and proton binding at internal and external carboxylate sites of proteins that use proton transfer for function. An inter-helical carboxylate-hydroxyl hydrogen-bond motif is present at functionally important sites of membrane proteins from archaea to the brain. External carboxylate-containing H-bond clusters are observed at putative proton-binding sites of protonation-coupled model proteins, raising the question of similar functionality in spike protein S.

Frontiers in Multiscale Modeling of Photoreceptor Proteins

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.

MerMAIDs: a family of metagenomically discovered marine anion-conducting and intensely desensitizing channelrhodopsins

Channelrhodopsins (ChRs) are algal light-gated ion channels widely used as optogenetic tools for manipulating neuronal activity. ChRs desensitize under continuous bright-light illumination, resulting in a significant decline of photocurrents. Here we describe a metagenomically identified family of phylogenetically distinct anion-conducting ChRs (designated MerMAIDs). MerMAIDs almost completely desensitize during continuous illumination due to accumulation of a late non-conducting photointermediate that disrupts the ion permeation pathway. MerMAID desensitization can be fully explained by a single photocycle in which a long-lived desensitized state follows the short-lived conducting state. A conserved cysteine is the critical factor in desensitization, as its mutation results in recovery of large stationary photocurrents. The rapid desensitization of MerMAIDs enables their use as optogenetic silencers for transient suppression of individual action potentials without affecting subsequent spiking during continuous illumination. Our results could facilitate the development of optogenetic tools from metagenomic databases and enhance general understanding of ChR function.

Bacteriorhodopsin: Structural Insights Revealed Using X-Ray Lasers and Synchrotron Radiation

Directional transport of protons across an energy transducing membrane—proton pumping—is ubiquitous in biology. Bacteriorhodopsin (bR) is a light-driven proton pump that is activated by a buried all- trans retinal chromophore being photoisomerized to a 13- cis conformation. The mechanism by which photoisomerization initiates directional proton transport against a proton concentration gradient has been studied by a myriad of biochemical, biophysical, and structural techniques. X-ray free electron lasers (XFELs) have created new opportunities to probe the structural dynamics of bR at room temperature on timescales from femtoseconds to milliseconds using time-resolved serial femtosecond crystallography (TR-SFX). Wereview these recent developments and highlight where XFEL studies reveal new details concerning the structural mechanism of retinal photoisomerization and proton pumping. We also discuss the extent to which these insights were anticipated by earlier intermediate trapping studies using synchrotron radiation. TR-SFX will open up the field for dynamical studies of other proteins that are not naturally light-sensitive.

Reactions at Biomembrane Interfaces

Membranes surrounding the biological cell and its internal compartments host proteins that catalyze chemical reactions essential for the functioning of the cell. Rather than being a passive structural matrix that holds membrane-embedded proteins in place, the membrane can largely shape the conformational energy landscape of membrane proteins and impact the energetics of their chemical reaction. Here, we highlight the challenges in understanding how lipids impact the conformational energy landscape of macromolecular membrane complexes whose functioning involves chemical reactions including proton transfer. We review here advances in our understanding of how chemical reactions occur at membrane interfaces gleaned with both theoretical and experimental advances using simple protein systems as guides. Our perspective is that of bridging experiments with theory to understand general physicochemical principles of membrane reactions, with a long term goal of furthering our understanding of the role of the lipids on the functioning of complex macromolecular assemblies at the membrane interface.

Interhelical interactions between D92 and C218 in the cytoplasmic domain regulate proton uptake upon N-decay in the proton transport of Acetabularia rhodopsin II

Acetabularia rhodopsin II (ARII or Ace2), an outward light-driven algal proton pump found in the giant unicellular marine alga Acetabularia acetabulum, has a unique property in the cytoplasmic (CP) side of its channel. The X-ray crystal structure of ARII in a dark state suggested the formation of an interhelical hydrogen bond between C218^ARII and D92^ARII, an internal proton donor to the Schiff base (Wada et al., 2011). In this report, we investigated the photocycles of two mutants at position C218^ARII: C218A^ARII which disrupts the interaction with D92^ARII, and C218S^ARII which potentially forms a stronger hydrogen bond. Both mutants exhibited slower photocycles compared to the wild-type pump. Together with several kinetic changes of the photoproducts in the first half of the photocycle, these replacements led to specific retardation of the N-to-O transition in the second half of the photocycle. In addition, measurements of the flash-induced proton uptake and release using a pH-sensitive indium-tin oxide electrode revealed a concomitant delay in the proton uptake. These observations strongly suggest the importance of a native weak hydrogen bond between C218^ARII and D92^ARII for proper proton translocation in the CP channel during N-decay. A putative role for the D92^ARII-C218^ARII interhelical hydrogen bond in the function of ARII is discussed.

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.

The photochemistry of sodium ion pump rhodopsin observed by watermarked femto- to submillisecond stimulated Raman spectroscopy

Krokinobacter rhodopsin 2 (KR2) is a recently discovered light-driven Na(+) pump that holds significant promise for application as a neural silencer in optogenetics. KR2 transports Na(+) (in NaCl solution) or H(+) (in larger cation solution, e.g. in CsCl) during its photocycle. Here, we investigate the photochemistry of KR2 with the recently developed watermarked, baseline-free femto- to submillisecond transient stimulated Raman spectroscopy (TSRS), which enables us to investigate retinal chromophore dynamics in real time with high spectral resolution over a large time range. We propose a new photocycle from femtoseconds to submilliseconds: J (formed in ∼200 fs) → K (∼3 ps) → K/L1 (∼20 ps) → K/L2 (∼30 ns) → L/M (∼20 μs). KR2 binds a Na(+) ion that is not transported on the extracellular side, of which the function is unclear. We demonstrate with TSRS that for the D102N mutant in NaCl (with Na(+) unbound, Na(+) transport) and for WT KR2 in CsCl (with Na(+) unbound, H(+) transport), the extracellular Na(+) binding significantly influences the intermediate K/L/M state equilibrium on the photocycle, while the identity of the transported ion, Na(+) or H(+), does not affect the photocycle. Our findings will contribute to further elucidation of the molecular mechanisms of KR2.

A delocalized proton-binding site within a membrane protein

The role of protein-bound water molecules in protein function and catalysis is an emerging topic. Here, we studied the solvation of an excess proton by protein-bound water molecules and the contribution of the surrounding amino acid residues at the proton release site of the membrane protein bacteriorhodopsin. It hosts an excess proton within a protein-bound water cluster, which is hydrogen bonded to several surrounding amino acids. Indicative of delocalization is a broad continuum absorbance experimentally observed by time-resolved Fourier transform infrared spectroscopy. In combination with site-directed mutagenesis, the involvement of several amino acids (especially Glu-194 and Glu-204) in the delocalization was elaborated. Details regarding the contributions of the glutamates and water molecules to the delocalization mode in biomolecular simulations are controversial. We carried out quantum mechanics/molecular mechanics (QM/MM) self-consistent charge density functional tight-binding simulations for all amino acids that have been experimentally shown to be involved in solvation of the excess proton, and systematically investigated the influence of the quantum box size. We compared calculated theoretical infrared spectra with experimental ones as a measure for the correct description of excess proton delocalization. A continuum absorbance can only be observed for small quantum boxes containing few amino acids and/or water molecules. Larger quantum boxes, including all experimentally shown involved amino acids, resulted in narrow absorbance bands, indicating protonation of a single binding site in contradiction to experimental results. We conclude that small quantum boxes seem to reproduce representative extreme cases of proton delocalization modes: proton delocalization only on water molecules or only between Glu-194 and Glu-204. Extending the experimental spectral region to lower wave numbers, a water-delocalized proton reproduces the observed continuum absorbance better than a glutamate-shared delocalized proton. However, a full agreement between QM simulations and experimental results on the delocalized excess proton will require a larger quantum box as well as more sophisticated QM/MM methods.

Crystal structure of a light-driven sodium pump

Recently, the first known light-driven sodium pumps, from the microbial rhodopsin family, were discovered. We have solved the structure of one of them, Krokinobacter eikastus rhodopsin 2 (KR2), in the monomeric blue state and in two pentameric red states, at resolutions of 1.45 Å and 2.2 and 2.8 Å, respectively. The structures reveal the ion-translocation pathway and show that the sodium ion is bound outside the protein at the oligomerization interface, that the ion-release cavity is capped by a unique N-terminal α-helix and that the ion-uptake cavity is unexpectedly large and open to the surface. Obstruction of the cavity with the mutation G263F imparts KR2 with the ability to pump potassium. These results pave the way for the understanding and rational design of cation pumps with new specific properties valuable for optogenetics.

Structure of archaerhodopsin-2 at 1.8 Å resolution

Archaerhodopsin-2 (aR2), the sole protein found in the claret membrane of Halorubrum sp. Aus-2, functions as a light-driven proton pump. In this study, structural analysis of aR2 was performed using a novel three-dimensional crystal prepared by the successive fusion of claret membranes. The crystal is made up of stacked membranes, in each of which aR2 trimers are arranged on a hexagonal lattice. This lattice structure resembles that found in the purple membrane of H. salinarum, except that lipid molecules trapped within the trimeric structure are not distributed with perfect threefold symmetry. Nonetheless, diffraction data at 1.8 Å resolution provide accurate structural information about functionally important residues. It is shown that two glutamates in the proton-release channel form a paired structure that is maintained by a low-barrier hydrogen bond. Although the structure of the proton-release pathway is highly conserved among proton-pumping archaeal rhodopsins, aR2 possesses the following peculiar structural features: (i) the motional freedom of the tryptophan residue that makes contact with the C13 methyl group of retinal is restricted, affecting the formation/decay kinetics of the L state, and (ii) the N-terminal polypeptide folds into an Ω-loop, which may play a role in organizing the higher-order structure.

Crystal structure of Cruxrhodopsin-3 from Haloarcula vallismortis

Cruxrhodopsin-3 (cR3), a retinylidene protein found in the claret membrane of Haloarcula vallismortis, functions as a light-driven proton pump. In this study, the membrane fusion method was applied to crystallize cR3 into a crystal belonging to space group P321. Diffraction data at 2.1 Å resolution show that cR3 forms a trimeric assembly with bacterioruberin bound to the crevice between neighboring subunits. Although the structure of the proton-release pathway is conserved among proton-pumping archaeal rhodopsins, cR3 possesses the following peculiar structural features: 1) The DE loop is long enough to interact with a neighboring subunit, strengthening the trimeric assembly; 2) Three positive charges are distributed at the cytoplasmic end of helix F, affecting the higher order structure of cR3; 3) The cytoplasmic vicinity of retinal is more rigid in cR3 than in bacteriorhodopsin, affecting the early reaction step in the proton-pumping cycle; 4) the cytoplasmic part of helix E is greatly bent, influencing the proton uptake process. Meanwhile, it was observed that the photobleaching of retinal, which scarcely occurred in the membrane state, became significant when the trimeric assembly of cR3 was dissociated into monomers in the presence of an excess amount of detergent. On the basis of these observations, we discuss structural factors affecting the photostabilities of ion-pumping rhodopsins.

Bacteriorhodopsin: Would the real structural intermediates please stand up?

BACKGROUND

Bacteriorhodopsin (bR) is the simplest known light driven proton pump and has been heavily studied using structural methods: eighty four X-ray diffraction, six electron diffraction and three NMR structures of bR are deposited within the protein data bank. Twenty one X-ray structures report light induced structural changes and changes induced by mutation, changes in pH, thermal annealing or X-ray induced photo-reduction have also been examined.

SCOPE OF REVIEW

We argue that light-induced structural changes that are replicated across several studies by independent research groups are those most likely to represent what is happening in reality. We present both internal distance matrix analyses that sort deposited bR structures into hierarchal trees, and difference Fourier analysis of deposited X-ray diffraction data.

MAJOR CONCLUSIONS

An internal distance matrix analysis separates most wild-type bR structures according to their different crystal forms, indicating how the protein's structure is influenced by crystallization conditions. A similar analysis clusters eleven studies of illuminated bR crystals as one branch of a hierarchal tree with reproducible movements of the extracellular portion of helix C towards helix G, and of the cytoplasmic portion of helix F away from helices A, B and G. All crystallographic data deposited for illuminated crystals show negative difference density on a water molecule (Wat402) that forms H-bonds to the retinal Schiff Base and two aspartate residues (Asp85, Asp212) in the bR resting state. Other recurring difference density features indicated reproducible side-chain, backbone and water molecule displacements. X-ray induced radiation damage also disorders Wat402 but acts via cleaving the head-groups of Asp85 and Asp212.

GENERAL SIGNIFICANCE

A remarkable level of agreement exists when deposited structures and crystallographic observations are viewed as a whole. From this agreement a unified picture of the structural mechanism of light-induced proton pumping by bR emerges. This article is part of a Special Issue entitled Structural biochemistry and biophysics of membrane proteins.

Channelrhodopsins: a bioinformatics perspective

Channelrhodopsins are microbial-type rhodopsins that function as light-gated cation channels. Understanding how the detailed architecture of the protein governs its dynamics and specificity for ions is important, because it has the potential to assist in designing site-directed channelrhodopsin mutants for specific neurobiology applications. Here we use bioinformatics methods to derive accurate alignments of channelrhodopsin sequences, assess the sequence conservation patterns and find conserved motifs in channelrhodopsins, and use homology modeling to construct three-dimensional structural models of channelrhodopsins. The analyses reveal that helices C and D of channelrhodopsins contain Cys, Ser, and Thr groups that can engage in both intra- and inter-helical hydrogen bonds. We propose that these polar groups participate in inter-helical hydrogen-bonding clusters important for the protein conformational dynamics and for the local water interactions. This article is part of a Special Issue entitled: Retinal Proteins - You can teach an old dog new tricks.

Exploring cavity dynamics in biomolecular systems

Background

The internal cavities of proteins are dynamic structures and their dynamics may be associated with conformational changes which are required for the functioning of the protein. In order to study the dynamics of these internal protein cavities, appropriate tools are required that allow rapid identification of the cavities as well as assessment of their time-dependent structures.

Results

In this paper, we present such a tool and give results that illustrate the applicability for the analysis of molecular dynamics trajectories. Our algorithm consists of a pre-processing step where the structure of the cavity is computed from the Voronoi diagram of the van der Waals spheres based on coordinate sets from the molecular dynamics trajectory. The pre-processing step is followed by an interactive stage, where the user can compute, select and visualize the dynamic cavities. Importantly, the tool we discuss here allows the user to analyze the time-dependent changes of the components of the cavity structure. An overview of the cavity dynamics is derived by rendering the dynamic cavities in a single image that gives the cavity surface colored according to its time-dependent dynamics.

Conclusion

The Voronoi-based approach used here enables the user to perform accurate computations of the geometry of the internal cavities in biomolecules. For the first time, it is possible to compute dynamic molecular paths that have a user-defined minimum constriction size. To illustrate the usefulness of the tool for understanding protein dynamics, we probe the dynamic structure of internal cavities in the bacteriorhodopsin proton pump.

Organization and dynamics of the N-terminal domain of chemokine receptor CXCR1 in reverse micelles: effect of graded hydration

Water plays a fundamental role in the folding, structure, dynamics, and function of proteins and peptides. The extracellular N-terminal domain of chemokine receptors is crucial in mediating binding affinity, receptor selectivity, and regulating function. The flexible N-terminal domain becomes ordered in membranes and membrane-mimetic assemblies, thereby indicating that the membrane could play an important role in regulating CXC chemokine receptor 1 (CXCR1) function. In view of the role of hydration in lipid-protein interactions in membranes, we explored the organization and dynamics of a 34-mer peptide of the CXCR1 N-terminal domain in reverse micelles by utilizing a combination of fluorescence-based approaches and circular dichroism spectroscopy. Our results show that the secondary structure adopted by the CXCR1 N-domain is critically dependent on hydration. The tryptophan residues of the CXCR1 N-domain experience motional restriction and exhibit red edge excitation shift (REES) upon incorporation in reverse micelles. REES and fluorescence lifetime exhibit reduction with increasing reverse micellar hydration. Time-resolved fluorescence anisotropy measurements reveal the effect of hydration on peptide rotational dynamics. Taken together, these results constitute the first report demonstrating modulation in the organization and dynamics of the N-terminal domain of a chemokine receptor in a membrane-like environment of varying hydration. We envisage that these results are relevant in the context of hydration in the function of G protein-coupled receptors.

Proton storage site in bacteriorhodopsin: new insights from quantum mechanics/molecular mechanics simulations of microscopic pK(a) and infrared spectra

Identifying the group that acts as the proton storage/loading site is a challenging but important problem for understanding the mechanism of proton pumping in biomolecular proton pumps, such as bacteriorhodopsin (bR) and cytochrome c oxidase. Recent experimental studies of bR propelled the idea that the proton storage/release group (PRG) in bR is not an amino acid but a water cluster embedded in the protein. We argue that this idea is at odds with our knowledge of protein electrostatics, since invoking the water cluster as the PRG would require the protein to raise the pK(a) of a hydronium by almost 11 pK(a) units, which is difficult considering known cases of pK(a) shifts in proteins. Our recent quantum mechanics/molecular mechanics (QM/MM) simulations suggested an alternative "intermolecular proton bond" model in which the stored proton is shared between two conserved Glu residues (194 and 204). Here we show that this model leads to microscopic pK(a) values consistent with available experimental data and the functional requirement of a PRG. Extensive QM/MM simulations also show that, independent of a number of technical issues, such as the influence of QM region size, starting X-ray structure, and nuclear quantum effects, the "intermolecular proton bond" model is qualitatively consistent with available spectroscopic data. Potential of mean force calculations show explicitly that the stored proton strongly prefers the pair of Glu residues over the water cluster. The results and analyses help highlight the importance of considering protein electrostatics and provide arguments for why the "intermolecular proton bond" model is likely applicable to the PRG in biomolecular proton pumps in general.

Crystal structures of an O-like blue form and an anion-free yellow form of pharaonis halorhodopsin

Halorhodopsin from Natronomonas pharaonis (pHR) was previously crystallized into a monoclinic space group C2, and the structure of the chloride-bound purple form was determined. Here, we report the crystal structures of two chloride-free forms of pHR, that is, an O-like blue form and an M-like yellow form. When the C2 crystal was soaked in a chloride-free alkaline solution, the protein packing was largely altered and the yellow form containing all-trans retinal was generated. Upon neutralization, this yellow form was converted into the blue form. From structural comparison of the different forms of pHR, it was shown that the removal of a chloride ion from the primary binding site (site I), which is located between the retinal Schiff base and Thr126, is accompanied by such a deformation of helix C that the side chain of Thr126 moves toward helix G, leading to a significant shrinkage of site I. A large structural change is also induced in the chloride uptake pathway, where a flip motion of the side chain of Glu234 is accompanied by large movements of the surrounding aromatic residues. Irrespective of different charge distributions at the active site, there was no large difference in the structures of the yellow form and the blue form. It is shown that the yellow-to-purple transition is initiated by the entrance of one water and one HCl to the active site, where the proton and the chloride ion in HCl are transferred to the Schiff base and site I, respectively.

Role of Arg82 in the early steps of the bacteriorhodopsin proton-pumping cycle

Proton-transfer reactions in the bacteriorhodopsin light-driven proton pump are coupled with structural rearrangements of protein amino acids and internal water molecules. It is generally thought that the first proton-transfer step from retinal Schiff base to the nearby Asp85 is coupled with movement of the Arg82 side chain away from Asp85 and toward the extracellular proton release group. This movement of Arg82 likely triggers the release of the proton from the proton release group to the extracellular bulk. The exact timing of the movement of Arg82 and how this movement is coupled with proton transfer are still not understood in molecular detail. Here, we address these questions by computing the free energy for the movement of the Arg82 side chain. The calculations indicate that protonation of Asp85 leads to a fast reorientation of the Arg82 side chain toward the extracellular proton release group.

X-ray-radiation-induced changes in bacteriorhodopsin structure

Bacteriorhodopsin (bR) provides light-driven vectorial proton transport across a cell membrane. Creation of electrochemical potential at the membrane is a universal step in energy transformation in a cell. Published atomic crystallographic models of early intermediate states of bR show a significant difference between them, and conclusions about pumping mechanisms have been contradictory. Here, we present a quantitative high-resolution crystallographic study of conformational changes in bR induced by X-ray absorption. It is shown that X-ray doses that are usually accumulated during data collection for intermediate-state studies are sufficient to significantly alter the structure of the protein. X-ray-induced changes occur primarily in the active site of bR. Structural modeling showed that X-ray absorption triggers retinal isomerization accompanied by the disappearance of electron densities corresponding to the water molecule W402 bound to the Schiff base. It is demonstrated that these and other X-ray-induced changes may mimic functional conformational changes of bR leading to misinterpretation of the earlier obtained X-ray crystallographic structures of photointermediates.

Structural divergence and functional versatility of the rhodopsin superfamily

Seven-transmembrane-helix retinylidene proteins, which constitute the rhodopsin superfamily, have been discovered in diverse species, including Archaea, Eubacteria, fungi, algae and animals. Some members of this super-family were specialized to function as light-driven proton pumps, light-driven chloride pumps, photoisomerases, or light-gated ion channels, where the photochemical reactions are self-completed without interactions with other proteins. Other members evolved to acquire the ability to modulate soluble cytoplasmic or membrane-embedded signal transducers. During the last decade, high-resolution crystal structures were reported for ten members of the rhodopsin superfamily; viz., four proton pumps, two chloride pumps, two microbial photosensors and two visual pigments. Comparison of these structures provides us with a hint to elucidate the common structural motif that is utilized to stabilize their tertiary structures as well as unique architectures that are relevant to specific functions.

Helical membrane peptides to modulate cell function

In recent years there has been an abundance of research into the potential of helical peptides to influence cell function. These peptides have been used to achieve a variety of different outcomes from cell repair to cell death, depending upon the peptide sequence and the nature of its interactions with cell membranes and membrane proteins. In this critical review, we summarise several mechanisms by which helical peptides, acting as either transporters, inhibitors, agonists or antibiotics, can have significant effects on cell membranes and can radically affect the internal mechanisms of the cell. The various approaches to peptide design are discussed, including the role of naturally-occurring proteins in the design of these helical peptides and current breakthroughs in the use of non-natural (and therefore more stable) peptide scaffolds. Most importantly, the current successful applications of these peptides, and their potential uses in the field of medicine, are reviewed (131 references).

Time-resolved structural studies of protein reaction dynamics: a smorgasbord of X-ray approaches

Time-resolved structural studies of proteins have undergone several significant developments during the last decade. Recent developments using time-resolved X-ray methods, such as time-resolved Laue diffraction, low-temperature intermediate trapping, time-resolved wide-angle X-ray scattering and time-resolved X-ray absorption spectroscopy, are reviewed.

Proteins undergo conformational changes during their biological function. As such, a high-resolution structure of a protein’s resting conformation provides a starting point for elucidating its reaction mechanism, but provides no direct information concerning the protein’s conformational dynamics. Several X-ray methods have been developed to elucidate those conformational changes that occur during a protein’s reaction, including time-resolved Laue diffraction and intermediate trapping studies on three-dimensional protein crystals, and time-resolved wide-angle X-ray scattering and X-ray absorption studies on proteins in the solution phase. This review emphasizes the scope and limitations of these complementary experimental approaches when seeking to understand protein conformational dynamics. These methods are illustrated using a limited set of examples including myoglobin and haemoglobin in complex with carbon monoxide, the simple light-driven proton pump bacteriorhodopsin, and the superoxide scavenger superoxide reductase. In conclusion, likely future developments of these methods at synchrotron X-ray sources and the potential impact of emerging X-ray free-electron laser facilities are speculated upon.