[论文解读] Computational Modeling of Exciton-bath Hamiltonians for LH2 and LH3 Complexes of Purple Photosynthetic Bacteria at Room Temperature

Daniel Villanueva Montemayor, Eva Rivera+1 more|arXiv (Cornell University)|Jan 2, 2026

Photosynthetic Processes and Mechanisms被引用 0

一句话总结

本文将全原子分子动力学（MD）模拟与TD-DFT及激子-浴模型相结合，解释LH2与LH3在光谱上的差异，并为这些复合体构建一个共同的激子-浴哈密顿量.

ABSTRACT

Light harvesting 2 (LH2) complex is a primary component of the photosynthetic unit of purple bacteria that is responsible for harvesting and relaying excitons. The electronic absorption line shape of LH2 contains two major bands at 800 nm and 850 nm wavelength regions. Under low light condition, some species of purple bacteria replace LH2 with LH3, a variant form with almost the same structure as the former but with distinctively different spectral features. The major difference between the absorption line shapes of LH2 and LH3 is the shift of the 850 nm band of the former to a new 820 nm region. The microscopic origin of this difference has been subject to some theoretical/computational investigations. However, the genuine molecular level source of such difference is not clearly understood yet. This work reports a comprehensive computational study of LH2 and LH3 complexes so as to clarify different molecular level features of LH2 and LH3 complexes and to construct simple exciton-bath models with a common form. All-atomistic molecular dynamics (MD) simulations of both LH2 and LH3 complexes provide detailed molecular level structural differences of BChls in the two complexes, in particular, in their patterns of hydrogen bonding (HB) and torsional angles of the acetyl group. Time-dependent density functional theory calculation of the excitation energies of BChls for structures sampled from the MD simulations, suggests that the observed differences in HB and torsional angles cannot fully account for the experimentally observed spectral shift of LH3. Potential sources that can explain the actual spectral shift of LH3 are discussed, and their magnitudes are assessed through fitting of experimental line shapes.

研究动机与目标

• 澄清在分子水平上区分LH2与LH3光谱特性的因素。
• 开发可用于LH2和LH3的简单、共同的激子-浴模型。
• 评估氢键和乙酰基取向如何影响Qy激发能。
• 提供基于MD得到的几何和耦合参数，以用于激子动力学模拟。

提出的方法

• 在膜环境中的LH2和LH3的全原子MD模拟，以获得结构参数和动力学信息。
• 对来自MD快照的相邻蛋白残基的BChl的Qy激发能进行TD-DFT计算。
• 转移电荷近似（TrEsp）用于从MD数据在介电筛选因子f=0.55的条件下计算电子耦合J_ss'(n-m)。
• 包括H_e^0、H_e、H_b和H_eb的粗粒度激子-浴哈密顿量H_T。
• 从MD得到的经典能隙相关函数C_E,sn^cl(t)通过傅里叶变换获得浴谱密度J_s(n)(ω)。
• 讨论除了HB效应之外的来源（如乙酰基的旋转），它们可能解释LH3的光谱位移。

实验结果

研究问题

• RQ1 哪些分子特征区分LH2和LH3及其对吸收光谱的影响？
• RQ2 是否可以用MD信息参数的共同激子-浴哈密顿量充分描述LH2和LH3的动力学？
• RQ3 氢键、乙酰基取向及局部环境在B800/B850（B820）光谱位移中贡献到何种程度？
• RQ4 MD推导的几何参数和电子耦合如何影响LH2与LH3的组装线形？

主要发现

Parameter	LH2 (cm^-1)	LH3 (cm^-1)	rel. diff. %
J_{αβ}(0)	244.74 ± 4.97	238.74 ± 5.11	2.45
J_{αβ}(1)	139.75 ± 10.22	74.13 ± 10.69	46.96
J_{αα}(1)	-59.07 ± 1.01	-61.21 ± 1.01	3.62
J_{ββ}(1)	-29.27 ± 1.61	-21.54 ± 1.70	26.41
J_{βα}(1)	13.83 ± 0.17	13.52 ± 0.18	2.24
J_{αβ}(2)	12.66 ± 0.35	10.83 ± 0.40	14.45
J_{γγ}(1)	-24.46 ± 0.77	-23.45 ± 0.75	4.13
J_{αγ}(1)	28.29 ± 0.59	26.58 ± 0.52	6.04
J_{βγ}(1)	2.63 ± 0.61	3.31 ± 0.54	25.86
J_{γα}(0)	-12.31 ± 0.42	-12.17 ± 0.37	1.14
J_{γβ}(0)	-2.68 ± 0.83	0.88 ± 0.77	67.16
λ_{α}	23.9	27.1	13.39
λ_{β}	28.4	26.4	7.04
λ_{γ}	73.3	72.4	1.23

• MD采样显示LH2与LH3在三级结构和色素排列上几乎相同，只有方位角和一些角度存在较小差异。
• 来自相邻单体的α-和β-BChl之间最近邻耦合（J_αβ(1)）在LH3中大约比LH2下降了一半。
• 在MD轨迹上平均的电子耦合通常小于固定晶体结构估计，与最近的高水平计算一致。
• 谱密度J_s(n)(ω)显示γ-BChl（B800）具有更大的浴耦合和较多的低频成分，原因是周围环境更亲水。
• 对于LH2和LH3，可行的共同激子-浴模型结构是可行的，能级无序用高斯分布建模，分离耦合通过TrEsp方法处理。
• 需要对能级无序进行建模以再现实验组态线形；建议无序的标准差大致在~200 cm^-1量级，适用于B850/B820单元。