In an aim to probe the structure-function relationship of prostacyclin synthase (PGIS), resonance Raman (RR) spectroscopy and molecular dynamic (MD) simulation approaches have been exploited to characterize the heme conformation and heme-protein matrix interactions for human PGIS (hPGIS) and zebrafish PGIS (zPGIS) in the presence and absence of ligands. The high-frequency RR (1300-1700 cm(-1)) indicates that the heme group is in the ferric, six-coordinate, low-spin state for both resting and ligand-bound hPGIS/zPGIS. The low-frequency RR (300-500 cm(-1)) and MD simulation reveal a salient difference in propionate-protein matrix interactions between hPGIS and zPGIS, as evident by a predominant propionate bending vibration at 386 cm(-1) in resting hPGIS, but two vibrations near 370 and 387 cm(-1) in resting zPGIS. Upon binding of a substrate analogue (U46619, U51605, or U44069), both hPGIS and zPGIS induce a distinctive perturbation of the propionate-protein matrix interactions, resulting in similar Raman shifts to similar to 381 cm(-1). On the contrary, the bending vibration remains unchanged upon binding of inhibitor/ligand (minoxidil, clotrimazole, or miconazole), indicating that these inhibitors/ligands do not interfere with the propionate-protein matrix interactions. These results, together with subtle changes in vinyl bending modes, demonstrate drastically different RR shifts with heme conformational changes in both hPGIS and zPGIS upon different ligand bindings, suggesting that PGIS exhibits a ligand-specific heme conformational change to accommodate the substrate binding. This substrate-induced modulation of the heme conformation may confer high product fidelity upon PGIS catalysis.