Bacteriorhodopsin photocycle at cryogenic temperatures reveals distributed barriers of conformational substates

Four-dimensional microED of conformational dynamics in protein microcrystals on the femto-to-microsecond timescales

As structural determination of protein complexes approaches atomic resolution, there is an increasing focus on conformational dynamics. Here we conceptualize the combination of two techniques which have become established in recent years: microcrystal electron diffraction and ultrafast electron microscopy. We show that the extremely low dose of pulsed photoemission still enables microED due to the strength of the electron bunching from diffraction of the protein crystals. Indeed, ultrafast electron diffraction experiments on protein crystals have already been demonstrated to be effective in measuring intermolecular forces in protein microcrystals. We discuss difficulties that may arise in the acquisition and processing of data and the overall feasibility of the experiment, paying specific attention to dose and signal-to-noise ratio. In doing so, we outline a detailed workflow that may be effective in minimizing the dose on the specimen. A series of model systems that would be good candidates for initial experiments is provided.

Modulation of Light Energy Transfer from Chromophore to Protein in the Channelrhodopsin ReaChR

The function of photoreceptors relies on efficient transfer of absorbed light energy from the chromophore to the protein to drive conformational changes that ultimately generate an output signal. In retinal-binding proteins, mainly two mechanisms exist to store the photon energy after photoisomerization: 1) conformational distortion of the prosthetic group retinal, and 2) charge separation between the protonated retinal Schiff base (RSBH⁺) and its counterion complex. Accordingly, energy transfer to the protein is achieved by chromophore relaxation and/or reduction of the charge separation in the RSBH⁺-counterion complex. Combining FTIR and UV-Vis spectroscopy along with molecular dynamics simulations, we show here for the widely used, red-activatable Volvox carteri channelrhodopsin-1 derivate ReaChR that energy storage and transfer into the protein depends on the protonation state of glutamic acid E163 (Ci1), one of the counterions of the RSBH⁺. Ci1 retains a pK_a of 7.6 so that both its protonated and deprotonated forms equilibrate at physiological conditions. Protonation of Ci1 leads to a rigid hydrogen-bonding network in the active-site region. This stabilizes the distorted conformation of the retinal after photoactivation and decelerates energy transfer into the protein by impairing the release of the strain energy. In contrast, with deprotonated Ci1 or removal of the Ci1 glutamate side chain, the hydrogen-bonded system is less rigid, and energy transfer by chromophore relaxation is accelerated. Based on the hydrogen out-of-plane (HOOP) band decay kinetics, we determined the activation energy for these processes in dependence of the Ci1 protonation state.

Infrared Difference Spectroscopy of Proteins: From Bands to Bonds

Infrared difference spectroscopy probes vibrational changes of proteins upon their perturbation. Compared with other spectroscopic methods, it stands out by its sensitivity to the protonation state, H-bonding, and the conformation of different groups in proteins, including the peptide backbone, amino acid side chains, internal water molecules, or cofactors. In particular, the detection of protonation and H-bonding changes in a time-resolved manner, not easily obtained by other techniques, is one of the most successful applications of IR difference spectroscopy. The present review deals with the use of perturbations designed to specifically change the protein between two (or more) functionally relevant states, a strategy often referred to as reaction-induced IR difference spectroscopy. In the first half of this contribution, I review the technique of reaction-induced IR difference spectroscopy of proteins, with special emphasis given to the preparation of suitable samples and their characterization, strategies for the perturbation of proteins, and methodologies for time-resolved measurements (from nanoseconds to minutes). The second half of this contribution focuses on the spectral interpretation. It starts by reviewing how changes in H-bonding, medium polarity, and vibrational coupling affect vibrational frequencies, intensities, and bandwidths. It is followed by band assignments, a crucial aspect mostly performed with the help of isotopic labeling and site-directed mutagenesis, and complemented by integration and interpretation of the results in the context of the studied protein, an aspect increasingly supported by spectral calculations. Selected examples from the literature, predominately but not exclusively from retinal proteins, are used to illustrate the topics covered in this review.

