Crystallographic structures of the M and N intermediates of bacteriorhodopsin: assembly of a hydrogen-bonded chain of water molecules between Asp-96 and the retinal Schiff base

Determinants of hydrogen bond distances in proteins

Hydrogen bonds (H-bonds) between oxygen atoms, with the O-H bond donated to the acceptor O atom (Odonor-H⋯Oacceptor), are essential for stabilizing protein structures and facilitating enzymatic reactions. The dielectric and electrostatic environment of proteins, as well as structural constraints imposed by protein folding, influence the nature of H-bonds. In this study, we investigated how these factors affect H-bond distances in proteins. Analysis of 906 high-resolution protein structures (≤1.2 Å) from the Protein Data Bank revealed that H-bond distances for H-bonds with the same donor and acceptor groups are distributed around a value primarily determined by the pKa difference between these groups (ΔpKa) in water, with lower ΔpKa values leading to shorter distances. This correlation arises from enhanced electron redistribution from the H-bond acceptor to the donor in lower ΔpKa H-bonds, which increases the covalent character of the H-bond and decreases the H⋯Oacceptor distance. In contrast, H-bond distances are largely unaffected by whether the H-bond is buried in the protein interior or exposed to bulk water, as the strength of the electrostatic interaction between the donor and acceptor groups plays a minor role in determining distances. Furthermore, analysis of H-bonds in microbial rhodopsins using a quantum mechanical/molecular mechanical approach demonstrates that the protein environment primarily influences H-bond distances electrostatically by altering the ΔpKa of the H-bond, while structural constraints impose a secondary influence by altering Odonor-H⋯Oacceptor angles or H⋯Oacceptor distances without changing ΔpKa.

Weak Organic Acid Effect of Bacterial Light-Driven Proton-Pumping Rhodopsin

Microbial rhodopsins are photoreceptor proteins that utilize light to elicit various biological functions. The best-studied microbial rhodopsins are outward proton (H+)-pumping rhodopsins, which transport H+ from the cytoplasmic to the extracellular side. Recently, the weak organic acid (WOA) effect, specifically the enhancement of pumping activity by WOAs such as acetic acid and indole-3-acetic acid (IAA), was discovered in outward H+-pumping rhodopsins from fungi. However, it remains unclear whether the WOA effect exists in nonfungal H+-pumping rhodopsins. Here, we revealed that the H+-pumping activity of a bacterial outward H+ pump rhodopsin, PspR, from the rhizobacterium Pseudomonas putida, is also enhanced by extracellular acetic acid and IAA. Using transient absorption measurements on purified PspR protein, we found that extracellular WOAs accelerate cytoplasmic H+ uptake and extracellular H+ release from a protonated counterion during its photocycle. Furthermore, acetic acid applied on the cytoplasmic side has an inhibitory effect on the H+ pump activity of PspR, which is less significant for IAA and can be mitigated by increasing the H+ concentration or introducing a cytoplasmic donor residue. These findings on the WOA effect in a bacterial rhodopsin provide new insights into the physiological function of outward H+-pumping rhodopsins in bacteria, particularly in their interaction with plants.

Stretching vibrational frequencies and pKa differences in H-bond networks of protein environments

The experimentally measured stretching vibrational frequencies of O-D [νO-D(donor)] and C=O [νC=O(donor)] H-bond donor groups can provide valuable information about the H-bonds in proteins. Here, using a quantum mechanical/molecular mechanical approach, the relationship between these vibrational frequencies and the difference in pKa values between H-bond donor and acceptor groups [ΔpKa(donor … acceptor)] in bacteriorhodopsin and photoactive yellow protein environments was investigated. The results show that νO-D(donor) is correlated with ΔpKa(donor … acceptor), regardless of the specific protein environment. νC=O(donor) is also correlated with ΔpKa(donor … acceptor), although the correlation is weak because the C=O bond does not have a proton. Importantly, the shifts in νO-D(donor) and νC=O(donor) are not caused by changes in pKa(donor) alone, but rather by changes in ΔpKa(donor … acceptor). Specifically, a decrease in ΔpKa(donor … acceptor) can lead to proton release from the H-bond donor group toward the acceptor group, resulting in shifts in the vibrational frequencies of the protein environment. These findings suggest that changes in the stretching vibrational frequencies, in particular νO-D(donor), can be used to monitor proton transfer in protein environments.

Voltage Imaging with Engineered Proton-Pumping Rhodopsins: Insights from the Proton Transfer Pathway

Voltage imaging using genetically encoded voltage indicators (GEVIs) has taken the field of neuroscience by storm in the past decade. Its ability to create subcellular and network level readouts of electrical dynamics depends critically on the kinetics of the response to voltage of the indicator used. Engineered microbial rhodopsins form a GEVI subclass known for their high voltage sensitivity and fast response kinetics. Here we review the essential aspects of microbial rhodopsin photocycles that are critical to understanding the mechanisms of voltage sensitivity in these proteins and link them to insights from efforts to create faster, brighter and more sensitive microbial rhodopsin-based GEVIs.

Structural studies of bacteriorhodopsin in BC era

It marked half a century since the discovery of bacteriorhodopsin two years ago. On this occasion, I have revisited historically important diffraction studies of this membrane protein, based on my recollections. X-ray diffraction and electron diffraction, and electron microscopy, described the low-resolution structure of bacteriorhodopsin within the purple membrane. Neutron diffraction was effective to assign the helical regions in the primary structure with 7 rods revealed by low-resolution structure as well as to describe the retinal position. Substantial conformational changes upon light illumination were clarified by the structures of various photointermediates. Early trials of time-resolved studies were also introduced. Models for the mechanism of light-driven proton pump based on the low-resolution structural studies are also described. Significantly, they are not far from the today's understanding. I believe that the spirit of the early research scientists in this field and the essence of their studies, which constitute the foundations of the field, still actively fertilizes current membrane protein research.

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.

