Chlorophyll a Radical Ions: A Density Functional Study

A dimeric chlorophyll electron acceptor differentiates type I from type II photosynthetic reaction centers

This research addresses one of the most compelling issues in the field of photosynthesis, namely, the role of the accessory chlorophyll molecules in primary charge separation. Using a combination of empirical and computational methods, we demonstrate that the primary acceptor of photosystem (PS) I is a dimer of accessory and secondary chlorophyll molecules, Chl_2A and Chl_3A, with an asymmetric electron charge density distribution. The incorporation of highly coupled donors and acceptors in PS I allows for extensive delocalization that prolongs the lifetime of the charge-separated state, providing for high quantum efficiency. The discovery of this motif has widespread implications ranging from the evolution of naturally occurring reaction centers to the development of a new generation of highly efficient artificial photosynthetic systems.

Video abstract

Analysis of Photosynthetic Systems and Their Applications with Mathematical and Computational Models

In biological and life science applications, photosynthesis is an important process that involves the absorption and transformation of sunlight into chemical energy. During the photosynthesis process, the light photons are captured by the green chlorophyll pigments in their photosynthetic antennae and further funneled to the reaction center. One of the most important light harvesting complexes that are highly important in the study of photosynthesis is the membrane-attached Fenna–Matthews–Olson (FMO) complex found in the green sulfur bacteria. In this review, we discuss the mathematical formulations and computational modeling of some of the light harvesting complexes including FMO. The most recent research developments in the photosynthetic light harvesting complexes are thoroughly discussed. The theoretical background related to the spectral density, quantum coherence and density functional theory has been elaborated. Furthermore, details about the transfer and excitation of energy in different sites of the FMO complex along with other vital photosynthetic light harvesting complexes have also been provided. Finally, we conclude this review by providing the current and potential applications in environmental science, energy, health and medicine, where such mathematical and computational studies of the photosynthesis and the light harvesting complexes can be readily integrated.

Probing the Electronic Structure of Bacteriochlorophyll Radical Ions-A Theoretical Study of the Effect of Substituents on Hyperfine Parameters

In reaction centers (RCs) of photosynthesis, a light-induced charge separation takes place creating radical cations and anions of the participating cofactors. In photosynthetic bacteria, different bacteriochlorophylls (BChl) are involved in this process. Information about the electronic structure of the BChl radical cations and anions can be obtained by measuring the electron spin density distribution via the electron-nuclear hyperfine interaction using EPR and ENDOR techniques. In this communication, we report isotropic hyperfine coupling constants (hfcs) of the BChl b and g radical cations and anions, calculated by density functional theory, and compare them with the more common radical ions of BChl a and with available experimental data. The observed differences in the computed hyperfine data are discussed in view of a possible distinction between these species by EPR/ENDOR methods. In addition, ¹⁴ N nuclear quadrupole coupling constants (nqcs) computed for BChl a, b, g, and also for Chl a in their charge neutral, radical cation and radical anion states are presented. These nqcs are compared with experimental values obtained by ESEEM spectroscopy on several different radical ions.

Optical absorption and electronic spectra of chlorophylls a and b

We report optical and electronic properties of the two main chlorophylls in green plants, namely, chlorophylls a and b. We estimate the electric moments of these molecules and study absorption spectra of the chlorophylls.

Interactions of copper(ii) and zinc(ii) with chlorophyll: insights from density functional theory studies

The most favored reaction of chlorophyll is computed to be substitution for Cu²⁺ and peripheral chelation for Zn²⁺.

Fine tuning of copper(II)-chlorophyll interactions in organic media. Metalation versus oxidation of the macrocycle

The nature of chlorophyll interactions with copper(II) ions varies considerably in organic solvents, depending on the dominant coordinative form. Besides formation of the metallo tetrapyrrolic complex, Cu(II) ions can cause oxidation of the pigment, reversible or irreversible, which can lead to the destruction of the macrocyclic structure. All these reaction types can be distinguished within a quite narrow range of reaction conditions. The ability to form new metallo derivatives in either metalation or transmetalation reactions is obviously limited by the concentration of the potential oxidant, but can be secured below this level via suitable composition of the reaction system. The decisive factor in the selection of a specific reaction pathway is the presence of a potential ligand that can affect the reactivity of Cu(II) for example by shifting its redox potential. Spectroscopic and electrochemical studies were performed in order to determine the predominant species of Cu(II) in methanol, nitromethane and acetonitrile in the presence of chloride and acetate ions, as well as to assign their appropriate oxidizing ability. This allowed us to estimate the boundary conditions for the electron transfer processes in chlorophyll-Cu(II) systems. Chlorophyll and its free base can undergo both types of electron transfer processes, however, they reveal different susceptibilities that make this class of ligands quite versatile markers in tuning the reactivity of metal ions in solutions.

Density-functional geometry optimization of the 150 000-atom photosystem-I trimer

