Structures of the excited states of phospholamban and shifts in their populations upon phosphorylation

Pathological mutations in the phospholamban cytoplasmic region affect its topology and dynamics modulating the extent of SERCA inhibition

Phospholamban (PLN) is a 52 amino acid regulin that allosterically modulates the activity of the sarco(endo)plasmic reticulum Ca2+-ATPase (SERCA) in the heart muscle. In its unphosphorylated form, PLN binds SERCA within its transmembrane (TM) domains, approximately 20 Å away from the Ca2+ binding site, reducing SERCA's apparent Ca2+ affinity (pKCa) and decreasing cardiac contractility. During the enzymatic cycle, the inhibitory TM domain of PLN remains anchored to SERCA, whereas its cytoplasmic region transiently binds the ATPase's headpiece. Phosphorylation of PLN at Ser16 by protein kinase A increases the affinity of its cytoplasmic domain to SERCA, weakening the TM interactions with the ATPase, reversing its inhibitory function, and augmenting muscle contractility. How the structural changes caused by pathological mutations in the PLN cytoplasmic region are transmitted to its inhibitory TM domain is still unclear. Using solid-state NMR spectroscopy and activity assays, we analyzed the structural and functional effects of a series of mutations and their phosphorylated forms located in the PLN cytoplasmic region and linked to dilated cardiomyopathy. We found that these missense mutations affect the overall topology and dynamics of PLN and ultimately modulate its inhibitory potency. Also, the changes in the TM tilt angle and cytoplasmic dynamics of PLN caused by these mutations correlate well with the extent of SERCA inhibition. Our study unveils new molecular determinants for designing variants of PLN that outcompete endogenous PLN to regulate SERCA in a tunable manner.

The role of structural dynamics in the thermal adaptation of hyperthermophilic enzymes

Proteins from hyperthermophilic organisms are evolutionary optimised to adopt functional structures and dynamics under conditions in which their mesophilic homologues are generally inactive or unfolded. Understanding the nature of such adaptation is of crucial interest to clarify the underlying mechanisms of biological activity in proteins. Here we measured NMR residual dipolar couplings of a hyperthermophilic acylphosphatase enzyme at 80°C and used these data to generate an accurate structural ensemble representative of its native state. The resulting energy landscape was compared to that obtained for a human homologue at 37°C, and additional NMR experiments were carried out to probe fast (¹⁵N relaxation) and slow (H/D exchange) backbone dynamics, collectively sampling fluctuations of the two proteins ranging from the nanosecond to the millisecond timescale. The results identified key differences in the strategies for protein-protein and protein-ligand interactions of the two enzymes at the respective physiological temperatures. These include the dynamical behaviour of a β-strand involved in the protection against aberrant protein aggregation and concerted motions of loops involved in substrate binding and catalysis. Taken together these results elucidate the structure-dynamics-function relationship associated with the strategies of thermal adaptation of protein molecules.

Structures of PKA-phospholamban complexes reveal a mechanism of familial dilated cardiomyopathy

Several mutations identified in phospholamban (PLN) have been linked to familial dilated cardiomyopathy (DCM) and heart failure, yet the underlying molecular mechanism remains controversial. PLN interacts with sarco/endoplasmic reticulum Ca²⁺-ATPase (SERCA) and regulates calcium uptake, which is modulated by the protein kinase A (PKA)-dependent phosphorylation of PLN during the fight-or-flight response. Here, we present the crystal structures of the catalytic domain of mouse PKA in complex with wild-type and DCM-mutant PLNs. Our structures, combined with the results from other biophysical and biochemical assays, reveal a common disease mechanism: the mutations in PLN reduce its phosphorylation level by changing its conformation and weakening its interactions with PKA. In addition, we demonstrate that another more ubiquitous SERCA-regulatory peptide, called another-regulin (ALN), shares a similar mechanism mediated by PKA in regulating SERCA activity.

Structural basis for allosteric control of the SERCA-Phospholamban membrane complex by Ca2+ and phosphorylation