Kinetic study on the molecular mechanism of light-driven inward proton transport by schizorhodopsins

Schizorhodopsins (SzRs) are light-driven inward proton pumping membrane proteins. A H⁺ is released to the cytoplasmic solvent from the chromophore, retinal Schiff base (RSB), after light absorption, and then another H⁺ is bound to the RSB at the end of photocyclic reaction. However, the mechanistic detail of H⁺ transfers in SzR is almost unknown. Here we studied the deuterium isotope effect and the temperature dependence of the reaction rate constants of elementary steps in the photocycles of SzRs. The former indicated that deprotonation and reprotonation of RSB is mainly accomplished by H⁺ hopping between heavy atoms with similar H⁺ affinity. Furthermore, the temperature dependence of the rate constants revealed that most of H⁺ transfer events have a high entropy barrier. In contrast, the activation enthalpy and entropy of extremely thermostable SzR (MsSzR) are significantly higher than other types of SzRs (SzR1 and MtSzR) suggesting that its highly thermostable structure is optimized with at the cost of slower reaction rates at ambient temperatures.

Common Mechanism of Activated Catalysis in P-loop Fold Nucleoside Triphosphatases-United in Diversity

To clarify the obscure hydrolysis mechanism of ubiquitous P-loop-fold nucleoside triphosphatases (Walker NTPases), we analysed the structures of 3136 catalytic sites with bound Mg-NTP complexes or their analogues. Our results are presented in two articles; here, in the second of them, we elucidated whether the Walker A and Walker B sequence motifs-common to all P-loop NTPases-could be directly involved in catalysis. We found that the hydrogen bonds (H-bonds) between the strictly conserved, Mg-coordinating Ser/Thr of the Walker A motif ([Ser/Thr]WA) and aspartate of the Walker B motif (AspWB) are particularly short (even as short as 2.4 ångströms) in the structures with bound transition state (TS) analogues. Given that a short H-bond implies parity in the pKa values of the H-bond partners, we suggest that, in response to the interactions of a P-loop NTPase with its cognate activating partner, a proton relocates from [Ser/Thr]WA to AspWB. The resulting anionic [Ser/Thr]WA alkoxide withdraws a proton from the catalytic water molecule, and the nascent hydroxyl attacks the gamma phosphate of NTP. When the gamma-phosphate breaks away, the trapped proton at AspWB passes by the Grotthuss relay via [Ser/Thr]WA to beta-phosphate and compensates for its developing negative charge that is thought to be responsible for the activation barrier of hydrolysis.

Proton transfer reactions in donor site mutants of ESR, a retinal protein from Exiguobacterium sibiricum

Light-driven proton transport by microbial retinal proteins such as archaeal bacteriorhodopsin involves carboxylic residues as internal proton donors to the catalytic center which is a retinal Schiff base (SB). The proton donor, Asp96 in bacteriorhodopsin, supplies a proton to the transiently deprotonated Schiff base during the photochemical cycle. Subsequent proton uptake resets the protonated state of the donor. This two step process became a distinctive signature of retinal based proton pumps. Similar steps are observed also in many natural variants of bacterial proteorhodopsins and xanthorhodopsins where glutamic acid residues serve as a proton donor. Recently, however, an exception to this rule was found. A retinal protein from Exiguobacterium sibiricum, ESR, contains a Lys residue in place of Asp or Glu, which facilitates proton transfer from the bulk to the SB. Lys96 can be functionally replaced with the more common donor residues, Asp or Glu. Proton transfer to the SB in the mutants containing these replacements (K96E and K96D/A47T) is much faster than in the proteins lacking the proton donor (K96A and similar mutants), and in the case of K96D/A47T, comparable with that in the wild type, indicating that carboxylic residues can replace Lys96 as proton donors in ESR. We show here that there are important differences in the functioning of these residues in ESR from the way Asp96 functions in bacteriorhodopsin. Reprotonation of the SB and proton uptake from the bulk occur almost simultaneously during the M to N transition (as in the wild type ESR at neutral pH), whereas in bacteriorhodopsin these two steps are well separated in time and occur during the M to N and N to O transitions, respectively.

Microsolvation Effects in the Spectral Tuning of Heliorhodopsin

Heliorhodopsins (HeR) are a new category of heptahelical transmembrane photoactive proteins with a covalently linked all-trans retinal. The protonated Schiff base (PSB) nitrogen in the retinal is stabilized by a negatively charged counterion. It is well-known that stronger or weaker electrostatic interactions with the counterion cause a significant spectral blue- or red-shift, respectively, in both microbial and animal rhodopsins. In HeR, however, while Glu107 acts as the counterion, mutations of this residue are not directly correlated with a spectral shift. A molecular dynamics analysis revealed that a water cluster pocket produces a microsolvation effect on the Schiff base, compensating to various extents the replacement of the native counterion. Using a combination of molecular dynamics and quantum mechanical/molecular mechanics (QM/MM), we study this microsolvation effect on the electronic absorption of the retinylidene Schiff base chromophore of HeR.

Hydroxide Ion Mechanism for Long-Range Proton Pumping in the Third Proton Transfer of Bacteriorhodopsin

