Yin Bo-Wen, Wang Jie-Lei, Xue Pu-Jie, Zhang Teng-Shuo, Xie Bin-Bin, Shen Lin, Fang Wei-Hai
Hangzhou Institute of Advanced Studies, Zhejiang Normal University, 1108 Gengwen Road, Hangzhou 311231 Zhejiang, P. R. China.
Key Laboratory of Theoretical and Computational Photochemistry of Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, P. R. China.
J Chem Inf Model. 2023 Aug 14;63(15):4679-4690. doi: 10.1021/acs.jcim.3c00640. Epub 2023 Jul 25.
The contradictory behaviors in light harvesting and non-photochemical quenching make xanthophyll lutein the most attractive functional molecule in photosynthesis. Despite several theoretical simulations on the spectral properties and excited-state dynamics, the atomic-level photophysical mechanisms need to be further studied and established, especially for an accurate description of geometric and electronic structures of conical intersections for the lowest several electronic states of lutein. In the present work, semiempirical OM2/MRCI and multi-configurational restricted active space self-consistent field methods were performed to optimize the minima and conical intersections in and between the 1Ag, 2Ag, 1Bu, and 1Bu states. Meanwhile, the relative energies were refined by MS-CASPT2(10,8)/6-31G*, which can reproduce correct electronic state properties as those in the spectroscopic experiments. Based on the above calculation results, we proposed a possible excited-state relaxation mechanism for lutein from its initially populated 1Bu state. Once excited to the optically bright 1Bu state, the system will propagate along the key reaction coordinate, i.e., the stretching vibration of the conjugated carbon chain. During this period of time, the 1Bu state will participate in and forms a resonance state between the 1Bu and 1Bu states. Later, the system will rapidly hop to the 2Ag state via the 1Bu/2Ag conical intersection. Finally, the lutein molecule will survive in the 2Ag state for a relatively long time before it internally converts to the ground state directly or via a twisted S/S conical intersection. Notably, though the photophysical picture may be very different in solvents and proteins, the current theoretical study proposed a promising calculation protocol and also provided many valuable mechanistic insights for lutein and similar carotenoids.
光捕获和非光化学猝灭中的矛盾行为使叶黄素成为光合作用中最具吸引力的功能分子。尽管对其光谱性质和激发态动力学进行了若干理论模拟,但原子水平的光物理机制仍需进一步研究和确立,特别是要准确描述叶黄素最低几个电子态的锥形交叉点的几何和电子结构。在本工作中,采用半经验OM2/MRCI和多组态限制活性空间自洽场方法优化了¹Ag、²Ag、¹Bu和¹¹Bu态及其之间的极小值和锥形交叉点。同时,用MS-CASPT2(10,8)/6-31G*对相对能量进行了精细计算,该方法能够重现与光谱实验中一致的正确电子态性质。基于上述计算结果,我们提出了叶黄素从初始占据的¹¹Bu态开始的一种可能的激发态弛豫机制。一旦被激发到光学明亮的¹¹Bu态,系统将沿着关键反应坐标,即共轭碳链的拉伸振动进行传播。在此期间,¹¹Bu态将参与并在¹¹Bu和¹¹Bu态之间形成一个共振态。随后,系统将通过¹¹Bu/²Ag锥形交叉点迅速跃迁到²Ag态。最后,叶黄素分子将在²Ag态存活相对较长时间,然后直接或通过扭曲的S/S锥形交叉点内转换回基态。值得注意的是,尽管在溶剂和蛋白质中的光物理情况可能非常不同,但目前的理论研究提出了一种有前景的计算方案,也为叶黄素和类似类胡萝卜素提供了许多有价值的机理见解。