Abstract 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 A0A/A0B ligands (Met or Asn in the 688/668 position) on the redox potential of chlorophylls A0A/A0B and phylloquinones A1A/A1B. A combination of MD and SC approaches was used to estimate reorganization energies λ of the primary (λ1) and secondary (λ2) charge separation reactions, which were found to be independent of the active branch of electron transfer; in PS I from the wild type, λ1 was estimated to be 390±20mV, while λ2 was estimated to be higher at 445±15mV. MD and QC approaches were used to describe the effect of substituting Met688PsaA/Met668PsaB by Asn688PsaA/Asn668PsaB 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 A0A/A0B and phylloquinones A1A/A1B, which may explain the experimentally observed slowdown of secondary electron transfer in the M688NPsaA variant. This article is part of a Special Issue entitled: Photosynthesis research for sustainability: Keys to produce clean energy.