We conducted research on the adsorption of carmustine (BCNU) on the surface of pristine, N-doped, and BN-doped Hexa-peri-hexabenzocoronene (HBC) nanographene using density functional theory calculations. We used the B3LYP functional in conjunction with 6-31+G(d) and 6-311+G(d,p) basis sets to perform geometrical optimization calculations, vibrational frequencies, and natural bond orbitals (NBO) analysis. The dispersion correction term of Grimme-D3 was incorporated to account for the dispersion interactions. We found that the geometry of the nanographenes remains unchanged after the drug molecule adsorption. The calculated adsorption energies (Eads) were −17.1 kcal.mol−1 for HBC-BCNU, −17.0 kcal.mol−1 for NHBC-BCNU, and −17.9 kcal.mol−1 for BNHBC-BCNU when using the 6-31+G(d) basis set. The corresponding values calculated employing the 6-311+G(d,p) basis set were −17.9, −17.3, and −19.1 kcal.mol−1, respectively. Analysis of frontier molecular orbitals (HOMO and LUMO) showed that the electronic properties of NHBC were significantly affected by the presence of BCNU. We used NBO analysis to investigate the distribution of partial atomic charges in different systems. Finally, using the transition state equation and the smaller basis set, we found that the required time for BCNU to release from HBC, NHBC, and BNHBC is approximately 3.3, 1.8, and 13.0 s, respectively. The obtained values are 13.0, 4.6, and 96.0 s, respectively, for the higher level of calculations. The recovery time for all nanographenes is within a reasonable range, while NHBC exhibits the highest sensitivity for drug detection.