Regiospecific replacement of a single phosphate (PO2) by a dithiophosphate (PS2) group in an RNA can dramatically increase its binding affinity for a target protein. Thus, complexes between antithrombin and anti-VEGF RNA aptamers with single dithiophosphate moieties and thrombin and VEGF, respectively, display equilibrium dissociation constants K-D of ca. 1 pM, 1000-fold tighter than the native RNA complexes (ca. 1 nM). Inspection of crystal structures of the native and PS2-RNA aptamer:thrombin complexes reveals an RNA-induced fit in the latter. This leads to a close approach between PS2 and the phenyl ring edge of Phe-232 that is surrounded by pairs of lysines and arginines. To better understand the origins of the tighter binding and individual contributions to the interaction energy, we carried out QM calculations with phosphate- and dithiophosphate-benzene and dimethyl phosphate- and dimethyl dithiophosphate-benzene model systems. These calculations demonstrate that the dithiophosphate-benzene interaction is much stronger than the corresponding interaction with phosphate. QM/MM calculations with the full complexes confirmed this finding and support the hypothesis that the electric field generated by basic residues surrounding Phe-232 is key to the polarization of the PS2 moiety. Thus, disparate polarization and dispersion energies between the PO2 and PS2 complexes contribute critically to the difference in binding affinity. By comparison, easier desolvation of the dithiophosphate group compared to phosphate does not contribute decisively to the observed difference in binding affinity. Favorable polarization and dispersion energies may be a general feature of the dramatic affinity gains seen for complexes between RNAs carrying dithiophosphate groups and their binding proteins.