Mike Reppert

Site Selective and Single Complex Laser-Based Spectroscopies: A Window on Excited State Electronic Structure, Excitation Energy Transfer, and Electron–Phonon Coupling of Selected Photosynthetic Complexes

Site Selective and Single Complex Laser-Based Spectroscopies: A Window on Excited State Electronic Structure, Excitation Energy Transfer, and Electron Phonon Coupling of Selected Photosynthetic Complexes Ryszard Jankowiak,* Mike Reppert, Valter Zazubovich, J€org Pieper, and Tonu Reinot Department of Chemistry, Kansas State University, Manhattan, Kansas 66506, United States Department of Physics, Concordia University, Montreal H4B1R6 Quebec, Canada Max-Volmer-Laboratories for Biophysical Chemistry, Technical University of Berlin, Germany Institute of Physics, University of Tartu, Riia 142, 51014 Tartu, Estonia

On the Shape of the Phonon Spectral Density in Photosynthetic Complexes

We provide a critical assessment of typical phonon spectral densities, J(ω), used to describe linear and nonlinear optical spectra in photosynthetic complexes. Evaluation is based on a more careful comparison to experiment than has been provided in the past. J(ω) describes the frequency-dependent coupling of the system to the bath and is an important component in calculations of excitation energy transfer times. On the basis of the shape of experimental J(ω) obtained for several photosynthetic complexes, we argue that the shape of J(ω) strongly depends on the pigment-protein complex. We show that many densities (especially the Drude-Lorentz/constant damping Brownian oscillator) display qualitatively wrong behavior when compared to experiment. Because of divergence of J(ω) at zero frequency, the Brownian oscillator cannot fit a single-site spectrum correctly. It is proposed that a log-normal distribution can be used to fit experimental data and exhibits desired attributes for a physically meaningful phonon J(ω), in contrast to several commonly used spectral densities which exhibit low-frequency behavior in qualitative disagreement with experiment. We anticipate that the log-normal J(ω) function proposed in this work will be further tested in theoretical modeling of both time- and frequency-domain data.

Electrostatic frequency shifts in amide I vibrational spectra: Direct parameterization against experiment

The interpretation of protein amide I infrared spectra has been greatly assisted by the observation that the vibrational frequency of a peptide unit reports on its local electrostatic environment. However, the interpretation of spectra remains largely qualitative due to a lack of direct quantitative connections between computational models and experimental data. Here, we present an empirical parameterization of an electrostatic amide I frequency map derived from the infrared absorption spectra of 28 dipeptides. The observed frequency shifts are analyzed in terms of the local electrostatic potential, field, and field gradient, evaluated at sites near the amide bond in molecular dynamics simulations. We find that the frequency shifts observed in experiment correlate very well with the electric field in the direction of the C=O bond evaluated at the position of the amide oxygen atom. A linear best-fit mapping between observed frequencies and electric field yield sample standard deviations of 2.8 and 3.7 cm(-1) for the CHARMM27 and OPLS-AA force fields, respectively, and maximum deviations (within our data set) of 9 cm(-1). These results are discussed in the broader context of amide I vibrational models and the effort to produce quantitative agreement between simulated and experimental absorption spectra.

Identifying Residual Structure in Intrinsically Disordered Systems: A 2D IR Spectroscopic Study of the GVGXPGVG Peptide

The peptide amide-I vibration of a proline turn encodes information on the turn structure. In this study, FTIR, two-dimensional IR spectroscopy and molecular dynamics simulations were employed to characterize the varying turn conformations that exist in the GVGX(L)PGVG family of disordered peptides. This analysis revealed that changing the size of the side chain at the X amino acid site from Gly to Ala to Val substantially alters the conformation of the peptide. To quantify this effect, proline peak shifts and intensity changes were compared to a structure-based spectroscopic model. These simulated spectra were used to assign the population of type-II β turns, bulged turns, and irregular β turns for each peptide. Of particular interest was the Val variant commonly found in the protein elastin, which contained a 25% population of irregular β turns containing two peptide hydrogen bonds to the proline C═O.