Phospholamban (PLN) is a mini-membrane protein that directly controls the cardiac Ca²⁺-transport response to β-adrenergic stimulation, thus modulating cardiac output during the fight-or-flight response. In the sarcoplasmic reticulum membrane, PLN binds to the sarco(endo)plasmic reticulum Ca²⁺-ATPase (SERCA), keeping this enzyme's function within a narrow physiological window. PLN phosphorylation by cAMP-dependent protein kinase A or increase in Ca²⁺ concentration reverses the inhibitory effects through an unknown mechanism. Using oriented-sample solid-state NMR spectroscopy and replica-averaged NMR-restrained structural refinement, we reveal that phosphorylation of PLN's cytoplasmic regulatory domain signals the disruption of several inhibitory contacts at the transmembrane binding interface of the SERCA-PLN complex that are propagated to the enzyme's active site, augmenting Ca²⁺ transport. Our findings address long-standing questions about SERCA regulation, epitomizing a signal transduction mechanism operated by posttranslationally modified bitopic membrane proteins.

Atomistic Structure and Dynamics of the Ca2+-ATPase Bound to Phosphorylated Phospholamban

Sarcoplasmic reticulum Ca²⁺-ATPase (SERCA) and phospholamban (PLB) are essential components of the cardiac Ca²⁺ transport machinery. PLB phosphorylation at residue Ser16 (pSer16) enhances SERCA activity in the heart via an unknown structural mechanism. Here, we report a fully atomistic model of SERCA bound to phosphorylated PLB and study its structural dynamics on the microsecond time scale using all-atom molecular dynamics simulations in an explicit lipid bilayer and water environment. The unstructured N-terminal phosphorylation domain of PLB samples different orientations and covers a broad area of the cytosolic domain of SERCA but forms a stable complex mediated by pSer16 interactions with a binding site formed by SERCA residues Arg324/Lys328. PLB phosphorylation does not affect the interaction between the transmembrane regions of the two proteins; however, pSer16 stabilizes a disordered structure of the N-terminal phosphorylation domain that releases key inhibitory contacts between SERCA and PLB. We found that PLB phosphorylation is sufficient to guide the structural transitions of the cytosolic headpiece that are required to produce a competent structure of SERCA. We conclude that PLB phosphorylation serves as an allosteric molecular switch that releases inhibitory contacts and strings together the catalytic elements required for SERCA activation. This atomistic model represents a vivid atomic-resolution visualization of SERCA bound to phosphorylated PLB and provides previously inaccessible insights into the structural mechanism by which PLB phosphorylation releases SERCA inhibition in the heart.

Probing membrane protein ground and conformationally excited states using dipolar- and J-coupling mediated MAS solid state NMR experiments

The intrinsic conformational plasticity of membrane proteins directly influences the magnitude of the orientational-dependent NMR interactions such as dipolar couplings (DC) and chemical shift anisotropy (CSA). As a result, the conventional cross-polarization (CP)-based techniques mainly capture the more rigid regions of membrane proteins, while the most dynamic regions are essentially invisible. Nonetheless, dynamic regions can be detected using experiments in which polarization transfer takes place via J-coupling interactions. Here, we review our recent efforts to develop single and dual acquisition pulse sequences with either 1H or 13C detection that utilize both DC and J-coupling mediated transfer to detect both rigid and mobile regions of membrane proteins in native-like lipid environments. We show the application of these new methods for studying the conformational equilibrium of a single-pass membrane protein, phospholamban, which regulates the calcium transport across the sarcoplasmic reticulum (SR) membrane by interacting with the SR Ca2+-ATPase. We anticipate that these methods will be ideal to portray the complex dynamics of membrane proteins in their native environments.

Determination of the conformational states of strychnine in solution using NMR residual dipolar couplings in a tensor-free approach

Small molecules with rotatable bonds can occupy different conformational states in solution as a consequence of their thermal fluctuations. The accurate determination of the structures of such states, as well as of their statistical weights, has been challenging because of the technical difficulties in extracting information from experimental measurements, which are normally averaged over the conformational space available. Here, to achieve this objective, we present an approach based on a recently proposed tensor-free method for incorporating NMR residual dipolar couplings as structural restraints in replica-averaged molecular dynamics simulations. This approach enables the information provided by the experimental data to be used in the spirit of the maximum entropy principle to determine the structural ensembles of small molecules. Furthermore, in order to enhance the sampling of the conformational space we incorporated the metadynamics method in the simulations. We illustrate the method in the case of strychnine, determining the three major conformational states of this small molecule and their associated occupation probabilities.

