1
|
Saraceno P, Sardar S, Caferri R, Camargo FVA, Dall'Osto L, D'Andrea C, Bassi R, Cupellini L, Cerullo G, Mennucci B. Probing the Effect of Mutations on Light Harvesting in CP29 by Transient Absorption and First-Principles Simulations. J Phys Chem Lett 2024; 15:6398-6408. [PMID: 38861672 DOI: 10.1021/acs.jpclett.4c01040] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2024]
Abstract
Natural light harvesting is exceptionally efficient thanks to the local energy funnel created within light-harvesting complexes (LHCs). To understand the design principles underlying energy transport in LHCs, ultrafast spectroscopy is often complemented by mutational studies that introduce perturbations into the excitonic structure of the natural complexes. However, such studies may fall short of identifying all excitation energy transfer (EET) pathways and their changes upon mutation. Here, we show that a synergistic combination of first-principles calculations and ultrafast spectroscopy can give unprecedented insight into the EET pathways occurring within LHCs. We measured the transient absorption spectra of the minor CP29 complex of plants and of two mutants, systematically mapping the kinetic components seen in experiments to the simulated exciton dynamics. With our combined strategy, we show that EET in CP29 is surprisingly robust to the changes in the exciton states induced by mutations, explaining the versatility of plant LHCs.
Collapse
Affiliation(s)
- Piermarco Saraceno
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy
| | - Samim Sardar
- Center for Nano Science and Technology @PoliMi, Istituto Italiano di Tecnologia, Via Rubattino 81, 20134 Milan, Italy
| | - Roberto Caferri
- Dipartimento di Biotecnologie, Università di Verona, Strada Le Grazie 15, 37134 Verona, Italy
| | - Franco V A Camargo
- Istituto di Fotonica e Nanotecnologie, Consiglio Nazionale delle Ricerche, Piazza L. da Vinci 32, 20133 Milano, Italy
| | - Luca Dall'Osto
- Dipartimento di Biotecnologie, Università di Verona, Strada Le Grazie 15, 37134 Verona, Italy
| | - Cosimo D'Andrea
- Center for Nano Science and Technology @PoliMi, Istituto Italiano di Tecnologia, Via Rubattino 81, 20134 Milan, Italy
- Dipartimento di Fisica, Politecnico di Milano, Piazza L. da Vinci 32, 20133 Milano, Italy
| | - Roberto Bassi
- Dipartimento di Biotecnologie, Università di Verona, Strada Le Grazie 15, 37134 Verona, Italy
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy
| | - Giulio Cerullo
- Istituto di Fotonica e Nanotecnologie, Consiglio Nazionale delle Ricerche, Piazza L. da Vinci 32, 20133 Milano, Italy
- Dipartimento di Fisica, Politecnico di Milano, Piazza L. da Vinci 32, 20133 Milano, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, 56124 Pisa, Italy
| |
Collapse
|
2
|
Betti E, Saraceno P, Cignoni E, Cupellini L, Mennucci B. Insights into Energy Transfer in Light-Harvesting Complex II Through Machine-Learning Assisted Simulations. J Phys Chem B 2024; 128:5188-5200. [PMID: 38761151 DOI: 10.1021/acs.jpcb.4c01494] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/20/2024]
Abstract
Light-harvesting complex II (LHCII) is the major antenna of higher plants. Energy transfer processes taking place inside its aggregate of chlorophylls have been experimentally investigated with time-resolved techniques, but a complete understanding of the most relevant energy transfer pathways and relative characteristic times remains elusive. Theoretical models to disentangle experimental data in LHCII have long been challenged by the large size and complex nature of the system. Here, we show that a fully first-principles approach combining molecular dynamics and machine learning can be successfully used to reproduce transient absorption spectra and characterize the EET pathways and the involved times.
Collapse
Affiliation(s)
- Elena Betti
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Piermarco Saraceno
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Edoardo Cignoni
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Benedetta Mennucci
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| |
Collapse
|
3
|
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
|
4
|
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] [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
|
5
|
Saraceno P, Sláma V, Cupellini L. First-principles simulation of excitation energy transfer and transient absorption spectroscopy in the CP29 light-harvesting complex. J Chem Phys 2023; 159:184112. [PMID: 37962444 DOI: 10.1063/5.0170295] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2023] [Accepted: 10/24/2023] [Indexed: 11/15/2023] Open
Abstract
The dynamics of delocalized excitons in light-harvesting complexes (LHCs) can be investigated using different experimental techniques, and transient absorption (TA) spectroscopy is one of the most valuable methods for this purpose. A careful interpretation of TA spectra is essential for the clarification of excitation energy transfer (EET) processes occurring during light-harvesting. However, even in the simplest LHCs, a physical model is needed to interpret transient spectra as the number of EET processes occurring at the same time is very large to be disentangled from measurements alone. Physical EET models are commonly built by fittings of the microscopic exciton Hamiltonians and exciton-vibrational parameters, an approach that can lead to biases. Here, we present a first-principles strategy to simulate EET and transient absorption spectra in LHCs, combining molecular dynamics and accurate multiscale quantum chemical calculations to obtain an independent estimate of the excitonic structure of the complex. The microscopic parameters thus obtained are then used in EET simulations to obtain the population dynamics and the related spectroscopic signature. We apply this approach to the CP29 minor antenna complex of plants for which we follow the EET dynamics and transient spectra after excitation in the chlorophyll b region. Our calculations reproduce all the main features observed in the transient absorption spectra and provide independent insight on the excited-state dynamics of CP29. The approach presented here lays the groundwork for the accurate simulation of EET and unbiased interpretation of transient spectra in multichromophoric systems.
Collapse
Affiliation(s)
- Piermarco Saraceno
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, Via G. Moruzzi 13, 56124 Pisa, Italy
| | - Vladislav Sláma
- 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
| |
Collapse
|
6
|
Ozaydin B, Curutchet C. Unraveling the role of thermal fluctuations on the exciton structure of the cryptophyte PC612 and PC645 photosynthetic antenna complexes. Front Mol Biosci 2023; 10:1268278. [PMID: 37790875 PMCID: PMC10544999 DOI: 10.3389/fmolb.2023.1268278] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Accepted: 08/30/2023] [Indexed: 10/05/2023] Open
Abstract
Protein scaffolds play a crucial role in tuning the light harvesting properties of photosynthetic pigment-protein complexes, influencing pigment-protein and pigment-pigment excitonic interactions. Here, we investigate the influence of thermal dynamic effects on the protein tuning mechanisms of phycocyanin PC645 and PC612 antenna complexes of cryptophyte algae, featuring closed or open quaternary structures. We employ a dual molecular dynamics (MD) strategy that combines extensive classical MD simulations with multiple short Born-Oppenheimer quantum/molecular mechanical (QM/MM) simulations to accurately account for both static and dynamic disorder effects. Additionally, we compare the results with an alternative protocol based on multiple QM/MM geometry optimizations of the pigments. Subsequently, we employ polarizable QM/MM calculations using time-dependent density functional theory (TD-DFT) to compute the excited states, and we adopt the full cumulant expansion (FCE) formalism to describe the absorption and circular dichroism spectra. Our findings indicate that thermal effects have only minor impacts on the energy ladder in PC612, despite its remarkable flexibility owing to an open quaternary structure. In striking contrast, thermal effects significantly influence the properties of PC645 due to the absence of a hydrogen bond controlling the twist of ring D in PCB β82 bilins, as well as the larger impact of fluctuations on the excited states of MBV pigments, which possess a higher conjugation length compared to other bilin types. Overall, the dual MD protocol combined with the FCE formalism yields excellent spectral properties for PC612 and PC645, and the resultant excitonic Hamiltonians pave the way for future investigations concerning the implications of open and closed quaternary structures on phycocyanin light harvesting properties.
Collapse
Affiliation(s)
- Beste Ozaydin
- Departament de Farmàcia i Tecnologia Farmacèutica, i Fisicoquímica, Facultat de Farmàcia i Ciències de l’Alimentació, Universitat de Barcelona (UB), Barcelona, Spain
- Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona (UB), Barcelona, Spain
| | - Carles Curutchet
- Departament de Farmàcia i Tecnologia Farmacèutica, i Fisicoquímica, Facultat de Farmàcia i Ciències de l’Alimentació, Universitat de Barcelona (UB), Barcelona, Spain
- Institut de Química Teòrica i Computacional (IQTCUB), Universitat de Barcelona (UB), Barcelona, Spain
| |
Collapse
|
7
|
Cupellini L, Qian P, Nguyen-Phan TC, Gardiner AT, Cogdell RJ. Quantum chemical elucidation of a sevenfold symmetric bacterial antenna complex. PHOTOSYNTHESIS RESEARCH 2023; 156:75-87. [PMID: 35672557 PMCID: PMC10070313 DOI: 10.1007/s11120-022-00925-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/09/2022] [Accepted: 05/12/2022] [Indexed: 06/15/2023]
Abstract
The light-harvesting complex 2 (LH2) of purple bacteria is one of the most studied photosynthetic antenna complexes. Its symmetric structure and ring-like bacteriochlorophyll arrangement make it an ideal system for theoreticians and spectroscopists. LH2 complexes from most bacterial species are thought to have eightfold or ninefold symmetry, but recently a sevenfold symmetric LH2 structure from the bacterium Mch. purpuratum was solved by Cryo-Electron microscopy. This LH2 also possesses unique near-infrared absorption and circular dichroism (CD) spectral properties. Here we use an atomistic strategy to elucidate the spectral properties of Mch. purpuratum LH2 and understand the differences with the most commonly studied LH2 from Rbl. acidophilus. Our strategy exploits a combination of molecular dynamics simulations, multiscale polarizable quantum mechanics/molecular mechanics calculations, and lineshape simulations. Our calculations reveal that the spectral properties of LH2 complexes are tuned by site energies and exciton couplings, which in turn depend on the structural fluctuations of the bacteriochlorophylls. Our strategy proves effective in reproducing the absorption and CD spectra of the two LH2 complexes, and in uncovering the origin of their differences. This work proves that it is possible to obtain insight into the spectral tuning strategies of purple bacteria by quantitatively simulating the spectral properties of their antenna complexes.
Collapse
Affiliation(s)
- Lorenzo Cupellini
- Department of Chemistry and Industrial Chemistry, University of Pisa, 56124, Pisa, Italy.
| | - Pu Qian
- Materials and Structure Analysis, Thermofisher Scientific, Achtseweg Nordic 5, 5651 GTC, Eindhoven, The Netherlands
| | - Tu C Nguyen-Phan
- Institute of Molecular, Cell and Systems Biology, University of Glasgow, Glasgow, G12 8QQ, UK
| | - Alastair T Gardiner
- Laboratory of Anoxygenic Phototrophs, Centre Algatech, Novohradská 237 - Opatovický mlýn, 379 01, Třeboň, Czech Republic
| | - Richard J Cogdell
- Institute of Molecular, Cell and Systems Biology, University of Glasgow, Glasgow, G12 8QQ, UK
| |
Collapse
|
8
|
Hu YY, Liu XL, Yao HD, Jiang YL, Li K, Chen MQ, Wang P, Zhang JP. PEG effects on excitonic properties of LH2 from Rhodobacter sphaeroides 2.4.1 in different environments. Chem Phys Lett 2023. [DOI: 10.1016/j.cplett.2023.140477] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/01/2023]
|
9
|
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
|
10
|
Nöthling JA, Mancal T, Kruger T. Accuracy of approximate methods for the calculation of absorption-type linear spectra with a complex system-bath coupling. J Chem Phys 2022; 157:095103. [DOI: 10.1063/5.0100977] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The accuracy of approximate methods for calculating linear optical spectra depends on many variables. In this study, we fix most of these parameters to typical values found in photosynthetic light-harvesting complexes of plants and determine the accuracy of approximate spectra with respect to exact calculation as a function of the energy gap and interpigment coupling in a pigment dimer. We use a spectral density with the first eight intramolecular modes of chlorophyll a and include inhomogeneous disorder for the calculation of spectra. We compare the accuracy of absorption, linear dichroism, and circular dichroism spectra calculated using the Full Cumulant Expansion (FCE), coherent time-dependent Redfield (ctR), and time-independent Redfield and modified Redfield methods. As a reference we use spectra calculated with the Exact Stochastic Path Integral Evaluation method. We find the FCE method to be the most accurate for the calculation of all spectra. The ctR method performs well for the qualitative calculation of absorption and linear dichroism spectra when pigments are moderately coupled (∼15 cm-1), but ctR spectra may differ significantly from exact spectra when strong interpigment coupling (>100 cm-1) is present. The dependence of the quality of Redfield and modified Redfield spectra on molecular parameters is similar, and these methods almost always perform worse than ctR, especially when the interpigment coupling is strong or the excitonic energy gap is small (for a given coupling). The accuracy of approximate spectra is not affected by resonance with intramolecular modes for typical system-bath coupling and disorder values found in plant light-harvesting complexes.
Collapse
Affiliation(s)
| | - Tomas Mancal
- Faculty of Mathematics and Physics, Charles University Faculty of Mathematics and Physics, Czech Republic
| | | |
Collapse
|
11
|
Chatterjee S, Nam Y, Salimi A, Lee JY. Monitoring early-stage β-amyloid dimer aggregation by histidine site-specific two-dimensional infrared spectroscopy in a simulation study. Phys Chem Chem Phys 2022; 24:18691-18702. [PMID: 35899740 DOI: 10.1039/d2cp02479a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Monitoring early-stage β-amyloid (Aβ) dimerization is a formidable challenge for understanding neurological diseases. We compared β-sheet formation and histidine site-specific two-dimensional infrared (2D IR) spectroscopic signatures of Aβ dimers with different histidine states (δ; Nδ1-H, ε; Nε2-H, or π; both protonated). Molecular dynamics (MD) simulations revealed that β-sheet formation is favored for the δδδ:δδδ and πππ:πππ tautomeric isomers showing strong couplings and frequent contacts between the central hydrophobic core and C-terminus compared with the εεε:εεε isomer. Characteristic blue-shifts in the 2D IR central bands were observed upon monomer-dimer transformation. The εεε:εεε dimer exhibited larger frequency shifts than δδδ:δδδ and πππ:πππ implying that the red-shift may have a correlation with Nδ1-H(δ) protonation. Our results support the tautomerization/protonation hypothesis that attributes Aβ misfolding to histidine tautomers as a possible primary initiator for Aβ aggregation and facilitates the application of histidine site-specific 2D IR spectroscopy for studying early-stage Aβ self-assembly.
Collapse
Affiliation(s)
| | - Yeonsig Nam
- Department of Chemistry, Sungkyunkwan University, Suwon 440-746, Korea. .,Department of Chemistry, University of California, Irvine, California 92697-2025, USA
| | - Abbas Salimi
- Department of Chemistry, Sungkyunkwan University, Suwon 440-746, Korea.
| | - Jin Yong Lee
- Department of Chemistry, Sungkyunkwan University, Suwon 440-746, Korea.
| |
Collapse
|
12
|
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
|
13
|
Cignoni E, Slama V, Cupellini L, Mennucci B. The atomistic modeling of light-harvesting complexes from the physical models to the computational protocol. J Chem Phys 2022; 156:120901. [DOI: 10.1063/5.0086275] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
The function of light-harvesting complexes is determined by a complex network of dynamic interactions among all the different components: the aggregate of pigments, the protein, and the surrounding environment. Complete and reliable predictions on these types of composite systems can be only achieved with an atomistic description. In the last few decades, there have been important advances in the atomistic modeling of light-harvesting complexes. These advances have involved both the completeness of the physical models and the accuracy and effectiveness of the computational protocols. In this Perspective, we present an overview of the main theoretical and computational breakthroughs attained so far in the field, with particular focus on the important role played by the protein and its dynamics. We then discuss the open problems in their accurate modeling that still need to be addressed. To illustrate an effective computational workflow for the modeling of light harvesting complexes, we take as an example the plant antenna complex CP29 and its H111N mutant.
Collapse
Affiliation(s)
- Edoardo Cignoni
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Vladislav Slama
- 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
|
14
|
Braver Y, Valkunas L, Gelzinis A. Quantum-Classical Approach for Calculations of Absorption and Fluorescence: Principles and Applications. J Chem Theory Comput 2021; 17:7157-7168. [PMID: 34618457 PMCID: PMC8719324 DOI: 10.1021/acs.jctc.1c00777] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2021] [Indexed: 01/20/2023]
Abstract
Absorption and fluorescence spectroscopy techniques provide a wealth of information on molecular systems. The simulations of such experiments remain challenging, however, despite the efforts put into developing the underlying theory. An attractive method of simulating the behavior of molecular systems is provided by the quantum-classical theory─it enables one to keep track of the state of the bath explicitly, which is needed for accurate calculations of fluorescence spectra. Unfortunately, until now there have been relatively few works that apply quantum-classical methods for modeling spectroscopic data. In this work, we seek to provide a framework for the calculations of absorption and fluorescence lineshapes of molecular systems using the methods based on the quantum-classical Liouville equation, namely, the forward-backward trajectory solution (FBTS) and the non-Hamiltonian variant of the Poisson bracket mapping equation (PBME-nH). We perform calculations on a molecular dimer and the photosynthetic Fenna-Matthews-Olson complex. We find that in the case of absorption, the FBTS outperforms PBME-nH, consistently yielding highly accurate results. We next demonstrate that for fluorescence calculations, the method of choice is a hybrid approach, which we call PBME-nH-Jeff, that utilizes the effective coupling theory [Gelzinis, A.; J. Chem. Phys. 2020, 152, 051103] to estimate the excited state equilibrium density operator. Thus, we find that FBTS and PBME-nH-Jeff are excellent candidates for simulating, respectively, absorption and fluorescence spectra of real molecular systems.
Collapse
Affiliation(s)
- Yakov Braver
- Institute
of Chemical Physics, Faculty of Physics, Vilnius University, Saulėtekio Avenue 9-III, LT-10222 Vilnius, Lithuania
- Department
of Molecular Compound Physics, Center for
Physical Sciences and Technology, Saulėtekio Avenue 3, LT-10257 Vilnius, Lithuania
| | - Leonas Valkunas
- Institute
of Chemical Physics, Faculty of Physics, Vilnius University, Saulėtekio Avenue 9-III, LT-10222 Vilnius, Lithuania
- Department
of Molecular Compound Physics, Center for
Physical Sciences and Technology, Saulėtekio Avenue 3, LT-10257 Vilnius, Lithuania
| | - Andrius Gelzinis
- Institute
of Chemical Physics, Faculty of Physics, Vilnius University, Saulėtekio Avenue 9-III, LT-10222 Vilnius, Lithuania
- Department
of Molecular Compound Physics, Center for
Physical Sciences and Technology, Saulėtekio Avenue 3, LT-10257 Vilnius, Lithuania
| |
Collapse
|
15
|
Green JA, Asha H, Santoro F, Improta R. Excitonic Model for Strongly Coupled Multichromophoric Systems: The Electronic Circular Dichroism Spectra of Guanine Quadruplexes as Test Cases. J Chem Theory Comput 2021; 17:405-415. [PMID: 33378185 DOI: 10.1021/acs.jctc.0c01100] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
We here propose a general and flexible approach, based on fragment diabatization, which incorporates charge transfer states and significantly increases the reliability of excitonic Hamiltonians for systems where the chromophores are very close. This model (FrDEx) is used to compute the electronic circular dichroism and absorption spectra of two prototype guanine-rich DNA sequences folded in quadruple helices (GQs), i.e., a fragment of the human telomeric sequence (Tel21, antiparallel), and (TGGGGT)4 (TG4T, parallel). Calculations on different subsets of Tel21 and TG4T, from dimers to tetramers, show that FrDEx provides spectra close to the reference full quantum mechanical (QM) ones (obtained with time-dependent density functional theory), with significant improvements with respect to "standard" excitonic Hamiltonians. Furthermore, these tests enable the most cost-effective procedure for the whole GQ to be determined. FrDEx spectra of Tel21 and TG4T are also in good agreement with the QM and experimental ones and give access to interesting insights into the chemical-physical effects modulating the spectral signals. FrDEx could be profitably used to investigate many other biological and nanotechnological materials, from DNA to (opto)electronic polymers.
Collapse
Affiliation(s)
- James A Green
- Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), via Mezzocannone 16, I-80136 Napoli, Italy
| | - Haritha Asha
- Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), via Mezzocannone 16, I-80136 Napoli, Italy
| | - Fabrizio Santoro
- Consiglio Nazionale delle Ricerche, Istituto di Chimica dei Composti Organo Metallici (ICCOM-CNR), SS di Pisa, Area della Ricerca, via G. Moruzzi 1, I-56124 Pisa, Italy
| | - Roberto Improta
- Consiglio Nazionale delle Ricerche, Istituto di Biostrutture e Bioimmagini (IBB-CNR), via Mezzocannone 16, I-80136 Napoli, Italy
| |
Collapse
|
16
|
Cupellini L, Lipparini F, Cao J. Absorption and Circular Dichroism Spectra of Molecular Aggregates With the Full Cumulant Expansion. J Phys Chem B 2020; 124:8610-8617. [PMID: 32901476 PMCID: PMC7901647 DOI: 10.1021/acs.jpcb.0c05180] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
The exciton Hamiltonian of multichromophoric aggregates can be probed by spectroscopic
techniques such as linear absorption and circular dichroism. To compare calculated
Hamiltonians to experiments, a lineshape theory is needed, which takes into account the
coupling of the excitons with inter- and intramolecular vibrations. This coupling is
normally introduced in a perturbative way through the cumulant expansion formalism and
further approximated by assuming a Markovian exciton dynamics, for example with the
modified Redfield theory. Here, we present the implementation of the full cumulant
expansion (FCE) formalism (142, 2015, 09410625747060) to
efficiently compute absorption and circular dichroism spectra of molecular aggregates
beyond the Markov approximation, without restrictions on the form of
exciton–phonon coupling. By employing the LH2 system of purple bacteria as a
challenging test case, we compare the FCE lineshapes with the Markovian lineshapes
obtained with the modified Redfield theory, showing that the latter presents a less
satisfying agreement with experiments. The FCE approach instead accurately describes the
lineshapes, especially in the vibronic sideband of the B800 peak. We envision that the
FCE approach will become a valuable tool for accurately comparing model exciton
Hamiltonians with optical spectroscopy experiments.
Collapse
Affiliation(s)
- Lorenzo Cupellini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, University of Pisa, via G. Moruzzi 13, 56124 Pisa, Italy
| | - Jianshu Cao
- Department of Chemistry, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, United States
| |
Collapse
|