Coherent two-dimensional infrared spectroscopy: Quantitative analysis of protein secondary structure in solution

We present a method to quantitatively determine the secondary structure composition of globular proteins using coherent two-dimensional infrared (2DIR) spectroscopy of backbone amide I vibrations (1550-1720 cm(-1)). Sixteen proteins with known crystal structures were used to construct a library of 2DIR spectra, and the fraction of residues in α-helix, β-sheet, and unassigned conformations was determined by singular value decomposition (SVD) of the measured two-dimensional spectra. The method was benchmarked by removing each individual protein from the set and comparing the composition extracted from 2DIR against the composition determined from the crystal structures. To highlight the increased structural content extracted from 2DIR spectra a similar analysis was also carried out using conventional infrared absorption of the proteins in the library.

Insight into the electronic structure of the CP47 antenna protein complex of photosystem II: hole burning and fluorescence study.

We report low temperature (T) optical spectra of the isolated CP47 antenna complex from Photosystem II (PSII) with a low-T fluorescence emission maximum near 695 nm and not, as previously reported, at 690-693 nm. The latter emission is suggested to result from three distinct bands: a lowest-state emission band near 695 nm (labeled F1) originating from the lowest-energy excitonic state A1 of intact complexes (located near 693 nm and characterized by very weak oscillator strength) as well as emission peaks near 691 nm (FT1) and 685 nm (FT2) originating from subpopulations of partly destabilized complexes. The observation of the F1 emission is in excellent agreement with the 695 nm emission observed in intact PSII cores and thylakoid membranes. We argue that the band near 684 nm previously observed in singlet-minus-triplet spectra originates from a subpopulation of partially destabilized complexes with lowest-energy traps located near 684 nm in absorption (referred to as AT2) giving rise to FT2 emission. It is demonstrated that varying contributions from the F1, FT1, and FT2 emission bands led to different maxima of fluorescence spectra reported in the literature. The fluorescence spectra are consistent with the zero-phonon hole action spectra obtained in absorption mode, the profiles of the nonresonantly burned holes as a function of fluence, as well as the fluorescence line-narrowed spectra obtained for the Q(y) band. The lowest Q(y) state in absorption band (A1) is characterized by an electron-phonon coupling with the Huang-Rhys factor S of approximately 1 and an inhomogeneous width of approximately 180 cm(-1). The mean phonon frequency of the A1 band is 20 cm(-1). In contrast to previous observations, intact isolated CP47 reveals negligible contribution from the triplet-bottleneck hole, i.e., the AT2 trap. It has been shown that Chls in intact CP47 are connected via efficient excitation energy transfer to the A1 trap near 693 nm and that the position of the fluorescence maximum depends on the burn fluence. That is, the 695 nm fluorescence maximum shifts blue with increasing fluence, in agreement with nonresonant hole burned spectra. The above findings provide important constraints and parameters for future excitonic calculations, which in turn should offer new insight into the excitonic structure and composition of low-energy absorption traps.

Direct observation of intermolecular interactions mediated by hydrogen bonding

Although intermolecular interactions are ubiquitous in physicochemical phenomena, their dynamics have proven difficult to observe directly, and most experiments rely on indirect measurements. Using broadband two-dimensional infrared spectroscopy (2DIR), we have measured the influence of hydrogen bonding on the intermolecular vibrational coupling between dimerized N-methylacetamide molecules. In addition to strong intramolecular coupling between N-H and C=O oscillators, cross-peaks in the broadband 2DIR spectrum appearing upon dimerization reveal strong intermolecular coupling that changes the character of the vibrations. In addition, dimerization changes the effects of intramolecular coupling, resulting in Fermi resonances between high and low-frequency modes. These results illustrate how hydrogen bonding influences the interplay of inter- and intramolecular vibrations, giving rise to correlated nuclear motions and significant changes in the vibrational structure of the amide group. These observations have direct impact on modeling and interpreting the IR spectra of proteins. In addition, they illustrate a general approach to direct molecular characterization of intermolecular interactions.

