Molecular dynamics (MD) simulations are performed to investigate how ∑3〈110〉(121) symmetric tilt grain boundary (GB) affects point defects and defect clusters in He-doped α-iron at 300K in picosecond time scales. Molecular statics calculations are also performed and show that the formation energy is reduced in the GB, and the GB acts as a good sink for point defects, especially for interstitial He and self-interstitial atoms (SIAs). It is observed that the average size of HenVm (m>n) clusters becomes smaller in the GB-containing Fe system, where m and n represent the number of vacancies and He atoms in the cluster, respectively. It is also found that the number of HenV (n=2, 3) clusters in the GB region decreases, while the number of the HeV clusters increases. The GBs loaded with substitutional or interstitial helium atoms are found to facilitate the growth of helium clusters in the GB region.