Interpretation of the resonance Raman spectra of linear tetrapyrroles based on DFT calculations

Vibrational Spectroscopy of Phytochromes

Phytochromes are biological photoswitches that translate light into physiological functions. Spectroscopic techniques are essential tools for molecular research into these photoreceptors. This review is directed at summarizing how resonance Raman and IR spectroscopy contributed to an understanding of the structure, dynamics, and reaction mechanism of phytochromes, outlining the substantial experimental and theoretical challenges and describing the strategies to master them. It is shown that the potential of the various vibrational spectroscopic techniques can be most efficiently exploited using integral approaches via a combination of theoretical methods as well as other experimental techniques.

Redshifting and Blueshifting of β82 Chromophores in the Phycocyanin Hexamer of Porphyridium purpureum Phycobilisomes Due to Linker Proteins

Phycobilisomes in cyanobacteria and red algae are large protein complexes that absorb light and transfer energy for use in photosynthesis. The light energy absorbed by chromophores binding to phycobiliproteins in the peripheral rods can be funneled to the core through chromophores at very high efficiency. The molecular mechanism of excitation energy transfer within a phycobilisome is an example of a higher and unique function in a living organism. However, the mechanism underlying the high efficiency remains unclear. Thus, this study was carried out as a step to resolve this mechanism theoretically. The three-dimensional structure of phycobilisomes containing the linker proteins of the red alga Porphyridium purpureum was determined by cryoelectron microscopy at 2.82 Å resolution in 2020. Using these data, the absorption wavelength of each β82 chromophore in the phycocyanin hexamer located next to the core was calculated using quantum chemical treatment, considering the electric effect from its surrounding phycocyanin proteins and two linker proteins. In addition to unaffected chromophores, chromophores that were redshifted and blueshifted under the electrical influence of the two linker proteins were found. Namely, the chromophore serving as the energy sink in the rod was determined.

Excitation energy transfer and vibronic coherence in intact phycobilisomes

The phycobilisome is an oligomeric chromoprotein complex that serves as the principal mid-visible light-harvesting system in cyanobacteria. Here we report the observation of excitation-energy-transfer pathways involving delocalized optical excitations of the bilin (linear tetrapyrrole) chromophores in intact phycobilisomes isolated from Fremyella diplosiphon. By using broadband multidimensional electronic spectroscopy with 6.7-fs laser pulses, we are able to follow the progress of excitation energy from the phycoerythrin disks at the ends of the phycobilisome's rods to the C-phycocyanin disks along their length in <600 fs. Oscillation maps show that coherent wavepacket motions prominently involving the hydrogen out-of-plane vibrations of the bilins mediate non-adiabatic relaxation of a manifold of vibronic exciton states. However, the charge-transfer character of the bilins in the allophycocyanin-containing segments localizes the excitations in the core of the phycobilisome, yielding a kinetic bottleneck that enables photoregulatory mechanisms to operate efficiently on the >10-ps timescale.

On the Role of the Conserved Histidine at the Chromophore Isomerization Site in Phytochromes

Phytochromes are sensory photoreceptors that use light to drive protein structural changes, which in turn trigger physiological reaction cascades. The process starts with a double-bond photoisomerization of the linear methine-bridged tetrapyrrole chromophore in the photosensory core module. The molecular mechanism of the photoconversion depends on the structural and electrostatic properties of the chromophore environment, which are highly conserved in related phytochromes. However, the specific role of individual amino acids is yet not clear. A histidine in the vicinity of the isomerization site is highly conserved and almost invariant among all phytochromes. The present study aimed at analyzing its role by taking advantage of a myxobacterial phytochrome SaBphP1 from Stigmatella aurantiaca, where this histidine is naturally substituted with a threonine (Thr289), and comparing it to its normal, His-containing counterpart from the same organism SaBphP2 (His275). We have carried out a detailed resonance Raman and IR spectroscopic investigation of the wild-type proteins and their respective His- or Thr-substituted variants (SaBphP1-T289H and SaBphP2-H275T) using the well-characterized prototypical phytochrome Agp1 from Agrobacterium fabrum as a reference. The overall mechanism of the photoconversion is insensitive toward the His substitution. However, the chromophore geometry at the isomerization site appears to be affected, with a slightly stronger twist of ring D in the presence of Thr, which is sufficient to cause different light absorption properties in SaBphP1 and SaBphP2. Furthermore, the presence of His allows for multiple hydrogen-bonding interactions with the ring D carbonyl which may be the origin for the geometric differences of the C-D methine bridge compared to the Thr-containing variants. Other structural and mechanistic differences are independent of the presence of His. The most striking finding is the protonation of the ring C propionate in the Pfr states of SaBphP2, which is common among bathy phytochromes but so far has not been reported in prototypical phytochromes.