In bacteriorhodopsin, representative light-driven proton pump, five proton transfers yield vectorial active proton translocation, resulting in a proton gradient in microbes. Third proton transfer occurs from Asp96 to the Schiff base on the photocycle, which is expected to be a long-range proton transfer via the Grotthuss mechanism through internal water molecules. Here, large-scale quantum molecular dynamics simulations are performed for the third proton transfer, where all the atoms (~50000 atoms) are treated quantum-mechanically. The simulations demonstrate that two reaction paths exist along the water wire, namely, via hydronium and via hydroxide ions. The free energy analysis confirms that the path via hydroxide ions is considerably favorable and consistent with the observed lifetime of the transient water wire. Therefore, the proposed hydroxide ion mechanism, as in the first proton transfer, is responsible for the third long-range proton transfer.

Correlation between C═O Stretching Vibrational Frequency and pKa Shift of Carboxylic Acids

Identifying the pK_a values of aspartic acid (Asp) and glutamic acid (Glu) in active sites is essential for understanding enzyme reaction mechanisms. In this study, we investigated the correlation between the C=O stretching vibrational frequency (ν_C=O) of protonated carboxylic acids and the pK_a values using density functional theory calculations. In unsaturated carboxylic acids (e.g., benzoic acid analogues), ν_C=O decreases as the pK_a increases (the negative correlation), whereas in saturated carboxylic acids (e.g., acetic acid analogues, Asp, and Glu), ν_C=O increases as the pK_a increases (the positive correlation) as long as the structure of the H-bond network around the acid is identical. The negative/positive correlation between ν_C=O and pK_a can be rationalized by the presence or absence of the C=C double bond. The pK_a shift was estimated from the ν_C=O shift of Asp and Glu in proteins on the basis of the negative correlation derived from benzoic acids. The previous estimations should be revisited by using the positive correlation derived in this study, as demonstrated by quantum mechanical/molecular mechanical calculations of ν_C=O and electrostatic calculations of pK_a on a key Asp85 in the proton-transfer pathway of bacteriorhodopsin.

Archaeal Lipids Regulating the Trimeric Structure Dynamics of Bacteriorhodopsin for Efficient Proton Release and Uptake

S-TGA-1 and PGP-Me are native archaeal lipids associated with the bacteriorhodopsin (bR) trimer and contribute to protein stabilization and native dynamics for proton transfer. However, little is known about the underlying molecular mechanism of how these lipids regulate bR trimerization and efficient photocycling. Here, we explored the specific binding of S-TGA-1 and PGP-Me with the bR trimer and elucidated how specific interactions modulate the bR trimeric structure and proton release and uptake using long-term atomistic molecular dynamic simulations. Our results showed that S-TGA-1 and PGP-Me are essential for stabilizing the bR trimer and maintaining the coherent conformational dynamics necessary for proton transfer. The specific binding of S-TGA-1 with W80 and K129 regulates proton release on the extracellular surface by forming a “Glu-shared” model. The interaction of PGP-Me with K40 ensures proton uptake by accommodating the conformation of the helices to recruit enough water molecules on the cytoplasmic side. The present study results could fill in the theoretical gaps of studies on the functional role of archaeal lipids and could provide a reference for other membrane proteins containing similar archaeal lipids.

Crystallographic Studies of Rhodopsins: Structure and Dynamics

Crystal structures have provided detailed insight in the architecture of rhodopsin photoreceptors. Of particular interest are the protein-chromophore interactions that govern the light-induced retinal isomerization and ultimately induce the large structural changes important for the various biological functions of the family. The reaction intermediates occurring along the rhodopsin photocycle have vastly differing lifetimes, from hundreds of femtoseconds to milliseconds. Detailed insight at high spatial and temporal resolution can be obtained by time-resolved crystallography using pump-probe approaches at X-ray free-electron lasers. Alternatively, cryotrapping approaches can be used. Both the approaches are described, including illumination and sample delivery. The importance of appropriate photoexcitation avoiding multiphoton absorption is stressed.

Dynamic Coupling of Tyrosine 185 with the Bacteriorhodopsin Photocycle, as Revealed by Chemical Shifts, Assisted AF-QM/MM Calculations and Molecular Dynamic Simulations

Aromatic residues are highly conserved in microbial photoreceptors and play crucial roles in the dynamic regulation of receptor functions. However, little is known about the dynamic mechanism of the functional role of those highly conserved aromatic residues during the receptor photocycle. Tyrosine 185 (Y185) is one of the highly conserved aromatic residues within the retinal binding pocket of bacteriorhodopsin (bR). In this study, we explored the molecular mechanism of its dynamic coupling with the bR photocycle by automated fragmentation quantum mechanics/molecular mechanics (AF-QM/MM) calculations and molecular dynamic (MD) simulations based on chemical shifts obtained by 2D solid-state NMR correlation experiments. We observed that Y185 plays a significant role in regulating the retinal cis–trans thermal equilibrium, stabilizing the pentagonal H-bond network, participating in the orientation switch of Schiff Base (SB) nitrogen, and opening the F42 gate by interacting with the retinal and several key residues along the proton translocation channel. Our findings provide a detailed molecular mechanism of the dynamic couplings of Y185 and the bR photocycle from a structural perspective. The method used in this paper may be applied to the study of other microbial photoreceptors.

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.

Interactive Interface for Graph-Based Analyses of Dynamic H-Bond Networks: Application to Spike Protein S

