An Allosteric Mechanism Inferred from Molecular Dynamics Simulations on Phospholamban Pentamer in Lipid Membranes

Abstract

Phospholamban functions as a regulator of Ca2+ concentration of cardiac muscle cells by triggering the bioactivity of sarcoplasmic reticulum Ca2+-ATPase. In order to understand its dynamic mechanism in the environment of bilayer surroundings, we performed long time-scale molecular dynamic simulations based on the high-resolution NMR structure of phospholamban pentamer. It was observed from the molecular dynamics trajectory analyses that the conformational transitions between the “bellflower” and “pinwheel” modes were detected for phospholamban. Particularly, the two modes became quite similar to each other after phospholamban was phosphorylated at Ser16. Based on these findings, an allosteric mechanism was proposed to elucidate the dynamic process of phospholamban interacting with Ca2+-ATPase.

About the Authors:

Peng Lian

Affiliations College of Life Science and Biotechnology and Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, Shanghai, China, Key Laboratory of Systems Biomedicine, Ministry of Education, Shanghai Center for Systems Biomedicine, Shanghai Jiao Tong University, Shanghai, China

Dong-Qing Wei

* E-mail: [email protected] (DQW); [email protected] (JFW)

Affiliations College of Life Science and Biotechnology and Key Laboratory of Microbial Metabolism, Shanghai Jiao Tong University, Shanghai, China, Gordon Life Science Institute, San Diego, California, United States of America

Jing-Fang Wang

* E-mail: [email protected] (DQW); [email protected] (JFW)

Affiliations Key Laboratory of Systems Biomedicine, Ministry of Education, Shanghai Center for Systems Biomedicine, Shanghai Jiao Tong University, Shanghai, China, Shanghai Center for Bioinformation and Technology, Shanghai, China

Kuo-Chen Chou

Affiliation: Gordon Life Science Institute, San Diego, California, United States of America

Introduction

Phospholamban (PLN) is a 52-residue integral membrane protein that functions as a regulator of Ca2+ concentration of cardiac muscle cells by triggering the bioactivity of sarcoplasmic reticulum Ca2+-ATPase (SERCA) [1], [2]. Binding of PLN to SERCA can inhibit the SERCA Ca2+ pump activity so as to decrease the Ca2+ concentration within the cell, further leading the reduction of the relaxation rate of heart muscle cells [3], [4]. However, the inhibition will lost after PLN is phosphorylated at position 16 or 17 by cAMP- or calmodulin-dependent protein kinases, resulting in the dissolution of PLN/SERCA complex or an altered interaction between the two proteins [5], [6]. Based on a link between PLN mutations and heart failure in humans, it is found that the alteration of PLN inhibitory function can lead to degenerative cardiomyopathy [7], [8]. Thus, for its regulation of heart contractions in heart muscle cells, PLN has been of considerable interest as a potential target for the treatment of degenerative cardiac diseases.

With an aim of understanding the inhibitory mechanism of PLN, many studies have begun with the structural behaviors of PLN in the membrane or other solvents without SERCA. As reported, more than 75% of PLN in the membranes adopt the pentameric form [9], [10]. However, most experimental evidences support the fact that the inhibition of SERCA is primarily involved in the monomeric form of PLN rather than the pentameric form [11], [12], [13]. Thus, most of the theoretical and experimental studies of PLN are focused on the monomeric form of PLN [14], [15], [16], [17], [18], [19], [20]. It is noted that the monomeric PLN has three distinct structural domains: a short cytoplasmic (CP) helices (also called domain Ia or cytoplasmic domain, residues 1–16), a hinge with a β-turn type III conformation (residues 17–22), as well as a long hydrophobic transmembrane (TM) helix that is composed of domain Ib (residues 23–30) and domain II (residues 31–52) [21], [22]. Further investigations with multidimensional solid-state NMR and hybrid solution NMR give an indication that the CP and TM helices adopt angles of 93–102° and 22–24° respectively with respect to the lipid normal [17], [22]. Additionally, the pentameric PLN is reported to be able to form an ion channel for Ca2+ and Cl- [23], [24]. Also, the pentameric form is considered to be capable of storing the monomeric form, revealing a mechanism for the cell to control inhibition of Ca2+-ATPase [25].

By now, the bellflower and pinwheel models are of general acceptance to be the structural models for pentameric PLN. The former is a high-resolution NMR structure determined by James J. Chou and his co-workers with PDB entry 1zll [26], [27]. In the bellflower model, PLN pentamer shows a pore-forming coiled-coil structure with the TM helices remarkably bending away from the channel pore near the cytoplasmic side. The latter is a theoretical model (PDB entry 1xnu) obtained from the Forster resonance energy transfer measurements [28]. In this model, TM helices are less bended compared with the bellflower model, and the cytoplasmic helices are located in the plane of the membranes.

Although several MD studies of partial or full length PLN in different environments have been reported [15], [18], [19], [20], [29], [30], [31], most of these studies are only focused on the monomeric PLN. So far only two papers [30], [31] are involved with the pentameric PLN. However, the one published by Kim et al. [30] only used the bellflower models, while the other [31] did not consider the phosphorylated PLN at all. In view of this, the present study was initiated to use all-atom molecular dynamics (MD) simulations to study the conformational dynamics of both bellflower and pinwheel models of PLN pentamer as well as their phosphorylated forms in POPC bilayer surroundings, so as to gain the atomistic insights into the internal motion [32], [33] of PLN pentamer and its functions. Similar approaches have been successfully used to investigate various problems related to protein structure [34], [35], [36], [37], [38], [39], protein-protein interactions and protein-drug interactions [40], [41], [42], [43], as well as drug design [44], [45], [46], [47].