Functional roles of the hexamer structure of C-phycocyanin revealed by calculation of absorption wavelength

Cyanophyta-phycocyanin (C-PC) is the main constituent of the rod of phycobilisome (PBS), which is a highly ordered and large peripheral light-harvesting protein complex present on the cytoplasmic side of the thylakoid membrane in cyanobacteria and red algae. The C-PC monomer comprises two chains, α- and β-subunits, and aggregates to form ring-shaped trimers (αβ)3 with rotational symmetry. The ring-shaped trimer (αβ)3 is a structural block unit (SBU) that forms the rod of PBS. Two (αβ)3 SBUs are arranged in a face-to-face manner to form an (αβ)6 -hexamer. In this study, the electronic states of three phycocyanobilins, α84, β84, and β155 in C-phycocyanin, constituting the rod of the PBS, were calculated for both the trimer and hexamer models by considering the effect of the electrostatic field of protein moieties and water molecules. For the hexamer, the absorption wavelengths of α84, β84, and β155 were similar to those obtained experimentally; however, for the trimer, only the absorption wavelength of β155 shifted toward a shorter-wavelength. The nature of the hexamer structure as a hierarchical structure is revealed by considering the calculated absorption wavelength and energy transfer.

Protochromic absorption changes in the two-cysteine photocycle of a blue/orange cyanobacteriochrome

Cyanobacteriochromes (CBCRs) are phytochrome-related photosensors with diverse spectral sensitivities spanning the entire visible spectrum. They covalently bind bilin chromophores via conserved cysteine residues and undergo 15Z/15E bilin photoisomerization upon light illumination. CBCR subfamilies absorbing violet-blue light use an additional cysteine residue to form a second bilin-thiol adduct in a two-Cys photocycle. However, the process of second thiol adduct formation is incompletely understood, especially the involvement of the bilin protonation state. Here, we focused on the Oscil6304_2705 protein from the cyanobacterium Oscillatoria acuminata PCC 6304, which photoconverts between a blue-absorbing 15Z state ( ^15Z Pb) and orange-absorbing 15E state ( ^15E Po). pH titration analysis revealed that ^15Z Pb was stable over a wide pH range, suggesting that bilin protonation is stabilized by a second thiol adduct. As revealed by resonance Raman spectroscopy, ^15E Po harbored protonated bilin at both acidic and neutral pH, but readily converted to a deprotonated green-absorbing 15Z state ( ^15Z Pg) at alkaline pH. Site-directed mutagenesis revealed that the conserved Asp-71 and His-102 residues are required for second thiol adduct formation in ^15Z Pb and bilin protonation in ^15E Po, respectively. An Oscil6304_2705 variant lacking the second cysteine residue, Cys-73, photoconverted between deprotonated ^15Z Pg and protonated ^15E Pr, similarly to the protochromic photocycle of the green/red CBCR subfamily. Time-resolved spectroscopy revealed ^15Z Pg formation as an intermediate in the ^15E Pr-to- ^15Z Pg conversion with a significant solvent-isotope effect, suggesting the sequential occurrence of 15EP-to-15Z photoisomerization, deprotonation, and second thiol adduct formation. Our findings uncover the details of protochromic absorption changes underlying the two-Cys photocycle of violet-blue-absorbing CBCR subfamilies.

Identification of the Deprotonated Pyrrole Nitrogen of the Bilin-Based Photoreceptor by Raman Spectroscopy with an Advanced Computational Analysis

Phytochrome and cyanobacteriochrome utilize a linear methine-bridged tetrapyrrole (bilin) to control numerous biological processes. They show a reversible photoconversion between two spectrally distinct states. This photocycle is initiated by a C═C double-bond photoisomerization of the bilin followed by its thermal relaxations with transient and/or stationary changes in the protonation state of the pyrrole moiety. However, it has never been identified which of the four pyrrole nitrogen atoms is deprotonated. Here, we report a resonance Raman spectroscopic study on cyanobacteriochrome RcaE, which has been proposed to contain a deprotonated bilin for its green-absorbing 15 Z state. The observed Raman spectra were well reproduced by a simulated structure whose bilin B ring is deprotonated, with the aid of molecular dynamics and quantum mechanics/molecular mechanics calculations. The results revealed that the deprotonation of B and C rings has the distinct effect on the overall bilin structure, which will be relevant to the color tuning and photoconversion mechanisms of the phytochrome superfamily. Furthermore, this study documents the ability of vibrational spectroscopy combined with the advanced spectral analysis to visualize a proton of a cofactor molecule embedded in a protein moiety.

