Molecular dynamics study of the primary charge separation reactions in Photosystem I: Effect of the replacement of the axial ligands to the electron acceptor A₀

Molecular dynamics (MD) calculations, a semi-continuum (SC) approach, and quantum chemistry (QC) calculations were employed together to investigate the molecular mechanics of ultrafast charge separation reactions in Photosystem I (PS I) of Thermosynechococcus elongatus. A molecular model of PS I was developed with the aim to relate the atomic structure with electron transfer events in the two branches of cofactors. A structural flexibility map of PS I was constructed based on MD simulations, which demonstrated its rigid hydrophobic core and more flexible peripheral regions. The MD model permitted the study of atomic movements (dielectric polarization) in response to primary and secondary charge separations, while QC calculations were used to estimate the direct chemical effect of the A_0A/A_0B ligands (Met or Asn in the 688/668 position) on the redox potential of chlorophylls A_0A/A _0B and phylloquinones A_1A/A_1B. A combination of MD and SC approaches was used to estimate reorganization energies λ of the primary (λ₁) and secondary (λ₂) charge separation reactions, which were found to be independent of the active branch of electron transfer; in PS I from the wild type, λ₁ was estimated to be 390 ± 20 mV, while λ₂ was estimated to be higher at 445 ± 15 mV. MD and QC approaches were used to describe the effect of substituting Met688_PsaA/Met668_PsaB by Asn688_PsaA/Asn668_PsaB on the energetics of electron transfer. Unlike Met, which has limited degrees of freedom in the site, Asn was found to switch between two relatively stable conformations depending on cofactor charge. The introduction of Asn and its conformation flexibility significantly affected the reorganization energy of charge separation and the redox potentials of chlorophylls A_0A/A_0B and phylloquinones A_1A/A_1B, which may explain the experimentally observed slowdown of secondary electron transfer in the M688N_PsaA variant. This article is part of a Special Issue entitled: Photosynthesis research for sustainability: Keys to produce clean energy.