Materials and Methods

The initial coordinates for the bellflower and pinwheel models of PLN pentamer were taken from the RCSB Protein Data Bank [48] with PDB entries 1zll [26] and 1xnu [28], respectively. Based on the atomic coordinates of the bellflower structure (1zll), the following four initial systems were generated for MD simulation studies: (i) monomeric PLN as denoted by PLNbf-mono; (ii) pentameric PLN denoted by PLNbf-penta; (iii) phosphorylated PLN monomer denoted by pPLNbf-mono; (iv) phosphorylated PLN pentamer denoted by pPLNbf-penta. Likewise, based on the atomic coordinates of the pinwheel structure (1xnu), the following four initial systems were generated: (i) monomeric PLN as denoted by PLNpw-mono; (ii) pentameric PLN denoted by PLNpw-penta; (iii) phosphorylated PLN monomer denoted by pPLNpw-mono; (iv) phosphorylated PLN pentamer denoted by pPLNpw-penta.

Before performing MD simulations, the aforementioned eight systems were merged with a patch of a pre-equilibrated membrane with 256 POPC bilayers, and subsequently solvated in TIP3P water molecules [49]. Meanwhile, sodium ions were randomly placed to neutralize the simulated systems. Subsequently, 150,000-step energy minimization was performed on water molecules and ions for the monomeric systems using conjugate gradient and line search algorithm to remove energetically unfavorable contacts. Another 50,000 steps were added for the pentameric systems because of their larger sizes. In this procedure, both of the proteins and lipids were fixed. Then, all the systems were heated to 310 K over 500 ps by velocity rescaling, and equilibrated for 4 ns in the NPT ensemble with all non-hydrogen atoms of the protein and POPC membranes harmonically restrained. The equilibration was continued for another 10 ns with the harmonic restrain applied to the protein backbone only. Finally, the systems obtained at the end of the aforementioned equilibrations were used for the further unrestrained 40-ns MD simulations. As the simulated systems may be far too large to expect the statistical convergence in 40-ns MD simulations, we have done some repeats with different starting structures and starting vectors. Thus, for each simulation, at least 10 MD trajectories with different starting vectors were generated.

All the MD simulations involved in this study were performed by NADM 2.6 package [50] with the CHARMM27 force field parameters [51], periodic boundary conditions, and NPT ensemble. SHAKE algorithm with a tolerance of 10−6 Å was used to constrain all bonds involved with hydrogen atoms. The Nosé-Hoover Langevin piston pressure control and Langevin damping dynamics were employed to keep the simulated systems at a constant pressure of 1 atmosphere and a constant temperature of 310 K [52], [53]. The long-range electrostatic interactions were evaluated by particle mesh Ewald (PME) method [54] with the size of the grid at about 1Å. The van der Waals interactions were treated by using a cutoff of 12Å. All the simulations were performed with a time step of 1fs, and the coordinates were saved every 1ps.

Results and Discussion

Global properties

As a crucial criterion for the convergence measure of the protein systems concerned, the root mean square (RMS) deviation from the initial structures of the Cα atoms for both bellflower and pinwheel models of pentameric and monomeric PLN and pPLN systems were calculated. Both bellflower and pinwheel models of pentameric PLN and pPLN systems underwent various conformational changes during our 40-ns MD simulations. For the bellflower models, the overall RMS deviation value for PLN was around 7.52±0.68 Å, while that for pPLN was 7.50±1.21 Å. For the pinwheel model, the overall RMS deviation value for PLN was 8.29±1.45 Å, while that for pPLN was 9.67±1.61 Å. There was a difference of about 1 Å for the overall RMS deviations between the two models in PLN, indicating that the pinwheel model was more flexible, and involved with larger conformational changes. It is interesting to note that the phosphorylation at Ser16 remarkably enhanced this kind of difference to around 1.4 Å. In view of this, it is very likely that the phosphorylation could play a role in increasing the flexibility of both structural models.

However, such RMS differences are really small for a biomacromolecular system. One possible explanation for this observation is that most residues of PLN have little conformational changes. As aforementioned, PLN monomer has three distinct structural domains: a short cytoplasmic domain (domain Ia), a hinge with a β-turn type III conformation, and a long hydrophobic transmembrane domain (domain Ib and II). The long transmembrane domain, which occupies almost 60% residues of PLN monomer, is inserted into the lipid bilayer with little room for large amplitude fluctuations. Therefore, the transmembrane domain has little contributions to the RMS differences. In other words, the major contributions to the RMS differences were from the cytoplasmic and the β-turn, which contain much less residues than the transmembrane domain. To further confirm this, we calculated the RMS fluctuation values for each residue in both bellflower and pinwheel models of pentameric PLN and pPLN systems. The results are shown in Fig. 1A, from which we can see that for both the bellflower and pinwheel models, the RMS fluctuation values in the transmembrane regions were almost the same. However, the diversity was detected at the position of res.10–20 that was situated outside of the bilayers, belonging to the cytoplasmic domain where the phosphorylation occurred.

