Calculated vibrational properties of pigments in protein binding sites

FTIR difference spectroscopy is widely used to probe molecular bonding interactions of protein-bound electron transfer cofactors. The technique is particularly attractive because it provides information on both neutral and radical cofactor states. Such dual information is not easily obtainable using other techniques. Although FTIR difference spectroscopy has been used to study cofactors in biological protein complexes, in nearly all cases interpretation of the spectra has been purely qualitative. Virtually no computational work has been undertaken in an attempt to model the spectra. To address this problem we have developed the use of ONIOM (our own N-layered integrated molecular Orbital + Molecular mechanics package) (quantum mechanical:molecular mechanics) methods to calculate FTIR difference spectra associated with protein-bound cofactors. As a specific example showing the utility of the approach we have calculated isotope edited FTIR difference spectra associated with unlabeled and labeled ubiquinones in the Q(A) binding site in Rhodobacter sphaeroides photosynthetic reaction centers. The calculated spectra are in remarkable agreement with experiment. Such agreement cannot be obtained by considering ubiquinone molecules in the gas phase or in solution. A calculation including the protein environment is required. The ONIOM calculated spectra agree well with experiment but indicate a very different interpretation of the experimental data compared to that proposed previously. In particular the calculations do not predict that one of the carbonyl groups of Q(A) is very strongly hydrogen bonded. We show that a computational-based interpretation of FTIR difference spectra associated with protein-bound cofactors is now possible. This approach will be applicable to FTIR studies of many cofactor-containing proteins.