Interpretation of Amide I Difference Bands Observed during Protein Reactions Using Site-Directed Isotopically Labeled Bacteriorhodopsin as a Model System

De Novo Designed Heterochiral Blue Fluorescent Protein

Diversification of chain stereochemistry offers a tremendous increase in protein design space. We have designed a minimal fluorescent protein, pregnant with β-(1-azulenyl)-l-alanine in the hydrophobic core of a heterotactic protein scaffold, employing automated design tools such as automated repetitive simulated annealing molecular dynamics and IDeAS. The de novo designed heterochiral protein can be selectively excited at 342 nm, quite distant from the intrinsic fluorophore, and emits in the blue region. The structure and stability of the designed proteins were evaluated by established spectroscopic and calorimetric methods.

A Different hIAPP Polymorph Is Observed in Human Serum Than in Aqueous Buffer: Demonstration of a New Method for Studying Amyloid Fibril Structure Using Infrared Spectroscopy

There is enormous interest in measuring amyloid fibril structures, but most structural studies measure fibril formation in vitro using aqueous buffer. Ideally, one would like to measure fibril structure and mechanism under more physiological conditions. Toward this end, we have developed a method for studying amyloid fibril structure in human serum. Our approach uses isotope labeling, antibody depletion of the most abundant proteins (albumin and IgG), and infrared spectroscopy to measure aggregation in human serum with reduced protein content. Reducing the nonamyloid protein content enables the measurements by decreasing background signals but retains the full composition of salts, sugars, metal ions, etc. that are naturally present but usually missing from in vitro studies. We demonstrate the method by measuring the two-dimensional infrared (2D IR) spectra of isotopically labeled human islet amyloid polypeptide (hIAPP or amylin). We find that the fibril structure of hIAPP formed in serum differs from that formed via aggregation in aqueous buffer at residues Gly24 and Ala25, which reside in the putative "amyloidogenic core" or FGAIL region of the sequence. The spectra are consistent with extended parallel stacks of strands consistent with β-sheet-like structure, rather than a partially disordered loop that forms in aqueous buffer. These experiments provide a new method for using infrared spectroscopy to monitor the structure of proteins under physiological conditions and reveal the formation of a significantly different polymorph structure in the most important region of hIAPP.

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.

The arrestin-1 finger loop interacts with two distinct conformations of active rhodopsin

Signaling of the prototypical G protein-coupled receptor (GPCR) rhodopsin through its cognate G protein transducin (G_t) is quenched when arrestin binds to the activated receptor. Although the overall architecture of the rhodopsin/arrestin complex is known, many questions regarding its specificity remain unresolved. Here, using FTIR difference spectroscopy and a dual pH/peptide titration assay, we show that rhodopsin maintains certain flexibility upon binding the "finger loop" of visual arrestin (prepared as synthetic peptide ArrFL-1). We found that two distinct complexes can be stabilized depending on the protonation state of E3.49 in the conserved (D)ERY motif. Both complexes exhibit different interaction modes and affinities of ArrFL-1 binding. The plasticity of the receptor within the rhodopsin/ArrFL-1 complex stands in contrast to the complex with the C terminus of the G_t α-subunit (GαCT), which stabilizes only one specific substate out of the conformational ensemble. However, G_t α-subunit binding and both ArrFL-1-binding modes involve a direct interaction to conserved R3.50, as determined by site-directed mutagenesis. Our findings highlight the importance of receptor conformational flexibility and cytoplasmic proton uptake for modulation of rhodopsin signaling and thereby extend the picture provided by crystal structures of the rhodopsin/arrestin and rhodopsin/ArrFL-1 complexes. Furthermore, the two binding modes of ArrFL-1 identified here involve motifs of conserved amino acids, which indicates that our results may have elucidated a common modulation mechanism of class A GPCR-G protein/-arrestin signaling.

Tracking Reactive Water and Hydrogen-Bonding Networks in Photosynthetic Oxygen Evolution

