The role of phosphorylation on the structure and dynamics of phospholamban: a model from molecular simulations

Fat SIRAH: Coarse-Grained Phospholipids To Explore Membrane-Protein Dynamics

The capability to handle highly heterogeneous molecular assemblies in a consistent manner is among the greatest challenges faced when deriving simulation parameters. This is particularly the case for coarse-grained (CG) simulations in which chemical functional groups are lumped into effective interaction centers for which transferability between different chemical environments is not guaranteed. Here, we introduce the parametrization of a set of CG phospholipids compatible with the latest version of the SIRAH force field for proteins. The newly introduced lipid species include different acylic chain lengths and partial unsaturation, as well as polar and acidic head groups that show a very good reproduction of structural membrane determinants, such as areas per lipid, thickness, order parameter, etc., and their dependence with temperature. Simulation of membrane proteins showed unprecedented accuracy in the unbiased description of the thickness-dependent membrane-protein orientation in systems where this information is experimentally available (namely, the SarcoEndoplasmic Reticulum Calcium-SERCA-pump and its regulator Phospholamban). The interactions that lead to this faithful reproduction can be traced down to the single amino acid-lipid interaction level and show full agreement with biochemical data present in the literature. Finally, the present parametrization is implemented in the GROMACS and AMBER simulation packages facilitating its use by a wide portion of the biocomputing community.

Protein Interaction with Charged Macromolecules: From Model Polymers to Unfolded Proteins and Post-Translational Modifications

Interaction of proteins with charged macromolecules is involved in many processes in cells. Firstly, there are many naturally occurred charged polymers such as DNA and RNA, polyphosphates, sulfated glycosaminoglycans, etc., as well as pronouncedly charged proteins such as histones or actin. Electrostatic interactions are also important for "generic" proteins, which are not generally considered as polyanions or polycations. Finally, protein behavior can be altered due to post-translational modifications such as phosphorylation, sulfation, and glycation, which change a local charge of the protein region. Herein we review molecular modeling for the investigation of such interactions, from model polyanions and polycations to unfolded proteins. We will show that electrostatic interactions are ubiquitous, and molecular dynamics simulations provide an outstanding opportunity to look inside binding and reveal the contribution of electrostatic interactions. Since a molecular dynamics simulation is only a model, we will comprehensively consider its relationship with the experimental data.

Perturbations of Native Membrane Protein Structure in Alkyl Phosphocholine Detergents: A Critical Assessment of NMR and Biophysical Studies

Membrane proteins perform a host of vital cellular functions. Deciphering the molecular mechanisms whereby they fulfill these functions requires detailed biophysical and structural investigations. Detergents have proven pivotal to extract the protein from its native surroundings. Yet, they provide a milieu that departs significantly from that of the biological membrane, to the extent that the structure, the dynamics, and the interactions of membrane proteins in detergents may considerably vary, as compared to the native environment. Understanding the impact of detergents on membrane proteins is, therefore, crucial to assess the biological relevance of results obtained in detergents. Here, we review the strengths and weaknesses of alkyl phosphocholines (or foscholines), the most widely used detergent in solution-NMR studies of membrane proteins. While this class of detergents is often successful for membrane protein solubilization, a growing list of examples points to destabilizing and denaturing properties, in particular for α-helical membrane proteins. Our comprehensive analysis stresses the importance of stringent controls when working with this class of detergents and when analyzing the structure and dynamics of membrane proteins in alkyl phosphocholine detergents.

Self-assembling study of sarcolipin and its mutants in multiple molecular dynamic simulations

The Sarcolipin (SLN) is a single trans-membrane protein that can self-assembly to dimer and oligomer for playing importantphysiological function. In this work, we addressed the dimerization of wild type SLN (wSLN) and its mutants (mSLNs) - I17A and I20A, using both coarse-grained (CG) and atomistic (AT) molecular dynamics (MD) simulations. Our results demonstrated that wSLN homodimer assembled as a left-handed helical complex, while mSLNs heterodimers assembled as right-handed complexes. Analysis of residue-residue contacts map indicated that isoleucine (Ile)-leucione (Leu) zipper domain played an important role in dimerization. The potential of mean force (PMF) demonstrated that wSLN homodimer was more stable than mSLNs heterodimers. Meanwhile, the mSLNs heterodimers preferred right-handed rather than left-handed helix. AT-MD simulations for wSLN and mSLNs were also in line with CG-MD simulations. These results provided the insights for understanding the mechanisms of SLNs self-assembling. Proteins 2017; 85:1065-1077. © 2017 Wiley Periodicals, Inc.