[Figure omitted. See PDF.]

Figure 1. Global properties for the MD simulations.

(A) Time dependent RMS fluctuations of all the Cα atoms, (B) Definitions of the angle θ and distance d, (C) Time dependent angle θ value, and (D) Time dependent distance d value. The black, green, red, and blue curves stand for pinwheel PLN, bellflower PLN, pinwheel pPLN, and bellflower pPLN, respectively.

https://doi.org/10.1371/journal.pone.0018587.g001

Structure and dynamics of the cytoplasmic domain

In order to study the structural and dynamic diversity of the cytoplasmic domain, in Fig. 1B we defined an angle θ to describe the direction of the cytoplasmic domain to the POPC bilayers, and a distance d to describe the distance of the cytoplasmic domain with the POPC bilayers. Based on the 10 independent simulations for each simulated system (40-ns MD trajectories), time averages were used to compute expectation values for the angle θ and distance d. Following similar analyses by some theoretical studies [55], [56], bootstrap-style sampling-with-replacement approach was adopted to generate new re-sampled datasets for each simulation trajectory, each of which was used to estimate average values for the angle θ and distance d. The distribution of these average values over all re-sampled dataset for all the MD trajectories was subsequently employed to estimate the variability of the angle θ and distance d, and also to calculate variances and confidence intervals.

The angle θ was considered as the angle of the cytoplasmic domain with the lipid bilayer. The first and last residues in the cytoplasmic domain were used to determine a line, which was further used to determine the angle θ. With the development of simulation along time, the θ value for the bellflower model was about 41.19±8.46 degrees, while that for the pinwheel model was about 23.40±7.46 degrees (Fig. 1C). The differences in the θ value indicated that bellflower and pinwheel structures assumed quite different conformations in the cytoplasmic domain: the former was extended to the region away from the bilayers, while the latter bended to near the bilayers. However, after phosphorylation at Ser16, both the two models shared a similar θ value (about 30 degrees,), indicating that the two were with quite the same configurations due to the phosphorylation effect.

It can be seen from Fig. 1D that the d value for the bellflower model was about 30.00±1.19 Å, while that for the pinwheel model was about 33.72±2.40 Å. After phosphorylation at Ser16, the distance d for both models increased: for the bellflower structure, it was increased from 30 Å to 32 Å; for the pinwheel mode, it was increased from 32 Å to 37 Å. These evidences further indicate that phosphorylation at Ser16 could enhance the flexibility of PLN as mentioned above. This finding is also quite consistent with the NMR study in DPC micelles [14] and previous MD simulations [18], [19], [20], [30].

Structural transition in the cytoplasmic domain

The structural differences in the bellflower and pinwheel models are mainly located in the cytoplasmic domain. In general, the bellflower structure assumes an extended appearance; while the pinwheel model assumes a bent conformation having more interactions with the bilayers (Fig. 2). The cytoplasmic helices in the bellflower assembly are oriented at ∼20 degrees with respect to the bilayers, while the ones in the pinwheel assembly are ∼90 degrees. It is intriguing to have observed the structural transitions between the two models via our MD simulations. As shown in Fig. 3, for the bellflower model, although the initial structure was in the “extended” state (θ value is 57.43 degrees, and d value is 32.23 Å), after 21-ns simulations it became the “bent” state (θ value is 24.95 degrees, and d value is 31.86 Å) as shown in the pinwheel model. Likewise, for the pinwheel model, similar transition was detected after 19-ns simulations when it changed its original “bent” conformation (θ value is 16.91 degrees, and d value is 42.28 Å) to the “extended” one (θ value is 39.34 degrees, and d value is 30.57 Å). These findings were also supported by the experimental results that both extended and bent states of PLN have been readily detected in the micelles and lipid bilayers using electron paramagnetic resonance (EPR) spectroscopy [57], [58]. Further analyses of chemical shift perturbation for the PLN-SERCA complex showed clear evidences for the structural transition between the extended and bent states [57], [59].

[Figure omitted. See PDF.]

Figure 2. The crystal structures of the bellflower and pinwheel models for the pentameric PLN.

(A) In the bellflower model, PLN assumes “R”-state configuration, and (B) in the pinwheel model, PLN assumes “L”-state configuration.

https://doi.org/10.1371/journal.pone.0018587.g002

[Figure omitted. See PDF.]

Figure 3. Structural transitions between the “R” and “L” states.

In this snapshot figure, only one subunit of the pentamer is shown for the bellflower model (green) and pinwheel model (red), respectively.

https://doi.org/10.1371/journal.pone.0018587.g003

In the pPLN systems, the conformations from both the bellflower structure and pinwheel structure look almost the same. As shown in Fig.1C, after 20-ns MD simulations the angle θ values for both pPLN systems are 30.87±9.93 degrees (bellflower model) and 27.37±9.05 degrees (pinwheel model). One possible explanation for this observation is that phosphorylation can increase the interactions of PLN with the lipid bilayer, resulting in the fact that the cytoplasmic domain no longer moves as freely as the unphosphorylated structure. Thus, after phosphorylation at Ser16, both models have to adopt a similar configuration, which is also supported by the previous theoretical studies [19], [20].

