Computational Biology Journal

Volume 2013 (2013), Article ID 807592, 11 pages

http://dx.doi.org/10.1155/2013/807592

## Calculated Vibrational Properties of Ubisemiquinones

Department of Physics and Astronomy, Georgia State University, 29 Peachtree Center Avenue, Atlanta, GA 30303, USA

Received 15 October 2012; Accepted 27 November 2012

Academic Editor: Philip Crooke

Copyright © 2013 Hari P. Lamichhane and Gary Hastings. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Density functional theory has been used to calculate harmonic normal mode vibrational frequencies for unlabeled and isotope-labeled ubisemiquinones in both the gas phase and in several solvents. It is shown that four methoxy group conformations are likely to be present in solution at room temperature. Boltzmann weighted infrared and Raman spectra for the four conformers were calculated, and composite spectra that are the sum of the Boltzmann weighted spectra were produced. These composite spectra were compared to experimental FTIR and resonance Raman spectra, and it is shown that the calculated band frequencies, relative band intensities, and and isotope-induced band shifts are in excellent agreement with experiment. The calculations show that the C=O and C=C modes of ubisemiquinone strongly mix with methoxy methyl CH bending vibrations, and that the degree of mixing is altered upon isotope labeling, resulting in complicated changes in mode frequencies, intensities, and composition upon isotope labeling. Upon consideration of the calculated potential energy distributions of the normal modes of ubisemiquinone, and how they change upon isotope labeling, an explanation of some puzzling features in previously published Raman spectra is provided.

#### 1. Introduction

Ubiquinones (: 2,3-dimethoxy-5-methyl-6-polyprenyl-1,4-benzoquinones) play an important role in biological electron and proton transfer processes that occur in both respiration and photosynthesis [1]. In photosynthetic reaction centers from purple bacteria, two UQ molecules, called and , act as terminal electron acceptors [2]. In purple bacterial reaction centers (PBRCs) (see Abbreviations) from *Rhodobacter *(*Rb.*)* sphaeroides*, and are both ubiquinone-10 (UQ_{10}) molecules. and have very different functions; however, is an intermediary cofactor involved in transferring electrons from bacteriopheophytin to , while couples electron and proton transfer processes [3, 4]. The very different redox functions of and are testimony to the flexibility of UQs in biological processes. Since and are both UQ_{10} molecules, pigment-protein interactions must modulate the functional properties of UQ_{10} in PBRCs. Elucidation of these pigment-protein interactions is at the heart of much current research in photosynthesis [5, 6].

Fourier transform infrared (FTIR) difference spectroscopy (DS) is a sensitive molecular-level probe of pigment-protein interactions, and it is widely used to study both the neutral and reduced states of the quinones occupying the and binding sites in PBRCs [7]. Although and FTIR difference spectra have been obtained under a wide range of conditions for variously treated PBRC’s, these spectra continue to be difficult to interpret because many bands not associated with the quinone also contribute to the spectra. Reconstitution of PBRCs with isotopically labeled quinones, however, has allowed some separation of the contributions of the quinones from those of the protein to the spectra [7]. Nonetheless the hypothesized band assignments in the experimental spectra, particularly those assignments associated with the ubiquinone anion radical, are still ambiguous and have not been modeled computationally.

One basis for developing an understanding of bands in and FTIR DS is to first consider spectra of the relevant quinones in solution. Infrared (IR) absorption spectra [8, 9] and resonance Raman spectra [10] for ubisemiquinones in solution have been obtained. However, from a computational standpoint, even these simpler solution spectra are poorly understood. The work outlined in this paper is aimed at addressing this problem.

Few computational studies aimed at modeling the vibrational properties of ubisemiquinones () have been undertaken. The work that has been undertaken [11, 12] is limited in one way or another; for example, tail-less quinone models in only the gas phase were considered, using relatively low levels of theory. Previously it was claimed that the calculated normal modes and associated isotope-induced frequency shifts are in good agreement with experiment [12]. Isotope shifts do appear to agree with experiment. However, upon careful examination, it appears that the previously calculated normal modes (frequencies and intensities) are *not* in agreement with experimental spectra (see below). In the light of this finding we have used more robust computational methods to investigate the vibrational properties of ubisemiquinones in the gas phase and in solution.

In this paper we describe the simulation of FTIR and Raman spectra associated with labeled and unlabeled tail-containing ubisemiquinones in both the gas phase and in solvent.

#### 2. Materials and Methods

##### 2.1. Calculations

Molecular geometry optimizations and harmonic vibrational frequency calculations were performed using hybrid density functional theory (DFT) methods, employing the B3LYP functional and the basis set within Gaussian 03 [13]. is preferable to for calculations involving semiquinones [14]. For calculations including solvent, the integral equation formalism (IEF) [15–17] of the polarizable continuum model (PCM) [18, 19] was used. The PCM uses the united atom cavity approach. Cavity parameters used were (overlap index between interlocking spheres) and (minimum radius in Angstroms for overlapping spheres). Very similar spectra were calculated when a smaller number of added spheres were considered ( and ). The potential energy distribution (PED) (or total energy distribution) of normal modes was calculated using gar2ped [20].

