The alpha-helical propensity of the cytoplasmic domain of phospholamban: a molecular dynamics simulation of the effect of phosphorylation and mutation

Thirty years of molecular dynamics simulations on posttranslational modifications of proteins

Posttranslational modifications (PTMs) are an integral component to how cells respond to perturbation. While experimental advances have enabled improved PTM identification capabilities, the same throughput for characterizing how structural changes caused by PTMs equate to altered physiological function has not been maintained. In this Perspective, we cover the history of computational modeling and molecular dynamics simulations which have characterized the structural implications of PTMs. We distinguish results from different molecular dynamics studies based upon the timescales simulated and analysis approaches used for PTM characterization. Lastly, we offer insights into how opportunities for modern research efforts on in silico PTM characterization may proceed given current state-of-the-art computing capabilities and methodological advancements.

Molecular dynamics shows complex interplay and long-range effects of post-translational modifications in yeast protein interactions

Post-translational modifications (PTMs) play a vital, yet often overlooked role in the living cells through modulation of protein properties, such as localization and affinity towards their interactors, thereby enabling quick adaptation to changing environmental conditions. We have previously benchmarked a computational framework for the prediction of PTMs’ effects on the stability of protein-protein interactions, which has molecular dynamics simulations followed by free energy calculations at its core. In the present work, we apply this framework to publicly available data on Saccharomyces cerevisiae protein structures and PTM sites, identified in both normal and stress conditions. We predict proteome-wide effects of acetylations and phosphorylations on protein-protein interactions and find that acetylations more frequently have locally stabilizing roles in protein interactions, while the opposite is true for phosphorylations. However, the overall impact of PTMs on protein-protein interactions is more complex than a simple sum of local changes caused by the introduction of PTMs and adds to our understanding of PTM cross-talk. We further use the obtained data to calculate the conformational changes brought about by PTMs. Finally, conservation of the analyzed PTM residues in orthologues shows that some predictions for yeast proteins will be mirrored to other organisms, including human. This work, therefore, contributes to our overall understanding of the modulation of the cellular protein interaction networks in yeast and beyond.

Proteins are a diverse set of biological molecules responsible for numerous functions within cells, such as obtaining energy from food or transport of small molecules, and many processes rely on interactions of specific proteins. Moreover, a single protein may acquire different roles depending on cellular requirements and as a response to changes in the environment. A commonly used way to quickly change protein’s function or activity is by introducing small chemical modifications on specific locations within the protein. These modifications can cause the protein to interact in a more or less stable way with other proteins. We have previously developed a computational pipeline for predicting the effect of modifications on interactions of proteins, and in this work we apply it to all yeast proteins with known structures. We find differences in effects on the binding for different types of modifications. Importantly, we demonstrate that the modifications far from the interaction interface also significantly contribute to binding due to their impact on protein’s shape, which is often neglected by other methods. This work contributes to our understanding of the modulation of protein interactions in yeast due to modifications, while our widely applicable method will allow similar investigations in other organisms.

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.

Crystal structure of diamondback moth ryanodine receptor Repeat34 domain reveals insect-specific phosphorylation sites

Background

Ryanodine receptor (RyR), a calcium-release channel located in the sarcoplasmic reticulum membrane of muscles, is the target of insecticides used against a wide range of agricultural pests. Mammalian RyRs have been shown to be under the regulatory control of several kinases and phosphatases, but little is known about the regulation of insect RyRs by phosphorylation.

Results

Here we present the crystal structures of wild-type and phospho-mimetic RyR Repeat34 domain containing PKA phosphorylation sites from diamondback moth (DBM), a major lepidopteran pest of cruciferous vegetables. The structure has unique features, not seen in mammalian RyRs, including an additional α-helix near the phosphorylation loop. Using tandem mass spectrometry, we identify several PKA sites clustering in the phosphorylation loop and the newly identified α-helix. Bioinformatics analysis shows that this α-helix is only present in Lepidoptera, suggesting an insect-specific regulation. Interestingly, the specific phosphorylation pattern is temperature-dependent. The thermal stability of the DBM Repeat34 domain is significantly lower than that of the analogous domain in the three mammalian RyR isoforms, indicating a more dynamic domain structure that can be partially unfolded to facilitate the temperature-dependent phosphorylation. Docking the structure into the cryo-electron microscopy model of full-length RyR reveals that the interface between the Repeat34 and neighboring HD1 domain is more conserved than that of the phosphorylation loop region that might be involved in the interaction with SPRY3 domain. We also identify an insect-specific glycerol-binding pocket that could be potentially targeted by novel insecticides to fight the current resistance crisis.

