Molecular orbital calculations of the active site complex of two iron ferredoxins

The results presented here are the first part of a systematic theoretical study of some of the physical and biochemical properties of two iron ferredoxins obtained by the use of an Extended Hückel Self-Consistent Charge iteration method of molecular orbital calculations. In this initial study, attention is focused on the calculation of electronic energies as a function of molecular geometry and the nature of the bonding ligands at the active site in order to determine the most stable form of the active site complex. Included in the active site complex are two iron atoms, two acid labile sulfur atoms of unknown inorganic origin and four sulfur atoms presumably from nearby cysteine residues. Fifteen chemical-conformational variations of this basic active site complex were considered. Among these conformational variations of the sulfur ligands, Fe-Fe distances, bond lengths and angles and chemical variations such as the effect of axial ligands, disulfide bonds and added protons were included. Our results indicate that that with all reasonable variations of the ligands, the preferred molecular geometry about 4-coordinated Fe is tetrahedral rather than planar. The planar conformation is somewhat stabilized by the addition of axial ligands, but is still less favorable than the tetrahedral conformation. In this model, interactions between the two iron atoms occur automatically since they are both part of the same active site complex. Hence the absence of low temperature paramagnetism in the oxidized state is readily explained. Preliminary investigations of the reduced state with one additional electron indicate that the odd electron is delocalized, as observed in both ESR and ENDOR. Its presence apparently substantially destabilizes all of the molecular orbital energies in accord with the observation that only one electron can be added to these proteins without decomposing them.

---

Zusammenfassung Die vorliegenden Ergebnisse stellen den ersten Teil einer systematischen theoretischen Untersuchung einiger physikalischer und biochemischer Eigenschaften des Fe₂-Ferredoxins mit Hilfe einer selbstkonsistenten erweiterten Hückelmethode dar. In diesem ersten Teil wird besonders die elektronische Energie in Abhängigkeit von der Molekülgeometrie und der Art der Liganden an der aktiven Stelle untersucht, um die stabilste Form des Komplexes der aktiven Stelle zu finden. Dieser Komplex enthält zwei Eisenatome, zwei saure, nicht fest gebundene Schwefelatome unbekannter anorganischer Herkunft sowie vier Schwefelatome, vermutlich aus nahegelegenen Cysteinresten. Bei den 15 untersuchten Konformationen der Schwefelliganden wurden Änderungen der Fe-Fe-Abstände, Bindungslängen und -winkel sowie chemische Variationen wie die Einwirkung der axialen Liganden, der Disulfidbindungen und zusätzlicher Protonen berücksichtigt. Es ergibt sich, daß das Eisen vierfach tetraedrisch koordiniert ist und nicht planar, selbst bei der Addition axialer Liganden. In diesem Modell ergeben sich Wechselwirkungen zwischen den beiden Eisenatomen zwanglos, da sie zum gleichen aktiven Komplex gehören. Dadurch wird das Fehlen eines Tieftemperatur-Paramagnetismus im oxydierten Zustand hinreichend erklärt. Eine vorläufige Untersuchung des reduzierten Zustandes mit einem zusätzlichen Elektron zeigt, daß dieses delokalisiert ist, genau wie es mit ESR und ENDOR beobachtet wird. seine Gegenwart destabilisiert offenbar alle MO-Energien, was in Übereinstimmung mit der Beobachtung steht, daß nur ein Elektron von diesen Proteinen addiert werden kann, ohne sie zu zerstören.