Protonation of the Biliverdin IXα Chromophore in the Red and Far-Red Photoactive States of a Bacteriophytochrome

The tetrapyrrole chromophore biliverdin IXα (BV) in the bacteriophytochrome from Deinococcus radiodurans (DrBphP) is usually assumed to be fully protonated, but this assumption has not been systematically validated by experiments or extensive computations. Here, we use force field molecular dynamics simulations and quantum mechanics/molecular mechanics calculations with density functional theory and XMCQDPT2 methods to investigate the effect of the five most probable protonation forms of BV on structural stability, binding pocket interactions, and absorption spectra in the two photochromic states of DrBphP. While agreement with X-ray structural data and measured UV/vis spectra suggest that in both states the protonated form of the chromophore dominates, we also find that a minor population with a deprotonated D-ring could contribute to the red-shifted tail in the absorption spectra.

Determination of the protonation preferences of bilin pigments in cryptophyte antenna complexes

The light-harvesting mechanisms of cryptophyte antenna complexes have attracted considerable attention due to their ability to exhibit maximal photosynthetic activity under very low-light conditions and to display several colors, as well as the observation of vibronic coherent features in their two-dimensional electronic spectra. However, detailed investigations on the interplay between the protein environment and their light-harvesting properties are hampered by the uncertainty related to the protonation state of the underlying bilin pigments. Here we study the protonation preferences of four types of bilin pigments including 15,16-dihydrobiliverdin (DBV), phycoerythrobilin (PEB), phycocyanobilin (PCB) and mesobiliverdin (MBV), which are found in phycoerythrin PE545 and phycocyanin PC577, PC612, PC630 and PC645 complexes. We apply quantum chemical calculations coupled to continuum solvation calculations to predict the intrinsic acidity of bilins in aqueous solution, and then combine molecular dynamics simulations with empirical pKa estimates to investigate the impact of the local protein environment on the acidity of the pigments. We also report measurements of the absorption spectra of the five complexes in a wide range of pH in order to validate our simulations and investigate possible changes in the light harvesting properties of the complexes in the range of physiological pH found in the lumen (pH ∼ 5-7). The results suggest a pKa > 7 for DBV and MBV pigments in the α polypeptide chains of PE545 and PC630/PC645 complexes, which are not coordinated to a negatively charged amino acid. For the other PEB, DBV and PCB pigments, which interact with a Glu or Asp side chain, higher pKa values (pKa > 8) are estimated. Overall, the results support a preferential population of the fully protonated state for bilins in cryptophyte complexes under physiological conditions regardless of the specific type of pigment and local protein environment.

Photoconversion mechanism of the second GAF domain of cyanobacteriochrome AnPixJ and the cofactor structure of its green-absorbing state

Cyanobacteriochromes are members of the phytochrome superfamily. In contrast to classical phytochromes, these small photosensors display a considerable variability of electronic absorption maxima. We have studied the light-induced conversions of the second GAF domain of AnPixJ, AnPixJg2, a phycocyanobilin-binding protein from the cyanobacterium Anabaena PCC 7120, using low-temperature resonance Raman spectroscopy combined with molecular dynamics simulations. AnPixJg2 is formed biosynthetically as a red-absorbing form (Pr) and can be photoconverted into a green-absorbing form (Pg). Forward and backward phototransformations involve the same reaction sequences and intermediates of similar cofactor structures as the corresponding processes in canonical phytochromes, including a transient cofactor deprotonation. Whereas the cofactor of the Pr state shows far-reaching similarities to the Pr states of classical phytochromes, the Pg form displays significant upshifts of the methine bridge stretching frequencies concomitant to the hypsochromically shifted absorption maximum. However, the cofactor in Pg is protonated and adopts a conformation very similar to the Pfr state of classical phytochromes. The spectral differences are probably related to an increased solvent accessibility of the chromophore which may reduce the π-electron delocalization in the phycocyanobilin and thus raise the energies of the first electronic transition and the methine bridge stretching modes. Molecular dynamics simulations suggest that the Z → E photoisomerization of the chromophore at the C-D methine bridge alters the interactions with the nearby Trp90 which in turn may act as a gate, allowing the influx of water molecules into the chromophore pocket. Such a mechanism of color tuning AnPixJg2 is unique among the cyanobacteriochromes studied so far.

