Ubiquinone forms an integral part of the electron transport chain in cellular respiration and photosynthesis across a vast number of organisms. Prior experimental results have shown that the photosynthetic reaction center (RC) from Rhodobacter sphaeroides is only fully functional with a limited set of methoxy-bearing quinones, suggesting that specific interactions with this substituent are required to drive electron transport and the formation of quinol. The nature of these interactions has yet to be determined. Through parameterization of a CHARMM-compatible quinone force field and subsequent molecular dynamics simulations of the quinone-bound RC, we have investigated and characterized the protein interactions with the quinones in the QA and QB sites using both equilibrium simulation and thermodynamic integration. In particular, we identify a specific interaction between the 2-methoxy group of ubiquinone in the QB site and the amide nitrogen of GlyL225 that we implicate in locking the orientation of the 2-methoxy group, thereby tuning the redox potential difference between the quinones occupying the QA and QB sites. Disruption of this interaction leads to weaker binding in a ubiquinone analog that lacks a 2-methoxy group, a finding supported by reverse electron transfer EPR experiments of the
QA−QB− biradical and competitive binding assays.