Accordingly, our theoretical studies supported the allosteric model to elucidate the interactions between PLN and SERCA, as illustrated in Fig. 4. According to the allosteric model, there is a dynamic equilibrium between the “extended” and “bent” conformations for PLN in the free form. The major difference between the extended and bent states for PLN is located on the cytoplasmic domain. For the extended state of PLN (bellflower model), the cytoplasmic domain is in a dynamically disordered state in which the central part of this domain is unfolded. For the bent state of PLN (pinwheel model), this domain employs an ordered state. Both models can interact with SERCA, forming transition and active states respectively. Only the bellflower structure with the “extended” conformation can form a thermodynamically more stable complex (active state) with SERCA so as to inhibit the activity of Ca2+ pump. This is because the bent state of PLN is not propitious for the SERCA binding by the reason that the cytoplasmic domain in this state is much close to the lipid, far away from SERCA. Owing to the disordered structure, the extended state of PLN can be extended above the lipid surface and interact with SERCA.

[Figure omitted. See PDF.]

Figure 4. A proposed allosteric model for the regulation mechanism of SERCA by PLN.

In this figure, the membrane and PLN are represented by rectangles and colored in gray and blue, respectively. SERCA is shown in blue boxes. The arrows in black show the formation of SERCA-PLN complex. The gray arrows are for the phosphorylation and dephosphorylation. Reactions shown with red arrows are the processes we found in this study. Phosphoric group is represented by the red dotted circle.

https://doi.org/10.1371/journal.pone.0018587.g004

However the phosphorylation at Ser16 reduces the structural difference between the PLN monomer in bellflower model and that in pinwheel model, revealing a structural transition between the ordered and disordered states in the cytoplasmic domain. This order-to-disorder or disorder-to-order transition can influence the conformational equilibrium between the bent and extended states of PLN monomers. Specifically, phosphorylation can shift the conformational equilibrium towards the extended state via a cooperative effect [60]. Furthermore, the phosphorylation at Ser16 can also inactive the SERCA-PLN complex, independent of it bding in either the transition state or in the active state, by disrupting the interactions between SERCA and PLN. According to the experimental and theoretical studies, a crucial hydrogen bond formed by PLN (Gln26) and SERCA (Arg234) maybe lost. In this case, the activity of Ca2+ pump is recovered. Interestingly, the dephosphorylation of Ser16 may drive PLN to the original equilibrium again so as to make the regulation cycle move on. Using the allosteric model, our findings show for the first time how the two different conformations of PLN monomer are produced. The results of this study also reveal the mechanism of SERCA Ca2+ pump's regulation by PLN.

Acknowledgments

The authors are very much indebted to Prof. Zhi-Feng Shao, Dr. Daniel M. Czajkowsky, as well as the reviewer for their constructive comments, which are very helpful for strengthening the presentation of this paper.

Author Contributions

Conceived and designed the experiments: J-FW. Performed the experiments: PL. Analyzed the data: PL D-QW. Contributed reagents/materials/analysis tools: D-QW K-CC. Wrote the paper: J-FW K-CC.

Citation: Lian P, Wei D-Q, Wang J-F, Chou K-C (2011) An Allosteric Mechanism Inferred from Molecular Dynamics Simulations on Phospholamban Pentamer in Lipid Membranes. PLoS ONE6(4): e18587. https://doi.org/10.1371/journal.pone.0018587

References

1. Zamoon J, Mascioni A, Thomas DD, Veglia G (2003) NMR solution structure and topological orientation of monomeric phospholamban in dodecylphosphocholine micelles. Biophys J 85: 2589–2598.J. ZamoonA. MascioniDD ThomasG. Veglia2003NMR solution structure and topological orientation of monomeric phospholamban in dodecylphosphocholine micelles.Biophys J8525892598

2. Post RL, Hegyvary C, Kume S (1972) Activation by adenosine triphosphate in the phosphorylation kinetics of sodium and potassium ion transport adenosine triphosphatase. J Biol Chem 247: 6530–6540.RL PostC. HegyvaryS. Kume1972Activation by adenosine triphosphate in the phosphorylation kinetics of sodium and potassium ion transport adenosine triphosphatase.J Biol Chem24765306540

3. Simmerman HK, Jones LR (1998) Phospholamban: protein structure, mechanism of action, and role in cardiac function. Physiol Rev 78: 921–947.HK SimmermanLR Jones1998Phospholamban: protein structure, mechanism of action, and role in cardiac function.Physiol Rev78921947

4. Traaseth NJ, Ha KN, Verardi R, Shi L, Buffy JJ, et al. (2008) Structural and dynamic basis of phospholamban and sarcolipin inhibition of Ca2+-ATPase. Biochemistry 47: 3–13.NJ TraasethKN HaR. VerardiL. ShiJJ Buffy2008Structural and dynamic basis of phospholamban and sarcolipin inhibition of Ca2+-ATPase.Biochemistry47313

5. Colyer J (1993) Control of the calcium pump of cardiac sarcoplasmic reticulum. A specific role for the pentameric structure of phospholamban? Cardiovasc Res 27: 1766–1771.J. Colyer1993Control of the calcium pump of cardiac sarcoplasmic reticulum. A specific role for the pentameric structure of phospholamban?Cardiovasc Res2717661771

