Excitation wavelength-dependent electron-phonon and electron-vibrational coupling in the CP29 antenna complex of green plants

Electron-phonon and electron-vibrational coupling strengths of a weakly (excitonically) coupled chlorophyll a S1-->S0 transition of the CP29 antenna complex of plant photosystem II were studied by difference fluorescence-line-narrowing spectroscopy at 4.5 K. A strong, almost linear increase of the electron-phonon coupling strength toward longer wavelengths was observed, with Huang-Rhys factors Sph increasing from 0.41+/-0.05 at 680 nm to about 0.66+/-0.07 at 688 nm. The former and latter wavelengths are located close to the peak and on the red edge of the inhomogeneous site distribution function, respectively. The experimentally obtained wavelength dependence of Sph may originate either from an alteration of the electron-phonon coupling strength by the local environment of the fluorescing chromophore and/or from the presence of two isoforms of CP29, which are characterized by different coupling strengths to the protein environment. The one-phonon profile peaks at omegam=22 cm(-1) and is described by an asymmetric function composed of a Gaussian low-energy wing and a Lorentzian high-energy tail with half-widths at half-maximum of 10+/-1 and 60+/-10 cm(-1), respectively. Thirty-nine individual vibrational modes between 90 and 1665 cm(-1) were resolved, and their Huang-Rhys factors were determined, which fall in the range between 0.0004 and 0.032. The broad feature present in the overlap region of phonon and vibrational modes at about 90 cm(-1) is characterized by S=0.048. An integral value of vibrational coupling strengths Svib=0.36+/-0.05 was determined, which is similar to that observed earlier for the trimeric LHC II complex.     

Low-energy chlorophyll states in the CP43 antenna protein complex: simulation of various optical spectra. II

