Different Types of Vibrations Interacting with Electronic Excitations in Phycoerythrin 545 and Fenna-Matthews-Olson Antenna Systems

Protein Effects on the Excitation Energies and Exciton Dynamics of the CP24 Antenna Complex

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.

Intermolecular resonance energy transfer between two lutein pigments in light-harvesting complex II studied by frenkel exciton models

The energy transfer pathways in light-harvesting complex II are complicated and the discovery of the energy transfer between the two luteins revealed an unelucidated important role of carotenoids in the energy flow. This energy transfer between the two S2 states of luteins was for the first time investigated using Frenkel exciton models, using a hybrid scheme of molecular mechanics and quantum mechanics. The results show the energy flow between the two luteins under the Förster resonance energy transfer mechanism. The energy transfer caused by energy level resonance occurs in configurations with small energy gaps. This energy transfer pathway is particularly sensitive to conformation. Moreover, according to the statistical characteristics of the data of the energy gaps and coupling values between LUTs, we proposed stochastic exciton Hamiltonian models to facilitate clarification of the energy transfer among pigments in antenna complexes.

Recent progress in atomistic modeling of light-harvesting complexes: a mini review

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.

Spectral densities and absorption spectra of the core antenna complex CP43 from photosystem II

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.

Role of Pigment-Protein Coupling in the Energy Transport Dynamics in the Fenna-Matthews-Olson Complex

The role of pigment-protein coupling in the dynamics of photosynthetic energy transport in chromophoric complexes has not been fully understood. The excitation energy transfer in the photosynthetic system is tremendously efficient. In particular, we investigate the excitation energy transport in the Fenna-Matthews-Olson (FMO) complex. The exciton dynamics and excitation energy transfer (EET) depend on the interaction between the excited chromophores and their environment. Most theoretical models believe that all bacteriochlorophyll-a (BChla) sites are surrounded by the same local protein environment, which is contradicted by the structural analysis of the FMO complex. Based on different values of pigment-protein coupling for different sites, measured in the adiabatic limit, we have theoretically investigated the effect of the heterogeneous local protein environment on the EET process. By the realistic and site-dependent model of the system-bath couplings, the results show that this interaction may have a critical value for the coherent energy-transfer process. Furthermore, we verify that the two transport pathways are coherent and stable to the important parameter reorganization energy of environmental interactions. The quantum dynamical simulations show that the correlation fluctuation keeps the oscillation of the coherent excitation on a long timescale. In addition, due to the inhomogeneous pigment-protein coupling, different BChl sites have asymmetric excitation oscillation timescales.

Time-dependent atomistic simulations of the CP29 light-harvesting complex

Light harvesting as the first step in photosynthesis is of prime importance for life on earth. For a theoretical description of photochemical processes during light harvesting, spectral densities are key quantities. They serve as input functions for modeling the excitation energy transfer dynamics and spectroscopic properties. Herein, a recently developed procedure is applied to determine the spectral densities of the pigments in the minor antenna complex CP29 of photosystem II, which has recently gained attention because of its active role in non-photochemical quenching processes in higher plants. To this end, the density functional-based tight binding (DFTB) method has been employed to enable simulation of the ground state dynamics in a quantum-mechanics/molecular mechanics (QM/MM) scheme for each chlorophyll pigment. Subsequently, the time-dependent extension of the long-range corrected DFTB approach has been used to obtain the excitation energy fluctuations along the ground-state trajectories also in a QM/MM setting. From these results, the spectral densities have been determined and compared for different force fields and to spectral densities from other light-harvesting complexes. In addition, time-dependent and time-independent excitonic Hamiltonians of the system have been constructed and applied to the determination of absorption spectra as well as exciton dynamics.

Theoretical Study of the Spectral Differences of the Fenna-Matthews-Olson Protein from Different Species and Their Mutants

