The B3LYP/6-31G* level of theory was used for the optimization of nonclassical platinum drugs with planar aromatic heterocycle ligands in addition to their monoaqua species, diaqua species and monofunctional adducts. Single point calculations were performed on optimized geometries using the MP2/6-31G** method, and the Lanl2dz pseudo potential for the Pt atom remained constant in all calculations. It can be easily seen that the trans effect can influence both geometrical structures and bond dissociating energies (BDEs). On the basis of our calculation, we get that the long-lived monofunctional adducts with chloride ligand may be in existence, because they possessed higher stability energies, easily dissociated chloride compared with the monoaqua species dissociating chloride and that replacement of the first chloride by bases became exothermic in solution. Our calculated results also demonstrate that the strongest H-bonds appear in the complexes of q-wG and t-wG using different methods. In comparison with adenine complexes, the corresponding guanine complexes possess larger interaction energies as well as higher stability energies either corrected by basis set superposition error (bsse) of Boys-Bernardi counterpoise method or uncorrected both in gas phase and in solution with one exception that the stability energy of q-ClA complex in solution is larger than that of q-ClG complex. Finally, the nature of bond was analyzed in terms of partial charges distribution based on NBO population.