The [FeFe]-hydrogenases in the green alga Chlamydomonas reinhardtii utilize photogenerated electrons to reduce protons into hydrogen gas. The electrons are supplied from photosystem I and transferred to the [FeFe]hydrogenase through specific hydrogenase-ferredoxin association. To understand how structural and kinetic factors control the association better, we used Brownian dynamics simulation methods to simulate the charge-transfer complex formation between both native and in silico mutants of the [FeFe]-hydrogenase HYDA2 and the [2Fe2S]-ferredoxin FDX1 from C. reinhardtii. The chances in binding free energy between different HYDA2 mutants and the native FDX1 were calculated by the free-energy perturbation method. Within the limits of our current models, we found that two HYDA2 mutations, T99K(H) and D102K(H), led to lower binding free energies and higher association rate with FDX1 and are thus promising targets for improving hydrogen production rates in engineered organisms.