Based on recently published initial experimental results on the intercalation of a class of broad spectrum antiparasitic compounds, we present a purely theoretical approach for determining if these compounds may preferentially intercalate with guanosine/cytosine (GC)-rich or adenosine/thymidine (TA)-rich regions of DNA. The predictive model presented herein is based upon utilization of density functional theory (DFT) to determine a priori how the best intercalator may energetically and sterically interact with each of the nucleoside base pairs. A potential new method using electrostatic potential maps (EPMs) to visually select the best poses is introduced and compared to the existing brute-force center of mass (COM) approach. The EPM and COM predictions are in agreement with each other, but the EPM method is potentially much more efficient. We report that 4-azatryptantrin, the best intercalator, is predicted to favor π-stacking with GC over that of TA by approximately 2-4 kcal/mol. This represents a significant difference if one takes into account the Boltzmann distribution at physiological temperature. This theoretical method will be utilized to guide future experimental studies on the elucidation of possible mechanism(s) for the action of these antiparasitic compounds at the molecular level.