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.

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.

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.

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.

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.

Quantum chemical elucidation of a sevenfold symmetric bacterial antenna complex.

The light-harvesting complex 2 (LH2) of purple bacteria is one of the most studied photosynthetic antenna complexes. Its symmetric structure and ring-like bacteriochlorophyll arrangement make it an ideal system for theoreticians and spectroscopists. LH2 complexes from most bacterial species are thought to have eightfold or ninefold symmetry, but recently a sevenfold symmetric LH2 structure from the bacterium Mch. purpuratum was solved by Cryo-Electron microscopy. This LH2 also possesses unique near-infrared absorption and circular dichroism (CD) spectral properties. Here we use an atomistic strategy to elucidate the spectral properties of Mch. purpuratum LH2 and understand the differences with the most commonly studied LH2 from Rbl. acidophilus. Our strategy exploits a combination of molecular dynamics simulations, multiscale polarizable quantum mechanics/molecular mechanics calculations, and lineshape simulations. Our calculations reveal that the spectral properties of LH2 complexes are tuned by site energies and exciton couplings, which in turn depend on the structural fluctuations of the bacteriochlorophylls. Our strategy proves effective in reproducing the absorption and CD spectra of the two LH2 complexes, and in uncovering the origin of their differences. This work proves that it is possible to obtain insight into the spectral tuning strategies of purple bacteria by quantitatively simulating the spectral properties of their antenna complexes.

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.

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.

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.

Multiscale Models for Light-Driven Processes.

Multiscale models combining quantum mechanical and classical descriptions are a very popular strategy to simulate properties and processes of complex systems. Many alternative formulations have been developed, and they are now available in all of the most widely used quantum chemistry packages. Their application to the study of light-driven processes, however, is more recent, and some methodological and numerical problems have yet to be solved. This is especially the case for the polarizable formulation of these models, the recent advances in which we review here. Specifically, we identify and describe the most important specificities that the polarizable formulation introduces into both the simulation of excited-state dynamics and the modeling of excitation energy and electron transfer processes.

A Unique and Stable Polyproline I Helix Sorted out from Conformational Equilibrium by Solvent Polarity.

Polyproline I helical structures are often considered as the hidden face of their most famous geminal sibling, Polyproline II, as PPI is generally spotted only within a conformational equilibrium. We designed and synthesized a stable Polyproline I structure exploiting the striking tendency of (S)-indoline-2-carboxylic acid to drive the peptide bond conformation toward the cis amide isomer, when dissolved in polar solvents. The cooperative effect of only four amino acidic units is sufficient to form a preferential structure in solution. We shed light on this rare secondary structure with a thorough analysis of the spectroscopic and chiroptical properties of the tetramer, supported by X-ray crystallography and computational studies.

The atomistic modeling of light-harvesting complexes from the physical models to the computational protocol.

The function of light-harvesting complexes is determined by a complex network of dynamic interactions among all the different components: the aggregate of pigments, the protein, and the surrounding environment. Complete and reliable predictions on these types of composite systems can be only achieved with an atomistic description. In the last few decades, there have been important advances in the atomistic modeling of light-harvesting complexes. These advances have involved both the completeness of the physical models and the accuracy and effectiveness of the computational protocols. In this Perspective, we present an overview of the main theoretical and computational breakthroughs attained so far in the field, with particular focus on the important role played by the protein and its dynamics. We then discuss the open problems in their accurate modeling that still need to be addressed. To illustrate an effective computational workflow for the modeling of light harvesting complexes, we take as an example the plant antenna complex CP29 and its H111N mutant.

Quantum Chemical Modeling of the Photoinduced Activity of Multichromophoric Biosystems.

Multichromophoric biosystems represent a broad family with very diverse members, ranging from light-harvesting pigment-protein complexes to nucleic acids. The former are designed to capture, harvest, efficiently transport, and transform energy from sunlight for photosynthesis, while the latter should dissipate the absorbed radiation as quickly as possible to prevent photodamages and corruption of the carried genetic information. Because of the unique electronic and structural characteristics, the modeling of their photoinduced activity is a real challenge. Numerous approaches have been devised building on the theoretical development achieved for single chromophores and on model Hamiltonians that capture the essential features of the system. Still, a question remains: is a general strategy for the accurate modeling of multichromophoric systems possible? By using a quantum chemical point of view, here we review the advancements developed so far highlighting differences and similarities with the single chromophore treatment. Finally, we outline the important limitations and challenges that still need to be tackled to reach a complete and accurate picture of their photoinduced properties and dynamics.

Molecular Mechanisms of Activation in the Orange Carotenoid Protein Revealed by Molecular Dynamics.

Light-harvesting in photosynthesis is accompanied by photoprotective processes. In cyanobacteria, the photoprotective role is played by a specialized complex, the orange carotenoid protein, which is activated by strong blue-green light. This photoactivation involves a unique series of structural changes which terminate with an opening of the complex into two separate domains, one of which acts as a quencher for the light-harvesting complexes. Many experimental studies have tried to reveal the molecular mechanisms through which the energy absorbed by the carotenoid finally leads to the large conformational change of the complex. Here, for the first time, these mechanisms are revealed by simulating at the atomistic level the whole dynamics of the complex through an effective combination of enhanced sampling techniques. On the basis of our findings, we can conclude that the carotenoid does not act as a spring that, releasing its internal strain, induces the dissociation, as was previously proposed, but as a "latch" locking together the two domains. The photochemically triggered displacement of the carotenoid breaks this balance, allowing the complex to dissociate.

