Ferritin (Fn) proteins or their isolated subunits can be used as biomolecular templates for the selectively heterogeneous nucleation and growth of nanoparticles, in particular of iron oxyhydroxides. To shed light on the atomistic mechanisms of ferritin-promoted mineralization, in this study we perform molecular dynamics simulations to investigate the anchoring sites for Fe(III) clusters on Fn subunit assemblies using models of goethite and ferrihydrite nanoparticles. For this aim, we develop and parametrize a classical force field for Fe(III) oxyhydroxides based on reference density functional theory calculations. We then reveal that stable Fn-nanoparticle contacts are formed not only via negatively charged amino acid residues (glutamic and aspartic acid) but also, in a similar amount, via positively charged (lysine and arginine) and neutral (histidine) residues. A large majority of the anchoring sites are situated at the inner side of protein cages, consistent with the natural iron storage function of ferritin in many organisms. A slightly different distribution of anchoring sites is observed on heavy (H) and light (L) Fn subunits, with the former offering a larger amount of negative and neutral sites than the latter. This finding is exploited to develop a Fn mineralization protocol in which immobilized Fn subunits are first loaded with Fe2+ ions in a long "activation" step before starting their oxidation to Fe3+. This leads to the formation of very dense and uniform iron oxide films, especially when H subunits are employed.