Monte Carlo Modeling of Spectral Diffusion Employing Multiwell Protein Energy Landscapes: Application to Pigment-Protein Complexes Involved in Photosynthesis

Identification of Residues Potentially Involved in Optical Shifts in the Water-Soluble Chlorophyll a-Binding Protein through Molecular Dynamics Simulations

Reversible light and thermally induced spectral shifts are universally observed in a wide variety of pigment-protein complexes at temperatures ranging from cryogenic to ambient. In this paper, we employed large-scale molecular dynamics (MD) simulations of a prototypical pigment-protein complex to better understand these shifts at a molecular scale. Although multiple mechanisms have been proposed over the years, no verification of these proposals via MD simulations has thus far been performed; our work represents the first step in this direction. From simulations of the water-soluble chlorophyll-binding protein complex, we determined that rearrangements of long hydrogen bonds were unlikely to be the origin of the multiwell landscape features necessary to explain observed spectral shifts. We also assessed small motions of amino acid residues and identified side chain rotations of some of these residues as likely candidates for the origin of relevant multiwell landscape features. The protein free-energy landscapes associated with side chain rotations feature energy barriers of around 1100-1600 cm-1, in agreement with optical spectroscopy results, with the most promising residue type associated with experimental signatures being serine, which possesses a symmetric triple-well landscape and moment of inertia of a relevant magnitude.

Evidence of Simultaneous Spectral Hole Burning Involving Two Tiers of the Protein Energy Landscape in Cytochrome b6f

Cytochrome b₆f, with one chlorophyll molecule per protein monomer, is a simple model system whose studies can help achieve a better understanding of nonphotochemical spectral hole burning (NPHB) and single-complex spectroscopy results obtained in more complicated photosynthetic chlorophyll-protein complexes. We are reporting new data and proposing an alternative explanation for spectral dynamics that was recently observed in cytochrome b₆f using NPHB. The relevant distribution of the tunneling parameter λ is a superposition of two components that are nearly degenerate in terms of the resultant NPHB yield and represent two tiers of the energy landscape responsible for NPHB. These two components likely burn competitively; we present the first demonstration of modeling a competitive NPHB process. Similar values of the NPHB yield result from distinctly different combinations of barrier heights, shifts along the generalized coordinate d, and/or masses of the entities involved in conformational changes m, with md² parameter different by a factor of 2.7. Consequently, in cytochrome b₆f, the first (at least) 10 h of fixed-temperature recovery preferentially probe different components of the barrier- and λ-distributions encoded into the spectral holes than thermocycling experiments. Both components most likely represent dynamics of the protein and not of the surrounding buffer/glycerol glass.

Probing Energy Landscapes of Cytochrome b6f with Spectral Hole Burning: Effects of Deuterated Solvent and Detergent

In non-photochemical spectral hole burning (NPHB) and spectral hole recovery experiments, cytochrome b₆f protein exhibits behavior that is almost independent of the deuteration of the buffer/glycerol glassy matrix containing the protein, apart from some differences in heat dissipation. On the other hand, strong dependence of the hole burning properties on sample preparation procedures was observed and attributed to a large increase of the electron-phonon coupling and shortening of the excited-state lifetime occurring when n-dodecyl β-d-maltoside (DM) is used as a detergent instead of n-octyl β-d-glucopyranoside (OGP). The data was analyzed assuming that the tunneling parameter distribution or barrier distribution probed by NPHB and encoded into the spectral holes contains contributions from two nonidentical components with accidentally degenerate excited state λ-distributions. Both components likely reflect protein dynamics, although with some small probability one of them (with larger md²) may still represent the dynamics involving specifically the -OH groups of the water/glycerol solvent. Single proton tunneling in the water/glycerol solvent environment or in the protein can be safely excluded as the origin of observed NPHB and hole recovery dynamics. The intensity dependence of the hole growth kinetics in deuterated samples likely reflects differences in heat dissipation between protonated and deuterated samples. These differences are most probably due to the higher interface thermal resistivity between (still protonated) protein and deuterated water/glycerol outside environment.