In oxygenic photosynthesis, photosystem II (PSII) converts water to molecular oxygen through four photodriven oxidation events at a Mn₄CaO₅ cluster. A tyrosine, YZ (Y161 in the D1 polypeptide), transfers oxidizing equivalents from an oxidized, primary chlorophyll donor to the metal center. Calcium or its analogue, strontium, is required for activity. The Mn₄CaO₅ cluster and YZ are predicted to be hydrogen bonded in a water-containing network, which involves amide carbonyl groups, amino acid side chains, and water. This hydrogen-bonded network includes amino acid residues in intrinsic and extrinsic subunits. One of the extrinsic subunits, PsbO, is intrinsically disordered. This extensive (35 Å) network may be essential in facilitating proton release from substrate water. While it is known that some proteins employ internal water molecules to catalyze reactions, there are relatively few methods that can be used to study the role of water. In this Account, we review spectroscopic evidence from our group supporting the conclusion that the PSII hydrogen-bonding network is dynamic and that water in the network plays a direct role in catalysis. Two approaches, transient electron paramagnetic resonance (EPR) and reaction-induced FT-IR (RIFT-IR) spectroscopies, were used. The EPR experiments focused on the decay kinetics of YZ• via recombination at 190 K and the solvent isotope, pH, and calcium dependence of these kinetics. The RIFT-IR experiments focused on shifts in amide carbonyl frequencies, induced by photo-oxidation of the metal cluster, and on the isotope-based assignment of bands to internal, small protonated water clusters at 190, 263, and 283 K. To conduct these experiments, PSII was prepared in selected steps along the catalytic pathway, the S_n state cycle (n = 0-4). This cycle ultimately generates oxygen. In the EPR studies, S-state dependent changes were observed in the YZ• lifetime and in its solvent isotope effect. The YZ• lifetime depended on the presence of calcium at pH 7.5, but not at pH 6.0, suggesting a two-donor model for PCET. At pH 6.0 or 7.5, barium and ammonia both slowed the rate of YZ• recombination, consistent with disruption of the hydrogen-bonding network. In the RIFT-IR studies of the S state transitions, infrared bands associated with the transient protonation and deprotonation of internal waters were identified by D₂O and H₂¹⁸O labeling. The infrared bands of these protonated water clusters, W_n⁺ (or nH₂O(H₃O)⁺, n = 5-6), exhibited flash dependence and were produced during the S₁ to S₂ and S₃ to S₀ transitions. Calcium dependence was observed at pH 7.5, but not at pH 6.0. S-state induced shifts were observed in amide C═O frequencies during the S₁ to S₂ transition and attributed to alterations in hydrogen bonding, based on ammonia sensitivity. In addition, isotope editing of the extrinsic subunit, PsbO, established that amide vibrational bands of this lumenal subunit respond to the S state transitions and that PsbO is a structural template for the reaction center. Taken together, these spectroscopic results support the hypothesis that proton transfer networks, extending from YZ to PsbO, play a functional and dynamic role in photosynthetic oxygen evolution.

Redox-linked conformational control of proton-coupled electron transfer: Y122 in the ribonucleotide reductase β2 subunit

Tyrosyl radicals play essential roles in biological proton-coupled electron transfer (PCET) reactions. Ribonucleotide reductase (RNR) catalyzes the reduction of ribonucleotides and is vital in DNA replication in all organisms. Class Ia RNRs consist of α2 and β2 homodimeric subunits. In class Ia RNR, such as the E. coli enzyme, an essential tyrosyl radical (Y122O(•))-diferric cofactor is located in β2. Although Y122O(•) is extremely stable in free β2, Y122O(•) is highly reactive in the quaternary substrate-α2β2 complex and serves as a radical initiator in catalytic PCET between β2 and α2. In this report, we investigate the structural interactions that control the reactivity of Y122O(•) in a model system, isolated E. coli β2. Y122O(•) was reduced with hydroxyurea (HU), a radical scavenger that quenches the radical in a clinically relevant reaction. In the difference FT-IR spectrum, associated with this PCET reaction, amide I (CO) and amide II (CN/NH) bands were observed. Specific (13)C-labeling of the tyrosine C1 carbon assigned a component of these bands to the Y122-T123 amide bond. Comparison to density functional calculations on a model dipeptide, tyrosine-threonine, and structural modeling demonstrated that PCET is associated with a Y122 rotation and a 7.2 Å translation of the Y122 phenolic oxygen. To test for the functional consequences of this structural change, a proton inventory defined the origin of the large solvent isotope effect (SIE = 16.7 ± 1.0 at 25 °C) on this reaction. These data suggest that the one-electron, HU-mediated reduction of Y122O(•) is associated with two, rate-limiting (full or partial) proton transfer reactions. One is attributable to HU oxidation (SIE = 11.9, net H atom transfer), and the other is attributable to coupled, hydrogen-bonding changes in the Y122O(•)-diferric cofactor (SIE = 1.4). These results illustrate the importance of redox-linked changes to backbone and ring dihedral angles in high potential PCET and provide evidence for rate-limiting, redox-linked hydrogen-bonding interactions between Y122O(•) and the iron cluster.

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.