Transmembrane dynamics of the Thr-5 phosphorylated sarcolipin pentameric channel

Sarcolipin (SLN), an important membrane protein expressed in the sarcoplasmic reticulum (SR), regulates muscle contractions in cardiac and skeletal muscle. The phosphorylation at amino acid Thr5 of the SLN protein modulates the amount of Ca(2+) that passes through the SR. Using molecular dynamics simulation, we evaluated the phosphorylation at Thr5 of pentameric SLN (phospho-SLN) channel's energy barrier and pore characteristics by calculating the potential of mean force (PMF) along the channel pore and determining the diffusion coefficient. The results indicate that pentameric phospho-SLN promotes penetration of monovalent and divalent ions through the channel. The analysis of PMF, pore radius and diffusion coefficient indicates that Leu21 is the hydrophobic gate of the pentameric SLN channel. In the channel, water molecules near the Leu21 pore demonstrated a clear hydrated-dehydrated transition; however, the mutation of Leu21 to an Alanine (L21A) destroyed the hydrated-dehydrated transitions. These water-dynamic behaviors and PMF confirm that Leu21 is the key residue that regulates the ion permeability of the pentameric SLN channel. These results provide the structural-basis insights and molecular-dynamic information that are needed to understand the regulatory mechanisms of ion permeability in the pentameric SLN channel.

Molecular Dynamics Simulations of the Human Glucose Transporter GLUT1

Glucose transporters (GLUTs) provide a pathway for glucose transport across membranes. Human GLUTs are implicated in devastating diseases such as heart disease, hyper- and hypo-glycemia, type 2 diabetes and cancer. The human GLUT1 has been recently crystalized in the inward-facing open conformation. However, there is no other structural information for other conformations. The X-ray structures of E. coli Xylose permease (XylE), a glucose transporter homolog, are available in multiple conformations with and without the substrates D-xylose and D-glucose. XylE has high sequence homology to human GLUT1 and key residues in the sugar-binding pocket are conserved. Here we construct a homology model for human GLUT1 based on the available XylE crystal structure in the partially occluded outward-facing conformation. A long unbiased all atom molecular dynamics simulation starting from the model can capture a new fully opened outward-facing conformation. Our investigation of molecular interactions at the interface between the transmembrane (TM) domains and the intracellular helices (ICH) domain in the outward- and inward-facing conformation supports that the ICH domain likely stabilizes the outward-facing conformation in GLUT1. Furthermore, inducing a conformational transition, our simulations manifest a global asymmetric rocker switch motion and detailed molecular interactions between the substrate and residues through the water-filled selective pore along a pathway from the extracellular to the intracellular side. The results presented here are consistent with previously published biochemical, mutagenesis and functional studies. Together, this study shed light on the structure and functional relationships of GLUT1 in multiple conformational states.

The p.Cys169Tyr variant of connexin 26 is not a polymorphism

Mutations in the GJB2 gene, which encodes the gap junction protein connexin 26 (Cx26), are the primary cause of hereditary prelingual hearing impairment. Here, the p.Cys169Tyr missense mutation of Cx26 (Cx26C169Y), previously classified as a polymorphism, has been identified as causative of severe hearing loss in two Qatari families. We have analyzed the effect of this mutation using a combination of confocal immunofluorescence microscopy and molecular dynamics simulations. At the cellular level, our results show that the mutant protein fails to form junctional channels in HeLa transfectants despite being correctly targeted to the plasma membrane. At the molecular level, this effect can be accounted for by disruption of the disulfide bridge that Cys169 forms with Cys64 in the wild-type structure (Cx26WT). The lack of the disulfide bridge in the Cx26C169Y protein causes a spatial rearrangement of two important residues, Asn176 and Thr177. In the Cx26WT protein, these residues play a crucial role in the intra-molecular interactions that permit the formation of an intercellular channel by the head-to-head docking of two opposing hemichannels resident in the plasma membrane of adjacent cells. Our results elucidate the molecular pathogenesis of hereditary hearing loss due to the connexin mutation and facilitate the understanding of its role in both healthy and affected individuals.

