Molecular determinants of thapsigargin binding by SERCA Ca2+-ATPase: a computational docking study

Thapsigargin (TG) is a potent and commonly used inhibitor of the ion transport activity of sarco/endoplasmic reticulum Ca2+-ATPases (SERCA). Based on the recently published crystal structures of rabbit muscle SERCA1a in the Ca2+/E1 (E1) and TG/E2 (E2) conformations, we performed computational docking studies to characterize the molecular interactions that govern binding of TG and TG-analogs by the enzyme. Using the program GOLD (genetic optimization for ligand docking) in combination with the scoring function ChemScore, TG was docked into the binding site of the E1 and E2 conformations of SERCA1a. The docking results revealed a consensus ligand-binding mode consistent with the crystal structure and showed that hydrophobic interactions are the primary driving force of TG binding by SERCA. Moreover, it was shown that the conformational changes accompanying the E2 to E1 transition in the enzyme likely displace TG from its favored orientation in the binding site, thereby substantially reducing its binding affinity. This finding illustrates on the molecular level how TG may exert its inhibitory effect in binding tightly to the E2 form and preventing it from converting into its E1 form, a requirement for catalytic function. We also docked 9 TG analogs into the E2 conformation of the enzyme. Eight of the analogs adopted a binding mode very similar to that of TG, whereas one compound preferred a different orientation in the binding site. Analysis of the predicted binding affinities showed a good correlation with the experimentally observed inhibitory potencies of the analogs. Docking was also performed with several modeled mutants of SERCA1a, whose phenylalanine residue in position 256 (Phe256) had been modified. The experimentally observed declines in TG sensitivity in most of the Phe256 mutants was qualitatively accounted for and appears, at least in part, be due to a slightly altered TG-binding mode.