We perform pQCD analysis of the existing DIS data for charged leptons with account of corrections up to the NNLO. The parton distributions, value of strong coupling constant, and high-twist terms are extracted and their stability with respect to account of the NNLO corrections is analyzed. All the quantities are generally stable within their experimental errors. Obtained value of the strong coupling constant is α NNLO Perturbative method is a powerful tool of the modern quantum field theory. In particular, analysis of the deep-inelastic-scattering (DIS) data in terms of perturbative QCD (pQCD) allows for quantitative description of this process and extraction of parton distribution functions (PDFs) together with value of the strong coupling constant α s , which can be used for calculation of the cross sections for other processes with hadronic beams and targets and check of self-consistency of Standard Model. The common practice for analysis of such kind is to take into account only the O(α s ) (or next-to-leading-order (NLO)) corrections to the DIS cross sections since the next-to-next-to-leading-order (NNLO) corrections has not been completely calculated yet. Meanwhile the value of α s is rather large at the values of transferred momentum Q typical for existing DIS data and the higher-order (HO) correction may have impact on the value of extracted quantities. Particular important is to know the HO PDFs that is motivated by the need to calculate precise value of the Higgs boson production cross section since these calculations require account of the HO QCD corrections 1 and corresponding HO PDFs as consistent input. Account of the HO corrections is also important for reducing the total uncertainty in the value of α s determined from the existing DIS data since the theoretical error due to HO corrections dominates the error in value of α s obtained in the NLO approximation 2 . Besides, account of HO corrections is necessary for clarification of nature of the 1/Q 2 terms, which are apparent in the NLO analysis of DIS data. These terms correspond to the high-twist (HT) contributions, but also can be simulated by the HO contributions since they are proportional to factors [α s (Q)] n