Electronic structures of Zn 2+ and Cd 2+ thiolate clusters found in metallothioneins (MT) have been obtained using density functional theory. We have found that the inherent asymmetry of cluster architectures gives rise to seven distinct metal sites. Whereas the non-strained bond lengths of such tetrathiolate complexes are found to be 2.60Å and 2.39Å for Cd-S and Zn-S, in the MT clusters four characteristic terminal and bridging bonds are observed with average lengths 2.55Å (Cd-S t ); 2.35Å (Zn-S t ); 2.62Å (Cd-S b ); and 2.42Å (Zn-S b ). For each stoichiometry of Zn 2+ and Cd 2+ , all possible isomers have been characterized and ranked according to relative free energy and metal ion selectivity. The most stable distribution at low Cd 2+ concentration is computed to be Zn 4 + CdZn 2 , whereas at 2 : 1 Cd 2+ : Zn 2+ concentration, only heteroclusters are thermodynamically stable, explaining experimental data. The presence of two different clusters in MTs must and can be rationalized already in their intrinsic differences. The results indicate that the asymmetry allows for Zn 2+ transfer to various molecular targets having different thresholds for Zn 2+ binding, while maintaining detoxification sites.