This study proposes an innovative approach to estimate relative biological effectiveness (RBE) of fast neutrons using the Geant4 toolkit. The Geant4-DNA version cannot track heavy ions below 0.5 MeV/nucleon. In order to explore the impact of this issue, secondary particles are simulated instead of the primary low-energy neutrons. The Evaluated Nuclear Data File library is used to determine the cross sections for the elastic and inelastic interactions of neutrons with water and to find the contribution of each secondary particle spectrum. Two strategies are investigated in order to find the best possible approach and results. The first one takes into account only light particles, protons produced from elastic scattering, and α particles from inelastic scattering. Geantino particles are shot instead of heavy ions; hence all heavy ions are considered in the simulations, though their physical effects on DNA not. The second strategy takes into account all the heavy and light ions, although heavy ions cannot be tracked down to very low energies (E<0.5 MeV/nucleon). Our model is based on the combination of an atomic resolution DNA geometrical model and a Monte Carlo simulation toolkit for tracking particles. The atomic coordinates of the DNA double helix are extracted from the Protein Data Bank. Since secondary particle spectra are used instead of simulating the interaction of neutrons explicitly, this method reduces the computation times dramatically. Double-strand break induction is used as the end point for the estimation of the RBE of fast neutrons. ^{60}Co γ rays are used as the reference radiation quality. Both strategies succeed in reproducing the behavior of the RBE_{max} as a function of the incident neutron energy ranging from 0.1 to 14 MeV, including the position of its peak. A comparison of the behavior of the two strategies shows that for neutrons with energies less than 0.7 MeV, the effect of heavy ions would not be very significant, but above 0.7 MeV, heavy ions have an important role in neutron RBE.