Dynamic hydrogen-bond networks are key determinants of protein conformational dynamics. In the case of macromolecular protein complexes, which can have a large number of hydrogen bonds giving rise to extensive hydrogen-bond networks, efficient algorithms are required to analyze interactions that could be important for the dynamics and biological function of the complex. We present here a highly efficient, standalone interface designed for analyses of dynamical hydrogen-bond networks of biomolecules and macromolecular complexes. To facilitate a comprehensive description of protein dynamics, the interface includes analyses of hydrophobic interactions. We illustrate the usefulness and workflow of the interface by dissecting the dynamics of the ectodomain of SARS-CoV-2 protein S in its closed conformation. We find that protein S contains numerous local clusters of dynamic hydrogen bonds and identify hydrogen bonds that are sampled persistently. The receptor binding domain of the spike protein hosts only a handful of persistent hydrogen-bond clusters, suggesting structural plasticity. Our data analysis interface is released here for open use.

Photoactivated voltage imaging in tissue with an archaerhodopsin-derived reporter

Photoactivated genetically encoded voltage indicators (GEVIs) have the potential to enable optically sectioned voltage imaging at the intersection of a photoactivation beam and an imaging beam. We developed a pooled high-throughput screen to identify archaerhodopsin mutants with enhanced photoactivation. After screening ~105 cells, we identified a novel GEVI, NovArch, whose one-photon near-infrared fluorescence is reversibly enhanced by weak one-photon blue or two-photon near-infrared excitation. Because the photoactivation leads to fluorescent signals catalytically rather than stoichiometrically, high fluorescence signals, optical sectioning, and high time resolution are achieved simultaneously at modest blue or two-photon laser power. We demonstrate applications of the combined molecular and optical tools to optical mapping of membrane voltage in distal dendrites in acute mouse brain slices and in spontaneously active neurons in vivo.

Catalytic Reactions and Energy Conservation in the Cytochrome bc1 and b6f Complexes of Energy-Transducing Membranes

This review focuses on key components of respiratory and photosynthetic energy-transduction systems: the cytochrome bc1 and b6f (Cytbc1/b6f) membranous multisubunit homodimeric complexes. These remarkable molecular machines catalyze electron transfer from membranous quinones to water-soluble electron carriers (such as cytochromes c or plastocyanin), coupling electron flow to proton translocation across the energy-transducing membrane and contributing to the generation of a transmembrane electrochemical potential gradient, which powers cellular metabolism in the majority of living organisms. Cytsbc1/b6f share many similarities but also have significant differences. While decades of research have provided extensive knowledge on these enzymes, several important aspects of their molecular mechanisms remain to be elucidated. We summarize a broad range of structural, mechanistic, and physiological aspects required for function of Cytbc1/b6f, combining textbook fundamentals with new intriguing concepts that have emerged from more recent studies. The discussion covers but is not limited to (i) mechanisms of energy-conserving bifurcation of electron pathway and energy-wasting superoxide generation at the quinol oxidation site, (ii) the mechanism by which semiquinone is stabilized at the quinone reduction site, (iii) interactions with substrates and specific inhibitors, (iv) intermonomer electron transfer and the role of a dimeric complex, and (v) higher levels of organization and regulation that involve Cytsbc1/b6f. In addressing these topics, we point out existing uncertainties and controversies, which, as suggested, will drive further research in this field.

Proton uptake mechanism in bacteriorhodopsin captured by serial synchrotron crystallography

Conformational dynamics are essential for proteins to function. We adapted time-resolved serial crystallography developed at x-ray lasers to visualize protein motions using synchrotrons. We recorded the structural changes in the light-driven proton-pump bacteriorhodopsin over 200 milliseconds in time. The snapshot from the first 5 milliseconds after photoactivation shows structural changes associated with proton release at a quality comparable to that of previous x-ray laser experiments. From 10 to 15 milliseconds onwards, we observe large additional structural rearrangements up to 9 angstroms on the cytoplasmic side. Rotation of leucine-93 and phenylalanine-219 opens a hydrophobic barrier, leading to the formation of a water chain connecting the intracellular aspartic acid-96 with the retinal Schiff base. The formation of this proton wire recharges the membrane pump with a proton for the next cycle.

Light-induced effects in glycine aqueous solution studied by Fourier transform infrared-emission spectroscopy and ultraviolet-visible spectroscopy

Although amino acids are insensitive to visible light, as is generally accepted, we show that particular light-matter interaction can break this obviousness. Using sensitive (FT-IR)-technique in a combination of a broadband visible light source, we registered emission spectra of glycine in the range 2500-500 cm^-1. Sensitivity of the infrared emission spectrum to the exciting power -induced changes in the glycine structure was demonstrated experimentally.Vibrational spectra of glycine displayed the prominent spectral features of CH₂, COO-, COOH, NH⁺₃ groups in the "fingerprint region". Simultaneous appearance of ionised COO- and unionised COOH forms of glycine in solution at neutral pH clearly indicated that visible light induces the partial protonation of COO- groups; if so, visible light irradiation should lead to occurrence of dimers or dimeric hydrogen - bonded structures. Spectroscopic and microscopic evidence of visible light-mediated formation of aggregates and nucleus in aqueous solution was presented.Electronic absorption/emission spectra of glycine in aqueous solution were primarily characterized in the near ultraviolet-visible region (240-600 nm). Negligible absorption near 270 nm was observed for a 1.0 M solution and dramatically enhanced with its "aging". Moreover, an extension of the absorption edge into the region above 400 nm could be seen. Due to the visible light irradiation, we observed modification of electronic structure or occurrence of additional species causing changes in absorption of glycine amino acid. For "aged" solution, it was shown that excitation spectra corresponding to the different emission wavelengths were entirely different, at that each excitation-spectral band had a characteristic emission band.Communicated by Ramaswamy H. Sarma.

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.

The Role of Proton Transport in Gating Current in a Voltage Gated Ion Channel, as Shown by Quantum Calculations

