Theoretical description of photoinduced electron transfer in donor-acceptor supramolecular complexes based on carbon buckybowls.

Herein, we explore, from a theoretical perspective, the nonradiative photoinduced processes (charge separation and energy transfer) within a family of donor-acceptor supramolecular complexes based on the electron-donor truxene-tetrathiafulvalene (truxTTF) derivative and a series of curved fullerene fragments (buckybowls) of different shapes and sizes (C30H12, C32H12, and C38H14) as electron acceptors that successfully combine with truxTTF via non-covalent interactions. The resulting supramolecular complexes (truxTTF·C30H12, truxTTF·C32H12, and truxTTF·C38H14) undergo charge-separation processes upon photoexcitation through charge-transfer states involving the donor and acceptor units. Despite the not so different size of the buckybowls, they present noticeable differences in the charge-separation efficiency owing to a complex decay post-photoexcitation mechanism involving several low-lying excited states of different natures (local and charge-transfer excitations), all closely spaced in energy. In this intricate scenario, we have adopted a theoretical approach combining electronic structure calculations at (time-dependent) density functional theory, a multistate multifragment diabatization method, the Marcus-Levitch-Jortner semiclassical rate expression, and a kinetic model to estimate the charge separation rate constants of the supramolecular heterodimers. Our outcomes highlight that the efficiency of the photoinduced charge-separation process increases with the extension of the buckybowl backbone. The supramolecular heterodimer with the largest buckybowl (truxTTF·C38H14) displays multiple and efficient electron-transfer pathways, providing a global photoinduced charge separation in the ultrafast time scale in line with the experimental findings. The study reported indicates that modifications in the shape and size of buckybowl systems can give rise to attractive novel acceptors for potential photovoltaic applications.

Probing the Effect of Mutations on Light Harvesting in CP29 by Transient Absorption and First-Principles Simulations.

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.

Insights into Energy Transfer in Light-Harvesting Complex II Through Machine-Learning Assisted Simulations.

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.

Predicting Solvatochromism of Chromophores in Proteins through QM/MM and Machine Learning.

Solvatochromism occurs in both homogeneous solvents and more complex biological environments, such as proteins. While in both cases the solvatochromic effects report on the surroundings of the chromophore, their interpretation in proteins becomes more complicated not only because of structural effects induced by the protein pocket but also because the protein environment is highly anisotropic. This is particularly evident for highly conjugated and flexible molecules such as carotenoids, whose excitation energy is strongly dependent on both the geometry and the electrostatics of the environment. Here, we introduce a machine learning (ML) strategy trained on quantum mechanics/molecular mechanics calculations of geometrical and electrochromic contributions to carotenoids' excitation energies. We employ this strategy to compare solvatochromism in protein and solvent environments. Despite the important specifities of the protein, ML models trained on solvents can faithfully predict excitation energies in the protein environment, demonstrating the robustness of the chosen descriptors.

Electronic Excited States from Physically Constrained Machine Learning.

Data-driven techniques are increasingly used to replace electronic-structure calculations of matter. In this context, a relevant question is whether machine learning (ML) should be applied directly to predict the desired properties or combined explicitly with physically grounded operations. We present an example of an integrated modeling approach in which a symmetry-adapted ML model of an effective Hamiltonian is trained to reproduce electronic excitations from a quantum-mechanical calculation. The resulting model can make predictions for molecules that are much larger and more complex than those on which it is trained and allows for dramatic computational savings by indirectly targeting the outputs of well-converged calculations while using a parametrization corresponding to a minimal atom-centered basis. These results emphasize the merits of intertwining data-driven techniques with physical approximations, improving the transferability and interpretability of ML models without affecting their accuracy and computational efficiency and providing a blueprint for developing ML-augmented electronic-structure methods.

The OpenMMPol library for polarizable QM/MM calculations of properties and dynamics.

