Accurate prediction of mutation-induced frequency shifts in chlorophyll proteins with a simple electrostatic model

Controlling Vibronic Coupling in Chlorophyll Proteins: The Effects of Excitonic Delocalization and Vibrational Localization

Vibrational-electronic (vibronic) coupling plays a critical role in excitation energy transfer in molecular aggregates and pigment-protein complexes (PPCs). But the interplay between excitonic delocalization and vibronic interactions is complex, often leaving even qualitative questions as to what conceptual framework (e.g., Redfield versus Förster theory) should be used to interpret experimental results. To shed light on this issue, we report here on the interplay between excitonic delocalization and vibronic coupling in site-directed mutants of the water-soluble chlorophyll protein (WSCP), as reflected in 77 K fluorescence spectra. Experimentally, we find that in PPCs where excitonic delocalization is disrupted (either by mutagenesis or heterodimer formation), the relative intensity of the vibrational sideband (VSB) in fluorescence spectra is suppressed by up to 37% compared to that of the native protein. Numerical simulations reveal that this effect results from the localization of high-frequency vibrations in the coupled system; while excitonic delocalization suppresses the purely electronic transition due to H-aggregate-like dipole-dipole interference, high-frequency vibrations are unaffected, leading to a relative enhancement of the VSB. By comparing VSB intensities of PPCs both in the presence and absence of excitonic delocalization, we extract a set of "local" Huang-Rhys (HR) factors for Chl a in WSCP. More generally, our results suggest a significant role for geometric effects in controlling energy-transfer rates (which depend sensitively on absorption/fluorescence line shapes) in molecular aggregates and PPCs.

De novo design of proteins housing excitonically coupled chlorophyll special pairs

Natural photosystems couple light harvesting to charge separation using a 'special pair' of chlorophyll molecules that accepts excitation energy from the antenna and initiates an electron-transfer cascade. To investigate the photophysics of special pairs independently of the complexities of native photosynthetic proteins, and as a first step toward creating synthetic photosystems for new energy conversion technologies, we designed C2-symmetric proteins that hold two chlorophyll molecules in closely juxtaposed arrangements. X-ray crystallography confirmed that one designed protein binds two chlorophylls in the same orientation as native special pairs, whereas a second designed protein positions them in a previously unseen geometry. Spectroscopy revealed that the chlorophylls are excitonically coupled, and fluorescence lifetime imaging demonstrated energy transfer. The cryo-electron microscopy structure of a designed 24-chlorophyll octahedral nanocage with a special pair on each edge closely matched the design model. The results suggest that the de novo design of artificial photosynthetic systems is within reach of current computational methods.

Impact of Peripheral Hydrogen Bond on Electronic Properties of the Primary Acceptor Chlorophyll in the Reaction Center of Photosystem I

Photosystem I (PS I) is a photosynthetic pigment-protein complex that absorbs light and uses the absorbed energy to initiate electron transfer. Electron transfer has been shown to occur concurrently along two (A- and B-) branches of reaction center (RC) cofactors. The electron transfer chain originates from a special pair of chlorophyll a molecules (P700), followed by two chlorophylls and one phylloquinone in each branch (denoted as A-1, A0, A1, respectively), converging in a single iron-sulfur complex Fx. While there is a consensus that the ultimate electron donor-acceptor pair is P700+A0-, the involvement of A-1 in electron transfer, as well as the mechanism of the very first step in the charge separation sequence, has been under debate. To resolve this question, multiple groups have targeted electron transfer cofactors by site-directed mutations. In this work, the peripheral hydrogen bonds to keto groups of A0 chlorophylls have been disrupted by mutagenesis. Four mutants were generated: PsaA-Y692F; PsaB-Y667F; PsaB-Y667A; and a double mutant PsaA-Y692F/PsaB-Y667F. Contrary to expectations, but in agreement with density functional theory modeling, the removal of the hydrogen bond by Tyr → Phe substitution was found to have a negligible effect on redox potentials and optical absorption spectra of respective chlorophylls. In contrast, Tyr → Ala substitution was shown to have a fatal effect on the PS I function. It is thus inferred that PsaA-Y692 and PsaB-Y667 residues have primarily structural significance, and their ability to coordinate respective chlorophylls in electron transfer via hydrogen bond plays a minor role.

