Recent crystal structures have revealed that regulatory subunit RIalpha of PKA undergoes a dramatic conformational change upon complex formation with the catalytic subunit. Molecular dynamics studies were initiated to elucidate the contributions of intrinsic conformational flexibility and interactions with the catalytic subunit in formation and stabilization of the complex. Simulations of a single RIalpha nucleotide binding domain (NBD), missing cAMP, showed that its C helix spontaneously occupies two distinct conformations: either packed against the nucleotide binding domain as in its cAMP bound structure, or extended into an intermediate form resembling that of the holoenzyme structure. C helix extension was not seen in a simulation of either RIalpha NBD. In a model complex containing both NBDs and the catalytic subunit, well-conserved residues at the interface between the NBDs in the cAMP bound form were found to stabilize the complex through contacts with the catalytic subunit. The model structure is consistent with available experimental data.