Proton-coupled sugar transport in the prototypical major facilitator superfamily protein XylE

The major facilitator superfamily (MFS) is the largest collection of structurally related membrane proteins that transport a wide array of substrates. The proton-coupled sugar transporter XylE is the first member of the MFS that has been structurally characterized in multiple transporting conformations, including both the outward and inward-facing states. Here we report the crystal structure of XylE in a new inward-facing open conformation, allowing us to visualize the rocker-switch movement of the N-domain against the C-domain during the transport cycle. Using molecular dynamics simulation, and functional transport assays, we describe the movement of XylE that facilitates sugar translocation across a lipid membrane and identify the likely candidate proton-coupling residues as the conserved Asp27 and Arg133. This study addresses the structural basis for proton-coupled substrate transport and release mechanism for the sugar porter family of proteins.

Glucose transporters are a medically important class of membrane proteins often deregulated in diseases such as Type 2 diabetes. Here, Wisedchaisri et al. report the crystal structure of XylE in an inward-facing open conformation to provide a general mechanism of substrate transport for the sugar porter family of proteins.

Molecular dynamics simulations highlight structural and functional alterations in deafness-related M34T mutation of connexin 26

Mutations of the GJB2 gene encoding the connexin 26 (Cx26) gap junction protein, which is widely expressed in the inner ear, are the primary cause of hereditary non-syndromic hearing loss in several populations. The deafness–associated single amino acid substitution of methionine 34 (M34) in the first transmembrane helix (TM1) with a threonine (T) ensues in the production of mutant Cx26M34T channels that are correctly synthesized and assembled in the plasma membrane. However, mutant channels overexpressed in HeLa cells retain only 11% of the wild type unitary conductance. Here we extend and rationalize those findings by comparing wild type Cx26 (Cx26WT) and Cx26M34T mutant channels in silico, using molecular dynamics simulations. Our results indicate that the quaternary structure of the Cx26M34T hemichannel is altered at the level of the pore funnel due to the disruption of the hydrophobic interaction between M34 and tryptophan 3 (W3) in the N–terminal helix (NTH). Our simulations also show that external force stimuli applied to the NTHs can detach them from the inner wall of the pore more readily in the mutant than in the wild type hemichannel. These structural alterations significantly increase the free energy barrier encountered by permeating ions, correspondingly decreasing the unitary conductance of the Cx26M34T hemichannel. Our results accord with the proposal that the mutant resides most of the time in a low conductance state. However, the small displacement of the NTHs in our Cx26M34T hemichannel model is not compatible with the formation of a pore plug as in the related Cx26M34A mutant.

Calcium pumps: why so many?

Ca(2+)-ATPases (pumps) are key to the regulation of Ca(2+) in eukaryotic cells: nine are known today, belonging to three multigene families. The three endo(sarco)plasmic reticulum (SERCA) and the four plasma membrane (PMCA) pumps have been known for decades, the two Secretory Pathway Ca(2+) ATPase (SPCA) pumps have only become known recently. The number of pump isoforms is further increased by alternative splicing processes. The three pump types share the basic features of the catalytic mechanism, but differ in a number of properties related to tissue distribution, regulation, and role in the cellular homeostasis of Ca(2+). The molecular understanding of the function of all pumps has received great impetus from the solution of the three-dimensional (3D) structure of one of them, the SERCA pump. This landmark structural advance has been accompanied by the emergence and rapid expansion of the area of pump malfunction. Most of the pump defects described so far are genetic and produce subtler, often tissue and isoform specific, disturbances that affect individual components of the Ca(2+)-controlling and/or processing machinery, compellingly indicating a specialized role for each Ca(2+) pump type and/or isoform.

Revised AMBER parameters for bioorganic phosphates

We report AMBER force field parameters for biological simulations involving phosphorylation of serine, threonine or tyrosine. The initial parameters used RESP fitting for the atomic partial charges and standard values for all other parameters such as Lennard-Jones coefficients. These were refined with the aid of a thermodynamic cycle consisting of experimentally determined pKa values, solvation energies from molecular dynamics free energy simulations, and gas phase basicities from QM calculations. A polarization energy term was included to account for the charge density change between the gas-phase and solution, and solvation free energies were determined using thermodynamic integration. Parameter adjustment is required to obtain consistent thermodynamic results with better balanced electrostatic interactions between water and the phosphate oxygens. To achieve this we modified the phosphate oxygen radii. A thermodynamically consistent parameter set can be derived for monoanions and requires an increase of the van der Waals phosphate oxygen radii of approximately 0.09 Å. Larger, residue-specific radii appear to be needed for dianions. The revised parameters developed here should be of particular interest for environments where simulations of multiple protonation states may be of interest.

