In the present study, nucleophilic properties of adenine and guanine are examined by means of density functional theory. H+ is used as a model electrophile. Two modes of H+ attack on the bases are considered: on the neutral molecule and on the anion. Solvent effects are modeled by means of polarizable continuum model. Regioselectivity of attack is studied by analyzing two contributions. The first one is the energetic ordering of the tautomers. The second is the relative inherent reactivity of nucleophilic sites in the bases. Atomic softnesses calculated by means of charge sensitivity analysis are employed for this purpose. The most reactive sites in various tautomers are identified on the ground of Li–Evans model. For adenine, it is demonstrated that both in basic and in neutral pH N7 atom possesses the most nucleophilic character. In polar solvents, N7 substitution is also most favored energetically. In basic pH and nonpolar solvents as well as in the gas phase, N9 substitution is slightly more probable. For guanine, a mixture of N7- and N9-substituted products can be expected in basic pH. In neutral pH, inherent reactivity and energy trends are opposite to each other; therefore, the substitution does not occur. Experimentally observed products of reactions with various electrophiles and in various conditions confirm the results obtained in this study.