Iron-sulfur proteins are an important class of electron carriers in a wide variety of biological reactions. Determining the intrinsic contribution of the metal site to the redox potential is crucial in understanding how the protein environment influences the overall redox properties of the Fe-S proteins. Here we combine density functional theory and coupled cluster methods with photodetachment spectroscopy to study the electronic structures and gas-phase redox potentials of the [Fe(SCH3)(4)](2-/-/0) and [Fe(SCH3)(3)](-/0) analogues of the rubredoxin redox site. The calculations show that oxidations of [Fe(SCH3)(4)](2-) and [Fe(SCH3)(4)](-) involve mainly the Fe 3d and S 3p orbitals, respectively. The calculated adiabatic and vertical detachment energies are in good agreement with the experiment for [Fe(SCH3)(3)](-) and [Fe(SCH3)(4)](-). The current results further confirm the "inverted level scheme" for the high-spin [1Fe] systems. The redox couple, [Fe(SCH3)(4)](-/2), which is the one found in rubredoxin, but cannot be accessed experimentally in the gas phase, was investigated using a thermodynamic cycle that relates it to the [Fe(SCH3)(3)](-/0) couple and the ligand association reaction, [Fe(SCH3)(3)](0/-) + SCH3- --> [Fe(SCH3)(4)](-/2-). The calculated reduction energy of [Fe(SCH3)(4)](-) (1.7 eV) compares well with the value (1.6 eV) estimated from the calculated bond energies and the experimental detachment energy of [Fe(SCH3)(3)](-). Thus, this thermodynamic cycle method can be used to estimate metal-ligand bonding energies and determine intrinsic reduction potentials from photodetachment experiments when the reduced forms are not stable in the gas phase.