[논문 리뷰] Computational Modeling of Exciton-bath Hamiltonians for LH2 and LH3 Complexes of Purple Photosynthetic Bacteria at Room Temperature

Daniel Villanueva Montemayor, Eva Rivera외 1명|arXiv (Cornell University)|2026. 01. 02.

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-유래 기하학적 매개변수와 결합 매개변수를 제공하여 엑시온 다이내믹스 시뮬레이션에 정보를 제공한다.

제안 방법

• 막 환경에서의 모든 원자 MD 시뮬레이션을 통해 구조 매개변수와 역동성을 얻는다.
• MD 스냅샷에서 근접한 단백질 잔기로부터 BChl의 Qy 들뜬 에너지를 TD-DFT로 계산한다.
• MD 데이터로부터 전자 커플링 J_ss'(n-m)을 유도하는 안정화 스크리닝 요인 f=0.55를 사용하는 전이 전하 근사(TrEsp) 방법.
• H_T를 포함하는 거시적 엑시온-배스 해밀토니안(H_e^0, 불규칙성 H_e, H_b, H_eb)
• MD로부터 얻은 에너지 격차의 고전적 상관 함수 C_E,sn^cl(t)로부터 배스 스펙트럼 밀도 J_s(n)(ω)를 푸리에 변환으로 구한다.
• LH3의 스펙트럴 이동을 설명할 수 있는 HB 효과를 넘는 소스(예: 아세틸기 회전)에 대한 논의

실험 결과

연구 질문

• RQ1 LH2와 LH3를 차별화하는 분자적 특징은 무엇이며 이들이 흡수 스펙트럼에 어떻게 영향을 미치는가?
• RQ2 MD 정보를 반영한 매개변수로 LH2와 LH3의 동역학을 충분히 설명하는 공통의 엑시온-배스 해밀토니안이 가능한가?
• RQ3 수소 결합, 아세틸기 방향성, 국소 환경이 B800/B850(B820) 스펙트럴 이동에 얼마나 기여하는가?
• RQ4 MD에서 얻은 기하학적 매개변수와 전자 커플링이 LH2와 LH3의 앙상블 선모양에 어떤 영향을 미치는가?

주요 결과

매개변수	LH2 (cm^-1)	LH3 (cm^-1)	상대 차이 %
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의 3차 구조 및 색소 배치가 거의 동일하며 방위 및 일부 각도에서 작은 차이가 나타난다.
• 인접 프로토머의 α- 및 β-BChl 간 최근접 커플링(J_αβ(1))은 LH3에서 LH2에 비해 약 절반으로 감소한다.
• MD 궤적에서 평균된 전자 커플링은 일반적으로 고정 결정 구조 추정값보다 작으며 최근의 고수준 계산과 일치한다.
• 스펙트럴 밀도 J_s(n)(ω)는 γ-BChl(B800)이 더 큰 배스 커플링과 더 많은 저주파 구성분을 보여 주며, 이는 더 친수성 환경 때문
• 공통 엑시온-배스 모델 구조가 LH2와 LH3에 대해 가능하며, 사이트 에너지의 불규칙성은 가우시안으로 모델링되고 역방향 커플링은 TrEsp 방법으로 처리한다.
• 사이트 에너지의 무작위성은 실험적 앙상블 선형모양을 재현하는 데 필요하다; 표준 편차는 B850/B820 단위에서 약 200 cm^-1 수준으로 제시된다.