The CP43 protein complex of the core antenna of higher plant photosystem II (PSII) has two quasidegenerate "red" absorption states. It has been shown in the accompanying paper I (Dang, N. C., et al. J. Phys. Chem. B 2008, 112, 9921.) that the site distribution functions (SDFs) of red-states A and B are uncorrelated and the narrow holes are burned in subpopulations of chlorophylls (Chls) from states A and B that are the lowest-energy pigments in their particular CP43 complexes and cannot further transfer energy downhill. In this work, we present the results of a series of Monte Carlo simulations using the 3.0-A structure of the PSII core complex from cyanobacteria (Loll, B., et al. Nature 2005, 303, 1040.) to model absorption, emission, persistent, and transient hole burned (HB) spectra. At the current structural resolution, we found calculated site energies (obtained from INDO/S calculations) to be only suggestive because their values are different for the two monomers of CP43 in the PS II dimer. As a result, to probe the excitonic structure, a simple fitting procedure was employed to optimize Chl site energies from various starting values corresponding to different A/B pigment combinations to provide simultaneously good fits to several types of optical spectra. It is demonstrated that the shape of the calculated absorption, emission, and transient/persistent hole-burned spectra is consistent with experimental data and our model for excitation energy transfer between two quasi-degenerate lowest-E states (A and B) with uncorrelated SDFs discussed in paper I. Calculations revealed that absorption changes observed near 670 nm in the non-line-narrowed persistent HB spectra (assigned to photoconversion involving Chl-protein hydrogen-bonding by Hughes (Biochemistry 2006, 45, 12345.) are most likely the result of nonphotochemical hole-burning (NPHB) accompanied by the redistribution of oscillator strength due to modified excitonic interactions. We argue that a unique redistribution of oscillator strength during the NPHB process helps to assign Chls contributing to the low-energy states. It is demonstrated that the 4.2 K asymmetric triplet-bottleneck (transient) hole is mostly contributed to by both A and B states, with the hole profile described by a subensemble of pigments, which are the lowest-energy pigments (B s- and A s-type) in their complexes. The same lowest-energy Chls contribute to the observed fluorescence spectra. On the basis of our excitonic calculations, the best Chl candidates that contribute to the low-energy A and B states are Chl 44 and Chl 37, respectively.     

Red antenna states of photosystem I from cyanobacteria Synechocystis PCC 6803 and Thermosynechococcus elongatus: single-complex spectroscopy and spectral hole-burning study

Hole-burning and single photosynthetic complex spectroscopy were used to study the excitonic structure and excitation energy-transfer processes of cyanobacterial trimeric Photosystem I (PS I) complexes from Synechocystis PCC 6803 and Thermosynechococcus elongatus at low temperatures. It was shown that individual PS I complexes of Synechocystis PCC 6803 (which have two red antenna states, i.e., C706 and C714) reveal only a broad structureless fluorescence band with a maximum near 720 nm, indicating strong electron-phonon coupling for the lowest energy C714 red state. The absence of zero-phonon lines (ZPLs) belonging to the C706 red state in the emission spectra of individual PS I complexes from Synechocystis PCC 6803 suggests that the C706 and C714 red antenna states of Synechocystis PCC 6803 are connected by efficient energy transfer with a characteristic transfer time of approximately 5 ps. This finding is in agreement with spectral hole-burning data obtained for bulk samples of Synechocystis PCC 6803. The importance of comparing the results of ensemble (spectral hole burning) and single-complex measurements was demonstrated. The presence of narrow ZPLs near 710 nm in addition to the broad fluorescence band at approximately 730 nm in Thermosynechococcus elongatus (Jelezko et al. J. Phys. Chem. B 2000, 104, 8093-8096) has been confirmed. We also demonstrate that high-quality samples obtained by dissolving crystals of PS I of Thermosynechococcus elongatus exhibit stronger absorption in the red antenna region than any samples studied so far by us and other groups.     

Assignment and extracting dynamics from experimentally and theoretically obtained spectroscopic hamiltonians in the complex spectral and classically chaotic regions

An analysis of existing algebraic multiresonance spectroscopic Hamiltonians, derived by fitting to experimental data or from classical canonical or quantum Van Vleck perturbation theory, allows without any significant further classical or quantum calculation the assignment of quantum numbers and motions to states observed in spectra that were previously thought to be irregular or just unassignable. In such cases, inspection of the amplitude and phase of eigenfunctions previously calculated in the Hamiltonians derivation process but now transformed to a reduced dimension semiclassical action-angle representation reveals extremely simple albeit unfamiliar topologies that give quantum numbers by simply counting nodes and phase advances. The topology allows these simple wave functions to be sorted into dynamically different excitation ladders or classes of states which are associated with different regions of phase space. The rungs of these ladders or the states in the classes intersperse in energy causing the spectral complexity. No experimental procedure allows such sorting. The success of the work stems from (1) the qualitative insights of nonlinear dynamics, (2) the conversion of the quantum problem in full dimension to a semiclassical one in reduced dimension by use of a canonical transform that takes advantage of the polyad and other constants of the motion, and (3) the judicious choice of the reduced angle variables to reflect rational ratio resonance frequency conditions. This leads to localization of those semiclassical wave functions, which are affected by the particular resonance. In reverse, the localized appearance of the reduced dimension wave function reveals which resonances govern it and makes sorting visually simple. The success of the work also stems from (4) the revealing use of plots of phase advances as well as the usual densities of the eigenstates for sorting and assignment purposes. Even in classically chaotic regions, organizing trajectories, which correspond to averages over regional phase space structures that need not be computed, can easily be drawn as the structure about which eigenfunction localization takes place. The organizing trajectories when transformed back to the full dimensional configuration space reveal the internal molecular motions. The complexity of the usual quantum stationary and propagated wave functions and associated classical trajectories forbids most often such assignments and sorting. This procedure brings the ability to interpret complex vibrational spectra to a degree previously thought possible only for lower excitations. The new methodology replaces and extends the computationally more difficult parts of a procedure used by the authors that was applied successfully to several molecules during the past few years. The new methodology is applied to DCO, CHBrClF, and the bending of acetylene.     

Assignment and extraction of dynamics of a small molecule with a complex vibrational spectrum: thiophosgene

The dispersed fluorescence spectrum of the ground electronic state of thiophosgene, SCCl2, is analyzed in a very complex region of vibrational excitation, 7000-9000 cm(-1). The final result is that most of the inferred excited vibrational levels are assigned in terms of approximate constants of the motion. Furthermore, each level is associated with a rung on a ladder of quantum states on the basis of common reduced dimension fundamental motions. The resulting ladders cannot be identified by any experimental means, and it is the interspersing in energy of their rungs that makes the spectrum complex even after the process of level separation into polyads. Van Vleck perturbation theory is used to create polyad constants of the motion and a spectroscopic Hamiltonian from a potential fitted to experimental data. The eigen functions of this spectroscopic Hamiltonian are rewritten as semiclassical wave functions and transformed to a representation that allows us to analyze and assign the spectra with no other work other than to utilize concepts from nonlinear dynamics.     

Optical dynamics of the reaction center of photosystem II. A hole-burning and photon-echo study

Photon-echo and transient hole-burning experiments on the P-band of the reaction center of photosystem II are reported which show that the P-band exhibits a large homogeneous width. A photon-echo signal with a decay time of 500 ps is also observed in the same spectral region and assigned to extraneous chlorophyll. These observations are discussed with reference to the bacterial reaction center and photosystem I.     

Immobilization of light-harvesting chlorophyll a/b complex (LHCIIb) studied by surface plasmon field-enhanced fluorescence spectroscopy

The major light-harvesting chlorophyll a/ b complex (LHCIIb) of the photosynthetic apparatus in green plants can be viewed as a protein scaffold binding and positioning a large number of pigment molecules that engage in rapid excitation energy transfer. This property makes LHCIIb potentially interesting as a light harvester (or a model thereof) in photoelectronic applications. Such applications would require the immobilization of LHCIIb (or similar dye-protein complexes) on a solid surface. In this work, the immobilization of recombinant LHCIIb is tested and optimized on functionalized gold surfaces via a histidine 6 tag (His tag) in the protein moiety. Immobilization efficiency and kinetics are analyzed by using surface plasmon resonance (SPR) and surface plasmon field-enhanced fluorescence spectroscopy (SPFS). The latter was also used to assess the integrity of immobilized LHCIIb by recording Chl b-sensitized Chl a emission spectra. Since His tags have been included in a substantial number of recombinant proteins, the immobilization technique developed here for LHCIIb presumably can be extended to a large range of other membrane and water-soluble proteins.     

Mass spectrometric analysis of the interactions between CP12, a chloroplast protein, and metal ions: a possible regulatory role within a PRK/GAPDH/CP12 complex

The small chloroplast protein CP12 plays the role of a protein linker in the assembly process of a PRK/GAPDH/CP12 complex that is involved in CO2 assimilation in photosynthetic organisms. The redox state of CP12 regulates its role as a protein linker. Only the oxidized protein, with two disulfide bonds, is active in complex formation. Several observations indicating that CP12 might bind a metal ion led us to screen the binding of different metal ions on oxidized or reduced CP12 using non-covalent electrospray ionization mass spectrometry (ESI-MS) experiments. The oxidized protein bound specifically Cu2+ and Ni2+ (Kd of 26+/-1 microM and 11+/-1 microM, respectively); other cations such as Fe2+ and Zn2+ did not bind, while cations such as Cd2+ formed non-specific adducts to CP12. Similar results were obtained for metal ions on screening with the reduced CP12. Interestingly, the present results suggest that Cu2+ catalyzes the re-formation of the disulfide bonds of the reduced CP12, leading to recovery of the fully oxidized CP12 that is then able to bind a Cu2+ ion. Finally the high similarity between CP12 and copper chaperones from Arabidopsis thaliana, as judged by hydrophobic cluster analysis, provides additional evidence for the relevance of metal binding for the in vivo situation. The findings that CP12 is able to bind a metal ion, and that Cu2+ catalyzes the oxidation of the thiol groups of CP12, are new characteristics of this protein that may prove to be important in the regulation of the assembly process of the PRK/GAPDH/CP12 complex.     

Excitation energy transfer in the photosystem II core antenna complex CP43 studied by femtosecond visible/visible and visible/mid-infrared pump probe spectroscopy

Excitation energy transfer in the Photosystem II core antenna complex CP43 has been investigated by vis/vis and vis/mid-IR pump-probe spectroscopy with the aim of understanding the relation between the dynamics of energy transfer and the structural arrangement of individual chlorophyll molecules within the protein. Energy transfer was found to occur on time scales of 250 fs, 2-4 ps, and 10-12 ps. The vis/mid-IR difference spectra show that the excitation is initially distributed over chlorophylls located in environments with different polarity, since two 9-keto C=O stretching bleachings, at 1691 and 1677 cm-1, are observable at early delay times. Positive signals in the initial difference spectra around 1750 and 1720 cm-1 indicate the presence of a charge transfer state between strongly interacting chlorophylls. We conclude, both from the spectral behavior in the visible when the annihilation processes are increased and from the vis/mid-IR data, that there are two pigments (one absorbing around 670 nm and one at 683 nm) which are not connected to the other pigments on a time scale faster than 10-20 ps. Since, in the IR, on a 10 ps time scale the population of the 1691 cm-1 mode almost disappears, while the 1677 cm-1 mode is still significantly populated, we can conclude that at least some of the red absorbing pigments are located in a polar environment, possibly forming H-bonds with the surrounding protein.     

Uranyl nitrate complex extraction into TBP/dodecane organic solutions: a molecular dynamics study

Liquid-liquid extraction of uranyl is studied by conducting atomistic molecular dynamics simulation using quantum chemistry calibrated force fields via restrained electrostatic potential fitting of atomic forces. The simulations depict the migration of uranyl nitrate complexes from the aqueous-organic interface into the tri-n-butyl phosphate (TBP)/dodecane organic phase, in the form of UO(2)(NO(3))(2)·H(2)O·2TBP and UO(2)(NO(3))(2)·3TBP. The migration process is characterized by the gradual breaking of all the hydrogen bonds between the complex and the water molecules at the interface. Moreover, our simulation results suggest that the experimentally observed complex UO(2)(NO(3))(2)·2TBP is formed after the migration of the aforementioned complexes into the organic phase by means of a reorganization of the nitrate binding mode from mono to bidentate which removes the excess oxygen atoms bound to uranyl.     

EPR study of electron spin polarization in the photoexcited triplet state of chlorophyll a and b

The photoexcited triplet states of chlorophyll à and b are studied by the EPR method at ≈85 K using modulated light excitation. Both compounds show anomalous EPR line intensities and transient kinetics, indicating electron spin polarization (ESP) in the photoexcited triplet state. EPR studies, using Mg-tetraphenyl porphyrin (MgTPP) dissolved in n-octane show that ESP occurs also in that solvent. It is shown that the zero field splitting (ZFS) parameters of MgTPP depend strongly on the solvent. From the analysis of the data for chlorophyll a and b we evaluate: (1) the population rate constants (k_p); (2) the ratio between the population rate constants (A_p) (p = x, y, z) and, (3) the spin lattice relaxation rate W. In both chlorophylls the in-plane component, x, is predominantly populated and depopulated. The ZFS parameters have been also determined for the above compounds.     

Chlorophyll excitation energies and structural stability of the CP47 antenna of photosystem II: a case study in the first-principles simulation of light-harvesting complexes

Natural photosynthesis relies on light harvesting and excitation energy transfer by specialized pigment–protein complexes. Their structure and the electronic properties of the embedded chromophores define the mechanisms of energy transfer. An important example of a pigment–protein complex is CP47, one of the integral antennae of the oxygen-evolving photosystem II (PSII) that is responsible for efficient excitation energy transfer to the PSII reaction center. The charge-transfer excitation induced among coupled reaction center chromophores resolves into charge separation that initiates the electron transfer cascade driving oxygenic photosynthesis. Mapping the distribution of site energies among the 16 chlorophyll molecules of CP47 is essential for understanding excitation energy transfer and overall antenna function. In this work, we demonstrate a multiscale quantum mechanics/molecular mechanics (QM/MM) approach utilizing full time-dependent density functional theory with modern range-separated functionals to compute for the first time the excitation energies of all CP47 chlorophylls in a complete membrane-embedded cyanobacterial PSII dimer. The results quantify the electrostatic effect of the protein on the site energies of CP47 chlorophylls, providing a high-level quantum chemical excitation profile of CP47 within a complete computational model of “near-native” cyanobacterial PSII. The ranking of site energies and the identity of the most red-shifted chlorophylls (B3, followed by B1) differ from previous hypotheses in the literature and provide an alternative basis for evaluating past approaches and semiempirically fitted sets. Given that a lot of experimental studies on CP47 and other light-harvesting complexes utilize extracted samples, we employ molecular dynamics simulations of isolated CP47 to identify which parts of the polypeptide are most destabilized and which pigments are most perturbed when the antenna complex is extracted from PSII. We demonstrate that large parts of the isolated complex rapidly refold to non-native conformations and that certain pigments (such as chlorophyll B1 and β-carotene h1) are so destabilized that they are probably lost upon extraction of CP47 from PSII. The results suggest that the properties of isolated CP47 are not representative of the native complexed antenna. The insights obtained from CP47 are generalizable, with important implications for the information content of experimental studies on biological light-harvesting antenna systems.

Advanced QM/MM simulations explore the excited states of a photosynthetic light-harvesting antenna in its physiologically complexed state and model the consequences of extraction on conformational and electronic properties.     

Influence of the host on vibronic relaxation: a study of free-base porphin in a series of n-alkanes by photochemical hole-burning

Photochemical hole-burning is used to determine the relaxation times of vibronic bands of the S₁ ← S₀ transition to free-base porphin in different substitutional sites of n-hexane, n-heptane, n-octane and n-decane at 1.6 K. The vibronic relaxation depends strongly on site and host. A correlation between the n-alkane chain length and the vibronic relaxation time is observed.     

Isolation and characterization of a photosystem I-associated antenna (LHC I) and a photosystem I—core complex from the chlorophyll c-containing alga Pleurochloris meiringensis (Xanthophyceae)

A photosystem (PS) I holocomplex was isolated from Pleurochloris meiringensis Vischer (Xanthophyceae) using sucrose density centrifugation. This complex exhibited a fluorescence emission maximum at 715 nm, which is in accordance with the long wavelength emission of whole cells. The complex was further dissociated into a core complex and a light-harvesting protein (LHC I). The core protein contains mainly Chl a and β-carotene, is 8.25 times enriched in P700 and has its main emission maximum at 715 nm. Therefore, the longest wavelength emission of P. meiringensis is due to the PS I core, which is in contrast to higher plants. The LHC I differs from LHC II with regard to its polypeptide pattern as well as its spectral properties. The arrangement of antennae is discussed in relation to the regulation of energy transfer between the photosystems.     

Mechanism of dihydroneopterin aldolase: a molecular dynamics study of the apo enzyme and its product complex

Dihydroneopterin aldolase (DHNA), an enzyme in the pathway that generates folic acid in bacteria, is investigated by a series of molecular dynamics simulations in its free form and complexed with its product, 6-hydroxymethyl-7,8-dihydropterin (HP). The active sites in DHNA are formed at the interface between pairs of protomers in this octameric protein. On the basis of root-mean-square deviation and root-mean-square fluctuation analyses of the trajectories, which take advantage of the presence of eight active sites, flexible regions of the apo protein surrounding the active site are identified and, upon binding HP, show that the active site is rigidified. Specific residues, associated with binding and the catalytic mechanism of DHNA, are associated with these flexible regions, and their interactions with HP account for most of the binding energy. A Principal Component Analysis shows rigidification of DHNA upon HP binding and that only a few modes of motion capture most of the atomic fluctuations in both apo and HP-bound forms. HP is pushed out of the active site in a series of simulations with different restrained positions between HP and DHNA to obtain a view of the exit pathway and energetic barrier to product release. The chosen pathway leads to a minimal disturbance of the system and provides a barrier consistent with the experimentally determined rate of product release. An analysis of the various components that contribute to the exit path energy and entropy provides insight into the energy-entropy compensation for product release.     

Structure and dynamics of a dihydrogen/hydride ansa molybdenocene complex

In contrast to [Cp(2)MoH(3)](+), which is a thermally stable trihydride complex, the ansa-bridged analogue [(eta-C(5)H(4))(2)CMe(2)MoH(H(2))](+) (1) is a thermally labile dihydrogen/hydride complex. Partial deuteration of the hydride ligands allows observation of J(H)(-)(D) = 11.9 Hz in 1-d(1) and 9.9 Hz in 1-d(2) (245 K), indicative of a dihydrogen/hydride structure. There is a slight preference for deuterium to concentrate in the dihydrogen ligand. A rapid dynamic process interchanges the hydride and dihydrogen moieties in complex 1. Low temperature (1)H NMR spectra of 1 give a single hydride resonance, which broadens at very low temperature due to rapid dipole-dipole relaxation (T(1) = 23 ms (750 MHz, 175 K) for the hydride resonance in 1). Low temperature (1)H NMR spectra of 1-d(2) allow the observation of decoalescence at 180 K into two resonances. The bound dihydrogen ligand exhibits hindered rotation with DeltaG(150) = 7.4 kcal/mol, but H atom exchange is still rapid at all accessible temperatures (down to 130 K). Density functional calculations confirm the dihydrogen/hydride structure as the ground state for the molecule and give estimates for the energy of two hydrogen exchange processes in good agreement with experiment. The presence of the C ansa bridge is shown to decrease the ability of the metallocene fragment to donate to the hydrogens, thus stabilizing the (eta(2)-H(2)) unit and modulating the barrier to H(2) rotation.     

Structure and dynamics of the homologous series of alanine peptides: a joint molecular dynamics/NMR study

The phi,psi backbone angle distribution of small homopolymeric model peptides is investigated by a joint molecular dynamics (MD) simulation and heteronuclear NMR study. Combining the accuracy of the measured scalar coupling constants and the atomistic detail of the all-atom MD simulations with explicit solvent, the thermal populations of the peptide conformational states are determined with an uncertainty of <5 %. Trialanine samples mainly ( approximately 90%) a poly-l-proline II helix-like structure, some ( approximately 10%) beta extended structure, but no alphaR helical conformations. No significant change in the distribution of conformers is observed with increasing chain length (Ala(3) to Ala(7)). Trivaline samples all three major conformations significantly. Triglycine samples the four corner regions of the Ramachandran space and exists in a slow conformational equilibrium between the cis and trans conformation of peptide bonds. The backbone angle distribution was also studied for the segment Ala3 surrounded by either three or eight amino acids on both N- and C-termini from a sequence derived from the protein hen egg white lysozyme. While the conformational distribution of the central three alanine residues in the 9mer is similar to that for the small peptides Ala(3)-Ala(7), major differences are found for the 19mer, which significantly (30-40%) samples alphaR helical stuctures.     

Isolation and protein chemical characterization of the B806-866 antenna complex of the green thermophilic bacterium Chloroflexus aurantiacus

The B806-866 antenna complex was isolated from cytoplasmic membranes of the green thermophilic bacterium Chloroflexus aurantiacus. The membranes were treated with 7 M urea at 50 degrees C, the B806-866 antenna complex was solubilized with a mixture of Noni-fjdet P-40 (octylphenoxypolyethoxyethanol (Sigma)) and sodium dodecylsulphate (2:1) and isolated by sucrose density gradient centrifugation. This antenna complex was characterized by reversed-phase chromatography (fast polypeptide and polynucleotide liquid chromatography), amino acid and sequence analyses. The B806-866 antenna of Chloroflexus aurantiacus consists of two polypeptides: the B806-866-alpha and B806-866-beta polypeptides in an apparent stoichiometric ratio of 1:1, which may be equivalent to the structural elementary unit found in the antenna systems of many species of Rhodospirillaceae.     

Unraveling complex three-body photodissociation dynamics of dimethyl sulfoxide: a femtosecond time-resolved spectroscopic study

Photodissociation of dimethyl sulfoxide at 200 nm has been studied using femtosecond time-resolved spectroscopy. The temporal evolutions of the initial state, intermediates, and products (CH3 and SO) were measured by means of fs pump-probe mass-selected multiphoton ionization and laser-induced fluorescence. Femtosecond time-resolved photofragment translational spectroscopy was also employed to measure the CH3 product kinetic energy distributions as a function of reaction time. The ionization experiments revealed that there are at least three major CH3 product components, whereas the fluorescence experiments indicated that two SO product components are present. The combination of experimental and theoretical results suggested a complex multichannel mechanism involving both concerted and stepwise three-body dissociation pathways.