Calculated normal mode vibrational frequencies presented here were scaled by 0.9808. Such a scale factor is standard for calculations using the specified functional and basis set and was derived by comparing the frequencies of bands in experimental and calculated spectra. Such a scaling of the calculated frequencies is undertaken only to facilitate a comparison between calculated and experimental spectra. We are primarily interested in vibrational frequency changes that occur upon isotope labeling, and these frequency differences are accurately calculated without scaling [14, 21].

#### 3. Results

##### 3.1. UQ Structure and Numbering

Figure 1 shows a geometry-optimized model with the atom numbering scheme displayed. UQ has two carbonyl groups (C_{1}=O_{18} and C_{4}=O_{15}), two methoxy groups (C_{3}–O_{16}–CH_{3} and C_{2}–O_{17}–CH_{3}), a methyl group at C_{5}, and an isoprene unit at C_{6}. In our calculations we used UQ models with only a single isoprene unit. As outlined previously [21], the calculated vibrational properties of UQ_{1} (or ) are very similar to that found for (with ). Also shown in Figure 1 are relevant internal coordinates of . The normal modes will be expressed in terms of contributions from these internal coordinates. Of particular interest in this paper are the coordinates R3, R9, R4, and R10 which are due to , , , and stretching vibrations, respectively. The methoxy methyl CH bending vibrations (coordinates *δ*C8 and *δ*C9) are also of considerable relevance in this paper, as they strongly couple to the vibrations (see below). This was not found to occur for neutral UQ [21].

##### 3.2. Calculated Structure of Ubisemiquinone_{1} ()

Previously we showed that neutral UQ_{1} can adopt at least eight different methoxy group conformations at room temperature [21]. To establish which conformations may be present for , single-point energy calculations were undertaken for methoxy group dihedral angles that were stepped in 10° increments. That is, structures with fixed methoxy group dihedral angles were geometry optimized.

A contour plot of energy versus the C_{2} and C_{3} dihedral angles is shown in Figure 2, which indicates that there are four low-energy conformations, each with C_{2} and C_{3} dihedral angles close to ±120°. The four conformers are labeled A, B, E, and F in Figure 2. These four conformers are similar to the neutral UQ_{1} conformers labeled J, L, I, and K, respectively, that were described previously [21].

Following single-point energy calculations, the four conformations were further geometry optimized (energy minimized) without constraining the dihedral angles. Calculations were undertaken for the four conformations in the gas phase and in several solvents that have dielectric constants spanning a wide range (2.2–78).

Calculated bond lengths, the C_{6}–C_{10}–C_{11} bond angle, and methoxy group dihedral angles for the various conformers in the gas phase and CCl_{4} are listed in Table 1. Similar trends in the listed data are calculated for the conformers in other solvents (data not shown). Data for in the binding site is also listed in Table 1.

The data presented in Table 1 demonstrates that all four conformers in solvent are within 0.45 kcal/mol in energy (*kT* at 298 K is ~0.59 kcal/mol), so all four conformers would be expected to be present to some degree in solvent at room temperature. The orientation of the methoxy groups of the four geometry-optimized conformers (in CCl_{4}) as well as the calculated dihedral angles are shown in the insets in Figure 2. The corresponding dihedral angles for the four conformers in different solvents are similar (data not shown).

The hydrocarbon chain (isoprene unit) attached at C_{6} makes a distinct kink at C_{10}. The C_{6}–C_{10}–C_{11} angle is close to 113° for all four conformers (Table 1). This angle is also similar to that found for occupying the binding site, respectively (Table 1).

##### 3.3. Calculated Vibrational Frequencies of

Figure 3(a) shows calculated IR spectra for the four conformers in CCl_{4}, in the 1530–1425 cm^{−1} region. This spectral region is chosen because it is the region where the main and modes of lie, and it is therefore the region generally focused upon in FTIR studies of in solution [8–10]. The spectra of the conformers in Figure 3(a) have been scaled by the appropriate Boltzmann factors, which were calculated based on the relative energies of the four conformations (Table 1). A composite spectrum which is the sum of the four Boltzmann weighted spectra is also shown in Figure 3(a). The corresponding calculated composite spectra for in various solvents are presented in Figures 3(b) and 3(c).

In the composite spectra an intense band is observed at 1500–1478 cm^{−1}, depending on the solvent. The frequency of this absorption band decreases, and the intensity increases, as the dielectric constant of the solvent increases. The frequency changes as a function of dielectric constant are outlined in the inset in Figure 3(a), which demonstrates that the band frequency is strongly solvent dependant only for solvents with dielectric constant ranging from ~1 to 20. Similar results have been found for PCM calculations of small neutral ketones in nonprotic solvents [22].

