Following previous works on that topic, we consider Euclidean hadronic matrix elements in position space of two spatially separated local currents on the lattice, in order to extract the x dependence of parton distribution functions (PDFs). The corresponding approach is often referred to by the term . In this work we will consider valence quark PDFs of an unpolarized proton. We adapt the previously established formalism to our choice of operators. The calculation of the two-current matrix elements requires the evaluation of four-point functions. The corresponding calculation is carried out on a nf=2+1 gauge ensemble with lattice spacing a=0.0856 fm and pseudoscalar masses mπ=355 MeV and mK=441 MeV. The four-point functions have been evaluated in a previous project. The lattice data is converted to the MS¯ scheme at a scale μ=2 GeV and improved with respect to lattice artifacts. We use a common model as fit ansatz for the lattice data in order to extract the PDFs. Published by the American Physical Society 2024