1
|
Sarngadharan P, Holtkamp Y, Kleinekathöfer U. Protein Effects on the Excitation Energies and Exciton Dynamics of the CP24 Antenna Complex. J Phys Chem B 2024; 128:5201-5217. [PMID: 38756003 PMCID: PMC11145653 DOI: 10.1021/acs.jpcb.4c01637] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2024] [Revised: 04/29/2024] [Accepted: 05/06/2024] [Indexed: 05/18/2024]
Abstract
In this study, the site energy fluctuations, energy transfer dynamics, and some spectroscopic properties of the minor light-harvesting complex CP24 in a membrane environment were determined. For this purpose, a 3 μs-long classical molecular dynamics simulation was performed for the CP24 complex. Furthermore, using the density functional tight binding/molecular mechanics molecular dynamics (DFTB/MM MD) approach, we performed excited state calculations for the chlorophyll a and chlorophyll b molecules in the complex starting from five different positions of the MD trajectory. During the extended simulations, we observed variations in the site energies of the different sets as a result of the fluctuating protein environment. In particular, a water coordination to Chl-b 608 occurred only after about 1 μs in the simulations, demonstrating dynamic changes in the environment of this pigment. From the classical and the DFTB/MM MD simulations, spectral densities and the (time-dependent) Hamiltonian of the complex were determined. Based on these results, three independent strongly coupled chlorophyll clusters were revealed within the complex. In addition, absorption and fluorescence spectra were determined together with the exciton relaxation dynamics, which reasonably well agrees with experimental time scales.
Collapse
Affiliation(s)
- Pooja Sarngadharan
- School of Science, Constructor
University, Campus Ring
1, 28759 Bremen, Germany
| | - Yannick Holtkamp
- School of Science, Constructor
University, Campus Ring
1, 28759 Bremen, Germany
| | | |
Collapse
|
2
|
Ahad S, Lin C, Reppert M. PigmentHunter: A point-and-click application for automated chlorophyll-protein simulations. J Chem Phys 2024; 160:154111. [PMID: 38639311 DOI: 10.1063/5.0198443] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2024] [Accepted: 03/21/2024] [Indexed: 04/20/2024] Open
Abstract
Chlorophyll proteins (CPs) are the workhorses of biological photosynthesis, working together to absorb solar energy, transfer it to chemically active reaction centers, and control the charge-separation process that drives its storage as chemical energy. Yet predicting CP optical and electronic properties remains a serious challenge, driven by the computational difficulty of treating large, electronically coupled molecular pigments embedded in a dynamically structured protein environment. To address this challenge, we introduce here an analysis tool called PigmentHunter, which automates the process of preparing CP structures for molecular dynamics (MD), running short MD simulations on the nanoHUB.org science gateway, and then using electrostatic and steric analysis routines to predict optical absorption, fluorescence, and circular dichroism spectra within a Frenkel exciton model. Inter-pigment couplings are evaluated using point-dipole or transition-charge coupling models, while site energies can be estimated using both electrostatic and ring-deformation approaches. The package is built in a Jupyter Notebook environment, with a point-and-click interface that can be used either to manually prepare individual structures or to batch-process many structures at once. We illustrate PigmentHunter's capabilities with example simulations on spectral line shapes in the light harvesting 2 complex, site energies in the Fenna-Matthews-Olson protein, and ring deformation in photosystems I and II.
Collapse
Affiliation(s)
- S Ahad
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - C Lin
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - M Reppert
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| |
Collapse
|
3
|
Cherepanov DA, Milanovsky GE, Neverov KV, Obukhov YN, Maleeva YV, Aybush AV, Kritsky MS, Nadtochenko VA. Exciton interactions of chlorophyll tetramer in water-soluble chlorophyll-binding protein BoWSCP. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2024; 309:123847. [PMID: 38217986 DOI: 10.1016/j.saa.2024.123847] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/31/2023] [Revised: 01/02/2024] [Accepted: 01/03/2024] [Indexed: 01/15/2024]
Abstract
The exciton interaction of four chlorophyll a (Chl a) molecules in a symmetrical tetrameric complex of the water-soluble chlorophyll-binding protein BoWSCP was analyzed in the pH range of 3-11. Exciton splitting ΔE = 232 ± 2 cm-1 of the Qy band of Chl a into two subcomponents with relative intensities of 78.1 ± 0.7 % and 21.9 ± 0.7 % was determined by a joint decomposition of the absorption and circular dichroism spectra into Gaussian functions. The exciton coupling parameters were calculated based on the BoWSCP atomic structure in three approximations: the point dipole model, the distributed atomic monopoles, and direct ab initio calculations in the TDDFT/PCM approximation. The Coulomb interactions of monomers were calculated within the continuum model using three values of optical permittivity. The models based on the properties of free Chl a in solution suffer from significant errors both in estimating the absolute value of the exciton interaction and in the relative intensity of exciton transitions. Calculations within the TDDFT/PCM approximation reproduce the experimentally determined parameters of the exciton splitting and the relative intensities of the exciton bands. The following factors of pigment-protein and pigment-pigment interactions were examined: deviation of the macrocycle geometry from the planar conformation of free Chl; the formation of hydrogen bonds between the macrocycle and water molecules; the overlap of wave functions of monomers at close distances. The most significant factor is the geometrical deformation of the porphyrin macrocycle, which leads to an increase in the dipole moment of Chl monomer from 5.5 to 6.9 D and to a rotation of the dipole moment by 15° towards the cyclopentane ring. The contributions of resonant charge-transfer states to the wave functions of the Chl dimer were determined and the transition dipole moments of the symmetric and antisymmetric charge-transfer states were estimated.
Collapse
Affiliation(s)
- D A Cherepanov
- N.N. Semenov Federal Research Center for Chemical Physics, Russian Academy of Sciences, 119991 Moscow, Kosygina str., 4, Russian Federation; A.N. Belozersky Institute Of Physico-Chemical Biology, Moscow State University, 119992 Moscow, Leninskye gory, 1b.40, Russian Federation.
| | - G E Milanovsky
- A.N. Belozersky Institute Of Physico-Chemical Biology, Moscow State University, 119992 Moscow, Leninskye gory, 1b.40, Russian Federation
| | - K V Neverov
- A.N. Bach Institute of Biochemistry, Federal Research Center "Fundamentals of Biotechnology", Russian Academy of Sciences", 119071 Moscow, Leninsky prospect, 33b.2, Russian Federation; Faculty of Biology, Moscow State University, 119234 Moscow, Leninskye gory, 1b.12, Russian Federation
| | - Yu N Obukhov
- A.N. Bach Institute of Biochemistry, Federal Research Center "Fundamentals of Biotechnology", Russian Academy of Sciences", 119071 Moscow, Leninsky prospect, 33b.2, Russian Federation
| | - Yu V Maleeva
- Faculty of Biology, Moscow State University, 119234 Moscow, Leninskye gory, 1b.12, Russian Federation
| | - A V Aybush
- N.N. Semenov Federal Research Center for Chemical Physics, Russian Academy of Sciences, 119991 Moscow, Kosygina str., 4, Russian Federation
| | - M S Kritsky
- A.N. Bach Institute of Biochemistry, Federal Research Center "Fundamentals of Biotechnology", Russian Academy of Sciences", 119071 Moscow, Leninsky prospect, 33b.2, Russian Federation
| | - V A Nadtochenko
- N.N. Semenov Federal Research Center for Chemical Physics, Russian Academy of Sciences, 119991 Moscow, Kosygina str., 4, Russian Federation; Department of Chemistry, Moscow State University, 119991 Moscow, Leninskye gory, 1b.3, Russian Federation
| |
Collapse
|
4
|
Gemeinhardt FG, Lahav Y, Schapiro I, Noy D, Müh F, Lindorfer D, Renger T. Short-Range Effects in the Special Pair of Photosystem II Reaction Centers: The Nonconservative Nature of Circular Dichroism. J Phys Chem Lett 2023; 14:11758-11767. [PMID: 38117270 PMCID: PMC10758115 DOI: 10.1021/acs.jpclett.3c02693] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2023] [Revised: 11/29/2023] [Accepted: 12/04/2023] [Indexed: 12/21/2023]
Abstract
Photosystem II reaction centers extract electrons from water, providing the basis of oxygenic life on earth. Among the light-sensitive pigments of the reaction center, a central chlorophyll a dimer, known as the special pair, so far has escaped a complete theoretical characterization of its excited state properties. The close proximity of the special pair pigments gives rise to short-range effects that comprise a coupling between local and charge transfer (CT) excited states as well as other intermolecular quantum effects. Using a multiscale simulation and a diabatization technique, we show that the coupling to CT states is responsible for 45% of the excitonic coupling in the special pair. The other short-range effects cause a nonconservative nature of the circular dichroism spectrum of the reaction center by effectively rotating the electric transition dipole moments of the special pair pigments inverting and strongly enhancing their intrinsic rotational strength.
Collapse
Affiliation(s)
- Felix G. Gemeinhardt
- Institut
für Theoretische Physik, Johannes
Kepler Universität Linz, Altenberger Strasse 69, 4040 Linz, Austria
| | - Yigal Lahav
- Fritz
Haber Center for Molecular Dynamics Research, Institute of Chemistry, Hebrew University of Jerusalem, 9190401 Jerusalem, Israel
- MIGAL
- Galilee Research Institute, S. Industrial Zone, 1101602 Kiryat Shmona, Israel
| | - Igor Schapiro
- Fritz
Haber Center for Molecular Dynamics Research, Institute of Chemistry, Hebrew University of Jerusalem, 9190401 Jerusalem, Israel
| | - Dror Noy
- MIGAL
- Galilee Research Institute, S. Industrial Zone, 1101602 Kiryat Shmona, Israel
- Faculty
of Sciences and Technology, Tel-Hai Academic
College, 1220800 Upper Galilee, Israel
| | - Frank Müh
- Institut
für Theoretische Physik, Johannes
Kepler Universität Linz, Altenberger Strasse 69, 4040 Linz, Austria
| | - Dominik Lindorfer
- Institut
für Theoretische Physik, Johannes
Kepler Universität Linz, Altenberger Strasse 69, 4040 Linz, Austria
| | - Thomas Renger
- Institut
für Theoretische Physik, Johannes
Kepler Universität Linz, Altenberger Strasse 69, 4040 Linz, Austria
| |
Collapse
|
5
|
Ress L, Malý P, Landgraf JB, Lindorfer D, Hofer M, Selby J, Lambert C, Renger T, Brixner T. Time-resolved circular dichroism of excitonic systems: theory and experiment on an exemplary squaraine polymer. Chem Sci 2023; 14:9328-9349. [PMID: 37712031 PMCID: PMC10498725 DOI: 10.1039/d3sc01674a] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Accepted: 07/28/2023] [Indexed: 09/16/2023] Open
Abstract
Experimental and theoretical foundations for femtosecond time-resolved circular dichroism (TRCD) spectroscopy of excitonic systems are presented. In this method, the system is pumped with linearly polarized light and the signal is defined as the difference between the transient absorption spectrum probed with left and with right circularly polarized light. We present a new experimental setup with a polarization grating as key element to generate circularly polarized pulses. Herein the positive (negative) first order of the diffracted light is left-(right-)circularly polarized and serves as a probe pulse in a TRCD experiment. The grating is capable of transferring ultrashort broadband pulses ranging from 470 nm to 720 nm into two separate beams with opposite ellipticity. By applying a specific chopping scheme we can switch between left and right circular polarizations and detect transient absorption (TA) and TRCD spectra on a shot-to-shot basis simultaneously. We perform experiments on a squaraine polymer, investigating excitonic dynamics, and we develop a general theory for TRCD experiments of excitonically coupled systems that we then apply to describe the experimental data in this particular example. At a magic angle of 54.7° between the pump-pulse polarization and the propagation direction of the probe pulse, the TRCD and TA signals become particularly simple to analyze, since the orientational average over random orientations of complexes factorizes into that of the interaction with the pump and the probe pulse, and the intrinsic electric quadrupole contributions to the TRCD signal average to zero for isotropic samples. Application of exciton theory to linear absorption and to linear circular dichroism spectra of squaraine polymers reveals the presence of two fractions of polymer conformations, a dominant helical conformation with close interpigment distances that are suggested to lead to short-range contributions to site energy shifts and excitonic couplings of the squaraine molecules, and a fraction of unfolded random coils. Theory demonstrates that TRCD spectra of selectively excited helices can resolve state populations that are practically invisible in TA spectroscopy due to the small dipole strength of these states. A qualitative interpretation of TRCD and TA spectra in the spectral window investigated experimentally is offered. The 1 ps time component found in these spectra is related to the slow part of exciton relaxation obtained between states of the helix in the low-energy half of the exciton manifold. The dominant 140 ps time constant reflects the decay of excited states to the electronic ground state.
Collapse
Affiliation(s)
- Lea Ress
- Institut für Physikalische und Theoretische Chemie, Universität Würzburg Am Hubland 97074 Würzburg Germany
| | - Pavel Malý
- Institut für Physikalische und Theoretische Chemie, Universität Würzburg Am Hubland 97074 Würzburg Germany
- Faculty of Mathematics and Physics, Charles University Ke Karlovu 5 121 16 Praha 2 Czech Republic
| | - Jann B Landgraf
- Institut für Physikalische und Theoretische Chemie, Universität Würzburg Am Hubland 97074 Würzburg Germany
- Freiburg Center for Interactive Materials and Bioinspired Technologies (FIT), Universität Freiburg Georges-Köhler-Allee 105 79110 Freiburg Germany
| | - Dominik Lindorfer
- Institut für Theoretische Physik, Johannes Kepler Universität Linz Altenberger Str. 69 4040 Linz Austria
| | - Michael Hofer
- Institut für Theoretische Physik, Johannes Kepler Universität Linz Altenberger Str. 69 4040 Linz Austria
| | - Joshua Selby
- Institut für Organische Chemie, Universität Würzburg Am Hubland 97074 Würzburg Germany
| | - Christoph Lambert
- Institut für Organische Chemie, Universität Würzburg Am Hubland 97074 Würzburg Germany
| | - Thomas Renger
- Institut für Theoretische Physik, Johannes Kepler Universität Linz Altenberger Str. 69 4040 Linz Austria
| | - Tobias Brixner
- Institut für Physikalische und Theoretische Chemie, Universität Würzburg Am Hubland 97074 Würzburg Germany
- Center for Nanosystems Chemistry (CNC), Universität Würzburg Theodor-Boveri-Weg 97074 Würzburg Germany
| |
Collapse
|
6
|
Kim Y, Mitchell Z, Lawrence J, Morozov D, Savikhin S, Slipchenko LV. Predicting Mutation-Induced Changes in the Electronic Properties of Photosynthetic Proteins from First Principles: The Fenna-Matthews-Olson Complex Example. J Phys Chem Lett 2023; 14:7038-7044. [PMID: 37524046 DOI: 10.1021/acs.jpclett.3c01461] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/02/2023]
Abstract
Multiscale molecular modeling is utilized to predict optical absorption and circular dichroism spectra of two single-point mutants of the Fenna-Matthews-Olson photosynthetic pigment-protein complex. The modeling approach combines classical molecular dynamics simulations with structural refinement of photosynthetic pigments and calculations of their excited states in a polarizable protein environment. The only experimental input to the modeling protocol is the X-ray structure of the wild-type protein. The first-principles modeling reproduces changes in the experimental optical spectra of the considered mutants, Y16F and Q198V. Interestingly, the Q198V mutation has a negligible effect on the electronic properties of the targeted bacteriochlorophyll a pigment. Instead, the electronic properties of several other pigments respond to this mutation. The molecular modeling demonstrates that a single-point mutation can induce long-range effects on the protein structure, while extensive structural changes near a pigment do not necessarily lead to significant changes in the electronic properties of that pigment.
Collapse
Affiliation(s)
- Yongbin Kim
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| | - Zach Mitchell
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
| | - Jack Lawrence
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| | - Dmitry Morozov
- Nanoscience Center and Department of Chemistry, University of Jyväskylä, P.O. Box 35, 40014 Jyväskylä, Finland
| | - Sergei Savikhin
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
| | - Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| |
Collapse
|
7
|
Delgado F, Enríquez M. Quantum Entanglement and State-Transference in Fenna-Matthews-Olson Complexes: A Post-Experimental Simulation Analysis in the Computational Biology Domain. Int J Mol Sci 2023; 24:10862. [PMID: 37446061 DOI: 10.3390/ijms241310862] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2023] [Revised: 06/13/2023] [Accepted: 06/27/2023] [Indexed: 07/15/2023] Open
Abstract
Fenna-Mathews-Olson complexes participate in the photosynthetic process of Sulfur Green Bacteria. These biological subsystems exhibit quantum features which possibly are responsible for their high efficiency; the latter may comprise multipartite entanglement and the apparent tunnelling of the initial quantum state. At first, to study these aspects, a multidisciplinary approach including experimental biology, spectroscopy, physics, and math modelling is required. Then, a global computer modelling analysis is achieved in the computational biology domain. The current work implements the Hierarchical Equations of Motion to numerically solve the open quantum system problem regarding this complex. The time-evolved states obtained with this method are then analysed under several measures of entanglement, some of them already proposed in the literature. However, for the first time, the maximum overlap with respect to the closest separable state is employed. This authentic multipartite entanglement measure provides information on the correlations, not only based on the system bipartitions as in the usual analysis. Our study has led us to note a different view of FMO multipartite entanglement as tiny contributions to the global entanglement suggested by other more basic measurements. Additionally, in another related trend, the initial state, considered as a Förster Resonance Energy Transfer, is tracked using a novel approach, considering how it could be followed under the fidelity measure on all possible permutations of the FMO subsystems through its dynamical evolution by observing the tunnelling in the most probable locations. Both analyses demanded significant computational work, making for a clear example of the complexity required in computational biology.
Collapse
Affiliation(s)
- Francisco Delgado
- School of Engineering and Sciences, Tecnologico de Monterrey, Atizapan 52926, Mexico
| | - Marco Enríquez
- School of Engineering and Sciences, Tecnologico de Monterrey, Santa Fe 01389, Mexico
| |
Collapse
|
8
|
Reiter S, Kiss FL, Hauer J, de Vivie-Riedle R. Thermal site energy fluctuations in photosystem I: new insights from MD/QM/MM calculations. Chem Sci 2023; 14:3117-3131. [PMID: 36970098 PMCID: PMC10034153 DOI: 10.1039/d2sc06160k] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2022] [Accepted: 02/03/2023] [Indexed: 02/09/2023] Open
Abstract
Cyanobacterial photosystem I (PSI) is one of the most efficient photosynthetic machineries found in nature. Due to the large scale and complexity of the system, the energy transfer mechanism from the antenna complex to the reaction center is still not fully understood. A central element is the accurate evaluation of the individual chlorophyll excitation energies (site energies). Such an evaluation must include a detailed treatment of site specific environmental influences on structural and electrostatic properties, but also their evolution in the temporal domain, because of the dynamic nature of the energy transfer process. In this work, we calculate the site energies of all 96 chlorophylls in a membrane-embedded model of PSI. The employed hybrid QM/MM approach using the multireference DFT/MRCI method in the QM region allows to obtain accurate site energies under explicit consideration of the natural environment. We identify energy traps and barriers in the antenna complex and discuss their implications for energy transfer to the reaction center. Going beyond previous studies, our model also accounts for the molecular dynamics of the full trimeric PSI complex. Via statistical analysis we show that the thermal fluctuations of single chlorophylls prevent the formation of a single prominent energy funnel within the antenna complex. These findings are also supported by a dipole exciton model. We conclude that energy transfer pathways may form only transiently at physiological temperatures, as thermal fluctuations overcome energy barriers. The set of site energies provided in this work sets the stage for theoretical and experimental studies on the highly efficient energy transfer mechanisms in PSI.
Collapse
Affiliation(s)
- Sebastian Reiter
- Department of Chemistry, Ludwig-Maximilians-Universität München Butenandtstr. 11 81377 Munich Germany
| | - Ferdinand L Kiss
- Department of Chemistry, Ludwig-Maximilians-Universität München Butenandtstr. 11 81377 Munich Germany
| | - Jürgen Hauer
- Department of Chemistry, Technical University of Munich Lichtenbergstr. 4, Garching 85747 Germany
| | - Regina de Vivie-Riedle
- Department of Chemistry, Ludwig-Maximilians-Universität München Butenandtstr. 11 81377 Munich Germany
| |
Collapse
|
9
|
Abstract
Biological pigment-protein complexes (PPCs) exhibit a remarkable ability to tune the optical properties of biological excitons (bioexcitons) through specific pigment-protein interactions. While such fine-tuning allows natural systems (e.g., photosynthetic proteins) to carry out their native functions with near-optimal performance, native function itself is often suboptimal for applications such as biofuel production or quantum technology development. This perspective offers a look at near-term prospects for the rational reoptimization of PPC bioexcitons for new functions using site-directed mutagenesis. The primary focus is on the "structure-spectrum" challenge of understanding the relationships between structural features and spectroscopic properties. While recent examples demonstrate that site-directed mutagenesis can be used to tune nearly all key bioexciton parameters (e.g., site energies, interpigment couplings, and electronic-vibrational interactions), critical challenges remain before we achieve truly rational design of bioexciton properties.
Collapse
Affiliation(s)
- Mike Reppert
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, United States
| |
Collapse
|
10
|
Saito K, Mitsuhashi K, Tamura H, Ishikita H. Quantum mechanical analysis of excitation energy transfer couplings in photosystem II. Biophys J 2023; 122:470-483. [PMID: 36609140 PMCID: PMC9941724 DOI: 10.1016/j.bpj.2023.01.002] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2022] [Revised: 11/21/2022] [Accepted: 01/04/2023] [Indexed: 01/07/2023] Open
Abstract
We evaluated excitation energy transfer (EET) coupling (J) between all pairs of chlorophylls (Chls) and pheophytins (Pheos) in the protein environment of photosystem II based on the time-dependent density functional theory with a quantum mechanical/molecular mechanics approach. In the reaction center, the EET coupling between Chls PD1 and PD2 is weaker (|J(PD1/PD2)| = 79 cm-1), irrespective of a short edge-to-edge distance of 3.6 Å (Mg-to-Mg distance of 8.1 Å), than the couplings between PD1 and the accessory ChlD1 (|J(PD1/ChlD2)| = 104 cm-1) and between PD2 and ChlD2 (|J(PD2/ChlD1)| = 101 cm-1), suggesting that PD1 and PD2 are two monomeric Chls rather than a "special pair". There exist strongly coupled Chl pairs (|J| > ∼100 cm-1) in the CP47 and CP43 core antennas, which may be candidates for the red-shifted Chls observed in spectroscopic studies. In CP47 and CP43, Chls ligated to CP47-His26 and CP43-His56, which are located in the middle layer of the thylakoid membrane, play a role in the "hub" that mediates the EET from the lumenal to stromal layers. In the stromal layer, Chls ligated to CP47-His466, CP43-His441, and CP43-His444 mediate the EET from CP47 to ChlD2/PheoD2 and from CP43 to ChlD1/PheoD1 in the reaction center. Thus, the excitation energy from both CP47 and CP43 can always be utilized for the charge-separation reaction in the reaction center.
Collapse
Affiliation(s)
- Keisuke Saito
- Department of Applied Chemistry, The University of Tokyo, Bunkyo-ku, Tokyo, Japan; Research Center for Advanced Science and Technology, The University of Tokyo, Meguro-ku, Tokyo, Japan.
| | - Koji Mitsuhashi
- Department of Applied Chemistry, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
| | - Hiroyuki Tamura
- Department of Applied Chemistry, The University of Tokyo, Bunkyo-ku, Tokyo, Japan; Research Center for Advanced Science and Technology, The University of Tokyo, Meguro-ku, Tokyo, Japan.
| | - Hiroshi Ishikita
- Department of Applied Chemistry, The University of Tokyo, Bunkyo-ku, Tokyo, Japan; Research Center for Advanced Science and Technology, The University of Tokyo, Meguro-ku, Tokyo, Japan.
| |
Collapse
|
11
|
Cignoni E, Cupellini L, Mennucci B. Machine Learning Exciton Hamiltonians in Light-Harvesting Complexes. J Chem Theory Comput 2023; 19:965-977. [PMID: 36701385 PMCID: PMC9933434 DOI: 10.1021/acs.jctc.2c01044] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
We propose a machine learning (ML)-based strategy for an inexpensive calculation of excitonic properties of light-harvesting complexes (LHCs). The strategy uses classical molecular dynamics simulations of LHCs in their natural environment in combination with ML prediction of the excitonic Hamiltonian of the embedded aggregate of pigments. The proposed ML model can reproduce the effects of geometrical fluctuations together with those due to electrostatic and polarization interactions between the pigments and the protein. The training is performed on the chlorophylls of the major LHC of plants, but we demonstrate that the model is able to extrapolate well beyond the initial training set. Moreover, the accuracy in predicting the effects of the environment is tested on the simulation of the small changes observed in the absorption spectra of the wild-type and a mutant of a minor LHC.
Collapse
|
12
|
Erić V, Li X, Dsouza L, Frehan SK, Huijser A, Holzwarth AR, Buda F, Sevink GJA, de Groot HJM, Jansen TLC. Manifestation of Hydrogen Bonding and Exciton Delocalization on the Absorption and Two-Dimensional Electronic Spectra of Chlorosomes. J Phys Chem B 2023; 127:1097-1109. [PMID: 36696537 PMCID: PMC9923760 DOI: 10.1021/acs.jpcb.2c07143] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Chlorosomes are supramolecular aggregates that contain thousands of bacteriochlorophyll molecules. They perform the most efficient ultrafast excitation energy transfer of all natural light-harvesting complexes. Their broad absorption band optimizes light capture. In this study, we identify the microscopic sources of the disorder causing the spectral width and reveal how it affects the excited state properties and the optical response of the system. We combine molecular dynamics, quantum chemical calculations, and response function calculations to achieve this goal. The predicted linear and two-dimensional electronic spectra are found to compare well with experimental data reproducing all key spectral features. Our analysis of the microscopic model reveals the interplay of static and dynamic disorder from the molecular perspective. We find that hydrogen bonding motifs are essential for a correct description of the spectral line shape. Furthermore, we find that exciton delocalization over tens to hundreds of molecules is consistent with the two-dimensional electronic spectra.
Collapse
Affiliation(s)
- Vesna Erić
- University
of Groningen, Zernike Institute
for Advanced Materials, 9747
AG Groningen, The Netherlands
| | - Xinmeng Li
- Department
of Chemistry and Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, Sem Sælands vei 26, 0315 Oslo, Norway
| | - Lolita Dsouza
- Leiden
Institute of Chemistry, Leiden University, Einsteinweg 55, 2300 RA Leiden, The Netherlands
| | - Sean K. Frehan
- MESA+
Institute for Nanotechnology, University
of Twente, Drienerlolaan 5, 7522 NB Enschede, The Netherlands
| | - Annemarie Huijser
- MESA+
Institute for Nanotechnology, University
of Twente, Drienerlolaan 5, 7522 NB Enschede, The Netherlands
| | - Alfred R. Holzwarth
- Department
of Biophysical Chemistry, Max Planck Institute
for Chemical Energy Conversion, Stiftstraße 34-36, 45470 Mülheim, Germany
| | - Francesco Buda
- Leiden
Institute of Chemistry, Leiden University, Einsteinweg 55, 2300 RA Leiden, The Netherlands
| | - G. J. Agur Sevink
- Leiden
Institute of Chemistry, Leiden University, Einsteinweg 55, 2300 RA Leiden, The Netherlands
| | - Huub J. M. de Groot
- Leiden
Institute of Chemistry, Leiden University, Einsteinweg 55, 2300 RA Leiden, The Netherlands
| | - Thomas L. C. Jansen
- University
of Groningen, Zernike Institute
for Advanced Materials, 9747
AG Groningen, The Netherlands,
| |
Collapse
|
13
|
Sarngadharan P, Maity S, Kleinekathöfer U. Spectral densities and absorption spectra of the core antenna complex CP43 from photosystem II. J Chem Phys 2022; 156:215101. [DOI: 10.1063/5.0091005] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Besides absorbing light, the core antenna complex CP43 of photosystem II is of great importance in transferring excitation energy from the antenna complexes to the reaction center. Excitation energies, spectral densities, and linear absorption spectra of the complex have been evaluated by a multiscale approach. In this scheme, quantum mechanics/molecular mechanics molecular dynamics simulations are performed employing the parameterized density functional tight binding (DFTB) while the time-dependent long-range-corrected DFTB scheme is applied for the excited state calculations. The obtained average spectral density of the CP43 complex shows a very good agreement with experimental results. Moreover, the excitonic Hamiltonian of the system along with the computed site-dependent spectral densities was used to determine the linear absorption. While a Redfield-like approximation has severe shortcomings in dealing with the CP43 complex due to quasi-degenerate states, the non-Markovian full second-order cumulant expansion formalism is able to overcome the drawbacks. Linear absorption spectra were obtained, which show a good agreement with the experimental counterparts at different temperatures. This study once more emphasizes that by combining diverse techniques from the areas of molecular dynamics simulations, quantum chemistry, and open quantum systems, it is possible to obtain first-principle results for photosynthetic complexes, which are in accord with experimental findings.
Collapse
Affiliation(s)
- Pooja Sarngadharan
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Sayan Maity
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Ulrich Kleinekathöfer
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| |
Collapse
|
14
|
Cignoni E, Cupellini L, Mennucci B. A fast method for electronic couplings in embedded multichromophoric systems. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:304004. [PMID: 35552268 DOI: 10.1088/1361-648x/ac6f3c] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/17/2022] [Accepted: 05/12/2022] [Indexed: 06/15/2023]
Abstract
Electronic couplings are key to understanding exciton delocalization and transport in natural and artificial light harvesting processes. We develop a method to compute couplings in multichromophoric aggregates embedded in complex environments without running expensive quantum chemical calculations. We use a transition charge approximation to represent the quantum mechanical transition densities of the chromophores and an atomistic and polarizable classical model to describe the environment atoms. We extend our framework to estimate transition charges directly from the chromophore geometry, i.e., bypassing completely the quantum mechanical calculations using a regression approach. The method allows to rapidly compute accurate couplings for a large number of geometries along molecular dynamics trajectories.
Collapse
Affiliation(s)
- Edoardo Cignoni
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124, Pisa, Italy
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124, Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124, Pisa, Italy
| |
Collapse
|
15
|
Friedl C, Fedorov DG, Renger T. Towards a quantitative description of excitonic couplings in photosynthetic pigment-protein complexes: quantum chemistry driven multiscale approaches. Phys Chem Chem Phys 2022; 24:5014-5038. [PMID: 35142765 PMCID: PMC8865841 DOI: 10.1039/d1cp03566e] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
A structure-based quantitative calculation of excitonic couplings between photosynthetic pigments has to describe the dynamical polarization of the protein/solvent environment of the pigments, giving rise to reaction field and screening effects. Here, this challenging problem is approached by combining the fragment molecular orbital (FMO) method with the polarizable continuum model (PCM). The method is applied to compute excitonic couplings between chlorophyll a (Chl a) pigments of the water-soluble chlorophyll-binding protein (WSCP). By calibrating the vacuum dipole strength of the 0–0 transition of the Chl a chromophores according to experimental data, an excellent agreement between calculated and experimental linear absorption and circular dichroism spectra of WSCP is obtained. The effect of the mutual polarization of the pigment ground states is calculated to be very small. The simple Poisson-Transition-charge-from-Electrostatic-potential (Poisson-TrEsp) method is found to accurately describe the screening part of the excitonic coupling, obtained with FMO/PCM. Taking into account that the reaction field effects of the latter method can be described by a scalar constant leads to an improvement of Poisson-TrEsp that is expected to provide the basis for simple and realistic calculations of optical spectra and energy transfer in photosynthetic light-harvesting complexes. In addition, we present an expression for the estimation of Huang–Rhys factors of high-frequency pigment vibrations from experimental fluorescence line-narrowing spectra that takes into account the redistribution of oscillator strength by the interpigment excitonic coupling. Application to WSCP results in corrected Huang–Rhys factors that are less than one third of the original values obtained by the standard electronic two-state analysis that neglects the above redistribution. These factors are important for the estimation of the dipole strength of the 0–0 transition of the chromophores and for the development of calculation schemes for the spectral density of the exciton-vibrational coupling. The importance of reaction field and screening effects on the excitonic couplings is demonstrated, and from quantum-chemical calculations a single scaling factor is derived that can be used to improve simple models based on the Poisson equation.![]()
Collapse
Affiliation(s)
- Christian Friedl
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria.
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba, 305-8568, Japan.
| | - Thomas Renger
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria.
| |
Collapse
|
16
|
Sirohiwal A, Pantazis DA. Electrostatic profiling of photosynthetic pigments: implications for directed spectral tuning. Phys Chem Chem Phys 2021; 23:24677-24684. [PMID: 34708851 DOI: 10.1039/d1cp02580e] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Photosynthetic pigment-protein complexes harvest solar energy with a high quantum efficiency. Protein scaffolds are known to tune the spectral properties of embedded pigments principally through structured electrostatic environments. Although the physical nature of electrostatic tuning is straightforward, the precise spatial principles of electrostatic preorganization remain poorly explored for different protein matrices and incompletely characterized with respect to the intrinsic properties of different photosynthetic pigments. In this work, we study the electronic structure features associated with the lowest excited state of a series of eight naturally occurring (bacterio)chlorophylls and pheophytins to describe the precise topological differences in electrostatic potentials and hence determine intrinsic differences in the expected mode and impact of electrostatic tuning. The difference electrostatic potentials between the ground and first excited states are used as fingerprints. Both the spatial profile and the propensity for spectral tuning are found to be unique for each pigment, indicating spatially and directionally distinct modes of electrostatic tuning. The results define a specific partitioning of the protein matrix around each pigment as an aid to identify regions with a maximal impact on spectral tuning and have direct implications for dimensionality reduction in protein design and engineering. Thus, a quantum mechanical basis is provided for understanding, predicting, and ultimately designing sequence-modified or pigment-exchanged biological systems, as suggested for selected examples of pigment-reconstituted proteins.
Collapse
Affiliation(s)
- Abhishek Sirohiwal
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Dimitrios A Pantazis
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| |
Collapse
|
17
|
Abstract
Numerous linear and non-linear spectroscopic techniques have been developed to elucidate structural and functional information of complex systems ranging from natural systems, such as proteins and light-harvesting systems, to synthetic systems, such as solar cell materials and light-emitting diodes. The obtained experimental data can be challenging to interpret due to the complexity and potential overlapping spectral signatures. Therefore, computational spectroscopy plays a crucial role in the interpretation and understanding of spectral observables of complex systems. Computational modeling of various spectroscopic techniques has seen significant developments in the past decade, when it comes to the systems that can be addressed, the size and complexity of the sample types, the accuracy of the methods, and the spectroscopic techniques that can be addressed. In this Perspective, I will review the computational spectroscopy methods that have been developed and applied for infrared and visible spectroscopies in the condensed phase. I will discuss some of the questions that this has allowed answering. Finally, I will discuss current and future challenges and how these may be addressed.
Collapse
Affiliation(s)
- Thomas L C Jansen
- Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| |
Collapse
|
18
|
Srivastava A, Ahad S, Wat JH, Reppert M. Accurate prediction of mutation-induced frequency shifts in chlorophyll proteins with a simple electrostatic model. J Chem Phys 2021; 155:151102. [PMID: 34686046 DOI: 10.1063/5.0064567] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Photosynthetic pigment-protein complexes control local chlorophyll (Chl) transition frequencies through a variety of electrostatic and steric forces. Site-directed mutations can modify this local spectroscopic tuning, providing critical insight into native photosynthetic functions and offering the tantalizing prospect of creating rationally designed Chl proteins with customized optical properties. Unfortunately, at present, no proven methods exist for reliably predicting mutation-induced frequency shifts in advance, limiting the method's utility for quantitative applications. Here, we address this challenge by constructing a series of point mutants in the water-soluble chlorophyll protein of Lepidium virginicum and using them to test the reliability of a simple computational protocol for mutation-induced site energy shifts. The protocol uses molecular dynamics to prepare mutant protein structures and the charge density coupling model of Adolphs et al. [Photosynth. Res. 95, 197-209 (2008)] for site energy prediction; a graphical interface that implements the protocol automatically is published online at http://nanohub.org/tools/pigmenthunter. With the exception of a single outlier (presumably due to unexpected structural changes), we find that the calculated frequency shifts match the experiment remarkably well, with an average error of 1.6 nm over a 9 nm spread in wavelengths. We anticipate that the accuracy of the method can be improved in the future with more advanced sampling of mutant protein structures.
Collapse
Affiliation(s)
- Amit Srivastava
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - Safa Ahad
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - Jacob H Wat
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - Mike Reppert
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| |
Collapse
|
19
|
Huai Z, Tong Z, Mei Y, Mo Y. Theoretical Study of the Spectral Differences of the Fenna-Matthews-Olson Protein from Different Species and Their Mutants. J Phys Chem B 2021; 125:8313-8324. [PMID: 34314175 DOI: 10.1021/acs.jpcb.1c01686] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The structural basis for the spectral differences between the Fenna-Matthews-Olson (FMO) proteins from Chlorobaculum tepidum (C. tepidum) and Prosthecochloris aestuarii 2K (P. aestuarii) is yet to be fully understood. Mutation-induced perturbation to the exciton structure and the optical spectra of the complex provide a suitable means to investigate the critical role played by the protein scaffold. In this work, we have performed quantum-mechanics/molecular-mechanics calculations over the molecular dynamics simulation trajectories with the polarized protein-specific charge scheme for both wild-type FMOs and two mutants. Our result reveals that a single-point mutation in the vicinity of BChl 6, namely, W183F of C. tepidum, significantly affects the absorption spectrum, resulting in a switch of the absorption spectrum from type 2, for which the 806 nm band is more pronounced than the 815 nm band, to type 1, for which the 815 nm band is pronounced. Our observations agree with the single-point mutation experiments reported by Saer et al. (Biochim. Biophys. Acta, Bioenerg. 2017, 1858, 288-296) and Khmelnitskiy et al. (J. Phys. Chem. Lett. 2018, 9, 3378-3386). In contrast, the absorption spectrum of the P. aestuarii experiences the opposite transition (from type 1 to type 2) upon the same mutation. Furthermore, by comparing the contributions of individual pigments to the spectra in the wild type and its mutant, we find that a single-point mutation near BChl 6 not only induces changes in excitation energy of BChl 6 per se but also affects the excitonic structures of the neighboring BChls 5 and 7 through strong interpigment electronic couplings, resulting in a significant change in the absorption spectra.
Collapse
Affiliation(s)
- Zhe Huai
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| | - Zhengqing Tong
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China.,Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| | - Yan Mo
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China.,Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| |
Collapse
|
20
|
Gelin MF, Borrelli R. Simulation of Nonlinear Femtosecond Signals at Finite Temperature via a Thermo Field Dynamics-Tensor Train Method: General Theory and Application to Time- and Frequency-Resolved Fluorescence of the Fenna-Matthews-Olson Complex. J Chem Theory Comput 2021; 17:4316-4331. [PMID: 34076412 DOI: 10.1021/acs.jctc.1c00158] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Addressing needs of contemporary nonlinear femtosecond optical spectroscopy, we have developed a fully quantum, numerically accurate wave function-based approach for the calculation of third-order spectroscopic signals of polyatomic molecules and molecular aggregates at finite temperature. The systems are described by multimode nonadiabatic vibronic-coupling Hamiltonians, in which diagonal terms are treated in harmonic approximation, while off-diagonal interstate couplings are assumed to be coordinate independent. The approach is based on the Thermo Field Dynamics (TFD) representation of quantum mechanics and tensor-train (TT) machinery for efficient numerical simulation of quantum evolution of systems with many degrees of freedom. The developed TFD-TT approach is applied to the calculation of time- and frequency-resolved fluorescence spectra of the Fenna-Matthews-Olson (FMO) antenna complex at room temperature taking into account finite time-frequency resolution in fluorescence detection, orientational averaging, and static disorder.
Collapse
Affiliation(s)
- Maxim F Gelin
- School of Sciences, Hangzhou Dianzi University, Hangzhou 310018, China
| | | |
Collapse
|
21
|
Renger T. Semiclassical Modified Redfield and Generalized Förster Theories of Exciton Relaxation/Transfer in Light-Harvesting Complexes: The Quest for the Principle of Detailed Balance. J Phys Chem B 2021; 125:6406-6416. [PMID: 34126008 PMCID: PMC8237266 DOI: 10.1021/acs.jpcb.1c01479] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
A conceptual problem
of transfer theories that use a semiclassical
description of the electron-vibrational coupling is the neglect of
the correlation between momenta and coordinates of nuclei. In the
Redfield theory of exciton relaxation, this neglect leads to a violation
of the principle of detailed balance; equal “uphill”
and “downhill” transfer rate constants are obtained.
Here, we investigate how this result depends on nuclear reorganization
effects, neglected in Redfield but taken into account in the modified
Redfield theory. These reorganization effects, resulting from a partial
localization of excited states, are found to promote a preferential
“downhill” relaxation of excitation energy. However,
for realistic spectral densities of light-harvesting antennae in photosynthesis,
the reorganization effects are too small to compensate for the missing
coordinate–momentum uncertainty. For weaker excitonic couplings
as they occur between domains of strongly coupled pigments, we find
the principle of detailed balance to be fulfilled in a semiclassical
variant of the generalized Förster theory. A qualitatively
correct description of the transfer is obtained with this theory at
a significantly lower computational cost as with the quantum generalized
Förster theory. Larger deviations between the two theories
are expected for large energy gaps as they occur in complexes with
chemically different pigments.
Collapse
Affiliation(s)
- Thomas Renger
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria
| |
Collapse
|
22
|
Kimura A, Kitoh-Nishioka H, Shigeta Y, Itoh S. Comparison between the Light-Harvesting Mechanisms of Type-I Photosynthetic Reaction Centers of Heliobacteria and Photosystem I: Pigment Site Energy Distribution and Exciton State. J Phys Chem B 2021; 125:3727-3738. [DOI: 10.1021/acs.jpcb.0c09400] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Affiliation(s)
- Akihiro Kimura
- Department of Physics, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan
| | - Hirotaka Kitoh-Nishioka
- JST, PRESTO, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8571, Japan
- Graduate School of System Informatics, Kobe University, Kobe 657-8501, Japan
| | - Yasuteru Shigeta
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba, Ibaraki 305-8571, Japan
| | - Shigeru Itoh
- Department of Physics, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan
| |
Collapse
|
23
|
Reinot T, Khmelnitskiy A, Kell A, Jassas M, Jankowiak R. Exciton Lifetime Distributions and Population Dynamics in the FMO Protein Complex from Prosthecochloris aestuarii. ACS OMEGA 2021; 6:5990-6008. [PMID: 33681637 PMCID: PMC7931385 DOI: 10.1021/acsomega.1c00286] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/20/2021] [Accepted: 02/09/2021] [Indexed: 06/12/2023]
Abstract
Significant protein rearrangement upon excitation and energy transfer in Fenna-Matthews-Olson protein of Prosthecochloris aestuarii results in a modified energy landscape, which induces more changes in pigment site energies than predicted by the "standard" hole-burning theory. The energy changes are elucidated by simulations while investigating the effects of site-dependent disorder, both static (site-energy distribution widths) and dynamic (spectral density shapes). The resulting optimized site energies and their fluctuations are consistent with relative differences observed in inhomogeneous widths calculated by recent molecular dynamic simulations. Two sets of different spectral densities reveal how their shapes affect the population dynamics and distribution of exciton lifetimes. Calculations revealed the wavelength-dependent distributions of exciton lifetimes (T 1) in the femtosecond to picosecond time frame. We suggest that the calculated multimodal and asymmetric wavelength-dependent T 1 distributions offer more insight into the interpretation of resonant hole-burned (HB) spectra, kinetic traces in two-dimensional (2D) electronic spectroscopy experiments, and widely used global analyses in fitting data from transient absorption experiments.
Collapse
Affiliation(s)
- Tonu Reinot
- Department
of Chemistry, Department of Physics, Kansas State University, Manhattan, Kansas 66506, United States
| | - Anton Khmelnitskiy
- Department
of Chemistry, Department of Physics, Kansas State University, Manhattan, Kansas 66506, United States
| | - Adam Kell
- Department
of Chemistry, Department of Physics, Kansas State University, Manhattan, Kansas 66506, United States
| | - Mahboobe Jassas
- Department
of Chemistry, Department of Physics, Kansas State University, Manhattan, Kansas 66506, United States
| | - Ryszard Jankowiak
- Department
of Chemistry, Department of Physics, Kansas State University, Manhattan, Kansas 66506, United States
| |
Collapse
|
24
|
Sirohiwal A, Neese F, Pantazis DA. 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. Chem Sci 2021; 12:4463-4476. [PMID: 34163712 PMCID: PMC8179452 DOI: 10.1039/d0sc06616h] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Abstract
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.
Collapse
Affiliation(s)
- Abhishek Sirohiwal
- Max-Planck-Institut für Kohlenforschung Kaiser-Wilhelm-Platz 1 45470 Mülheim an der Ruhr Germany.,Fakultät für Chemie und Biochemie, Ruhr-Universität Bochum 44780 Bochum Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung Kaiser-Wilhelm-Platz 1 45470 Mülheim an der Ruhr Germany
| | - Dimitrios A Pantazis
- Max-Planck-Institut für Kohlenforschung Kaiser-Wilhelm-Platz 1 45470 Mülheim an der Ruhr Germany
| |
Collapse
|
25
|
Rose PA, Krich JJ. Efficient numerical method for predicting nonlinear optical spectroscopies of open systems. J Chem Phys 2021; 154:034108. [PMID: 33499622 DOI: 10.1063/5.0024104] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Abstract
Nonlinear optical spectroscopies are powerful tools for probing quantum dynamics in molecular and nanoscale systems. While intuition about ultrafast spectroscopies is often built by considering impulsive optical pulses, actual experiments have finite-duration pulses, which can be important for interpreting and predicting experimental results. We present a new freely available open source method for spectroscopic modeling, called Ultrafast Ultrafast (UF2) spectroscopy, which enables computationally efficient and convenient prediction of nonlinear spectra, such as treatment of arbitrary finite duration pulse shapes. UF2 is a Fourier-based method that requires diagonalization of the Liouvillian propagator of the system density matrix. We also present a Runge-Kutta-Euler (RKE) direct propagation method. We include open system dynamics in the secular Redfield, full Redfield, and Lindblad formalisms with Markovian baths. For non-Markovian systems, the degrees of freedom corresponding to memory effects are brought into the system and treated nonperturbatively. We analyze the computational complexity of the algorithms and demonstrate numerically that, including the cost of diagonalizing the propagator, UF2 is 20-200 times faster than the direct propagation method for secular Redfield models with arbitrary Hilbert space dimension; it is similarly faster for full Redfield models at least up to system dimensions where the propagator requires more than 20 GB to store; and for Lindblad models, it is faster up to Hilbert space dimension near 100 with speedups for small systems by factors of over 500. UF2 and RKE are part of a larger open source Ultrafast Software Suite, which includes tools for automatic generation and calculation of Feynman diagrams.
Collapse
Affiliation(s)
- Peter A Rose
- Department of Physics, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
| | - Jacob J Krich
- Department of Physics, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
| |
Collapse
|
26
|
Chaillet M, Lengauer F, Adolphs J, Müh F, Fokas AS, Cole DJ, Chin AW, Renger T. Static Disorder in Excitation Energies of the Fenna-Matthews-Olson Protein: Structure-Based Theory Meets Experiment. J Phys Chem Lett 2020; 11:10306-10314. [PMID: 33227205 PMCID: PMC7751012 DOI: 10.1021/acs.jpclett.0c03123] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2020] [Accepted: 11/17/2020] [Indexed: 06/11/2023]
Abstract
Inhomogeneous broadening of optical lines of the Fenna-Matthews-Olson (FMO) light-harvesting protein is investigated by combining a Monte Carlo sampling of low-energy conformational substates of the protein with a quantum chemical/electrostatic calculation of local transition energies (site energies) of the pigments. The good agreement between the optical spectra calculated for the inhomogeneous ensemble and the experimental data demonstrates that electrostatics is the dominant contributor to static disorder in site energies. Rotamers of polar amino acid side chains are found to cause bimodal distribution functions of site energy shifts, which can be probed by hole burning and single-molecule spectroscopy. When summing over the large number of contributions, the resulting distribution functions of the site energies become Gaussians, and the correlations in site energy fluctuations at different sites practically average to zero. These results demonstrate that static disorder in the FMO protein is in the realm of the central limit theorem of statistics.
Collapse
Affiliation(s)
- Marten
L. Chaillet
- Bijvoet
Centre for Biomolecular Research, University
of Utrecht, Heidelberglaan 8, 3584 CS Utrecht, The Netherlands
| | - Florian Lengauer
- Institute
of Theoretical Physics, Johannes Kepler
University Linz, Altenberger Str. 69, 4040 Linz, Austria
| | - Julian Adolphs
- Leibniz
Institute for Agricultural Engineering and Bioeconomy, Max-Eyth-Allee 100, 14469 Potsdam, Germany
| | - Frank Müh
- Institute
of Theoretical Physics, Johannes Kepler
University Linz, Altenberger Str. 69, 4040 Linz, Austria
| | - Alexander S. Fokas
- TCM
Group, Cavendish Laboratory, 19 J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Daniel J. Cole
- School
of Natural and Environmental Sciences, Newcastle
University, Newcastle
upon Tyne NE1 7RU, United
Kingdom
| | - Alex W. Chin
- Centre
National de la Recherce Scientifique, Institute des Nanosciences de
Paris, Sorbonne Université, Paris, France
| | - Thomas Renger
- Institute
of Theoretical Physics, Johannes Kepler
University Linz, Altenberger Str. 69, 4040 Linz, Austria
| |
Collapse
|
27
|
Klinger A, Lindorfer D, Müh F, Renger T. Normal mode analysis of spectral density of FMO trimers: Intra- and intermonomer energy transfer. J Chem Phys 2020; 153:215103. [PMID: 33291900 DOI: 10.1063/5.0027994] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
The intermolecular contribution to the spectral density of the exciton-vibrational coupling of the homotrimeric Fenna-Matthews-Olson (FMO) light-harvesting protein of green sulfur bacteria P. aestuarii is analyzed by combining a normal mode analysis of the protein with the charge density coupling method for the calculation of local transition energies of the pigments. Correlations in site energy fluctuations across the whole FMO trimer are found at low vibrational frequencies. Including, additionally, the high-frequency intrapigment part of the spectral density, extracted from line-narrowing spectra, we study intra- and intermonomer exciton transfer. Whereas the intrapigment part of the spectral density is important for fast intramonomer exciton relaxation, the intermolecular contributions (due to pigment-environment coupling) determine the intermonomer exciton transfer. Neither the variations of the local Huang-Rhys factors nor the correlations in site energy fluctuations have a critical influence on energy transfer. At room temperature, the intermonomer transfer in the FMO protein occurs on a 10 ps time scale, whereas intramonomer exciton equilibration is roughly two orders of magnitude faster. At cryogenic temperatures, intermonomer transfer limits the lifetimes of the lowest exciton band. The lifetimes are found to increase between 20 ps in the center of this band up to 100 ps toward lower energies, which is in very good agreement with the estimates from hole burning data. Interestingly, exciton delocalization in the FMO monomers is found to slow down intermonomer energy transfer, at both physiological and cryogenic temperatures.
Collapse
Affiliation(s)
- Alexander Klinger
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria
| | - Dominik Lindorfer
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria
| | - Frank Müh
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria
| | - Thomas Renger
- Institut für Theoretische Physik, Johannes Kepler Universität Linz, Altenberger Str. 69, 4040 Linz, Austria
| |
Collapse
|
28
|
Suzuki Y, Watanabe H, Okiyama Y, Ebina K, Tanaka S. Comparative study on model parameter evaluations for the energy transfer dynamics in Fenna–Matthews–Olson complex. Chem Phys 2020. [DOI: 10.1016/j.chemphys.2020.110903] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
29
|
Kunsel T, Löhner A, Mayo JJ, Köhler J, Jansen TLC, Knoester J. Unraveling intra-aggregate structural disorder using single-molecule spectroscopy. J Chem Phys 2020; 153:134304. [PMID: 33032400 DOI: 10.1063/5.0023551] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Structural disorder within self-assembled molecular aggregates may have strong effects on their optical functionality. Such disorder, however, is hard to explore using standard ensemble measurements. In this paper, we report on the characterization of intra-aggregate structural disorder through a linewidth analysis of fluorescence excitation experiments on individual zinc-chlorin (ZnChl) nanotubular molecular aggregates. Recent experiments suggest an anomaly in the linewidths of the two absorption bands that dominate the spectra: the higher-energy bands on average show a smaller linewidth than the lower-energy bands. This anomaly is explored in this paper by analyzing and modeling the correlation of the two linewidths for each aggregate. We exploit a Frenkel exciton model to show that the experimentally observed correlation of linewidths and other statistical properties of the single-aggregate spectra can be explained from small variations of the molecular orientations within individual aggregates.
Collapse
Affiliation(s)
- T Kunsel
- University of Groningen, Zernike Institute for Advanced Materials, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - A Löhner
- Spectroscopy of Soft Matter, University of Bayreuth, Universitätsstraße 30, 94557 Bayreuth, Germany
| | - J J Mayo
- University of Groningen, Zernike Institute for Advanced Materials, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - J Köhler
- Spectroscopy of Soft Matter, University of Bayreuth, Universitätsstraße 30, 94557 Bayreuth, Germany
| | - T L C Jansen
- University of Groningen, Zernike Institute for Advanced Materials, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - J Knoester
- University of Groningen, Zernike Institute for Advanced Materials, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| |
Collapse
|
30
|
Analysis of Photosynthetic Systems and Their Applications with Mathematical and Computational Models. APPLIED SCIENCES-BASEL 2020. [DOI: 10.3390/app10196821] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
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.
Collapse
|
31
|
Chen MS, Zuehlsdorff TJ, Morawietz T, Isborn CM, Markland TE. Exploiting Machine Learning to Efficiently Predict Multidimensional Optical Spectra in Complex Environments. J Phys Chem Lett 2020; 11:7559-7568. [PMID: 32808797 DOI: 10.1021/acs.jpclett.0c02168] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
The excited-state dynamics of chromophores in complex environments determine a range of vital biological and energy capture processes. Time-resolved, multidimensional optical spectroscopies provide a key tool to investigate these processes. Although theory has the potential to decode these spectra in terms of the electronic and atomistic dynamics, the need for large numbers of excited-state electronic structure calculations severely limits first-principles predictions of multidimensional optical spectra for chromophores in the condensed phase. Here, we leverage the locality of chromophore excitations to develop machine learning models to predict the excited-state energy gap of chromophores in complex environments for efficiently constructing linear and multidimensional optical spectra. By analyzing the performance of these models, which span a hierarchy of physical approximations, across a range of chromophore-environment interaction strengths, we provide strategies for the construction of machine learning models that greatly accelerate the calculation of multidimensional optical spectra from first principles.
Collapse
Affiliation(s)
- Michael S Chen
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Tim J Zuehlsdorff
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Tobias Morawietz
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Christine M Isborn
- Chemistry and Chemical Biology, University of California Merced, Merced, California 95343, United States
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| |
Collapse
|
32
|
Yan Y, Liu Y, Xing T, Shi Q. Theoretical study of excitation energy transfer and nonlinear spectroscopy of photosynthetic light‐harvesting complexes using the nonperturbative reduced dynamics method. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2020. [DOI: 10.1002/wcms.1498] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Yaming Yan
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species Institute of Chemistry, Chinese Academy of Sciences Beijing China
- University of Chinese Academy of Sciences Beijing China
- Physical Science Laboratory Huairou National Comprehensive Science Center Beijing China
| | - Yanying Liu
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species Institute of Chemistry, Chinese Academy of Sciences Beijing China
- University of Chinese Academy of Sciences Beijing China
- Physical Science Laboratory Huairou National Comprehensive Science Center Beijing China
| | - Tao Xing
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species Institute of Chemistry, Chinese Academy of Sciences Beijing China
- University of Chinese Academy of Sciences Beijing China
- Physical Science Laboratory Huairou National Comprehensive Science Center Beijing China
| | - Qiang Shi
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species Institute of Chemistry, Chinese Academy of Sciences Beijing China
- University of Chinese Academy of Sciences Beijing China
- Physical Science Laboratory Huairou National Comprehensive Science Center Beijing China
| |
Collapse
|
33
|
Sánchez Muñoz C, Schlawin F. Photon Correlation Spectroscopy as a Witness for Quantum Coherence. PHYSICAL REVIEW LETTERS 2020; 124:203601. [PMID: 32501097 DOI: 10.1103/physrevlett.124.203601] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/10/2020] [Accepted: 04/28/2020] [Indexed: 06/11/2023]
Abstract
The development of spectroscopic techniques able to detect and verify quantum coherence is a goal of increasing importance given the rapid progress of new quantum technologies, the advances in the field of quantum thermodynamics, and the emergence of new questions in chemistry and biology regarding the possible relevance of quantum coherence in biochemical processes. Ideally, these tools should be able to detect and verify the presence of quantum coherence in both the transient dynamics and the steady state of driven-dissipative systems, such as light-harvesting complexes driven by thermal photons in natural conditions. This requirement poses a challenge for standard laser spectroscopy methods. Here, we propose photon correlation measurements as a new tool to analyze quantum dynamics in molecular aggregates in driven-dissipative situations. We show that the photon correlation statistics of the light emitted in several models of molecular aggregates can signal the presence of coherent dynamics. Deviations from the counting statistics of independent emitters constitute a direct fingerprint of quantum coherence in the steady state. Furthermore, the analysis of frequency resolved photon correlations can signal the presence of coherent dynamics even in the absence of steady state coherence, providing direct spectroscopic access to the much sought-after site energies in molecular aggregates.
Collapse
Affiliation(s)
- Carlos Sánchez Muñoz
- Clarendon Laboratory, University of Oxford, Parks Road, Oxford OX1 3PU, United Kingdom
| | - Frank Schlawin
- Clarendon Laboratory, University of Oxford, Parks Road, Oxford OX1 3PU, United Kingdom
| |
Collapse
|
34
|
Kramer T, Rodríguez M. Effect of disorder and polarization sequences on two-dimensional spectra of light-harvesting complexes. PHOTOSYNTHESIS RESEARCH 2020; 144:147-154. [PMID: 31872335 DOI: 10.1007/s11120-019-00699-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/23/2019] [Accepted: 12/10/2019] [Indexed: 06/10/2023]
Abstract
Two-dimensional electronic spectra (2DES) provide unique ways to track the energy transfer dynamics in light-harvesting complexes. The interpretation of the peaks and structures found in experimentally recorded 2DES is often not straightforward, since several processes are imaged simultaneously. The choice of specific pulse polarization sequences helps to disentangle the sometimes convoluted spectra, but brings along other disturbances. We show by detailed theoretical calculations how 2DES of the Fenna-Matthews-Olson complex are affected by rotational and conformational disorder of the chromophores.
Collapse
|
35
|
Tong Z, Huai Z, Mei Y, Mo Y. Reproducing the low-temperature excitation energy transfer dynamics of phycoerythrin 545 light-harvesting complex with a structure-based model Hamiltonian. J Chem Phys 2020; 152:135101. [DOI: 10.1063/1.5135999] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Zhengqing Tong
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| | - Zhe Huai
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| | - Yan Mo
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| |
Collapse
|
36
|
Cao J, Cogdell RJ, Coker DF, Duan HG, Hauer J, Kleinekathöfer U, Jansen TLC, Mančal T, Miller RJD, Ogilvie JP, Prokhorenko VI, Renger T, Tan HS, Tempelaar R, Thorwart M, Thyrhaug E, Westenhoff S, Zigmantas D. Quantum biology revisited. SCIENCE ADVANCES 2020; 6:eaaz4888. [PMID: 32284982 PMCID: PMC7124948 DOI: 10.1126/sciadv.aaz4888] [Citation(s) in RCA: 175] [Impact Index Per Article: 43.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/12/2019] [Accepted: 01/06/2020] [Indexed: 05/21/2023]
Abstract
Photosynthesis is a highly optimized process from which valuable lessons can be learned about the operating principles in nature. Its primary steps involve energy transport operating near theoretical quantum limits in efficiency. Recently, extensive research was motivated by the hypothesis that nature used quantum coherences to direct energy transfer. This body of work, a cornerstone for the field of quantum biology, rests on the interpretation of small-amplitude oscillations in two-dimensional electronic spectra of photosynthetic complexes. This Review discusses recent work reexamining these claims and demonstrates that interexciton coherences are too short lived to have any functional significance in photosynthetic energy transfer. Instead, the observed long-lived coherences originate from impulsively excited vibrations, generally observed in femtosecond spectroscopy. These efforts, collectively, lead to a more detailed understanding of the quantum aspects of dissipation. Nature, rather than trying to avoid dissipation, exploits it via engineering of exciton-bath interaction to create efficient energy flow.
Collapse
Affiliation(s)
- Jianshu Cao
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA
| | - Richard J. Cogdell
- Institute of Molecular, Cell and Systems Biology, College of Medical, Veterinary and Life Science, University of Glasgow, Glasgow G12 8QQ, UK
| | - David F. Coker
- Department of Chemistry, Boston University, 590 Commonwealth Avenue, Boston, MA 02215, USA
| | - Hong-Guang Duan
- Atomically Resolved Dynamics Department, Max Planck Institute for the Structure and Dynamics of Matter, 22761 Hamburg, Germany
- I. Institut für Theoretische Physik, Universität Hamburg, Jungiusstrasse 9, 20355 Hamburg, Germany
- The Hamburg Centre for Ultrafast Imaging, Universität Hamburg, 22761 Hamburg, Germany
| | - Jürgen Hauer
- Technische Universität München, Dynamische Spektroskopien, Fakultät für Chemie, Lichtenbergstr. 4, 85748 Garching, Germany, and Photonics Institute, TU Wien, 1040 Vienna, Austria
| | - Ulrich Kleinekathöfer
- Department of Physics and Earth Science, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Thomas L. C. Jansen
- Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, Netherlands
| | - Tomáš Mančal
- Faculty of Mathematics and Physics, Charles University, Ke Karlovu 5, CZ-12116 Prague 2, Czech Republic
| | - R. J. Dwayne Miller
- Atomically Resolved Dynamics Department, Max Planck Institute for the Structure and Dynamics of Matter, 22761 Hamburg, Germany
- The Hamburg Centre for Ultrafast Imaging, Universität Hamburg, 22761 Hamburg, Germany
- Departments of Chemistry and Physics, University of Toronto, Toronto, ON M5S 3H6, Canada
| | | | - Valentyn I. Prokhorenko
- Atomically Resolved Dynamics Department, Max Planck Institute for the Structure and Dynamics of Matter, 22761 Hamburg, Germany
| | - Thomas Renger
- Institute of Theoretical Physics, Department of Theoretical Biophysics, Johannes Kepler University Linz, Altenberger Str. 69, 4040 Linz, Austria
| | - Howe-Siang Tan
- Division of Chemistry and Biological Chemistry, School of Physical and Mathematical Sciences, Nanyang Technological University, 21 Nanyang Link, Singapore 637371, Singapore
| | - Roel Tempelaar
- Department of Chemistry, Columbia University, 3000 Broadway, New York, NY 10027, USA
| | - Michael Thorwart
- I. Institut für Theoretische Physik, Universität Hamburg, Jungiusstrasse 9, 20355 Hamburg, Germany
- The Hamburg Centre for Ultrafast Imaging, Universität Hamburg, 22761 Hamburg, Germany
| | - Erling Thyrhaug
- Technische Universität München, Dynamische Spektroskopien, Fakultät für Chemie, Lichtenbergstr. 4, 85748 Garching, Germany, and Photonics Institute, TU Wien, 1040 Vienna, Austria
| | - Sebastian Westenhoff
- Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg 40530, Sweden
| | | |
Collapse
|
37
|
Kim Y, Morozov D, Stadnytskyi V, Savikhin S, Slipchenko LV. Predictive First-Principles Modeling of a Photosynthetic Antenna Protein: The Fenna-Matthews-Olson Complex. J Phys Chem Lett 2020; 11:1636-1643. [PMID: 32013435 DOI: 10.1021/acs.jpclett.9b03486] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
High efficiency of light harvesting in photosynthetic pigment-protein complexes is governed by evolutionary-perfected protein-assisted tuning of individual pigment properties and interpigment interactions. Due to the large number of spectrally overlapping pigments in a typical photosynthetic complex, experimental methods often fail to unambiguously identify individual chromophore properties. Here, we report a first-principles-based modeling protocol capable of predicting properties of pigments in protein environment to a high precision. The technique was applied to successfully uncover electronic properties of the Fenna-Matthews-Olson (FMO) pigment-protein complex. Each of the three subunits of the FMO complex contains eight strongly coupled bacteriochlorophyll a (BChl a) pigments. The excitonic structure of FMO can be described by an electronic Hamiltonian containing excitation (site) energies of BChl a pigments and electronic couplings between them. Several such Hamiltonians have been developed in the past based on the information from various spectroscopic measurements of FMO; however, fine details of the excitonic structure and energy transfer in FMO, especially assignments of short-lived high-energy sites, remain elusive. Utilizing polarizable embedding quantum mechanics/molecular mechanics with the effective fragment potentials, we computed the electronic Hamiltonian of FMO that is in general agreement with previously reported empirical Hamiltonians and quantitatively reproduces experimental absorption and circular dichroism spectra of the FMO protein. The developed computational protocol is sufficiently simple and can be utilized for predictive modeling of other wild-type and mutated photosynthetic pigment-protein complexes.
Collapse
Affiliation(s)
- Yongbin Kim
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| | - Dmitry Morozov
- Nanoscience Center and Department of Chemistry, University of Jyväskylä, P.O. Box 35, Jyväskylä 40014, Finland
| | - Valentyn Stadnytskyi
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
- Laboratory of Chemical Physics, National Institute of Diabetes, Digestion and Kidney Diseases, National Institutes of Health, 5 Memorial Drive, Bethesda, Maryland 20892, United States
| | - Sergei Savikhin
- Department of Physics and Astronomy, Purdue University, 525 Northwestern Avenue, West Lafayette, Indiana 47907, United States
| | - Lyudmila V Slipchenko
- Department of Chemistry, Purdue University, 560 Oval Drive, West Lafayette, Indiana 47907, United States
| |
Collapse
|
38
|
Kell A, Khmelnitskiy AY, Reinot T, Jankowiak R. On uncorrelated inter-monomer Förster energy transfer in Fenna-Matthews-Olson complexes. J R Soc Interface 2020; 16:20180882. [PMID: 30958204 DOI: 10.1098/rsif.2018.0882] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
The Fenna-Matthews-Olson (FMO) light-harvesting antenna protein of green sulfur bacteria is a long-studied pigment-protein complex which funnels energy from the chlorosome to the reaction centre where photochemistry takes place. The structure of the FMO protein from Chlorobaculum tepidum is known as a homotrimeric complex containing eight bacteriochlorophyll a per monomer. Owing to this structure FMO has strong intra-monomer and weak inter-monomer electronic coupling constants. While long-lived (sub-picosecond) coherences within a monomer have been a prevalent topic of study over the past decade, various experimental evidence supports the presence of subsequent inter-monomer energy transfer on a picosecond time scale. The latter has been neglected by most authors in recent years by considering only sub-picosecond time scales or assuming that the inter-monomer coupling between low-energy states is too weak to warrant consideration of the entire trimer. However, Förster theory predicts that energy transfer of the order of picoseconds is possible even for very weak (less than 5 cm-1) electronic coupling between chromophores. This work reviews experimental data (with a focus on emission and hole-burned spectra) and simulations of exciton dynamics which demonstrate inter-monomer energy transfer. It is shown that the lowest energy 825 nm absorbance band cannot be properly described by a single excitonic state. The energy transfer through FMO is modelled by generalized Förster theory using a non-Markovian, reduced density matrix approach to describe the electronic structure. The disorder-averaged inter-monomer transfer time across the 825 nm band is about 27 ps. While only isolated FMO proteins are presented, the presence of inter-monomer energy transfer in the context of the overall photosystem is also briefly discussed.
Collapse
Affiliation(s)
- Adam Kell
- 1 Department of Chemistry, Kansas State University , Manhattan, KS , USA
| | | | - Tonu Reinot
- 1 Department of Chemistry, Kansas State University , Manhattan, KS , USA
| | - Ryszard Jankowiak
- 1 Department of Chemistry, Kansas State University , Manhattan, KS , USA.,2 Department of Physics, Kansas State University , Manhattan, KS , USA
| |
Collapse
|
39
|
Müh F, Zouni A. Structural basis of light-harvesting in the photosystem II core complex. Protein Sci 2020; 29:1090-1119. [PMID: 32067287 PMCID: PMC7184784 DOI: 10.1002/pro.3841] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/26/2020] [Revised: 02/06/2020] [Accepted: 02/06/2020] [Indexed: 12/20/2022]
Abstract
Photosystem II (PSII) is a membrane-spanning, multi-subunit pigment-protein complex responsible for the oxidation of water and the reduction of plastoquinone in oxygenic photosynthesis. In the present review, the recent explosive increase in available structural information about the PSII core complex based on X-ray crystallography and cryo-electron microscopy is described at a level of detail that is suitable for a future structure-based analysis of light-harvesting processes. This description includes a proposal for a consistent numbering scheme of protein-bound pigment cofactors across species. The structural survey is complemented by an overview of the state of affairs in structure-based modeling of excitation energy transfer in the PSII core complex with emphasis on electrostatic computations, optical properties of the reaction center, the assignment of long-wavelength chlorophylls, and energy trapping mechanisms.
Collapse
Affiliation(s)
- Frank Müh
- Department of Theoretical Biophysics, Institute for Theoretical Physics, Johannes Kepler University Linz, Linz, Austria
| | - Athina Zouni
- Humboldt-Universität zu Berlin, Institute for Biology, Biophysics of Photosynthesis, Berlin, Germany
| |
Collapse
|
40
|
Kaliakin DS, Nakata H, Kim Y, Chen Q, Fedorov DG, Slipchenko LV. FMOxFMO: Elucidating Excitonic Interactions in the Fenna-Matthews-Olson Complex with the Fragment Molecular Orbital Method. J Chem Theory Comput 2020; 16:1175-1187. [PMID: 31841349 DOI: 10.1021/acs.jctc.9b00621] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In order to study Förster resonance energy transfer (FRET), the fragment molecular orbital (FMO) method is extended to compute electronic couplings between local excitations via the excited state transition density model, enabling efficient calculations of nonlocal excitations in a large molecular system and overcoming the previous limitation of being able to compute only local excitations. The results of these simple but accurate models are validated against full quantum calculations without fragmentation. The developed method is applied to a very important photosynthetic pigment-protein complex, the Fenna-Matthews-Olson complex (FMOc), that is responsible for the energy transfer from a chlorosome to the reaction center in the green sulfur bacteria. Absorption and circular dichroism spectra of FMOc are simulated, and the role of the molecular environment on the excitations is revealed.
Collapse
Affiliation(s)
- Danil S Kaliakin
- Department of Chemistry , Purdue University , 560 Oval Drive , West Lafayette , Indiana 47907 , United States
| | - Hiroya Nakata
- Research Institute for Advanced Materials and Devices , Kyocera , 5-3 Hikaridai-3 , Seika-cho Soraku-gun, Kyoto 619-0237 , Japan
| | - Yongbin Kim
- Department of Chemistry , Purdue University , 560 Oval Drive , West Lafayette , Indiana 47907 , United States
| | - Qifeng Chen
- Department of Chemistry , Purdue University , 560 Oval Drive , West Lafayette , Indiana 47907 , United States
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat) , National Institute of Advanced Industrial Science and Technology (AIST) , Central 2, Umezono 1-1-1 , Tsukuba 305-8568 , Japan
| | - Lyudmila V Slipchenko
- Department of Chemistry , Purdue University , 560 Oval Drive , West Lafayette , Indiana 47907 , United States
| |
Collapse
|
41
|
Saller MAC, Kelly A, Richardson JO. Improved population operators for multi-state nonadiabatic dynamics with the mixed quantum-classical mapping approach. Faraday Discuss 2020; 221:150-167. [DOI: 10.1039/c9fd00050j] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Application to the 7-state Frenkel-exciton Hamiltonian for the Fenna–Matthews–Olson complex shows that using a different representation of the electronic population operators can drastically improve the accuracy of the quasiclassical mapping approach without increasing the computational effort.
Collapse
Affiliation(s)
| | - Aaron Kelly
- Department of Chemistry
- Dalhousie University
- Halifax
- Canada
| | | |
Collapse
|
42
|
Paulikat M, Mata RA, Gelabert R. A high-throughput computational approach to UV-Vis spectra in protein mutants. Phys Chem Chem Phys 2019; 21:20678-20692. [PMID: 31508628 DOI: 10.1039/c9cp03908b] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
In this work we present a high-throughput approach to the computation of absorption UV-Vis spectra tailored to mutagenesis studies. The scheme makes use of a single molecular dynamics trajectory of a reference (non-mutated) species. The shifts in absorption energy caused by a residue mutation are evaluated by building an effective potential of the environment and computing a correction term based on perturbation theory. The sampling is only performed in the phase space of the initial protein. We analyze the robustness of the method by comparing different approximations for the effective potential, the sampling of mutant residue geometries and observing the impact in the prediction of both bathocromic and hypsochromic shifts. As a test subject, we consider a red fluorescent protein variant with potential biotechnological applications.
Collapse
Affiliation(s)
- Mirko Paulikat
- Institute of Physical Chemistry, University of Goettingen, Tammannstraße 6, D-37077 Göttingen, Germany.
| | - Ricardo A Mata
- Institute of Physical Chemistry, University of Goettingen, Tammannstraße 6, D-37077 Göttingen, Germany.
| | - Ricard Gelabert
- Departament de Química, Universitat Autònoma de Barcelona, 08193 Bellaterra, Barcelona, Spain.
| |
Collapse
|
43
|
Chin A, Mangaud E, Chevet V, Atabek O, Desouter-Lecomte M. Visualising the role of non-perturbative environment dynamics in the dissipative generation of coherent electronic motion. Chem Phys 2019. [DOI: 10.1016/j.chemphys.2019.110392] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
|
44
|
Oh SA, Coker DF, Hutchinson DAW. Optimization of energy transport in the Fenna-Matthews-Olson complex via site-varying pigment-protein interactions. J Chem Phys 2019; 150:085102. [DOI: 10.1063/1.5048058] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023] Open
Affiliation(s)
- S. A. Oh
- Dodd-Walls Centre for Photonic and Quantum Technologies, Department of Physics, University of Otago, Dunedin 9016, New Zealand
| | - D. F. Coker
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | - D. A. W. Hutchinson
- Dodd-Walls Centre for Photonic and Quantum Technologies, Department of Physics, University of Otago, Dunedin 9016, New Zealand
- Centre for Quantum Technologies, National University of Singapore, 3 Science Drive 2, Singapore 117543, Singapore
| |
Collapse
|
45
|
Lee SH, Matula AJ, Hu G, Troiano JL, Karpovich CJ, Crabtree RH, Batista VS, Brudvig GW. Strongly Coupled Phenazine-Porphyrin Dyads: Light-Harvesting Molecular Assemblies with Broad Absorption Coverage. ACS APPLIED MATERIALS & INTERFACES 2019; 11:8000-8008. [PMID: 30698407 DOI: 10.1021/acsami.8b20996] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
The development of light-harvesting architectures with broad absorption coverage in the visible region continues to be an important research area in the field of artificial photosynthesis. Here, we introduce a new class of ethynyl-linked panchromatic dyads composed of dibenzophenazines coupled ortho and meta to tetrapyrroles with an anchoring group that can be grafted onto metal oxide surfaces. Quantum chemical calculations and photophysical measurements of the synthesized materials reveal that both of the dibenzophenazine dyads absorb broadly from 300 to 636 nm and exhibit absorption bands different from those of the constituent chromophore units. Moreover, the different points of attachment of dibenzophenazines to tetrapyrroles give different absorption profiles which computations suggest result from differences in the planarity of the two dyads. Applicability of the dyads in artificial photosynthesis systems was assessed by their incorporation and characterization of their performance in dye-sensitized solar cells.
Collapse
Affiliation(s)
- Shin Hee Lee
- Department of Chemistry, and Yale Energy Sciences Institute , Yale University , New Haven , Connecticut 06520-8107 , United States
| | - Adam J Matula
- Department of Chemistry, and Yale Energy Sciences Institute , Yale University , New Haven , Connecticut 06520-8107 , United States
| | - Gongfang Hu
- Department of Chemistry, and Yale Energy Sciences Institute , Yale University , New Haven , Connecticut 06520-8107 , United States
| | - Jennifer L Troiano
- Department of Chemistry, and Yale Energy Sciences Institute , Yale University , New Haven , Connecticut 06520-8107 , United States
| | - Christopher J Karpovich
- Department of Chemistry, and Yale Energy Sciences Institute , Yale University , New Haven , Connecticut 06520-8107 , United States
| | - Robert H Crabtree
- Department of Chemistry, and Yale Energy Sciences Institute , Yale University , New Haven , Connecticut 06520-8107 , United States
| | - Victor S Batista
- Department of Chemistry, and Yale Energy Sciences Institute , Yale University , New Haven , Connecticut 06520-8107 , United States
| | - Gary W Brudvig
- Department of Chemistry, and Yale Energy Sciences Institute , Yale University , New Haven , Connecticut 06520-8107 , United States
| |
Collapse
|
46
|
Tong Z, Huai Z, Mei Y, Mo Y. Influence of the Protein Environment on the Electronic Excitation of Chromophores in the Phycoerythrin 545 Light–Harvesting Complex: A Combined MD-QM/MM Method with Polarized Protein–Specific Charge Scheme. J Phys Chem B 2019; 123:2040-2049. [DOI: 10.1021/acs.jpcb.8b11764] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Zhengqing Tong
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Zhe Huai
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| | - Yan Mo
- State Key Laboratory of Precision Spectroscopy, School of Physics and Materials Science, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| |
Collapse
|
47
|
Hsieh ST, Zhang L, Ye DW, Huang X, Cheng YC. A theoretical study on the dynamics of light harvesting in the dimeric photosystem II core complex: regulation and robustness of energy transfer pathways. Faraday Discuss 2019; 216:94-115. [PMID: 31016302 DOI: 10.1039/c8fd00205c] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023]
Abstract
Here we present our theoretical investigations into the light reaction in the dimeric photosystem II (PSII) core complex. An effective model for excitation energy transfer (EET) and primary charge separation (CS) in the PSII core complex was developed, with model parameters constructed based on molecular dynamics (MD) simulation data. Compared to experimental results, we demonstrated that this model faithfully reproduces the absorption spectra of the RC and core light-harvesting complexes (CP43 and CP47) as well as the full EET dynamics among the chromophores in the PSII core complex. We then applied master equation simulations and network analysis to investigate detailed EET plus CS dynamics in the system, allowing us to identify key EET pathways and produce a coarse-grained cluster model for the light reaction in the dimeric PSII core complex. We show that non-equilibrium energy transfer channels play important roles in the efficient light harvesting process and that multiple EET pathways exist between subunits of PSII to ensure the robustness of light harvesting in the system. Furthermore, we revealed that inter-monomer energy transfer dominated by the coupling between the two CLA625 molecules enables efficient energy exchange between two CP47s in the dimeric PSII core complex, which leads to significant energy pooling in the CP47 domain during the light reaction. Our study provides a blueprint for the design of light harvesting in the PSII core and show that a structure-based approach using molecular dynamics simulations and quantum chemistry calculations can be effectively utilized to elucidate the dynamics of light harvesting in complex photosynthetic systems.
Collapse
Affiliation(s)
- Shou-Ting Hsieh
- Department of Chemistry, National Taiwan University, Taipei City, Taiwan.
| | - Lu Zhang
- State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences Fuzhou, Fujian CN 350002, China
| | - De-Wei Ye
- Department of Chemistry, National Taiwan University, Taipei City, Taiwan.
| | - Xuhui Huang
- Department of Chemistry, Institute for Advance Study and School of Science, Hong Kong University of Science and Technology, Kowloon, Hong Kong.
| | - Yuan-Chung Cheng
- Department of Chemistry, National Taiwan University, Taipei City, Taiwan.
| |
Collapse
|
48
|
Kramer T, Noack M, Reinefeld A, Rodríguez M, Zelinskyy Y. Efficient calculation of open quantum system dynamics and time-resolved spectroscopy with distributed memory HEOM (DM-HEOM). J Comput Chem 2018; 39:1779-1794. [DOI: 10.1002/jcc.25354] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2018] [Revised: 04/20/2018] [Accepted: 04/22/2018] [Indexed: 12/17/2022]
Affiliation(s)
- Tobias Kramer
- Zuse Institute Berlin (ZIB), Takustr. 7; 14195 Berlin Germany
- Department of Physics; Harvard University, 17 Oxford Street; Cambridge Massachusetts 02138
| | - Matthias Noack
- Zuse Institute Berlin (ZIB), Takustr. 7; 14195 Berlin Germany
| | | | - Mirta Rodríguez
- Zuse Institute Berlin (ZIB), Takustr. 7; 14195 Berlin Germany
| | | |
Collapse
|
49
|
Collette F, Renger T, Müh F, Schmidt am Busch M. Red/Green Color Tuning of Visual Rhodopsins: Electrostatic Theory Provides a Quantitative Explanation. J Phys Chem B 2018; 122:4828-4837. [DOI: 10.1021/acs.jpcb.8b02702] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Florimond Collette
- Institut für Theoretische Physik, Johannes Kepler University Linz, Altenberger Strasse 69, 4040 Linz, Austria
| | - Thomas Renger
- Institut für Theoretische Physik, Johannes Kepler University Linz, Altenberger Strasse 69, 4040 Linz, Austria
| | - Frank Müh
- Institut für Theoretische Physik, Johannes Kepler University Linz, Altenberger Strasse 69, 4040 Linz, Austria
| | - Marcel Schmidt am Busch
- Institut für Theoretische Physik, Johannes Kepler University Linz, Altenberger Strasse 69, 4040 Linz, Austria
| |
Collapse
|
50
|
Montemayor D, Rivera E, Jang SJ. Computational Modeling of Exciton-Bath Hamiltonians for Light Harvesting 2 and Light Harvesting 3 Complexes of Purple Photosynthetic Bacteria at Room Temperature. J Phys Chem B 2018. [PMID: 29533664 DOI: 10.1021/acs.jpcb.8b00358] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Light harvesting 2 (LH2) complex is the primary component of the photosynthetic unit of purple bacteria that is responsible for harvesting and relaying excitons. The electronic absorption line shape of LH2 contains two major bands at 800 and 850 nm wavelength regions. Under low light conditions, some species of purple bacteria replace LH2 with light harvesting 3 (LH3), a variant form with almost the same structure as the former but with distinctively different spectral features. The major difference between the absorption line shapes of LH2 and LH3 is the shift of the 850 nm band of the former to a new 820 nm region. The microscopic origin of this difference has been the subject of some theoretical/computational investigations. However, the genuine molecular level source of such a difference is not clearly understood yet. This work reports a comprehensive computational study of LH2 and LH3 complexes so as to clarify different molecular level features of LH2 and LH3 complexes and to construct simple exciton-bath models with a common form. All-atomistic molecular dynamics simulations of both LH2 and LH3 complexes provide detailed molecular level structural differences of bacteriochlorophylls (BChls) in the two complexes, in particular, in their patterns of hydrogen bonding (HB) and torsional angles of the acetyl group. Time-dependent density functional theory calculation of the excitation energies of BChls for structures sampled from the MD simulations suggests that the observed differences in the HB and torsional angles cannot fully account for the experimentally observed spectral shift of LH3. Potential sources that can explain the actual spectral shift of LH3 are discussed, and their magnitudes are assessed through fitting of experimental line shapes. These results demonstrate the feasibility of developing simple exciton-bath models for both LH2 and LH3, which can be employed for large-scale exciton quantum dynamics in their aggregates.
Collapse
Affiliation(s)
- Daniel Montemayor
- Department of Chemistry and Biochemistry, Queens College , City University of New York , 65-30 Kissena Boulevard , Queens , New York 11367 , United States.,PhD Programs in Chemistry and Physics, and Initiative for the Theoretical Sciences, Graduate Center , City University of New York , 365 Fifth Avenue , New York , New York 10016 , United States
| | - Eva Rivera
- Department of Chemistry and Biochemistry, Queens College , City University of New York , 65-30 Kissena Boulevard , Queens , New York 11367 , United States.,PhD Programs in Chemistry and Physics, and Initiative for the Theoretical Sciences, Graduate Center , City University of New York , 365 Fifth Avenue , New York , New York 10016 , United States
| | - Seogjoo J Jang
- Department of Chemistry and Biochemistry, Queens College , City University of New York , 65-30 Kissena Boulevard , Queens , New York 11367 , United States.,PhD Programs in Chemistry and Physics, and Initiative for the Theoretical Sciences, Graduate Center , City University of New York , 365 Fifth Avenue , New York , New York 10016 , United States
| |
Collapse
|