We use numerical N-body hydrodynamical simulations with varying PopIII stellar models to investigate the possibility of detecting first star signatures with observations of high-redshift damped Lya absorbers (DLAs). The simulations include atomic and molecular cooling, star formation, energy feedback and metal spreading due to the evolution of stars with a range of masses and metallicities. Different initial mass functions (IMFs) and corresponding metaldependent yields and lifetimes are adopted to model primordial stellar populations. The DLAs in the simulations are selected according to either the local gas temperature (temperature selected) or the host mass (mass selected). We find that 3 per cent (40 per cent) of mass (temperature)-selected high-z (z = 5.5) DLAs retain signatures of pollution from PopIII stars, independent of the first star model. Such DLAs have low halo mass (< 109.6 M-circle dot), metallicity (< 10(-3) Z(circle dot)) and star formation rate (< 10(-1.5) M-circle dot yr(-1)). Metal abundance ratios of DLAs imprinted in the spectra of quasi-stellar object can be useful as tools to infer the properties of the polluting stellar generation and to constrain the first star mass ranges. Comparing the abundance ratios derived from our simulations to those observed in DLAs at z = 5, we find that most of these DLAs are consistent within errors with PopII star dominated enrichment and strongly disfavour the pollution pattern of very massive first stars (i.e. 100-500 M-circle dot). However, some of them could still result from the pollution of first stars in the mass range [0.1, 100] M-circle dot. In particular, we find that the abundance ratios from SDSS J1202+ 3235 are consistent with those expected from PopIII enrichment dominated by massive (but not extreme) first stars.