The Eighth Bacteriochlorophyll Completes the Excitation Energy Funnel in the FMO Protein

Non-adiabatic molecular dynamics simulations provide new insights into the exciton transfer in the Fenna-Matthews-Olson complex

A trajectory surface hopping approach, which uses machine learning to speed up the most time-consuming steps, has been adopted to investigate the exciton transfer in light-harvesting systems. The present neural networks achieve high accuracy in predicting both Coulomb couplings and excitation energies. The latter are predicted taking into account the environment of the pigments. Direct simulation of exciton dynamics through light-harvesting complexes on significant time scales is usually challenging due to the coupled motion of nuclear and electronic degrees of freedom in these rather large systems containing several relatively large pigments. In the present approach, however, we are able to evaluate a statistically significant number of non-adiabatic molecular dynamics trajectories with respect to exciton delocalization and exciton paths. The formalism is applied to the Fenna-Matthews-Olson complex of green sulfur bacteria, which transfers energy from the light-harvesting chlorosome to the reaction center with astonishing efficiency. The system has been studied experimentally and theoretically for decades. In total, we were able to simulate non-adiabatically more than 30 ns, sampling also the relevant space of parameters within their uncertainty. Our simulations show that the driving force supplied by the energy landscape resulting from electrostatic tuning is sufficient to funnel the energy towards site 3, from where it can be transferred to the reaction center. However, the high efficiency of transfer within a picosecond timescale can be attributed to the rather unusual properties of the BChl a molecules, resulting in very low inner and outer-sphere reorganization energies, not matched by any other organic molecule, e.g., used in organic electronics. A comparison with electron and exciton transfer in organic materials is particularly illuminating, suggesting a mechanism to explain the comparably high transfer efficiency.

Exciton dynamics from the mapping approach to surface hopping: comparison with Förster and Redfield theories

We compare the recently introduced multi-state mapping approach to surface hopping (MASH) with the Förster and Redfield theories of excitation energy transfer. Whereas Förster theory relies on weak coupling between chromophores, and Redfield theory assumes the electronic excitations to be weakly coupled to fast chromophore vibrations, MASH is free from any perturbative or Markovian approximations. We illustrate this with an example application to the rate of energy transfer in a Frenkel-exciton dimer, showing that MASH interpolates correctly between the opposing regimes in which the Förster and Redfield results are reliable. We then compare the three methods for a realistic model of the Fenna-Matthews-Olson complex with a structured vibrational spectral density and static disorder in the excitation energies. In this case there are no exact results for comparison so we use MASH to assess the validity of Förster and Redfield theories. We find that Förster theory is the more accurate of the two on the picosecond timescale, as has been shown previously for a simpler model of this particular light-harvesting complex. We also explore various ways to sample the initial electronic state in MASH and find that they all give very similar results for exciton dynamics.

Incoherent Energy Transfer between the Baseplate and FMO Protein Explored at Ideal Geometries

The Förster resonance energy transfer (FRET) between the Fenna-Matthews-Olson (FMO) protein complex and the chlorosomal baseplate (CBP) is investigated by using an idealized model. This simplified model is based on crystal structure and molecular dynamics conformations. Some of the further input, such as the transition dipole moments, was extracted from earlier molecular-level simulations. The resulting model mimics the effects of the relative position between the CBP and the FMO complex on the corresponding FRET efficiency under ideal conditions, involving about 1.3 billion FRET calculations per investigated model. In this idealized model and employing some approximations, it is found that FRET efficiency is almost completely independent of the FMO trimer orientation (displacement, distance, and rotation), despite FMO and CBP being highly structured complexes. Even removing individual FMO BChl triples will only reduce the FRET efficiency by up to 8.6%. An FMO containing only the least efficient BChl triple will retain about 25% of the FRET efficiency of a full FMO complex. In addition to its proposed function as an energetic funnel, FMO is thus identified to act as a highly robust spatial funnel for CBP excitation harvesting, independent of the mutual CBP-FMO orientation.

