Adsorption of propranolol enantiomers on naturally chiral copper (Cu(3,1,17)S) and achiral copper (Cu(100)) surfaces were studied by molecular dynamics simulation to unravel the features of adsorbate-adsorbent enantioselectivity. Adsorption of S- and R-propranolol on Cu(3,1,17)S terraces (with 100 plane) leads mainly to endo- and exo-conformers, respectively. Simulated pair correlation function (g(r)) and mean square displacement (MSD) were analyzed to identify adsorption sites of enantiomers on Cu(3,1,17)S substrate surface, and their simulated binding energies were used to access the adsorption strength. According to (g(r)), R-propranolol adsorbs via naphtyl group while S-propranolol mainly adsorbs through chain group. R-enantiomer binds more tightly to the chiral substrate surface than S-enantiomer as indicated by a higher simulated binding energy by 2.74kJmol−1 per molecule. The difference in binding energies of propranolol enantiomers on naturally chiral Cu(3,1,17)S is almost six times larger than on the achiral Cu(100) surface, which substantiates the appreciably strong specific enantioselective adsorption on the former surface.