Read one of the shapefiles with the {sf} package and get familiar enough working with geometries in that system, then come back for help in combining to sf data frames, and eventually a for-loop or vectorized function recipe to sweep them all into a single object. The hard part will probably in the first step to sweep in both the point geometries and the associated concentrations.