Structural basis for relief of phospholamban-mediated inhibition of the sarcoplasmic reticulum Ca2+-ATPase at saturating Ca2+ conditions

Sarcoplasmic reticulum Ca²⁺-ATPase (SERCA) is critical for cardiac Ca²⁺ transport. Reversal of phospholamban (PLB)-mediated SERCA inhibition by saturating Ca²⁺ conditions operates as a physiological rheostat to reactivate SERCA function in the absence of PLB phosphorylation. Here, we performed extensive atomistic molecular dynamics simulations to probe the structural mechanism of this process. Simulation of the inhibitory complex at superphysiological Ca²⁺ concentrations ([Ca²⁺] = 10 mm) revealed that Ca²⁺ ions interact primarily with SERCA and the lipid headgroups, but not with PLB's cytosolic domain or the cytosolic side of the SERCA-PLB interface. At this [Ca²⁺], a single Ca²⁺ ion was translocated from the cytosol to the transmembrane transport sites. We used this Ca²⁺-bound complex as an initial structure to simulate the effects of saturating Ca²⁺ at physiological conditions ([Ca²⁺]_total ≈ 400 μm). At these conditions, ∼30% of the Ca²⁺-bound complexes exhibited structural features consistent with an inhibited state. However, in ∼70% of the Ca²⁺-bound complexes, Ca²⁺ moved to transport site I, recruited Glu⁷⁷¹ and Asp⁸⁰⁰, and disrupted key inhibitory contacts involving the conserved PLB residue Asn³⁴ Structural analysis showed that Ca²⁺ induces only local changes in interresidue inhibitory interactions, but does not induce repositioning or changes in PLB structural dynamics. Upon relief of SERCA inhibition, Ca²⁺ binding produced a site I configuration sufficient for subsequent SERCA activation. We propose that at saturating [Ca²⁺] and in the absence of PLB phosphorylation, binding of a single Ca²⁺ ion in the transport sites rapidly shifts the equilibrium toward a noninhibited SERCA-PLB complex.

Effects of the Arg9Cys and Arg25Cys mutations on phospholamban's conformational equilibrium in membrane bilayers

Approximately, 70% of the Ca2+ ion transport into the sarcoplasmic reticulum is catalyzed by the sarcoplasmic reticulum Ca2+-ATPase (SERCA), whose activity is endogenously regulated by phospholamban (PLN). PLN comprises a TM inhibitory region and a cytoplasmic regulatory region that harbors a consensus sequence for cAMP-dependent protein kinase (PKA). The inhibitory region binds the ATPase, reducing its apparent Ca2+ binding affinity. β-adrenergic stimulation activates PKA, which phosphorylates PLN at Ser 16, reversing its inhibitory function. Mutations and post-translational modifications of PLN may lead to dilated cardiomyopathy (DCM) and heart failure. PLN's cytoplasmic region interconverts between a membrane-associated T state and a membrane-detached R state. The importance of these structural transitions on SERCA regulation is emerging, but the effects of natural occurring mutations and their relevance to the progression of heart disease are unclear. Here we use solid-state NMR spectroscopy to investigate the structural dynamics of two lethal PLN mutations, R9C and R25C, which lead to DCM. We found that the R25C mutant enhances the dynamics of PLN and shifts the conformational equilibrium toward the R state confirmation, whereas the R9C mutant drives the amphipathic cytoplasmic domain toward the membrane-associate state, enriching the T state population. The changes in membrane interactions caused by these mutations may explain the aberrant regulation of SERCA.

Electrostatic Interactions in Protein Structure, Folding, Binding, and Condensation