We present a new library designed to provide a simple and straightforward way to implement QM/AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Applications) and other polarizable QM/MM (Molecular Mechanics) methods based on induced point dipoles. The library, herein referred to as OpenMMPol, is free and open-sourced and is engineered to address the increasing demand for accurate and efficient QM/MM simulations. OpenMMPol is specifically designed to allow polarizable QM/MM calculations of ground state energies and gradients and excitation properties. Key features of OpenMMPol include a modular architecture facilitating extensibility, parallel computing capabilities for enhanced performance on modern cluster architectures, a user-friendly interface for intuitive implementation, and a simple and flexible structure for providing input data. To show the capabilities offered by the library, we present an interface with PySCF to perform QM/AMOEBA molecular dynamics, geometry optimization, and excited-state calculation based on (time-dependent) density functional theory.

Structural and quantum chemical basis for OCP-mediated quenching of phycobilisomes.

Cyanobacteria use large antenna complexes called phycobilisomes (PBSs) for light harvesting. However, intense light triggers non-photochemical quenching, where the orange carotenoid protein (OCP) binds to PBS, dissipating excess energy as heat. The mechanism of efficiently transferring energy from phycocyanobilins in PBS to canthaxanthin in OCP remains insufficiently understood. Using cryo-electron microscopy, we unveiled the OCP-PBS complex structure at 1.6- to 2.1-angstrom resolution, showcasing its inherent flexibility. Using multiscale quantum chemistry, we disclosed the quenching mechanism. Identifying key protein residues, we clarified how canthaxanthin's transition dipole moment in its lowest-energy dark state becomes large enough for efficient energy transfer from phycocyanobilins. Our energy transfer model offers a detailed understanding of the atomic determinants of light harvesting regulation and antenna architecture in cyanobacteria.

The nature of carotenoid S* state and its role in the nonphotochemical quenching of plants.

In plants, light-harvesting complexes serve as antennas to collect and transfer the absorbed energy to reaction centers, but also regulate energy transport by dissipating the excitation energy of chlorophylls. This process, known as nonphotochemical quenching, seems to be activated by conformational changes within the light-harvesting complex, but the quenching mechanisms remain elusive. Recent spectroscopic measurements suggest the carotenoid S* dark state as the quencher of chlorophylls' excitation. By investigating lutein embedded in different conformations of CP29 (a minor antenna in plants) via nonadiabatic excited state dynamics simulations, we reveal that different conformations of the complex differently stabilize the lutein s-trans conformer with respect to the dominant s-cis one. We show that the s-trans conformer presents the spectroscopic signatures of the S* state and rationalize its ability to accept energy from the closest excited chlorophylls, providing thus a relationship between the complex's conformation and the nonphotochemical quenching.

Deciphering Photoreceptors Through Atomistic Modeling from Light Absorption to Conformational Response.

In this review, we discuss the successes and challenges of the atomistic modeling of photoreceptors. Throughout our presentation, we integrate explanations of the primary methodological approaches, ranging from quantum mechanical descriptions to classical enhanced sampling methods, all while providing illustrative examples of their practical application to specific systems. To enhance the effectiveness of our analysis, our primary focus has been directed towards the examination of applications across three distinct photoreceptors. These include an example of Blue Light-Using Flavin (BLUF) domains, a bacteriophytochrome, and the orange carotenoid protein (OCP) employed by cyanobacteria for photoprotection. Particular emphasis will be placed on the pivotal role played by the protein matrix in fine-tuning the initial photochemical event within the embedded chromophore. Furthermore, we will investigate how this localized perturbation initiates a cascade of events propagating from the binding pocket throughout the entire protein structure, thanks to the intricate network of interactions between the chromophore and the protein.

Ultrafast excited-state dynamics of Luteins in the major light-harvesting complex LHCII.

Carotenoid pigments are known to present a functional versatility when bound to light-harvesting complexes. This versatility originates from a strong correlation between a complex electronic structure and a flexible geometry that is easily tunable by the surrounding protein environment. Here, we investigated how the different L1 and L2 sites of the major trimeric light-harvesting complex (LHCII) of green plants tune the electronic structure of the two embedded luteins, and how this reflects on their ultrafast dynamics upon excitation. By combining molecular dynamics and quantum mechanics/molecular mechanics calculations, we found that the two luteins feature a different conformation around the second dihedral angle in the lumenal side. The s-cis preference of the lutein in site L2 allows for a more planar geometry of the π -conjugated backbone, which results in an increased degree of delocalization and a reduced excitation energy, explaining the experimentally observed red shift. Despite these remarkable differences, according to surface hopping simulations the two luteins present analogous ultrafast dynamics upon excitation: the bright S 2 state quickly decays (in ∼ 50 fs) to the dark intermediate S x , eventually ending up in the S 1 state. Furthermore, by employing two different theoretical approaches (i.e., Förster theory and an excitonic version of surface hopping), we investigated the experimentally debated energy transfer between the two luteins. With both approaches, no evident energy transfer was observed in the ultrafast timescale.

