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 this type of composite systems can be only achieved with an atomistic description. In the last 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 LH complexes, we take as an example the plant antenna complex CP29 and its H111N mutant.

Optical spectra of the LHCII antenna complex are simulated using a fully atomistic first-principles exciton model, derived from multiscale quantum mechanical calculations of the pigments and a polarizable classical model for the environment.

Force-field parameters of the first row transition metals together with a few additional common elements such as those from the second (Rh, Ru) and third (Hg, Pt) rows of elements in ligated forms were determined based on the density functional theory calculations. Bonding characteristics were determined by averaging metal-ligand force constants in optimal geometries from several chosen complexes of each metal in the most common oxidation numbers and structural arrangements. Parameters of Lennard-Jones potential were determined based on a supermolecular model. Our determined molecular mechanical parameters are compared with presently available parameters published by other groups. We performed two different kinds of testing in order to demonstrate the reliability of these parameters in the case of ligated metallo complexes. First, the nonbonding potential was constructed for an additional set of 19 larger systems containing common complexes with organic molecules. The second test compares the Pt-O and Pt-H radial distribution functions for cisplatin in a box of TIP3P water with lately published studies.

Light-harvesting complexes (LHCs) are pigment-protein complexes whose main function is to capture sunlight and transfer the energy to reaction centers of photosystems. In response to varying light conditions, LH complexes also play photoregulation and photoprotection roles. In algae and mosses, a sub-family of LHCs, light-harvesting complex stress-related (LHCSR), is responsible for photoprotective quenching. Despite their functional and evolutionary importance, no direct structural information on LHCSRs is available that can explain their unique properties. In this work, we propose a structural model of LHCSR1 from the moss P. patens, obtained through an integrated computational strategy that combines homology modeling, molecular dynamics, and multiscale quantum chemical calculations. The model is validated by reproducing the spectral properties of LHCSR1. Our model reveals the structural specificity of LHCSR1, as compared with the CP29 LH complex, and poses the basis for understanding photoprotective quenching in mosses.

No abstract