Objective: The scientific community has considered internal dosimetry by the Monte Carlo method the gold standard. However, there is a trade-off between simulation processing time and the statistical quality of the results that makes it a challenge to obtain accurate absorbed dose values in some situations, such as dose estimation in organs affected by cross-irradiation or limited computing power. Variance reduction techniques are used to reduce computational processing time without impairing the statistical quality of the results, such as tracking energy cutoff, secondary particle production threshold, and parallelism of different types of emissions from radionuclides. Approach: In this work, GATE Monte Carlo code and its variance reduction techniques were evaluated to calculate S values of organs from the International Commission on Radiological Protection (ICRP) report 110 male phantom for the lutetium-177, iodine-131, yttrium-90, and radium-223 radionuclides. The results are compared with the data from the OpenDose collaboration. Main Results: A cutoff of 5 MeV for local electron deposition and 2.0 mm of secondary particle production range resulted in a computational efficiency increase of 7.9 and 1.05 times, respectively. Simulation of ICRP 107 spectra-based source proved to be about 5 times more efficient when compared to a decay simulation using \verb|G4RadioactiveDecay| (Geant4-based radioactive decay processes). TLE (Track Length Estimator) and seTLE (Split Exponential Track Length Estimator) techniques were used to calculate the absorbed dose of photon emissions, resulting in computational efficiency up to 29.4 and 62.5 times higher when compared to traditional simulations, respectively. In particular, the seTLE technique accelerates the simulation time by up to 1426 times, achieving a statistical uncertainty of 10\% in volumes affected by cross-irradiation. Significance: The variance reduction techniques used in this work drastically reduced the simulation time and maintained the statistical quality of the calculated absorbed dose values, proving the feasibility of the use of the Monte Carlo method in internal dosimetry under challenging situations and making it viable for clinical routine or web applications.