Pharmaceutical cocrystals are currently gaining interest among the scientific community, due to their great potential for providing novel crystalline forms with superior properties such as solubility, dissolution rate, bioavailability, and stability. Robust computational tools are valuable tools in the rationalization of cocrystal formation, by providing insight into the intermolecular interactions of multicomponent molecular solids. In this study, various computational techniques based on charge density analysis were implemented to assess structural and energetical perspectives of the interactions responsible for the formation and stability of entacapone-theophylline-water (ETP-THP-water, 1:1:1). Significant non-covalent interactions (NCIs) were identified and evaluated by Hirshfeld surface analysis and density functional theory (DFT) computations, and three-dimensional networks (energy vector diagrams, lattice energy frameworks) were constructed, outlining the crucial stabilizing role of water and the dominance of π-π stacking interactions in the cocrystal. Furthermore, thermal dehydration studies confirmed the strong binding of water molecules in the crystal lattice, as expressed by the high activation energy.
