Direct solution of renormalization group equations of QCD in x-space: NLO implementations at leading twist