Charged and polar groups, through forming ion pairs, hydrogen bonds, and other less specific electrostatic interactions, impart important properties to proteins. Modulation of the charges on the amino acids, e.g., by pH and by phosphorylation and dephosphorylation, have significant effects such as protein denaturation and switch-like response of signal transduction networks. This review aims to present a unifying theme among the various effects of protein charges and polar groups. Simple models will be used to illustrate basic ideas about electrostatic interactions in proteins, and these ideas in turn will be used to elucidate the roles of electrostatic interactions in protein structure, folding, binding, condensation, and related biological functions. In particular, we will examine how charged side chains are spatially distributed in various types of proteins and how electrostatic interactions affect thermodynamic and kinetic properties of proteins. Our hope is to capture both important historical developments and recent experimental and theoretical advances in quantifying electrostatic contributions of proteins.

Accurate Determination of Conformational Transitions in Oligomeric Membrane Proteins

The structural dynamics governing collective motions in oligomeric membrane proteins play key roles in vital biomolecular processes at cellular membranes. In this study, we present a structural refinement approach that combines solid-state NMR experiments and molecular simulations to accurately describe concerted conformational transitions identifying the overall structural, dynamical, and topological states of oligomeric membrane proteins. The accuracy of the structural ensembles generated with this method is shown to reach the statistical error limit, and is further demonstrated by correctly reproducing orthogonal NMR data. We demonstrate the accuracy of this approach by characterising the pentameric state of phospholamban, a key player in the regulation of calcium uptake in the sarcoplasmic reticulum, and by probing its dynamical activation upon phosphorylation. Our results underline the importance of using an ensemble approach to characterise the conformational transitions that are often responsible for the biological function of oligomeric membrane protein states.

Atomic-level mechanisms for phospholamban regulation of the calcium pump

We performed protein pKa calculations and molecular dynamics (MD) simulations of the calcium pump (sarcoplasmic reticulum Ca(2+)-ATPase (SERCA)) in complex with phospholamban (PLB). X-ray crystallography studies have suggested that PLB locks SERCA in a low-Ca(2+)-affinity E2 state that is incompatible with metal-ion binding, thereby blocking the conversion toward a high-Ca(2+)-affinity E1 state. Estimation of pKa values of the acidic residues in the transport sites indicates that at normal intracellular pH (7.1-7.2), PLB-bound SERCA populates an E1 state that is deprotonated at residues E309 and D800 yet protonated at residue E771. We performed three independent microsecond-long MD simulations to evaluate the structural dynamics of SERCA-PLB in a solution containing 100 mM K(+) and 3 mM Mg(2+). Principal component analysis showed that PLB-bound SERCA lies exclusively along the structural ensemble of the E1 state. We found that the transport sites of PLB-bound SERCA are completely exposed to the cytosol and that K(+) ions bind transiently (≤5 ns) and nonspecifically (nine different positions) to the two transport sites, with a total occupancy time of K(+) in the transport sites of 80%. We propose that PLB binding to SERCA populates a novel (to our knowledge) E1 intermediate, E1⋅H(+)771. This intermediate serves as a kinetic trap that controls headpiece dynamics and depresses the structural transitions necessary for Ca(2+)-dependent activation of SERCA. We conclude that PLB-mediated regulation of SERCA activity in the heart results from biochemical and structural transitions that occur primarily in the E1 state of the pump.

E2 superfamily of ubiquitin-conjugating enzymes: constitutively active or activated through phosphorylation in the catalytic cleft

Protein phosphorylation is a modification that offers a dynamic and reversible mechanism to regulate the majority of cellular processes. Numerous diseases are associated with aberrant regulation of phosphorylation-induced switches. Phosphorylation is emerging as a mechanism to modulate ubiquitination by regulating key enzymes in this pathway. The molecular mechanisms underpinning how phosphorylation regulates ubiquitinating enzymes, however, are elusive. Here, we show the high conservation of a functional site in E2 ubiquitin-conjugating enzymes. In catalytically active E2s, this site contains aspartate or a phosphorylatable serine and we refer to it as the conserved E2 serine/aspartate (CES/D) site. Molecular simulations of substrate-bound and -unbound forms of wild type, mutant and phosphorylated E2s, provide atomistic insight into the role of the CES/D residue for optimal E2 activity. Both the size and charge of the side group at the site play a central role in aligning the substrate lysine toward E2 catalytic cysteine to control ubiquitination efficiency. The CES/D site contributes to the fingerprint of the E2 superfamily. We propose that E2 enzymes can be divided into constitutively active or regulated families. E2s characterized by an aspartate at the CES/D site signify constitutively active E2s, whereas those containing a serine can be regulated by phosphorylation.