Water-Polymer Coupling Induces a Dynamical Transition in Microgels

The long debated protein dynamical transition was recently found also in nonbiological macromolecules, such as poly- N-isopropylacrylamide (PNIPAM) microgels. Here, by using atomistic molecular dynamics simulations, we report a description of the molecular origin of the dynamical transition in these systems. We show that PNIPAM and water dynamics below the dynamical transition temperature T _d are dominated by methyl group rotations and hydrogen bonding, respectively. By comparing with bulk water, we unambiguously identify PNIPAM-water hydrogen bonding as mainly responsible for the occurrence of the transition. The observed phenomenology thus crucially depends on the water-macromolecule coupling, being relevant to a wide class of hydrated systems, independently from the biological function.

Far infrared spectroscopy of hydrogen bonding collective motions in complex molecular systems

Far infrared spectroscopy as a tool for the study of inter and intramolecular interactions in complex molecular structures.

Light-induced structural changes during early photo-intermediates of the eubacterial Cl−pump Fulvimarina rhodopsin observed by FTIR difference spectroscopy

Fulvimarina pelagirhodopsin (FR) is a member of inward eubacterial light-activated Cl⁻translocating rhodopsins (ClR) that were found recently in marine bacteria.

Light-driven Na(+) pump from Gillisia limnaea: a high-affinity Na(+) binding site is formed transiently in the photocycle

A group of microbial retinal proteins most closely related to the proton pump xanthorhodopsin has a novel sequence motif and a novel function. Instead of, or in addition to, proton transport, they perform light-driven sodium ion transport, as reported for one representative of this group (KR2) from Krokinobacter. In this paper, we examine a similar protein, GLR from Gillisia limnaea, expressed in Escherichia coli, which shares some properties with KR2 but transports only Na(+). The absorption spectrum of GLR is insensitive to Na(+) at concentrations of ≤3 M. However, very low concentrations of Na(+) cause profound differences in the decay and rise time of photocycle intermediates, consistent with a switch from a "Na(+)-independent" to a "Na(+)-dependent" photocycle (or photocycle branch) at ∼60 μM Na(+). The rates of photocycle steps in the latter, but not the former, are linearly dependent on Na(+) concentration. This suggests that a high-affinity Na(+) binding site is created transiently after photoexcitation, and entry of Na(+) from the bulk to this site redirects the course of events in the remainder of the cycle. A greater concentration of Na(+) is needed for switching the reaction path at lower pH. The data suggest therefore competition between H(+) and Na(+) to determine the two alternative pathways. The idea that a Na(+) binding site can be created at the Schiff base counterion is supported by the finding that upon perturbation of this region in the D251E mutant, Na(+) binds without photoexcitation. Binding of Na(+) to the mutant shifts the chromophore maximum to the red like that of H(+), which occurs in the photocycle of the wild type.

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.

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.

Probing specific molecular processes and intermediates by time-resolved Fourier transform infrared spectroscopy: application to the bacteriorhodopsin photocycle

We present a general approach for probing the kinetics of specific molecular processes in proteins by time-resolved Fourier transform infrared (IR) spectroscopy. Using bacteriorhodopsin (bR) as a model we demonstrate that by appropriately monitoring some selected IR bands it is possible obtaining the kinetics of the most important events occurring in the photocycle, namely changes in the chromophore and the protein backbone conformation, and changes in the protonation state of the key residues implicated in the proton transfers. Besides confirming widely accepted views of the bR photocycle, our analysis also sheds light into some disputed issues: the degree of retinal torsion in the L intermediate to respect the ground state; the possibility of a proton transfer from Asp85 to Asp212; the relationship between the protonation/deprotonation of Asp85 and the proton release complex; and the timing of the protein backbone dynamics. By providing a direct way to estimate the kinetics of photocycle intermediates the present approach opens new prospects for a robust quantitative kinetic analysis of the bR photocycle, which could also benefit the study of other proteins involved in photosynthesis, in phototaxis, or in respiratory chains.