Probing exciton dynamics with spectral selectivity through the use of quantum entangled photons

Quantum light is increasingly recognized as a promising resource for developing optical measurement techniques. Particular attention has been paid to enhancing the precision of the measurements beyond classical techniques by using nonclassical correlations between quantum entangled photons. Recent advances in the quantum optics technology have made it possible to manipulate spectral and temporal properties of entangled photons, and photon correlations can facilitate the extraction of matter information with relatively simple optical systems compared to conventional schemes. In these respects, the applications of entangled photons to time-resolved spectroscopy can open new avenues for unambiguously extracting information on dynamical processes in complex molecular and materials systems. Here, we propose time-resolved spectroscopy in which specific signal contributions are selectively enhanced by harnessing nonclassical correlations of entangled photons. The entanglement time characterizes the mutual delay between an entangled twin and determines the spectral distribution of photon correlations. The entanglement time plays a dual role as the knob for controlling the accessible time region of dynamical processes and the degrees of spectral selectivity. In this sense, the role of the entanglement time is substantially equivalent to the temporal width of the classical laser pulse. The results demonstrate that the application of quantum entangled photons to time-resolved spectroscopy leads to monitoring dynamical processes in complex molecular and materials systems by selectively extracting desired signal contributions from congested spectra. We anticipate that more elaborately engineered photon states would broaden the availability of quantum light spectroscopy.

Predicting Mutation-Induced Changes in the Electronic Properties of Photosynthetic Proteins from First Principles: The Fenna-Matthews-Olson Complex Example

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.

Quantum Entanglement and State-Transference in Fenna-Matthews-Olson Complexes: A Post-Experimental Simulation Analysis in the Computational Biology Domain

Fenna-Mathews-Olson complexes participate in the photosynthetic process of Sulfur Green Bacteria. These biological subsystems exhibit quantum features which possibly are responsible for their high efficiency; the latter may comprise multipartite entanglement and the apparent tunnelling of the initial quantum state. At first, to study these aspects, a multidisciplinary approach including experimental biology, spectroscopy, physics, and math modelling is required. Then, a global computer modelling analysis is achieved in the computational biology domain. The current work implements the Hierarchical Equations of Motion to numerically solve the open quantum system problem regarding this complex. The time-evolved states obtained with this method are then analysed under several measures of entanglement, some of them already proposed in the literature. However, for the first time, the maximum overlap with respect to the closest separable state is employed. This authentic multipartite entanglement measure provides information on the correlations, not only based on the system bipartitions as in the usual analysis. Our study has led us to note a different view of FMO multipartite entanglement as tiny contributions to the global entanglement suggested by other more basic measurements. Additionally, in another related trend, the initial state, considered as a Förster Resonance Energy Transfer, is tracked using a novel approach, considering how it could be followed under the fidelity measure on all possible permutations of the FMO subsystems through its dynamical evolution by observing the tunnelling in the most probable locations. Both analyses demanded significant computational work, making for a clear example of the complexity required in computational biology.

Effects of Heterogeneous Protein Environment on Excitation Energy Transfer Dynamics in the Fenna-Matthews-Olson Complex

The Fenna-Matthews-Olson (FMO) complex of green sulfur bacteria has been serving as a prototypical light-harvesting protein for studying excitation energy transfer (EET) dynamics in photosynthesis. The most widely used Frenkel exciton model for FMO complex assumes that each excited bacteriochlorophyll site couples to an identical and isolated harmonic bath, which does not account for the heterogeneous local protein environment. To better describe the realistic environment, we propose to use the recently developed multistate harmonic (MSH) model, which contains a globally shared bath that couples to the different pigment sites according to the atomistic quantum mechanics/molecular mechanics simulations with explicit protein scaffold and solvent. In this work, the effects of heterogeneous protein environment on EET in FMO complexes from Prosthecochloris aestuarii and Chlorobium tepidum, specifically including realistic spectral density, site-dependent reorganization energies, and system-bath couplings are investigated. Semiclassical and mixed quantum-classical mapping dynamics were applied to obtain the nonadiabatic EET dynamics in several models ranging from the Frenkel exciton model to the MSH model and their variants. The MSH model with realistic spectral density and site-dependent system-bath couplings displays slower EET dynamics than the Frenkel exciton model. Our comparative study shows that larger average reorganization energy, heterogeneity in spectral densities, and low-frequency modes could facilitate energy dissipation, which is insensitive to the static disorder in reorganization energies. The effects of the spectral densities and system-bath couplings along with the MSH model can be used to optimize EET dynamics for artificial light-harvesting systems.