Vibrational spectroscopic analyses of unique yellow feather pigments (spheniscins) in penguins

Many animals extract, synthesize and refine chemicals for colour display, where a range of compounds and structures can produce a diverse colour palette. Feather colours, for example, span the visible spectrum and mostly result from pigments in five chemical classes (carotenoids, melanins, porphyrins, psittacofulvins and metal oxides). However, the pigment that generates the yellow colour of penguin feathers appears to represent a sixth, poorly characterized class of feather pigments. This pigment class, here termed 'spheniscin', is displayed by half of the living penguin genera; the larger and richer colour displays of the pigment are highly attractive. Using Raman and mid-infrared spectroscopies, we analysed yellow feathers from two penguin species (king penguin, Aptenodytes patagonicus; macaroni penguin, Eudyptes chrysolophus) to further characterize spheniscin pigments. The Raman spectrum of spheniscin is distinct from spectra of other feather pigments and exhibits 17 distinctive spectral bands between 300 and 1700 cm(-1). Spectral bands from the yellow pigment are assigned to aromatically bound carbon atoms, and to skeletal modes in an aromatic, heterocyclic ring. It has been suggested that the penguin pigment is a pterin compound; Raman spectra from yellow penguin feathers are broadly consistent with previously reported pterin spectra, although we have not matched it to any known compound. Raman spectroscopy can provide a rapid and non-destructive method for surveying the distribution of different classes of feather pigments in the avian family tree, and for correlating the chemistry of spheniscin with compounds analysed elsewhere. We suggest that the sixth class of feather pigments may have evolved in a stem-lineage penguin and endowed modern penguins with a costly plumage trait that appears to be chemically unique among birds.

Quantitative investigations of quantum coherence for a light-harvesting protein at conditions simulating photosynthesis

Recent measurements using two-dimensional electronic spectroscopy (2D ES) have shown that the initial dynamic response of photosynthetic proteins can involve quantum coherence. We show how electronic coherence can be differentiated from vibrational coherence in 2D ES. On that basis we conclude that both electronic and vibrational coherences are observed in the phycobiliprotein light-harvesting complex PC645 from Chroomonas sp. CCMP270 at ambient temperature. These light-harvesting antenna proteins of the cryptophyte algae are suspended in the lumen, where the pH drops significantly under sustained illumination by sunlight. Here we measured 2D ES of PC645 at increasing levels of acidity to determine if the change in pH affects the quantum coherence; quantitative analysis reveals that the dynamics are insensitive to the pH change.

Structure of the chromophore binding pocket in the Pr state of plant phytochrome phyA

A homology structural model was generated for plant phytochrome phyA utilizing the crystal structure of the sensory module of cyanobacterial phytochrome Cph1 (Cph1Δ2). As chromophores, either the native phytochromobilin cofactor (PΦB) or phycocyanobilin (PCB), the natural cofactor in Cph1, was incorporated. These homology models were further optimized by molecular dynamics (MD) simulations revealing a satisfying overall agreement with the crystal structure of Cph1Δ2. Notable differences in the PΦB adduct of phyA result from a restructuring of the small helical segment α(7) that leads to displacements of a few amino acids away from the cofactor. This repositioning of residues also include aspartate 218 such that, instead of its carbonyl function as in Cph1Δ2, an additional water molecule forms hydrogen bonds with the ring B and C NH groups. To validate the phyA structural model in the chromophore binding pocket, Raman spectra of the cofactor were calculated by means of the quantum mechanics/molecular mechanics (QM/MM) hybrid methodology and compared with the experimental resonance Raman (RR) spectra. The satisfactory overall agreement between calculated and experimental spectra is taken as an indication for the good quality of the structural model. Moreover, the methine bridge stretching modes and the effects of isotopic labeling at selected positions of the chromophore are very well reproduced to allow confirming even details of the methine bridge geometry as predicted by the homology model. Specifically, it is demonstrated that the experimental RR spectra are consistent with a torsional angle of ring D with respect to ring C that is distinctly higher for phyA-PCB (45°) and phyA-PΦB (42°) than for Cph1Δ2 (30°). Raman spectra calculated from different points of the MD trajectory display variations of the mode frequencies and intensities reflecting the structural fluctuations from snapshot to snapshot. The snapshot spectrum of the lowest energy structure and the sum of all snapshot spectra afford an equally good description of the experimental data. Particularly large variations between the snapshots are noted for the N-H in-plane bending mode of the pyrrole rings B and C, which reflect alterations of the hydrogen bond interactions brought about by fluctuations of water molecules in the cofactor cavity. This overestimation of the water molecule mobility is a consequence of the deficiency of the current QM/MM methodology that, due to the lack of appropriate protein force fields, cannot adequately account for the electrostatics in the cofactor pocket.