Differential stability of TATA box binding proteins from archaea with different optimal growth temperatures

The TATA box binding protein (TBP) is involved in promoter recognition, the first step of transcription initiation. TBP is universally conserved and essential in archaea and eukaryotes. In archaea, TBPs have to be stable and to function in species that cover an extremely wide range of optimal growth temperatures (OGTs), from below 0 degrees C to more than 100 degrees C. Thus, the archaeal TBP family is ideally suited to study the evolutionary adaptation of proteins to an extremely wide range of temperatures. We characterized the thermostability of one mesophilic and one thermophilic TBP by infrared spectroscopy. Transition temperatures (T(m)s) of thermal unfolding have been determined using TBPs from Methanosarcina mazei (OGT 37 degrees C) and from Methanothermobacter thermautotrophicus (OGT 65 degrees C). Furthermore, the influence of protein and salt concentration on thermostability has been characterized. Together with previous studies, our results reveal that the T(m)s of archaeal TBPs are closely correlated with the OGTs of the respective species. Noteworthy, this is also true for the TBP from M. mazei representing the first characterized TBP from a mesophilic archaeon. In contrast, the only characterized eukaryotic TBP of the mesophilic plant Arabidopsis thaliana has a T(m) more than 40 degrees C above the OGT.

Site-specific difference 2D-IR spectroscopy of bacteriorhodopsin

We demonstrate the extension of the principle of difference Fourier transform infrared (FTIR) spectroscopy to difference 2D-IR spectroscopy. To this end, we measure difference 2D-IR spectra of the protein bacteriorhodopsin in its early J- and K-intermediates. By comparing with the static 2D-IR spectrum of the protonated Schiff base of all-trans retinal, we demonstrate that the 2D-IR spectrum of the all-trans retinal chromophore in bacteriorhodopsin can be measured with the background from the remainder of the protein completely suppressed. We discuss several models to interpret the detailed line shape of the difference 2D-IR spectrum.

Site-specific relaxation kinetics of a tryptophan zipper hairpin peptide using temperature-jump IR spectroscopy and isotopic labeling

Two antiparallel beta-strands connected by a turn make beta-hairpins an ideal model system to analyze the interactions and dynamics of beta-sheets. Site-specific conformational dynamics were studied by temperature-jump IR spectroscopy and isotopic labeling in a model based on the tryptophan zipper peptide, Trpzip2, developed by Cochran et al. (Proc. Natl. Acad. Sci. U.S.A. 2001, 98, 5578). The modified Trpzip2C peptides have nearly identical equilibrium spectral behavior as Trpzip2 showing that they also form well-characterized beta-hairpin conformations in aqueous solution. Selective introduction of 13C=O groups on opposite strands lead to distinguishable cross-strand coupling of the labeled residues as monitored in the amide I' band. These frequency patterns reflect theoretical predictions, and the coupled 13C=O band loses intensity with increase in temperature and unfolding of the hairpin. Thermal relaxation kinetics were analyzed for unlabeled and cross-strand isotopically labeled variants. T-jumps of approximately 10 degrees C induce relaxation times of a few microseconds that decrease with increase of the peptide temperature. Differences in kinetic behavior for the loss of beta-strand and gain of disordered structure can be used to distinguish localized structure dynamics by comparison of nonlabeled and labeled amide I' components. Analysis of the data supports multistate dynamic and equilibrium behavior, but because of this process it is not possible to clearly define a folding and unfolding rate. Nonetheless, site-specific relaxation kinetics could be seen to be consistent with a hydrophobic collapse hypothesis for hairpin folding.

