Density functional theory calculations have been performed at the B3LYP/6-31++G** level to investigate the effect of the isomorphously substituted ZSM-5 zeolite with trivalent elements B, Al, and Ga on the hydrogen atom of the Brönsted acid site via nuclear quadrupole resonance (NQR) and nuclear magnetic resonance (NMR) parameters. ZSM-5 clusters were modeled as ((H) 3 SiO) 3 -Si-O(H)-M-(OSi(H) 3 ) 3 (M = B, Al, and Ga). Variations of isotropic chemical shifts, δ iso , of 1 H atom in the Brönsted acid site as well as quadrupole coupling constant, C Q , and asymmetry parameter, η Q , of 2 H atom in the Brönsted acid site are well correlated with the trend in acidic strength of isomorphously substituted H-ZSM-5, estimated from the experimental results: B-ZSM-5 < Ga-ZSM-5 < Al-ZSM-5.The C Q value of 2 H atom in the Brönsted acid site decreases in the order: B-ZSM-5, Ga-ZSM-5, Al-ZSM-5, respectively; eventually, their η Q values of 2 H atom in the Brönsted acid site increase concluding models with high C Q have a lower Brönsted acidity than those with low C Q . In addition, the calculated chemical shifts of the Brönsted proton increase in the sequence: B-ZSM-5 < Ga-ZSM-5 < Al-ZSM-5. All the results show a high consistency among NQR and NMR parameters and relative acidity of the isomorphously substituted ZSM-5 zeolite.