Search

Photoreceptors containing the light-oxygen-voltage (LOV) domain elicit biological responses upon excitation of their flavin mononucleotide (FMN) chromophore by blue light. The mechanism and kinetics of dark-state recovery are not well understood. Here we incorporated the non-canonical amino acid p-cyanophenylalanine (CNF) by genetic code expansion technology at 45 positions of the bacterial transcription factor EL222. Screening of light-induced changes in infrared (IR) absorption frequency, electric field and hydration of the nitrile groups identified residues CNF31 and CNF35 as reporters of monomer/oligomer and caged/decaged equilibria, respectively. Time-resolved multi-probe UV/visible and IR spectroscopy experiments of the lit-to-dark transition revealed four dynamical events. Predominantly, rearrangements around the A'α helix interface (CNF31 and CNF35) precede FMN-cysteinyl adduct scission, folding of α-helices (amide bands), and relaxation of residue CNF151. This study illustrates the importance of characterizing all parts of a protein and suggests a key role for the N-terminal A'α extension of the LOV domain in controlling EL222 photocycle length., (© 2023 The Authors. Protein Science published by Wiley Periodicals LLC on behalf of The Protein Society.)

The catalytic activity of proteins is a function of structural changes. Very often these are as minute as protonation changes, hydrogen bonding changes, and amino acid side chain reorientations. To resolve these, a methodology is afforded that not only provides the molecular sensitivity but allows for tracing the sequence of these hierarchical reactions at the same time. This feature article showcases results from time-resolved IR spectroscopy on channelrhodopsin (ChR), light-oxygen-voltage (LOV) domain protein, and cryptochrome (CRY). All three proteins are activated by blue light, but their biological role is drastically different. Channelrhodopsin is a transmembrane retinylidene protein which represents the first light-activated ion channel of its kind and which is involved in primitive vision (phototaxis) of algae. LOV and CRY are flavin-binding proteins acting as photoreceptors in a variety of signal transduction mechanisms in all kingdoms of life. Beyond their biological relevance, these proteins are employed in exciting optogenetic applications. We show here how IR difference absorption resolves crucial structural changes of the protein after photonic activation of the chromophore. Time-resolved techniques are introduced that cover the time range from nanoseconds to minutes along with some technical considerations. Finally, we provide an outlook toward novel experimental approaches that are currently developed in our laboratories or are just in our minds ("Gedankenexperimente"). We believe that some of them have the potential to provide new science.

Sensory rhodopsin II (SRII) is the primary light sensor in the photophobic reaction of the halobacterium Natronomonas pharaonis. Photoactivation of SRII results in a movement of helices F and G of this seven-helical transmembrane protein. This conformational change is conveyed to the transducer protein (HtrII). Global changes in the protein backbone have been monitored by IR difference spectroscopy by recording frequency shifts in the amide bands. Here we investigate local structural changes by judiciously inserting thiocyanides at different locations of SRII. These vibrational Stark probes absorb in a frequency range devoid of any protein vibrations and respond to local changes in the dielectric, electrostatics, and hydrogen bonding. As a proof of principle, we demonstrate the use of Stark probes to test the conformational changes occurring in SRII 12 ms after photoexcitation and later. Thus, a methodology is provided to trace local conformational changes in membrane proteins by a minimal invasive probe at the high temporal resolution inherent to IR spectroscopy.

The discovery of channelrhodopsins introduced a new class of light-gated ion channels, which when genetically encoded in host cells resulted in the development of optogenetics. Channelrhodopsin-2 from Chlamydomonas reinhardtii, CrChR2, is the most widely used optogenetic tool in neuroscience. To explore the connection between the gating mechanism and the influx and efflux of water molecules in CrChR2, we have integrated light-induced time-resolved infrared spectroscopy and electrophysiology. Cross-correlation analysis revealed that ion conductance tallies with peptide backbone amide I vibrational changes at 1,665(-) and 1,648(+) cm(-1). These two bands report on the hydration of transmembrane α-helices as concluded from vibrational coupling experiments. Lifetime distribution analysis shows that water influx proceeded in two temporally separated steps with time constants of 10 μs (30%) and 200 μs (70%), the latter phase concurrent with the start of ion conductance. Water efflux and the cessation of the ion conductance are synchronized as well, with a time constant of 10 ms. The temporal correlation between ion conductance and hydration of helices holds for fast (E123T) and slow (D156E) variants of CrChR2, strengthening its functional significance.

