Atomic motions and energetics for a phosphate transfer reaction catalyzed by the cAMP-dependent protein kinase are calculated by plane-wave density functional theory, starting from structures of proteins crystallized in both the reactant conformation (RC) and the transition-state conformation (TC). In TC, we calculate that the reactants and products are nearly isoenergetic with a 20-kJ/mol barrier, whereas phosphate transfer is unfavorable by 120 kJ/mol in the RC, with an even higher barrier. With the protein in TC, the motions involved in reaction are small, with only P␥ and the catalytic proton moving >0.5 Å. Examination of the structures reveals that in the RC the active site cleft is not completely closed and there is insufficient space for the phosphorylated serine residue in the product state. Together, these observations imply that the phosphate transfer reaction occurs rapidly and reversibly in a particular conformation of the protein, and that the reaction can be gated by changes of a few tenths of an angstrom in the catalytic site.protein dynamics ͉ quantum chemistry ͉ reaction kinetics P rotein kinases regulate many biological processes by transferring a phosphate group from adenosine triphosphate (ATP) to the side chains of particular serine, threonine, or tyrosine residues. The bulky, charged phosphate group alters the conformation and function of the target protein (1, 2). Different kinases recognize different primary sequence motifs surrounding the residue to be phosphorylated, in a highly regulated fashion (3-6). Structural studies have revealed several conformational changes, such as closing of the active-site cleft, the packing of the activation loop, and rotation of the C-helix, which are often implicated in controlling the activity of protein kinases (2). The reasons for such control are clear, but no answer has been provided to such questions as these: ''How closed is closed?'' or ''Is this particular conformation of the activation loop 'good enough' for phosphorylation to occur?'' Quantum chemistry is required to objectively answer these questions.The extent of conformational heterogeneity in a covalent protein reaction was first quantified in a series of experiments monitoring the temperature-dependent rebinding of CO to myoglobin after flash photolysis (7). Agmon and Hopfield created a concise phenomenological model describing this situation, using transition-state theory to describe the vibrational reaction, and a diffusive coordinate that describes the protein conformation and modulates the reaction barrier to the vibrational transition (8, 9). Moving beyond the phenomenological model requires specifying both the conformational heterogeneity and the sensitivity of the reaction barrier to this heterogeneity. Careful structural analysis (10) and quantum chemistry calculations (11) showed that the reaction barrier heterogeneity is indeed a reasonable consequence of observed structural heterogeneity at the myoglobin active site. A reevaluation of a wide variety of myoglobin data also shows ...