6. Oxenoid K, Rice AJ, Chou JJ (2007) Comparing the structure and dynamics of phospholamban pentamer in its unphosphorylated and pseudo-phosphorylated states. Protein Sci 16: 1977–1983.K. OxenoidAJ RiceJJ Chou2007Comparing the structure and dynamics of phospholamban pentamer in its unphosphorylated and pseudo-phosphorylated states.Protein Sci1619771983

7. Haghighi K, Gregory KN, Kranias EG (2004) Sarcoplasmic reticulum Ca-ATPase-phospholamban interactions and dilated cardiomyopathy. Biochem Biophys Res Commun 322: 1214–1222.K. HaghighiKN GregoryEG Kranias2004Sarcoplasmic reticulum Ca-ATPase-phospholamban interactions and dilated cardiomyopathy.Biochem Biophys Res Commun32212141222

8. Schmitt JP, Kamisago M, Asahi M, Li GH, Ahmad F, et al. (2003) Dilated cardiomyopathy and heart failure caused by a mutation in phospholamban. Science 299: 1410–1413.JP SchmittM. KamisagoM. AsahiGH LiF. Ahmad2003Dilated cardiomyopathy and heart failure caused by a mutation in phospholamban.Science29914101413

9. Cornea RL, Autry JM, Chen Z, Jones LR (2000) Reexamination of the role of the leucine/isoleucine zipper residues of phospholamban in inhibition of the Ca2+ pump of cardiac sarcoplasmic reticulum. J Biol Chem 275: 41487–41494.RL CorneaJM AutryZ. ChenLR Jones2000Reexamination of the role of the leucine/isoleucine zipper residues of phospholamban in inhibition of the Ca2+ pump of cardiac sarcoplasmic reticulum.J Biol Chem2754148741494

10. Kimura Y, Kurzydlowski K, Tada M, MacLennan DH (1997) Phospholamban inhibitory function is activated by depolymerization. J Biol Chem 272: 15061–15064.Y. KimuraK. KurzydlowskiM. TadaDH MacLennan1997Phospholamban inhibitory function is activated by depolymerization.J Biol Chem2721506115064

11. Autry JM, Jones LR (1997) Functional Co-expression of the canine cardiac Ca2+ pump and phospholamban in Spodoptera frugiperda (Sf21) cells reveals new insights on ATPase regulation. J Biol Chem 272: 15872–15880.JM AutryLR Jones1997Functional Co-expression of the canine cardiac Ca2+ pump and phospholamban in Spodoptera frugiperda (Sf21) cells reveals new insights on ATPase regulation.J Biol Chem2721587215880

12. Thomas DD, Reddy LG, Karim CB, Li M, Cornea R, et al. (1998) Direct spectroscopic detection of molecular dynamics and interactions of the calcium pump and phospholamban. Ann N Y Acad Sci 853: 186–194.DD ThomasLG ReddyCB KarimM. LiR. Cornea1998Direct spectroscopic detection of molecular dynamics and interactions of the calcium pump and phospholamban.Ann N Y Acad Sci853186194

13. Li M, Reddy LG, Bennett R, Silva ND Jr, Jones LR, et al. (1999) A fluorescence energy transfer method for analyzing protein oligomeric structure: application to phospholamban. Biophys J 76: 2587–2599.M. LiLG ReddyR. BennettND Silva JrLR Jones1999A fluorescence energy transfer method for analyzing protein oligomeric structure: application to phospholamban.Biophys J7625872599

14. Metcalfe EE, Traaseth NJ, Veglia G (2005) Serine 16 phosphorylation induces an order-to-disorder transition in monomeric phospholamban. Biochemistry 44: 4386–4396.EE MetcalfeNJ TraasethG. Veglia2005Serine 16 phosphorylation induces an order-to-disorder transition in monomeric phospholamban.Biochemistry4443864396

15. Houndonougbo Y, Kuczera K, Jas GS (2005) Structure and dynamics of phospholamban in solution and in membrane bilayer: computer simulations. Biochemistry 44: 1780–1792.Y. HoundonougboK. KuczeraGS Jas2005Structure and dynamics of phospholamban in solution and in membrane bilayer: computer simulations.Biochemistry4417801792

16. Mascioni A, Karim C, Zamoon J, Thomas DD, Veglia G (2002) Solid-state NMR and rigid body molecular dynamics to determine domain orientations of monomeric phospholamban. J Am Chem Soc 124: 9392–9393.A. MascioniC. KarimJ. ZamoonDD ThomasG. Veglia2002Solid-state NMR and rigid body molecular dynamics to determine domain orientations of monomeric phospholamban.J Am Chem Soc12493929393

17. Traaseth NJ, Buffy JJ, Zamoon J, Veglia G (2006) Structural dynamics and topology of phospholamban in oriented lipid bilayers using multidimensional solid-state NMR. Biochemistry 45: 13827–13834.NJ TraasethJJ BuffyJ. ZamoonG. Veglia2006Structural dynamics and topology of phospholamban in oriented lipid bilayers using multidimensional solid-state NMR.Biochemistry451382713834

18. Paterlini MG, Thomas DD (2005) The alpha-helical propensity of the cytoplasmic domain of phospholamban: a molecular dynamics simulation of the effect of phosphorylation and mutation. Biophys J 88: 3243–3251.MG PaterliniDD Thomas2005The alpha-helical propensity of the cytoplasmic domain of phospholamban: a molecular dynamics simulation of the effect of phosphorylation and mutation.Biophys J8832433251

