This Article Full Text Full Text (PDF) All Versions of this Article: biophysj.104.054460v1 88/5/3243    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Paterlini, M. G. Articles by Thomas, D. D. Search for Related Content PubMed PubMed Citation Articles by Paterlini, M. G. Articles by Thomas, D. D.

The -Helical Propensity of the Cytoplasmic Domain of Phospholamban: A Molecular Dynamics Simulation of the Effect of Phosphorylation and Mutation

M. Germana Paterlini ^* and David D. Thomas 

^* Certusoft, Inc., Minneapolis, Minnesota; and Department of Biochemistry, Molecular Biology, and Biophysics, University of Minnesota, Minneapolis, Minnesota

Correspondence: Address reprint requests to M. Germana Paterlini, Certusoft, Inc., 7900 W. 78th St., Suite 165, Minneapolis, MN 55439. Tel.: 952-921-0351; E-mail: germana{at}certusoft.com.

We have used molecular dynamics simulations to investigate the effect of phosphorylation and mutation on the cytoplasmic domain of phospholamban (PLB), a 52-residue protein that regulates the calcium pump in cardiac muscle. Simulations were carried out in explicit water systems at 300 K for three peptides spanning the first 25 residues of PLB: wild-type (PLB_1–25), PLB_1–25 phosphorylated at Ser16 and PLB_1–25 with the R9C mutation, which is known to cause human heart disease. The unphosphorylated peptide maintains a helical conformation from 3 to 15 throughout a 26-ns simulation, in agreement with spectroscopic data. Comparison with simulations of a fourth peptide truncated at Pro21 showed the importance of the region from 17 to 21 in preventing local unfolding of the helix. The results suggest that residues 11–16 are more likely to unfold when specific capping motifs are not present. It is proposed that protein kinase A exploits the intrinsic flexibility of the 11–21 region when binding PLB. In agreement with available CD and NMR data, the simulations show a decrease in the helical content upon phosphorylation. The phosphorylated peptide is characterized by helix spanning residues 3–11, followed by a turn that optimizes the salt-bridge interaction between the side chains of the phosphorylated Ser-16 and Arg-13. Replacing Arg-9 with Cys results in unfolding of the helix from C9 and an overall decrease of the helical conformation. The simulations show that initiation of unfolding is due to increased solvent accessibility of the backbone atoms near the smaller Cys. It is proposed that the loss of inhibitory potency upon Ser-16 phosphorylation or R9C mutation of PLB is due to a similar mechanism, in which the partial unfolding of the cytoplasmic helix of PLB results in a conformation that interacts with the cytoplasmic domain of the calcium pump to relieve its inhibition.

This article has been cited by other articles:

				L. M. Espinoza-Fonseca, D. Kast, and D. D. Thomas Molecular Dynamics Simulations Reveal a Disorder-to-Order Transition on Phosphorylation of Smooth Muscle Myosin Biophys. J., September 15, 2007; 93(6): 2083 - 2090. [Abstract] [Full Text] [PDF]