Calculation of iron-sulfur clusters in FeS protein active sites by broken symmetry density functional method