Low-temperature FTIR study of multiple K intermediates in the photocycles of bacteriorhodopsin and xanthorhodopsin

Low-temperature FTIR spectroscopy of bacteriorhodopsin and xanthorhodopsin was used to elucidate the number of K-like bathochromic states, their sequence, and their contributions to the photoequilibrium mixtures created by illumination at 80-180 K. We conclude that in bacteriorhodopsin the photocycle includes three distinct K-like states in the sequence bR (hv)--> I* --> J --> K(0) --> K(E) --> L --> ..., and similarly in xanthorhodopsin. K(0) is the main fraction in the mixture at 77 K that is formed from J. K(0) becomes thermally unstable above approximately 50 K in both proteins. At 77 K, both J-to-K(0) and K(0)-to-K(E) transitions occur and, contrarily to long-standing belief, cryogenic trapping at 77 K does not produce a pure K state but a mixture of the two states, K(0) and K(E), with contributions from K(E) of approximately 15 and approximately 10% in the two retinal proteins, respectively. Raising the temperature leads to increasing conversion of K(0) to K(E), and the two states coexist (without contamination from non-K-like states) in the 80-140 K range in bacteriorhodopsin, and in the 80-190 K range in xanthorhodopsin. Temperature perturbation experiments in these regions of coexistence revealed that, in spite of the observation of apparently stable mixtures of K(0) and K(E), the two states are not in thermally controlled equilibrium. The K(0)-to-K(E) transition is unidirectional, and the partial transformation to K(E) is due to distributed kinetics, which governs the photocycle dynamics at temperatures below approximately 245 K (Dioumaev and Lanyi, Biochemistry 2008, 47, 11125-11133). From spectral deconvolution, we conclude that the K(E) state, which is increasingly present at higher temperatures, is the same intermediate that is detected by time-resolved FTIR prior to its decay, on a time scale of hundreds of nanoseconds at ambient temperature (Dioumaev and Braiman, J. Phys. Chem. B 1997, 101, 1655-1662), into the K(L) state. We were unable to trap the latter separately from K(E) at low temperature, due to the slow distributed kinetics and the increasingly faster overlapping formation of the L state. Formation of the two consecutive K-like states in both proteins is accompanied by distortion of two different weakly bound water molecules: one in K(0), the other in K(E). The first, well-documented in bacteriorhodopsin at 77 K where K(0) dominates, was assigned to water 401 in bacteriorhodopsin. The other water molecule, whose participation has not been described previously, is disturbed on the next step of the photocycle, in K(E), in both proteins. In bacteriorhodopsin, the most likely candidate is water 407. However, unlike bacteriorhodopsin, the crystal structure of xanthorhodopsin lacks homologous weakly bound water molecules.

Photoreactions and structural changes of anabaena sensory rhodopsin

Anabaena sensory rhodopsin (ASR) is an archaeal-type rhodopsin found in eubacteria. The gene encoding ASR forms a single operon with ASRT (ASR transducer) which is a 14 kDa soluble protein, suggesting that ASR functions as a photochromic sensor by activating the soluble transducer. This article reviews the detailed photoreaction processes of ASR, which were studied by low-temperature Fourier-transform infrared (FTIR) and UV-visible spectroscopy. The former research reveals that the retinal isomerization is similar to bacteriorhodopsin (BR), but the hydrogen-bonding network around the Schiff base and cytoplasmic region is different. The latter study shows the stable photoproduct of the all-trans form is 100% 13-cis, and that of the 13-cis form is 100% all-trans. These results suggest that the structural changes of ASR in the cytoplasmic domain play important roles in the activation of the transducer protein, and photochromic reaction is optimized for its sensor function.