Over two-thirds of a century ago, Hodgkin and Huxley proposed the existence of voltage gated ion channels (VGICs) to carry Na⁺ and K⁺ ions across the cell membrane to create the nerve impulse, in response to depolarization of the membrane. The channels have multiple physiological roles, and play a central role in a wide variety of diseases when they malfunction. The first channel structure was found by MacKinnon and coworkers in 1998. Subsequently, the structure of a number of VGICs was determined in the open (ion conducting) state. This type of channel consists of four voltage sensing domains (VSDs), each formed from four transmembrane (TM) segments, plus a pore domain through which ions move. Understanding the gating mechanism (how the channel opens and closes) requires structures. One TM segment (S4) has an arginine in every third position, with one such segment per domain. It is usually assumed that these arginines are all ionized, and in the resting state are held toward the intracellular side of the membrane by voltage across the membrane. They are assumed to move outward (extracellular direction) when released by depolarization of this voltage, producing a capacitive gating current and opening the channel. We suggest alternate interpretations of the evidence that led to these models. Measured gating current is the total charge displacement of all atoms in the VSD; we propose that the prime, but not sole, contributor is proton motion, not displacement of the charges on the arginines of S4. It is known that the VSD can conduct protons. Quantum calculations on the K_v1.2 potassium channel VSD show how; the key is the amphoteric nature of the arginine side chain, which allows it to transfer a proton. This appears to be the first time the arginine side chain has had its amphoteric character considered. We have calculated one such proton transfer in detail: this proton starts from a tyrosine that can ionize, transferring to the NE of the third arginine on S4; that arginine’s NH then transfers a proton to a glutamate. The backbone remains static. A mutation predicted to affect the proton transfer has been qualitatively confirmed experimentally, from the change in the gating current-voltage curve. The total charge displacement in going from a normal closed potential of −70 mV across the membrane to 0 mV (open), is calculated to be approximately consistent with measured values, although the error limits on the calculation require caution in interpretation.

Allosteric Effects of the Proton Donor on the Microbial Proton Pump Proteorhodopsin

Proteorhodopsin (PR) is a microbial proton pump that is ubiquitous in marine environments and may play an important role in the oceanic carbon cycle. Photoisomerization of the retinal chromophore in PR leads to a series of proton transfers between specific acidic amino acid residues and the Schiff base of retinal, culminating in a proton motive force to facilitate ATP synthesis. The proton donor in a similar retinal protein, bacteriorhodopsin, acts as a latch to allow the influx of bulk water. However, it is unclear if the proton donor in PR, E108, utilizes the same latch mechanism to become internally hydrated. Here, we used molecular dynamics simulations to model the changes in internal hydration of the blue variant of PR during photoactivation with the proton donor in protonated and deprotonated states. We find that there is a stark contrast in the levels of internal hydration of the cytoplasmic half of PR based on the protonation state of E108. Instead of a latch mechanism, deprotonation of E108 acts as a gate, taking advantage of a nearby polar residue (S61) to promote the formation of a stable water wire from bulk cytoplasm to the retinal-binding pocket over hundreds of nanoseconds. No large-scale conformational changes occur in PR over the microsecond timescale. This subtle yet clear difference in the effect of deprotonation of the proton donor in PR may help explain why the photointermediates that involve the proton donor (i.e., M and N states) have timescales that are orders of magnitude different from the archaeal proton pump, bacteriorhodopsin. In general, our study highlights the importance of understanding how structural fluctuations lead to differences in the way that retinal proteins accomplish the same task.

Catalysis of Ground State cis[Formula: see text] trans Isomerization of Bacteriorhodopsin's Retinal Chromophore by a Hydrogen-Bond Network

For the photocycle of the membrane protein bacteriorhodopsin to proceed efficiently, the thermal 13-cis to all-trans back-isomerization of the retinal chromophore must return the protein to its resting state on a time-scale of milliseconds. Here, we report on quantum mechanical/molecular mechanical energy calculations examining the structural and energetic determinants of the retinal cis-trans isomerization in the protein environment. The results suggest that a hydrogen-bonded network consisting of the retinal Schiff base, active site amino acid residues, and water molecules can stabilize the twisted retinal, thus reducing the intrinsic energy cost of the cis-trans thermal isomerization barrier.

Bacteriorhodopsin-like channelrhodopsins: Alternative mechanism for control of cation conductance

The recently discovered cation-conducting channelrhodopsins in cryptophyte algae are far more homologous to haloarchaeal rhodopsins, in particular the proton pump bacteriorhodopsin (BR), than to earlier known channelrhodopsins. They uniquely retain the two carboxylate residues that define the vectorial proton path in BR in which Asp-85 and Asp-96 serve as acceptor and donor, respectively, of the photoactive site Schiff base (SB) proton. Here we analyze laser flash-induced photocurrents and photochemical conversions in Guillardia theta cation channelrhodopsin 2 (GtCCR2) and its mutants. Our results reveal a model in which the GtCCR2 retinylidene SB chromophore rapidly deprotonates to the Asp-85 homolog, as in BR. Opening of the cytoplasmic channel to cations in GtCCR2 requires the Asp-96 homolog to be unprotonated, as has been proposed for the BR cytoplasmic channel for protons. However, reprotonation of the GtCCR2 SB occurs not from the Asp-96 homolog, but by proton return from the earlier protonated acceptor, preventing vectorial proton translocation across the membrane. In GtCCR2, deprotonation of the Asp-96 homolog is required for cation channel opening and occurs >10-fold faster than reprotonation of the SB, which temporally correlates with channel closing. Hence in GtCCR2, cation channel gating is tightly coupled to intramolecular proton transfers involving the same residues that define the vectorial proton path in BR.

