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

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.

NH Stretching Frequencies of Intramolecularly Hydrogen-Bonded Systems: An Experimental and Theoretical Study

The vibrational NH stretching transitions in secondary amines with intramolecular NH···O hydrogen bonds were investigated by experimental and theoretical methods, considering a large number of compounds and covering a wide range of stretching wavenumbers. The assignment of the NH stretching transitions in the experimental IR spectra was, in several instances, supported by measurement of the corresponding ND wavenumbers and by correlation with the observed NH proton chemical shifts. The observed wavenumbers were correlated with theoretical wavenumbers predicted with B3LYP density functional theory, using the basis sets 6-311++G(d,p) and 6-31G(d) and considering the harmonic as well as the anharmonic VPT2 approximation. Excellent correlations were established between observed wavenumbers and calculated harmonic values. However, the correlations were non-linear, in contrast to the results of previous investigations of the corresponding OH···O systems. The anharmonic VPT2 wavenumbers were found to be linearly related to the corresponding harmonic values. The results provide correlation equations for the prediction of NH stretching bands on the basis of standard B3LYP/6-311++G(d,p) and B3LYP/6-31G(d) harmonic analyses, with standard deviations close to 38 cm⁻¹. This is significant because the full anharmonic VPT2 analysis tends to be impractical for large molecules, requiring orders of magnitude more computing time than the harmonic analysis.

Ultrafast proton release reaction and primary photochemistry of phycocyanobilin in solution observed with fs-time-resolved mid-IR and UV/Vis spectroscopy

Deactivation processes of photoexcited (λ_ex = 580 nm) phycocyanobilin (PCB) in methanol were investigated by means of UV/Vis and mid-IR femtosecond (fs) transient absorption (TA) as well as static fluorescence spectroscopy, supported by density-functional-theory calculations of three relevant ground state conformers, PCB_A, PCB_B and PCB_C, their relative electronic state energies and normal mode vibrational analysis. UV/Vis fs-TA reveals time constants of 2.0, 18 and 67 ps, describing decay of PCB_B*, of PCB_A* and thermal re-equilibration of PCB_A, PCB_B and PCB_C, respectively, in line with the model by Dietzek et al. (Chem Phys Lett 515:163, 2011) and predecessors. Significant substantiation and extension of this model is achieved first via mid-IR fs-TA, i.e. identification of molecular structures and their dynamics, with time constants of 2.6, 21 and 40 ps, respectively. Second, transient IR continuum absorption (CA) is observed in the region above 1755 cm^-1 (CA1) and between 1550 and 1450 cm^-1 (CA2), indicative for the IR absorption of highly polarizable protons in hydrogen bonding networks (X-H^…Y). This allows to characterize chromophore protonation/deprotonation processes, associated with the electronic and structural dynamics, on a molecular level. The PCB photocycle is suggested to be closed via a long living (> 1 ns), PCB_C-like (i.e. deprotonated), fluorescent species.

Vibrational spectra of solid cis- and trans-2-thioxohexahydroquinazolin-4(1H)-one and theoretical calculations towards the interpretation of its thermal reactivity

FT-Raman and FT-IR spectra ofcisandtrans2-thioxohexahydroquinazolin-4(1H)-one are reported. Both compounds are dimers in the solid phase, withC₂symmetry. This work contributes to the knowledge of data which are rather scarce for quinazolinones.

Comparison of resonance assisted and charge assisted effects in strengthening of hydrogen bonds in dipyrrins

This paper deals with the study of two types of hydrogen bonding: a quasi-aromatic hydrogen bonding in dipyrromethene and the ionic one in dipyrromethane. The study focuses on two phenomena-the proton transfer process and tautomeric equilibrium. Metric parameters and spectroscopic assignments have been calculated; this allowed a further comparison of spectral features calculated with four methods (Car-Parrinello molecular dynamics (CPMD), ab initio, density functional theory (DFT), and numerical calculation of anharmonic vibrational levels via a solution of the corresponding 1D Schrödinger equation). A significant dynamics of the bridged proton and bent vibration of pyrrole fragments in dipyrromethane have been exposed by the CPMD calculations. The prevailing of the ionic effect over the π-electronic coupling in the strengthening of the hydrogen bonding has been shown on the basis of the calculated structural, electron-topological, and spectral data as well as potential energy surface (PES). The analysis of the aromaticity and electronic state of pyrrole and chelate moieties depending on the tautomeric equilibrium by the quantum theory of atoms in molecules (QTAIM) method was conducted. The principle divergence in the behavior of aromaticity of the chelate chains in the analyzed compounds was demonstrated.

Refined syntheses of hydrodipyrrin precursors to chlorin and bacteriochlorin building blocks