The calculated composite spectrum for in CCl_{4} (Figure 3(a)) displays an intense band at 1493 cm^{−1}. Lower intensity peaks are observed at 1483 and 1450 cm^{−1}. Table 2(b) lists the frequencies, IR intensities, Raman activities, and potential energy distributions for the normal modes that contribute to the bands in the spectra of in CCl_{4}. Similar results are calculated for conformers A, B, and E (data not shown), as would be expected given the similarity in the spectra of the conformers in Figure 3(a). For comparison, Table 2(a) also lists data for in the gas phase.

For in CCl_{4}, the band at 1493 cm^{−1} (Figure 3(a)) is due to two intense normal modes at ~1491 and ~1495 cm^{−1}. The 1491 cm^{−1} normal mode is due predominantly to stretching [R9(56%)] while the 1491 cm^{−1} normal mode is due predominantly to stretching [R3(46%)]. For all four conformers, the and groups vibrate separately at a similar frequency with similar intensity. This is also observed for in other solvents (not shown). This behavior is different from that found in calculations for neutral UQ_{1}, however, where most of the intensity is in only one of the C=O modes [21].

For in the gas phase the most intense band is calculated at 1500 cm^{−1}. In gas phase calculations, however, this band is due to the out-of-phase vibration of both groups [R3(29%)–R9(27%)] (Table 2(a)). In gas phase calculations, the in-phase vibration of both groups is found at 1495 cm^{−1}, and it is approximately a factor of seven lower in intensity than the out-of-phase vibration (Table 2(a)). In gas phase calculations the in-phase vibration is very strongly Raman active while the out-of-phase vibration is not. In contrast, for calculations in CCl_{4}, both the and vibrations are strongly Raman active.

In all spectra in Figure 3 a weak band is found at 1522 cm^{−1}. This band is due to an out-of-phase vibration of the groups of the quinone ring (R4–R10). Given the antisymmetric nature of the vibration it is very weakly Raman active. The in-phase vibration of the groups of the quinone ring occurs at 1607 cm^{−1} and is IR silent but very strongly Raman active (Table 2).

A relatively intense band is found at 1450–1456 cm^{−1} in all of the spectra in Figure 3. This band is due predominantly to CH bending vibrations of both methoxy methyl groups (*δ*C8 and *δ*C9) (Table 2(b)). Given that the mode is due to CH bending vibrations of the methoxy methyl groups it is not surprising that the precise frequency of this normal mode can vary by as much as 5 cm^{−1} among the four conformers (Figure 3(a)).

The calculated spectra of the four conformers are similar (Figure 3(a)). We also find that the spectra are very similar for isotope-labeled versions of the conformers (not shown). For this reason we will consider only the Boltzmann weighted composite spectra below. In addition, we will consider spectra for in CCl_{4}, noting that similar results and conclusions hold for in other solvents.

Figure 4 shows calculated IR (*left*) and Raman (*right*) spectra for unlabeled, ^{13}C, and ^{18}O isotope-labeled in the gas phase (a) and CCl_{4} (b). The normal modes (frequencies, intensities, Raman activities, and PEDs) that give rise to the bands in the spectra in Figure 4 are also listed in Table 2.

As discussed above, for unlabeled in CCl_{4} the 1493 cm^{−1} band (IR spectrum) is due to separate and vibrations. Upon ^{13}C labeling the 1493 cm^{−1} band appears to downshift 39 cm^{−1} to 1454 cm^{−1} (Figure 4(b)). Such a downshift is expected for a band that is due to groups. Table 2(b) indicates that the 1454 cm^{−1} band in the spectrum of ^{13}C labeled in CCl_{4} is due to a stretching vibration mixed with CH methyl bending vibrations (associated with both methoxy methyl groups). A very low-intensity normal mode at 1458 cm^{−1} also contributes to the 1454 cm^{−1} band in the IR spectrum of ^{13}C-labeled in CCl_{4}. This 1458 cm^{−1} mode is due to a stretching vibration mixed with CH methoxy methyl bending vibrations (Table 2(b)). So both groups give rise to intense normal modes for unlabeled in CCl_{4}. However, upon ^{13}C labeling, only one intense mode is found while the other is considerably weaker. Similar ^{13}C isotope-induced changes are found for calculations of in the gas phase (Table 2(a)).

The band at 1449 cm^{−1} in the IR spectrum for unlabeled in CCl_{4} is due predominantly to CH bending vibrations of both methoxy methyl groups. Upon ^{13}C labeling the 1449 cm^{−1} band downshifts from 18 cm^{−1} to 1431 cm^{−1} (Figure 4(b)). The 1431 cm^{−1} mode is due to the out-of-phase vibration of both groups [–R3(24%) + R9(11%)] coupled to a C_{3} methoxy methyl bending vibration [*δ*C8(32%)].