Role of conformational sampling of Ser16 and Thr17-phosphorylated phospholamban in interactions with SERCA

Phosphorylation of phospholamban (PLB) at Ser16 and/ or Thr17 is believed to release its inhibitory effect on sarcoplasmic reticulum calcium ATPase. Ser16 phosphorylation of PLB has been suggested to cause a conformational change that alters the interaction between the enzyme and protein. Using computer simulations, the conformational sampling of Ser16 phosphorylated PLB in implicit membrane environment is compared here with the unphosphorylated PLB system to investigate these conformational changes. The results suggest that conformational changes in the cytoplasmic domain of PLB upon phosphorylation at Ser16 increase the likelihood of unfavorable interactions with SERCA in the E2 state prompting a conformational switch of SERCA from E2 to E1. Phosphorylation of PLB at Thr17 on the other hand does not appear to affect interactions with SERCA significantly suggesting that the mechanism of releasing the inhibitory effect is different between Thr17 phosphorylated and Ser16 phosphorylated PLB.

An allosteric mechanism inferred from molecular dynamics simulations on phospholamban pentamer in lipid membranes

Phospholamban functions as a regulator of Ca(2+) concentration of cardiac muscle cells by triggering the bioactivity of sarcoplasmic reticulum Ca(2+)-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 Ca(2+)-ATPase.

Effect of membrane thickness on conformational sampling of phospholamban from computer simulations

The conformational sampling of monomeric, membrane-bound phospholamban is described from computer simulations. Phospholamban (PLB) plays a key role as a regulator of sarcoplasmic reticulum calcium ATPase. An implicit membrane model is used in conjunction with replica exchange molecular dynamics simulations to reach mus-ms timescales. The implicit membrane model was also used to study the effect of different membrane thicknesses by scaling the low-dielectric region. The conformational sampling with the membrane model mimicking dipalmitoylphosphatidylcholine bilayers is in good agreement overall with experimental measurements, but consists of a wide variety of different conformations including structures not described previously. The conformational ensemble shifts significantly in the presence of thinner or thicker membranes. This has implications for the structure and dynamics of PLB in physiological membranes and offers what we believe to be a new interpretation of previous experimental measurements of PLB in detergents and microsomal membrane.

Phospholamban modulates the functional coupling between nucleotide domains in Ca-ATPase oligomeric complexes in cardiac sarcoplasmic reticulum

Oligomeric interactions between Ca-ATPase polypeptide chains and their modulation by phospholamban (PLB) were measured in native cardiac sarcoplasmic reticulum (SR) microsomes. Progressive modification of Lys(514) with fluorescein 5-isothiocyanate (FITC), which physically blocks access to the nucleotide binding site by ATP, demonstrates that Ca-ATPase active sites function independently of one another prior to the phosphorylation of PLB. However, upon cAMP-dependent protein kinase (PKA) phosphorylation of PLB, a second-order dependence between residual enzyme activity and the fraction of active sites is observed, consistent with a dimeric functional complex. Complementary distance measurements were made using FITC or 5-iodoacetamidofluorescein (IAF) bound to Cys(674) within the N- or P-domains, respectively, to detect structural coupling within oligomeric complexes. Accompanying the phosphorylation of PLB, neighboring Ca-ATPase polypeptide chains exhibit a 4 +/- 2 A decrease in the proximity between FITC sites within the N-domain and a 9 +/- 3 A increase in the proximity between IAF sites within P-domains. Thus, the phosphorylation of PLB induces spatial rearrangements between the N- and P-domain elements of proximal Ca-ATPase polypeptide chains which restore functional interactions between neighboring polypeptide chains and, in turn, result in increased rates of catalytic turnover. These results are interpreted in terms of a structural model, calculated through optimization of shape complementarity, desolvation, and electrostatic energies, which suggests a dimeric arrangement of Ca-ATPase polypeptide chains through the proximal association of N-domains that accommodates interaction with PLB. We suggest that the phosphorylation of PLB acts to release constraints involving interdomain subunit interactions that enhance catalytically important N-domain motions.

