The evaluation of skin doses during manipulation of radioactive sources can be a critical issue for which the most accurate calculation strategies available should be used. The aim of this work was to compare the results of the analytical approach used in VARSKIN with the simulation of radiation transport and interaction by Monte Carlo calculations in GAMOS (GEANT4-based Architecture for Medicine-Oriented Simulations), and to provide an accurate and versatile tool for the evaluation of skin doses from radionuclide sources of any realistic shape (e.g. cylindrical, parallelepiped), even in the presence of multiple interposed absorber layers. A set of 20 radionuclides (pure β, β-γ, Auger and γ emitters) from among the most frequently employed in nuclear medicine and laboratory practices were selected for comparison. We studied a point-like and a cylindrical source, in the presence of varying thicknesses of absorbing layers. We found a general agreement for most nuclides when the source was directly in contact with skin or in the presence of a thin layer of absorbing material. However, when the thickness of the absorber increased, significant differences were found for several nuclides. In these cases, the proposed method based on a dedicated Monte Carlo simulation could give more accurate results in a reasonable time, which could optimise accuracy when assessing skin doses in routine as well as incidental exposure scenarios.