A Synthetic Oxygen Sensor for Plants Based on Animal Hypoxia Signaling.

Due to the involvement of oxygen in many essential metabolic reactions, all living organisms have developed molecular systems that allow adaptive physiological and metabolic transitions depending on oxygen availability. In mammals, the expression of hypoxia-response genes is controlled by the heterodimeric Hypoxia-Inducible Factor. The activity of this transcriptional regulator is linked mainly to the oxygen-dependent hydroxylation of conserved proline residues in its α-subunit, carried out by prolyl-hydroxylases, and subsequent ubiquitination via the E3 ligase von Hippel-Lindau tumor suppressor, which targets Hypoxia-Inducible Factor-α to the proteasome. By exploiting bioengineered versions of this mammalian oxygen sensor, we designed and optimized a synthetic device that drives gene expression in an oxygen-dependent fashion in plants. Transient assays in Arabidopsis (Arabidopsis thaliana) mesophyll protoplasts indicated that a combination of the yeast Gal4/upstream activating sequence system and the mammalian oxygen sensor machinery can be used effectively to engineer a modular, oxygen-inducible transcriptional regulator. This synthetic device also was shown to be selectively controlled by oxygen in whole plants when its components were expressed stably in Arabidopsis seedlings. We envision the exploitation of our genetically encoded controllers to generate plants able to switch gene expression selectively depending on oxygen availability, thereby providing a proof of concept for the potential of synthetic biology to assist agricultural practices in environments with variable oxygen provision.

Exciton properties and optical spectra of light harvesting complex II from a fully atomistic description.

We present a fully atomistic simulation of linear optical spectra (absorption, fluorescence and circular dichroism) of the Light Harvesting Complex II (LHCII) trimer using a hybrid approach, which couples a quantum chemical description of the chlorophylls with a classical model for the protein and the external environment (membrane and water). The classical model uses a polarizable Molecular Mechanics force field, thus allowing mutual polarization effects in the calculations of the excitonic properties. The investigation is performed both on the crystal structure and on structures generated by a µs long classical molecular dynamics simulation of the complex within a solvated membrane. The results show that this integrated approach not only provides a good description of the excitonic properties and optical spectra without the need for additional refinements of the excitonic parameters, but it also allows an atomistic investigation of the relative importance of electronic, structural and environment effects in determining the optical spectra.

Elucidating the role of structural fluctuations, and intermolecular and vibronic interactions in the spectroscopic response of a bacteriophytochrome.

We present the first comprehensive multiscale computational investigation of Resonance Raman, absorption and Circular Dichroism spectra of the resting state of the Deinococcus radiodurans phytochrome. The spectra are simulated in all their components, namely the energy position and the lineshapes of both the far-red and the blue bands. To achieve such a goal, we have combined a 4.5 µs MD simulation of the solvated dimeric phytochrome with a hybrid quantum mechanics/molecular mechanics (QM/MM) model, which accounts for both electrostatic and mutual polarization effects between the QM and the MM subsystems. A good agreement with experiments is found for all the three spectra. Moreover, we find a transient H-bond network within the binding pocket of the biliverdin chromophore that, unexpectedly, does not significantly affect the spectra. In parallel, we characterize the vibrations that are more strongly coupled to the biliverdin excitation, confirming the important role of the hydrogen-out-of-plane mode of its vinyl C-H together with the expected C[double bond, length as m-dash]C stretching of the double bond involved in the photoisomerization.

Polarizable embedding QM/MM: the future gold standard for complex (bio)systems?

Nowadays, hybrid QM/MM approaches are widely used to study (supra)molecular systems embedded in complex biological matrices. However, in their common formulation, mutual interactions between the quantum and classical parts are neglected. To go beyond such a picture, a polarizable embedding can be used. In this perspective, we focus on the induced point dipole formulation of polarizable QM/MM approaches and we show how efficient and linear scaling implementations have allowed their application to the modeling of complex biosystems. In particular, we discuss their use in the prediction of spectroscopies and in molecular dynamics simulations, including Born-Oppenheimer dynamics, enhanced sampling techniques and nonadiabatic descriptions. We finally suggest the theoretical and computational developments that still need to be achieved to overcome the limitations which have prevented so far larger diffusion of these methods.

Electronic couplings for photo-induced processes from subsystem time-dependent density-functional theory: The role of the diabatization.

Subsystem time-dependent density-functional theory (sTDDFT) making use of approximate non-additive kinetic energy (NAKE) functionals is known to be capable of describing excitation energy transfer processes in a variety of applications. Here, we show that sTDDFT, especially when combined with projection-based embedding (PbE), can be employed for the entire range of photo-induced electronic couplings essential for modeling photophysical properties of complex chemical and biological systems and therefore represents a complete toolbox for this class of problems. This means that it is capable of capturing the interaction/coupling associated with local- and charge-transfer (CT) excitons. However, this requires the choice of a reasonable diabatic basis. We therefore propose different diabatization strategies of the virtual orbital space in PbE-sTDDFT and show how CT excitations can be included in sTDDFT using NAKE functionals via a phenomenological approach. Finally, these electronic couplings are compared to couplings from a multistate fragment excitation difference (FED)-fragment charge difference (FCD) diabatization procedure. We show that both procedures, multistate FED-FCD and sTDDFT (with the right diabatization procedure chosen), lead to an overall good agreement for the electronic couplings, despite differences in their general diabatization strategy. We conclude that the entire range of photo-induced electronic couplings can be obtained using sTDDFT (with the right diabatization procedure chosen) in a black-box manner.