Structural Changes in the Cytoplasmic Domain of Phospholamban by Phosphorylation at Ser16: A Molecular Dynamics Study

Abstract

Phospholamban is a 52-residue integral membrane protein that regulates the activity of the sarcoplasmic reticulum calcium pump in cardiac muscle. Its inhibitory action is relieved when phospholamban is phosphorylated at Ser16 by cAMP-dependent protein kinase. To computationally explore all possible conformations of the phosphorylated form, and thereby to understand the structural effects of phosphorylation, replica-exchange molecular dynamics (REMD) was applied to the cytoplasmic domain that includes Ser16. The simulations showed that (i) without phosphorylation, the region from Lys3 to Ser16 takes all α-helical conformations; (ii) when phosphorylated, the α-helix is partially unwound in the C-terminal part (from Ser10 to Ala15) resulting in less extended conformations; (iii) the phosphate at Ser16 forms salt bridges with Arg9, Arg13, and/or Arg14; and (iv) the salt bridges with Arg13 and Arg14 distort the α-helix and induce unwinding of the C-terminal part. We then applied conventional all-atom molecular dynamics simulations to the full-length phospholamban in the phospholipid bilayer. The results were consistent with those obtained with REMD simulations, suggesting that the transmembrane part of phospholamban and the lipid bilayer itself have only minor effects on the conformational changes in the cytoplasmic domain. The distortions caused by the salt bridges involving the phosphate at Ser16 readily explain the relief of the inhibitory effect of phospholamban by phosphorylation, as they will substantially reduce the population of all helical conformations, which are presumably required for the binding to the calcium pump. This will also be the mechanism for releasing the phosphorylated phospholamban from kinase.