1
|
Maity S, Daskalakis V, Jansen TLC, Kleinekathöfer U. Electric Field Susceptibility of Chlorophyll c Leads to Unexpected Excitation Dynamics in the Major Light-Harvesting Complex of Diatoms. J Phys Chem Lett 2024; 15:2499-2510. [PMID: 38410961 PMCID: PMC10926154 DOI: 10.1021/acs.jpclett.3c03241] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2023] [Revised: 02/04/2024] [Accepted: 02/12/2024] [Indexed: 02/28/2024]
Abstract
Diatoms are one of the most abundant photosynthetic organisms on earth and contribute largely to atmospheric oxygen production. They contain fucoxanthin and chlorophyll-a/c binding proteins (FCPs) as light-harvesting complexes with a remarkable adaptation to the fluctuating light on ocean surfaces. To understand the basis of the photosynthetic process in diatoms, the excitation energy funneling within FCPs must be probed. A state-of-the-art multiscale analysis within a quantum mechanics/molecular mechanics framework has been employed. To this end, the chlorophyll (Chl) excitation energies within the FCP complex from the diatom Phaeodactylum tricornutum have been determined. The Chl-c excitation energies were found to be 5-fold more susceptible to electric fields than those of Chl-a pigments and thus are significantly lower in FCP than in organic solvents. This finding challenges the general belief that the excitation energy of Chl-c is always higher than that of Chl-a in FCP proteins and reveals that Chl-c molecules are much more sensitive to electric fields within protein scaffolds than in Chl-a pigments. The analysis of the linear absorption spectrum and the two-dimensional electronic spectra of the FCP complex strongly supports these findings and allows us to study the excitation transfer within the FCP complex.
Collapse
Affiliation(s)
- Sayan Maity
- School
of Science, Constructor University, Campus Ring 1, 28759 Bremen, Germany
| | - Vangelis Daskalakis
- Department
of Chemical Engineering, School of Engineering,
University of Patras, Patras 26504, Greece
| | - Thomas L. C. Jansen
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh 4, 9747 AG Groningen, Netherlands
| | | |
Collapse
|
2
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
3
|
Kubař T, Elstner M, Cui Q. Hybrid Quantum Mechanical/Molecular Mechanical Methods For Studying Energy Transduction in Biomolecular Machines. Annu Rev Biophys 2023; 52:525-551. [PMID: 36791746 PMCID: PMC10810093 DOI: 10.1146/annurev-biophys-111622-091140] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023]
Abstract
Hybrid quantum mechanical/molecular mechanical (QM/MM) methods have become indispensable tools for the study of biomolecules. In this article, we briefly review the basic methodological details of QM/MM approaches and discuss their applications to various energy transduction problems in biomolecular machines, such as long-range proton transports, fast electron transfers, and mechanochemical coupling. We highlight the particular importance for these applications of balancing computational efficiency and accuracy. Using several recent examples, we illustrate the value and limitations of QM/MM methodologies for both ground and excited states, as well as strategies for calibrating them in specific applications. We conclude with brief comments on several areas that can benefit from further efforts to make QM/MM analyses more quantitative and applicable to increasingly complex biological problems.
Collapse
Affiliation(s)
- T Kubař
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany;
| | - M Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Karlsruhe, Germany;
- Institute of Biological Interfaces (IBG-2), Karlsruhe Institute of Technology, Karlsruhe, Germany;
| | - Q Cui
- Department of Chemistry, Boston University, Boston, Massachusetts, USA;
- Department of Physics, Boston University, Boston, Massachusetts, USA
- Department of Biomedical Engineering, Boston University, Boston, Massachusetts, USA
| |
Collapse
|
4
|
Kaiser A, Daoud RE, Aquilante F, Kühn O, De Vico L, Bokarev SI. A Multiconfigurational Wave Function Implementation of the Frenkel Exciton Model for Molecular Aggregates. J Chem Theory Comput 2023; 19:2918-2928. [PMID: 37115036 DOI: 10.1021/acs.jctc.3c00185] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/29/2023]
Abstract
We present an implementation of the Frenkel exciton model into the OpenMolcas program package enabling calculations of collective electronic excited states of molecular aggregates based on a multiconfigurational wave function description of the individual monomers. The computational protocol avoids using diabatization schemes and, thus, supermolecule calculations. Additionally, the use of the Cholesky decomposition of the two-electron integrals entering pair interactions enhances the efficiency of the computational scheme. The application of the method is exemplified for two test systems, that is, a formaldehyde oxime and a bacteriochlorophyll-like dimer. For the sake of comparison with the dipole approximation, we restrict our considerations to situations where intermonomer exchange can be neglected. The protocol is expected to be beneficial for aggregates composed of molecules with extended π systems, unpaired electrons such as radicals or transition metal centers, where it should outperform widely used methods based on time-dependent density functional theory.
Collapse
Affiliation(s)
- Andy Kaiser
- Institut für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock. Germany
| | - Razan E Daoud
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Francesco Aquilante
- Theory and Simulation of Materials (THEOS) and National Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Oliver Kühn
- Institut für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock. Germany
| | - Luca De Vico
- Dipartimento di Biotecnologie, Chimica e Farmacia, Università degli Studi di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Sergey I Bokarev
- Institut für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock. Germany
- Chemistry Department, Technical University of Munich, Lichtenbergstr. 4, 85748 Garching, Germany
| |
Collapse
|
5
|
Maity S, Kleinekathöfer U. Recent progress in atomistic modeling of light-harvesting complexes: a mini review. Photosynth Res 2023; 156:147-162. [PMID: 36207489 PMCID: PMC10070314 DOI: 10.1007/s11120-022-00969-w] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/30/2022] [Accepted: 09/22/2022] [Indexed: 06/16/2023]
Abstract
In this mini review, we focus on recent advances in the atomistic modeling of biological light-harvesting (LH) complexes. Because of their size and sophisticated electronic structures, multiscale methods are required to investigate the dynamical and spectroscopic properties of such complexes. The excitation energies, in this context also known as site energies, excitonic couplings, and spectral densities are key quantities which usually need to be extracted to be able to determine the exciton dynamics and spectroscopic properties. The recently developed multiscale approach based on the numerically efficient density functional tight-binding framework followed by excited state calculations has been shown to be superior to the scheme based on pure classical molecular dynamics simulations. The enhanced approach, which improves the description of the internal vibrational dynamics of the pigment molecules, yields spectral densities in good agreement with the experimental counterparts for various bacterial and plant LH systems. Here, we provide a brief overview of those results and described the theoretical foundation of the multiscale protocol.
Collapse
Affiliation(s)
- 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
|
6
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
7
|
Manicum AE, Louis H, Mathias GE, Agwamba EC, Malan FP, Unimuke TO, Nzondomyo WJ, Sithole SA, Biswas S, Prince S. Single crystal investigation, spectroscopic, DFT studies, and in-silico molecular docking of the anticancer activities of acetylacetone coordinated Re(I) tricarbonyl complexes. Inorganica Chim Acta 2023; 546:121335. [DOI: 10.1016/j.ica.2022.121335] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
|
8
|
Feighan O, Manby FR, Bourne-Worster S. An efficient protocol for excited states of large biochromophores. J Chem Phys 2023; 158:024107. [PMID: 36641400 DOI: 10.1063/5.0132417] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Abstract
Efficient energy transport in photosynthetic antenna is a long-standing source of inspiration for artificial light harvesting materials. However, characterizing the excited states of the constituent chromophores poses a considerable challenge to mainstream quantum chemical and semiempirical excited state methods due to their size and complexity and the accuracy required to describe small but functionally important changes in their properties. In this paper, we explore an alternative approach to calculating the excited states of large biochromophores, exemplified by a specific method for calculating the Qy transition of bacteriochlorophyll a, which we name Chl-xTB. Using a diagonally dominant approximation to the Casida equation and a bespoke parameterization scheme, Chl-xTB can match time-dependent density functional theory's accuracy and semiempirical speed for calculating the potential energy surfaces and absorption spectra of chlorophylls. We demonstrate that Chl-xTB (and other prospective realizations of our protocol) can be integrated into multiscale models, including concurrent excitonic and point-charge embedding frameworks, enabling the analysis of biochromophore networks in a native environment. We exploit this capability to probe the low-frequency spectral densities of excitonic energies and interchromophore interactions in the light harvesting antenna protein LH2 (light harvesting complex 2). The impact of low-frequency protein motion on interchromophore coupling and exciton transport has routinely been ignored due to the prohibitive costs of including it in simulations. Our results provide a more rigorous basis for continued use of this approximation by demonstrating that exciton transition energies are unaffected by low-frequency vibrational coupling to exciton interaction energies.
Collapse
Affiliation(s)
- Oliver Feighan
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Frederick R Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Susannah Bourne-Worster
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| |
Collapse
|
9
|
Mondragón-Solórzano G, Sandoval-Lira J, Nochebuena J, Cisneros GA, Barroso-Flores J. Electronic Structure Effects Related to the Origin of the Remarkable Near-Infrared Absorption of Blastochloris viridis' Light Harvesting 1-Reaction Center Complex. J Chem Theory Comput 2022; 18:4555-4564. [PMID: 35767461 PMCID: PMC10408377 DOI: 10.1021/acs.jctc.2c00497] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Various photosynthetic organisms have evolved to absorb light in different regions of the visible light spectrum, thus adapting to the various lighting conditions available on Earth. While most of these autotrophic organisms absorb wavelengths around the 700-800 nm region, some are capable of red-shifted absorptions above this range, but none as remarkably as Blastochloris viridis whose main absorption is observed at 1015 nm, approximately 220 nm (0.34 eV) lower in energy than their main constituent pigments, BChl-b, whose main absorption is observed at 795 nm. The structure of its light harvesting 1-reaction center was recently elucidated by cryo-EM; however, the electronic structure details behind this red-shifted absorption remain unattended. We used hybrid quantum mechanics/molecular mechanics (QM/MM) calculations to optimize one of the active centers and performed classical molecular dynamics (MD) simulations to sample conformations beyond the optimized structure. We did excited state calculations with the time-dependent density functional theory method at the CAM-B3LYP/cc-pVDZ level of theory. We reproduced the near IR absorption by sequentially modifying the number of components involved in our systems using representative structures from the calculated MD ensemble. Natural transition orbital analysis reveals the participation of the BChl-b fragments to the main transition in the native structure and the structures obtained from the QM/MM and MD simulations. H-bonding pigment-protein interactions play a role on the conformation stabilization and orientation; however, the bacteriochlorin ring conformations and the exciton delocalization are the most relevant factors to explain the red-shifting phenomenon.
Collapse
Affiliation(s)
- Gustavo Mondragón-Solórzano
- Centro Conjunto de Investigación en Química Sustentable UAEM-UNAM. Carretera Toluca-Atlacomulco Km. 14.5, Unidad San Cayetano. Toluca de Lerdo 50200, México
- Instituto de Química. Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, CDMX 04510, México
| | - Jacinto Sandoval-Lira
- Centro Conjunto de Investigación en Química Sustentable UAEM-UNAM. Carretera Toluca-Atlacomulco Km. 14.5, Unidad San Cayetano. Toluca de Lerdo 50200, México
- Instituto de Química. Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, CDMX 04510, México
- Departamento de Ingeniería Ambiental, Instituto Tecnológico Superior de San Martín Texmelucan, TecNM, Camino a la Barranca de Pesos, C.P. 74120 San Martín Texmelucan, Puebla, México
| | - Jorge Nochebuena
- Department of Physics, University of Texas at Dallas, Richardson, Texas 75801, United States
| | - G Andrés Cisneros
- Department of Physics, University of Texas at Dallas, Richardson, Texas 75801, United States
- Department of Chemistry and Biochemistry, University of Texas at Dallas, Richardson, Texas 75801, United States
| | - Joaquín Barroso-Flores
- Centro Conjunto de Investigación en Química Sustentable UAEM-UNAM. Carretera Toluca-Atlacomulco Km. 14.5, Unidad San Cayetano. Toluca de Lerdo 50200, México
- Instituto de Química. Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, CDMX 04510, México
| |
Collapse
|
10
|
Cho KH, Rhee YM. Computational elucidations on the role of vibrations in energy transfer processes of photosynthetic complexes. Phys Chem Chem Phys 2021; 23:26623-26639. [PMID: 34842245 DOI: 10.1039/d1cp04615b] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Coupling between pigment excitations and nuclear movements in photosynthetic complexes is known to modulate the excitation energy transfer (EET) efficiencies. Toward providing microscopic information, researchers often apply simulation techniques and investigate how vibrations are involved in EET processes. Here, reports on such roles of nuclear movements are discussed from a theory perspective. While vibrations naturally present random thermal fluctuations that can affect energy transferring characteristics, they can also be intertwined with exciton structures and create more specific non-adiabatic energy transfer pathways. For reliable simulations, a bath model that accurately mimics a given molecular system is required. Methods for obtaining such a model in combination with quantum chemical electronic structure calculations and molecular dynamics trajectory simulations are discussed. Various quantum dynamics simulation tools that can handle pigment-to-pigment energy transfers together with their vibrational characters are also touched on. Behaviors of molecular vibrations often deviate from ideality, especially when all-atom details are included, which practically forces us to treat them classically. We conclude this perspective by considering some recent reports that suggest that classical descriptions of bath effects with all-atom details may still produce valuable information for analyzing sophisticated contributions by vibrations to EET processes.
Collapse
Affiliation(s)
- Kwang Hyun Cho
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Korea.
| |
Collapse
|
11
|
Abstract
Theoretical prediction of electronic absorption spectra without input from experiments is no easy feat, as it requires addressing all of the factors that affect line shapes. In practice, however, the methodologies are limited to treat these ingredients only to a certain extent. Here, we present a multiscale protocol that addresses the temperature, solvent, and nuclear quantum effects as well as anharmonicity and the reconstruction of the final spectra from individual transitions. First, quantum mechanics/molecular mechanics (QM/MM) molecular dynamics is conducted to obtain trajectories of solute-solvent configurations, from which the corresponding quantum-corrected ensembles are generated through the generalized smoothed trajectory analysis (GSTA). The optical spectra of the ensembles are then produced by calculating vertical transitions using time-dependent density-functional theory (TDDFT) with implicit solvation. To obtain the final spectral shapes, the stick spectra from TDDFT are convoluted with Gaussian kernels where the half-widths are determined by a statistically motivated strategy. We have tested our method by calculating the UV-vis spectra of a recently discovered acridine photocatalyst in two redox states. Vibronic progressions and broadenings due to the finite lifetime of the excited states are not included in the methodology yet. Nuclear quantization affects the relative peak intensities and widths, which is necessary to reproduce the experimental spectrum. We have also found that using only the optimized geometry of each molecule works surprisingly well if a proper empirical broadening factor is applied. This is explained by the rigidity of the conjugated chromophore moieties of the selected molecules, which are mainly responsible for the excitations in the spectra. In contrast, we have also shown that other parts of the molecules are flexible enough to feature anharmonicities that impair the use of other techniques such as Wigner sampling.
Collapse
Affiliation(s)
- Péter P. Fehér
- Institute
of Organic Chemistry, Research Centre for
Natural Sciences, Magyar tudósok krt. 2, 1117 Budapest, Hungary
| | - Ádám Madarász
- Institute
of Organic Chemistry, Research Centre for
Natural Sciences, Magyar tudósok krt. 2, 1117 Budapest, Hungary
| | - András Stirling
- Institute
of Organic Chemistry, Research Centre for
Natural Sciences, Magyar tudósok krt. 2, 1117 Budapest, Hungary
- Department
of Chemistry, Eszterházy Károly
University, Leányka
u. 6, 3300 Eger, Hungary
| |
Collapse
|
12
|
Hashemi Z, Leppert L. Assessment of the Ab Initio Bethe-Salpeter Equation Approach for the Low-Lying Excitation Energies of Bacteriochlorophylls and Chlorophylls. J Phys Chem A 2021; 125:2163-2172. [PMID: 33656894 PMCID: PMC8028335 DOI: 10.1021/acs.jpca.1c01240] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Bacteriochlorophyll and chlorophyll molecules are crucial building blocks of the photosynthetic apparatus in bacteria, algae, and plants. Embedded in transmembrane protein complexes, they are responsible for the primary processes of photosynthesis: excitation energy and charge transfer. Here, we use ab initio many-body perturbation theory within the GW approximation and Bethe-Salpeter equation (BSE) approach to calculate the electronic structure and optical excitations of bacteriochlorophylls a, b, c, d, and e and chlorophylls a and b. We systematically study the effects of the structure, basis set size, partial self-consistency in GW, and the underlying exchange-correlation approximation and compare our calculations with results from time-dependent density functional theory, multireference RASPT2, and experimental literature results. We find that optical excitations calculated with GW+BSE are in excellent agreement with experimental data, with an average deviation of less than 100 meV for the first three bright excitations of the entire family of (bacterio)chlorophylls. Contrary to state-of-the-art time-dependent density functional theory (TDDFT) with an optimally tuned range-separated hybrid functional, this accuracy is achieved in a parameter-free approach. Moreover, GW+BSE predicts the energy differences between the low-energy excitations correctly and eliminates spurious charge transfer states that TDDFT with (semi)local approximations is known to produce. Our study provides accurate reference results and highlights the potential of the GW+BSE approach for the simulation of larger pigment complexes.
Collapse
Affiliation(s)
- Zohreh Hashemi
- Institute of Physics, University of Bayreuth, Bayreuth 95440, Germany
| | - Linn Leppert
- Institute of Physics, University of Bayreuth, Bayreuth 95440, Germany.,MESA+ Institute for Nanotechnology, University of Twente, 7500 AE Enschede, The Netherlands
| |
Collapse
|
13
|
Abstract
In photosynthetic complexes, tuning of chlorophyll light-absorption spectra by the protein environment is crucial to their efficiency and robustness. Recombinant type II water soluble chlorophyll-binding proteins from Brassicaceae (WSCPs) are useful for studying spectral tuning mechanisms due to their symmetric homotetramer structure, and the ability to rigorously modify the chlorophyll's protein surroundings. Our previous comparison of the crystal structures of two WSCP homologues suggested that protein-induced chlorophyll ring deformation is the predominant spectral tuning mechanism. Here, we implement a more rigorous analysis based on hybrid quantum mechanics and molecular mechanics calculations to quantify the relative contributions of geometrical and electrostatic factors to the absorption spectra of WSCP-chlorophyll complexes. We show that when considering conformational dynamics, geometry distortions such as chlorophyll ring deformation accounts for about one-third of the spectral shift, whereas the direct polarization of the electron density accounts for the remaining two-thirds. From a practical perspective, protein electrostatics is easier to manipulate than chlorophyll conformations, thus, it may be more readily implemented in designing artificial protein-chlorophyll complexes.
Collapse
Affiliation(s)
- Yigal Lahav
- Fritz Haber Centre for Molecular Dynamics Research, Institute of Chemistry, Hebrew University of Jerusalem, Israel.
| | | | | |
Collapse
|
14
|
Sirohiwal A, Neese F, Pantazis DA. How Can We Predict Accurate Electrochromic Shifts for Biochromophores? A Case Study on the Photosynthetic Reaction Center. J Chem Theory Comput 2021; 17:1858-1873. [PMID: 33566610 PMCID: PMC8023663 DOI: 10.1021/acs.jctc.0c01152] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2020] [Indexed: 01/28/2023]
Abstract
Protein-embedded chromophores are responsible for light harvesting, excitation energy transfer, and charge separation in photosynthesis. A critical part of the photosynthetic apparatus are reaction centers (RCs), which comprise groups of (bacterio)chlorophyll and (bacterio)pheophytin molecules that transform the excitation energy derived from light absorption into charge separation. The lowest excitation energies of individual pigments (site energies) are key for understanding photosynthetic systems, and form a prime target for quantum chemistry. A major theoretical challenge is to accurately describe the electrochromic (Stark) shifts in site energies produced by the inhomogeneous electric field of the protein matrix. Here, we present large-scale quantum mechanics/molecular mechanics calculations of electrochromic shifts for the RC chromophores of photosystem II (PSII) using various quantum chemical methods evaluated against the domain-based local pair natural orbital (DLPNO) implementation of the similarity-transformed equation of motion coupled cluster theory with single and double excitations (STEOM-CCSD). We show that certain range-separated density functionals (ωΒ97, ωΒ97X-V, ωΒ2PLYP, and LC-BLYP) correctly reproduce RC site energy shifts with time-dependent density functional theory (TD-DFT). The popular CAM-B3LYP functional underestimates the shifts and is not recommended. Global hybrid functionals are too insensitive to the environment and should be avoided, while nonhybrid functionals are strictly nonapplicable. Among the applicable approximate coupled cluster methods, the canonical versions of CC2 and ADC(2) were found to deviate significantly from the reference results both for the description of the lowest excited state and for the electrochromic shifts. By contrast, their spin-component-scaled (SCS) and particularly the scale-opposite-spin (SOS) variants compare well with the reference DLPNO-STEOM-CCSD and the best range-separated DFT methods. The emergence of RC excitation asymmetry is discussed in terms of intrinsic and protein electrostatic potentials. In addition, we evaluate a minimal structural scaffold of PSII, the D1-D2-CytB559 RC complex often employed in experimental studies, and show that it would have the same site energy distribution of RC chromophores as the full PSII supercomplex, but only under the unlikely conditions that the core protein organization and cofactor arrangement remain identical to those of the intact enzyme.
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
|
15
|
Osella S. Artificial Photosynthesis: Is Computation Ready for the Challenge Ahead? Nanomaterials (Basel) 2021; 11:nano11020299. [PMID: 33498961 PMCID: PMC7911014 DOI: 10.3390/nano11020299] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/07/2021] [Revised: 01/19/2021] [Accepted: 01/20/2021] [Indexed: 12/13/2022]
Abstract
A tremendous effort is currently devoted to the generation of novel hybrid materials with enhanced electronic properties for the creation of artificial photosynthetic systems. This compelling and challenging problem is well-defined from an experimental point of view, as the design of such materials relies on combining organic materials or metals with biological systems like light harvesting and redox-active proteins. Such hybrid systems can be used, e.g., as bio-sensors, bio-fuel cells, biohybrid photoelectrochemical cells, and nanostructured photoelectronic devices. Despite these efforts, the main bottleneck is the formation of efficient interfaces between the biological and the organic/metal counterparts for efficient electron transfer (ET). It is within this aspect that computation can make the difference and improve the current understanding of the mechanisms underneath the interface formation and the charge transfer efficiency. Yet, the systems considered (i.e., light harvesting protein, self-assembly monolayer and surface assembly) are more and more complex, reaching (and often passing) the limit of current computation power. In this review, recent developments in computational methods for studying complex interfaces for artificial photosynthesis will be provided and selected cases discussed, to assess the inherent ability of computation to leave a mark in this field of research.
Collapse
Affiliation(s)
- Silvio Osella
- Chemical and Biological Systems Simulation Lab, Center of New Technologies, University of Warsaw, Banacha 2C, 02-097 Warsaw, Poland
| |
Collapse
|
16
|
Maity S, Bold BM, Prajapati JD, Sokolov M, Kubař T, Elstner M, Kleinekathöfer U. DFTB/MM Molecular Dynamics Simulations of the FMO Light-Harvesting Complex. J Phys Chem Lett 2020; 11:8660-8667. [PMID: 32991176 DOI: 10.1021/acs.jpclett.0c02526] [Citation(s) in RCA: 33] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Because of the size of light-harvesting complexes and the involvement of electronic degrees of freedom, computationally these systems need to be treated with a combined quantum-classical description. To this end, Born-Oppenheimer molecular dynamics simulations have been employed in a quantum mechanics/molecular mechanics (QM/MM) fashion for the ground state followed by excitation energy calculations again in a QM/MM scheme for the Fenna-Matthews-Olson (FMO) complex. The self-consistent-charge density functional tight-binding (DFTB) method electrostatically coupled to a classical description of the environment was applied to perform the ground-state dynamics. Subsequently, long-range-corrected time-dependent DFTB calculations were performed to determine the excitation energy fluctuations of the individual bacteriochlorophyll a molecules. The spectral densities obtained using this approach show an excellent agreement with experimental findings. In addition, the fluctuating site energies and couplings were used to estimate the exciton transfer dynamics.
Collapse
Affiliation(s)
- Sayan Maity
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Beatrix M Bold
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | | | - Monja Sokolov
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Tomáš Kubař
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Marcus Elstner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
- Institute of Biological Interfaces, Karlsruhe Institute of Technology, Kaiserstrasse 12, 76131 Karlsruhe, Germany
| | - Ulrich Kleinekathöfer
- Department of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| |
Collapse
|
17
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
18
|
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] [What about the content of this article? (0)] [Affiliation(s)] [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
|
19
|
Bold BM, Sokolov M, Maity S, Wanko M, Dohmen PM, Kranz JJ, Kleinekathöfer U, Höfener S, Elstner M. Benchmark and performance of long-range corrected time-dependent density functional tight binding (LC-TD-DFTB) on rhodopsins and light-harvesting complexes. Phys Chem Chem Phys 2020; 22:10500-10518. [PMID: 31950960 DOI: 10.1039/c9cp05753f] [Citation(s) in RCA: 34] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The chromophores of rhodopsins (Rh) and light-harvesting (LH) complexes still represent a major challenge for a quantum chemical description due to their size and complex electronic structure. Since gradient corrected and hybrid density functional approaches have been shown to fail for these systems, only range-separated functionals seem to be a promising alternative to the more time consuming post-Hartree-Fock approaches. For extended sampling of optical properties, however, even more approximate approaches are required. Recently, a long-range corrected (LC) functional has been implemented into the efficient density functional tight binding (DFTB) method, allowing to sample the excited states properties of chromophores embedded into proteins using quantum mechanical/molecular mechanical (QM/MM) with the time-dependent (TD) DFTB approach. In the present study, we assess the accuracy of LC-TD-DFT and LC-TD-DFTB for rhodopsins (bacteriorhodopsin (bR) and pharaonis phoborhodopsin (ppR)) and LH complexes (light-harvesting complex II (LH2) and Fenna-Matthews-Olson (FMO) complex). This benchmark study shows the improved description of the color tuning parameters compared to standard DFT functionals. In general, LC-TD-DFTB can exhibit a similar performance as the corresponding LC functionals, allowing a reliable description of excited states properties at significantly reduced cost. The two chromophores investigated here pose complementary challenges: while huge sensitivity to external field perturbation (color tuning) and charge transfer excitations are characteristic for the retinal chromophore, the multi-chromophoric character of the LH complexes emphasizes a correct description of inter-chromophore couplings, giving less importance to color tuning. None of the investigated functionals masters both systems simultaneously with satisfactory accuracy. LC-TD-DFTB, at the current stage, although showing a systematic improvement compared to TD-DFTB cannot be recommended for studying color tuning in retinal proteins, similar to some of the LC-DFT functionals, because the response to external fields is still too weak. For sampling of LH-spectra, however, LC-TD-DFTB is a viable tool, allowing to efficiently sample absorption energies, as shown for three different LH complexes. As the calculations indicate, geometry optimization may overestimate the importance of local minima, which may be averaged over when using trajectories. Fast quantum chemical approaches therefore may allow for a direct sampling of spectra in the near future.
Collapse
Affiliation(s)
- Beatrix M Bold
- Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Kaiserstrasse 12, 76131 Karlsruhe, Germany.
| | | | | | | | | | | | | | | | | |
Collapse
|
20
|
Saito S, Higashi M, Fleming GR. Site-Dependent Fluctuations Optimize Electronic Energy Transfer in the Fenna–Matthews–Olson Protein. J Phys Chem B 2019; 123:9762-9772. [DOI: 10.1021/acs.jpcb.9b07456] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Affiliation(s)
- Shinji Saito
- Institute for Molecular Science, Myodaiji, Okazaki, Aichi 444-8585, Japan
- The Graduate University for Advanced Studies, Myodaiji, Okazaki, Aichi 444-8585, Japan
| | - Masahiro Higashi
- Department of Molecular Engineering, Graduate School of Engineering, Kyoto University, Kyoto daigaku-katsura, Nishikyo-ku, Kyoto 615-8510, Kyoto, Japan
- JST, PRESTO, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| | - Graham R. Fleming
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Molecular Biophysics and Integrated Bioimaging Division, Lawrence Berkeley National Lab, Berkeley, California 94720, United States
- Kavli Energy Nanosciences Institute at Berkeley, Berkeley, California 94720, United States
| |
Collapse
|
21
|
|
22
|
Cupellini L, Bondanza M, Nottoli M, Mennucci B. Successes & challenges in the atomistic modeling of light-harvesting and its photoregulation. Biochim Biophys Acta Bioenerg 2019; 1861:148049. [PMID: 31386831 DOI: 10.1016/j.bbabio.2019.07.004] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/10/2019] [Revised: 07/10/2019] [Accepted: 07/15/2019] [Indexed: 12/21/2022]
Abstract
Light-harvesting is a crucial step of photosynthesis. Its mechanisms and related energetics have been revealed by a combination of experimental investigations and theoretical modeling. The success of theoretical modeling is largely due to the application of atomistic descriptions combining quantum chemistry, classical models and molecular dynamics techniques. Besides the important achievements obtained so far, a complete and quantitative understanding of how the many different light-harvesting complexes exploit their structural specificity is still missing. Moreover, many questions remain unanswered regarding the mechanisms through which light-harvesting is regulated in response to variable light conditions. Here we show that, in both fields, a major role will be played once more by atomistic descriptions, possibly generalized to tackle the numerous time and space scales on which the regulation takes place: going from the ultrafast electronic excitation of the multichromophoric aggregate, through the subsequent conformational changes in the embedding protein, up to the interaction between proteins.
Collapse
Affiliation(s)
- Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy
| | - Mattia Bondanza
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy
| | - Michele Nottoli
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, Pisa 56124, Italy.
| |
Collapse
|
23
|
Affiliation(s)
- André Anda
- Chemical and Quantum Physics, School of Science, RMIT University, Melbourne, VIC 3001, Australia
- ARC Centre of Excellence in Exciton Science, School of Science, RMIT University, Melbourne, VIC 3001, Australia
| | - Thorsten Hansen
- Department of Chemistry, Copenhagen University, Universitetsparken 5, DK-2100, Copenhagen Ø, Denmark
| | - Luca De Vico
- Department of Biotechnologies, Chemistry and Pharmacy, University of Siena, Via Aldo Moro 2, I-53100, Siena, Italy
| |
Collapse
|
24
|
De Vico L, Anda A, Osipov VA, Madsen AØ, Hansen T. Macrocycle ring deformation as the secondary design principle for light-harvesting complexes. Proc Natl Acad Sci U S A 2018; 115:E9051-7. [PMID: 30194231 DOI: 10.1073/pnas.1719355115] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Natural light-harvesting is performed by pigment-protein complexes, which collect and funnel the solar energy at the start of photosynthesis. The identity and arrangement of pigments largely define the absorption spectrum of the antenna complex, which is further regulated by a palette of structural factors. Small alterations are induced by pigment-protein interactions. In light-harvesting systems 2 and 3 from Rhodoblastus acidophilus, the pigments are arranged identically, yet the former has an absorption peak at 850 nm that is blue-shifted to 820 nm in the latter. While the shift has previously been attributed to the removal of hydrogen bonds, which brings changes in the acetyl moiety of the bacteriochlorophyll, recent work has shown that other mechanisms are also present. Using computational and modeling tools on the corresponding crystal structures, we reach a different conclusion: The most critical factor for the shift is the curvature of the macrocycle ring. The bending of the planar part of the pigment is identified as the second-most important design principle for the function of pigment-protein complexes-a finding that can inspire the design of novel artificial systems.
Collapse
|
25
|
Kozlov MI, Poddubnyy VV, Glebov IO. Calculation of the excited states properties of LH1 complex of Thermochromatium tepidum. J Comput Chem 2018; 39:1599-1606. [PMID: 29701316 DOI: 10.1002/jcc.25232] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/12/2017] [Revised: 03/28/2018] [Accepted: 04/03/2018] [Indexed: 11/09/2022]
Abstract
Calculation of the excited states properties of pigment complexes is one of the key problems in the photosynthesis research. The excited states of LH1 complex of Thermochromatium tepidum were studied by means of the high-precision quantum chemistry methods. The influence of different parameters of the calculation procedure was examined. The optimal scheme of calculation was chosen by comparison of calculated results with the experimental data on absorption, electronic and magnetic circular dichroism spectra. The high importance of the account of the second excited states of bacteriochlorophylls and of site heterogeneity was shown. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Maxim I Kozlov
- Department of Chemistry, Lomonosov Moscow State University, Moscow, Russia, 119991
| | - Vladimir V Poddubnyy
- Department of Chemistry, Lomonosov Moscow State University, Moscow, Russia, 119991
| | - Ilya O Glebov
- Department of Chemistry, Lomonosov Moscow State University, Moscow, Russia, 119991
| |
Collapse
|
26
|
Kim CW, Choi B, Rhee YM. Excited state energy fluctuations in the Fenna-Matthews-Olson complex from molecular dynamics simulations with interpolated chromophore potentials. Phys Chem Chem Phys 2018; 20:3310-3319. [PMID: 29186231 DOI: 10.1039/c7cp06303b] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
We analyze the environment-induced fluctuation of pigment excitation energies in the Fenna-Matthews-Olson (FMO) complex from various perspectives, by employing an interpolation-based all-atom potential energy model for describing realistic pigment vibrations. We conduct molecular dynamics simulations on a 100 ns timescale, which is an extent that can enclose the effect of static disorder, and demonstrate its timescale separation from fast dynamic disorder. We extract the spectral densities of the complex by considering both the site and the exciton bases. We show that exciton delocalization reduces the effective environmental fluctuation and rationalize this aspect based on a model of fluctuating molecular aggregates. We also obtained the spectral density of the lowest exciton state under low temperature conditions and show that it reasonably well reproduces the experimental result. Finally, by additionally performing non-equilibrium excited state trajectory simulations, we show that the system lies well within the linear response regime after photo-absorption and that the pigments do not visit anharmonic regions of the potential surface to a significant extent. This indicates that methodologies based on harmonic bath models are indeed reasonable approaches for describing the excited state dynamics of the FMO complex.
Collapse
Affiliation(s)
- Chang Woo Kim
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang 37673, Korea
| | | | | |
Collapse
|
27
|
Abstract
In this paper, we introduce two simple quantum dynamics methods. One is based on the popular surface-hopping method, and the other is based on rescaling of the propagation on the bath ground-state potential surface. The first method is special, as it avoids specific feedback from the simulated quantum system to the bath and can be applied for precalculated classical trajectories. It is based on the equipartition theorem to determine if hops between different potential energy surfaces are allowed. By comparing with the formally exact Hierarchical Equations Of Motion approach for four model systems we find that the method generally approximates the quantum dynamics toward thermal equilibrium very well. The second method is based on rescaling of the nonadiabatic coupling and also neglect the effect of the state of the quantum system on the bath. By the nature of the approximations, they cannot reproduce the effect of bath relaxation following excitation. However, the methods are both computationally more tractable than the conventional fewest switches surface hopping, and we foresee that the methods will be powerful for simulations of quantum dynamics in systems with complex bath dynamics, where the system-bath coupling is not too strong compared to the thermal energy.
Collapse
Affiliation(s)
- Thomas L. C. Jansen
- Zernike Institute for Advanced
Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| |
Collapse
|
28
|
Abstract
In this review, we discuss our recent work on modelling biological pigment-protein complexes, such as the Fenna-Matthews-Olson complex and light-harvesting complex-II, to explain their electronic energy transport properties. In particular, we highlight how a network-based analysis approach, where the light-absorbing pigments are treated as a network of interconnected nodes, can provide a qualitative picture of quantum dynamic energy transport. With this in mind, we demonstrate how other properties such as robustness to environmental changes can be assessed in a simple and computationally tractable manner. Such analyses could prove useful for the design of artificial energy transport networks such as those which might find application in solar cells.
Collapse
|
29
|
Anda A, De Vico L, Hansen T. Intermolecular Modes between LH2 Bacteriochlorophylls and Protein Residues: The Effect on the Excitation Energies. J Phys Chem B 2017; 121:5499-5508. [PMID: 28485594 DOI: 10.1021/acs.jpcb.7b02071] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Light-harvesting system 2 (LH2) executes the primary processes of photosynthesis in purple bacteria; photon absorption, and energy transportation to the reaction center. A detailed mechanistic insight into these operations is obscured by the complexity of the light-harvesting systems, particularly by the chromophore-environment interaction. In this work, we focus on the effects of the protein residues that are ligated to the bacteriochlorophylls (BChls) and construct potential energy surfaces of the ground and first optically excited state for the various BChl-residue systems where we in each case consider two degrees of freedom in the intermolecular region. We find that the excitation energies are only slightly affected by the considered modes. In addition, we see that axial ligands and hydrogen-bonded residues have opposite effects on both excitation energies and oscillator strengths by comparing to the isolated BChls. Our results indicate that only a small part of the chromophore-environment interaction can be associated with the intermolecular region between a BChl and an adjacent residue, but that it may be possible to selectively raise or lower the excitation energy at the axial and planar residue positions, respectively.
Collapse
Affiliation(s)
- André Anda
- Department of Chemistry, H. C. Ørsted Institute, University of Copenhagen , Universitetsparken 5, DK-2100 Copenhagen, Denmark
| | - Luca De Vico
- Department of Chemistry, H. C. Ørsted Institute, University of Copenhagen , Universitetsparken 5, DK-2100 Copenhagen, Denmark.,Department of Biotechnology, Chemistry and Pharmacy, University of Siena , via Aldo Moro 2, 53100 Siena, Italy
| | - Thorsten Hansen
- Department of Chemistry, H. C. Ørsted Institute, University of Copenhagen , Universitetsparken 5, DK-2100 Copenhagen, Denmark
| |
Collapse
|
30
|
Varsano D, Caprasecca S, Coccia E. Theoretical description of protein field effects on electronic excitations of biological chromophores. J Phys Condens Matter 2017; 29:013002. [PMID: 27830666 DOI: 10.1088/0953-8984/29/1/013002] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Photoinitiated phenomena play a crucial role in many living organisms. Plants, algae, and bacteria absorb sunlight to perform photosynthesis, and convert water and carbon dioxide into molecular oxygen and carbohydrates, thus forming the basis for life on Earth. The vision of vertebrates is accomplished in the eye by a protein called rhodopsin, which upon photon absorption performs an ultrafast isomerisation of the retinal chromophore, triggering the signal cascade. Many other biological functions start with the photoexcitation of a protein-embedded pigment, followed by complex processes comprising, for example, electron or excitation energy transfer in photosynthetic complexes. The optical properties of chromophores in living systems are strongly dependent on the interaction with the surrounding environment (nearby protein residues, membrane, water), and the complexity of such interplay is, in most cases, at the origin of the functional diversity of the photoactive proteins. The specific interactions with the environment often lead to a significant shift of the chromophore excitation energies, compared with their absorption in solution or gas phase. The investigation of the optical response of chromophores is generally not straightforward, from both experimental and theoretical standpoints; this is due to the difficulty in understanding diverse behaviours and effects, occurring at different scales, with a single technique. In particular, the role played by ab initio calculations in assisting and guiding experiments, as well as in understanding the physics of photoactive proteins, is fundamental. At the same time, owing to the large size of the systems, more approximate strategies which take into account the environmental effects on the absorption spectra are also of paramount importance. Here we review the recent advances in the first-principle description of electronic and optical properties of biological chromophores embedded in a protein environment. We show their applications on paradigmatic systems, such as the light-harvesting complexes, rhodopsin and green fluorescent protein, emphasising the theoretical frameworks which are of common use in solid state physics, and emerging as promising tools for biomolecular systems.
Collapse
Affiliation(s)
- Daniele Varsano
- S3 Center, CNR Institute of Nanoscience, Via Campi 213/A, 41125 Modena, Italy
| | | | | |
Collapse
|
31
|
Sisto A, Stross C, van der Kamp MW, O’Connor M, McIntosh-Smith S, Johnson GT, Hohenstein EG, Manby FR, Glowacki DR, Martinez TJ. Atomistic non-adiabatic dynamics of the LH2 complex with a GPU-accelerated ab initio exciton model. Phys Chem Chem Phys 2017; 19:14924-14936. [DOI: 10.1039/c7cp00492c] [Citation(s) in RCA: 50] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
We present GPU-accelerated ab initio molecular dynamics simulations of nonadiabatic dynamics in the LH2 complex in full atomistic detail.
Collapse
Affiliation(s)
- Aaron Sisto
- PULSE Institute and Department of Chemistry
- Stanford University
- Stanford
- USA
- SLAC National Accelerator Laboratory
| | - Clem Stross
- School of Chemistry
- University of Bristol
- Bristol
- UK
| | | | - Michael O’Connor
- School of Chemistry
- University of Bristol
- Bristol
- UK
- Department of Computer Science
| | | | - Graham T. Johnson
- California Institute for Quantitative Biosciences (QB3)
- University of California
- San Francisco
- USA
- Department of Bioengineering and Therapeutic Sciences
| | | | | | - David R. Glowacki
- School of Chemistry
- University of Bristol
- Bristol
- UK
- Department of Computer Science
| | - Todd J. Martinez
- PULSE Institute and Department of Chemistry
- Stanford University
- Stanford
- USA
- SLAC National Accelerator Laboratory
| |
Collapse
|
32
|
Morgan SE, Cole DJ, Chin AW. Nonlinear network model analysis of vibrational energy transfer and localisation in the Fenna-Matthews-Olson complex. Sci Rep 2016; 6:36703. [PMID: 27827409 PMCID: PMC5101523 DOI: 10.1038/srep36703] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2016] [Accepted: 10/19/2016] [Indexed: 11/10/2022] Open
Abstract
Collective protein modes are expected to be important for facilitating energy transfer in the Fenna-Matthews-Olson (FMO) complex of photosynthetic green sulphur bacteria, however to date little work has focussed on the microscopic details of these vibrations. The nonlinear network model (NNM) provides a computationally inexpensive approach to studying vibrational modes at the microscopic level in large protein structures, whilst incorporating anharmonicity in the inter-residue interactions which can influence protein dynamics. We apply the NNM to the entire trimeric FMO complex and find evidence for the existence of nonlinear discrete breather modes. These modes tend to transfer energy to the highly connected core pigments, potentially opening up alternative excitation energy transfer routes through their influence on pigment properties. Incorporating localised modes based on these discrete breathers in the optical spectra calculations for FMO using ab initio site energies and excitonic couplings can substantially improve their agreement with experimental results.
Collapse
Affiliation(s)
- Sarah E Morgan
- Theory of Condensed Matter Group, Physics Department, University of Cambridge, CB3 0HE, United Kingdom
| | - Daniel J Cole
- Theory of Condensed Matter Group, Physics Department, University of Cambridge, CB3 0HE, United Kingdom.,School of Chemistry, Newcastle University, Newcastle upon Tyne, NE1 7RU, United Kingdom
| | - Alex W Chin
- Theory of Condensed Matter Group, Physics Department, University of Cambridge, CB3 0HE, United Kingdom
| |
Collapse
|
33
|
Cupellini L, Jurinovich S, Campetella M, Caprasecca S, Guido CA, Kelly SM, Gardiner AT, Cogdell R, Mennucci B. An Ab Initio Description of the Excitonic Properties of LH2 and Their Temperature Dependence. J Phys Chem B 2016; 120:11348-11359. [PMID: 27791372 DOI: 10.1021/acs.jpcb.6b06585] [Citation(s) in RCA: 52] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
The spectroscopic properties of light-harvesting (LH) antennae in photosyntehtic organisms represent a fingerprint that is unique for each specific pigment-protein complex. Because of that, spectroscopic observations are generally combined with structural data from X-ray crystallography to obtain an indirect representation of the excitonic properties of the system. Here, an alternative strategy is presented which goes beyond this empirical approach and introduces an ab initio computational description of both structural and electronic properties and their dependence on the temperature. The strategy is applied to the peripheral light-harvesting antenna complex (LH2) present in purple bacteria. By comparing this model with the one based on the crystal structure, a detailed, molecular level explanation of the absorption and circular dichroism (CD) spectra and their temperature dependence is achieved. The agreement obtained with the experiments at both low and room temperature lays the groundwork for an atomistic understanding of the excitation dynamics in the LH2 system.
Collapse
Affiliation(s)
- Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa , Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Sandro Jurinovich
- Dipartimento di Chimica e Chimica Industriale, University of Pisa , Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Marco Campetella
- Dipartimento di Chimica e Chimica Industriale, University of Pisa , Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Stefano Caprasecca
- Dipartimento di Chimica e Chimica Industriale, University of Pisa , Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Ciro A Guido
- Dipartimento di Chimica e Chimica Industriale, University of Pisa , Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Sharon M Kelly
- Life Sciences Biomolecular Sci, Joseph Black Building, University of Glasgow , Glasgow G12 8QQ, Scotland
| | - Alastair T Gardiner
- Glasgow Biomedical Research Centre, Institute of Molecular Cell and Systems Biology, University of Glasgow , 126 University Place, Glasgow G12 8TA, Scotland
| | - Richard Cogdell
- Glasgow Biomedical Research Centre, Institute of Molecular Cell and Systems Biology, University of Glasgow , 126 University Place, Glasgow G12 8TA, Scotland
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, University of Pisa , Via G. Moruzzi 13, I-56124 Pisa, Italy
| |
Collapse
|
34
|
Abstract
Obtaining the exciton dynamics of large photosynthetic complexes by using mixed quantum mechanics/molecular mechanics (QM/MM) is computationally demanding. We propose a machine learning technique, multi-layer perceptrons, as a tool to reduce the time required to compute excited state energies. With this approach we predict time-dependent density functional theory (TDDFT) excited state energies of bacteriochlorophylls in the Fenna-Matthews-Olson (FMO) complex. Additionally we compute spectral densities and exciton populations from the predictions. Different methods to determine multi-layer perceptron training sets are introduced, leading to several initial data selections. In addition, we compute spectral densities and exciton populations. Once multi-layer perceptrons are trained, predicting excited state energies was found to be significantly faster than the corresponding QM/MM calculations. We showed that multi-layer perceptrons can successfully reproduce the energies of QM/MM calculations to a high degree of accuracy with prediction errors contained within 0.01 eV (0.5%). Spectral densities and exciton dynamics are also in agreement with the TDDFT results. The acceleration and accurate prediction of dynamics strongly encourage the combination of machine learning techniques with ab initio methods.
Collapse
Affiliation(s)
- Florian Häse
- Department of Chemistry and Chemical Biology , Harvard University , Cambridge , 02138 , USA . ; Tel: +1-617-384-8188
- Physik-Department T38 , Technische Universität München , Garching , 85748 , Germany
| | - Stéphanie Valleau
- Department of Chemistry and Chemical Biology , Harvard University , Cambridge , 02138 , USA . ; Tel: +1-617-384-8188
| | - Edward Pyzer-Knapp
- Department of Chemistry and Chemical Biology , Harvard University , Cambridge , 02138 , USA . ; Tel: +1-617-384-8188
| | - Alán Aspuru-Guzik
- Department of Chemistry and Chemical Biology , Harvard University , Cambridge , 02138 , USA . ; Tel: +1-617-384-8188
| |
Collapse
|
35
|
Affiliation(s)
- Albrecht Goez
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster , Münster, Germany
| | - Johannes Neugebauer
- Theoretische Organische Chemie, Organisch-Chemisches Institut and Center for Multiscale Theory and Computation, Westfälische Wilhelms-Universität Münster , Münster, Germany
| |
Collapse
|
36
|
Affiliation(s)
- Mi Kyung Lee
- Department of Chemistry, Boston University, Boston, Massachusetts 02215;
| | - Pengfei Huo
- Department of Chemistry, University of Rochester, Rochester, New York 14627;
| | - David F. Coker
- Department of Chemistry, Boston University, Boston, Massachusetts 02215;
| |
Collapse
|
37
|
Abstract
The design of optimal light-harvesting (supra)molecular systems and materials is one of the most challenging frontiers of science. Theoretical methods and computational models play a fundamental role in this difficult task, as they allow the establishment of structural blueprints inspired by natural photosynthetic organisms that can be applied to the design of novel artificial light-harvesting devices. Among theoretical strategies, the application of quantum chemical tools represents an important reality that has already reached an evident degree of maturity, although it still has to show its real potentials. This Review presents an overview of the state of the art of this strategy, showing the actual fields of applicability but also indicating its current limitations, which need to be solved in future developments.
Collapse
Affiliation(s)
- Carles Curutchet
- Departament de Fisicoquímica, Facultat de Farmàcia, Universitat de Barcelona , Av. Joan XXIII s/n, 08028 Barcelona, Spain
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, University of Pisa , via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
38
|
Anda A, Hansen T, De Vico L. Multireference Excitation Energies for Bacteriochlorophylls A within Light Harvesting System 2. J Chem Theory Comput 2016; 12:1305-13. [DOI: 10.1021/acs.jctc.5b01104] [Citation(s) in RCA: 28] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- André Anda
- Department of Chemistry, Copenhagen University, Universitetsparken 5, DK-2100, Copenhagen Ø, Denmark
| | - Thorsten Hansen
- Department of Chemistry, Copenhagen University, Universitetsparken 5, DK-2100, Copenhagen Ø, Denmark
| | - Luca De Vico
- Department of Chemistry, Copenhagen University, Universitetsparken 5, DK-2100, Copenhagen Ø, Denmark
| |
Collapse
|
39
|
MacGowan SA, Senge MO. Contribution of bacteriochlorophyll conformation to the distribution of site-energies in the FMO protein. Biochim Biophys Acta 2016; 1857:427-42. [PMID: 26851682 DOI: 10.1016/j.bbabio.2016.02.001] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/04/2015] [Revised: 01/30/2016] [Accepted: 02/01/2016] [Indexed: 10/22/2022]
Abstract
The structural data for the Fenna-Matthews-Olson (FMO) protein indicate that the bacteriochlorophylls (BChls) display a significant degree of conformational heterogeneity of their peripheral substituents and the protein-induced nonplanar skeletal deformations of the tetrapyrrole macrocycle. As electronic properties of chromophores are altered by such differences, a conformational effect may influence the site-energies of specific pigments and thus play a role in mediating the excitation energy transfer dynamics, but this has not yet been established. The difficulty of assessing this question is shown to be partly the result of the inability of the sequential truncation approach usually employed to account for interactions between the conformations of the macrocycle and its substituents and an alternative approach is suggested. By assigning the BChl atoms to meaningful atom groups and performing all possible permutations of partial optimizations in a full-factorial design, where each group is either frozen in the crystal geometry or optimized in vacuo, followed by excited state calculations on each resulting structure (PM6//ZIndo/S), the specific effects of the conformations of each BChl component as well as mutual interactions between the molecular fragments on the site-energy can be delineated. This factorial relaxation procedure gives different estimates of the macrocycle conformational perturbation than the approach of sequentially truncating the BChl periphery. The results were evaluated in the context of published site-energies for the FMO pigments from three species to identify how conformational effects contribute to their distribution and instances of cross-species conservation and functional divergence of the BChl nonplanarity conformational contribution are described.
Collapse
Affiliation(s)
- Stuart A MacGowan
- School of Chemistry, SFI Tetrapyrrole Laboratory, Trinity Biomedical Sciences Institute, 152-160 Pearse Street, Trinity College Dublin, The University of Dublin, Dublin 2, Ireland
| | - Mathias O Senge
- School of Chemistry, SFI Tetrapyrrole Laboratory, Trinity Biomedical Sciences Institute, 152-160 Pearse Street, Trinity College Dublin, The University of Dublin, Dublin 2, Ireland.
| |
Collapse
|
40
|
Jia X, Mei Y, Zhang JZH, Mo Y. Hybrid QM/MM study of FMO complex with polarized protein-specific charge. Sci Rep 2015; 5:17096. [PMID: 26611739 PMCID: PMC4661465 DOI: 10.1038/srep17096] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2015] [Accepted: 10/26/2015] [Indexed: 12/20/2022] Open
Abstract
The Fenna-Matthews-Olson (FMO) light-harvesting complex is now one of the primary model systems for the study of excitation energy transfer (EET). However, the mechanism of the EET in this system is still controversial. In this work, molecular dynamics simulations and the electrostatic-embedding quantum-mechanics/molecular-mechanics single-point calculations have been employed to predict the energy transfer pathways utilizing the polarized protein-specific charge (PPC), which provides a more realistic description of Coulomb interaction potential in the protein than conventional mean-field charge scheme. The recently discovered eighth pigment has also been included in this study. Comparing with the conventional mean-field charges, more stable structures of FMO complex were found under PPC scheme during molecular dynamic simulation. Based on the electronic structure calculations, an exciton model was constructed to consider the couplings during excitation. The results show that pigments 3 and 4 dominate the lowest exciton levels whereas the highest exciton level are mainly constituted of pigments 1 and 6. This observation agrees well with the assumption based on the spatial distribution of the pigments. Moreover, the obtained spectral density in this study gives a reliable description of the diverse local environment embedding each pigment.
Collapse
Affiliation(s)
- Xiangyu Jia
- State Key Laboratory of Precision Spectroscopy and Department of Physics and Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy and Department of Physics and Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
| | - John Z H Zhang
- State Key Laboratory of Precision Spectroscopy and Department of Physics and Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China.,Department of Chemistry, New York University, New York, NY 10003
| | - Yan Mo
- State Key Laboratory of Precision Spectroscopy and Department of Physics and Institute of Theoretical and Computational Science, East China Normal University, Shanghai 200062, China
| |
Collapse
|
41
|
Steinmann C, Kongsted J. Electronic Energy Transfer in Polarizable Heterogeneous Environments: A Systematic Investigation of Different Quantum Chemical Approaches. J Chem Theory Comput 2015; 11:4283-93. [DOI: 10.1021/acs.jctc.5b00470] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Casper Steinmann
- Department of Physics, Chemistry,
and Pharmacy, University of Southern Denmark, DK-5230 Odense
M, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry,
and Pharmacy, University of Southern Denmark, DK-5230 Odense
M, Denmark
| |
Collapse
|
42
|
Chandrasekaran S, Aghtar M, Valleau S, Aspuru-Guzik A, Kleinekathöfer U. Influence of Force Fields and Quantum Chemistry Approach on Spectral Densities of BChl a in Solution and in FMO Proteins. J Phys Chem B 2015; 119:9995-10004. [DOI: 10.1021/acs.jpcb.5b03654] [Citation(s) in RCA: 75] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
| | - Mortaza Aghtar
- Department
of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| | - Stéphanie Valleau
- Department
of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Alán Aspuru-Guzik
- Department
of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Ulrich Kleinekathöfer
- Department
of Physics and Earth Sciences, Jacobs University Bremen, Campus Ring 1, 28759 Bremen, Germany
| |
Collapse
|
43
|
Jurinovich S, Viani L, Curutchet C, Mennucci B. Limits and potentials of quantum chemical methods in modelling photosynthetic antennae. Phys Chem Chem Phys 2015; 17:30783-92. [DOI: 10.1039/c5cp00986c] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Abstract
A critical overview of quantum chemical approaches to simulate the light-harvesting process in photosynthetic antennae is presented together with a perspective on the developments that need to be introduced to reach a quantitative predictive power.
Collapse
Affiliation(s)
- Sandro Jurinovich
- Dipartimento di Chimica e Chimica Industriale
- University of Pisa
- 56124 Pisa
- Italy
| | - Lucas Viani
- Dipartimento di Chimica e Chimica Industriale
- University of Pisa
- 56124 Pisa
- Italy
- Institute for Fluid Dynamics
| | - Carles Curutchet
- Departament de Fisicoquímica
- Facultat de Farmàcia
- Universitat de Barcelona
- 08028 Barcelona
- Spain
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale
- University of Pisa
- 56124 Pisa
- Italy
| |
Collapse
|
44
|
Higashi M, Kosugi T, Hayashi S, Saito S. Theoretical study on excited states of bacteriochlorophyll a in solutions with density functional assessment. J Phys Chem B 2014; 118:10906-18. [PMID: 25153487 DOI: 10.1021/jp507259g] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The excited-state properties of bacteriochlorophyll (BChl) a in triethylamine, 1-propanol, and methanol are investigated with the time-dependent density functional theory by using the quantum mechanical and molecular mechanical reweighting free energy self-consistant field method. It is found that no prevalent density functionals can reproduce the experimental excited-state properties, i.e., the absorption and reorganization energies, of BChl a in the solutions. The parameter μ in the range-separated hybrid functional is therefore optimized to reproduce the differences of the absorption energies in the solutions. We examine the origin of the differences of the absorption energies in the solutions and find that sensitive balance between contributions of structural changes and solute-solvent interactions determines the differences. The accurate description of the excitation with the density functional with the adjusted parameter is therefore essential to the understanding of the excited-state properties of BChl a in proteins and also the mechanism of the photosynthetic systems.
Collapse
Affiliation(s)
- Masahiro Higashi
- Department of Chemistry, Biology and Marine Science, University of the Ryukyus , 1 Senbaru, Nishihara, Okinawa 903-0213, Japan
| | | | | | | |
Collapse
|
45
|
Affiliation(s)
- Jugal Kumawat
- National Chemical Laboratory (NCL), Dr. Homi Bhabha Road, Pashan, Pune 411008, Maharashtra, India
| | - Virendra Kumar Gupta
- Reliance Technology Group-Hazira, Reliance Industries Ltd., Hazira Complex, Surat, Gujarat 394 510, India
| | - Kumar Vanka
- National Chemical Laboratory (NCL), Dr. Homi Bhabha Road, Pashan, Pune 411008, Maharashtra, India
| |
Collapse
|
46
|
Jurinovich S, Curutchet C, Mennucci B. The Fenna-Matthews-Olson protein revisited: a fully polarizable (TD)DFT/MM description. Chemphyschem 2014; 15:3194-204. [PMID: 25080315 DOI: 10.1002/cphc.201402244] [Citation(s) in RCA: 61] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2014] [Indexed: 01/21/2023]
Abstract
We report a combined molecular dynamics and quantum mechanics (QM)/molecular mechanics (MM) analysis of the excitonic properties of the Fenna-Matthews-Olson (FMO) protein by using a polarizable MM model combined with a time-dependent density functional theory description. Overall, our results indicate that structural fluctuations, electrostatic interactions, and short-range quantum effects can significantly modulate the model Hamiltonian parameters (site energies and couplings). We find that the specific interactions with the axial ligand and the hydrogen-bonded residues are responsible for the energy ladder, with their effects being mainly due to electrostatic interactions, but with short-range quantum contributions that are not negligible. In addition, a striking modulation of the screening effects experienced by the BChl pairs, due to the heterogeneous polarizability of the FMO and solvent environment, is observed. Finally, we find that the exciton model gives a reliable description of the delocalized excited states in the complex.
Collapse
Affiliation(s)
- Sandro Jurinovich
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, Via Risorgimento 35, I-56126 Pisa (Italy)
| | | | | |
Collapse
|