  Determination of proton parton distribution functions is present under the
dynamical parton model assumption by applying DGLAP equations with GLR-MQ-ZRS
corrections. We provide two data sets, referred as IMParton16, which are from
two different nonperturbative inputs. One is the naive input of three valence
quarks and the other is the input of three valence quarks with
flavor-asymmetric sea components. Basically, both data sets are compatible with
the experimental measurements at high scale ($Q^2>2$ GeV$^2$). Furthermore, our
analysis shows that the input with flavor-asymmetric sea components better
reproduce the structure functions at high $Q^2$. Generally, the obtained parton
distribution functions, especially the gluon distribution function, are the
good options of inputs for simulations of high energy scattering processes. The
analysis is performed under the fixed-flavor number scheme for $n_f=$ 3, 4, 5.
Both data sets start from very low scales around 0.07 GeV$^2$, where the
nonperturbative input is directly connected to the simple picture of quark
model. These results may shed some lights on the origin of the parton
distributions observed at high $Q^2$.