Functional and shunt states of bacteriorhodopsin resolved by 250 GHz dynamic nuclear polarization-enhanced solid-state NMR

Observation and structural studies of reaction intermediates of proteins are challenging because of the mixtures of states usually present at low concentrations. Here, we use a 250 GHz gyrotron (cyclotron resonance maser) and cryogenic temperatures to perform high-frequency dynamic nuclear polarization (DNP) NMR experiments that enhance sensitivity in magic-angle spinning NMR spectra of cryo-trapped photocycle intermediates of bacteriorhodopsin (bR) by a factor of approximately 90. Multidimensional spectroscopy of U-(13)C,(15)N-labeled samples resolved coexisting states and allowed chemical shift assignments in the retinylidene chromophore for several intermediates not observed previously. The correlation spectra reveal unexpected heterogeneity in dark-adapted bR, distortion in the K state, and, most importantly, 4 discrete L substates. Thermal relaxation of the mixture of L's showed that 3 of these substates revert to bR(568) and that only the 1 substate with both the strongest counterion and a fully relaxed 13-cis bond is functional. These definitive observations of functional and shunt states in the bR photocycle provide a preview of the mechanistic insights that will be accessible in membrane proteins via sensitivity-enhanced DNP NMR. These observations would have not been possible absent the signal enhancement available from DNP.

Infrared monitoring of interlayer water in stacks of purple membranes

The thermodynamic behavior of films of hydrated purple membranes from Halobacterium salinarum and the water confined in it was studied by Fourier transform infrared spectroscopy in the 180-280 K range. Unlike bulk water, water in the thin layers sandwiched between the biological membranes does not freeze at 273 K but will be supercooled to approximately 256 K. The melting point is unaffected, leading to hysteresis between 250 and 273 K. In its heating branch, a gradually increasing light-scattering by ice is observed with rate-limiting kinetics of tens of minutes. Infrared (IR) spectra decomposition provided extinction coefficients for the confined water vibrational bands and their changes upon freezing. Because of the hysteresis, at any given temperature in the 255-270 K range, the interbilayer water could be either liquid or frozen, depending on thermal history. We find that this difference affects the dynamics of the bacteriorhodopsin photocycle in the hysteresis range: the decay of the M and N states and the redistribution between them are different depending on whether or not the water was initially precooled to below the freezing point. However, freezing of interbilayer water does block the M to N transition. Unlike the water, the purple membrane lipids do not undergo any IR-detectable phase transition in the 180-280 K range.

High resolution approach to the native state ensemble kinetics and thermodynamics

Many biologically interesting functions such as allosteric switching or protein-ligand binding are determined by the kinetics and mechanisms of transitions between various conformational substates of the native basin of globular proteins. To advance our understanding of these processes, we constructed a two-dimensional free energy surface (FES) of the native basin of a small globular protein, Trp-cage. The corresponding order parameters were defined using two native substructures of Trp-cage. These calculations were based on extensive explicit water all-atom molecular dynamics simulations. Using the obtained two-dimensional FES, we studied the transition kinetics between two Trp-cage conformations, finding that switching process shows a borderline behavior between diffusive and weakly-activated dynamics. The transition is well-characterized kinetically as a biexponential process. We also introduced a new one-dimensional reaction coordinate for the conformational transition, finding reasonable qualitative agreement with the two-dimensional kinetics results. We investigated the distribution of all the 38 native nuclear magnetic resonance structures on the obtained FES, analyzing interactions that stabilize specific low-energy conformations. Finally, we constructed a FES for the same system but with simple dielectric model of water instead of explicit water, finding that the results were surprisingly similar in a small region centered on the native conformations. The dissimilarities between the explicit and implicit model on the larger-scale point to the important role of water in mediating interactions between amino acid residues.

Switch from conventional to distributed kinetics in the bacteriorhodopsin photocycle

