Molecular Docking Analysis of Heparin–Diclofenac Complexes: Insights into Enhanced Cox Enzyme Inhibition for Pain Management

The aim of this study was to investigate the molecular interactions of heparin, diclofenac, and their supramolecular complexes with cyclooxygenase enzymes (COX-1 and COX-2) using computational docking techniques. Diclofenac is a widely used nonsteroidal anti-inflammatory drug (NSAID) that inhibits COX isoforms, whereas heparin is a polyanionic glycosaminoglycan with established anticoagulant and emerging anti-inflammatory properties. Supramolecular association between these agents may modulate their physicochemical behavior and target engagement. Molecular modeling, dual-drug docking, and molecular dynamics (MD) simulations were employed to characterize the interactions of heparin, diclofenac, and pre-formed heparin-diclofenac complexes with COX-1 and COX-2. Geometry optimization and lipophilicity (logP) estimates were obtained using HyperChem, while protein-ligand docking was performed in HEX using crystallographic COX structures from the Protein Data Bank. Docking poses were analyzed in Chimera, and selected complexes were refined through short MD simulations. Pre-formed heparin-diclofenac assemblies exhibited markedly enhanced docking scores toward both COX isoforms compared with single ligands. Binding orientation strongly influenced affinity: for COX-1, the heparin-diclofenac configuration yielded the most favorable interaction, whereas for COX-2 the diclofenac-heparin configuration was preferred. Both assemblies adopted binding modes distinct from free diclofenac, suggesting cooperative electrostatic and hydrophobic contacts at the enzyme surface. Supramolecular complexation also altered calculated logP values relative to the individual compounds. MD simulations supported the relative stability of the top-ranked complex-COX assemblies. These findings indicate that heparin-diclofenac assemblies may enhance and reorganize predicted COX interactions in a configuration-dependent manner and illustrate the utility of dual-drug docking for modeling potential synergistic effects. Such insights may inform the design of localized or topical formulations, potentially incorporating non-anticoagulant heparin derivatives, to achieve effective COX inhibition with reduced systemic exposure. However, the results rely on simplified heparin fragments, legacy docking tools, and short MD simulations, and should therefore be interpreted qualitatively. Experimental studies will be essential to confirm whether such supramolecular assemblies form under physiological conditions and whether they influence COX inhibition in vivo.