Origin of the 701-nm fluorescence emission of the Lhca2 subunit of higher plant photosystem I

Photosystem I of higher plants is characterized by red-shifted spectral forms deriving from chlorophyll chromophores. Each of the four Lhca1 to -4 subunits exhibits a specific fluorescence emission spectrum, peaking at 688, 701, 725, and 733 nm, respectively. Recent analysis revealed the role of chlorophyll-chlorophyll interactions of the red forms in Lhca3 and Lhca4, whereas the basis for the fluorescence emission at 701 nm in Lhca2 is not yet clear. We report a detailed characterization of the Lhca2 subunit using molecular biology, biochemistry, and spectroscopy and show that the 701-nm emission form originates from a broad absorption band at 690 nm. Spectroscopy on recombinant mutant proteins assesses that this band represents the low energy form of an excitonic interaction involving two chlorophyll a molecules bound to sites A5 and B5, the same protein domains previously identified for Lhca3 and Lhca4. The resulting emission is, however, substantially shifted to higher energies. These results are discussed on the basis of the structural information that recently became available from x-ray crystallography (Ben Shem, A., Frolow, F., and Nelson, N. (2003) Nature 426, 630-635). We suggest that, within the Lhca subfamily, spectroscopic properties of chromophores are modulated by the strength of the excitonic coupling between the chromophores A5 and B5, thus yielding fluorescence emission spanning a large wavelength interval. It is concluded that the interchromophore distance rather than the transition energy of the individual chromophores or the orientation of transition vectors represents the critical factor in determining the excitonic coupling in Lhca pigment-protein complexes.     

Light harvesting in photosystem I: modeling based on the 2.5-A structure of photosystem I from Synechococcus elongatus

The structure of photosystem I from the thermophilic cyanobacterium Synechococcus elongatus has been recently resolved by x-ray crystallography to 2.5-A resolution. Besides the reaction center, photosystem I consists also of a core antenna containing 90 chlorophyll and 22 carotenoid molecules. It is their function to harvest solar energy and to transfer this energy to the reaction center (RC) where the excitation energy is converted into a charge separated state. Methods of steady-state optical spectroscopy such as absorption, linear, and circular dichroism have been applied to obtain information on the spectral properties of the complex, whereas transient absorption and fluorescence studies reported in the literature provide information on the dynamics of the excitation energy transfer. On the basis of the structure, the spectral properties and the energy transfer kinetics are simultaneously modeled by application of excitonic coupling theory to reveal relationships between structure and function. A spectral assignment of the 96 chlorophylls is suggested that allows us to reproduce both optical spectra and transfer and emission spectra and lifetimes of the photosystem I complex from S. elongatus. The model calculation allowed to study the influence of the following parameters on the excited state dynamics: the orientation factor, the heterogeneous site energies, the modifications arising from excitonic coupling (redistribution of oscillator strength, energetic splitting, reorientation of transition dipoles), and presence or absence of the linker cluster chlorophylls between antenna and reaction center. For the F?rster radius and the intrinsic primary charge separation rate, the following values have been obtained: R(0) = 7.8 nm and k(CS) = 0.9 ps(-1). Variations of these parameters indicate that the excited state dynamics is neither pure trap limited, nor pure transfer (to-the-trap) limited but seems to be rather balanced.     

How proteins trigger excitation energy transfer in the FMO complex of green sulfur bacteria

A simple electrostatic method for the calculation of optical transition energies of pigments in protein environments is presented and applied to the Fenna-Matthews-Olson (FMO) complex of Prosthecochloris aestuarii and Chlorobium tepidum. The method, for the first time, allows us to reach agreement between experimental optical spectra and calculations based on transition energies of pigments that are calculated in large part independently, rather than fitted to the spectra. In this way it becomes possible to understand the molecular mechanism allowing the protein to trigger excitation energy transfer reactions. The relative shift in excitation energies of the seven bacteriochlorophyll-a pigments of the FMO complex of P. aestuarii and C. tepidum are obtained from calculations of electrochromic shifts due to charged amino acids, assuming a standard protonation pattern of the protein, and by taking into account the three different ligand types of the pigments. The calculations provide an explanation of some of the earlier results for the transition energies obtained from fits of optical spectra. In addition, those earlier fits are verified here by using a more advanced theory of optical spectra, a genetic algorithm, and excitonic couplings obtained from electrostatic calculations that take into account the influence of the dielectric protein environment. The two independent calculations of site energies strongly favor one of the two possible orientations of the FMO trimer relative to the photosynthetic membrane, which were identified by electron microscopic studies and linear dichroism experiments. Efficient transfer of excitation energy to the reaction center requires bacteriochlorophylls 3 and 4 to be the linker pigments. The temporal and spatial transfer of excitation energy through the FMO complex is calculated to proceed along two branches, with transfer times that differ by an order of magnitude.     