First-principles simulation of excitation energy transfer and transient absorption spectroscopy in the CP29 light-harvesting complex.

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.

How orange carotenoid protein controls the excited state dynamics of canthaxanthin.

Orange Carotenoid Protein (OCP) is a ketocarotenoid-binding protein essential for photoprotection in cyanobacteria. The main steps of the photoactivated conversion which converts OCP from its resting state to the active one have been extensively investigated. However, the initial photochemical event in the ketocarotenoid which triggers the large structural changes finally leading to the active state is still not understood. Here we employ QM/MM surface hopping nonadiabatic dynamics to investigate the excited-state decay of canthaxanthin in OCP, both in the ultrafast S2 to S1 internal conversion and the slower decay leading back to the ground state. For the former step we show the involvement of an additional excited state, which in the literature has been often named the SX state, and we characterize its nature. For the latter step, we reveal an excited state decay characterized by multiple timescales, which are related to the ground-state conformational heterogeneity of the ketocarotenoid. We assigned the slowly decaying population to the so-called S* state. Finally, we identify a minor decay pathway involving double-bond photoisomerization, which could be the initial trigger to photoactivation of OCP.

Enabling Stereochemical Communication and Stimuli-Responsive Chiroptical Properties in Biphenyl-Capped Cyclodextrins.

Chiroptical materials are gaining increasing interest due to their innovative character and their applications in optoelectronics and data encryption technologies. Fully harnessing the potential of building blocks from the "chiral pool", such as native cyclodextrins (CDs), as they often lack chromophores suitable for the construction of materials with significant chiroptical properties. Here, we present the synthesis and characterization of a two-level molecular stack consisting of a point-chiral element (CD) and an axially chiral element (biphenyl), capable of effectively translating the overall stereochemical information contained in CDs into stimuli-responsive chiroptical properties. α- and ß-permethylated CDs were efficiently capped with two different 2,2'-difunctionalized 1,1'-biphenyl units. In CD derivatives containing the rigid 2,2'-dihydroxy-1,1'-biphenyl cap, two intramolecular hydrogen bonds act synergistically as stereoselective actuators, enabling effective communication between the two levels and the transfer of nonchromophoric stereochemical information from the cyclic-oligosaccharide to the atropoisomeric cap. The chiroptical properties can be finely tuned by external stimuli such as temperature and solvent. The way chirality is transferred from the CD platform to the biphenyl cap was revealed thanks to crystallographic and computational analyses, together with electronic circular dichroism (ECD) studies.

Origin of Low-Lying Red States in the Lhca4 Light-Harvesting Complex of Photosystem I.

The antenna complexes of Photosystem I present low-lying states visible as red-shifted and broadened absorption and fluorescence bands. Among these, Lhca4 has the most evident features of these "red" states, with a fluorescence band shifted by more than 25 nm from typical LHC emission. This signal arises from a mixing of exciton and charge-transfer (CT) states within the excitonically coupled a603-a609 chlorophyll (Chl) dimer. Here we combine molecular dynamics, multiscale quantum chemical calculations, and spectral simulations to uncover the molecular mechanism for the formation and tuning of exciton-CT interactions in Lhca4. We show that the coupling between exciton and CT states is extremely sensitive to tiny variations in the Chl dimer arrangement, explaining both the red-shifted bands and the switch between conformations with blue and red emission observed in single-molecule spectroscopy. Finally, we show that mutating the axial ligand of a603 diminishes the exciton-CT coupling, removing any red-state fingerprint.

Tidying up the conformational ensemble of a disordered peptide by computational prediction of spectroscopic fingerprints.