Bromo-substituted hydrodipyrrins are valuable precursors to synthetic bromo-chlorins and bromo-bacteriochlorins, which in turn are versatile substrates for derivatization in pursuit of diverse molecular designs. 8-bromo-2,3-dihydro-1-(1,1-dimethoxymethyl)-3,3-dimethyldipyrrin (1) is a crucial precursor in the rational synthesis of the bacteriochlorin building block 3,13-dibromo-8,8,18,18-tetramethylbacteriochlorin ( H2BC-Br3Br13) . 8-bromo-2,3,4,5-tetrahydro-1,3,3-trimethyldipyrrin (2) is a crucial precursor in the rational synthesis of the analogous 3,13-disubstituted chlorin building block (e.g. H2C-Br3M10Br13 ). The routes to 1 and 2 share a common precursor, namely 4-bromo-2-(2-nitroethyl)-1-N-tosylpyrrole (6-Ts), which is derived from pyrrole-2-carboxaldehyde. The prior seven-step synthesis of 1 from pyrrole-2-carboxaldehyde has limited access to H2BC-Br3Br13 given the large excesses of materials, extensive reliance on column chromatography, and low overall yield (1.4%). Refined procedures for synthesis of the common precursor 6-Ts as well as 1 and 2 afford the advantages of (1) diminished consumption of solvents and reagents, (2) limited or no use of chlorinated solvents, (3) limited or no chromatography, and (4) improved yields of most steps. Streamlined procedures enable the final two or three transformations to be performed without purification of intermediates. The new procedures facilitate the expedient preparation of 1 and 2 at the multigram scale.

Light-induced activation of bacterial phytochrome Agp1 monitored by static and time-resolved FTIR spectroscopy

Phytochromes, which regulate many biological processes in plants, bacteria, and fungi, can exist in two stable states, Pr and Pfr, that can be interconverted by light, via a number of intermediates such as meta-Rc. Herein we employ FTIR spectroscopy to study the Pr-to-Pfr conversion of the bacteriophytochrome Agp1 from Agrobacterium tumefaciens. Static FTIR Pfr/Pr and meta-Rc/Pr difference spectra are disentangled in terms of cofactor and protein structural changes. Guided by DFT calculations on cofactor models, the chromophore conformational changes can be grouped into structural adjustments of the cofactor-protein interactions localized in the C-D dipyrrole moiety, that is, the photoisomerisation site, and in the A-B dipyrrole moiety including the protein attachment site. Whereas changes at the C and D rings appear to be largely completed in the meta-Rc state, the structural changes in the A-B unit occur during the transition from meta-Rc to Pfr, concomitant with the main protein structural changes, as demonstrated by static and time-resolved FTIR difference spectroscopy. We employ this technique to monitor, for the first time, the dynamics of the photocycle of phytochrome on the millisecond timescale. By extending the studies to genetically engineered protein variants of Agp1, we further demonstrate that H250 and D197 as well as the PHY domain are essential for formation of the Pfr state. Based on the IR spectroscopic and available crystallographic data we discuss the role of critical amino acid residues for the protein-cofactor interactions during the photoinduced reaction cycle.

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.

The chromophore structures of the Pr States in plant and bacterial phytochromes

The resonance Raman spectra of the Pr state of the N-terminal 65-kDa fragment of plant phytochrome phyA have been measured and analyzed in terms of the configuration and conformation of the tetrapyrroles methine bridges. Spectra were obtained from phyA adducts reconstituted with the natural chromophore phytochromobilin as well as phycocyanobilin and its isotopomers labeled at the terminal methine bridges through (13)C/(12)C and D/H substitution. Upon comparing the resonance Raman spectra of the various phyA adducts, it was possible to identify the bands that originate from normal modes dominated by the stretching coordinates of the terminal methine bridges A-B and C-D. Quantum chemical calculations of the isolated tetrapyrroles reveal that these modes are sensitive indicators for the methine bridge configuration and conformation. For all phyA adducts, the experimental spectra of Pr including this marker band region are well reproduced by the calculated spectra obtained for the ZZZasa configuration. In contrast, there are substantial discrepancies between the experimental spectra and the spectra calculated for the ZZZssa configuration, which has been previously shown to be the chromophore geometry in the Pr state of the bacterial, biliverdin-binding phytochrome from Deinococcus radiodurans (Wagner, J. R., J. S. Brunzelle, K. T. Forest, R. D. Vierstra. 2005. Nature. 438:325-331). The results of this work, therefore, suggest that plant and bacterial (biliverdin-binding) phytochromes exhibit different structures in the parent state although the mechanism of the photoinduced reaction cycle may be quite similar.

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.