We show that the QCD factorization approach for B-meson decays to charmless hadronic two-body final states can be extended to include electromagnetic corrections. The presence of electrically charged final-state particles complicates the framework. Nevertheless, the factorization formula takes the same form as in QCD alone, with appropriate generalizations of the definitions of light-cone distribution amplitudes and form factors to include QED effects. More precisely, we factorize QED effects above the strong interaction scale ΛQCD for the non-radiative matrix elements $$ \left\langle {M}_1{M}_2\left|{Q}_i\right|\overline{B}\right\rangle $$
M
1
M
2
Q
i
B
¯
of the current-current operators from the effective weak interactions. The rates of the branching fractions for the infrared-finite observables $$ \overline{B}\to {M}_1{M}_2\left(\gamma \right) $$
B
¯
→
M
1
M
2
γ
with photons of maximal energy ∆E ≪ ΛQCD is then obtained by multiplying with the soft-photon exponentiation factors. We provide first estimates for the various electromagnetic corrections, and in particular quantify their impact on the πK ratios and sum rules that are often used as diagnostics of New Physics.