19. Pantano S, Carafoli E (2007) The role of phosphorylation on the structure and dynamics of phospholamban: a model from molecular simulations. Proteins 66: 930–940.S. PantanoE. Carafoli2007The role of phosphorylation on the structure and dynamics of phospholamban: a model from molecular simulations.Proteins66930940

20. Sugita Y, Miyashita N, Yoda T, Ikeguchi M, Toyoshima C (2006) Structural changes in the cytoplasmic domain of phospholamban by phosphorylation at Ser16: a molecular dynamics study. Biochemistry 45: 11752–11761.Y. SugitaN. MiyashitaT. YodaM. IkeguchiC. Toyoshima2006Structural changes in the cytoplasmic domain of phospholamban by phosphorylation at Ser16: a molecular dynamics study.Biochemistry451175211761

21. Metcalfe EE, Zamoon J, Thomas DD, Veglia G (2004) (1)H/(15)N heteronuclear NMR spectroscopy shows four dynamic domains for phospholamban reconstituted in dodecylphosphocholine micelles. Biophys J 87: 1205–1214.EE MetcalfeJ. ZamoonDD ThomasG. Veglia2004(1)H/(15)N heteronuclear NMR spectroscopy shows four dynamic domains for phospholamban reconstituted in dodecylphosphocholine micelles.Biophys J8712051214

22. Traaseth NJ, Shi L, Verardi R, Mullen DG, Barany G, et al. (2009) Structure and topology of monomeric phospholamban in lipid membranes determined by a hybrid solution and solid-state NMR approach. Proc Natl Acad Sci USA 106: 10165–10170.NJ TraasethL. ShiR. VerardiDG MullenG. Barany2009Structure and topology of monomeric phospholamban in lipid membranes determined by a hybrid solution and solid-state NMR approach.Proc Natl Acad Sci USA1061016510170

23. Kovacs RJ, Nelson MT, Simmerman HK, Jones LR (1988) Phospholamban forms Ca2+-selective channels in lipid bilayers. J Biol Chem 263: 18364–18368.RJ KovacsMT NelsonHK SimmermanLR Jones1988Phospholamban forms Ca2+-selective channels in lipid bilayers.J Biol Chem2631836418368

24. Decrouy A, Juteau M, Proteau S, Teijiera J, Rousseau E (1996) Biochemical regulation of sarcoplasmic reticulum Cl- channel from human atrial myocytes: involvement of phospholamban. J Mol Cell Cardiol 28: 767–780.A. DecrouyM. JuteauS. ProteauJ. TeijieraE. Rousseau1996Biochemical regulation of sarcoplasmic reticulum Cl- channel from human atrial myocytes: involvement of phospholamban.J Mol Cell Cardiol28767780

25. MacLennan DH, Kranias EG (2003) Phospholamban: a crucial regulator of cardiac contractility. Nat Rev Mol Cell Biol 4: 566–577.DH MacLennanEG Kranias2003Phospholamban: a crucial regulator of cardiac contractility.Nat Rev Mol Cell Biol4566577

26. Oxenoid K, Chou JJ (2005) The structure of phospholamban pentamer reveals a channel-like architecture in membranes. Proc Natl Acad Sci USA 102: 10870–10875.K. OxenoidJJ Chou2005The structure of phospholamban pentamer reveals a channel-like architecture in membranes.Proc Natl Acad Sci USA1021087010875

27. Potluri S, Yan AK, Chou JJ, Donald BR, Bailey-Kellogg C (2006) Structure determination of symmetric homo-oligomers by a complete search of symmetry configuration space, using NMR restraints and van der Waals packing. Proteins 65: 203–219.S. PotluriAK YanJJ ChouBR DonaldC. Bailey-Kellogg2006Structure determination of symmetric homo-oligomers by a complete search of symmetry configuration space, using NMR restraints and van der Waals packing.Proteins65203219

28. Robia SL, Flohr NC, Thomas DD (2005) Phospholamban pentamer quaternary conformation determined by in-gel fluorescence anisotropy. Biochemistry 44: 4302–4311.SL RobiaNC FlohrDD Thomas2005Phospholamban pentamer quaternary conformation determined by in-gel fluorescence anisotropy.Biochemistry4443024311

29. Sayadi M, Tanizaki S, Feig M Effect of membrane thickness on conformational sampling of phospholamban from computer simulations. Biophys J 98: 805–814.Sayadi M, Tanizaki S, Feig M Effect of membrane thickness on conformational sampling of phospholamban from computer simulations.Biophys J98805814

30. Kim T, Lee J, Im W (2009) Molecular dynamics studies on structure and dynamics of phospholamban monomer and pentamer in membranes. Proteins 76: 86–98.T. KimJ. LeeW. Im2009Molecular dynamics studies on structure and dynamics of phospholamban monomer and pentamer in membranes.Proteins768698

31. Maffeo C, Aksimentiev A (2009) Structure, dynamics, and ion conductance of the phospholamban pentamer. Biophys J 96: 4853–4865.C. MaffeoA. Aksimentiev2009Structure, dynamics, and ion conductance of the phospholamban pentamer.Biophys J9648534865

32. Chou KC (1988) Review: Low-frequency collective motion in biomacromolecules and its biological functions. Biophys Chem 30: 3–48.KC Chou1988Review: Low-frequency collective motion in biomacromolecules and its biological functions.Biophys Chem30348