Infrared spectral marker bands characterizing a transient water wire inside a hydrophobic membrane protein

Proton conduction along protein-bound "water wires" is an essential feature in membrane proteins. Here, we analyze in detail a transient water wire, which conducts protons via a hydrophobic barrier within a membrane protein to create a proton gradient. It is formed only for a millisecond out of three water molecules distributed at inactive positions in a polar environment in the ground state. The movement into a hydrophobic environment causes characteristic shifts of the water bands reflecting their different chemical properties. These band shifts are identified by time-resolved Fourier Transform Infrared difference spectroscopy and analyzed by biomolecular Quantum Mechanical/Molecular Mechanical simulations. A non-hydrogen bonded ("dangling") O-H stretching vibration band and a broad continuum absorbance caused by a combined vibration along the water wire are identified as characteristic marker bands of such water wires in a hydrophobic environment. The results provide a basic understanding of water wires in hydrophobic environments.

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.

Altered backbone and side-chain interactions result in route heterogeneity during the folding of interleukin-1β (IL-1β)

Deletion of the β-bulge trigger-loop results in both a switch in the preferred folding route, from the functional loop packing folding route to barrel closure, as well as conversion of the agonist activity of IL-1β into antagonist activity. Conversely, circular permutations of IL-1β conserve the functional folding route as well as the agonist activity. These two extremes in the folding-functional interplay beg the question of whether mutations in IL-1β would result in changes in the populations of heterogeneous folding routes and the signaling activity. A series of topologically equivalent water-mediated β-strand bridging interactions within the pseudosymmetric β-trefoil fold of IL-1β highlight the backbone water interactions that stabilize the secondary and tertiary structure of the protein. Additionally, conserved aromatic residues lining the central cavity appear to be essential for both stability and folding. Here, we probe these protein backbone-water molecule and side chain-side chain interactions and the role they play in the folding mechanism of this geometrically stressed molecule. We used folding simulations with structure-based models, as well as a series of folding kinetic experiments to examine the effects of the F42W core mutation on the folding landscape of IL-1β. This mutation alters water-mediated backbone interactions essential for maintaining the trefoil fold. Our results clearly indicate that this perturbation in the primary structure alters a structural water interaction and consequently modulates the population of folding routes accessed during folding and signaling activity.

Stable closure of the cytoplasmic half-channel is required for efficient proton transport at physiological membrane potentials in the bacteriorhodopsin catalytic cycle

The bacteriorhodopsin (BR) Asp96Gly/Phe171Cys/Phe219Leu triple mutant has been shown to translocate protons 66% as efficiently as the wild-type protein. Light-dependent ATP synthesis in haloarchaeal cells expressing the triple mutant is 85% that of the wild-type BR expressing cells. Therefore, the functional activity of BR seems to be largely preserved in the triple mutant despite the observations that its ground-state structure resembles that of the wild-type M state (i.e., the so-called cytoplasmically open state) and that the mutant shows no significant structural changes during its photocycle, in sharp contrast to what occurs in the wild-type protein in which a large structural opening and closing occurs on the cytoplasmic side. To resolve the contradiction between the apparent functional robustness of the triple mutant and the presumed importance of the opening and closing that occurs in the wild-type protein, we conducted additional experiments to compare the behavior of wild-type and mutant proteins under different operational loads. Specifically, we characterized the ability of the two proteins to generate light-driven proton currents against a range of membrane potentials. The wild-type protein showed maximal conductance between −150 and −50 mV, whereas the mutant showed maximal conductance at membrane potentials >+50 mV. Molecular dynamics (MD) simulations of the triple mutant were also conducted to characterize structural changes in the protein and in solvent accessibility that might help to functionally contextualize the current–voltage data. These simulations revealed that the cytoplasmic half-channel of the triple mutant is constitutively open and dynamically exchanges water with the bulk. Collectively, the data and simulations help to explain why this mutant BR does not mediate photosynthetic growth of haloarchaeal cells, and they suggest that the structural closing observed in the wild-type protein likely plays a key role in minimizing substrate back flow in the face of electrochemical driving forces present at physiological membrane potentials.

Role of Thr218 in the light-driven anion pump halorhodopsin from Natronomonas pharaonis

Halorhodopsin (HR) is an inward-directed light-driven halogen ion pump, and NpHR is a HR from Natronomonas pharaonis. Unphotolyzed NpHR binds halogen ion in the vicinity of the Schiff base, which links retinal to Lys256. This halogen ion is transported during the photocycle. We made various mutants of Thr218, which is located one half-turn up from the Schiff base to the cytoplasm (CP) channel, and analyzed the photocycle using a sequential irreversible model. Four photochemically defined intermediates (P(i), i = 1-4) were adequate to describe the photocycle. The third component, P₃, was a quasi-equilibrium complex between the N and O intermediates, where a N ↔ O + Cl⁻ equilibrium was attained. The K(d,N↔O) values of this equilibrium for various mutants were determined, and the value of Thr (wild type) was the highest. The partial molar volume differences between N and O, ΔV(N→O), were estimated from the pressure dependence of K(d,N↔O). A comparison between K(d,N↔O) and ΔV(N→O) led to the conclusion that water entry by the F-helix opening at O may occur, which may increase K(d,N↔O). For some mutants, however, large ΔV(N→O) values were found, whereas the K(d,N↔O) values were small. This suggests that the special coordination of a water molecule with the OH group of Thr is necessary for the increase in K(d,N↔O). Mutants with a small K(d,N↔O) showed low pumping activities in the presence of inside negative membrane potential, while the mutant activities were not different in the absence of membrane potential. The effect of the mutation on the pumping activities is discussed.