The structural basis for the spectral differences between the Fenna-Matthews-Olson (FMO) proteins from Chlorobaculum tepidum (C. tepidum) and Prosthecochloris aestuarii 2K (P. aestuarii) is yet to be fully understood. Mutation-induced perturbation to the exciton structure and the optical spectra of the complex provide a suitable means to investigate the critical role played by the protein scaffold. In this work, we have performed quantum-mechanics/molecular-mechanics calculations over the molecular dynamics simulation trajectories with the polarized protein-specific charge scheme for both wild-type FMOs and two mutants. Our result reveals that a single-point mutation in the vicinity of BChl 6, namely, W183F of C. tepidum, significantly affects the absorption spectrum, resulting in a switch of the absorption spectrum from type 2, for which the 806 nm band is more pronounced than the 815 nm band, to type 1, for which the 815 nm band is pronounced. Our observations agree with the single-point mutation experiments reported by Saer et al. (Biochim. Biophys. Acta, Bioenerg. 2017, 1858, 288-296) and Khmelnitskiy et al. (J. Phys. Chem. Lett. 2018, 9, 3378-3386). In contrast, the absorption spectrum of the P. aestuarii experiences the opposite transition (from type 1 to type 2) upon the same mutation. Furthermore, by comparing the contributions of individual pigments to the spectra in the wild type and its mutant, we find that a single-point mutation near BChl 6 not only induces changes in excitation energy of BChl 6 per se but also affects the excitonic structures of the neighboring BChls 5 and 7 through strong interpigment electronic couplings, resulting in a significant change in the absorption spectra.

DFTB/MM Molecular Dynamics Simulations of the FMO Light-Harvesting Complex

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.

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

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

Designing and understanding light-harvesting devices with machine learning

Understanding the fundamental processes of light-harvesting is crucial to the development of clean energy materials and devices. Biological organisms have evolved complex metabolic mechanisms to efficiently convert sunlight into chemical energy. Unraveling the secrets of this conversion has inspired the design of clean energy technologies, including solar cells and photocatalytic water splitting. Describing the emergence of macroscopic properties from microscopic processes poses the challenge to bridge length and time scales of several orders of magnitude. Machine learning experiences increased popularity as a tool to bridge the gap between multi-level theoretical models and Edisonian trial-and-error approaches. Machine learning offers opportunities to gain detailed scientific insights into the underlying principles governing light-harvesting phenomena and can accelerate the fabrication of light-harvesting devices.

Determination of the protonation preferences of bilin pigments in cryptophyte antenna complexes

The light-harvesting mechanisms of cryptophyte antenna complexes have attracted considerable attention due to their ability to exhibit maximal photosynthetic activity under very low-light conditions and to display several colors, as well as the observation of vibronic coherent features in their two-dimensional electronic spectra. However, detailed investigations on the interplay between the protein environment and their light-harvesting properties are hampered by the uncertainty related to the protonation state of the underlying bilin pigments. Here we study the protonation preferences of four types of bilin pigments including 15,16-dihydrobiliverdin (DBV), phycoerythrobilin (PEB), phycocyanobilin (PCB) and mesobiliverdin (MBV), which are found in phycoerythrin PE545 and phycocyanin PC577, PC612, PC630 and PC645 complexes. We apply quantum chemical calculations coupled to continuum solvation calculations to predict the intrinsic acidity of bilins in aqueous solution, and then combine molecular dynamics simulations with empirical pKa estimates to investigate the impact of the local protein environment on the acidity of the pigments. We also report measurements of the absorption spectra of the five complexes in a wide range of pH in order to validate our simulations and investigate possible changes in the light harvesting properties of the complexes in the range of physiological pH found in the lumen (pH ∼ 5-7). The results suggest a pKa > 7 for DBV and MBV pigments in the α polypeptide chains of PE545 and PC630/PC645 complexes, which are not coordinated to a negatively charged amino acid. For the other PEB, DBV and PCB pigments, which interact with a Glu or Asp side chain, higher pKa values (pKa > 8) are estimated. Overall, the results support a preferential population of the fully protonated state for bilins in cryptophyte complexes under physiological conditions regardless of the specific type of pigment and local protein environment.

Vibration-Assisted and Vibration-Hampered Excitonic Quantum Transport

