Enantioselectivity of the propranolol on β-cyclodextrin was simulated by molecular modeling. Monte Carlo (MC) docking and molecular dynamics (MD) simulations were applied to investigate the molecular mechanism of enantioselective difference of both enantiomeric complexes. An energetic analysis of MC docking simulations coupled to the MD simulations successfully explains the experimental elution order of propranolol enantiomers. Molecular dynamics simulations indicate that average energy difference between the enantiomeric complexes, frequently used as a measure of chiral recognition, depends on the length of the simulation time. We found that, only in case of much longer MD simulations, noticeable chiral separation was observed.