Chromophore structure of cyanobacterial phytochrome Cph1 in the Pr state: reconciling structural and spectroscopic data by QM/MM calculations

A quantum mechanics (QM)/molecular mechanics (MM) hybrid method was applied to the Pr state of the cyanobacterial phytochrome Cph1 to calculate the Raman spectra of the bound PCB cofactor. Two QM/MM models were derived from the atomic coordinates of the crystal structure. The models differed in the protonation site of His(260) in the chromophore-binding pocket such that either the delta-nitrogen (M-HSD) or the epsilon-nitrogen (M-HSE) carried a hydrogen. The optimized structures of the two models display small differences specifically in the orientation of His(260) with respect to the PCB cofactor and the hydrogen bond network at the cofactor-binding site. For both models, the calculated Raman spectra of the cofactor reveal a good overall agreement with the experimental resonance Raman (RR) spectra obtained from Cph1 in the crystalline state and in solution, including Cph1 adducts with isotopically labeled PCB. However, a distinctly better reproduction of important details in the experimental spectra is provided by the M-HSD model, which therefore may represent an improved structure of the cofactor site. Thus, QM/MM calculations of chromoproteins may allow for refining crystal structure models in the chromophore-binding pocket guided by the comparison with experimental RR spectra. Analysis of the calculated and experimental spectra also allowed us to identify and assign the modes that sensitively respond to chromophore-protein interactions. The most pronounced effect was noted for the stretching mode of the methine bridge A-B adjacent to the covalent attachment site of PCB. Due a distinct narrowing of the A-B methine bridge bond angle, this mode undergoes a large frequency upshift as compared with the spectrum obtained by QM calculations for the chromophore in vacuo. This protein-induced distortion of the PCB geometry is the main origin of a previous erroneous interpretation of the RR spectra based on QM calculations of the isolated cofactor.

Analysis and prediction of absorption band shapes, fluorescence band shapes, resonance Raman intensities, and excitation profiles using the time-dependent theory of electronic spectroscopy

A general method for the simulation of absorption (ABS) and fluorescence band shapes, resonance-Raman (rR) spectra, and excitation profiles based on the time-dependent theory of Heller is discussed. The following improvements to Heller's theory have been made: (a) derivation of new recurrence relations for the time-dependent wave packet overlap in the case of frequency changes between the ground and electronically excited states, (b) a new series expansion that gives insight into the nature of Savin's preresonance approximation, (c) incorporation of inhomogeneous broadening effects into the formalism at no additional computational cost, and (d) derivation of a new and simple short-time dynamics based equation for the Stokes shift that remains valid in the case of partially resolved vibrational structure. Our implementation of the time-dependent theory for the fitting of experimental spectra and the simulation of model spectra as well as the quantum mechanical calculation of the model parameters is discussed. The implementation covers all electronic structure approaches which are able to deliver ground- and excited-state energies and transition dipole moments. The technique becomes highly efficient if analytic gradients for the excited-state surface are available. In this case, the computational cost for the simultaneous prediction of ABS, fluorescence, and rR spectra is equal to that of a single excited-state geometry optimization step while the limitations of the short-time dynamics approximation are completely avoided. As a test case we discuss the well-known case of the strongly allowed 1 (1)A(g) --> 1 (1)B(u) transition in 1,3,5 trans-hexatriene in detail using method ranging from simple single-reference treatments to elaborate multireference electronic structure approaches. At the highest computational level, the computed spectra show the best agreement that has so far been obtained with quantum chemical methods for this problem.

Assignments of the Pfr-Pr FTIR difference spectrum of cyanobacterial phytochrome Cph1 using 15N and 13C isotopically labeled phycocyanobilin chromophore