The interplay between excitons and vibrations is considered to be a key factor in determining the exciton-transfer properties in light-harvesting complexes. Here we study this interplay theoretically in a model for exciton transport, composed of two chromophores coupled to an exciton source and sink in the presence of vibrations. We consider two cases that show qualitatively distinct transport features. In the first, the vibrations are global and affect the two chromophores simultaneously. In the second case, the vibrations are localized on each chromophore. For global vibrations, the current exhibits antiresonances as a function of the chromophore energy difference, which are due to exciton-polaron interference. For local vibrations, on the contrary, the currents show tunneling resonances at multiples of the vibration energy. Counterintuitively, both effects increase with increasing temperature. Our results demonstrate that an environment can either assist or hamper exciton transport and is in accord with the current understanding of energy transfer in natural exciton-transfer complexes.

Coherent electronic and nuclear dynamics in a rhodamine heterodimer-DNA supramolecular complex

Elucidating the role of quantum coherences in energy migration within biological and artificial multichromophoric antenna systems is the subject of an intense debate. It is also a practical matter because of the decisive implications for understanding the biological processes and engineering artificial materials for solar energy harvesting. A supramolecular rhodamine heterodimer on a DNA scaffold was suitably engineered to mimic the basic donor-acceptor unit of light-harvesting antennas. Ultrafast 2D electronic spectroscopic measurements allowed identifying clear features attributable to a coherent superposition of dimer electronic and vibrational states contributing to the coherent electronic charge beating between the donor and the acceptor. The frequency of electronic charge beating is found to be 970 cm^-1 (34 fs) and can be observed for 150 fs. Through the support of high level ab initio TD-DFT computations of the entire dimer, we established that the vibrational modes preferentially optically accessed do not drive subsequent coupling between the electronic states on the 600 fs of the experiment. It was thereby possible to characterize the time scales of the early time femtosecond dynamics of the electronic coherence built by the optical excitation in a large rigid supramolecular system at a room temperature in solution.

How fine-tuned for energy transfer is the environmental noise produced by proteins around biological chromophores?

Analysis of intermolecular motions of pigment–protein complexes shows no significant difference in influence of local environment despite different biological functions.

Simulation of Quantum Dynamics of Excitonic Systems at Finite Temperature: an efficient method based on Thermo Field Dynamics

Quantum electron-vibrational dynamics in molecular systems at finite temperature is described using an approach based on Thermo Field Dynamics theory. This formulation treats temperature effects in the Hilbert space without introducing the Liouville space. The solution of Thermo Field Dynamics equations with a novel technique for the propagation of Tensor Trains (Matrix Product States) is implemented and discussed. The methodology is applied to the study of the exciton dynamics in the Fenna-Mathews-Olsen complex using a realistic structured spectral density to model the electron-phonon interaction. The results of the simulations highlight the effect of specific vibrational modes on the exciton dynamics and energy transfer process, as well as call for careful modeling of electron-phonon couplings.

Relation between Vibrational Dephasing Time and Energy Gap Fluctuations

Dephasing processes are present in basically all applications in which quantum mechanics plays a role. These applications certainly include excitation energy and charge transfer in biological systems. In a previous study, we have analyzed the vibrational dephasing time as a function of energy gap fluctuation for a large set of molecular simulations. In that investigation, individual molecular subunits were the focus of the calculations. The set of studied molecules included bacteriochlorophylls in Fenna-Matthews-Olson and light-harvesting system 2 complexes as well as bilins in PE545 aggregates. The present work extends this study to entire complexes, including the respective intermolecular couplings. Again, it can be concluded that a universal and inverse proportionality exists between dephasing time and variance of the excitonic energy gap fluctuations, whereas the respective proportionality constants can be rationalized using the energy gap autocorrelation functions. Furthermore, these findings can be extended to the gaps between higher-lying neighboring excitonic states.

Constructing an Interpolated Potential Energy Surface of a Large Molecule: A Case Study with Bacteriochlorophyll a Model in the Fenna-Matthews-Olson Complex