33. Chou KC (1989) Low-frequency resonance and cooperativity of hemoglobin. Trends Biochem Sci 14: 212.KC Chou1989Low-frequency resonance and cooperativity of hemoglobin.Trends Biochem Sci14212

34. Wang JF, Wei DQ, Li L, Zheng SY, Li YX, et al. (2007) 3D structure modeling of cytochrome P450 2C19 and its implication for personalized drug design. Biochem Biophys Res Commun 355: 513-519. Erratum in: Biochem Biophys Res Commun 2007, 357: 330. Biochem Biophys Res Commun 2009, 384: 399.JF WangDQ WeiL. LiSY ZhengYX Li20073D structure modeling of cytochrome P450 2C19 and its implication for personalized drug design. Biochem Biophys Res Commun 355: 513-519. Erratum in: Biochem Biophys Res Commun 2007, 357: 330.Biochem Biophys Res Commun 2009,384399

35. Wang JF, Wei DQ, Lin Y, Wang YH, Du HL, et al. (2007) Insights from modeling the 3D structure of NAD(P)H-dependent D-xylose reductase of Pichia stipitis and its binding interactions with NAD and NADP. Biochem Biophys Res Commun 359: 323–329.JF WangDQ WeiY. LinYH WangHL Du2007Insights from modeling the 3D structure of NAD(P)H-dependent D-xylose reductase of Pichia stipitis and its binding interactions with NAD and NADP.Biochem Biophys Res Commun359323329

36. Wang JF, Wei DQ, Chen C, Li Y, Chou KC (2008) Molecular modeling of two CYP2C19 SNPs and its implications for personalized drug design. Protein Pept Lett 15: 27–32.JF WangDQ WeiC. ChenY. LiKC Chou2008Molecular modeling of two CYP2C19 SNPs and its implications for personalized drug design.Protein Pept Lett152732

37. Wang JF, Chou KC (2010) Insights from studying the mutation-induced allostery in the M2 proton channel by molecular dynamics. Protein Eng Des Sel 23: 663–666.JF WangKC Chou2010Insights from studying the mutation-induced allostery in the M2 proton channel by molecular dynamics.Protein Eng Des Sel23663666

38. Wang JF, Chou KC (2009) Insight into the molecular switch mechanism of human Rab5a from molecular dynamics simulations. Biochem Biophys Res Commun 390: 608–612.JF WangKC Chou2009Insight into the molecular switch mechanism of human Rab5a from molecular dynamics simulations.Biochem Biophys Res Commun390608612

39. Wang JF, Gong K, Wei DQ, Li YX (2009) Structure flexibility and interactions of PTP1B's S-loop. Interdiscip Sci 1: 214–219.JF WangK. GongDQ WeiYX Li2009Structure flexibility and interactions of PTP1B's S-loop.Interdiscip Sci1214219

40. Wang JF, Wei DQ, Chou KC (2009) Insights from investigating the interactions of adamantane-based drugs with the M2 proton channel from the H1N1 swine virus. Biochem Biophys Res Commun 388: 413–417.JF WangDQ WeiKC Chou2009Insights from investigating the interactions of adamantane-based drugs with the M2 proton channel from the H1N1 swine virus.Biochem Biophys Res Commun388413417

41. Wang JF, Yan JY, Wei DQ, Chou KC (2009) Binding of CYP2C9 with diverse drugs and its implications for metabolic mechanism. Med Chem 5: 263–270.JF WangJY YanDQ WeiKC Chou2009Binding of CYP2C9 with diverse drugs and its implications for metabolic mechanism.Med Chem5263270

42. Wang JF, Wei DQ (2009) Role of structural bioinformatics and traditional Chinese medicine databases in pharmacogenomics. Pharmacogenomics 10: 1213–1215.JF WangDQ Wei2009Role of structural bioinformatics and traditional Chinese medicine databases in pharmacogenomics.Pharmacogenomics1012131215

43. Wang Y, Wei DQ, Wang JF (2010) Molecular dynamics studies on T1 lipase: insight into a double-flap mechanism. J Chem Inf Model 50: 875–878.Y. WangDQ WeiJF Wang2010Molecular dynamics studies on T1 lipase: insight into a double-flap mechanism.J Chem Inf Model50875878

44. Wang JF, Wei DQ, Chou KC (2008) Drug candidates from traditional chinese medicines. Curr Top Med Chem 8: 1656–1665.JF WangDQ WeiKC Chou2008Drug candidates from traditional chinese medicines.Curr Top Med Chem816561665

45. Gong K, Li L, Wang JF, Cheng F, Wei DQ, et al. (2009) Binding mechanism of H5N1 influenza virus neuraminidase with ligands and its implication for drug design. Med Chem 5: 242–249.K. GongL. LiJF WangF. ChengDQ Wei2009Binding mechanism of H5N1 influenza virus neuraminidase with ligands and its implication for drug design.Med Chem5242249

46. Gu RX, Gu H, Xie ZY, Wang JF, Arias HR, et al. (2009) Possible drug candidates for Alzheimer's disease deduced from studying their binding interactions with alpha7 nicotinic acetylcholine receptor. Med Chem 5: 250–262.RX GuH. GuZY XieJF WangHR Arias2009Possible drug candidates for Alzheimer's disease deduced from studying their binding interactions with alpha7 nicotinic acetylcholine receptor.Med Chem5250262

