Protein kinases are essential for the regulation of cellular growth and metabolism. Since their dysfunction leads to debilitating diseases, they represent key targets for pharmaceutical research. The rational design of kinase inhibitors requires an understanding of the determinants of ligand binding to these proteins. In the present study, a theoretical model based on continuum electrostatics and a surface-area-dependent nonpolar term is used to calculate binding affinities of balanol derivatives, H-series inhibitors, and ATP analogues toward the catalytic subunit of cAMP-dependent protein kinase (cAPK or protein kinase A). The calculations reproduce most of the experimental trends and provide insight into the driving forces responsible for binding. Nonpolar interactions are found to govern protein-ligand affinity. Hydrogen bonds represent a negligible contribution, because hydrogen bond formation in the complex requires the desolvation of the interacting partners. However, the binding affinity is decreased if hydrogen-bonding groups of the ligand remain unsatisfied in the complex. The disposition of hydrogen-bonding groups in the ligand is therefore crucial for binding specificity. These observations should be valuable guides in the design of potent and specific kinase inhibitors.