Constructing a reliable potential energy surface (PES) is a key step toward computationally studying the chemical dynamics of any molecular system. The interpolation scheme is a useful tool that can closely follow the accuracy of quantum chemical means at a dramatically reduced computational cost. However, applying interpolation to building a PES of a large molecule is not a straightforward black-box approach, as it frequently encounters practical difficulties associated with its large dimensionality. Here, we present detailed courses of applying interpolation toward building a PES of a large chromophore molecule. We take the example of S₀ and S₁ electronic states of bacteriochlorophyll a (BChla) molecules in the Fenna-Matthews-Olson light harvesting complex. With a reduced model molecule that bears BChla's main π-conjugated ring, various practical approaches are designed for improving the PES quality in a stable manner and for fine-tuning the final surface such that the surface can be adopted for long time molecular dynamics simulations. Combined with parallel implementation, we show that interpolated mechanics/molecular mechanics (IM/MM) simulations of the entire complex in the nanosecond time scale can be conducted readily without any practical issues. With 1500 interpolation data points for each chromophore unit, the PES error relative to the reference quantum chemical calculation is found to be ∼0.15 eV in the thermally accessible region of the conformational space, together with ∼0.01 eV error in S₀ - S₁ transition energies. The performance issue related to the use of a large interpolation database within the framework of our parallel routines is also discussed.

Modeling Electronic-Nuclear Interactions for Excitation Energy Transfer Processes in Light-Harvesting Complexes

An accurate approach for computing intermolecular and intrachromophore contributions to spectral densities to describe the electronic-nuclear interactions relevant for modeling excitation energy transfer processes in light harvesting systems is presented. The approach is based on molecular dynamics (MD) calculations of classical correlation functions of long-range contributions to excitation energy fluctuations and a separate harmonic analysis and single-point gradient quantum calculations for electron-intrachromophore vibrational couplings. A simple model is also presented that enables detailed analysis of the shortcomings of standard MD-based excitation energy fluctuation correlation function approaches. The method introduced here avoids these problems, and its reliability is demonstrated in accurate predictions for bacteriochlorophyll molecules in the Fenna-Matthews-Olson pigment-protein complex, where excellent agreement with experimental spectral densities is found. This efficient approach can provide instantaneous spectral densities for treating the influence of fluctuations in environmental dissipation on fast electronic relaxation.

Relation between Dephasing Time and Energy Gap Fluctuations in Biomolecular Systems

Excitation energy and charge transfer are fundamental processes in biological systems. Because of their quantum nature, the effect of dephasing on these processes is of interest especially when trying to understand their efficiency. Moreover, recent experiments have shown quantum coherences in such systems. As a first step toward a better understanding, we studied the relationship between dephasing time and energy gap fluctuations of the individual molecular subunits. A larger set of molecular simulations has been investigated to shed light on this dependence. This set includes bacterio-chlorophylls in Fenna-Matthews-Olson complexes, the PE545 aggregate, the LH2 complexes, DNA, photolyase, and cryptochromes. For the individual molecular subunits of these aggregates it has been confirmed quantitatively that an inverse proportionality exists between dephasing time and average gap energy fluctuation. However, for entire complexes including the respective intermolecular couplings, such a relation still needs to be verified.

Quantum Chemical Studies of Light Harvesting

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.

Chebyshev Expansion Applied to Dissipative Quantum Systems

To determine the dynamics of a molecular aggregate under the influence of a strongly time-dependent perturbation within a dissipative environment is still, in general, a challenge. The time-dependent perturbation might be, for example, due to external fields or explicitly treated fluctuations within the environment. Methods to calculate the dynamics in these cases do exist though some of these approaches assume that the corresponding correlation functions can be written as a weighted sum of exponentials. One such theory is the hierarchical equations of motion approach. If the environment, however, is described by a complex spectral density or if its temperature is low, these approaches become very inefficient. Therefore, we propose a scheme based on a Chebyshev decomposition of the bath correlation functions and detail the respective quantum master equations within second-order perturbation theory in the environmental coupling. Similar approaches have recently been proposed for systems coupled to Fermionic reservoirs. The proposed scheme is tested for a simple two-level system and compared to existing results. Furthermore, the advantages and disadvantages of the present Chebyshev approach are discussed.

Open quantum system parameters for light harvesting complexes from molecular dynamics

We elucidate the difference between various parameter extraction methods and demonstrate sensitivity to molecular dynamics equilibration.

Developing accurate molecular mechanics force fields for conjugated molecular systems

A rapid method to parameterize the intramolecular component of classical force fields is proposed and applied to a molecular semiconductor, oligomers of conjugated polymers and a biological chromophore.

Limits and potentials of quantum chemical methods in modelling photosynthetic antennae

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.