The control of nanoparticle shape offers promise for improving catalytic activity and selectivity through optimization of the structure of the catalytically active site. Here, we have employed density functional theory calculations with a correction for the long-range interactions (DFT-D2) to investigate the effect of adsorption of the amino acid cysteine on the {001}, {011}, {100}, and {111} surfaces of mackinawite, which are commonly found in FeS nanoparticles. We have calculated the surface energies and adsorption energies for all the surfaces considered, and compared the surface energies of the pure and adsorbed systems. Based on the calculated surface energies, we have simulated the thermodynamic crystal morphology of the pure and cysteine-modified FeS nanoparticles using Wulff's construction. The strength of cysteine adsorption is found to be related to the stability of different surfaces, where it adsorbs most strongly onto the least stable FeS{111} surface via bidentate Fe–S and Fe–N chemical bonds and most weakly onto the most stable FeS{001} surface via hydrogen-bonded interactions; the adsorption energy decreases in the order {111} > {100} > {011} > {001}. We demonstrate that the stronger binding of the cysteine to the {011}, {100}, and {111} surfaces rather than to the {001} facet results in shape modulation of the FeS nanoparticles, with the reactive surfaces more expressed in the thermodynamic crystal morphology compared to the unmodified FeS crystals. Information regarding the structural parameters, electronic structures and vibrational frequency assignments of the cysteine–FeS complexes is also presented.