Other than the normal modes just discussed, stretching vibrations (R9 and R3) contribute to at least 6 other modes in the 1500–1400 cm^{−1} region for ^{13}C-labeled . Similar results are found for ^{13}C-labeled in the gas phase.

In the IR spectrum of unlabeled in CCl_{4} the weak band at 1523 cm^{−1} is due predominantly to an out-of-phase vibration (R4–R10). The in-phase vibration (R4 + R10) occurs at 1607 cm^{−1}, with negligible IR intensity but high Raman activity (Figure 4(b)). The in-phase mode downshifts 57 cm^{−1} to 1550 cm^{−1} upon ^{13}C labeling with little change in the mode composition (Table 2(b)).

In the unlabeled species the relatively pure modes are found near 1493 cm^{−1} (in CCl_{4}). Upon ^{18}O labeling these modes are expected to downshift from ~40 cm^{−1} to ~1453 cm^{−1}. In the unlabeled species the methoxy-methyl bending mode is found at 1450 cm^{−1}. So upon ^{18}O labeling the modes and methyl bending modes will be similar in frequency and are therefore expected to strongly mix. Table 2(b) shows that upon ^{18}O labeling the stretching vibrations do mix extensively with methyl bending vibrations, and that the modes are distributed amongst at least five different mixed modes. From the IR spectra in Figure 4(b), one could argue that the 1493 cm^{−1} band downshifts from 18 cm^{−1} to 1475 cm^{−1} upon ^{18}O labeling. Such an ^{18}O isotope-induced frequency shift for semiquinones is in line with experimental observations [8, 23] (see below).

#### 4. Discussion

The calculated changes in frequency, intensity, and mode composition upon isotope labeling of ubisemiquinone are considerably more complex than those found for the neutral species [21]. In spite of this, however, the calculated data allow a clear and detailed interpretation of bands in experimental Raman and FTIR spectra of ubisemiquinone. The calculated IR and Raman spectra presented in Figure 4(b) correspond very well to experimental spectra [8–10].

##### 4.1. Modeling Isotope-Induced Bandshifts Observed in Resonance Raman Spectra

Resonance Raman spectra of unlabeled and ^{13}C-labeled in the and binding sites in purple bacterial reaction centers have been obtained [10]. For comparison, resonance Raman spectra of unlabeled in solution were also obtained [10]. For both *in vivo* and *in vitro* cases an intense Raman band was observed near 1608 cm^{−1}, with weaker bands observed near 1523 and 1488 cm^{−1}. We note that the calculated Raman spectrum for unlabeled (Figure 4(b)) looks very similar to the experimental spectrum.

The ~1608 cm^{−1} band was assigned to a mode, weakly coupled to a mode, because it downshifted to 49–58 cm^{−1} upon ^{13}C labeling [10] (Table 3). The ~1488 cm^{−1} band was assigned to a mode because it downshifted to ~28 cm^{−1} upon ^{13}C labeling [10] (Table 3). The ~1523 cm^{−1} band apparently disappears upon ^{13}C labeling. Although not suggested in the original manuscript, it is possible that the 1523 cm^{−1} band (of in the binding site) downshifts from ~67 cm^{−1} to 1456 cm^{−1} upon ^{13}C labeling, and is masked by the band at 1456 cm^{−1} (Table 3).

In our calculations the and stretching vibrations couple to give in-phase and out-of-phase vibrations. For in CCl_{4} the in-phase vibration is at 1607 cm^{−1} (Figure 4(b)). This is 84 cm^{−1} higher in frequency than the out-of-phase vibration (at 1523 cm^{−1}). Unlike the out-of-phase vibration, the in-phase vibration is not coupled with methyl bending and carbonyl stretching modes (Table 2(b)). The in-phase vibration has negligible IR intensity but huge Raman activity. It is basically unaffected by ^{18}O labeling, but downshifts to 57 cm^{−1} upon ^{13}C labeling. The intensely Raman active band calculated at 1607 cm^{−1} for clearly corresponds to the band observed at ~1608 cm^{−1} experimentally [10].

The out-of-phase vibrational mode at 1523 cm^{−1} mixes with other modes upon ^{13}C labeling and is not easily identifiable. A strongly Raman active mode of ^{13}C-labeled is calculated at 1458 cm^{−1} (Table 2(b)). The out-of-phase vibration contributes 17% to the PED of this mode [R10(6%)–R4(11%)]. The in-phase-coupled vibration of both groups [R3(15%) + R9(7%)] also contributes to this mode.