Below 195 K, the bacteriorhodopsin photocycle could not be adequately described with exponential kinetics [Dioumaev, A. K., and Lanyi, J. K. (2007) Proc. Natl. Acad. Sci. U.S.A. 104, 9621-9626] but required distributed kinetics, previously found in hemoglobin and myoglobin at temperatures below the vitrification point of the surrounding solvent. The aim of this study is to determine which factors cause the switch from this low-temperature regime to the conventional kinetics observed at ambient temperature. The photocycle was monitored by time-resolved FTIR between 180 and 280 K, using the D96N mutant. Depending on the temperature, decay and temporal redistribution of two or three intermediates (L, M, and N) were observed. Above approximately 245 K, an abrupt change in the kinetic behavior of the photocycle takes place. It does not affect the intermediates present but greatly accelerates their decay. Below approximately 240 K, a kinetic pattern with partial decay that cannot be explained by conventional kinetics, but suggesting distributed kinetics, was dominant, while above approximately 250 K, there were no significant deviations from exponential behavior. The approximately 245 K critical point is >/=10 K below the freezing point of interbilayer water, and we were unable to correlate it with any FTIR-detectable transition of the lipids. Therefore, we attribute the change from distributed to conventional kinetics to a thermodynamic phase transition in the protein. Most probably, it is related to the freezing and thawing of internal fluctuations of the protein, known as the dynamic phase transition, although in bacteriorhodopsin the latter is usually believed to take place at least 15 K below the observed critical temperature of approximately 245 K.

Observations on rate theory for rugged energy landscapes

The potential energy profile for many complex reactions of proteins, such as folding or allosteric conformational change, involves many different scales of molecular motion along the reaction coordinate. Although it is natural to model the dynamics of motion along such rugged energy landscapes as diffusional (the Smoluchowski equation; SE), problems arise because the frictional forces generated by the molecular surround are typically not strong enough to justify the use of the SE. Here, we discuss the fundamental theory behind the SE and note that it may be justified through a master equation when reduced to its continuum limit. However, the SE cannot be used for rough energy landscapes, where the continuum limit is ill defined. Instead, we suggest that one should use a mean first passage time expression derived from a master equation, and show how this approach can be used to glean information about the underlying dynamics of barrier crossing. We note that the potential profile in the SE is that of the microbarriers between conformational substates, and that there is a temperature-dependent, effective friction associated with the long residence time in the microwells that populate the rough landscape. The number of recrossings of the overall barrier is temperature-dependent, governed by the microbarriers and not by the effective friction. We derive an explicit expression for the mean number of recrossings and its temperature dependence. Finally, we note that the mean first passage time can be used as a departure point for measuring the roughness of the landscape.

Energy transformations early in the bacteriorhodopsin photocycle revealed by DNP-enhanced solid-state NMR

By exploiting dynamic nuclear polarization (DNP) at 90 K, we observe the first NMR spectrum of the K intermediate in the ion-motive photocycle of bacteriorhodopsin. The intermediate is identified by its reversion to the resting state of the protein in red light and by its thermal decay to the L intermediate. The (15)N chemical shift of the Schiff base in K indicates that contact has been lost with its counterion. Under these circumstances, the visible absorption of K is expected to be more red-shifted than is observed and this suggests torsion around single bonds of the retinylidene chromophore. This is in contrast to the development of a strong counterion interaction and double bond torsion in L. Thus, photon energy is stored in electrostatic modes in K and is transferred to torsional modes in L. This transfer is facilitated by the reduction in bond alternation that occurs with the initial loss of the counterion interaction, and is driven by the attraction of the Schiff base to a new counterion. Nevertheless, the process appears to be difficult, as judged by the multiple L substates, with weaker counterion interactions, that are trapped at lower temperatures. The double-bond torsion ultimately developed in the first half of the photocycle is probably responsible for enforcing vectoriality in the pump by causing a decisive switch in the connectivity of the active site once the Schiff base and its counterion are neutralized by proton transfer.