The reversible red and far-red light-induced transitions of cyanobacterial phytochrome Cph1 from Synechocystis PCC 6803 were investigated by Fourier transform infrared (FTIR) difference spectroscopy. High-quality light-induced Pfr-Pr difference FTIR spectra were recorded for the 58 kDa N-terminal domain of Cph1 by repetitive photochemical cycling and signal averaging. The Pfr-Pr difference spectra in H(2)O and D(2)O were very similar to those previously reported for full-length 85 kDa Cph1.(1) Published assignments were extended by analysis of the effects of (13)C and (15)N isotope substitutions at selected sites in the phycocyanobilin chromophore and by (15)N global labeling of the protein. The Pfr-Pr difference spectra were dominated by an amide I peak/trough at 1653 cm(-1)(+)/1631 cm(-1)(-) and a smaller amide II band at 1554 cm(-1). Labeling effects allowed specific chromophore assignments for the C(1)=O (1736 cm(-1)(-)/1724 cm(-1)(+)) and C(19)=O (1704 cm(-1)(-)) carbonyl vibrations, C=C vibrations at 1589 cm(-1)(+), and bands at 1537(-), 1512(+), 1491(-), 1163(+), 1151(-), 1134(+), 1109(-), and 1072(-) cm(-1) that must involve chromophore C-N bonds. A variety of additional changes were insensitive to isotope labeling of the chromophore. Effects of (15)N labeling of the protein were used to tentatively assign some of these to specific amino acid changes. Those insensitive to (15)N labeling included a protonated aspartic or glutamic acid at 1734 cm(-1)(-)/1722 cm(-1)(+) and a cysteine at 2575 cm(-1)(+)/2557 cm(-1)(-). Bands sensitive to (15)N protein labeling at 1487 cm(-1)(+)/1502 cm(-1)(-) might arise from trytophan and bands at 1261 cm(-1)(+)/1244 cm(-1)(-) and 1107 cm(-1)(-)/1095 cm(-1)(+) might arise from a histidine environment or protonation change. These assignments are discussed in light of the 15Z-E photoisomerization model of phototransformation and the associated protein conformational changes.

Calculation of vibrational spectra of linear tetrapyrroles. 3. Hydrogen-bonded hexamethylpyrromethene dimers

The structure and vibrational spectra of hexamethylpyrromethene (HMPM) have been investigated by X-ray crystallography, IR and Raman spectroscopies, and density functional theory calculations. HMPM crystallizes in the form of dimers, which are held together by bifurcated N-H(...N)(2) hydrogen bonds, involving one intramolecular and one intermolecular N-H...N interaction. The monomers are essentially planar, and the mean planes of the monomers lie approximately perpendicular to one another, so that the four N atoms in the dimer form a distorted tetrahedron. The structure of the HMPM dimer is well-reproduced by B3LYP/6-31G calculations. A comparison of the calculated geometry of the dimer with that of the monomer reveals only small changes in the N-H...N entity and the methine bridge angles upon dimerization. These are a result of weakening of the intramolecular N-H...N hydrogen bond and the formation of a more linear N-H...N intermolecular hydrogen bond. Using an empirical relation between the shift of the N-H stretching frequency of pyrrole and the enthalpy of adduct formation with bases [Nozari, M. S.; Drago, R. S. J. Am. Chem. Soc. 1970, 92, 7086-7090], estimates of the strength of the intra- and intermolecular hydrogen bonds are obtained. IR and Raman spectroscopies of HMPM and its isotopomers deuterated at the pyrrolic nitrogen atom and at the methine bridge reveal that the molecule is monomeric in nonpolar organic solvents but dimeric in a solid Ar matrix and in KBr pellets. The matrix IR spectra show a splitting of vibrational modes for the dimer, particularly those involving the N-H coordinates. Due to intrinsic deficiencies of the B3LYP/6-31G approximation, a satisfactory reproduction of these modes of the monomeric and dimeric HMPM requires specific adjustments of the NH scaling factors for the calculated force constants and, in the case of the NH out-of-plane modes of HMPM dimers, also of intra- and intermolecular coupling constants. This parametrization does not significantly affect the other calculated modes, which in general reveal a very good agreement with the experimental data.

Quantum mechanics/molecular mechanics calculation of the Raman spectra of the phycocyanobilin chromophore in alpha-C-phycocyanin