Clearly, the calculated out-of-phase vibrational mode at 1523 cm^{−1} can be associated with the band observed at ~1521 cm^{−1} in resonance Raman spectra of in solution [10]. We suggest that the mode calculated at 1523 cm^{−1} forms part of a new mode that appears at 1458 cm^{−1} upon ^{13}C labeling (Table 3). In phase vibrations also contribute to the 1458 cm^{−1} mode. Our calculated data therefore provides an explanation as to why the ~1521 cm^{−1} resonance Raman band that is observed experimentally is not identified in spectra of ^{13}C-labeled [10]. Upon ^{13}C labeling the mode mixes with modes (and methyl bending modes) to become a new mode that is not distinctly identifiable as a mode.

Bands at 1486/1489 cm^{−1} in resonance Raman spectra of in the binding site downshift 30/27 cm^{−1} upon ^{13}C labeling of (Table 3), respectively. They were therefore associated with modes coupled to modes. Computationally, we find two modes at 1495 and 1491 cm^{−1}. Both modes are Raman active with the 1491 cm^{−1} mode displaying the greater activity (Table 2(b)). These modes give rise to the 1492 cm^{−1} band in the calculated Raman spectrum (Figure 3(b)), which appears to downshift 19/35 cm^{−1} to 1473/1457 cm^{−1} upon ^{13}C-labeling. The 1473/1457 cm^{−1} band in the calculated Raman spectrum for ^{13}C labeled is dominated by a mode at 1474/1458 cm^{−1}, respectively. The 1474 cm^{−1} mode and to a lesser degree the 1458 cm^{−1} mode are due predominantly to methyl CH bending vibrations of both methoxy groups coupled to a vibration. Notice that the coupling of the vibration is to methoxy methyl CH bending vibrations, *not * ring vibrations, as was originally proposed based on the experimental spectra.

##### 4.2. Modeling Isotope-Induced Bandshifts Observed in FTIR Spectra

Electrochemically generated FTIR difference spectra of UQ in various solvents have been obtained [8]. For in acetonitrile, THF, or dichloromethane an intense FTIR absorption band was observed at 1483–1488 cm^{−1}. The observation of predominantly a single intense band in experimental FTIR spectra of unlabeled and in solution is in line with our calculated IR spectra, which are dominated by an intense band at 1478–1493 cm^{−1} for in a variety of solvents (Figures 3(b) and 3(c)).

In experimental FTIR difference spectra for in dichloromethane, a band is observed at 1483 cm^{−1}, which downshifts to 41 cm^{−1} upon ^{13}C labeling (Table 3) [8]. From Figure 3(b) it can be seen that upon ^{13}C labeling the 1493 cm^{−1} band downshifts from 39 cm^{−1} to 1454 cm^{−1}. The calculated result therefore agrees very well with the experimental observation.

Experimentally, for in dichloromethane, it is also observed that the 1483 cm^{−1} band downshifts from 15 cm^{−1} to 1468 cm^{−1} upon ^{18}O labeling. From the calculated IR spectra in Figure 4(b), the most obvious suggestion is that the 1493 cm^{−1} band (of unlabeled ) downshifts from 18 cm^{−1} to 1475 cm^{−1} upon ^{18}O labeling. The calculated PEDs in Table 2(b) indicate a complicated situation: the 1495 and 1491 cm^{−1} modes of unlabeled are due to the groups (R3 and R9, resp.). Upon ^{18}O labeling modes appear at 1476 [R9(23%)] and 1465 cm^{−1} [R9(27%)]. Thus, the 1491 cm^{−1} mode in the unlabeled species appears to split and downshift to 15 and 26 cm^{−1} upon ^{18}O labeling (Table 3). The former is in excellent agreement with experiment [8]. Upon ^{18}O labeling modes also appear at 1481 [R3(14%)] and 1483 cm^{−1} [R3(7%)]. Thus, the 1495 cm^{−1} mode in the unlabeled species also appears to split and downshift to 14 and 12 cm^{−1} upon ^{18}O labeling (Table 3). Again, these conclusions are in good agreement with experiment [8]. It is the plethora of mixed modes that appear upon ^{18}O labeling that give rise to the broad band with a peak near 1475 cm^{−1} in the calculated spectrum (Figure 3(b)). Unfortunately FTIR spectra for ^{18}O-labeled have never been presented. Only the observed shifts upon labeling were presented.

From electrochemically generated FTIR difference spectra of ^{13}C-labeled in various solvents [8] a band was observed at 1412 cm^{−1}. It was suggested that this band was due to a vibration that was downshifted to 71 cm^{−1} from 1483 cm^{−1} in the unlabeled species. Neither the calculated data presented here nor the resonance Raman data presented previously support this hypothesis.

##### 4.3. Experimental and FTIR DS.