47. Guo X, Wang JF, Zhu Y, Wei DQ (2010) Recent Progress on Computer-Aided Inhibitor Design of H5N1 Influenza A Virus. Curr Comput Aided Drug Des 6: 139–146.X. GuoJF WangY. ZhuDQ Wei2010Recent Progress on Computer-Aided Inhibitor Design of H5N1 Influenza A Virus.Curr Comput Aided Drug Des6139146

48. Berman HM, Battistuz T, Bhat TN, Bluhm WF, Bourne PE, et al. (2002) The Protein Data Bank. Acta Crystallogr D Biol Crystallogr 58: 899–907.HM BermanT. BattistuzTN BhatWF BluhmPE Bourne2002The Protein Data Bank.Acta Crystallogr D Biol Crystallogr58899907

49. Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klei ML (1983) Comparison of simple potential functions for simulating liquid water. J Chem Phys 79: 926–935.WL JorgensenJ. ChandrasekharJD MaduraRW ImpeyML Klei1983Comparison of simple potential functions for simulating liquid water.J Chem Phys79926935

50. Phillips JC, Braun R, Wang W, Gumbart J, Tajkhorshid E, et al. (2005) Scalable molecular dynamics with NAMD. J Comput Chem 26: 1781–1802.JC PhillipsR. BraunW. WangJ. GumbartE. Tajkhorshid2005Scalable molecular dynamics with NAMD.J Comput Chem2617811802

51. Feller SE, MacKerell AD (2000) Am improved empirical potential energy function for molecular simulations of phospholipids. J Phys Chem B 104: 7510–7515.SE FellerAD MacKerell2000Am improved empirical potential energy function for molecular simulations of phospholipids.J Phys Chem B10475107515

52. Nose S (1984) A unified formulation of the constant temperature molecular dynamics methods. J Chem Phys 81: 511–519.S. Nose1984A unified formulation of the constant temperature molecular dynamics methods.J Chem Phys81511519

53. Hoover WG (1985) Canonical dynamics: Equilibrium phase-space destributions. Phys Rev A 31: 1695–1697.WG Hoover1985Canonical dynamics: Equilibrium phase-space destributions.Phys Rev A3116951697

54. Darden T, York D, Pedersen L (1993) Particle Mesh Ewald: an N-log(N) method for Ewald sums in large systems. J Chem Phys 98: 10089–10092.T. DardenD. YorkL. Pedersen1993Particle Mesh Ewald: an N-log(N) method for Ewald sums in large systems.J Chem Phys981008910092

55. Chen AA, Pappu RV (2007) Quantitative characterization of ion pairing and cluster formation in strong 1∶1 electrolytes. J Phys Chem B 111: 6469–6478.AA ChenRV Pappu2007Quantitative characterization of ion pairing and cluster formation in strong 1∶1 electrolytes.J Phys Chem B11164696478

56. Lee SJ, Song Y, Baker NA (2008) Molecular dynamics simulaitons of asymmetric NaCl and KCl solutions separated by phosphatidylcholne bilayers: potential drops and structural changes induced by strong Na+-lipid interactions and finite size effect. Biophys J 94: 3565–3576.SJ LeeY. SongNA Baker2008Molecular dynamics simulaitons of asymmetric NaCl and KCl solutions separated by phosphatidylcholne bilayers: potential drops and structural changes induced by strong Na+-lipid interactions and finite size effect.Biophys J9435653576

57. Zamoon J, Nitu F, Karim C, Thomas DD, Veglia G (2005) Mapping the interaction surface of a membrane protein: unveiling the conformational switch of phospholamban in calcium pump regulation. Proc Natl Acad Sci USA 102: 4747–4752.J. ZamoonF. NituC. KarimDD ThomasG. Veglia2005Mapping the interaction surface of a membrane protein: unveiling the conformational switch of phospholamban in calcium pump regulation.Proc Natl Acad Sci USA10247474752

58. Karim CB, Kirby TL, Zhang Z, Nesmelov Y, Thomas DD (2004) Phospholamban structural dynamics in lipid bilayers probed by a spin label rigidly coupled to the peptide backbone. Proc Natl Acad Sci USA 101: 14437–14442.CB KarimTL KirbyZ. ZhangY. NesmelovDD Thomas2004Phospholamban structural dynamics in lipid bilayers probed by a spin label rigidly coupled to the peptide backbone.Proc Natl Acad Sci USA1011443714442

59. Karim CB, Zhang Z, Howard EC, Torgersen KD, Thomas DD (2006) Phosphorylation-dependent conformational switch in spin-labeled phospholamban bound to SERCA. J Mol Biol 358: 1032–1040.CB KarimZ. ZhangEC HowardKD TorgersenDD Thomas2006Phosphorylation-dependent conformational switch in spin-labeled phospholamban bound to SERCA.J Mol Biol35810321040

60. Traaseth NJ, Thomas , DD , Veglia G (2006) Effects of Ser16 phosphorylation on the allosteric transitions of phospholamban/Ca2+-ATPase complex. J Mol Biol 358: 1041–1050.NJ TraasethThomasDDG. Veglia2006Effects of Ser16 phosphorylation on the allosteric transitions of phospholamban/Ca2+-ATPase complex.J Mol Biol35810411050