Gloeobacter rhodopsin, limitation of proton pumping at high electrochemical load

We studied the photocurrents of a cyanobacterial rhodopsin Gloeobacter violaceus (GR) in Xenopus laevis oocytes and HEK-293 cells. This protein is a light-driven proton pump with striking similarities to marine proteorhodopsins, including the D121-H87 cluster of the retinal Schiff base counterion and a glutamate at position 132 that acts as a proton donor for chromophore reprotonation during the photocycle. Interestingly, at low extracellular pH(o) and negative voltage, the proton flux inverted and directed inward. Using electrophysiological measurements of wild-type and mutant GR, we demonstrate that the electrochemical gradient limits outward-directed proton pumping and converts it into a purely passive proton influx. This conclusion contradicts the contemporary paradigm that at low pH, proteorhodopsins actively transport H(+) into cells. We identified E132 and S77 as key residues that allow inward directed diffusion. Substitution of E132 with aspartate or S77 with either alanine or cysteine abolished the inward-directed current almost completely. The proton influx is likely caused by the pK(a) of E132 in GR, which is lower than that of other microbial ion pumping rhodopsins. The advantage of such a low pK(a) is an acceleration of the photocycle and high pump turnover at high light intensities.

The role of protein-bound water molecules in microbial rhodopsins

Protein-bound internal water molecules are essential features of the structure and function of microbial rhodopsins. Besides structural stabilization, they act as proton conductors and even proton storage sites. Currently, the most understood model system exhibiting such features is bacteriorhodopsin (bR). During the last 20 years, the importance of water molecules for proton transport has been revealed through this protein. It has been shown that water molecules are as essential as amino acids for proton transport and biological function. In this review, we present an overview of the historical development of this research on bR. We furthermore summarize the recently discovered protein-bound water features associated with proton transport. Specifically, we discuss a pentameric water/amino acid arrangement close to the protonated Schiff base as central proton-binding site, a protonated water cluster as proton storage site at the proton-release site, and a transient linear water chain at the proton uptake site. We highlight how protein conformational changes reposition or reorient internal water molecules, thereby guiding proton transport. Last, we compare the water positions in bR with those in other microbial rhodopsins to elucidate how protein-bound water molecules guide the function of microbial rhodopsins. This article is part of a Special Issue entitled: Retinal Proteins - You can teach an old dog new tricks.

Molecular mechanisms for generating transmembrane proton gradients

Membrane proteins use the energy of light or high energy substrates to build a transmembrane proton gradient through a series of reactions leading to proton release into the lower pH compartment (P-side) and proton uptake from the higher pH compartment (N-side). This review considers how the proton affinity of the substrates, cofactors and amino acids are modified in four proteins to drive proton transfers. Bacterial reaction centers (RCs) and photosystem II (PSII) carry out redox chemistry with the species to be oxidized on the P-side while reduction occurs on the N-side of the membrane. Terminal redox cofactors are used which have pKas that are strongly dependent on their redox state, so that protons are lost on oxidation and gained on reduction. Bacteriorhodopsin is a true proton pump. Light activation triggers trans to cis isomerization of a bound retinal. Strong electrostatic interactions within clusters of amino acids are modified by the conformational changes initiated by retinal motion leading to changes in proton affinity, driving transmembrane proton transfer. Cytochrome c oxidase (CcO) catalyzes the reduction of O2 to water. The protons needed for chemistry are bound from the N-side. The reduction chemistry also drives proton pumping from N- to P-side. Overall, in CcO the uptake of 4 electrons to reduce O2 transports 8 charges across the membrane, with each reduction fully coupled to removal of two protons from the N-side, the delivery of one for chemistry and transport of the other to the P-side.

Deprotonation of D96 in bacteriorhodopsin opens the proton uptake pathway

Despite extensive investigation, the precise mechanism controlling the opening of the cytoplasmic proton uptake pathway in bacteriorhodopsin (bR) has remained a mystery. From an analysis of the X-ray structure of the D96G/F171C/F219L triple mutant of bR and 60 independent molecular dynamics simulations of bR photointermediates, we report that the deprotonation of D96, a key residue in proton transfer reactions, serves two roles that occur sequentially. First, D96 donates a proton to the Schiff base. Subsequently, the deprotonation of D96 serves to "unlatch" the cytoplasmic side. The latching function of D96 appears to be remarkably robust, functioning to open hydration channels in all photointermediate structures. These results suggest that the protonation state of D96 may be the critical biophysical cue controlling the opening and closing of the cytoplasmic half-channel in bR. We suspect that this protonation-switch mechanism could also be utilized in other proton pumps to minimize backflow and reinforce directionality.

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.

Photocycle of Exiguobacterium sibiricum rhodopsin characterized by low-temperature trapping in the IR and time-resolved studies in the visible

The photocycle of the retinal protein from Exiguobacterium sibiricum, which differs from bacteriorhodopsin in both its primary donor and acceptor, is characterized by visible and infrared spectroscopy. At pH above pKa ~6.5, we find a bacteriorhodopsin-like photocycle, which originates from excitation of the all-trans retinal chromophore with K-, L-, M-, and N-like intermediates. At pH below pKa ~6.5, the M state, which reflects Schiff base deprotonation during proton pumping, is not accumulated. However, using the infrared band at ~1760 cm(-1) as a marker for transient protonation of the primary acceptor, we find that Schiff base deprotonation must have occurred at pH not only above but also below the pKa ~6.5. Thus, the M state is formed but not accumulated for kinetic reasons. Further, chromophore reisomerization from the 13-cis to the all-trans conformation occurs very late in the photocycle. The strongly red-shifted states that dominate the second half of the cycle are produced before the reisomerization step, and by this criterion, they are not O-like but rather N-like states. The assignment of photocycle intermediates enables reevaluation of the photocycle; its specific features are discussed in relation to the general mechanism of proton transport in retinal proteins.