The influence of cholesterol on membrane protein structure, function, and dynamics studied by molecular dynamics simulations

The plasma membrane, which encapsulates human cells, is composed of a complex mixture of lipids and embedded proteins. Emerging knowledge points towards the lipids as having a regulating role in protein function. Furthermore, insight from protein crystallography has revealed several different types of lipids intimately bound to membrane proteins and peptides, hereby possibly pointing to a site of action for the observed regulation. Cholesterol is among the lipid membrane constituents most often observed to be co-crystallized with membrane proteins, and the cholesterol levels in cell membranes have been found to play an essential role in health and disease. Remarkably little is known about the mechanism of lipid regulation of membrane protein function in health as well as in disease. Herein, we review molecular dynamics simulation studies aimed at investigating the effect of cholesterol on membrane protein and peptide properties. This article is part of a Special Issue entitled: Lipid-protein interactions.

Structural dynamics and conformational equilibria of SERCA regulatory proteins in membranes by solid-state NMR restrained simulations

Solid-state NMR spectroscopy is emerging as a powerful approach to determine structure, topology, and conformational dynamics of membrane proteins at the atomic level. Conformational dynamics are often inferred and quantified from the motional averaging of the NMR parameters. However, the nature of these motions is difficult to envision based only on spectroscopic data. Here, we utilized restrained molecular dynamics simulations to probe the structural dynamics, topology and conformational transitions of regulatory membrane proteins of the calcium ATPase SERCA, namely sarcolipin and phospholamban, in explicit lipid bilayers. Specifically, we employed oriented solid-state NMR data, such as dipolar couplings and chemical shift anisotropy measured in lipid bicelles, to refine the conformational ensemble of these proteins in lipid membranes. The samplings accurately reproduced the orientations of transmembrane helices and showed a significant degree of convergence with all of the NMR parameters. Unlike the unrestrained simulations, the resulting sarcolipin structures are in agreement with distances and angles for hydrogen bonds in ideal helices. In the case of phospholamban, the restrained ensemble sampled the conformational interconversion between T (helical) and R (unfolded) states for the cytoplasmic region that could not be observed using standard structural refinements with the same experimental data set. This study underscores the importance of implementing NMR data in molecular dynamics protocols to better describe the conformational landscapes of membrane proteins embedded in realistic lipid membranes.

Phospholamban phosphorylation, mutation, and structural dynamics: a biophysical approach to understanding and treating cardiomyopathy

We review the recent development of novel biochemical and spectroscopic methods to determine the site-specific phosphorylation, expression, mutation, and structural dynamics of phospholamban (PLB), in relation to its function (inhibition of the cardiac calcium pump, SERCA2a), with specific focus on cardiac physiology, pathology, and therapy. In the cardiomyocyte, SERCA2a actively transports Ca²⁺ into the sarcoplasmic reticulum (SR) during relaxation (diastole) to create the concentration gradient that drives the passive efflux of Ca²⁺ required for cardiac contraction (systole). Unphosphorylated PLB (U-PLB) inhibits SERCA2a, but phosphorylation at S16 and/or T17 (producing P-PLB) changes the structure of PLB to relieve SERCA2a inhibition. Because insufficient SERCA2a activity is a hallmark of heart failure, SERCA2a activation, by gene therapy (Andino et al. 2008; Fish et al. 2013; Hoshijima et al. 2002; Jessup et al. 2011) or drug therapy (Ferrandi et al. 2013; Huang 2013; Khan et al. 2009; Rocchetti et al. 2008; Zhang et al. 2012), is a widely sought goal for treatment of heart failure. This review describes rational approaches to this goal. Novel biophysical assays, using site-directed labeling and high-resolution spectroscopy, have been developed to resolve the structural states of SERCA2a-PLB complexes in vitro and in living cells. Novel biochemical assays, using synthetic standards and multidimensional immunofluorescence, have been developed to quantitate PLB expression and phosphorylation states in cells and human tissues. The biochemical and biophysical properties of U-PLB, P-PLB, and mutant PLB will ultimately resolve the mechanisms of loss of inhibition and gain of inhibition to guide therapeutic development. These assays will be powerful tools for investigating human tissue samples from the Sydney Heart Bank, for the purpose of analyzing and diagnosing specific disorders.

