In the heart, the hERG voltage-gated potassium channel mediates the IKr current, which is crucial for the duration of cardiac action potential. Undesired block of the channel may prolong the QT interval with increased risk of malignant ventricular arrhythmia called torsades de pointes. Although the molecular determinants of hERG block are intensively studied, stereoselectivity has been poorly investigated. Levo-(S)-bupivacaine was the first drug reported to have higher affinity for hERG than its enantiomer. This study aims at understanding the principles underlying the stereoselectivity of bupivacaine block with the help of molecular modeling. Putative binding modes of levo-(S)- and dextro-(R)-bupivacaine inside an open form model of hERG channel were predicted by docking simulations, allowing a clear depiction of ligand-protein interactions. Estimated binding energies for both enantiomers to wild-type channel are in line with previously published electrophysiology measurements. These results may be considered as a confirmation at the molecular level of bupivacaine stereoselective binding towards hERG. Moreover this information lays the foundations for a structural guideline to filter out potentially cardiotoxic drug candidates in silico.