Conclusions

The crystal structures of the DBM Repeat34 domain reveals insect-specific temperature-dependent phosphorylation sites that may regulate insect ryanodine receptor function. It also reveals insect-specific structural features and a potential ligand-binding site that could be targeted in an effort to develop green pesticides with high species-specificity.

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.

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.

Exploring the association of rs10490924 polymorphism with age-related macular degeneration: An in silico approach

The polymorphism rs10490924 (A69S) in the age-related maculopathy susceptibility 2 (ARMS2) gene is highly associated with age-related macular degeneration, which is the leading cause of blindness among the elderly population. ARMS2 gene encodes a putative small (11 kDa) protein, which the function and localization of the ARMS2 protein remain under debate. For a better understanding of functional impacts of A69S mutation, we performed a detailed analysis of an ARMS2 sequence with a broad set of bioinformatics tools. In silico analysis was followed to predict the tertiary structure, putative binding site regions, and binding site residues. Also, the effects of this mutation on protein stability, aggregation propensity, and homodimerization were analyzed. Next, a molecular dynamic simulation was carried out to understand the dynamic behavior of wild-type, A69S, and phosphorylated A69S structures. The results showed alterations in the putative post-translational modification sites on the ARMS2 protein, due to the mutation. Furthermore, the stability of protein and putative homodimer conformations were affected by the mutation. Molecular dynamic simulation results revealed that A69S mutation enhances the rigidity of the ARMS2 structure and residue serine at position 69 is buried and may not be phosphorylated; however, phosphorylated serine enhances the flexibility of the ARMS2 structure. In conclusion, our study provides new insights into the deleterious effects of A69S mutation on the ARMS2 structure.

Exploring the association of rs10490924 polymorphism with age-related macular degeneration: An in silico approach

The polymorphism rs10490924 (A69S) in the age-related maculopathy susceptibility 2 (ARMS2) gene is highly associated with age-related macular degeneration, which is the leading cause of blindness among the elderly population. The ARMS2 gene encodes a putative small (11kDa) protein, which the function and localization of the ARMS2 protein remain under debate. For a better understanding of functional impacts of the A69S mutation, we performed a detailed analysis of the ARMS2 sequence with a broad set of bioinformatics tools. In silico analysis was followed to predict the tertiary structure, putative binding site regions, and binding site residues. Also, the effects of this mutation on protein stability, aggregation propensity, and homodimerization were analyzed. Next, a molecular dynamic simulation was carried out to understand the dynamic behavior of wild-type, A69S, and phosphorylated A69S structures. The results showed alterations in the putative post-translational modification sites on the ARMS2 protein, due to the mutation. Furthermore, the stability of protein and putative homodimer conformations were affected by the mutation. Molecular dynamic simulation results revealed that the A69S mutation enhances the rigidity of the ARMS2 structure and residue serine at position 69 is buried and may not be phosphorylated; however, phosphorylated serine enhances the flexibility of the ARMS2 structure. In conclusion, our study provides new insights into the deleterious effects of the A69S mutation on the ARMS2 structure.

A phosphorylation-motif for tuneable helix stabilisation in intrinsically disordered proteins - Lessons from the sodium proton exchanger 1 (NHE1)

Intrinsically disordered proteins (IDPs) are involved in many pivotal cellular processes including phosphorylation and signalling. The structural and functional effects of phosphorylation of IDPs remain poorly understood and difficult to predict. Thus, a need exists to identify motifs that confer phosphorylation-dependent perturbation of the local preferences for forming e.g. helical structures as well as motifs that do not. The disordered distal tail of the Na⁺/H⁺ exchanger 1 (NHE1) is six-times phosphorylated (S693, S723, S726, S771, T779, S785) by the mitogen activated protein kinase 2 (MAPK1, ERK2). Using NMR spectroscopy, we found that two out of those six phosphorylation sites had a stabilizing effect on transient helices. One of these was further investigated by circular dichroism and NMR spectroscopy as well as by molecular dynamic simulations, which confirmed the stabilizing effect and resulted in the identification of a short linear motif for helix stabilisation: [S/T]-P-{3}-[R/K] where [S/T] is the phosphorylation-site. By analysing IDP and phosphorylation site databases we found that the motif is significantly enriched around known phosphorylation sites, supporting a potential wider-spread role in phosphorylation-mediated regulation of intrinsically disordered proteins. The identification of such motifs is important for understanding the molecular mechanism of cellular signalling, and is crucial for the development of predictors for the structural effect of phosphorylation; a tool of relevance for understanding disease-promoting mutations that for example interfere with signalling for instance through constitutive active and often cancer-promoting signalling.