The most advanced structure prediction methods are powerless in exploring the conformational ensemble of disordered peptides and proteins and for this reason the "protein folding problem" remains unsolved. We present a novel methodology that enables the accurate prediction of spectroscopic fingerprints (circular dichroism, infrared, Raman, and Raman optical activity), and by this allows for "tidying up" the conformational ensembles of disordered peptides and disordered regions in proteins. This concept is elaborated for and applied to a dodecapeptide, whose spectroscopic fingerprint is measured and theoretically predicted by means of enhanced-sampling molecular dynamics coupled with quantum mechanical calculations. Following this approach, we demonstrate that peptides lacking a clear propensity for ordered secondary-structure motifs are not randomly, but only conditionally disordered. This means that their conformational landscape, or phase-space, can be well represented by a basis-set of conformers including about 10 to 100 structures. The implications of this finding have profound consequences both for the interpretation of experimental electronic and vibrational spectral features of peptides in solution and for the theoretical prediction of these features using accurate and computationally expensive techniques. The here-derived methods and conclusions are expected to fundamentally impact the rationalization of so-far elusive structure-spectra relationships for disordered peptides and proteins, towards improved and versatile structure prediction methods.

Integrated Computational Study of the Light-Activated Structure of the AppA BLUF Domain and Its Spectral Signatures.

We apply an integrated approach combining microsecond MD simulations and (polarizable) QM/MM calculations of NMR, FTIR, and UV-vis spectra to validate the structure of the light-activated form of the AppA photoreceptor, an example of blue light using flavin (BLUF) protein domain. The latter photoactivate through a proton-coupled electron transfer (PCET) that results in a tautomerization of a conserved glutamine residue in the active site, but this mechanism has never been spectroscopically proven for AppA, which has been always considered as an exception. Our simulations instead confirm that the spectral features observed upon AppA photoactivation are indeed directly connected to the tautomer form of glutamine as predicted by the PCET mechanism. In addition, we observe small but significant changes in the AppA structure, which are transmitted from the flavin binding pocket to the surface of the protein.

How the pH Controls Photoprotection in the Light-Harvesting Complex of Mosses.

In response to varying light conditions, light-harvesting complexes (LHCs) switch from a light-harvesting state to a quenched state to protect the photosynthetic organism from excessive light irradiation in a strategy known as nonphotochemical quenching (NPQ). NPQ is activated by an acidification of the thylakoid lumen, which is sensed directly or indirectly by the LHC, resulting in a conformational change of the complex that leads to the quenched state. The conformational changes responsible for NPQ activation and their connection to specific quenching mechanisms are still unknown. Here, we investigate the pH-triggered conformational changes in the light-harvesting complex stress-related (LHCSR) of mosses. By combining constant-pH molecular dynamics and enhanced sampling techniques, we find that the pH sensitivity of the complex is driven by the coupled protonation of three residues modulating the conformation of the short amphipathic helix placed at the lumen side of the embedding membrane. Combining these results with quantum mechanics/molecular mechanics calculations, we show that the quenching mechanism sensitive to the pH goes through a charge-transfer between a carotenoid and an excited chlorophyll, which is controlled by the protein conformation.

pH-Dependent Capping Interactions Induce Large-Scale Structural Transitions in i-Motifs.

We study here a DNA oligonucleotide having the ability to form two different i-motif structures whose relative stability depends on pH and temperature. The major species at neutral pH is stabilized by two C:C+ base pairs capped by two minor groove G:C:G:C tetrads. The high pH and thermal stability of this structure are mainly due to the favorable effect of the minor groove tetrads on their adjacent positively charged C:C+ base pairs. At pH 5, we observe a more elongated i-motif structure consisting of four C:C+ base pairs capped by two G:T:G:T tetrads. Molecular dynamics calculations show that the conformational transition between the two structures is driven by the protonation state of key cytosines. In spite of large conformational differences, the transition between the acidic and neutral structures can occur without unfolding of the i-motif. These results represent the first case of a conformational switch between two different i-motif structures and illustrate the dramatic pH-dependent plasticity of this fascinating DNA motif.