One-Shot Trajectory Learning of Open Quantum Systems Dynamics

Nonadiabatic quantum dynamics is important for understanding light-harvesting processes, but its propagation with traditional methods can be rather expensive. Here we present a one-shot trajectory learning approach that allows us to directly make an ultrafast prediction of the entire trajectory of the reduced density matrix for a new set of such simulation parameters as temperature and reorganization energy. The whole 10-ps-long propagation takes 70 ms as we demonstrate on the comparatively large quantum system, the Fenna-Matthews-Olsen (FMO) complex. Our approach also significantly reduces time and memory requirements for training.

Explaining the Efficiency of Photosynthesis: Quantum Uncertainty or Classical Vibrations?

Photosynthetic organisms are known to use a mechanism of vibrationally assisted exciton energy transfer to efficiently harvest energy from light. The importance of quantum effects in this mechanism is a long-standing topic of debate, which has traditionally focused on the role of excitonic coherences. Here, we address another recent claim: that the efficient energy transfer in the Fenna-Matthews-Olson complex relies on nuclear quantum uncertainty and would not function if the vibrations were classical. We present a counter-example to this claim, showing by trajectory-based simulations that a description in terms of quantum electrons and classical nuclei is indeed sufficient to describe the funneling of energy to the reaction center. We analyze and compare these findings to previous classical-nuclear approximations that predicted the absence of an energy funnel and conclude that the key difference and the reason for the discrepancy is the ability of the trajectories to properly account for Newton's third law.

Predicting the future of excitation energy transfer in light-harvesting complex with artificial intelligence-based quantum dynamics

Exploring excitation energy transfer (EET) in light-harvesting complexes (LHCs) is essential for understanding the natural processes and design of highly-efficient photovoltaic devices. LHCs are open systems, where quantum effects may play a crucial role for almost perfect utilization of solar energy. Simulation of energy transfer with inclusion of quantum effects can be done within the framework of dissipative quantum dynamics (QD), which are computationally expensive. Thus, artificial intelligence (AI) offers itself as a tool for reducing the computational cost. Here we suggest AI-QD approach using AI to directly predict QD as a function of time and other parameters such as temperature, reorganization energy, etc., completely circumventing the need of recursive step-wise dynamics propagation in contrast to the traditional QD and alternative, recursive AI-based QD approaches. Our trajectory-learning AI-QD approach is able to predict the correct asymptotic behavior of QD at infinite time. We demonstrate AI-QD on seven-sites Fenna–Matthews–Olson (FMO) complex.

Simulations of energy transfer in light-harvesting complexes are computationally very demanding. Here the authors apply an artificial intelligence quantum dissipative algorithm to study the excited state energy transfer dynamics in a light-harvesting complex.

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.

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.

Semiclassical Modified Redfield and Generalized Förster Theories of Exciton Relaxation/Transfer in Light-Harvesting Complexes: The Quest for the Principle of Detailed Balance

A conceptual problem of transfer theories that use a semiclassical description of the electron-vibrational coupling is the neglect of the correlation between momenta and coordinates of nuclei. In the Redfield theory of exciton relaxation, this neglect leads to a violation of the principle of detailed balance; equal “uphill” and “downhill” transfer rate constants are obtained. Here, we investigate how this result depends on nuclear reorganization effects, neglected in Redfield but taken into account in the modified Redfield theory. These reorganization effects, resulting from a partial localization of excited states, are found to promote a preferential “downhill” relaxation of excitation energy. However, for realistic spectral densities of light-harvesting antennae in photosynthesis, the reorganization effects are too small to compensate for the missing coordinate–momentum uncertainty. For weaker excitonic couplings as they occur between domains of strongly coupled pigments, we find the principle of detailed balance to be fulfilled in a semiclassical variant of the generalized Förster theory. A qualitatively correct description of the transfer is obtained with this theory at a significantly lower computational cost as with the quantum generalized Förster theory. Larger deviations between the two theories are expected for large energy gaps as they occur in complexes with chemically different pigments.