Self-guided Langevin dynamics study of regulatory interactions in NtrC

Multiple self-guided Langevin dynamics (SGLD) simulations were performed to examine structural and dynamical properties of the receiver domain of nitrogen regulatory protein C (NtrC(r)). SGLD and MD simulations of the phosphorylated active form structure suggest a mostly stable but broad structural ensemble of this protein. The finite difference Poisson-Boltzmann calculations of the pK(a) values of the active site residues suggest an increase in the pK(a) of His-84 on phosphorylation of Asp-54. In SGLD simulations of the phosphorylated active form with charged His-84, the average position of the regulatory helix alpha4 is found closer to the starting structure than in simulations with the neutral His-84. To model the transition pathway, the phosphate group was removed from the simulations. After 7 ns of simulations, the regulatory helix alpha4 was found approximately halfway between positions in the NMR structures of the active and inactive forms. Removal of the phosphate group stimulated loss of helix alpha4, suggesting that the pathway of conformational transition may involve partial unfolding mechanism. The study illustrates the potential utility of the SGLD method in studies of the coupling between ligand binding and conformational transitions.

Structure, dynamics, and ion conductance of the phospholamban pentamer

A 52-residue membrane protein, phospholamban (PLN) is an inhibitor of an adenosine-5'-triphosphate-driven calcium pump, the Ca2+-ATPase. Although the inhibition of Ca2+-ATPase involves PLN monomers, in a lipid bilayer membrane, PLN monomers form stable pentamers of unknown biological function. The recent NMR structure of a PLN pentamer depicts cytoplasmic helices extending normal to the bilayer in what is known as the bellflower conformation. The structure shows transmembrane helices forming a hydrophobic pore 4 A in diameter, which is reminiscent of earlier reports of possible ion conductance through PLN pentamers. However, recent FRET measurements suggested an alternative structure for the PLN pentamer, known as the pinwheel model, which features a narrower transmembrane pore and cytoplasmic helices that lie against the bilayer. Here, we report on structural dynamics and conductance properties of the PLN pentamers from all-atom (AA) and coarse-grained (CG) molecular dynamics simulations. Our AA simulations of the bellflower model demonstrate that in a lipid bilayer membrane or a detergent micelle, the cytoplasmic helices undergo large structural fluctuations, whereas the transmembrane pore shrinks and becomes asymmetric. Similar asymmetry of the transmembrane region was observed in the AA simulations of the pinwheel model; the cytoplasmic helices remained in contact with the bilayer. Using the CG approach, structural dynamics of both models were investigated on a microsecond timescale. The cytoplasmic helices of the CG bellflower model were observed to fall against the bilayer, whereas in the CG pinwheel model the conformation of the cytoplasmic helices remained stable. Using steered molecular dynamics simulations, we investigated the feasibility of ion conductance through the pore of the bellflower model. The resulting approximate potentials of mean force indicate that the PLN pentamer is unlikely to function as an ion channel.

A fully atomistic model of the Cx32 connexon

Connexins are plasma membrane proteins that associate in hexameric complexes to form channels named connexons. Two connexons in neighboring cells may dock to form a "gap junction" channel, i.e. an intercellular conduit that permits the direct exchange of solutes between the cytoplasm of adjacent cells and thus mediate cell-cell ion and metabolic signaling. The lack of high resolution data for connexon structures has hampered so far the study of the structure-function relationships that link molecular effects of disease-causing mutations with their observed phenotypes. Here we present a combination of modeling techniques and molecular dynamics (MD) to infer side chain positions starting from low resolution structures containing only C alpha atoms. We validated this procedure on the structure of the KcsA potassium channel, which is solved at atomic resolution. We then produced a fully atomistic model of a homotypic Cx32 connexon starting from a published model of the C alpha carbons arrangement for the connexin transmembrane helices, to which we added extracellular and cytoplasmic loops. To achieve structural relaxation within a realistic environment, we used MD simulations inserted in an explicit solvent-membrane context and we subsequently checked predictions of putative side chain positions and interactions in the Cx32 connexon against a vast body of experimental reports. Our results provide new mechanistic insights into the effects of numerous spontaneous mutations and their implication in connexin-related pathologies. This model constitutes a step forward towards a structurally detailed description of the gap junction architecture and provides a structural platform to plan new biochemical and biophysical experiments aimed at elucidating the structure of connexin channels and hemichannels.

