We provide a perspective of the induced dipole formulation of polarizable QM/MM, showing how efficient implementations will enable their application to the modeling of dynamics, spectroscopy, and reactivity in complex biosystems.

The spectroscopic properties of light-harvesting (LH) antennae in photosyntehtic organisms represent a fingerprint that is unique for each specific pigment-protein complex. Because of that, spectroscopic observations are generally combined with structural data from X-ray crystallography to obtain an indirect representation of the excitonic properties of the system. Here, an alternative strategy is presented which goes beyond this empirical approach and introduces an ab initio computational description of both structural and electronic properties and their dependence on the temperature. The strategy is applied to the peripheral light-harvesting antenna complex (LH2) present in purple bacteria. By comparing this model with the one based on the crystal structure, a detailed, molecular level explanation of the absorption and circular dichroism (CD) spectra and their temperature dependence is achieved. The agreement obtained with the experiments at both low and room temperature lays the groundwork for an atomistic understanding of the excitation dynamics in the LH2 system.

The LH2 antenna complexes of purple bacteria occur, depending on light conditions, in various different spectroscopic forms, with a similar structure but different absorption spectra. The differences are related to point changes in the primary amino acid sequence, but the molecular-level relationship between these changes and the resulting spectrum is still not well understood. We undertook a systematic quantum chemical analysis of all the main factors that contribute to the exciton structure, looking at how the environment modulates site energies and couplings in the B800-850 and B800-820 spectroscopic forms of LH2. A multiscale approach combining quantum chemistry and an atomistic classical embedding has been used where mutual polarization effects between the two parts are taken into account. We find that the loss of hydrogen bonds following amino acid changes can only explain a part of the observed blue-shift in the B850 band. The coupling of excitonic states to charge-transfer states, which is different in the two forms, contributes with a similar amount to the overall blue-shift.

Nonlinear electronic spectroscopies represent one of the most powerful techniques to study complex multichromophoric architectures. For these systems, in fact, linear spectra are too congested to be used to disentangle the many coupled vibroelectronic processes that are activated. By using a 2D approach, instead, a clear picture can be achieved, but only when the recorded spectra are combined with a proper interpretative model. So far, this has been almost always achieved through parametrized exciton Hamiltonians that necessarily introduce biases and/or arbitrary assumptions. In this study, a first-principles approach is presented that combines accurate quantum chemical descriptions with state-of-the-art models for the environment through the use of atomistic and polarizable embeddings. Slow and fast bath dynamics, along with exciton transport between the pigments, are included. This approach is applied to the 2DES spectroscopy of the Light-Harvesting 2 (LH2) complex of purple bacteria. Simulations are extended over the entire visible-near-infrared spectral region to cover both carotenoid and bacteriochlorophyll signals. Our results provide an accurate description of excitonic properties and relaxation pathways, and give an unprecedented insight into the interpretation of the spectral signatures of the measured 2D signals.

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.