We have established a quantum mechanics (QM)/molecular mechanics (MM) hybrid method for calculating the Raman spectra of protein-bound cofactors using the alpha-subunit of C-phycocyanin containing a phycocyanobilin (PCB) chromophore as a test case. The PCB cofactor was described with density functional theory, whereas the protein matrix was treated with the CHARMM force field. The Hessian matrix of the QM region was built by taking into account bonded and nonbonded interactions with the protein environment and projected onto the internal coordinate space. Force constants were scaled with a global set of scaling factors, and the Raman intensities were computed using a finite-field method combined with a fourth-order differentiation algorithm for the calculation of the polarizability derivatives. In general, the QM/MM results provided a substantially improved description of the experimental resonance Raman (RR) spectra of the protein-bound cofactor compared to QM calculations of isolated PCB models in vacuo. The results allow the assessment of the effect of the protein-cofactor interactions on the RR spectra and reveal the potential and limitations of QM calculations on isolated tetrapyrroles for determining the chromophore structures in the various species and states of phytochromes for which three-dimensional structures are not available.

The chromophore structural changes during the photocycle of phytochrome: a combined resonance Raman and quantum chemical approach

Phytochromes are sensory photoreceptors that, upon light irradiation, can be transformed between an inactive and an active state. The conversion is initiated by the photoisomerization of the cofactor, a linear methine-bridged tetrapyrrole, followed by conformational relaxations of the chromophore and the protein matrix that finally leads to the formation of the signaling state. To elucidate the underlying molecular processes, resonance Raman spectroscopy combined with quantum chemical calculations constitutes a powerful approach since it allows determination of the chromophore structure in the various states of phytochrome. On the basis of these studies, a molecular model for the photoinduced reaction cycle is derived.

Calculation of Vibrational Spectra of Linear Tetrapyrroles. 4. Methine Bridge C−H Out-of-Plane Modes

Quantum chemical force fields obtained by density functional theory (DFT) calculations systematically overestimate the frequencies of normal modes including ethylenic C-H out-of-plane (HOOP) coordinates. Compensation of this deviation requires a specific scaling factor for this type of coordinate that is distinctly lower than those applicable to out-of-plane coordinates in general. Such a specific scaling factor (0.900) has been optimized for the DFT(B3LYP) level of theory on the basis of vibrational analyses of training molecules including the HOOP coordinate. Thus, the root-mean-square deviation for the calculated frequencies of these modes is reduced from 16 to 8 cm(-1). Although Raman intensities are yet not reproduced in a satisfactory manner, implementation of the HOOP scaling factor into the set of global scaling factors determined previously (Magdo et al. J. Phys. Chem. A 1999, 103, 289-303) allows for a substantially improved reproduction of the experimental (resonance) Raman spectra of test molecules including linear methine-bridged tetrapyrroles. A very good agreement between calculated and experimental spectra is noted for the phycocyanobilin dimethylester dimer as well as for the protein-bound phycocyanobilin in the antenna pigment alpha-CPC. However, for the phycocyanobilin chromophore in the P(r) state of the plant photoreceptor phytochrome phyA, considerable deviations remain in the spectral range between 800 and 500 cm(-1), which are attributed to the effect of specific protein-chromophore interactions. The influence of the protein environment is not considered in the present calculations that refer to the molecule in vacuo.

Proton transfer in the photoreceptors phytochrome and photoactive yellow protein

Light-induced activation of the photoreceptors phytochrome and photoactive yellow protein (PYP) is accompanied by protonation changes of the respective chromophores and key residues in the protein moiety. For both systems, proton exchange with the external medium could be observed with pH electrode measurements and with UV-visible absorption spectroscopy using appropriate pH indicator dyes. From these signals, the stoichiometry of proton release and uptake, respectively, was determined by different calibration procedures which will be presented and discussed. Kinetic information on these processes is only available from time-resolved measurements with pH indicator dyes. Vibrational spectroscopy methods such as Fourier transform infrared spectroscopy and resonance Raman scattering provided information on the protonation state of individual functional groups suggesting that internal proton transfer processes are involved as well. Deuterium kinetic isotope effects that occurred in the Pr --> Pfr phototransformation of the bacteriophytochromes Cph1 and Agp1 were consistent with proton transfer reactions as rate-limiting steps. In contrast, the apparent rate constants in the photocycle of PYP exhibited only small kinetic isotope effects that could not be interpreted conclusively. Possible mechanisms of proton transfer in the activation of phytochrome and PYP will be discussed.

Determination of the Chromophore Structures in the Photoinduced Reaction Cycle of Phytochrome

