The free energy of alpha-helix dimerization of the transmembrane (TM) region of glycophorin A was estimated from a 125-ns molecular dynamics (MID) simulation in a membrane mimetic. The free energy profile was obtained by allowing the TM helical segments to diffuse reversibly along the reaction pathway. Partition of the potential of mean force into free energy components illuminates the critical steps of alpha-helix recognition and association. At large separations, the TM segments are pushed together by the solvent, allowing initial, but not necessarily native, interhelical interactions to occur. This early recognition stage precedes the formation of native contacts, which is accompanied by a tilt of the helices, characteristic of the dimeric structure. This step is primarily driven by the van der Waals helix-helix interactions. Free energy perturbation calculations of the L75A and I76A point mutations reveal a disruption in helix-helix association due to a loss of favorable dispersion interactions. Additional MID simulations of the native TM dimer and of a single alpha-helix confirm that, prior to association, individual alpha-helices are independently stable, in agreement with the "two-stage" model of integral membrane protein folding.