This Article Full Text Full Text (PDF) Supplement All Versions of this Article: biophysj.107.108878v1 93/6/1885    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Mroginski, M. A. Articles by Hildebrandt, P. Search for Related Content PubMed PubMed Citation Articles by Mroginski, M. A. Articles by Hildebrandt, P.

Quantum Mechanics/Molecular Mechanics Calculation of the Raman Spectra of the Phycocyanobilin Chromophore in -C-Phycocyanin

Maria Andrea Mroginski ^*, Franz Mark , Walter Thiel  and Peter Hildebrandt ^*

^* Technische Universität Berlin, Institut für Chemie, Max-Volmer-Laboratorium für Biophysikalische Chemie, Sekr. PC 14, D-10623 Berlin, Germany; Max-Planck-Institut für Bioanorganische Chemie, D-45470 Mülheim, Germany; and Max-Planck-Institut für Kohlenforschung, D-45470 Mülheim, Germany

Correspondence: Address reprint requests to Maria Andrea Mroginski, Tel.: 49-30-314-21584; Fax: 49-30-314-21122; E-mail: andrea.mroginski{at}tu-berlin.de.

We have established a quantum mechanics (QM)/molecular mechanics (MM) hybrid method for calculating the Raman spectra of protein-bound cofactors using the -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.