Vibrational spectra from atomic fluctuations in dynamics simulations. II. Solvent-induced frequency fluctuations at femtosecond time resolution

The midinfrared (MIR) spectra of molecules in polar solvents exhibit inhomogeneously broadened bands whose spectral positions are shifted as compared to the gas phase. The shifts are caused by interactions with structured solvation shells and the broadenings by fluctuations of these interactions. The MIR spectra can be calculated from hybrid molecular dynamics (MD) simulations, which treat the solute molecule by density functional theory and the solvent by molecular mechanics by the so-called instantaneous normal mode analysis (INMA) or by Fourier transforming the time correlation function (FTTCF) of the molecular dipole moment. In Paper I of this work [M. Schmitz and P. Tavan, J. Chem. Phys. 121, 12233 (2004)] we explored an alternative method based on generalized virial (GV) frequencies noting, however, that GV systematically underestimates frequencies. As shown by us these artifacts are caused by solvent-induced fluctuations of the (i) equilibrium geometry, (ii) force constants, and (iii) normal mode directions as well as by (iv) diagonal and (v) off-diagonal anharmonicities. Here we now show, by analyzing the time scales of fluctuations and sample MD trajectories of formaldehyde in the gas phase and in water, that all these sources of computational artifacts can be made visible by a Fourier analysis of the normal coordinates. Correspondingly, the error sources (i) and (iii)-(v) can be removed by bandpass filtering, as long as the spectral signatures of the respective effects are well separated from the fundamental band. Furthermore, the artifacts arising from effect (ii) can be strongly diminished by a time-resolved version of the GV approach (TF-GV). The TF-GV method then yields for each mode j a trajectory of the vibrational frequency omega(j)(tmid R:tau) at a time resolution tau>tau(j), which is only limited by the corresponding oscillation time tau(j)=2pi/omega(j) and, thus, is in the femtosecond range. A correlation analysis of these trajectories clearly separates the librational motions from the conformational dynamics of the solvation shells and yields the inhomogeneously broadened MIR spectra, if the theory of motional narrowing is properly included. The MIR spectrum of formaldehyde in solution obtained by TF-GV agrees very well with the FTTCF result, if one applies the so-called "harmonic approximation" quantum correction factor and a temperature scaling to the FTTCF intensities. Also for INMA an excellent agreement is achieved if one disregards a slight INMA overestimate of linewidths.

Vibrational spectra from atomic fluctuations in dynamics simulations. I. Theory, limitations, and a sample application

Hybrid molecular dynamics (MD) simulations, which combine density functional theory (DFT) descriptions of a molecule with a molecular mechanics (MM) modeling of its solvent environment, have opened the way towards accurate computations of solvation effects in the vibrational spectra of molecules. Recently, Wheeler et al. [ChemPhysChem 4, 382 (2002)] have suggested to compute these spectra from DFT/MM-MD trajectories by diagonalizing the covariance matrix of atomic fluctuations. This so-called principal mode analysis (PMA) allegedly can replace the well-established approaches, which are based on Fourier transform methods or on conventional normal mode analyses. By scrutinizing and revising the PMA approach we identify five conditions, which must be guaranteed if PMA is supposed to render exact vibrational frequencies. Besides specific choices of (a) coordinates and (b) coordinate systems, these conditions cover (c) a harmonic intramolecular potential, (d) a complete thermal equilibrium within the molecule, and (e) a molecular Hamiltonian independent of time. However, the PMA conditions [(c)-(d)] and [(c)-(e)] are generally violated in gas phase DFT-MD and liquid phase DFT/MM-MD trajectories, respectively. Based on a series of simple analytical model calculations and on the analysis of MD trajectories calculated for the formaldehyde molecule in the gas phase (DFT) and in liquid water (DFT/MM) we show that in both phases the violation of condition (d) can cause huge errors in PMA frequency computations, whereas the inevitable violations of conditions (c) and (e), the latter being generic to the liquid phase, imply systematic and sizable underestimates of the vibrational frequencies by PMA. We demonstrate that the huge errors, which are caused by an incomplete thermal equilibrium violating (d), can be avoided if one introduces mode-specific temperatures T(j) and calculates the frequencies from a "generalized virial" (GV) expression instead from PMA. Concerning ways to additionally remove the remaining errors, which GV still shares with PMA, we refer to Paper II of this work [M. Schmitz and P. Tavan, J. Chem. Phys. 121, 12247 (2004)].

