In this proceedings we discuss a prescription to extract the QCD
strong coupling constant at N^{3}LON3LO
precision in perturbative QCD using a combination of
{{O}}(\alpha_{S}^{3})O(αS3)
calculations in pQCD and estimations of the
{{O}}(\alpha_{S}^{4})O(αS4)
corrections from the data. The method is applied to a set of event shape
averages measured in experiments at the LEP, PETRA, PEP and TRISTAN
colliders. In our analysis we account for hadronization effects with
models from modern Monte Carlo event generators and analytic
hadronization models. We conclude that the precision of the
\alpha_{S}αS
extraction cannot be improved significantly only with pQCD predictions
of higher orders, and further progress in these studies requires a
significant advances in the studies and modeling of hadronization
process.