Characterization of the Gly45Asp variant of human cytochrome P450 1A1 using recombinant expression

Cytochrome P450 1A1 (CYP1A1) is a heme-containing enzyme involved in metabolism of xenobiotics. CYP1A1 containing a Gly45Asp substitution has not yet been characterized. Escherichia coli expressing the Gly45Asp variant, as well as the purified variant protein, had lower CYP (i.e., holoenzyme) contents than their wild-type (WT) equivalents. The purified variant protein had reduced heme contents compared with their WT equivalents. Enhanced supplementation of a heme precursor during culture did not increase CYP content in E. coli expressing the variant, but did for the WT. Substitution of Gly45 with other residues, especially those having large side chains, decreased CYP contents of E. coli expressing the variants to a considerable extent. A 3D structure of CYP1A1 indicates that Gly45, along with other residues of the PR region, interacts with Arg77 of β- strand 1-1, which indirectly interacts with heme. Substitution analyses suggest the importance of residues of the PR region and Arg77 in holoenzyme expression. E. coli membrane and mammalian microsomes expressing the Gly45Asp variant, as well as the purified variant protein, had reduced ethoxyresorufin O-dealkylation activities, compared with the WT equivalents. These findings suggest the Gly45Asp substitution results in a structural disturbance of CYP1A1, reducing its holoenzyme formation and catalytic activity.

Acute inotropic and lusitropic effects of cardiomyopathic R9C mutation of phospholamban

A naturally occurring R9C mutation of phospholamban (PLB) triggers cardiomyopathy and premature death by altering regulation of sarco/endoplasmic reticulum calcium-ATPase (SERCA). The goal of this study was to investigate the acute physiological consequences of the R9C-PLB mutation on cardiomyocyte calcium kinetics and contractility. We measured the physiological consequences of R9C-PLB mutation on calcium transients and sarcomere shortening in adult cardiomyocytes. In contrast to studies of chronic R9C-PLB expression in transgenic mice, we found that acute expression of R9C-PLB exerts a positively inotropic and lusitropic effect in cardiomyocytes. Importantly, R9C-PLB exhibited blunted sensitivity to frequency potentiation and β-adrenergic stimulation, two major physiological mechanisms for the regulation of cardiac performance. To identify the molecular mechanism of R9C pathology, we quantified the effect of R9C on PLB oligomerization and PLB-SERCA binding. FRET measurements in live cells revealed that R9C-PLB exhibited an increased propensity for oligomerization, and this was further increased by oxidative stress. The R9C also decreased PLB binding to SERCA and altered the structure of the PLB-SERCA regulatory complex. The structural change after oxidative modification of R9C-PLB was similar to that observed after PLB phosphorylation. We conclude that R9C mutation of PLB decreases SERCA inhibition by decreasing the amount of the regulatory complex and altering its conformation. This has an acute inotropic/lusitropic effect but yields negative consequences of impaired frequency potentiation and blunted β-adrenergic responsiveness. We envision a self-reinforcing mechanism beginning with phosphomimetic R9C-PLB oxidation and loss of SERCA inhibition, leading to impaired calcium regulation and heart failure.

Molecular dynamics simulation of the phosphorylation-induced conformational changes of a tau peptide fragment

Aggregation of the microtubule associated protein tau (MAPT) within neurons of the brain is the leading cause of tauopathies such as Alzheimer's disease. MAPT is a phospho-protein that is selectively phosphorylated by a number of kinases in vivo to perform its biological function. However, it may become pathogenically hyperphosphorylated, causing aggregation into paired helical filaments and neurofibrillary tangles. The phosphorylation induced conformational change on a peptide of MAPT (htau225-250) was investigated by performing molecular dynamics simulations with different phosphorylation patterns of the peptide (pThr231 and/or pSer235) in different simulation conditions to determine the effect of ionic strength and phosphate charge. All phosphorylation patterns were found to disrupt a nascent terminal β-sheet pattern (226VAVVR230 and 244QTAPVP249), replacing it with a range of structures. The double pThr231/pSer235 phosphorylation pattern at experimental ionic strength resulted in the best agreement with NMR structural characterization, with the observation of a transient α-helix (239AKSRLQT245). PPII helical conformations were only found sporadically throughout the simulations.

