Recent breakthrough experiments revealed the iso-reticular Zr-MOFs, MIL-140B and MIL-140C, as promising sorbents for the separation of C6 isomers. Interestingly while the ultra-small pore MIL-140B exhibited hexane isomer sorption hierarchy according to the normal boiling point order (n-C6 > 3MP (3-methyl pentane)), an uncommon shift in the elution order was observed in the larger pore MIL-140C. It was only speculated that the flexibility of the MOFs might be the origin of this intriguing behavior. Herein, flexible force field hybrid osmotic Monte Carlo combined with molecular dynamics simulations were carried out to unravel the microscopic mechanism of the adsorption and dynamics of both C6 isomers in MIL140B and MIL140C. Thermodynamically preferred adsorption of n-C6 over 3MP was predicted for MIL-140B and to a slightly less extent for MIL-140C. Interestingly while the mobility of n-C6 was found to remain higher than that of 3MP in the whole range of loading for MIL-140B, 3MP becomes more mobile than n-C6 at saturation in MIL-140C. This suggests that this kinetics order is most probably the origin of the inversion of the elution order observed experimentally for MIL-140C. The translational and rotational dynamics of the two guests in MIL-140B and MIL-140C was further understood in-depth.