Breaking the carboxyl rule: lysine 96 facilitates reprotonation of the Schiff base in the photocycle of a retinal protein from Exiguobacterium sibiricum

A lysine instead of the usual carboxyl group is in place of the internal proton donor to the retinal Schiff base in the light-driven proton pump of Exiguobacterium sibiricum (ESR). The involvement of this lysine in proton transfer is indicated by the finding that its substitution with alanine or other residues slows reprotonation of the Schiff base (decay of the M intermediate) by more than 2 orders of magnitude. In these mutants, the rate constant of the M decay linearly decreases with a decrease in proton concentration, as expected if reprotonation is limited by the uptake of a proton from the bulk. In wild type ESR, M decay is biphasic, and the rate constants are nearly pH-independent between pH 6 and 9. Proton uptake occurs after M formation but before M decay, which is especially evident in D2O and at high pH. Proton uptake is biphasic; the amplitude of the fast phase decreases with a pKa of 8.5 ± 0.3, which reflects the pKa of the donor during proton uptake. Similarly, the fraction of the faster component of M decay decreases and the slower one increases, with a pKa of 8.1 ± 0.2. The data therefore suggest that the reprotonation of the Schiff base in ESR is preceded by transient protonation of an initially unprotonated donor, which is probably the ε-amino group of Lys-96 or a water molecule in its vicinity, and it facilitates proton delivery from the bulk to the reaction center of the protein.

Effect of cosolvents on nano-confined water: a molecular dynamics study

We present results from atomistic molecular dynamics simulations to characterize the effects of cosolvents, such as urea and guanidinium (Gdm) salts, on the water confined in hydrophobic carbon nanotubes. We observed complete drying of the nanotube interiors of diameter ranging from 8 to 17 Å in urea. In contrast, the water population within nanotube cores smaller than 12 Å remains unaffected in GdmCl solution, whereas larger nanotube interiors become partially dehydrated with prevailing presence of stable Gdm(+)-Gdm(+) dimers. The molecular arrangement and the lifetime inside the nanotube were found to be characteristics of a particular cosolvent. In both urea and GdmCl solutions, preferential cosolvent intrusion resulting in nanotube dehydration is driven by the stronger dispersion interaction of cosolvent than water with the nanotube. The partial drying of the hydrophobic core is attributed to guanidinium's better hydration and weaker self-association propensity compared to urea, as well as to its moderate ion-pairing with strongly hydrated chloride ions. The Gdm(+) induced dehydration varies with the charge density of counter-ions, as the presence of high charge-density sulfate ions impedes penetration of guanidinium, and consequent dehydration of the nanotube. These findings provide important insights into the effect of cosolvents on the nano-confined water in a hydrophobic environment.

Conformational changes in the archaerhodopsin-3 proton pump: detection of conserved strongly hydrogen bonded water networks

Archaerhodopsin-3 (AR3) is a light-driven proton pump from Halorubrum sodomense, but little is known about its photocycle. Recent interest has focused on AR3 because of its ability to serve both as a high-performance, genetically-targetable optical silencer of neuronal activity and as a membrane voltage sensor. We examined light-activated structural changes of the protein, retinal chromophore, and internal water molecules during the photocycle of AR3. Low-temperature and rapid-scan time-resolved FTIR-difference spectroscopy revealed that conformational changes during formation of the K, M, and N photocycle intermediates are similar, although not identical, to bacteriorhodopsin (BR). Positive/negative bands in the region above 3,600 cm( - 1), which have previously been assigned to structural changes of weakly hydrogen bonded internal water molecules, were substantially different between AR3 and BR. This included the absence of positive bands recently associated with a chain of proton transporting water molecules in the cytoplasmic channel and a weakly hydrogen bonded water (W401), which is part of a hydrogen-bonded pentagonal cluster located near the retinal Schiff base. However, many of the broad IR continuum absorption changes below 3,000 cm( - 1) assigned to networks of water molecules involved in proton transport through cytoplasmic and extracellular portions in BR were very similar in AR3. This work and subsequent studies comparing BR and AR3 structural changes will help identify conserved elements in BR-like proton pumps as well as bioengineer AR3 to optimize neural silencing and voltage sensing.

Proton transfer via a transient linear water-molecule chain in a membrane protein

High-resolution protein ground-state structures of proton pumps and channels have revealed internal protein-bound water molecules. Their possible active involvement in protein function has recently come into focus. An illustration of the formation of a protonated protein-bound water cluster that is actively involved in proton transfer was described for the membrane protein bacteriorhodopsin (bR) [Garczarek F, Gerwert K (2006) Nature 439:109-112]. Here we show through a combination of time-resolved FTIR spectroscopy and molecular dynamics simulations that three protein-bound water molecules are rearranged by a protein conformational change that resulted in a transient Grotthuss-type proton-transfer chain extending through a hydrophobic protein region of bR. This transient linear water chain facilitates proton transfer at an intermediate conformation only, thereby directing proton transfer within the protein. The rearrangement of protein-bound water molecules that we describe, from inactive positions in the ground state to an active chain in an intermediate state, appears to be energetically favored relative to transient incorporation of water molecules from the bulk. Our discovery provides insight into proton-transfer mechanisms through hydrophobic core regions of ubiquitous membrane spanning proteins such as G-protein coupled receptors or cytochrome C oxidases.