Secondary structure, backbone dynamics, and structural topology of phospholamban and its phosphorylated and Arg9Cys-mutated forms in phospholipid bilayers utilizing 13C and 15N solid-state NMR spectroscopy

Phospholamban (PLB) is a membrane protein that regulates heart muscle relaxation rates via interactions with the sarcoplasmic reticulum Ca(2+) ATPase (SERCA). When PLB is phosphorylated or Arg9Cys (R9C) is mutated, inhibition of SERCA is relieved. (13)C and (15)N solid-state NMR spectroscopy is utilized to investigate conformational changes of PLB upon phosphorylation and R9C mutation. (13)C═O NMR spectra of the cytoplasmic domain reveal two α-helical structural components with population changes upon phosphorylation and R9C mutation. The appearance of an unstructured component is observed on domain Ib. (15)N NMR spectra indicate an increase in backbone dynamics of the cytoplasmic domain. Wild-type PLB (WT-PLB), Ser16-phosphorylated PLB (P-PLB), and R9C-mutated PLB (R9C-PLB) all have a very dynamic domain Ib, and the transmembrane domain has an immobile component. (15)N NMR spectra indicate that the cytoplasmic domain of R9C-PLB adopts an orientation similar to P-PLB and shifts away from the membrane surface. Domain Ib (Leu28) of P-PLB and R9C-PLB loses the alignment. The R9C-PLB adopts a conformation similar to P-PLB with a population shift to a more extended and disordered state. The NMR data suggest the more extended and disordered forms of PLB may relate to inhibition relief.

Effect of Ser392 phosphorylation on the structure and dynamics of the polybasic domain of ADP ribosylation factor nucleotide site opener protein: a molecular simulation study

ADP ribosylation factor nucleotide site opener (ARNO) as a guanine nucleotide exchange factor (GEF) activates small GTPases called ADP ribosylation factors (Arfs), which function as molecular switches and regulate a variety of cell biological events. ARNO directly interacts with the transmembrane a2-subunit isoform of the proton-pumping vacuolar ATPase in an acidification-dependent manner, and this interaction plays a crucial role in the regulation of the protein degradation pathway. A recent study reported specific interactions of a2N with the ARNO375-400 peptide corresponding to the polybasic (PB) domain of ARNO, which is a crucial regulatory element in the autoregulation and modulation of Arf-GEF activity. Interestingly, phosphorylation of Ser392 completely abolishes this interaction, and the experimental structure shows significant structural rearrangements. To investigate the effect of Ser392 phosphorylation on the structure and dynamics of the ARNO375-400 peptide, we employed all atom molecular dynamics (MD) simulations of the phosphorylated and unphosphorylated PB domain of the ARNO protein. A Hamiltonian-based replica exchange method called biasing potential replica exchange MD was used to enhance conformational sampling. Simulations predicted that the isolated PB domain is highly flexible, with the C-terminal region of the unphosphorylated state being unstable. In contrast, Ser392 phosphorylation increases the overall stability of the peptide. In agreement with experimental results, our simulations further support the hypothesis that phosphorylation induces disorder to order transitions and provide new insights into the structural dynamics of the PB domain. Phosphorylation of Ser392 appears to stabilize the C-terminal α-helix via formation of salt bridges between phospho-Ser392 and Arg390, Lys395, and Lys396.

Probing the interaction of Arg9Cys mutated phospholamban with phospholipid bilayers by solid-state NMR spectroscopy

Phospholamban (PLB) is a 52 amino acid integral membrane protein that interacts with the sarcoplasmic reticulum Ca(2+) ATPase (SERCA) and helps to regulate Ca(2+) flow. PLB inhibits SERCA impairing Ca(2+) translocation. The inhibition can be relieved upon phosphorylation of PLB. The Arg9 to Cys (R9C) mutation is a loss of function mutation with reduced inhibitory potency. The effect R9C PLB has on the membrane surface and the hydrophobic region dynamics was investigated by (31)P and (2)H solid-state NMR spectroscopy in multilamellar vesicles (MLVs). The (31)P NMR spectra indicate that, like the phosphorylated PLB (P-PLB), the mutated R9C-PLB protein has significantly less interaction with the lipid bilayer headgroup when compared to wild-type PLB (WT-PLB). Similar to P-PLB, R9C-PLB slightly decreases (31)P T1 values in the lipid headgroup region. (2)H SCD order parameters of (2)H nuclei along the lipid acyl chain decrease less dramatically for R9C-PLB and P-PLB when compared to WT-PLB. The results suggest that R9C-PLB interacts less with the membrane surface and hydrophobic region than WT-PLB. Detachment of the cytoplasmic domain of R9C-PLB from the membrane surface could be related to its loss of function.

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.