and FTIR DS have been obtained using PBRCs from *R. sphaeroides* [7, 24–27]. In FTIR DS three intense IR bands are observed near 1485, 1466, and 1449 cm^{−1} [24, 25]. On the basis of ^{18}O, ^{13}C, ^{13}C_{1}, and ^{13}C_{4} labeling the 1486/1466 cm^{−1} bands were assigned to / vibrations, respectively [24]. The modes were suggested to be considerably mixed. The origin of the 1449 cm^{−1} band was not considered.

Another group, which undertook identical labeling experiments [25], assigned the 1485 cm^{−1} band to a vibration, the 1466 cm^{−1} band to C_{4}=O vibration, and the 1449 cm^{−1} band to a vibration. All modes were suggested to be strongly mixed.

Resonance Raman spectra for in the binding site display a weak band at 1486 cm^{−1}, but no bands were apparent at 1466 and 1449 cm^{−1}. Of course it may simply be the case that the 1466 and 1449 cm^{−1} normal modes are Raman inactive.

Our calculated spectra for in solution poorly model observed FTIR bands of in the binding site. For in the binding site, the modes appear to be separated by 19 cm^{−1}. For calculations in solvent the two modes do appear to be distinct, although the separation of the modes is only 4 cm^{−1}. In gas phase calculations the two modes are coupled. The separation of modes of in the binding site is due to asymmetric interactions with the protein environment. Calculations of in solvent or in the gas phase cannot model these interactions. Calculations including effects of the protein environment are essential. Such calculations are underway in our lab.

In FTIR DS a single IR band is observed near 1479 cm^{−1}. It was suggested that this band was due to both modes of in the binding site [7, 26, 27]. It was also suggested that the 1479 cm^{−1} band downshifts 33/52 cm^{−1} upon ^{18}O/^{13}C labeling, respectively [7, 26, 27]. Such shifts are difficult to rationalize in view of the shifts calculated (15/37 cm^{−1}) and observed experimentally (15/27–41 cm^{−1}) for UQ in solution (Table 3). Additionally, there appears to be some inequivalence in the modes of UQ in the binding site when perturbed specifically at the C_{1} or C_{4} position [7, 26, 27]. It was suggested that this inequivalence is a result of specific protein interactions [7, 26, 27]. Again, calculations including effects of the protein environment appear to be necessary (essential) in order to accurately simulate the vibrational spectra of UQ in the binding site.

##### 4.4. Previous Calculations of Ubisemiquinones

DFT-based vibrational frequency calculations (using the BP86 functional) have been undertaken for 2,3-dimethoxy-1,4-benzoquinone and 2,3-dimethoxy-5,6-dimethyl-1,4-benzoquinone in the gas phase [12]. Comparison of calculated data for the two models showed that substituents at C_{5} and C_{6} are required in order to better model the properties of ubiquinones and ubisemiquinones. In the above study isotope shifts were calculated. However how the and modes couple with each other and with CH methoxy methyl bending vibrations was not considered. As we have shown above, the extent of mode mixing can be considerably altered upon labeling, making it difficult to identify how the different bands shift upon labeling. As we show here, the detailed PEDs are a crucial tool in the analysis of how calculated bands shift upon isotope labeling.

One problem with previous DFT calculations (in the gas phase) [12] is that for 2,3-dimethoxy-1,4-benzoquinone the modes were found at a higher frequency than the modes. For 2,3-dimethoxy-5,6-dimethyl-1,4-benzoquinone (in the gas phase) the modes were found at slightly lower frequency than the modes (3-4 cm^{−1}). However, from Raman experiments the out-of-phase mode is found to be ~32 cm^{−1} higher in frequency than the mode [10] (Table 3).

Furthermore, the antisymmetrically coupled mode (for 2,3-dimethoxy-5,6-dimethyl-1,4-benzoquinone) was calculated to be more than a factor of 26 times more intense than the mode [12]. This calculated result is not in line with experimental IR spectra [8].

Clearly, previous DFT calculations [12] poorly model the experimental Raman and IR spectra. In contrast, in our calculations for in CCl_{4}, the out-of-phase modes are 28–32 cm^{−1} higher in frequency than either of the mode (Table 2(b)). In gas phase calculations the out-of-phase mode is still 24 cm^{−1} higher in frequency than the antisymmetrically coupled mode (Table 2(a)). In addition, in gas phase calculations and in solvent, the intensity of antisymmetrically coupled mode is ~7.5 times more intense than the out-of-phase mode. These results are in excellent agreement with experimental IR and Raman spectra. The limitations in previous calculations are most likely related to the choice of functional and basis set, and the inadequacy of a UQ structural model that lacks an isoprene unit.

#### 5. Conclusions

We calculate that four conformers are likely present in solution at room temperature. Calculated IR spectra for all four conformers are similar. Calculated IR spectra of unlabeled and isotope-labeled in the gas phase and in solution show a similar band pattern, although in some cases there are differences in the composition of the modes that contribute to the bands in the spectra.