Structure of a low-population intermediate state in the release of an enzyme product

Enzymes can increase the rate of biomolecular reactions by several orders of magnitude. Although the steps of substrate capture and product release are essential in the enzymatic process, complete atomic-level descriptions of these steps are difficult to obtain because of the transient nature of the intermediate conformations, which makes them largely inaccessible to standard structure determination methods. We describe here the determination of the structure of a low-population intermediate in the product release process by human lysozyme through a combination of NMR spectroscopy and molecular dynamics simulations. We validate this structure by rationally designing two mutations, the first engineered to destabilise the intermediate and the second to stabilise it, thus slowing down or speeding up, respectively, product release. These results illustrate how product release by an enzyme can be facilitated by the presence of a metastable intermediate with transient weak interactions between the enzyme and product.

Conformational recognition of an intrinsically disordered protein

There is a growing interest in understanding the properties of intrinsically disordered proteins (IDPs); however, the characterization of these states remains an open challenge. IDPs appear to have functional roles that diverge from those of folded proteins and revolve around their ability to act as hubs for protein-protein interactions. To gain a better understanding of the modes of binding of IDPs, we combined statistical mechanics, calorimetry, and NMR spectroscopy to investigate the recognition and binding of a fragment from the disordered protein Gab2 by the growth factor receptor-bound protein 2 (Grb2), a key interaction for normal cell signaling and cancer development. Structural ensemble refinement by NMR chemical shifts, thermodynamics measurements, and analysis of point mutations indicated that the population of preexisting bound conformations in the free-state ensemble of Gab2 is an essential determinant for recognition and binding by Grb2. A key role was found for transient polyproline II (PPII) structures and extended conformations. Our findings are likely to have very general implications for the biological behavior of IDPs in light of the evidence that a large fraction of these proteins possess a specific propensity to form PPII and to adopt conformations that are more extended than the typical random-coil states.

A tensor-free method for the structural and dynamical refinement of proteins using residual dipolar couplings

Residual dipolar couplings (RDCs) are parameters measured in nuclear magnetic resonance spectroscopy that can provide exquisitely detailed information about the structure and dynamics of biological macromolecules. We describe here a method of using RDCs for the structural and dynamical refinement of proteins that is based on the observation that the RDC between two atomic nuclei depends directly on the angle ϑ between the internuclear vector and the external magnetic field. For every pair of nuclei for which an RDC is available experimentally, we introduce a structural restraint to minimize the deviation from the value of the angle ϑ derived from the measured RDC and that calculated in the refinement protocol. As each restraint involves only the calculation of the angle ϑ of the corresponding internuclear vector, the method does not require the definition of an overall alignment tensor to describe the preferred orientation of the protein with respect to the alignment medium. Application to the case of ubiquitin demonstrates that this method enables an accurate refinement of the structure and dynamics of this protein to be obtained.

New opportunities for tensor-free calculations of residual dipolar couplings for the study of protein dynamics

Residual dipolar couplings (RDCs) can provide exquisitely detailed information about the structure and dynamics of proteins. It is challenging, however, to extract such information from RDC measurements in conformationally heterogeneous states of proteins because of the complex relationship between RDCs and protein structures. To obtain new insights into this problem, we discuss methods of calculating the RDCs that do not require the definition of an alignment tensor. These methods can help in particular in the search of effective ways to use RDCs to characterise disordered or partially disordered states of proteins.

Current and emerging opportunities for molecular simulations in structure-based drug design

An overview of the current capabilities and limitations of molecular simulation of biomolecular complexes in the context of computer-aided drug design is provided. Steady improvements in computer hardware coupled with more refined representations of energetics are leading to a new appreciation of the driving forces of molecular recognition. Molecular simulations are poised to more frequently guide the interpretation of biophysical measurements of biomolecular complexes. Ligand design strategies emerge from detailed analyses of computed structural ensembles. The feasibility of routine applications to ligand optimization problems hinges upon successful extensive large scale validation studies and the development of protocols to intelligently automate computations.