Photosynthetic Light-Harvesting (Antenna) Complexes-Structures and Functions

Chlorophylls and bacteriochlorophylls, together with carotenoids, serve, noncovalently bound to specific apoproteins, as principal light-harvesting and energy-transforming pigments in photosynthetic organisms. In recent years, enormous progress has been achieved in the elucidation of structures and functions of light-harvesting (antenna) complexes, photosynthetic reaction centers and even entire photosystems. It is becoming increasingly clear that light-harvesting complexes not only serve to enlarge the absorption cross sections of the respective reaction centers but are vitally important in short- and long-term adaptation of the photosynthetic apparatus and regulation of the energy-transforming processes in response to external and internal conditions. Thus, the wide variety of structural diversity in photosynthetic antenna “designs” becomes conceivable. It is, however, common for LHCs to form trimeric (or multiples thereof) structures. We propose a simple, tentative explanation of the trimer issue, based on the 2D world created by photosynthetic membrane systems.

Time-Domain Line-Shape Analysis from 2D Spectroscopy to Precisely Determine Hamiltonian Parameters for a Photosynthetic Complex

Optical signals come from coherences between quantum states, with spectral line widths determined by the coherences' dephasing dynamics. Using a 2D electronic spectrometer, we observe weak coherence- and rephasing-time-domain signals persisting to 1 ps in the Fenna-Matthews-Olson complex at 77 K. These are coherences between the ground and excited states prepared after the complex interacts once or three times with light, rather than zero-quantum coherences that are more frequently investigated following two interactions. Here, we use these small but persistent signal components to isolate spectral contributions with narrowed peaks and reveal the system's eigenenergies.

Exciton Lifetime Distributions and Population Dynamics in the FMO Protein Complex from Prosthecochloris aestuarii

Significant protein rearrangement upon excitation and energy transfer in Fenna-Matthews-Olson protein of Prosthecochloris aestuarii results in a modified energy landscape, which induces more changes in pigment site energies than predicted by the "standard" hole-burning theory. The energy changes are elucidated by simulations while investigating the effects of site-dependent disorder, both static (site-energy distribution widths) and dynamic (spectral density shapes). The resulting optimized site energies and their fluctuations are consistent with relative differences observed in inhomogeneous widths calculated by recent molecular dynamic simulations. Two sets of different spectral densities reveal how their shapes affect the population dynamics and distribution of exciton lifetimes. Calculations revealed the wavelength-dependent distributions of exciton lifetimes (T ₁) in the femtosecond to picosecond time frame. We suggest that the calculated multimodal and asymmetric wavelength-dependent T ₁ distributions offer more insight into the interpretation of resonant hole-burned (HB) spectra, kinetic traces in two-dimensional (2D) electronic spectroscopy experiments, and widely used global analyses in fitting data from transient absorption experiments.

Static Disorder in Excitation Energies of the Fenna-Matthews-Olson Protein: Structure-Based Theory Meets Experiment

Inhomogeneous broadening of optical lines of the Fenna-Matthews-Olson (FMO) light-harvesting protein is investigated by combining a Monte Carlo sampling of low-energy conformational substates of the protein with a quantum chemical/electrostatic calculation of local transition energies (site energies) of the pigments. The good agreement between the optical spectra calculated for the inhomogeneous ensemble and the experimental data demonstrates that electrostatics is the dominant contributor to static disorder in site energies. Rotamers of polar amino acid side chains are found to cause bimodal distribution functions of site energy shifts, which can be probed by hole burning and single-molecule spectroscopy. When summing over the large number of contributions, the resulting distribution functions of the site energies become Gaussians, and the correlations in site energy fluctuations at different sites practically average to zero. These results demonstrate that static disorder in the FMO protein is in the realm of the central limit theorem of statistics.