Structural dynamics of muscle protein phosphorylation

We have used site-directed spectroscopic probes to detect structural changes, motions, and interactions due to phosphorylation of proteins involved in the regulation of muscle contraction and relaxation. Protein crystal structures provide static snapshots that provide clues to the conformations that are sampled dynamically by proteins in the cellular environment. Our site-directed spectroscopic experiments, combined with computational simulations, extend these studies into functional assemblies in solution, and reveal details of protein regions that are too dynamic or disordered for crystallographic approaches. Here, we discuss phosphorylation-mediated structural transitions in the smooth muscle myosin regulatory light chain, the striated muscle accessory protein myosin binding protein-C, and the cardiac membrane Ca(2+) pump modulator phospholamban. In each of these systems, phosphorylation near the N terminus of the regulatory protein relieves an inhibitory interaction between the phosphoprotein and its regulatory target. Several additional unifying themes emerge from our studies: (a) The effect of phosphorylation is not to change the affinity of the phosphoprotein for its regulated binding partner, but to change the structure of the bound complex without dissociation. (b) Phosphorylation induces transitions between order and dynamic disorder. (c) Structural states are only loosely coupled to phosphorylation; i.e., complete phosphorylation induces dramatic functional effects with only a partial shift in the equilibrium between ordered and disordered structural states. These studies, which offer atomic-resolution insight into the structural and functional dynamics of these phosphoproteins, were inspired in part by the ground-breaking work in this field by Michael and Kate Barany.

Lethal, hereditary mutants of phospholamban elude phosphorylation by protein kinase A

The sarcoplasmic reticulum calcium pump (SERCA) and its regulator, phospholamban, are essential components of cardiac contractility. Phospholamban modulates contractility by inhibiting SERCA, and this process is dynamically regulated by β-adrenergic stimulation and phosphorylation of phospholamban. Herein we reveal mechanistic insight into how four hereditary mutants of phospholamban, Arg(9) to Cys, Arg(9) to Leu, Arg(9) to His, and Arg(14) deletion, alter regulation of SERCA. Deletion of Arg(14) disrupts the protein kinase A recognition motif, which abrogates phospholamban phosphorylation and results in constitutive SERCA inhibition. Mutation of Arg(9) causes more complex changes in function, where hydrophobic substitutions such as cysteine and leucine eliminate both SERCA inhibition and phospholamban phosphorylation, whereas an aromatic substitution such as histidine selectively disrupts phosphorylation. We demonstrate that the role of Arg(9) in phospholamban function is multifaceted: it is important for inhibition of SERCA, it increases the efficiency of phosphorylation, and it is critical for protein kinase A recognition in the context of the phospholamban pentamer. Given the synergistic consequences on contractility, it is not surprising that the mutants cause lethal, hereditary dilated cardiomyopathy.

Protein Electrostatic Properties Predefining the Level of Surface Hydrophobicity Change upon Phosphorylation

We use explicit-solvent, molecular dynamics simulations to study the change in polar properties of a solvent-accessible surface for proteins undergoing phosphorylation. We analyze eight different pairs of proteins representing different structural classes in native and phosphorylated states and estimate the polarity of their surface using the molecular hydrophobicity potential approach. Whereas the phosphorylation-induced hydrophobicity change in the vicinity of phosphosites does not vary strongly among the studied proteins, the equivalent change for complete proteins covers a surprisingly wide range of effects including even an increase in the overall hydrophobicity in some cases. Importantly, the observed changes are strongly related to electrostatic properties of proteins, such as the net charge per residue, the distribution of charged side-chain contacts, and the isoelectric point. These features predefine the level of surface hydrophobicity change upon phosphorylation and may thus contribute to the phosphorylation-induced alteration of the interactions between a protein and its environment.

An acidic loop and cognate phosphorylation sites define a molecular switch that modulates ubiquitin charging activity in Cdc34-like enzymes