PigmentHunter: A point-and-click application for automated chlorophyll-protein simulations

Chlorophyll proteins (CPs) are the workhorses of biological photosynthesis, working together to absorb solar energy, transfer it to chemically active reaction centers, and control the charge-separation process that drives its storage as chemical energy. Yet predicting CP optical and electronic properties remains a serious challenge, driven by the computational difficulty of treating large, electronically coupled molecular pigments embedded in a dynamically structured protein environment. To address this challenge, we introduce here an analysis tool called PigmentHunter, which automates the process of preparing CP structures for molecular dynamics (MD), running short MD simulations on the nanoHUB.org science gateway, and then using electrostatic and steric analysis routines to predict optical absorption, fluorescence, and circular dichroism spectra within a Frenkel exciton model. Inter-pigment couplings are evaluated using point-dipole or transition-charge coupling models, while site energies can be estimated using both electrostatic and ring-deformation approaches. The package is built in a Jupyter Notebook environment, with a point-and-click interface that can be used either to manually prepare individual structures or to batch-process many structures at once. We illustrate PigmentHunter's capabilities with example simulations on spectral line shapes in the light harvesting 2 complex, site energies in the Fenna-Matthews-Olson protein, and ring deformation in photosystems I and II.

Detangling Solvatochromic Effects by the Effective Fragment Potential Method

Understanding molecular interactions in complex systems opens avenues for the efficient design of new materials with target properties. Energy decomposition methods provide a means to obtain a detailed picture of intermolecular interactions. This work introduces a molecular modeling approach for decomposing the solvatochromic shifts of the electronic excited states into the contributions of the individual molecular fragments of the environment surrounding the chromophore. The developed approach is implemented for the QM/EFP (quantum mechanics/effective fragment potential) model that provides a rigorous first-principles-based description of the electronic states of the chromophores in complex polarizable environments. On the example of two model systems, water pentamer and hydrated uracil, we show how the decomposition of the solvatochromic shifts into the contributions of individual solvent water molecules provides a detailed picture of the intermolecular interactions in the ground and excited states of these systems. The analysis also demonstrates the nonadditivity of solute-solvent interactions and the significant contribution of solute polarization to the total values of solvatochromic shifts.

Bioexcitons by Design: How Do We Get There?

Biological pigment-protein complexes (PPCs) exhibit a remarkable ability to tune the optical properties of biological excitons (bioexcitons) through specific pigment-protein interactions. While such fine-tuning allows natural systems (e.g., photosynthetic proteins) to carry out their native functions with near-optimal performance, native function itself is often suboptimal for applications such as biofuel production or quantum technology development. This perspective offers a look at near-term prospects for the rational reoptimization of PPC bioexcitons for new functions using site-directed mutagenesis. The primary focus is on the "structure-spectrum" challenge of understanding the relationships between structural features and spectroscopic properties. While recent examples demonstrate that site-directed mutagenesis can be used to tune nearly all key bioexciton parameters (e.g., site energies, interpigment couplings, and electronic-vibrational interactions), critical challenges remain before we achieve truly rational design of bioexciton properties.

Machine Learning Exciton Hamiltonians in Light-Harvesting Complexes

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.

Magnetophotoselection in the Investigation of Excitonically Coupled Chromophores: The Case of the Water-Soluble Chlorophyll Protein

A magnetophotoselection (MPS) investigation of the photoexcited triplet state of chlorophyll a both in a frozen organic solvent and in a protein environment, provided by the water-soluble chlorophyll protein (WSCP) of Lepidium virginicum, is reported. The MPS experiment combines the photoselection achieved by exciting with linearly polarized light with the magnetic selection of electron paramagnetic resonance (EPR) spectroscopy, allowing the determination of the relative orientation of the optical transition dipole moment and the zero-field splitting tensor axes in both environments. We demonstrate the robustness of the proposed methodology for a quantitative description of the excitonic interactions among pigments. The orientation of the optical transition dipole moments determined by the EPR analysis in WSCP, identified as an appropriate model system, are in excellent agreement with those calculated in the point-dipole approximation. In addition, MPS provides information on the electronic properties of the triplet state, localized on a single chlorophyll a pigment of the protein cluster, in terms of orientation of the zero-field splitting tensor axes in the molecular frame.