Exciton interaction in molecular beacons: a sensitive sensor for short range modifications of the nucleic acid structure

Molecular beacons are hairpin-shaped, single-stranded oligonucleotides constituting sensitive fluorescent DNA probes widely used to report the presence of specific nucleic acids. In its closed form the stem of the hairpin holds the fluorophore covalently attached to one end, close to the quencher, which is covalently attached to the other end. Here we report that in the closed form the fluorophore and the quencher form a ground state intramolecular heterodimer whose spectral properties can be described by exciton theory. Formation of the heterodimers was found to be poorly sensitive to the stem sequence, the respective positions of the dyes and the nature of the nucleic acid (DNA or RNA). The heterodimer allows strong coupling between the transition dipoles of the two chromophores, leading to dramatic changes in the absorption spectrum that are not compatible with a Förster-type fluorescence resonance energy transfer (FRET) mechanism. The excitonic heterodimer and its associated absorption spectrum are extremely sensitive to the orientation of and distance between the dyes. Accordingly, the application of molecular beacons can be extended to monitoring short range modifications of the stem structure. Moreover, the excitonic interaction was also found to operate for doubly end-labeled duplexes.     

VenusA206 Dimers Behave Coherently at Room Temperature

Fluorescent proteins (FPs) have revolutionized cell biology by allowing genetic tagging of specific proteins inside living cells. In conjunction with Förster’s resonance energy transfer (FRET) measurements, FP-tagged proteins can be used to study protein-protein interactions and estimate distances between tagged proteins. FRET is mediated by weak Coulombic dipole-dipole coupling of donor and acceptor fluorophores that behave independently, with energy hopping discretely and incoherently between fluorophores. Stronger dipole-dipole coupling can mediate excitonic coupling in which excitation energy is distributed near instantaneously between coherently interacting excited states that behave as a single quantum entity. The interpretation of FP energy transfer measurements to estimate separation often assumes that donors and acceptors are very weakly coupled and therefore use a FRET mechanism. This assumption is considered reasonable as close fluorophore proximity, typically associated with strong excitonic coupling, is limited by the FP β-barrel structure. Furthermore, physiological temperatures promote rapid vibrational dephasing associated with a rapid decoherence of fluorophore-excited states. Recently, FP dephasing times that are 50 times slower than traditional organic fluorophores have been measured, raising the possibility that evolution has shaped FPs to allow stronger than expected coupling under physiological conditions. In this study, we test if excitonic coupling between FPs is possible at physiological temperatures. FRET and excitonic coupling can be distinguished by monitoring spectral changes associated with fluorophore dimerization. The weak coupling mediating FRET should not cause a change in fluorophore absorption, whereas strong excitonic coupling causes Davydov splitting. Circular dichroism spectroscopy revealed Davydov splitting when the yellow FP Venus_A206 dimerizes, and a novel approach combining photon antibunching and fluorescence correlation spectroscopy was used to confirm that the two fluorophores in a Venus_A206 homodimer behave as a single-photon emitter. We conclude that excitonic coupling between Venus_A206 fluorophores is possible at physiological temperatures.     

Significance of the Excitonic Intensity Borrowing in the J-/H-aggregates of Bacteriochlorophylls/Chlorophylls

A quantitative analysis of the excitonic intensity borrowing for the J-/H-aggregates of the bacteriochlorophylls/chlorophylls (BChls/Chls) in specific, and of porphyrins in general, is presented. The analysis is based on the argument that the mixing between the two energetically well-separated bands, such as the Q and B bands of BChls/Chls, should be considered important if the aggregated system possesses an excitonic superstate. A remarkably simple explanation of the significance of the excitonic intensity borrowing is given: superhyperchromism is manifested by the mediation of interband coupling between the superstates in␣the two well-separated bands of such aggregates. A comprehensive discussion on the significance of superhyperchromism and on its size-dependence is provided in connection with its effects on the absorption spectra of the BChl/Chl J- and H-aggregates.     

Visualization of excitonic structure in the Fenna-Matthews-Olson photosynthetic complex by polarization-dependent two-dimensional electronic spectroscopy