Calculations show that upon isotope labeling the out-of-phase ring modes and modes of strongly couple with methyl C–H bending vibrations of the methoxy groups. This leads to complicated splitting of modes and unusual downshifts upon isotope labeling. Nonetheless by consideration of PEDs of the calculated normal modes, sense can be made of the isotope-induced shifts and intensity changes, and it is shown that the calculated data provide a rational and detailed interpretation of experimentally observed isotope-induced band shifts in experimental FTIR and Raman spectra of in solution.

#### Abbreviations

DFT: | Density functional theory |

DS: | Difference spectra |

FTIR: | Fourier transform infrared |

IR: | Infrared |

IEF: | Integral equation formalism |

PBRCs: | Purple bacterial reaction centers |

PCM: | Polarizable continuum model |

PED: | Potential energy distribution |

UQ: | Ubiquinone |

: | Ubisemiquinone. |

#### Acknowledgments

H. P. Lamichhane acknowledge support from a fellowship from the Molecular Basis of Disease Program at Georgia State University. G. Hastings acknowledges the support from Qatar National Research Fund.

#### References

- B. Trumpower,
*Function of Quinones in Energy Conserving Systems*, Academic Press, 1982. - B. Ke, “The bacterial photosynthetic reaction center: chemical composition and crystal structure,” in
*Photosynthesis: Photobiochemistry and Photobiophysics*, pp. 47–62, Kluwer Academic Publishers, Dordrecht, The Netherlands, 2001. View at Google Scholar - B. Ke, “The, “Stable” primary electron acceptor (${\text{Q}}_{\text{A}}$) of photosynthetic bacteria,” in
*Photosynthesis: Photobiochemistry and Photobiophysics*, pp. 101–110, Kluwer Academic Publishers, Dordrecht, The Netherlands, 2001. View at Google Scholar - B. Ke, “The secondary electron acceptor (${\text{Q}}_{\text{B}}$) of photosynthetic bacteria,” in
*Photobiochemistry and Photobiophysics*, pp. 111–128, Kluwer Academic Publishers, Dordrecht, The Netherlands, 2001. View at Google Scholar - N. Srinivasan and J. H. Golbeck, “Protein-cofactor interactions in bioenergetic complexes: the role of the A1A and A1B phylloquinones in Photosystem I,”
*Biochimica et Biophysica Acta*, vol. 1787, no. 9, pp. 1057–1088, 2009. View at Publisher · View at Google Scholar · View at Scopus - C. A. Wraight and M. R. Gunner, “The acceptor quinones of purple photosynthetic Bacteria-structure and spectroscopy,” in
*The Purple Photosynthetic Bacteria*, C. N. Hunter, F. Daldal, M. C. Thurnauer, and J. T. Beatty, Eds., pp. 379–405, Springer, 2009. View at Google Scholar - J. Breton and E. Nabedryk, “Protein-quinone interactions in the bacterial photosynthetic reaction center: light-induced FTIR difference spectroscopy of the quinone vibrations,”
*Biochimica et Biophysica Acta*, vol. 1275, no. 1-2, pp. 84–90, 1996. View at Publisher · View at Google Scholar · View at Scopus - M. Bauscher and W. Mäntele, “Electrochemical and infrared-spectroscopic characterization of redox reactions of p-quinones,”
*Journal of Physical Chemistry*, vol. 96, no. 26, pp. 11101–11108, 1992. View at Google Scholar · View at Scopus - M. Bauscher, E. Nabedryk, K. Bagley, J. Breton, and W. Mantele, “Investigation of models for photosynthetic electron acceptors. Infrared spectroelectrochemistry of ubiquinone and its anions,”
*FEBS Letters*, vol. 261, no. 1, pp. 191–195, 1990. View at Publisher · View at Google Scholar · View at Scopus - X. Zhao, T. Ogura, M. Okamura, and T. Kitagawa, “Observation of the resonance raman spectra of the semiquinones QA
*·*- and QB*·*- in photosynthetic reaction centers from*Rhodobacter sphaeroides*R26,”*Journal of the American Chemical Society*, vol. 119, pp. 5263–5264, 1997. View at Google Scholar - G. Balakrishnan, P. Mohandas, and S. Umapathy, “Ab initio studies on structure and vibrational spectra of ubiquinone and its radical anion,”
*Spectrochimica Acta A*, vol. 53, no. 10, pp. 1553–1561, 1997. View at Google Scholar · View at Scopus - M. Nonella, “A density functional investigation of model molecules for ubisemiquinone radical anions,”
*Journal of Physical Chemistry B*, vol. 102, no. 21, pp. 4217–4225, 1998. View at Google Scholar · View at Scopus - M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., 2004. View at Publisher · View at Google Scholar
- K. M. Bandaranayake, V. Sivakumar, R. Wang, and G. Hastings, “Modeling the A1 binding site in photosystem. I. Density functional theory for the calculation of “anion—neutral“ FTIR difference spectra of phylloquinone,”
*Vibrational Spectroscopy*, vol. 42, no. 1, pp. 78–87, 2006. View at Publisher · View at Google Scholar · View at Scopus - E. Cancès, C. Le Bris, B. Mennucci, and J. Tomasi, “Integral equation methods for molecular scale calculations in the liquid phase,”
*Mathematical Models and Methods in Applied Sciences*, vol. 9, no. 1, pp. 35–44, 1999. View at Google Scholar · View at Scopus - E. Cancès, B. Mennucci, and J. Tomasi, “A new integral equation formalism for the polarizable continuum model: theoretical background and applications to Isotropic and anisotropic dielectrics,”
*Journal of Chemical Physics*, vol. 107, no. 8, pp. 3032–3041, 1997. View at Google Scholar · View at Scopus - J. Tomasi, B. Mennucci, and E. Cancès, “The IEF version of the PCM solvation method: an overview of a new method addressed to study molecular solutes at the QM ab initio level,”
*Journal of Molecular Structure*, vol. 464, no. 1–3, pp. 211–226, 1999. View at Publisher · View at Google Scholar · View at Scopus - J. Tomasi, R. Cammi, B. Mennucci, C. Cappelli, and S. Corni, “Molecular properties in solution described with a continuum solvation model,”
*Physical Chemistry Chemical Physics*, vol. 4, no. 23, pp. 5697–5712, 2002. View at Publisher · View at Google Scholar · View at Scopus - J. Tomasi, B. Mennucci, and R. Cammi, “Quantum mechanical continuum solvation models,”
*Chemical Reviews*, vol. 105, no. 8, pp. 2999–3093, 2005. View at Publisher · View at Google Scholar · View at Scopus - J. M. L. Martin and C. Van Alsenoy,
*GAR2PED*, University of Antwerp, 1995. - H. Lamichhane, R. Wang, and G. Hastings, “Comparison of calculated and experimental FTIR spectra of specifically labeled ubiquinones,”
*Vibrational Spectroscopy*, vol. 55, no. 2, pp. 279–286, 2011. View at Publisher · View at Google Scholar · View at Scopus - C. Cappelli, C. O. Silva, and J. Tomasi, “Solvent effects on vibrational modes: ab-initio calculations, scaling and solvent functions with applications to the carbonyl stretch of dialkyl ketones,”
*Journal of Molecular Structure*, vol. 544, pp. 191–203, 2001. View at Publisher · View at Google Scholar · View at Scopus - G. Hastings, K. M. P. Bandaranayake, and E. Carrion, “Time-resolved FTlR difference spectroscopy in combination with specific isotope labeling for the study of A1, the secondary electron acceptor in photosystem 1,”
*Biophysical Journal*, vol. 94, no. 11, pp. 4383–4392, 2008. View at Publisher · View at Google Scholar · View at Scopus - J. Breton, J. R. Burie, C. Berthomieu, G. Berger, and E. Nabedry, “The binding sites of quinones in photosynthetic bacterial reaction centers investigated by light-induced FTIR difference spectroscopy: assignment of the ${\text{Q}}_{\text{A}}$ vibrations in
*Rhodobacter sphaeroides*using^{18}O- Or^{13}C-labeled ubiquinone and vitamin K_{1},”*Biochemistry*, vol. 33, no. 16, pp. 4953–4965, 1994. View at Google Scholar · View at Scopus - R. Brudler, H. J. M. De Groot, W. B. S. Van Liemt et al., “Asymmetric binding of the 1- and 4- $\text{C}=\text{O}$ groups of ${\text{Q}}_{\text{A}}$ in
*Rhodobacter sphaeroides*R26 reaction centres monitored by Fourier transform infra-red spectroscopy using site-specific isotopically labelled ubiquinone-10,”*EMBO Journal*, vol. 13, no. 23, pp. 5523–5530, 1994. View at Google Scholar · View at Scopus - J. Breton, C. Boullais, G. Berger, C. Mioskowski, and E. Nabedryk, “Binding sites of quinones in photosynthetic bacterial reaction centers investigated by light-induced FTIR difference spectroscopy: symmetry of the carbonyl interactions and close equivalence of the QB vibrations in
*Rhodobacter sphaeroides*and Rhodopseudomonas viridis probed by isotope labeling,”*Biochemistry*, vol. 34, no. 36, pp. 11606–11616, 1995. View at Google Scholar · View at Scopus - R. Brudler, H. J. M. De Groot, W. B. S. Van Liemt et al., “FTIR spectroscopy shows weak symmetric hydrogen bonding of the ${\text{Q}}_{\text{B}}$ carbonyl groups in
*Rhodobacter sphaeroides*R26 reaction centres,”*FEBS Letters*, vol. 370, no. 1-2, pp. 88–92, 1995. View at Publisher · View at Google Scholar · View at Scopus