E2 ubiquitin-conjugating enzymes are crucial mediators of protein ubiquitination, which strongly influence the ultimate fate of the target substrates. Recently, it has been shown that the activity of several enzymes of the ubiquitination pathway is finely tuned by phosphorylation, an ubiquitous mechanism for cellular regulation, which modulates protein conformation. In this contribution, we provide the first rationale, at the molecular level, of the regulatory mechanism mediated by casein kinase 2 (CK2) phosphorylation of E2 Cdc34-like enzymes. In particular, we identify two co-evolving signature elements in one of the larger families of E2 enzymes: an acidic insertion in β4α2 loop in the proximity of the catalytic cysteine and two conserved key serine residues within the catalytic domain, which are phosphorylated by CK2. Our investigations, using yeast Cdc34 as a model, through 2.5 µs molecular dynamics simulations and biochemical assays, define these two elements as an important phosphorylation-controlled switch that modulates opening and closing of the catalytic cleft. The mechanism relies on electrostatic repulsions between a conserved serine phosphorylated by CK2 and the acidic residues of the β4α2 loop, promoting E2 ubiquitin charging activity. Our investigation identifies a new and unexpected pivotal role for the acidic loop, providing the first evidence that this loop is crucial not only for downstream events related to ubiquitin chain assembly, but is also mandatory for the modulation of an upstream crucial step of the ubiquitin pathway: the ubiquitin charging in the E2 catalytic cleft.

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.

Phosphorylation and mutation of phospholamban alter physical interactions with the sarcoplasmic reticulum calcium pump

Phospholamban physically interacts with the sarcoplasmic reticulum calcium pump (SERCA) and regulates contractility of the heart in response to adrenergic stimuli. We studied this interaction using electron microscopy of 2D crystals of SERCA in complex with phospholamban. In earlier studies, phospholamban oligomers were found interspersed between SERCA dimer ribbons and a 3D model was constructed to show interactions with SERCA. In this study, we examined the oligomeric state of phospholamban and the effects of phosphorylation and mutation of phospholamban on the interaction with SERCA in the 2D crystals. On the basis of projection maps from negatively stained and frozen-hydrated crystals, phosphorylation of Ser16 selectively disordered the cytoplasmic domain of wild type phospholamban. This was not the case for a pentameric gain-of-function mutant (Lys27Ala), which retained inhibitory activity and remained ordered in the phosphorylated state. A partial loss-of-function mutation that altered the charge state of phospholamban (Arg14Ala) retained an ordered state, while a complete loss-of-function mutation (Asn34Ala) was also disordered. The functional state of phospholamban was correlated with an order-to-disorder transition of the phospholamban cytoplasmic domain in the 2D co-crystals. Furthermore, co-crystals of the gain-of-function mutant (Lys27Ala) facilitated data collection from frozen-hydrated crystals. An improved projection map was calculated to a resolution of 8 Å, which supports the pentamer as the oligomeric state of phospholamban in the crystals. The 2D co-crystals with SERCA require a functional pentameric form of phospholamban, which physically interacts with SERCA at an accessory site distinct from that used by the phospholamban monomer for the inhibitory association.

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.

Phosphorylation-induced structural changes in smooth muscle myosin regulatory light chain

We have performed complementary time-resolved fluorescence resonance energy transfer (TR-FRET) experiments and molecular dynamics (MD) simulations to elucidate structural changes in the phosphorylation domain (PD) of smooth muscle regulatory light chain (RLC) bound to myosin. PD is absent in crystal structures, leaving uncertainty about the mechanism of regulation. Donor-acceptor pairs of probes were attached to three site-directed di-Cys mutants of RLC, each having one Cys at position 129 in the C-terminal lobe and the other at position 2, 3, or 7 in the N-terminal PD. Labeled RLC was reconstituted onto myosin subfragment 1 (S1). TR-FRET resolved two simultaneously populated structural states of RLC, closed and open, in both unphosphorylated and phosphorylated biochemical states. All three FRET pairs show that phosphorylation shifts the equilibrium toward the open state, increasing its mol fraction by approximately 20%. MD simulations agree with experiments in remarkable detail, confirming the coexistence of two structural states, with phosphorylation shifting the system toward the more dynamic open structural state. This agreement between experiment and simulation validates the additional structural details provided by MD simulations: In the closed state, PD is bent onto the surface of the C-terminal lobe, stabilized by interdomain salt bridges. In the open state, PD is more helical and straight, resides farther from the C-terminal lobe, and is stabilized by an intradomain salt bridge. The result is a vivid atomic-resolution visualization of the first step in the molecular mechanism by which phosphorylation activates smooth muscle.

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.