Photosynthetic light-harvesting proceeds by the collection and highly efficient transfer of energy through a network of pigment-protein complexes. Interchromophore electronic couplings and interactions between pigments and the surrounding protein determine energy levels of excitonic states, and dictate the mechanism of energy flow. The excitonic structure (orientation of excitonic transition dipoles) of pigment-protein complexes is generally deduced indirectly from x-ray crystallography, in combination with predictions of transition energies and couplings in the chromophore site basis. We demonstrate that coarse-grained, excitonic, structural information in the form of projection angles between transition dipole moments can be obtained from the polarization-dependent, two-dimensional electronic spectroscopy of an isotropic sample, particularly when the nonrephasing or free polarization decay signal, rather than the photon echo signal, is considered. This method provides an experimental link between atomic and electronic structure, and accesses dynamical information with femtosecond time resolution. In an investigation of the Fenna-Matthews-Olson complex from green sulfur bacteria, the energy transfer connecting two particular exciton states in the protein was isolated as the primary contributor to a crosspeak in the nonrephasing two-dimensional spectrum at 400 femtoseconds under a specific sequence of polarized excitation pulses. The results suggest the possibility of designing experiments using combinations of tailored polarization sequences to separate and monitor individual relaxation pathways.     

Dihedral angles of tripeptides in solution directly determined by polarized Raman and FTIR spectroscopy

The amide I mode of the peptide linkage is highly delocalized in peptides and protein segments due to through-bond and through-space vibrationally coupling between adjacent peptide groups. J. Phys. Chem. B. 104:11316-11320) used coherent femtosecond infrared (IR) spectroscopy to determine the excitonic coupling energy and the orientational angle between the transition dipole moments of the interacting amide I modes of cationic tri-alanine in D(2)O. Recently, the same parameters were determined for all protonation states of tri-alanine by analyzing the amide I bands in the respective IR and isotropic Raman spectra (. J. Am. Chem. Soc. 119:1720-1726.). In both studies, the dihedral angles phi and psi were then obtained by utilizing the orientational dependence of the coupling energy obtained from ab initio calculations on tri-glycine in vacuo (. J. Raman Spectrosc. 29:81-86) to obtain an extended 3(1) helix-like structure for the tripeptide. In the present paper, a novel algorithm for the analysis of excitonic coupling between amide I modes is presented, which is based on the approach by Schweitzer-Stenner et al. but avoids the problematic use of results from ab initio calculations. Instead, the dihedral angles are directly determined from infrared and visible polarized Raman spectra. First, the interaction energy and the corresponding degree of wave-function mixing were obtained from the amide I profile in the isotropic Raman spectrum. Second, the depolarization ratios and the amide I profiles in the anisotropic Raman and IR-absorption spectra were used to determine the orientational angle between the peptide planes and the transition dipole moments, respectively. Finally, these two geometric parameters were utilized to determine the dihedral angles phi and psi between the interacting peptide groups. Stable extended conformations with dihedral angles in the beta-sheet region were obtained for all protonation states of tri-alanine, namely phi(+) = -126 degrees, psi(+) = 178 degrees; phi(+/-) = -110 degrees, psi(+/-) = 155 degrees; and phi(-) = -127 degrees, psi(-) = 165 degrees for the cationic, zwitterionic, and anionic state, respectively. These values reflect an extended beta-helix structure. Tri-glycine was found to be much more heterogeneous in that different extended conformers coexist in the cationic and zwitterionic state, which yield a noncoincidence between isotropic and anisotropic Raman scattering. Our study introduces vibrational spectroscopy as a suitable tool for the structure analysis of peptides in solution and tripeptides as suitable model systems for investigating the role of local interactions in determining the propensity of peptide segments for distinct secondary structure motifs.     

Excitation energy transfer in the LHC-II trimer: a model based on the new 2.72 Å structure

Energy transfer of the light harvesting complex LHC-II trimer, extracted from spinach, was studied in the Q(y) region at room temperature by femtosecond transient absorption spectroscopy. Configuration interaction exciton method [Linnanto et al. (1999) J Phys Chem B 103: 8739-8750] and 2.72 A structural information reported by Liu et al. was used to calculate spectroscopic properties and excitation energy transfer rates of the complex. Site energies of the pigments and coupling constants of pigment pairs in close contact were calculated by using a quantum chemical configuration interaction method. Gaussian random variation of the diagonal and off-diagonal exciton matrix elements was used to account for inhomogeneous broadening. Rate calculations included only the excitonic states initially excited and probed in the experiments. A kinetic model was used to simulate time and wavelength dependent absorption changes after excitation on the blue side of the Q(y) transition and compared to experimentally recorded rates. Analysis of excitonic wavefunctions allowed identification of pigments initially excited and probed into later. It was shown that excitation of the blue side of the Q(y) band of a single LHC-II complex results in energy transfer from chlorophyll b's of the lumenal side to chlorophyll a's located primarly on one of the monomers of the stromal side.     