Normal mode analysis of spectral density of FMO trimers: Intra- and intermonomer energy transfer

The intermolecular contribution to the spectral density of the exciton-vibrational coupling of the homotrimeric Fenna-Matthews-Olson (FMO) light-harvesting protein of green sulfur bacteria P. aestuarii is analyzed by combining a normal mode analysis of the protein with the charge density coupling method for the calculation of local transition energies of the pigments. Correlations in site energy fluctuations across the whole FMO trimer are found at low vibrational frequencies. Including, additionally, the high-frequency intrapigment part of the spectral density, extracted from line-narrowing spectra, we study intra- and intermonomer exciton transfer. Whereas the intrapigment part of the spectral density is important for fast intramonomer exciton relaxation, the intermolecular contributions (due to pigment-environment coupling) determine the intermonomer exciton transfer. Neither the variations of the local Huang-Rhys factors nor the correlations in site energy fluctuations have a critical influence on energy transfer. At room temperature, the intermonomer transfer in the FMO protein occurs on a 10 ps time scale, whereas intramonomer exciton equilibration is roughly two orders of magnitude faster. At cryogenic temperatures, intermonomer transfer limits the lifetimes of the lowest exciton band. The lifetimes are found to increase between 20 ps in the center of this band up to 100 ps toward lower energies, which is in very good agreement with the estimates from hole burning data. Interestingly, exciton delocalization in the FMO monomers is found to slow down intermonomer energy transfer, at both physiological and cryogenic temperatures.

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.

Importance of Intermonomer Couplings of the FMO Complex in Coherently Initiated Dynamics

In this work, the hierarchical equation of motion approach was employed to study the excitation transfer efficiency (ETE) from the Fenna-Matthews-Olson (FMO) complex to the reaction center core complex (RCC). Further, the realistic initial pure state (created by the interaction of femtosecond laser pulse with the isolated FMO complex) was used in the presence of an inhomogeneous protein environment. We observed the same amount of excitation transfer from the FMO to the RCC, in agreement with the experiment using femtosecond absorption spectroscopy. Simulation analysis shows weak intermonomer couplings are strongly affecting the dynamics of ETE, which implies that all the monomers of the FMO complex need to be considered as connected with the RCC pigment for a better estimate of the ETE.

Quantum biology revisited

Photosynthesis is a highly optimized process from which valuable lessons can be learned about the operating principles in nature. Its primary steps involve energy transport operating near theoretical quantum limits in efficiency. Recently, extensive research was motivated by the hypothesis that nature used quantum coherences to direct energy transfer. This body of work, a cornerstone for the field of quantum biology, rests on the interpretation of small-amplitude oscillations in two-dimensional electronic spectra of photosynthetic complexes. This Review discusses recent work reexamining these claims and demonstrates that interexciton coherences are too short lived to have any functional significance in photosynthetic energy transfer. Instead, the observed long-lived coherences originate from impulsively excited vibrations, generally observed in femtosecond spectroscopy. These efforts, collectively, lead to a more detailed understanding of the quantum aspects of dissipation. Nature, rather than trying to avoid dissipation, exploits it via engineering of exciton-bath interaction to create efficient energy flow.

Theoretical Study on the Effect of Environment on Excitation Energy Transfer in Photosynthetic Light-Harvesting Systems