Thermodynamic and structural basis of phosphorylation-induced disorder-to-order transition in the regulatory light chain of smooth muscle myosin

We have performed molecular dynamics simulations of the phosphorylation domain (PD) of the regulatory light chain (RLC) of smooth muscle myosin, to gain insight into the thermodynamic principles governing the phosphorylation-induced disorder-to-order transition. Simulations were performed in explicit water under near-physiological conditions, starting with an ideal alpha-helix. In the absence of phosphorylation, the helical periodicity of the peptide was disrupted at residues T9-K11, while phosphorylation significantly favored the helical periodicity, in agreement with experimental data. Using the MM/PBSA approach, we calculated a relative free energy of -7.1 kcal/mol for the disorder-to-order transition. A large enthalpic decrease was compensated by a large loss of conformational entropy, despite the small helical increase (no more than three residues) upon phosphorylation. Phosphorylation decreased the conformational dynamics of K and R side chains, especially R16, which forms a salt bridge with pS19. Mutation of R16 to A or E prevented this phosphorylation-dependent ordering. We propose that phosphorylation balances the enthalpy-entropy compensation of the disorder-to-order transition of RLC via short and long-range electrostatic interactions with positively charged residues of the phosphorylation domain. We suggest that this balance is necessary to induce a disorder-to-order conformational change through a subtle energy switching.

Synthesis of TOAC spin-labeled proteins and reconstitution in lipid membranes

A procedure is described for the synthetic incorporation into membrane proteins of the non-natural amino acid TOAC (2,2,6,6-tetramethyl-piperidine-1-oxyl-4-amino-4-carboxylic acid), which is coupled rigidly to the alpha-carbon, providing direct detection of peptide backbone dynamics by electron paramagnetic resonance (EPR). Also included is a protocol for the functional reconstitution of the spin-labeled protein in lipid vesicles. This protocol can be completed in 17 d.

Molecular dynamics simulations reveal a disorder-to-order transition on phosphorylation of smooth muscle myosin

We have performed molecular dynamics simulations of the phosphorylated (at S-19) and the unphosphorylated 25-residue N-terminal phosphorylation domain of the regulatory light chain (RLC) of smooth muscle myosin to provide insight into the structural basis of regulation. This domain does not appear in any crystal structure, so these simulations were combined with site-directed spin labeling to define its structure and dynamics. Simulations were carried out in explicit water at 310 K, starting with an ideal alpha-helix. In the absence of phosphorylation, large portions of the domain (residues S-2 to K-11 and R-16 through Y-21) were metastable throughout the simulation, undergoing rapid transitions among alpha-helix, pi-helix, and turn, whereas residues K-12 to Q-15 remained highly disordered, displaying a turn motif from 1 to 22.5 ns and a random coil pattern from 22.5 to 50 ns. Phosphorylation increased alpha-helical order dramatically in residues K-11 to A-17 but caused relatively little change in the immediate vicinity of the phosphorylation site (S-19). Phosphorylation also increased the overall dynamic stability, as evidenced by smaller temporal fluctuations in the root mean-square deviation. These results on the isolated phosphorylation domain, predicting a disorder-to-order transition induced by phosphorylation, are remarkably consistent with published experimental data involving site-directed spin labeling of the intact RLC bound to the two-headed heavy meromyosin. The simulations provide new insight into structural details not revealed by experiment, allowing us to propose a refined model for the mechanism by which phosphorylation affects the N-terminal domain of the RLC of smooth muscle myosin.

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

Phospholamban (PLB) is a small membrane protein that regulates the activity of the calcium ATP-ase in the cardiac, slow-twitch, and smooth muscle sarcoplasmic reticulum through the reversible phosphorylation of Ser16. We present here a comparative molecular dynamics study of unmodified and phosphorylated PLB immersed in a phospholipid membrane. The study has been performed under different ionic strength conditions, using the NMR structures of two PLB variants determined in mixed organic solvent and dodecylphosphocholine micelles. The simulations indicate that all PLB forms studied display a highly dynamic behavior of the N-terminal cytoplasmic moiety, with a decrease of its helical content in the phosphorylated forms. The cytoplasmic domain undergoes large collective motions sampling conformations parallel as well as perpendicular to the membrane surface in all the simulations. The transmembrane domain retains a tightly folded helical conformation with a small tilt with respect to the membrane plane probably induced by the presence of Asn30 and Asn34 within the hydrophobic environment. Furthermore, the phosphoric group on Ser16 establishes transient electrostatic interactions with the phospholipid heads. We propose a model in which phosphorylation diminishes the probability of interactions of PLB with residues near Lys400 in the SERCA pump, thus relieving its inhibition.

