Joining the global demand for the discovery of potent NSAIDs with minimized ulcerogenic effect, new pyrazole clubbed thiazole derivatives 5a-o were designed and synthesized. The new derivatives were initially evaluated for their analgesic activity. Eight compounds 5a, 5c, 5d, 5e, 5f, 5h, 5m, and 5o showed higher activity than Indomethacin (potency = 105–130 % vs. 100 %). Subsequently, they were picked for further evaluation of their anti-inflammatory activity, ulcerogenic liability as well as toxicological studies. Derivatives 5h and 5m showed a potential % edema inhibition after 3 h (79.39 % and 72.12 %, respectively), with a promising safety profile and low ulcer indices (3.80 and 3.20, respectively). The two compounds 5h and 5m were subjected to in vitro COX-1 and COX-2 inhibition assay. The candidate 5h showed nearly equipotent COX-1 inhibition (IC50 = 38.76 nM) compared to the non-selective reference drug Indomethacin (IC50 = 35.72 nM). Compound 5m expressed significant inhibitory activities and a higher COX-2 selectivity index (IC50 = 87.74 nM, SI = 2.05) in comparison with Indomethacin (SI = 0.52), with less selectivity than Celecoxib (SI = 8.31). Simulation docking studies were carried out to gain insights into the binding interaction of compounds 5h and 5m in the vicinity of COX-1 and COX-2 enzymes that illustrated the importance of pyrazole clubbed thiazole core in hydrogen bonding interactions. The thiazole motif of compounds 5h and 5m exhibited a well orientation toward COX-1 Arg120 key residue by hydrogen bonding interactions. Compound 5h revealed an additional arene-cation interaction with Arg120 that could rationalize its superior COX-1 inhibitory activity. Compounds 5h and 5m overlaid the co-crystallized ligand Celecoxib I differently in the active site of COX-2. Compound 5m showed an enhanced accommodation with binding energy of − 6.13 vs. − 1.70 kcal/mol of compounds 5h. The naphthalene ring of compound 5m adopted the Celecoxib I benzene sulfonamide region that is stabilized by hydrogen-arene interactions with the hydrophobic sidechains of the key residues Ser339 and Phe504. Further, the core structure of compound 5m, pyrazole clubbed thiazole, revealed deeper hydrophobic interactions with Ala513, Leu517 and Val509 residues. Finally, a sensitive and accurate UPLC-MS/MS method was developed for the simultaneous estimation of some selected promising pyrazole derivatives in rat plasma. Accordingly, compounds 5h and 5m were suggested to be promising potent analgesic and anti-inflammatory agents with improved safety profiles and a novel COX isozyme modulation activity.