We present a linear-scaling method based on the use of density-functional theory (DFT) for the system-wide optimization of x-ray structural coordinates and apply it to optimize the 150,000 atoms of the photosystem-I (PS-I) trimer. The method is based on repetitive applications of a multilevel ONIOM procedure using the PW916-31G(d) DFT calculations for the high level and PM3 for the lower level; this method treats all atoms in the structure equivalently, a structure in which the majority of the atoms can be considered as part of some internal "active site." To obtain a realistic single structure, some changes to the original protein model were necessary but these are kept to a minimum in order that the optimized structure most closely resembles the original x-ray one. Optimization has profound effects on the perceived electronic properties of the cofactors, with, e.g., optimization lowering the internal energy of the chlorophylls by on average 53 kcal mol(-1) and eliminates the enormous 115 kcal mol(-1) energy spread depicted by the original x-ray heavy-atom coordinates. A highly precise structure for PS-I results that is suitable for analysis of device function. Significant qualitative features of the structure are also improved such as correction of an error in the stereochemistry of one of the chlorophylls in the "special pair" of the reaction center, as well as the replacement of a water molecule with a metal cation in a critical region on the C3 axis. The method also reveals other unusual features of the structure, leading both to suggestions concerning device functionality and possible mutations between gene sequencing and x-ray structure determination. The optimization scheme is thus shown to augment the molecular modeling schemes that are currently used to add medium-resolution structural information to the raw scattering data in order to obtain atomically resolved structures. System-wide optimization is now a feasible process and its use within protein x-ray data refinement should be considered.

Quantum chemical simulation of excited states of chlorophylls, bacteriochlorophylls and their complexes

The present review describes the use of quantum chemical methods in estimation of structures and electronic transition energies of photosynthetic pigments in vacuum, in solution and imbedded in proteins. Monomeric Mg-porphyrins, chlorophylls and bacteriochlorophylls and their solvent 1:1 and 1:2 complexes were studied. Calculations were performed for Mg-porphyrin, Mg-chlorin, Mg-bacteriochlorin, mesochlorophyll a, chlorophylls a, b, c(1), c(2), c(3), d and bacteriochlorophylls a, b, c, d, e, f, g, h, plus several homologues. Geometries were optimised with PM3, PM3/CISD, PM5, ab initio HF (6-31G*/6-311G**) and density functional B3LYP (6-31G*/6-311G**) methods. Spectroscopic transition energies were calculated with ZINDO/S CIS, PM3 CIS, PM3 CISD, ab initio CIS, time-dependent HF and time-dependent B3LYP methods. Estimates for experimental transition energies were obtained from linear correlations of the calculated transition energies of 1:1 solvent complexes against experimentally recorded solution energies (scaling). According to the calculations in five-coordinated solvent complexes the magnesium atom lies out of the porphyrin plane, while in six-coordinated complexes the porphyrin is nearly planar. Charge densities on magnesium and nitrogen atoms were strongly dependent on the computational method deployed. Several dark states of low oscillator strength below the main Soret band were predicted for solvent complexes and chlorophylls and bacteriochlorophylls in protein environment. Such states, though not yet identified experimentally, might serve as intermediate states for excitation energy transfer in photosynthetic complexes. Q(y), Q(x) and Soret transition energies were found to depend on the orientation of the acetyl group and external pressure. A method to estimate site energies and dimeric interaction energies and to simulate absorption and CD spectra of photosynthetic complexes is described. Simulations for the light harvesting complexes Rhodospirillum molischianum, chlorosomes of Chlorobium tepidum and Chloroflexus aurantiacus, and LHC-II of Spinacia oleracea are presented as examples.

Density Functional Theory Calculations on the Dielectric Constant Dependence of the Oxidation Potential of Chlorophyll: Implication for the High Potential of P680 in Photosystem II†

The primary donor chlorophyll (Chl) of photosystem II (PSII), P680, has an extremely high oxidation redox potential (E(ox)) of approximately 1.2 V, which is essential for photosynthetic water oxidation. The mechanism for achieving a high potential such as that of P680 has been one of the central questions in photosynthesis research. Here, we have examined the dielectric constant (epsilon) dependence of the E(ox) of monomer Chl using density functional theory calculations with the polarizable continuum model. The calculated E(ox) of a model Chl compound exhibited a sharp increase with a decrease in epsilon in the relatively low epsilon region (epsilon < 5). In contrast, in the higher-epsilon region, E(ox) was rather insensitive to epsilon and converged to a constant value at very high epsilon values. This tendency in the high-epsilon region explains the experimental E(ox) values of isolated Chl a that have been observed in a relatively narrow range of 0.74-0.93 V. The E(ox) of Chl in an ideal hydrophobic protein was estimated to be approximately 1.4 V at an epsilon value of 2. This value indicates that Chl in a hydrophobic environment originally has a high E(ox) that is sufficient for oxidizing water (E(ox) = 0.88 V at pH 6). On the basis of the reported X-ray crystallographic structures, the protein environment of P680 in PSII was estimated to be more hydrophobic than that of the primary donors in bacterial reaction centers. It is therefore suggested that the low-dielectric environment around P680 is one of the major factors in its very high E(ox), and thus, introducing nonpolar amino acids into the binding pocket of P680 was an important process in the evolution of PSII.