Consequences of Counterion Mutation in Sensory Rhodopsin II ofNatronobacterium pharaonisfor Photoreaction and Receptor Activation: An FTIR Study

In many retinal proteins the proton transfer from the Schiff base to the counterion represents a functionally important step of the photoreaction. In the signaling state of sensory rhodopsin II from Natronobacterium pharaonis this transfer has already occurred, but in the counterion mutant Asp75Asn it is blocked during all steps of the photocycle. Therefore, the study of the molecular changes during the photoreaction of this mutant should provide a deeper understanding of the activation mechanism, and for this, we have applied time-resolved step-scan FTIR spectroscopy. The photoreaction is drastically altered; only red-shifted intermediates are formed with a chromophore strongly twisted around the 14-15 single bond. In addition, the photocycle is shortened by 2 orders of magnitude. Nevertheless, a transition involving only protein changes similar to that of the wild type is observed, which has been correlated with the formation of the signaling state. However, whereas in the wild type this transition occurs in the millisecond range, it is shortened to 200 micros in the mutant. The results are discussed with respect to the altered electrostatic interactions, role of proton transfer, the published 3D structure, and physiological activity.

Fourier transform IR spectroscopy study for new insights into molecular properties and activation mechanisms of visual pigment rhodopsin

Fourier transform IR (FTIR) spectroscopy has been successfully applied in recent years to examine the functional and structural properties of the membrane protein rhodopsin, a prototype G protein coupled receptor. Unlike UV-visible spectroscopy, FTIR spectroscopy is structurally sensitive. It may give us both global information about the conformation of the protein and very detailed information about the retinal chromophore and all other functional groups, even when these are not directly related to the chromophore. Furthermore, it can be successfully applied to the photointermediates of rhodopsin, including the active receptor species, metarhodopsin II, and its decay products, which is not expected presently or even in the near future from crystallographic approaches. In this review we show how FTIR spectroscopy has significantly contributed to the understanding of very different aspects of rhodopsin, comprising both structural properties and the mechanisms leading to receptor activation and deactivation.

Time-resolved FTIR studies of sensory rhodopsin II (NpSRII) from Natronobacterium pharaonis: implications for proton transport and receptor activation

The photocycle of the photophobic receptor from Natronobacterium pharaonis, NpSRII, is studied by static and time-resolved step-scan Fourier transform infrared spectroscopy. Both low-temperature static and time-resolved spectra resolve a K-like intermediate, and the corresponding spectra show little difference within the noise of the time-resolved data. As compared to intermediate K of bacteriorhodopsin, relatively large amide I bands indicate correspondingly larger distortions of the protein backbone. The time-resolved spectra identify an intermediate L-like state with surprisingly small additional molecular alterations. With the formation of intermediate M, the Schiff-base proton is transferred to the counterion Asp-75. This state is characterized by larger amide bands indicating larger distortions of the protein. We can identify a second M state that differs only in small-protein bands. Reisomerization of the chromophore to all-trans occurs with the formation of intermediate O. The accelerated decay of intermediate M caused by azide results in another red-shifted intermediate with a protonated Schiff base. The chromophore in this state, however, still has 13-cis geometry. Nevertheless, the reisomerization is still as slow as under the conditions without azide. The results are discussed with respect to mechanisms of the observed proton pumping and the possible roles of the intermediates in receptor activation.