The Fenna-Matthews-Olson (FMO) light-harvesting protein connects the outer antenna system (chlorosome/baseplate) with the reaction center complex in green sulfur bacteria. Since its first structure determination in the mid-70s, this pigment-protein complex has become an important model system to study excitation energy transfer. Recently, an additional bacteriochlorophyll a (the eighth) pigment was discovered in each subunit of this homotrimer. Our structure-based calculations of the optical properties of the FMO protein demonstrate that the eighth pigment is the linker to the baseplate, confirming recent suggestions from crystallographic studies.

Optical line shape theory is combined with a quantum-chemical/electrostatic calculation of the site energies of the 96 chlorophyll a pigments and their excitonic couplings to simulate optical spectra of photosystem I core complexes from Thermosynechococcus elongatus. The absorbance, linear dichroism and circular dichroism spectra, calculated on the basis of the 2.5 A crystal structure, match the experimental data semiquantitatively allowing for a detailed analysis of the pigment-protein interaction. The majority of site energies are determined by multiple interactions with a large number (>20) of amino acid residues, a result which demonstrates the importance of long-range electrostatic interactions. The low-energy exciton states of the antenna are found to be located at a nearest distance of about 25 A from the special pair of the reaction center. The intermediate pigments form a high-energy bridge, the site energies of which are stabilized by a particularly large number (>100) of amino acid residues. The concentration of low energy exciton states in the antenna is larger on the side of the A-branch of the reaction center, implying an asymmetric delivery of excitation energy to the latter. This asymmetry in light-harvesting may provide the key for understanding the asymmetric use of the two branches in primary electron transfer reactions. Experiments are suggested to check for this possibility.

We report a method for the structure-based calculation of the spectral density of the pigment–protein coupling in light-harvesting complexes that combines normal-mode analysis with the charge density coupling (CDC) and transition charge from electrostatic potential (TrEsp) methods for the computation of site energies and excitonic couplings, respectively. The method is applied to the Fenna–Matthews–Olson (FMO) protein in order to investigate the influence of the different parts of the spectral density as well as correlations among these contributions on the energy transfer dynamics and on the temperature-dependent decay of coherences. The fluctuations and correlations in excitonic couplings as well as the correlations between coupling and site energy fluctuations are found to be 1 order of magnitude smaller in amplitude than the site energy fluctuations. Despite considerable amplitudes of that part of the spectral density which contains correlations in site energy fluctuations, the effect of these correlations on the exciton population dynamics and dephasing of coherences is negligible. The inhomogeneous charge distribution of the protein, which causes variations in local pigment–protein coupling constants of the normal modes, is responsible for this effect. It is seen thereby that the same building principle that is used by nature to create an excitation energy funnel in the FMO protein also allows for efficient dissipation of the excitons’ excess energy.

Single-molecule studies that allow to compute pKa values, proton affinities (gas-phase acidity/basicity) and the electrostatic energy of solvation have been performed for a heterogeneous set of 26 organic compounds. Quantum mechanical density functional theory (DFT) using the Becke-half&half and B3LYP functionals on optimized molecular geometries have been carried out to investigate the energetics of gas-phase protonation. The electrostatic contribution to the solvation energies of protonated and deprotonated compounds were calculated by solving the Poisson equation using atomic charges generated by fitting the electrostatic potential derived from the molecular wave functions in vacuum. The combination of gas-phase and electrostatic solvation energies by means of the thermodynamic cycle enabled us to compute pKa values for the 26 compounds, which cover six distinct chemical groups (carboxylic acids, benzoic acids, phenols, imides, pyridines and imidazoles). The computational procedure for determining pKa values is accurate and transferable with a root-mean-square deviation of 0.53 and 0.57 pKa units and a maximum error of 1.0 pKa and 1.3 pKa units for Becke-half&half and B3LYP DFT functionals, respectively.