In recent years, two-dimensional (2D) electronic spectroscopy experiments prove that the excitation energy transfer (EET) in photosynthetic light-harvesting systems presents long-lived electronic quantum beating signals. After being discovered in the light-harvesting system, the quantum coherence effect has aroused widespread discussion. To illustrate the EET process in the Fenna-Matthews-Olson (FMO) and phycocyanin 645 (PC645) complex, the local protein environment is often thought to be the same; however, this is ambivalent to the practical structural analysis of the light-harvesting complex. By adopting the dissipaton equation of motion theory, we present the effect of a heterogeneous protein environment on the energy transfer process with accurate numerical results. We demonstrate that the energy transfer process relies on the local heterogeneous environment for the FMO complex. A similar good agreement is found for the PC645 complex. Furthermore, we discuss the optimal value of different chromophores in the excitation energy transfer process by controlling the environmental characteristics.

Predictive First-Principles Modeling of a Photosynthetic Antenna Protein: The Fenna-Matthews-Olson Complex

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.

FMOxFMO: Elucidating Excitonic Interactions in the Fenna-Matthews-Olson Complex with the Fragment Molecular Orbital Method

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.

Successes & challenges in the atomistic modeling of light-harvesting and its photoregulation

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.

The influence of quaternary structure on the stability of Fenna-Matthews-Olson (FMO) antenna complexes

The trimeric nature of the Fenna-Matthews-Olson (FMO) protein antenna complex from green sulfur phototrophic bacteria was investigated. Mutations were introduced into the protein at positions 142 and 198, which were chosen to destabilize the intra-trimer salt bridges between adjacent monomers. Strains bearing the mutations R142L, R198L, or their combination, exhibited altered optical absorption spectra of purified membranes and fluoresced more intensely than the wild type. In particular, the introduction of the R142L mutation resulted in slower culture growth rates, as well as an FMO complex that was not able to be isolated in appreciable quantities, while the R198L mutation yielded an FMO complex with increased sensitivity to sodium thiocyanate and Triton X-100 treatments. Native and denaturing PAGE experiments suggest that much of the FMO complexes in the mutant strains pool with the insoluble material upon membrane solubilization with n-dodecyl β-D-maltoside, a mild nonionic detergent. Taken together, our results suggest that the quaternary structure of the FMO complex, the homotrimer, is an important factor in the maintenance of the complex's tertiary structure.

Neutron and X-ray analysis of the Fenna-Matthews-Olson photosynthetic antenna complex from Prosthecochloris aestuarii

The Fenna-Matthews-Olson protein from Prosthecochloris aestuarii (PaFMO) has been crystallized in a new form that is amenable to high-resolution X-ray and neutron analysis. The crystals belonged to space group H3, with unit-cell parameters a = b = 83.64, c = 294.78 Å, and diffracted X-rays to ∼1.7 Å resolution at room temperature. Large PaFMO crystals grown to volumes of 0.3-0.5 mm³ diffracted neutrons to 2.2 Å resolution on the MaNDi neutron diffractometer at the Spallation Neutron Source. The resolution of the neutron data will allow direct determination of the positions of H atoms in the structure, which are believed to be fundamentally important in tuning the individual excitation energies of bacteriochlorophylls in this archetypal photosynthetic antenna complex. This is one of the largest unit-cell systems yet studied using neutron diffraction, and will allow the first high-resolution neutron analysis of a photosynthetic antenna complex.

Strongly Coupled Phenazine-Porphyrin Dyads: Light-Harvesting Molecular Assemblies with Broad Absorption Coverage

The development of light-harvesting architectures with broad absorption coverage in the visible region continues to be an important research area in the field of artificial photosynthesis. Here, we introduce a new class of ethynyl-linked panchromatic dyads composed of dibenzophenazines coupled ortho and meta to tetrapyrroles with an anchoring group that can be grafted onto metal oxide surfaces. Quantum chemical calculations and photophysical measurements of the synthesized materials reveal that both of the dibenzophenazine dyads absorb broadly from 300 to 636 nm and exhibit absorption bands different from those of the constituent chromophore units. Moreover, the different points of attachment of dibenzophenazines to tetrapyrroles give different absorption profiles which computations suggest result from differences in the planarity of the two dyads. Applicability of the dyads in artificial photosynthesis systems was assessed by their incorporation and characterization of their performance in dye-sensitized solar cells.

