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 .