Backbone dynamics determined by electron paramagnetic resonance to optimize solid-phase peptide synthesis of TOAC-labeled phospholamban

Electron paramagnetic resonance (EPR) was used to optimize the solid-phase peptide synthesis of a membrane-bound peptide labeled with TOAC (2,2,6,6-tetramethyl-piperidine-1-oxyl-4-amino-4-carboxylic acid). The incorporation of this paramagnetic amino acid results in a nitroxide spin label coupled rigidly to the alpha-carbon, providing direct detection of peptide backbone dynamics by EPR. We applied this approach to phospholamban, which regulates cardiac calcium transport. The synthesis of this amphipathic 52-amino-acid membrane peptide including TOAC is a challenge, especially in the addition of TOAC and the next several amino acids. Therefore, EPR of synthetic intermediates, reconstituted into lipid bilayers, was used to ensure complete coupling and 9-fluorenylmethoxycarbonyl (Fmoc) deprotection. The attachment of Fmoc-TOAC-OH leads to strong immobilization of the spin label, whereas Fmoc deprotection dramatically mobilizes it, producing an EPR spectral peak that is completely resolved from that observed before deprotection. Similarly, coupling of the next amino acid (Ser) restores the spin label to strong immobilization, giving a peak that is completely resolved from that of the preceding step. For several subsequent steps, the effect of coupling and deprotection is similar but less dramatic. Thus, the sensitivity and resolution of EPR provides a quantitative monitor of completion at each of these critical steps in peptide synthesis. Mass spectrometry, circular dichroism, and Edman degradation were used in concert with EPR to verify the chemistry and characterize the secondary structure. In conclusion, the application of conventional analytical methods in combination with EPR offers an improved approach to optimize the accurate synthesis of TOAC spin-labeled membrane peptides.

Sarcolipin and phospholamban as regulators of cardiac sarcoplasmic reticulum Ca2+ ATPase

The cardiac sarcoplasmic reticulum calcium ATPase (SERCA2a) plays a critical role in maintaining the intracellular calcium homeostasis during cardiac contraction and relaxation. It has been well documented over the years that altered expression and activity of SERCA2a can lead to systolic and diastolic dysfunction. The activity of SERCA2a is regulated by two structurally similar proteins, phospholamban (PLB) and sarcolipin (SLN). Although, the relevance of PLB has been extensively studied over the years, the role SLN in cardiac physiology is an emerging field of study. This review focuses on the advances in the understanding of the regulation of SERCA2a by SLN and PLB. In particular, it highlights the similarities and differences between the two proteins and their roles in cardiac patho-physiology.

Effects of Ser16 phosphorylation on the allosteric transitions of phospholamban/Ca(2+)-ATPase complex

Phosphorylation by protein kinase A and dephosphorylation by protein phosphatase 1 modulate the inhibitory activity of phospholamban (PLN), the endogenous regulator of the sarco(endo)plasmic reticulum calcium Ca(2+) ATPase (SERCA). This cyclic mechanism constitutes the driving force for calcium reuptake from the cytoplasm into the myocite lumen, regulating cardiac contractility. PLN undergoes a conformational transition between a relaxed (R) and tense (T) state, an equilibrium perturbed by the addition of SERCA. Here, we show that the single phosphoryl transfer at Ser16 induces a more pronounced conformational switch to the R state in phosphorylated PLN (pPLN). The binding affinity of PLN to SERCA is not affected (K(d) values for the transmembrane domains of pPLN and PLN are approximately 60 microM), supporting the hypothesis that phosphorylation at Ser16 does not dissociate PLN from SERCA. However, the binding surface and dynamics in domain Ib (residues 22-31) change substantially upon phosphorylation. Since PLN can be singly or doubly phosphorylated at Ser16 and Thr17, we propose that these sites remotely control the conformation of domain Ib. These findings constitute a paradigm for how post-translational modifications such as phosphorylation in the cytoplasmic portion of membrane proteins control intramembrane protein-protein interactions.