The three-dimensional structures of a number of [M(SR)(4)](n-) complexes, where M is a 3d transition metal and R is an alkyl or aryl group, have been analyzed using density functional theory (DFT). Special attention is paid to the Fe(II)/Fe(III) mimics of rubredoxin. The Fe(II) model complex [Fe(SCH(3))(4)](2-) has an equilibrium conformation with D2d symmetry. The DFT energy has been decomposed into contributions for ligand-ligand and metal-ligand interactions. The latter contribution is analyzed with the angular overlap model (AOM) and constitutes the dominant stereospecific interaction in the Fe(II) complex. The sulfur lone-pair electrons exert anisotropic pi interactions on the 3d(6) shell of Fe(II), which are controlled by the torsion angles, omega(i), for the rotations of the S(i)-C(beta) bonds around the Fe-S(i) axes. In contrast, the pi interactions acting on the high-spin 3d(5) shell of Fe(III) are isotropic. As a consequence, the stereochemistry of the Fe(III) complexes is determined by the Coulomb repulsions between the ligands and has S(4) symmetry. The electrostatic repulsions between the lone pairs of the sulfurs are an essential component of the ligand-ligand interaction. The lone-pair repulsions distort the 90 degree angle SFeS' angles (delta + delta(t)) and give rise to a correlation between delta and omega, which is confirmed by crystallographic data. Both the Fe(II) and Fe(III) complexes exhibit structural bistability due to the presence of low-lying equilibrium conformations with S(4) symmetry in which the complex can be trapped by the crystalline host.