The chromophore structures in the parent states Pr and Pfr as well as in the photocycle intermediate Lumi-R of oat phytochrome phyA are determined by comparison of the experimental resonance Raman spectra with calculated Raman spectra that have been obtained by density functional theory calculations (B3LYP) using scaled force fields. The spectra were calculated for various tetrapyrrole geometries including more than twenty different methine bridge isomers. For the parent states Pr and Pfr the best agreement in terms of vibrational frequencies, isotopic shifts, and Raman intensities was achieved with the ZZZasa and ZZEssa geometry, respectively. For the first intermediate Lumi-R, the chromophore geometry is concluded to be the ZZEasa configuration. These finding imply that the primary step of the photoactivation of phytochrome is the Z/E isomerization of the C-D methine bridge double bond, whereas the single bond remains in the anti conformation. The subsequent transition to the physiologically active state Pfr includes a (partial) single bond rotation of the A-B methine bridge.

Resonance Raman spectroscopic study of the tryptic 39-kDa fragment of phytochrome

The 39-kDa fragment of oat phytochrome phyA, obtained by tryptic digestion at the amino acids 65 and 425, was studied by resonance Raman spectroscopy. The parent state P(r) reveals far-reaching similarities with that of the native phytochrome implying that the structures of the tetrapyrrole chromophore and its immediate protein environment are not affected by the proteolysis. However, the resonance Raman spectrum of the final product of the P(r) phototransformation, denoted as P(bl), is more closely related to that of the P(fr) precursor of the native phytochrome, i.e. meta-R(C), rather than to that of P(fr) itself. The resonance Raman spectra indicate a high conformational flexibility of the chromophore in P(bl) so that, unlike in P(fr), the tetrapyrrole rings C and D adopt a largely coplanar conformation. The protein interactions with ring D of the chromophore, which in the native phytochrome stabilize the specific chromophore structure of P(fr), cannot be established in the 39-kDa fragment due to the lack of the major C-terminal part of the protein. These findings, furthermore, support the view that the meta-R(C)-->P(fr) transition is associated with a coupling of chromophore and protein structural changes that represent crucial events for the photoactivation of phytochrome.

Vibrational spectroscopy as a tool for probing protein function

Vibrational spectroscopy has become increasingly important as a tool for understanding the mechanisms of photosystem II, phytochrome and terminal oxidases. More general enzymatic or receptor systems have been studied, opening a new field of applications. Femtosecond infrared pump/probe studies of the important amide-I band seem to provide a basis for its molecular and structural interpretation.

Protonation state and structural changes of the tetrapyrrole chromophore during the Pr --> Pfr phototransformation of phytochrome: a resonance Raman spectroscopic study

The photoconversion of phytochrome (phytochrome A from Avena satina) from the inactive (Pr) to the physiologically active form (Pfr) was studied by near-infrared Fourier transform resonance Raman spectroscopy at cryogenic temperatures, which allow us to trap the intermediate states. Nondeuterated and deuterated buffer solutions were used to determine the effect of H/D exchange on the resonance Raman spectra. For the first time, reliable spectra of the "bleached" intermediates meta-R(A) and meta-R(C) were obtained. The vibrational bands in the region 1300-1700 cm(-)(1), which is particularly indicative of structural changes in tetrapyrroles, were assigned on the basis of recent calculations of the Raman spectra of the chromophore in C-phycocyanin and model compounds [Kneip, C., Hildebrandt, P., Németh, K., Mark, F., Schaffner, K. (1999) Chem. Phys. Lett. 311, 479-485]. The experimental resonance Raman spectra Pr are compatible with the Raman spectra calculated for the protonated ZZZasa configuration, which hence is suggested to be the chromophore structure in this parent state of phytochrome. Furthermore, marker bands could be identified that are of high diagnostic value for monitoring structural changes in individual parts of the chromophore. Specifically, it could be shown that not only in the parent states Pr and Pfr but also in all intermediates the chromophore is protonated at the pyrroleninic nitrogen. The spectral changes observed for lumi-R confirm the view that the photoreaction of Pr is a Z --> E isomerization of the CD methine bridge. The subsequent thermal decay reaction to meta-R(A) includes relaxations of the CD methine bridge double bond, whereas the formation of meta-R(C) is accompanied by structural adaptations of the pyrrole rings B and C in the protein pocket. The far-reaching similarities between the chromophores of meta-R(A) and Pfr suggest that in the step meta-R(A) --> Pfr the ultimate structural changes of the protein matrix occur.