Molecular complexes formed in the organic phase during solvent extraction may self-assemble as reverse micelles, and therefore induce a supramolecular organization of this phase. In most of the cases, water molecules play an essential role in the organization of this non polar medium. The aim of this work is to investigate the speciation of the aggregates formed in the organic phase during solvent extraction, and especially to assess their stability as a function of the number of water molecules included in their polar core. We have focused on malonamide extractants that have already been investigated experimentally. Different stoichiometries of reverse micelles in the organic phase have been studied by means of classical molecular dynamics simulations. Furthermore, umbrella-sampling molecular dynamics simulations have been used to calculate the equilibrium constant (K°) representing the association/dissociation pathways of water molecules in the aggregates and the corresponding reaction free energies (ΔrG°).