The enhancement of electrostriction caused by lowering the solvent dielectric constant leads to the decrease of activation energy in trypsin catalysis.

The contribution of electrostriction of the solvent to the stabilization of the negatively charged tetrahedral transition state of a trypsin-catalyzed reaction was probed by means of kinetic studies involving high-pressure and solvent dielectric constant. A good correlation was observed between the increased catalytic efficiency of trypsin and the decreased solvent dielectric constant. When the dielectric constant of the solvents was lowered by 4.68 units, the loss of activation energy and that of free energy of activation were 2.26 kJ/mol and 3.09 kJ/mol, respectively. The activation volume for k(cat) decreased significantly as the dielectric constant of the solvent decreased, indicating that the degree of electrostriction of the solvent around the charged tetrahedral transition state has been enhanced. These observations demonstrate that the increase in the catalytic efficiency of the trypsin reaction with decreasing dielectric constant resulted from the stabilization of electrostatic energy for the formation of an oxyanion hole, and this stabilization was caused by the increase of electrostricted water around the charged tetrahedral transition state. Therefore, we conclude that control of the solvent dielectric constant can stabilize the tetrahedral transition state, and this lowers the activation energy.     

Significance of molecular crowding in grana membranes of higher plants for light harvesting by photosystem II

Significance of molecular crowding in grana thylakoids of higher plants on photosystem II function was studied by 'titrating' the naturally high protein density by fusing unilamellar liposomes of the native lipid mixture with isolated grana membranes (BBY). The incorporation of lipids was monitored by equilibrium density gradient centrifugation and two-dimensional thin layer chromatography. The excitonic coupling between light-harvesting (LHC) II and photosystem (PS) II was analysed by chlorophyll a fluorescence spectroscopy. The fluorescence parameters Fv/Fm and Fo clearly depend on the protein density indicating the importance of molecular crowding for establishing an efficient excitonic protein network. In addition the strong dependency of Fo on the protein density reveals weak interactions between LHCII complexes which could be important for dynamic adjustment of the photosynthetic apparatus in higher plants.     

Fundamental considerations in the effect of molecular weight on the glass transition of the gelatin/cosolute system

Four molecular fractions of gelatin produced by alkaline hydrolysis of collagen were investigated in the presence of cosolute to record the mechanical properties of the glass transition in high-solid preparations. Dynamic oscillatory and stress relaxation moduli in shear were recorded from 40°C to temperatures as low as -60°C. The small-deformation behavior of these linear polymers was separated by the method of reduced variables into a basic function of time alone and a basic function of temperature alone. The former allowed the reduction of isothermal runs into a master curve covering 17 orders of magnitude in the time domain. The latter follows the passage from the rubbery plateau through the glass transition region to the glassy state seen in the variation of shift factor, a(T) , as a function of temperature. The mechanical glass transition temperature (T(g) ) is pinpointed at the operational threshold of the free volume theory and the predictions of the reaction rate theory. Additional insights into molecular dynamics are obtained via the coupling model of cooperativity, which introduces the concept of coupling constant or interaction strength of local segmental motions that govern structural relaxation at the vicinity of T(g) . The molecular weight of the four gelatin fractions appears to have a profound effect on the transition temperature or coupling constant of vitrified matrices, as does the protein chemistry in relation to that of amorphous synthetic polymers or gelling polysaccharides.     

The Effect of Inserted Gold Nanosphere on the Second Harmonic Generation (SHG) Enhancement Factor of Three-Layered Dielectric-Gold Nanoshell

