In the present study NBO method  contain the HOMO and the LUMO energies are calculated for 10 different heterocycles anticancer drug using B3LYP/6-31G(d,p). Frontier molecular orbitals (HOMO and LUMO) and Molecular Electrostatic Potential map of the compound was produced by using the π stacking of structures and anticancer activity of molecules. The NBO analysis was suggested that the molecular system contains π- π interaction, strong conjugative interactions and the molecule become more polarized owing to the movement of π-electron cloud from donor to acceptor. NBO, HOMO and LUMO energies, were investigated and Anticancer activity of Aromatic Heterocyclic compounds was investigated by NBO study and result was compared with our previous study about NICS and S-NICS of these 10 anticancer drug. the HOMO/LUMO gap of the heterocycle anticancer drug is significantly different from each other. The NBO method is used in both symmetric and asymmetric molecules and provides accurate information on the aromatics of the compound, especially the heterocyclic rings. It also provides accurate information in protected areas. Molecule 8 has the highest amount of HOMO and therefore aromaticity among the studied compounds which confirms the result of molecular orbital examination.

In this study, we have sujessted a theoritical aproach via computing of nucleus-independent chemical shifts (S-NICS) in view point of probes motions in a sphere of de-shielding and shielding spaces of rings of anticancer contain heterocyclic rings. Few research in theoretical of a statistical approach in NMR shielding and nucleus independent chemical shifts for study of anticancer drugs. S-NICS method is an accurate method for estimation the amount of aromaticity in the non-benzene rings such as heterocyclic rings which is a famous for organic chemical compounds as anti-cancer disease. Theoretical of a statistical method of NICS-nucleus independent chemical shift in small distance and are alternate in large distance in the center of heterocyclic rings. In this research, we have explored the statistical approaches through nucleus-independent chemical shifts-SNICS calculations in view point of Bq motions in the center of sphere spaces of heterocyclic rings of some anticancer drugs. The S-NICS method is an accurate computing method for estimation the aromaticity in the non-benzene rings such as heterocyclic rings which are standard molecules in organic anti-cancer compounds. Although NICS values for benzene and naphthalene and so on can be indicated as aromaticity criterion, for other molecules such as heterocyclic rings and their derivatives, S-NICS values are much more accurate compare to NICS index. Hetrocyclic compound with good π stacks are good cacdidate for anticancer medicines.  Precence of Nitrogen, oxygen and sulfur in these compounds make them similar biochemical parts such as nuclice acids and. Aromaticity of heterocyclic ring.   This similarity especially them π stacks give the the ability of interaction with DNAs and RNAs and the heterocyclic compounds.