Lowest Electronic States of the CP47 Antenna Protein Complex of Photosystem II: Simulation of Optical Spectra and Revised Structural Assignments†

In this work, we present simulated steady-state absorption, emission, and nonresonant hole burning (HB) spectra for the CP47 antenna complex of photosystem II (PS II) based on fits to recently refined experimental data (Neupane et al. J. Am. Chem. Soc. 2010, 132, 4214). Excitonic simulations are based on the 2.9 Å resolution structure of the PS II core from cyanobacteria (Guskov et al. Nat. Struct. Mol. Biol. 2009, 16, 334), and allow for preliminary assignment of the chlorophylls (Chls) contributing to the lowest excitonic states. The search for realistic site energies was guided by experimental constraints and aided by simple fitting algorithms. The following experimental constraints were used: (i) the oscillator strength of the lowest-energy state should be approximately ≤0.5 Chl equivalents; (ii) the excitonic structure must explain the experimentally observed red-shifted (∼695 nm) emission maximum; and (iii) the excitonic interactions of all states must properly describe the broad (non-line-narrowed, NLN) HB spectrum (including its antihole) whose shape is extremely sensitive to the excitonic structure of the complex, especially the lowest excitonic states. Importantly, our assignments differ significantly from those previously reported by Raszewski and Renger (J. Am. Chem. Soc. 2008, 130, 4431), due primarily to differences in the experimental data simulated. In particular, we find that the lowest state localized on Chl 526 possesses too high of an oscillator strength to fit low-temperature experimental data. Instead, we suggest that Chl 523 most strongly contributes to the lowest excitonic state, with Chl 526 contributing to the second excitonic state. Since the fits of nonresonant holes are more restrictive (in terms of possible site energies) than those of absorption and emission spectra, we suggest that fits of linear optical spectra along with HB spectra provide more realistic site energies.

Computational Amide I 2D IR Spectroscopy as a Probe of Protein Structure and Dynamics

Two-dimensional infrared spectroscopy of amide I vibrations is increasingly being used to study the structure and dynamics of proteins and peptides. Amide I, a primarily carbonyl stretching vibration of the protein backbone, provides information on secondary structures as a result of vibrational couplings and on hydrogen-bonding contacts when isotope labeling is used to isolate specific sites. In parallel with experiments, computational models of amide I spectra that use atomistic structures from molecular dynamics simulations have evolved to calculate experimental spectra. Mixed quantum-classical models use spectroscopic maps to translate the structural information into a quantum-mechanical Hamiltonian for the spectroscopically observed vibrations. This allows one to model the spectroscopy of large proteins, disordered states, and protein conformational dynamics. With improvements in amide I models, quantitative modeling of time-dependent structural ensembles and of direct feedback between experiments and simulations is possible. We review the advances in developing these models, their theoretical basis, and current and future applications.

Accurate modeling of fluorescence line narrowing difference spectra: Direct measurement of the single-site fluorescence spectrum

Accurate lineshape functions for modeling fluorescence line narrowing (FLN) difference spectra (DeltaFLN spectra) in the low-fluence limit are derived and examined in terms of the physical interpretation of various contributions, including photoproduct absorption and emission. While in agreement with the earlier results of Jaaniso [Proc. Est. Acad. Sci., Phys., Math. 34, 277 (1985)] and Funfschilling et al. [J. Lumin. 36, 85 (1986)], the derived formulas differ substantially from functions used recently [e.g., M. Ratsep et al., Chem. Phys. Lett. 479, 140 (2009)] to model DeltaFLN spectra. In contrast to traditional FLN spectra, it is demonstrated that for most physically reasonable parameters, the DeltaFLN spectrum reduces simply to the single-site fluorescence lineshape function. These results imply that direct measurement of a bulk-averaged single-site fluorescence lineshape function can be accomplished with no complicated extraction process or knowledge of any additional parameters such as site distribution function shape and width. We argue that previous analysis of DeltaFLN spectra obtained for many photosynthetic complexes led to strong artificial lowering of apparent electron-phonon coupling strength, especially on the high-energy side of the pigment site distribution function.