The tunable second harmonic generation (SHG) enhancement factor of gold-dielectric-gold three-layered nanoshells has been theoretically studied using the theory of quasi-static electrodynamics and plasmon hybridization. Because of the local surface plasmon resonance (LSPR)-induced local field effect, the SHG response corresponding to both fundamental frequency and second harmonic has been greatly enhanced. By changing the geometry parameters and local dielectric environment of the three-layered nanostructure, the intensity and shift of the SHG factor peaks could be fine tuned. As the radius of the inner gold sphere is increased, both the fundamental and the second harmonic SHG peaks from the anti-symmetric coupling between the outer bonding shell plasmon and the inner sphere plasmon decrease, whereas the SHG peaks from the symmetric coupling between the outer shell and the inner sphere get intense. These radius-dependent intensity changes of the SHG peaks also depend on the dielectric constant of the separate layer and outer surrounding. Thus, the number of SHG peak could be tuned from two to four. Furthermore, the wavelength gaps between the SHG peaks corresponding to anti-symmetric and symmetric coupling could be greatly reduced by increasing the thickness of the outer gold shell. Therefore, the nonmonotonous intensity change could be observed because of the switching of the SHG peaks. The corresponding physical origin has been illuminated by analyzing the plasmon hybridization and the polarization fields in the nanostructure.

     

Thermodynamics of a diffusional protein folding reaction

The folding reactions of several proteins are well described as diffusional barrier crossing processes, which suggests that they should be analyzed by Kramers' rate theory rather than by transition state theory. For the cold shock protein Bc-Csp from Bacillus caldolyticus, we measured stability and folding kinetics, as well as solvent viscosity as a function of temperature and denaturant concentration. Our analysis indicates that diffusional folding reactions can be treated by transition state theory, provided that the temperature and denaturant dependence of the solvent viscosity is properly accounted for, either at the level of the measured rate constants or of the calculated activation parameters. After viscosity correction the activation barriers for folding become less enthalpic and more entropic. The transition from an enthalpic to an entropic folding barrier with increasing temperature is, however, apparent in the data before and after this correction. It is a consequence of the negative activation heat capacity of refolding, which is independent of solvent viscosity. Bc-Csp and its mesophilic homolog Bs-CspB from Bacillus subtilis differ strongly in stability but show identical enthalpic and entropic barriers to refolding. The increased stability of Bc-Csp originates from additional enthalpic interactions that are established after passage through the activated state. As a consequence, the activation enthalpy of unfolding is increased relative to Bs-CspB.     

Kinetics of the lactate dehydrogenase reaction in high-viscosity media

The effect of the medium viscosity on kinetics parameters of lactate dehydrogenase reaction was studied. The viscosity increase results in a sharp decline in the catalytic rate for both the pyruvate reduction and lactate oxidation reactions. It is shown that the catalytic step and its associated conformational motions is the only step which is considerably retarded when the viscosity increases. The reaction is not sensitive to changes in the dielectric properties of the medium. An inverse power function observed between the rate constant and viscosity cannot be explained by the theory of absolute reaction rates. However, it can easily be interpreted on the basis of the Kramers theory dealing with the transition over the activation barrier as a diffusional motion in the field of random forces. The influence of the medium's viscosity on the kinetic parameters indicates the existence of strong coupling between the dynamics of the solvent and the conformational motions of the protein molecule, which are correlated with the catalytic step.     

Orientation of the chromophores in the reaction center of Rhodopseudomonas viridis. Comparison of low-temperature linear dichroism spectra with a model derived from X-ray crystallography

Linear dichroism (LD) and absorption (A) spectra of reaction centers from Rhodopseudomonas viridis included in the native chromatophores or reconstituted in planar aggregates have been recorded at 10 K. The samples were oriented in squeezed polyacrylamide gels and the primary donor P was in the reduced or (chemically) oxidized state. The LD spectra of reaction centers in these two states are in favor of a dimeric model of P in which excitonic coupling between the two non-parallel Q_Y transitions leads to a main transition at 990 nm (parallel to the membrane plane) and another one of smaller oscillator strength at 850 nm (tilted at approx. 60° out of the membrane plane). These assignments are in close agreement with the ones proposed in a previous LD study at room temperature (Paillotin, G., Verméglio, A. and Breton, J. (1979) Biochim. Biophys. Acta 545, 249–264). The main Q_X excitonic component of P has a broad absorption peaking at 620 nm and it corresponds to dipoles exhibiting the same orientation as those responsible for the 850 nm transition. On the basis of the present LD study and of CD data of chemically oxidized-minus-reduced reaction centers, we proposed that the minor Q_X excitonic component of P is oriented close to the membrane plane and absorbs around 660 nm. The two monomeric bacteriochlorophylls exhibit a positive LD for both their Q_Y transitions (unresolved at 834 nm) and their Q_X transitions (resolved at 600 and 607 nm), indicating that the planes of these molecules are only slightly tilted out of the membrane plane. The two bacteriopheophytins exhibit strong negative LD with identical LD/A values for their Q_Y transitions (resolved at 790 and 805 nm) and small positive LD for their Q_X transitions (resolved at 534 and 544 nm), demonstrating that these two molecules are strongly tilted out of the membrane plane with each of the Q_Y transitions tilted at approx. 50° out of that plane. A comparison of these LD data with the structural model derived from X-ray crystallography (Deisenhofer, J., Epp, O., Miki, K., Huber, R. and Michel, H. (1984) J. Mol. Biol. 180, 385–398) clearly suggests that a good agreement exists between the results of the two techniques under the following conditions: (i) the C-2 symmetry axis of the reaction center runs along the membrane normal; (ii) excitonic coupling is present only in the primary donor special pair; and (iii) the direction of the optical transitions of the monomeric bacteriochlorophylls and of the bacteriopheophytins is not significantly perturbed by the interactions among the pigments. In addition, a carotenoid is detected in the isolated reaction center with an orientation rather perpendicular to the C-2 symmetry axis. Finally, a comparison of these data with similar ones obtained on the bacteriochlorophyll a-containing reaction center of Rhodopseudomonas sphaeroides 241 points towards a geometrical arrangement of the chromophores which is indistinguishable from the one observed in the reaction center of Rps. viridis.     

