We study contributions of quark-squark-gluino interactions in Minimal Supersymmetric Standard Model (MSSM) to B → PP (PP = Kπ, ππ,KK) decays using QCD improved factorization method for the evaluation of the hadronic matrix elements and taking into account renormalization group running of the Wilson coefficients from SUSY scale (≈ mq̃) to the low energy scale (= mb) applicable to B decays. Using ...