Channelrhodopsins (ChRs) are light-gated cation channels. After blue-light excitation, the protein undergoes a photocycle with different intermediates. Here, we have recorded transient absorbance changes of ChR2 from Chlamydomonas reinhardtii in the visible and infrared regions with nanosecond time resolution, the latter being accomplished using tunable quantum cascade lasers. Because proton transfer reactions play a key role in channel gating, we determined vibrational as well as kinetic isotope effects (VIEs and KIEs) of carboxylic groups of various key aspartic and glutamic acid residues by monitoring their C=O stretching vibrations in H2O and in D2O. D156 exhibits a substantial KIE (>2) in its deprotonation and reprotonation, which substantiates its role as the internal proton donor to the retinal Schiff base. The unusual VIE of D156, upshifted from 1736 cm(-1) to 1738 cm(-1) in D2O, was scrutinized by studying the D156E variant. The C=O stretch of E156 shifted down by 8 cm(-1) in D2O, providing evidence for the accessibility of the carboxylic group. The C=O stretching band of E90 exhibits a VIE of 9 cm(-1) and a KIE of ∼2 for the de- and the reprotonation reactions during the lifetime of the late desensitized state. The KIE of 1 determined in the time range from 20 ns to 5 ms is incompatible with early deprotonation of E90., (Copyright © 2015 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

We examine the role of Lys-377, the only charged residue in helix XI, on the functional mechanism of the Na(+)-sugar melibiose symporter from Escherichia coli. Intrinsic fluorescence, FRET, and Fourier transform infrared difference spectroscopy reveal that replacement of Lys-377 with either Cys, Val, Arg, or Asp disables both Na(+) and melibiose binding. On the other hand, molecular dynamics simulations extending up to 200-330 ns reveal that Lys-377 (helix XI) interacts with the anionic side chains of two of the three putative ligands for cation binding (Asp-55 and Asp-59 in helix II). When Asp-59 is protonated during the simulations, Lys-377 preferentially interacts with Asp-55. Interestingly, when a Na(+) ion is positioned in the Asp-55-Asp-59 environment, Asp-124 in helix IV (a residue essential for melibiose binding) reorients and approximates the Asp-55-Asp-59 pair, and all three acidic side chains act as Na(+) ligands. Under these conditions, the side chain of Lys-377 interacts with the carboxylic moiety of these three Asp residues. These data highlight the crucial role of the Lys-377 residue in the spatial organization of the Na(+) binding site. Finally, the analysis of the second-site revertants of K377C reveals that mutation of Ile-22 (in helix I) preserves Na(+) binding, whereas that of melibiose is largely abolished according to spectroscopic measurements. This amino acid is located in the border of the sugar-binding site and might participate in sugar binding through apolar interactions., (© 2015 by The American Society for Biochemistry and Molecular Biology, Inc.)

Light-gated ion permeation by channelrhodopsin-2 (ChR2) relies on the photoisomerization of the retinal chromophore and the subsequent photocycle, leading to the formation (on-gating) and decay (off-gating) of the conductive state. Here, we have analyzed the photocycle of a fast-cycling ChR2 variant (E123T mutation, also known as ChETA), by time-resolved UV/vis, step-scan FT-IR, and tunable quantum cascade laser IR spectroscopies with nanosecond resolution. Pre-gating conformational changes rise with a half-life of 200 ns, silent to UV/vis but detected by IR spectroscopy. They involve changes in the peptide backbone and in the H-bond of the side chain of the critical residue D156. Thus, the P1(500) intermediate must be separated into early and late states. Light-adapted ChR2 contains a mixture of all-trans and 13-cis retinal in a 70:30 ratio which are both photoactive. Analysis of ethylenic and fingerprint vibrations of retinal provides evidence that the 13-cis photocycle recovers in 1 ms. This recovery is faster than channel off-gating and most of the proton transfer reactions, implying that the 13-cis photocycle is of minor functional relevance for ChR2.

Water plays an essential role in the structure and function of proteins, particularly in the less understood class of membrane proteins. As the first of its kind, channelrhodopsin is a light-gated cation channel and paved the way for the new and vibrant field of optogenetics, where nerve cells are activated by light. Still, the molecular mechanism of channelrhodopsin is not understood. Here, we applied time-resolved FT-IR difference spectroscopy to channelrhodopsin-1 from Chlamydomonas augustae. It is shown that the (conductive) P2(380) intermediate decays with τ ≈ 40 ms and 200 ms after pulsed excitation. The vibrational changes between the closed and the conductive states were analyzed in the X-H stretching region (X = O, S, N), comprising vibrational changes of water molecules, sulfhydryl groups of cysteine side chains and changes of the amide A of the protein backbone. The O-H stretching vibrations of "dangling" water molecules were detected in two different states of the protein using H2 (18)O exchange. Uncoupling experiments with a 1:1 mixture of H2O:D2O provided the natural uncoupled frequencies of the four O-H (and O-D) stretches of these water molecules, each with a very weakly hydrogen-bonded O-H group (3639 and 3628 cm(-1)) and with the other O-H group medium (3440 cm(-1)) to moderately strongly (3300 cm(-1)) hydrogen-bonded. Changes in amide A and thiol vibrations report on global and local changes, respectively, associated with the formation of the conductive state. Future studies will aim at assigning the respective cysteine group(s) and at localizing the "dangling" water molecules within the protein, providing a better understanding of their functional relevance in CaChR1.

Changes in vibrational circular dichroism (VCD) were recorded on-line during a chemical reaction. The chiral complex nickel-(-)-sparteine chloride was hydrolyzed to free (-)-sparteine base in a biphasic system of sodium hydroxide solution and chloroform (CHCl(3)). Infrared (IR) and VCD spectra were iteratively recorded after pumping a sample from the CHCl(3) phase through a lab-built VCD spectrometer equipped with a tunable mid-IR quantum cascade laser light source, which allows for VCD measurements even in the presence of strongly absorbing backgrounds. Time-dependent VCD spectra were analyzed by singular value decomposition and global exponential fitting. Spectral features corresponding to the complex and free (-)-sparteine could be clearly identified in the fitted amplitude spectrum, which was associated with an exponential decay with an apparent time constant of 127 min (t(½) = 88 min)., (© 2014 Wiley Periodicals, Inc.)

Channelrhodopsin-1 from Chlamydomonas augustae (CaChR1) is a light-activated cation channel, which is a promising optogenetic tool. We show by resonance Raman spectroscopy and retinal extraction followed by high pressure liquid chromatography (HPLC) that the isomeric ratio of all-trans to 13-cis of solubilized channelrhodopsin-1 is with 70:30 identical to channelrhodopsin-2 from Chlamydomonas reinhardtii (CrChR2). Critical frequency shifts in the retinal vibrations are identified in the Raman spectrum upon transition to the open (conductive P2(380)) state. Fourier transform infrared spectroscopy (FTIR) spectra indicate different structures of the open states in the two channelrhodopsins as reflected by the amide I bands and the protonation pattern of acidic amino acids., (Copyright © 2014 The Authors. Published by Elsevier B.V. All rights reserved.)

Monitoring the dynamics of protonation and protein backbone conformation changes during the function of a protein is an essential step towards understanding its mechanism. Protonation and conformational changes affect the vibration pattern of amino acid side chains and of the peptide bond, respectively, both of which can be probed by infrared (IR) difference spectroscopy. For proteins whose function can be repetitively and reproducibly triggered by light, it is possible to obtain infrared difference spectra with (sub)microsecond resolution over a broad spectral range using the step-scan Fourier transform infrared technique. With -10(2)-10(3) repetitions of the photoreaction, the minimum number to complete a scan at reasonable spectral resolution and bandwidth, the noise level in the absorption difference spectra can be as low as -10(-) (4), sufficient to follow the kinetics of protonation changes from a single amino acid. Lower noise levels can be accomplished by more data averaging and/or mathematical processing. The amount of protein required for optimal results is between 5-100 µg, depending on the sampling technique used. Regarding additional requirements, the protein needs to be first concentrated in a low ionic strength buffer and then dried to form a film. The protein film is hydrated prior to the experiment, either with little droplets of water or under controlled atmospheric humidity. The attained hydration level (g of water / g of protein) is gauged from an IR absorption spectrum. To showcase the technique, we studied the photocycle of the light-driven proton-pump bacteriorhodopsin in its native purple membrane environment, and of the light-gated ion channel channelrhodopsin-2 solubilized in detergent.

The new and vibrant field of optogenetics was founded by the seminal discovery of channelrhodopsin, the first light-gated cation channel. Despite the numerous applications that have revolutionised neurophysiology, the functional mechanism is far from understood on the molecular level. An arsenal of biophysical techniques has been established in the last decades of research on microbial rhodopsins. However, application of these techniques is hampered by the duration and the complexity of the photoreaction of channelrhodopsin compared with other microbial rhodopsins. A particular interest in resolving the molecular mechanism lies in the structural changes that lead to channel opening and closure. Here, we review the current structural and mechanistic knowledge that has been accomplished by integrating the static structure provided by X-ray crystallography and electron microscopy with time-resolved spectroscopic and electrophysiological techniques. The dynamical reactions of the chromophore are effectively coupled to structural changes of the protein, as shown by ultrafast spectroscopy. The hierarchical sequence of structural changes in the protein backbone that spans the time range from 10(-12)s to 10(-3)s prepares the channel to open and, consequently, cations can pass. Proton transfer reactions that are associated with channel gating have been resolved. In particular, glutamate 253 and aspartic acid 156 were identified as proton acceptor and donor to the retinal Schiff base. The reprotonation of the latter is the critical determinant for channel closure. The proton pathway that eventually leads to proton pumping is also discussed. This article is part of a Special Issue entitled: Retinal Proteins - You can teach an old dog new tricks., (Copyright © 2013 Elsevier B.V. All rights reserved.)

Vibrational circular dichroism (VCD) spectra of aqueous solutions of proline were recorded in the course of titrations from basic to acidic pH using a spectrometer equipped with a quantum cascade laser (QCL) as an infrared light source in the spectral range from 1320 to 1220 cm(-1). The pH-dependent spectra were analyzed by singular value decomposition and global fitting of a two-pK Henderson-Hasselbalch model. The analysis delivered relative fractions of the three different protonation species. Their agreement with the relative fractions obtained from performing the same analysis on pH-dependent Fourier transform infrared (FT-IR) and QCL-IR spectra validates the quantitative results from QCL-VCD. Global fitting of the pH-dependent VCD spectra of L-proline allowed for extraction of pure spectra corresponding to anionic, zwitterionic, and cationic L-proline. From a static experiment, only pure spectra of the zwitterion would be accessible in a straightforward way. A comparison to VCD spectra calculated for all three species led to assignment of vibrational modes that are characteristic for the respective protonation states. The study demonstrates the applicability of QCL-VCD both for quantitative evaluation and for qualitative interpretation of dynamic processes in aqueous solutions.

The discovery of the light-gated ion channel channelrhodopsin (ChR) set the stage for the novel field of optogenetics, where cellular processes are controlled by light. However, the underlying molecular mechanism of light-induced cation permeation in ChR2 remains unknown. Here, we have traced the structural changes of ChR2 by time-resolved FTIR spectroscopy, complemented by functional electrophysiological measurements. We have resolved the vibrational changes associated with the open states of the channel (P(2)(390) and P(3)(520)) and characterized several proton transfer events. Analysis of the amide I vibrations suggests a transient increase in hydration of transmembrane α-helices with a t(1/2) = 60 μs, which tallies with the onset of cation permeation. Aspartate 253 accepts the proton released by the Schiff base (t(1/2) = 10 μs), with the latter being reprotonated by aspartic acid 156 (t(1/2) = 2 ms). The internal proton acceptor and donor groups, corresponding to D212 and D115 in bacteriorhodopsin, are clearly different from other microbial rhodopsins, indicating that their spatial position in the protein was relocated during evolution. Previous conclusions on the involvement of glutamic acid 90 in channel opening are ruled out by demonstrating that E90 deprotonates exclusively in the nonconductive P(4)(480) state. Our results merge into a mechanistic proposal that relates the observed proton transfer reactions and the protein conformational changes to the gating of the cation channel.

The determination of protein conformational changes induced by the interaction of substrates with secondary transporters is an important step toward the elucidation of their transport mechanism. Since conformational changes in a protein alter its vibrational patterns, they can be detected with high sensitivity by infrared difference (IR(diff)) spectroscopy without the need for external probes. We describe a general procedure to obtain substrate-induced IR(diff) spectra by alternating perfusion of buffers over an attenuated total reflection (ATR) crystal containing an adhered film of a membrane protein reconstituted in lipids. As an example, we provide specific protocols to obtain melibiose and Na(+)-induced ATR-IR(diff) spectra of reconstituted melibiose permease, a sodium/melibiose co-transporter from E. coli. The presented methodology is applicable in principle to any membrane protein, provided that it can be purified and reconstituted in functional form, and appropriate substrates are available.

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

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.

The melibiose carrier from Escherichia coli (MelB) couples the accumulation of the disaccharide melibiose to the downhill entry of H(+), Na(+), or Li(+). In this work, substrate-induced FTIR difference spectroscopy was used in combination with fluorescence spectroscopy to quantitatively compare the conformational properties of MelB mutants, implicated previously in sodium binding, with those of a fully functional Cys-less MelB permease. The results first suggest that Asp55 and Asp59 are essential ligands for Na(+) binding. Secondly, though Asp124 is not essential for Na(+) binding, this acidic residue may play a critical role, possibly by its interaction with the bound cation, in the full Na(+)-induced conformational changes required for efficient coupling between the ion- and sugar-binding sites; this residue may also be a sugar ligand. Thirdly, Asp19 does not participate in Na(+) binding but it is a melibiose ligand. The location of these residues in two independent threading models of MelB is consistent with their proposed role.

Efficient retinal photoisomerization, signal transduction, and amplification contribute to single-photon electrical responses in vertebrates visual cells. However, spontaneous discrete electrical signals arising in the dark, with identical intensity and time profiles as those generated by genuine single photons (dark events), limit the potential capability of the rod visual system to discern single photons from thermal noise. It is accepted that the light and the thermal activation of the rod photoreceptor rhodopsin (Rho) triggers the light and the dark events, respectively. However the activation barrier for the dark events (80-110 kJ/mol) appears to be only half of the barrier for light-dependent activation of Rho (> or =180 kJ/mol). On the basis of these observations, it has been postulated that both processes should follow different pathways, but the molecular mechanism for the thermal activation process still remains an open question and subject of debate. Here, performing infrared difference spectroscopy measurements, we found that the -OH group of Thr118 from bovine Rho exhibits a slow but measurable hydrogen/deuterium exchange (HDX) under native conditions. Given the location of Thr118 in the X-ray structures, isolated from the aqueous phase and in steric contact with the buried retinal chromophore, we assume that a protein structural fluctuation must drive the retinal binding pocket (RBP) transiently open. We characterized the kinetics (rate and activation enthalpy) and thermodynamics (equilibrium constant and enthalpy) of this fluctuation from the global analysis of the HDX of Thr118-OH as a function of the temperature and pH. In parallel, using HPLC chromatography, we determined the kinetics of the thermal isomerization of the protonated 11-cis retinal in solution, as a model for retinal thermal isomerization in an open RBP. Finally, we propose a quantitative two-step model in which the dark activation of Rho is triggered by thermal isomerization of the retinal in a transiently opened RBP, which accurately reproduced both the experimental activation barrier and the rate of the dark events. We conclude that the absolute sensitivity threshold of our visual system is limited by structural fluctuations of the chromophore binding pocket rather than in the chromophore itself.

Attenuated total reflection infrared (ATR-IR) difference spectroscopy stands out because of its ability to provide information on the interaction of substrates with membrane proteins in their native lipid bilayer environment. We show how the study and interpretation of the structural changes in membrane proteins upon substrate binding is simplified by obtaining ATR-IR difference spectra with polarized light and then computing the difference spectra in the z and x,y directions, where structural and orientation changes give specific difference absorbance patterns. In combination with a maximum-entropy band-narrowing method and some simple spectroscopic rules, the present approach allows us to unambiguously identify changes in the tilt of some helices in the secondary transporter melibiose permease following melibiose binding in the presence of sodium, suggesting the formation of an occluded state during the transport mechanism of the substrates.

Overlapped bands often appear in applications of infrared spectroscopy, for instance in the analysis of the amide I band of proteins. Fourier self-deconvolution (FSD) is a popular band-narrowing mathematical method, allowing for the resolution of overlapped bands. The filter function used in FSD plays a significant role in the factor by which the deconvolved bands are actually narrowed (the effective narrowing), as well as in the final signal-to-noise degradation induced by FSD. Moreover, the filter function determines, to a good extent, the band-shape of the deconvolved bands. For instance, the intensity of the harmful side-lobule oscillations that appear in over-deconvolution depends importantly on the filter function used. In the present paper we characterized the resulting band shape, effective narrowing, and signal-to-noise degradation in infra-, self-, and over-deconvolution conditions for several filter functions: Triangle, Bessel, Hanning, Gaussian, Sinc2, and Triangle2. We also introduced and characterized new filters based on the modification of the Blackmann filter. Our conclusion is that the Bessel filter (in infra-, self-, and mild over-deconvolution), the newly introduced BL3 filter (in self- and mild/moderate over-deconvolution), and the Gaussian filter (in moderate/strong over-deconvolution) are the most suitable filter functions to be used in FSD.

Vectorial transport by pumps requires a switch in the accessibility (or affinity) of the ion binding site from the extracelullar to the cytoplamic side, or vice versa. In the proton-pump bacteriorhodopsin (bR) the nature of this switch mechanism is still controversial, although it is expected to occur during the transition between two M substates. Here, we characterized this transition by time-resolved Fourier transform infrared (FT-IR) spectroscopy under functional conditions, using a novel approach for the analysis of kinetic data: the regularized inversion of an eigenvalue problem. The use of IR spectroscopy allowed the simultaneous evaluation of the involvement of the protein backbone, retinal, amino acid side chains, and internal water molecules in the switch mechanism. We provide solid evidence that the switch is not associated with protein backbone conformational changes. On the other hand, changes in the retinal conformation (or in the orientation of the Schiff Base (SB) during the switch) are reasonably although not completely discarded. We found that the proton release group (PRG), a delocalized proton characterized by a broad continuum band in the infrared, deprotonates in the transition between two M substates. Vectorial proton transport is most likely guaranteed by the coupled proton affinity changes resulting from the PRG deprotonation, favoring an affinity-based over an accessibility-based switch mechanism.

Proline residues play a fundamental and subtle role in the dynamics, structure, and function in many membrane proteins. Temperature derivative spectroscopy and differential scanning calorimetry have been used to determine the effect of proline substitution in the structural stability of the active site and transmembrane arrangement of bacteriorhodopsin. We have analyzed the Pro-to-Ala mutation for the helix-embedded prolines Pro50, Pro91, and Pro186 in the native membrane environment. This information has been complemented with the analysis of the respective crystallographic structures by the FoldX force field. Differential scanning calorimetry allowed us to determine distorted membrane arrangement for P50A and P186A. The protein stability was severely affected for P186A and P91A. In the case of Pro91, a single point mutation is capable of strongly slowing down the conformational diffusion along the denaturation coordinate, becoming a barrier-free downhill process above 371 K. Temperature derivative spectroscopy, applied for first time to study thermal stability of proteins, has been used to monitor the stability of the active site of bacteriorhodopsin. The mutation of Pro91 and Pro186 showed the most striking effects on the retinal binding pocket. These residues are the Pro in closer contact to the active site (activation energies for retinal release of 60.1 and 76.8 kcal/mol, respectively, compared to 115.8 kcal/mol for WT). FoldX analysis of the protein crystal structures indicates that the Pro-to-Ala mutations have both local and long-range effects on the structural stability of residues involved in the architecture of the protein and the active site and in the proton pumping function. Thus, this study provides a complete overview of the substitution effect of helix-embedded prolines in the thermodynamic and dynamic stability of a membrane protein, also related to its structure and function.

We introduce a method for the estimation of the mean Lorentzian bandwidth of the component bands in a spectrum. The method is computationally simple, using only the module of the Fourier transform of the spectrum, and its first derivative. Moreover, the presented method does not require knowledge of the number of bands in the spectrum, their band positions, or their band areas. Furthermore, it works on spectra containing Lorentzian bands, as well as Gaussian and Voigtian bands. Therefore, the introduced method seems especially well suited for obtaining a representative Lorentzian width for highly overlapped bands, independent of their number and Lorentzian/Gaussian character. We describe how different experimental limitations (spectral truncation, offset error, presence of noise, etc.) may affect the performance of the method, and when required we propose effective alternatives to minimize their effects. Finally, we show the application of the method to an experimental spectrum: the amide I band of a dry film of the solubilized ADP/ATP carrier. The estimation of the mean Lorentzian width can allow, for instance, for a more objective selection of the deconvolution width in Fourier self-deconvolution, allowing for a more objective and reliable analysis of the amide I band of proteins. The mean Lorentzian width can also be useful to obtain an estimation of the homogenous broadening and vibrational relaxation of the amide I vibration of proteins, without requiring complex pump-probe experiments.

Analysis of infrared polarized absorbance spectra and linear dichroism spectra of reconstituted melibiose permease from Escherichia coli shows that the oriented structures correspond mainly to tilted transmembrane alpha-helices, forming an average angle of approximately 26 degrees with the membrane normal in substrate-free medium. Examination of the deconvoluted linear dichroism spectra in H(2)O and D(2)O makes apparent two populations of alpha-helices differing by their tilt angle (helix types I and II). Moreover, the average helical tilt angle significantly varies upon substrate binding: it is increased upon Na(+) binding, whereas it decreases upon subsequent melibiose binding in the presence of Na(+). In contrast, melibiose binding in the presence of H(+) causes virtually no change in the average tilt angle. The data also suggest that the two helix populations change their tilting and H/D exchange level in different ways depending on the bound substrate(s). Notably, cation binding essentially influences type I helices, whereas melibiose binding modifies the tilting of both helix populations.

We present time-resolved room-temperature infrared difference spectra for the bacteriorhodopsin (bR) photocycle at 8 cm (-1) spectral and 5 micros temporal resolution, from 4000 to 800 cm (-1). An in situ hydration method allowed for a controlled and stable sample hydration (92% relative humidity), largely improving the quality of the data without affecting the functionality of bR. Experiments in both H 2 (16)O and H 2 (18)O were conducted to assign bands to internal water molecules. Room-temperature difference spectra of the L and M intermediates minus the bR ground state (L-BR and M-BR, respectively) were comprehensively compared with their low-temperature counterparts. The room-temperature M-BR spectrum was almost identical to that obtained at 230 K, except for a continuum band. The continuum band contains water vibrations from this spectral comparison between H 2 (16)O and H 2 (18)O, and no continuum band at 230 K suggests that the protein/solvent dynamics are insufficient for deprotonation of the water cluster. On the other hand, an intense positive broadband in the low-temperature L-BR spectrum (170 K) assigned to the formation of a water cavity in the cytoplasmic domain is absent at room temperature. This water cavity, proposed to be an essential feature for the formation of L, seems now to be a low-temperature artifact caused by restricted protein dynamics at 170 K. The observed differences between low- and room-temperature FTIR spectra are further discussed in light of previously reported dynamic transitions in bR. Finally, we show that the kinetics of the transient heat relaxation of bR after photoexcitation proceeds as a thermal diffusion process, uncorrelated with the photocycle itself.

The behavior of the D115A mutant was analyzed by time-resolved UV-Vis and Fourier transformed infrared (FTIR) spectroscopies, aiming to clarify the role of Asp115 in the intra-protein signal transductions occurring during the bacteriorhodopsin photocycle. UV-Vis data on the D115A mutant show severely desynchronized photocycle kinetics. FTIR data show a poor transmission of the retinal isomerization to the chromoprotein, evidenced by strongly attenuated helical changes (amide I), the remarkable absence of environment alterations and protonation/deprotonation events related to Asp96 and direct Schiff base (SB) protonation form the bulk. This argues for the interactions of Asp115 with Leu87 (via water molecule) and Thr90 as key elements for the effective and vectorial proton path between Asp96 and the SB, in the cytoplasmic half of bacteriorhodopsin. The results strongly suggest the presence of a regulation motif enclosed in helices C and D (Thr90-Pro91/Asp115) which drives properly the dynamics of helix C through a set of interactions. It also supports the idea that intra-helical hydrogen bonding clusters in the buried regions of transmembrane proteins can be potential elements in intra-protein signal transduction.

Time-resolved spectroscopy is often used to monitor the relaxation processes (or reactions) of physical, chemical, and biochemical systems after some fast physical or chemical perturbation. Time-resolved spectra contain information about the relaxation kinetics, in the form of macroscopic time constants of decay and their decay associated spectra. In the present paper we show how the Bayesian maximum entropy inversion of the Laplace transform (MaxEnt-iLT) can provide a lifetime distribution without sign-restrictions (or two-dimensional (2D)-lifetime distribution), representing the most probable inference given the data. From the reconstructed (2D) lifetime distribution it is possible to obtain the number of exponentials decays, macroscopic rate constants, and exponential amplitudes (or their decay associated spectra) present in the data. More importantly, the obtained (2D) lifetime distribution is obtained free from pre-conditioned ideas about the number of exponential decays present in the data. In contrast to the standard regularized maximum entropy method, the Bayesian MaxEnt approach automatically estimates the regularization parameter, providing an unsupervised and more objective analysis. We also show that the regularization parameter can be automatically determined by the L-curve and generalized cross-validation methods, providing (2D) lifetime reconstructions relatively close to the Bayesian best inference. Finally, we propose the use of MaxEnt-iLT for a more objective discrimination between data-supported and data-unsupported quantitative kinetic models, which takes both the data and the analysis limitations into account. All these aspects are illustrated with realistic time-resolved Fourier transform infrared (FT-IR) synthetic spectra of the bacteriorhodopsin photocycle.

The number, position, area, and width of the bands in a lifetime distribution give the number of exponentials present in time-resolved data and their time constants, amplitudes, and heterogeneities. The maximum entropy inversion of the Laplace transform (MaxEnt-iLT) provides a lifetime distribution from time-resolved data, which is very helpful in the analysis of the relaxation of complex systems. In some applications both positive and negative values for the lifetime distribution amplitudes are physical, but most studies to date have focused on positive-constrained solutions. In this work, we first discuss optimal conditions to obtain a sign-unrestricted maximum entropy lifetime distribution, i.e., the selection of the entropy function and the regularization value. For the selection of the regularization value we compared four methods: the chi2 criterion and Bayesian inference (already used in sign-restricted MaxEnt-iLT), and the L-curve and the generalized cross-validation methods (not yet used in MaxEnt-iLT to our knowledge). Except for the frequently used chi2 criterion, these methods recommended similar regularization values, providing close to optimum solutions. However, even when an optimal entropy function and regularization value are used, a MaxEnt lifetime distribution will contain noise-induced errors, as well as systematic distortions induced by the entropy maximization (regularization-induced errors). We introduce the concept of the apparent resolution function in MaxEnt, which allows both the noise and regularization-induced errors to be estimated. We show the capability of this newly introduced concept in both synthetic and experimental time-resolved Fourier transform infrared (FT-IR) data from the bacteriorhodopsin photocycle.

The visual pigment rhodopsin has been extensively studied for the kinetics of its photointermediates by various spectroscopic methods. Unlike such archaeal retinal proteins as bacteriorhodopsin, visual rhodopsin does not thermally recover its dark state after photoexcitation, which precludes repeated excitation of a single sample and thereby complicates time-resolved experiments. Kinetic data on the late rhodopsin photointermediates have so far been available mainly from time-resolved ultraviolet (UV)-visible spectroscopy, but not from Fourier transform infrared (FTIR) spectroscopy. The latter has the advantage of being informative of structural changes of both chromophore and protein, but does not allow the highly reproducible, automated sample exchange procedures available to UV-visible spectroscopy. Using rapid-scan FTIR difference spectroscopy, we obtained time-resolved data sets that were analyzed by a maximum entropy inverse Laplace-transform. Covering the time range from 8 ms to 15 s at temperatures of 0 and -7 degrees C, the transitions from the Lumi to the Meta I and from the Meta I to the Meta II photoproduct states could be resolved. In the transition from Meta I to Meta II, our data reveal a partial deprotonation of the retinal Schiff base preceding the conformational change of the receptor protein to Meta II. The technique and the results are discussed in regard to its advantages as well as its limitations., (Copyright 2006 Wiley Periodicals, Inc.)

We present a method for the analysis of time-resolved spectroscopic data following first-order kinetics. The time traces at all the available spectroscopic channels (e.g., wavelength or wavenumber) are inverse Laplace transformed. The transformation is stabilized by the maximum entropy method generalized for solutions without sign-restriction. In this way, time-resolved spectra can be converted to lifetime-resolved spectra, where bands appear at coordinates corresponding to their spectroscopic maxima and time constant of appearance (negative amplitude) or disappearance (positive amplitude). From the lifetime-resolved spectra, the number of exponentially decaying components, their time constants, and their decay-associated spectra are readily available. Moreover, since bands are spread in two dimensions extra band-resolution is possible. We named this method of transforming time-resolved spectra into lifetime-resolved spectra multi-spectroscopic channel maximum entropy inversion of the Laplace transform (M-MaxEnt-iLT). The basis of M-MaxEnt-iLT is presented in detail and its properties and limitations are thoroughly discussed. We also show how the combination of M-MaxEnt-iLT with spectral smoothing or deconvolution can improve the appearance and/or band resolution of the obtained lifetime-resolved spectra.

Absorbance and difference infrared spectra are often acquired aiming to characterize protein structure and structural changes of proteins upon ligand binding, as well as for many other chemical and biochemical studies. Their analysis requires as a first step the identification of the component bands (number, position, and area) and as a second step their assignment. The first step of the analysis is challenged by the habitually strong band overlap in infrared spectra. Therefore, it is useful to make use of a mathematical method able to narrow the component bands to the extent to eliminate, or at least reduce, the band overlap. Additionally, to be of general applicability this method should permit negative values for the solution. We present a maximum entropy deconvolution approach for the band-narrowing of absorbance and difference spectra showing the required characteristics, which uses the generalized negative Burg-entropy (Itakura-Saito discrepancy) generalized for difference spectra. We present results on synthetic noisy absorbance and difference spectra, as well as on experimental infrared spectra from the membrane protein bacteriorhodopsin.

Fourier transform infrared difference spectroscopy has been used to obtain information about substrate-induced structural changes of the melibiose permease (MelB) from Escherichia coli reconstituted into liposomes. Binding of the cosubstrate Na(+) gives rise to several peaks in the amide I and II regions of the difference spectrum Na(+).MelB minus H(+).MelB, that denote the presence of conformational changes in all types of secondary structures (alpha-helices, beta-sheets, loops). In addition, peaks around 1400 and at 1740-1720 cm(-1) are indicative of changes in protonation/deprotonation or in environment of carboxylic groups. Binding of the cosubstrate Li(+) produces a difference spectrum that is also indicative of conformational changes, but that is at variance as compared to that induced by Na(+) binding. To analyze the following transport steps, the melibiose permease with either H(+), Na(+), or Li(+) bound was incubated with melibiose. The difference spectra obtained by subtracting the spectrum cation.MelB from the respective complex cation.melibiose.MelB were roughly similar among them, but different from those induced by cation binding, and more intense. Therefore, major conformational changes that are induced during melibiose binding/substrate translocation, like those denoted by intense peaks at 1668 and 1645 cm(-)(1), are similar for the three cotransporting cations. Changes in the protonation state and/or in the environment of given carboxylic residues were also induced by melibiose-MelB interaction in the presence of cations.

Estimation of the band parameters of overlapped bands often relies on curve-fitting. It has been demonstrated that curve-fitting provides the maximum likelihood estimation of band parameters under a series of assumptions. One of these assumptions is that the curve-fitting model is correct and any error in the data is random. Under real conditions, we have to acknowledge the unavoidable presence of errors in the model and systematic errors in the data. Here, we derive an expression for the estimation of how these errors affect the quality of the parameters obtained from curve-fitting. In addition, we derive theoretical expressions to quantify the extent to which different methods can improve the curve-fitting robustness to these errors. The methods considered are: (i) deterministic and (ii) probabilistic constraints in the band parameters, (iii) curve-fitting band-narrowed data, and (iv) building a more accurate model. The theoretical expressions obtained are tested in the curve-fitting of a synthetic noisy spectrum with either baseline or band shape errors, and in the curve-fitting of the experimental infrared amide I band of the membrane protein bacteriorhodopsin.

We present a general method for curve-fitting Fourier manipulated spectra, comprising apodized, smoothed, derivatised and deconvoluted spectra. The analytical expressions of Fourier manipulated bands in the spectral domain, needed for the curve-fitting, are usually very complex or do not even exist; hence an accurate curve-fit of Fourier manipulated spectra becomes unfeasible. Our strategy is to construct both the model and their derivatives in the Fourier domain, where they have simple and general expressions, and then Fourier transform them back to the spectral domain. The first benefit of this approach is the accurate curve-fitting of Fourier deconvoluted spectra, a main step in the secondary structure estimation of proteins by FTIR spectroscopy.

The ADP/ATP transporter shows a high instability when solubilized, making it difficult to obtain functional protein with sufficient purity for long-term spectroscopic studies. When solubilized in the detergent dodecyl maltoside the protein is in equilibrium between the so-called CATR and BA conformations and in a few hours it becomes nonfunctional, unable to bind either its inhibitors or its substrates. By Fourier transform infrared spectroscopy, we studied the structural changes involved in this denaturation process. To do so, the carboxyatractyloside-inhibited protein was used as a structural model for the protein in the CATR conformation and its spectrum was compared with that of the unliganded time-inactivated protein. From the difference spectra of the amide I, amide II, and amide A bands combined with dichroism spectra of the carboxyatractyloside-inhibited protein, we concluded that few structural differences exist between both states, affecting as few as 11 amino acids (3.5% of the protein); the structural changes consisted in the disappearance of large loop structure and the appearance of aggregated strands. We hypothesize that some mitochondrial loop (tentatively loop M1) shows a high tendency to aggregate, being responsible for the observed features. The functional consequences of this hypothesis are discussed.

The accessibility of Escherichia coli melibiose permease to aqueous solvent was studied following hydrogen-deuterium exchange kinetics monitored by attenuated total reflection-Fourier transform infrared spectroscopy under four distinct conditions where MelB forms different complexes with its substrates (H(+), Na(+), melibiose). Analysis of the amide II band upon (2)H(2)O exposure discloses a significant sugar protection of the protein against aqueous solvent, resulting in an 8% less exchange of the corresponding H(+)*melibiose*MelB complex compared with the protein in the absence of sugar. Investigation of the amide I exchange reveals clear substrate effects on beta-sheet accessibility, with the complex H(+)*melibiose*MelB being the most protected state against exchange, followed by Na(+)*melibiose*MelB. Although of smaller magnitude, similar changes in alpha-helices plus non-ordered structures are detected. Finally, no differences are observed when analyzing reverse turn structures. The results suggest that sugar binding induces a remarkable compactness of the carrier's structure, affecting mainly beta-sheet domains of the transporter, which, according to secondary structure predictions, may include cytoplasmic loops 4-5 and 10-11. A possible catalytic role of these two loops in the functioning of MelB is hypothesized.