We report a computation methodology, which leads to the ability to partition the Gibb's free energy for the complexation reaction of aromatic drug molecules with DNA. Using this approach, it is now possible to calculate the absolute values of the energy contributions of various physical factors to the DNA binding process, whose summation gives a value that is reasonably close to the experimentally measured Gibb's free energy of binding. Application of the methodology to binding of various aromatic drugs with DNA provides an answer to the question "What forces are the main contributors to the stabilization of aromatic ligand-DNA complexes?"