Covalently linked zinc chlorophyll dimers as a model of a chlorophyllous pair in photosynthetic reaction centers

A heterodimer, where zinc pyropheophorbide-a was linked with zinc pyropheophorbide-d through ethylene glycol diester, was prepared, as well as the corresponding homodimers. The synthetic dimers were complexed with methanol in benzene to give folded dimers by mutual Zn...O(Me)-H...O=C13(1) bonding. Such complexes had furthest red (Qy) absorption bands at longer wavelengths than the monomeric species. These red-shifts were ascribable to excitonic coupling of the Qy transition states in the chlorin pi-pi stacking conformer. In the heterodimeric system, a minor band was observed at the shorter wavelength side of the main Qy band. This observation can be explained by an additional contribution of Qy vibronic state to the exciton-coupled states. Based on the experimental results, a pair of chlorophyll(Chl)-d with Chl-a as well as a Chl-d homopair were proposed as dimers in reaction centers of Chl-d dominating cyanobacteria.     

Theoretical evidence for the coupling of winding to compression in the solution conformation of duplex DNA

Published data indicate linear variation of the duplex winding of DNA with logarithm of the ionic strength. We show that this observation is a natural prediction of a current poly-electrolyte theory if structural coupling between duplex rotation angle (rotation per base pair) and step height (axial translation per base pair) is assumed. We show also that the coupling must then be negative: the positive winding on augmentation of the ionic strength must be accompanied by a decrement of step height, or axial compression. Cation species specificity of winding, CD spectra, and cation binding affinity are discussed in terms of local variation of dielectric shielding, using current accurate data for dielectric constants of aqueous salt solutions. One need not invoke specific complexation of cations, for which there is no present evidence (in aqueous solution), to explain any of these effects. The analysis is limited to low and moderate ionic strengths. Experimental observations at very high ionic strengths are puzzling and possibly conflicting. It is concluded that a need exists for systematic extension of available duplex winding data throughout the region of moderate to very high ionic strengths. It would be useful to extend the region at the low end as well, in order to distinguish clearly the behavior predicted by our model from that predicted by a referee on the basis of an alternate model described herein.     

Calculation of pigment transition energies in the FMO protein

The Fenna-Matthews-Olson (FMO) protein of green sulfur bacteria represents an important model protein for the study of elementary pigment-protein couplings. We have previously used a simple approach [Adolphs and Renger (2006) Biophys J 91:2778-2797] to study the shift in local transition energies (site energies) of the FMO protein of Prosthecochloris aestuarii by charged amino acid residues, assuming a standard protonation pattern of the titratable groups. Recently, we have found strong evidence that besides the charged amino acids also the neutral charge density of the protein is important, by applying a combined quantum chemical/electrostatic approach [Müh et al. (2007) Proc Natl Acad Sci USA, in press]. Here, we extract the essential parts from this sophisticated method to obtain a relatively simple method again. It is shown that the main contribution to the site energy shifts is due to charge density coupling (CDC) between the pigments and their pigment, protein and water surroundings and that polarization effects for qualitative considerations can be approximated by screening the Coulomb coupling by an effective dielectric constant.