Coarse-grained representation of the quasi adiabatic propagator path integral for the treatment of non-Markovian long-time bath memory

The description of non-Markovian effects imposed by low frequency bath modes poses a persistent challenge for path integral based approaches like the iterative quasi-adiabatic propagator path integral (iQUAPI) method. We present a novel approximate method, termed mask assisted coarse graining of influence coefficients (MACGIC)-iQUAPI, that offers appealing computational savings due to substantial reduction of considered path segments for propagation. The method relies on an efficient path segment merging procedure via an intermediate coarse grained representation of Feynman-Vernon influence coefficients that exploits physical properties of system decoherence. The MACGIC-iQUAPI method allows us to access the regime of biological significant long-time bath memory on the order of hundred propagation time steps while retaining convergence to iQUAPI results. Numerical performance is demonstrated for a set of benchmark problems that cover bath assisted long range electron transfer, the transition from coherent to incoherent dynamics in a prototypical molecular dimer and excitation energy transfer in a 24-state model of the Fenna-Matthews-Olson trimer complex where in all cases excellent agreement with numerically exact reference data is obtained.

Excitonic Energy Landscape of the Y16F Mutant of the Chlorobium tepidum Fenna-Matthews-Olson (FMO) Complex: High Resolution Spectroscopic and Modeling Studies

We report high-resolution (low-temperature) absorption, emission, and nonresonant/resonant hole-burned (HB) spectra and results of excitonic calculations using a non-Markovian reduced density matrix theory (with an improved algorithm for parameter optimization in heterogeneous samples) obtained for the Y16F mutant of the Fenna-Matthews-Olson (FMO) trimer from the green sulfur bacterium Chlorobium tepidum. We show that the Y16F mutant is a mixture of FMO complexes with three independent low-energy traps (located near 817, 821, and 826 nm), in agreement with measured composite emission and HB spectra. Two of these traps belong to mutated FMO subpopulations characterized by significantly modified low-energy excitonic states. Hamiltonians for the two major subpopulations (Sub₈₂₁ and Sub₈₁₇) provide new insight into extensive changes induced by the single-point mutation in the vicinity of BChl 3 (where tyrosine Y16 was replaced with phenylalanine F16). The average decay time(s) from the higher exciton state(s) in the Y16F mutant depends on frequency and occurs on a picosecond time scale.

EXAT: EXcitonic analysis tool

We introduce EXcitonic Analysis Tool (EXAT), a program able to compute optical spectra of large excitonic systems directly from the output of quantum mechanical calculations performed with the popular Gaussian 16 package. The software is able to combine in an excitonic scheme the single-chromophore properties and exciton couplings to simulate energies, coefficients, and excitonic spectra (UV-vis, CD, and LD). The effect of the environment can also be included using a Polarizable Continuum Model. EXAT also presents a simple graphical user interface, which shows on-screen both site and exciton properties. To show the potential of the method, we report two applications on a a chiral perturbed BODIPY system and DNA G-quadruplexes, respectively. The program is available online at http://molecolab.dcci.unipi.it/tools/. © 2017 Wiley Periodicals, Inc.

Photosynthesis, pigment-protein complexes and electronic energy transport: simple models for complicated processes

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.

Nature does not rely on long-lived electronic quantum coherence for photosynthetic energy transfer

During the first steps of photosynthesis, the energy of impinging solar photons is transformed into electronic excitation energy of the light-harvesting biomolecular complexes. The subsequent energy transfer to the reaction center is commonly rationalized in terms of excitons moving on a grid of biomolecular chromophores on typical timescales [Formula: see text]100 fs. Today's understanding of the energy transfer includes the fact that the excitons are delocalized over a few neighboring sites, but the role of quantum coherence is considered as irrelevant for the transfer dynamics because it typically decays within a few tens of femtoseconds. This orthodox picture of incoherent energy transfer between clusters of a few pigments sharing delocalized excitons has been challenged by ultrafast optical spectroscopy experiments with the Fenna-Matthews-Olson protein, in which interference oscillatory signals up to 1.5 ps were reported and interpreted as direct evidence of exceptionally long-lived electronic quantum coherence. Here, we show that the optical 2D photon echo spectra of this complex at ambient temperature in aqueous solution do not provide evidence of any long-lived electronic quantum coherence, but confirm the orthodox view of rapidly decaying electronic quantum coherence on a timescale of 60 fs. Our results can be considered as generic and give no hint that electronic quantum coherence plays any biofunctional role in real photoactive biomolecular complexes. Because in this structurally well-defined protein the distances between bacteriochlorophylls are comparable to those of other light-harvesting complexes, we anticipate that this finding is general and directly applies to even larger photoactive biomolecular complexes.

