Parton distribution functions (PDFs) describe universal properties of bound states and allow us to calculate scattering amplitudes in processes with large momentum transfer. Calculating PDFs involves the evaluation of matrix elements with a Wilson line in a light-cone direction. In contrast to Monte Carlo methods in Euclidean spacetime, these matrix elements can be directly calculated in Minkowski-space using the Hamiltonian formalism. The necessary spatial- and time-evolution can be efficiently applied using established tensor network methods. We present PDFs in the Schwinger model calculated with matrix product states.