Tandem phosphorylation of Ser-911 and Thr-912 at the C terminus of yeast plasma membrane H+-ATPase leads to glucose-dependent activation

In recent years there has been growing interest in the post-translational regulation of P-type ATPases by protein kinase-mediated phosphorylation. Pma1 H(+)-ATPase, which is responsible for H(+)-dependent nutrient uptake in yeast (Saccharomyces cerevisiae), is one such example, displaying a rapid 5-10-fold increase in activity when carbon-starved cells are exposed to glucose. Activation has been linked to Ser/Thr phosphorylation in the C-terminal tail of the ATPase, but the specific phosphorylation sites have not previously been mapped. The present study has used nanoflow high pressure liquid chromatography coupled with electrospray electron transfer dissociation tandem mass spectrometry to identify Ser-911 and Thr-912 as two major phosphorylation sites that are clearly related to glucose activation. In carbon-starved cells with low Pma1 activity, peptide 896-918, which was derived from the C terminus upon Lys-C proteolysis, was found to be singly phosphorylated at Thr-912, whereas in glucose-metabolizing cells with high ATPase activity, the same peptide was doubly phosphorylated at Ser-911 and Thr-912. Reciprocal (14)N/(15)N metabolic labeling of cells was used to measure the relative phosphorylation levels at the two sites. The addition of glucose to carbon-starved cells led to a 3-fold reduction in the singly phosphorylated form and an 11-fold increase in the doubly phosphorylated form. These results point to a mechanism in which the stepwise phosphorylation of two tandemly positioned residues near the C terminus mediates glucose-dependent activation of the H(+)-ATPase.

Mechanism of reversal of phospholamban inhibition of the cardiac Ca2+-ATPase by protein kinase A and by anti-phospholamban monoclonal antibody 2D12

Our model of phospholamban (PLB) regulation of the cardiac Ca(2+)-ATPase in sarcoplasmic reticulum (SERCA2a) states that PLB binds to the Ca(2+)-free, E2 conformation of SERCA2a and blocks it from transitioning from E2 to E1, the Ca(2+)-bound state. PLB and Ca(2+) binding to SERCA2a are mutually exclusive, and PLB inhibition of SERCA2a is manifested as a decreased apparent affinity of SERCA2a for Ca(2+). Here we extend this model to explain the reversal of SERCA2a inhibition that occurs after phosphorylation of PLB at Ser(16) by protein kinase A (PKA) and after binding of the anti-PLB monoclonal antibody 2D12, which recognizes residues 7-13 of PLB. Site-specific cysteine variants of PLB were co-expressed with SERCA2a, and the effects of PKA phosphorylation and 2D12 on Ca(2+)-ATPase activity and cross-linking to SERCA2a were monitored. In Ca(2+)-ATPase assays, PKA phosphorylation and 2D12 partially and completely reversed SERCA2a inhibition by decreasing K(Ca) values for enzyme activation, respectively. In cross-linking assays, cross-linking of PKA-phosphorylated PLB to SERCA2a was inhibited at only two of eight sites when conducted in the absence of Ca(2+) favoring E2. However, at a subsaturating Ca(2+) concentration supporting some E1, cross-linking of phosphorylated PLB to SERCA2a was attenuated at all eight sites. K(Ca) values for cross-linking inhibition were decreased nearly 2-fold at all sites by PLB phosphorylation, demonstrating that phosphorylated PLB binds more weakly to SERCA2a than dephosphorylated PLB. In parallel assays, 2D12 blocked PLB cross-linking to SERCA2a at all eight sites regardless of Ca(2+) concentration. Our results demonstrate that 2D12 restores maximal Ca(2+)-ATPase activity by physically disrupting the binding interaction between PLB and SERCA2a. Phosphorylation of PLB by PKA weakens the binding interaction between PLB and SERCA2a (yielding more PLB-free SERCA2a molecules at intermediate Ca(2+) concentrations), only partially restoring Ca(2+) affinity and Ca(2+)-ATPase activity.