We extract polarized parton distribution functions (PPDFs), referred to as “KTA17,” together with the highly correlated strong coupling αs from recent and up-to-date g1 and g2 polarized structure functions world data at next-to-next-to-leading order in perturbative QCD. The stability and reliability of the results are ensured by including nonperturbative target mass corrections as well as higher-twist terms which are particularly important at the large-x region at low Q2. Their role in extracting the PPDFs in the nucleon is studied. Sum rules are discussed and compared with other results from the literature. This analysis is made by means of the Jacobi polynomials expansion technique to the DGLAP evolution. The uncertainties on the observables and on the PPDFs throughout this paper are computed using standard Hessian error propagation which served to provide a more realistic estimate of the PPDFs uncertainties.