Photosynthetic pigment-protein complexes as highly connected networks: implications for robust energy transport

Photosynthetic pigment-protein complexes (PPCs) are a vital component of the light-harvesting machinery of all plants and photosynthesizing bacteria, enabling efficient transport of the energy of absorbed light towards the reaction centre, where chemical energy storage is initiated. PPCs comprise a set of chromophore molecules, typically bacteriochlorophyll species, held in a well-defined arrangement by a protein scaffold; this relatively rigid distribution leads to a viewpoint in which the chromophore subsystem is treated as a network, where chromophores represent vertices and inter-chromophore electronic couplings represent edges. This graph-based view can then be used as a framework within which to interrogate the role of structural and electronic organization in PPCs. Here, we use this network-based viewpoint to compare excitation energy transfer (EET) dynamics in the light-harvesting complex II (LHC-II) system commonly found in higher plants and the Fenna-Matthews-Olson (FMO) complex found in green sulfur bacteria. The results of our simple network-based investigations clearly demonstrate the role of network connectivity and multiple EET pathways on the efficient and robust EET dynamics in these PPCs, and highlight a role for such considerations in the development of new artificial light-harvesting systems.

Probing the excitonic landscape of the Chlorobaculum tepidum Fenna-Matthews-Olson (FMO) complex: a mutagenesis approach

In this paper we report the steady-state optical properties of a series of site-directed mutants in the Fenna-Matthews-Olson (FMO) complex of Chlorobaculum tepidum, a photosynthetic green sulfur bacterium. The FMO antenna complex has historically been used as a model system for energy transfer due to the water-soluble nature of the protein, its stability at room temperature, as well as the availability of high-resolution structural data. Eight FMO mutants were constructed with changes in the environment of each of the bacteriochlorophyll a pigments found within each monomer of the homotrimeric FMO complex. Our results reveal multiple changes in low temperature absorption, as well as room temperature CD in each mutant compared to the wild-type FMO complex. These datasets were subsequently used to model the site energies of each pigment in the FMO complex by employing three different Hamiltonians from the literature. This enabled a basic approximation of the site energy shifts imparted on each pigment by the changed amino acid residue. These simulations suggest that, while the three Hamiltonians used in this work provide good fits to the wild-type FMO absorption spectrum, further efforts are required to obtain good fits to the mutant minus wild-type absorption difference spectra. This demonstrates that the use of FMO mutants can be a valuable tool to refine and iterate the current models of energy transfer in this system.

Multi-layer multi-configuration time-dependent Hartree (ML-MCTDH) approach to the correlated exciton-vibrational dynamics in the FMO complex

The coupled quantum dynamics of excitonic and vibrational degrees of freedom is investigated for high-dimensional models of the Fenna-Matthews-Olson complex. This includes a seven- and an eight-site model with 518 and 592 harmonic vibrational modes, respectively. The coupling between local electronic transitions and vibrations is described within the Huang-Rhys model using parameters that are obtained by discretization of an experimental spectral density. Different pathways of excitation energy flow are analyzed in terms of the reduced one-exciton density matrix, focussing on the role of vibrational and vibronic excitations. Distinct features due to both competing time scales of vibrational and exciton motion and vibronically assisted transfer are observed. The question of the effect of initial state preparation is addressed by comparing the case of an instantaneous Franck-Condon excitation at a single site with that of a laser field excitation.