The kinetics of charge recombination between the primary photoxidized donor (P(+)) and the secondary reduced quinone acceptor (Q(B)(-)) have been studied in reaction centers (RCs) from the purple photosynthetic bacterium Rhodobacter sphaeroides incorporated into lecithin vesicles containing large ubiquinone pools over the temperature range 275 K </= T </= 307 K. To account for the non-exponential kinetics of P(+) re-reduction observed following a flash, a new approach has been developed, based on the following assumptions: 1) the exchange of quinone between different vesicles is negligible; 2) the exchange of quinone between the Q(B) site of the RC and the quinone pool within each single vesicle is faster than the return of the electron from the primary reduced acceptor Q(A)(-) to P(+); 3) the size polydispersity of proteoliposomes and the distribution of quinone molecules among them result in a quinone concentration distribution function, P(Q). The first and second moments of P(Q) have been evaluated from the size distribution of proteoliposomes probed by quasi-elastic light scattering (mean radius, <R> = (50 +/- 15) nm). Following these premises, we describe the kinetics of P(+)Q(B)(-) recombination with a truncated cumulant expansion and relate it to P(Q) and to the free energy changes for Q(A)(-)Q(B) --> Q(A)Q(B)(-) electron transfer (DeltaG(AB)(o)) and for quinone binding (DeltaG(bind)(o)) at Q(B). The model accounts well for the temperature and quinone dependence of the charge recombination kinetics, yielding DeltaG(AB)(o) = -7.67 +/- 0.05 kJ mol(-1) and DeltaG(bind)(o) = -14.6 +/- 0.6 kJ mol(-1) at 298 K.