What the structure of a calcium pump tells us about its mechanism

Linking Biochemical and Structural States of SERCA: Achievements, Challenges, and New Opportunities

Sarcoendoplasmic reticulum calcium ATPase (SERCA), a member of the P-type ATPase family of ion and lipid pumps, is responsible for the active transport of Ca²⁺ from the cytoplasm into the sarcoplasmic reticulum lumen of muscle cells, into the endoplasmic reticulum (ER) of non-muscle cells. X-ray crystallography has proven to be an invaluable tool in understanding the structural changes of SERCA, and more than 70 SERCA crystal structures representing major biochemical states (defined by bound ligand) have been deposited in the Protein Data Bank. Consequently, SERCA is one of the best characterized components of the calcium transport machinery in the cell. Emerging approaches in the field, including spectroscopy and molecular simulation, now help integrate and interpret this rich structural information to understand the conformational transitions of SERCA that occur during activation, inhibition, and regulation. In this review, we provide an overview of the crystal structures of SERCA, focusing on identifying metrics that facilitate structure-based categorization of major steps along the catalytic cycle. We examine the integration of crystallographic data with different biophysical approaches and computational methods to link biochemical and structural states of SERCA that are populated in the cell. Finally, we discuss the challenges and new opportunities in the field, including structural elucidation of functionally important and novel regulatory complexes of SERCA, understanding the structural basis of functional divergence among homologous SERCA regulators, and bridging the gap between basic and translational research directed toward therapeutic modulation of SERCA.

Reduced SERCA Function Preferentially Affects Wnt Signaling by Retaining E-Cadherin in the Endoplasmic Reticulum

Calcium homeostasis in the lumen of the endoplasmic reticulum is required for correct processing and trafficking of transmembrane proteins, and defects in protein trafficking can impinge on cell signaling pathways. We show here that mutations in the endoplasmic reticulum calcium pump SERCA disrupt Wingless signaling by sequestering Armadillo/β-catenin away from the signaling pool. Armadillo remains bound to E-cadherin, which is retained in the endoplasmic reticulum when calcium levels there are reduced. Using hypomorphic and null SERCA alleles in combination with the loss of the plasma membrane calcium channel Orai allowed us to define three distinct thresholds of endoplasmic reticulum calcium. Wingless signaling is sensitive to even a small reduction, while Notch and Hippo signaling are disrupted at intermediate levels, and elimination of SERCA function results in apoptosis. These differential and opposing effects on three oncogenic signaling pathways may complicate the use of SERCA inhibitors as cancer therapeutics.

Suisse and Treisman describe genetic conditions that reduce calcium in the endoplasmic reticulum to three distinct extents. They find that Wnt signaling is more sensitive to changes in calcium levels than the Notch and Hippo pathways, potentially complicating the use of calcium pump inhibitors as cancer therapeutics.

Inhibition of the sarco/endoplasmic reticulum (ER) Ca2+-ATPase by thapsigargin analogs induces cell death via ER Ca2+ depletion and the unfolded protein response

Calcium (Ca²⁺) is a fundamental regulator of cell signaling and function. Thapsigargin (Tg) blocks the sarco/endoplasmic reticulum (ER) Ca²⁺-ATPase (SERCA), disrupts Ca²⁺ homeostasis, and causes cell death. However, the exact mechanisms whereby SERCA inhibition induces cell death are incompletely understood. Here, we report that low (0.1 μm) concentrations of Tg and Tg analogs with various long-chain substitutions at the O-8 position extensively inhibit SERCA1a-mediated Ca²⁺ transport. We also found that, in both prostate and breast cancer cells, exposure to Tg or Tg analogs for 1 day caused extensive drainage of the ER Ca²⁺ stores. This Ca²⁺ depletion was followed by markedly reduced cell proliferation rates and morphological changes that developed over 2-4 days and culminated in cell death. Interestingly, these changes were not accompanied by bulk increases in cytosolic Ca²⁺ levels. Moreover, knockdown of two key store-operated Ca²⁺ entry (SOCE) components, Orai1 and STIM1, did not reduce Tg cytotoxicity, indicating that SOCE and Ca²⁺ entry are not critical for Tg-induced cell death. However, we observed a correlation between the abilities of Tg and Tg analogs to deplete ER Ca²⁺ stores and their detrimental effects on cell viability. Furthermore, caspase activation and cell death were associated with a sustained unfolded protein response. We conclude that ER Ca²⁺ drainage and sustained unfolded protein response activation are key for initiation of apoptosis at low concentrations of Tg and Tg analogs, whereas high cytosolic Ca²⁺ levels and SOCE are not required.

High-salt intake negatively regulates fat deposition in mouse

High-salt (HS) intake contributes to hypertension and cardiopathy, but the effect of HS on fat deposition is controversial. Feed intake, fat mass, the percentage of abdominal fat, heat production, rate of oxygen consumption and the respiratory exchange ratio of mice on a HS diet were significantly decreased (P < 0.01 or 0.05) compared with mice on a normal-salt (NS) diet. An in vitro experiment with differentiating pre-adipocytes showed reduced fat deposition in the presence of high concentrations of NaCl (>0.05 M). Abdominal fat mRNA profiles and protein measurements showed that 5 known genes involved in lipolysis were up-regulated significantly and 9 genes related to lipogenesis were down-regulated in HS mice. Abundant genes and some proteins (ATP2a1, AGT, and ANGPTL4) related to calcium ion metabolism or the renin-angiotensin system (RAS) were differentially expressed between HS and NS mice. Of special interest, CREB1 phosphorylation (S133 and S142), a key factor involved in calcium signaling and other pathways, was up-regulated in HS mice. By IPA analysis, a network mediated by calcium was established providing the molecular mechanisms underlying the negative effect of HS on fat deposition.

Preexisting domain motions underlie protonation-dependent structural transitions of the P-type Ca2+-ATPase

We have performed microsecond molecular dynamics (MD) simulations to determine the mechanism for protonation-dependent structural transitions of the sarco/endoplasmic reticulum Ca2+-ATPase (SERCA), one of the most prominent members of the large P-type ATPase superfamily that transports ions across biological membranes. The release of two H+ from the transport sites activates SERCA by inducing a structural transition between low (E2) and high (E1) Ca2+-affinity states (E2-to-E1 transition), but the structural mechanism by which transport site deprotonation facilitates this transition is unknown. We performed microsecond all-atom MD simulations to determine the effects of transport site protonation on the structural dynamics of the E2 state in solution. We found that the protonated E2 state has structural characteristics that are similar to those observed in crystal structures of E2. Upon deprotonation, a single Na+ ion rapidly (<10 ns) binds to the transmembrane transport sites and induces a kink in M5, disrupts the M3-M5 interface, and increases the mobility of the M3/A-M3 linker. Principal component analysis showed that counter-rotation of the cytosolic N-A domains about the membrane normal axis, which is the primary motion driving the E2-to-E1 transition, is present in both protonated and deprotonated E2 states; however, protonation-dependent structural changes in the transmembrane domain control the hierarchical organization and amplitude of this motion. We propose that preexisting rigid-body domain motions underlie structural transitions of SERCA, where the functionally important directionality is preserved while transport site protonation controls the dominance and amplitude of motion to shift the equilibrium between the E1 and E2 states. We conclude that ligand-induced modulation of preexisting domain motions is likely a common theme in structural transitions of the P-type ATPase superfamily.

The role of skeletal-muscle-based thermogenic mechanisms in vertebrate endothermy

Thermogenesis is one of the most important homeostatic mechanisms that evolved during vertebrate evolution. Despite its importance for the survival of the organism, the mechanistic details behind various thermogenic processes remain incompletely understood. Although heat production from muscle has long been recognized as a thermogenic mechanism, whether muscle can produce heat independently of contraction remains controversial. Studies in birds and mammals suggest that skeletal muscle can be an important site of non-shivering thermogenesis (NST) and can be recruited during cold adaptation, although unequivocal evidence is lacking. Much research on thermogenesis during the last two decades has been focused on brown adipose tissue (BAT). These studies clearly implicate BAT as an important site of NST in mammals, in particular in newborns and rodents. However, BAT is either absent, as in birds and pigs, or is only a minor component, as in adult large mammals including humans, bringing into question the BAT-centric view of thermogenesis. This review focuses on the evolution and emergence of various thermogenic mechanisms in vertebrates from fish to man. A careful analysis of the existing data reveals that muscle was the earliest facultative thermogenic organ to emerge in vertebrates, long before the appearance of BAT in eutherian mammals. Additionally, these studies suggest that muscle-based thermogenesis is the dominant mechanism of heat production in many species including birds, marsupials, and certain mammals where BAT-mediated thermogenesis is absent or limited. We discuss the relevance of our recent findings showing that uncoupling of sarco(endo)plasmic reticulum Ca(2+)-ATPase (SERCA) by sarcolipin (SLN), resulting in futile cycling and increased heat production, could be the basis for NST in skeletal muscle. The overall goal of this review is to highlight the role of skeletal muscle as a thermogenic organ and provide a balanced view of thermogenesis in vertebrates.

Microsecond molecular dynamics simulations of Mg²⁺- and K⁺-bound E1 intermediate states of the calcium pump

We have performed microsecond molecular dynamics (MD) simulations to characterize the structural dynamics of cation-bound E1 intermediate states of the calcium pump (sarcoendoplasmic reticulum Ca²⁺-ATPase, SERCA) in atomic detail, including a lipid bilayer with aqueous solution on both sides. X-ray crystallography with 40 mM Mg²⁺ in the absence of Ca²⁺ has shown that SERCA adopts an E1 structure with transmembrane Ca²⁺-binding sites I and II exposed to the cytosol, stabilized by a single Mg²⁺ bound to a hybrid binding site I'. This Mg²⁺-bound E1 intermediate state, designated E1•Mg²⁺, is proposed to constitute a functional SERCA intermediate that catalyzes the transition from E2 to E1•2Ca²⁺ by facilitating H⁺/Ca²⁺ exchange. To test this hypothesis, we performed two independent MD simulations based on the E1•Mg²⁺ crystal structure, starting in the presence or absence of initially-bound Mg²⁺. Both simulations were performed for 1 µs in a solution containing 100 mM K⁺ and 5 mM Mg²⁺ in the absence of Ca²⁺, mimicking muscle cytosol during relaxation. In the presence of initially-bound Mg²⁺, SERCA site I' maintained Mg²⁺ binding during the entire MD trajectory, and the cytosolic headpiece maintained a semi-open structure. In the absence of initially-bound Mg²⁺, two K⁺ ions rapidly bound to sites I and I' and stayed loosely bound during most of the simulation, while the cytosolic headpiece shifted gradually to a more open structure. Thus MD simulations predict that both E1•Mg²⁺ and E•2K+ intermediate states of SERCA are populated in solution in the absence of Ca²⁺, with the more open 2K+-bound state being more abundant at physiological ion concentrations. We propose that the E1•2K⁺ state acts as a functional intermediate that facilitates the E2 to E1•2Ca²⁺ transition through two mechanisms: by pre-organizing transport sites for Ca²⁺ binding, and by partially opening the cytosolic headpiece prior to Ca²⁺ activation of nucleotide binding.

Hyperactivation of the human plasma membrane Ca2+ pump PMCA h4xb by mutation of Glu99 to Lys

The transport of calcium to the extracellular space carried out by plasma membrane Ca(2+) pumps (PMCAs) is essential for maintaining low Ca(2+) concentrations in the cytosol of eukaryotic cells. The activity of PMCAs is controlled by autoinhibition. Autoinhibition is relieved by the binding of Ca(2+)-calmodulin to the calmodulin-binding autoinhibitory sequence, which in the human PMCA is located in the C-terminal segment and results in a PMCA of high maximal velocity of transport and high affinity for Ca(2+). Autoinhibition involves the intramolecular interaction between the autoinhibitory domain and a not well defined region of the molecule near the catalytic site. Here we show that the fusion of GFP to the C terminus of the h4xb PMCA causes partial loss of autoinhibition by specifically increasing the Vmax. Mutation of residue Glu(99) to Lys in the cytosolic portion of the M1 transmembrane helix at the other end of the molecule brought the Vmax of the h4xb PMCA to near that of the calmodulin-activated enzyme without increasing the apparent affinity for Ca(2+). Altogether, the results suggest that the autoinhibitory interaction of the extreme C-terminal segment of the h4 PMCA is disturbed by changes of negatively charged residues of the N-terminal region. This would be consistent with a recently proposed model of an autoinhibited form of the plant ACA8 pump, although some differences are noted.

Calcium binding and allosteric signaling mechanisms for the sarcoplasmic reticulum Ca²+ ATPase

The sarcoplasmic reticulum Ca²⁺ ATPase (SERCA) is a membrane-bound pump that utilizes ATP to drive calcium ions from the myocyte cytosol against the higher calcium concentration in the sarcoplasmic reticulum. Conformational transitions associated with Ca²⁺-binding are important to its catalytic function. We have identified collective motions that partition SERCA crystallographic structures into multiple catalytically-distinct states using principal component analysis. Using Brownian dynamics simulations, we demonstrate the important contribution of surface-exposed, polar residues in the diffusional encounter of Ca²⁺. Molecular dynamics simulations indicate the role of Glu309 gating in binding Ca²⁺, as well as subsequent changes in the dynamics of SERCA's cytosolic domains. Together these data provide structural and dynamical insights into a multistep process involving Ca²⁺ binding and catalytic transitions.

Tracing cytoplasmic Ca(2+) ion and water access points in the Ca(2+)-ATPase

Sarco(endo)plasmic reticulum Ca(2+)-ATPase (SERCA) transports two Ca(2+) ions across the membrane of the sarco(endo)plasmic reticulum against the concentration gradient, harvesting the required energy by hydrolyzing one ATP molecule during each transport cycle. Although SERCA is one of the best structurally characterized membrane transporters, it is still largely unknown how the transported Ca(2+) ions reach their transmembrane binding sites in SERCA from the cytoplasmic side. Here, we performed extended all-atom molecular dynamics simulations of SERCA. The calculated electrostatic potential of the protein reveals a putative mechanism by which cations may be attracted to and bind to the Ca(2+)-free state of the transporter. Additional molecular dynamics simulations performed on a Ca(2+)-bound state of SERCA reveal a water-filled pathway that may be used by the Ca(2+) ions to reach their buried binding sites from the cytoplasm. Finally, several residues that are involved in attracting and guiding the cations toward the possible entry channel are identified. The results point to a single Ca(2+) entry site close to the kinked part of the first transmembrane helix, in a region loaded with negatively charged residues. From this point, a water pathway outlines a putative Ca(2+) translocation pathway toward the transmembrane ion-binding sites.

Kinetic and mesoscopic non-equilibrium description of the Ca(2+) pump: a comparison

We analyse the operation of the Ca(2+)-ATPase ion pump using a kinetic cycle diagram. Using the methodology of Hill, we obtain the cycle fluxes, entropy production and efficiency of the pump. We compare these results with a mesoscopic non-equilibrium description of the pump and show that the kinetic and mesoscopic pictures are in accordance with each other. This gives further support to the mesoscopic theory, which is less restricted and also can include the heat flux as a variable. We also show how motors can be characterised in terms of unidirectional backward fluxes. We proceed to show how the mesoscopic approach can be used to identify fast and slow steps of the model in terms of activation energies, and how this can be used to simplify the kinetic diagram.

Epigallocatechin-3-gallate has dual, independent effects on the cardiac sarcoplasmic reticulum/endoplasmic reticulum Ca2+ ATPase

We determined the effects of epigallocatechin-3-gallate (EGCG) and epicatechin (EC), on pump turnover and Ca2+ transport by the cardiac form of the sarcoplasmic/endoplasmic reticulum Ca2+-ATPase (SERCA). Fluorescence spectroscopy was used to directly measure SERCA ATPase activity and to measure Ca2+ uptake into cardiac sarcoplasmic reticulum (SR) vesicles and microsomes derived from human embryonic kidney (HEK) cells expressing human cardiac SERCA2a. We found that EGCG reduces the maximum velocity of Ca2+ uptake into cardiac SR vesicles and increases the Ca2+-sensitivity of uptake in a concentration-dependent manner. EC is less potent than EGCG in increasing the Ca2+-sensitivity of uptake and does not affect maximum uptake velocity. The EGCG-dependent reduction in Ca2+ uptake velocity is well correlated with direct inhibition of SERCA. The effect of EGCG on the Ca2+-sensitivity of Ca2+ uptake into cardiac SR vesicles is affected by the phosphorylation status of phospholamban (PLB). When cardiac SERCA2a is expressed in HEK cells without PLB, EGCG reduces the maximum velocity of Ca2+ uptake but does not affect the Ca2+-sensitivity of uptake into microsomes derived from these cells indicating that the effect of EGCG on Ca2+-sensitivity requires the presence of PLB. Our results show that EGCG has dual effects on SERCA function in cardiac SR vesicles: it directly affects SERCA by reducing maximum uptake velocity; it increases the Ca2+-sensitivity of Ca2+ uptake in a manner that appears to depend on the interaction between SERCA and PLB.

In and out of the cation pumps: P-type ATPase structure revisited

Active transport across membranes is a crucial requirement for life. P-type ATPases build up electrochemical gradients at the expense of ATP by forming and splitting a covalent phosphoenzyme intermediate, coupled to conformational changes in the transmembrane section where the ions are translocated. The marked increment during the last three years in the number of crystal structures of P-type ATPases has greatly improved our understanding of the similarities and differences of pumps with different ion specificities, since the structures of the Ca2+-ATPase, the Na+,K+-ATPase and the H+-ATPase can now be compared directly. Mechanisms for ion gating, charge neutralization and backflow prevention are starting to emerge from comparative structural analysis; and in combination with functional studies of mutated pumps this provides a framework for speculating on how the ions are bound and released as well as on how specificity is achieved.

What can be learned about the function of a single protein from its various X-ray structures: the example of the sarcoplasmic calcium pump

Improvements in the handling of membrane proteins for crystallization, combined with better synchrotron sources for X-ray diffraction analysis, are leading to clarification of the structural details of an ever increasing number of membrane transporters and receptors. Here we describe how this development has resulted in the elucidation at atomic resolution of a large number of structures of the sarcoplasmic Ca(2+)-ATPase (SERCA1a) present in skeletal muscle. The structures corresponding to the various intermediary states have been obtained after stabilization with structural analogues of ATP and of metal fluorides as mimicks of inorganic phosphate. From these results it is possible, in accordance with previous biochemical and molecular biology data, to give a detailed structural description of both ATP hydrolysis and Ca(2+) transport through the membrane, to serve as the starting point for a fuller understanding of the pump mechanism and, in future studies, on the regulatory role of this ubiquitous intracellular Ca(2+)-ATPase in cellular Ca(2+) metabolism in normal and pathological conditions.

Roles of transmembrane segment M1 of Na+,K+-ATPase and Ca2-ATPase, the gatekeeper and the pivot

In this review we summarize mutagenesis work on the structure-function relationship of transmembrane segment M1 in the Na+,K+-ATPase and the sarco(endo)plasmic reticulum Ca2+-ATPase. The original hypothesis that charged residues in the N-terminal part of M1 interact with the transported cations can be rejected. On the other hand hydrophobic residues in the middle part of M1 turned out to play crucial roles in Ca2+ interaction/occlusion in Ca2+-ATPase and K+ interaction/occlusion in Na+,K+-ATPase. Leu65 of the Ca2+-ATPase and Leu99 of the Na+,K+-ATPase, located at homologous positions in M1, function as gate-locking residues that restrict the mobility of the side chain of the cation binding/gating residue of transmembrane segment M4, Glu309/Glu329. A pivot formed between a pair of a glycine and a bulky residue in M1 and M3 seems critical to the opening of the extracytoplasmic gate in both the Ca2+-ATPase and the Na+,K+-ATPase.

The widely utilized brominated flame retardant tetrabromobisphenol A (TBBPA) is a potent inhibitor of the SERCA Ca2+ pump

TBBPA (tetrabromobisphenol A) is currently the most widely used type of BFR (brominated flame retardant) employed to reduce the combustibility of a large variety of electronic and other manufactured products. Recent studies have indicated that BFRs, including TBBPA, are bio-accumulating within animal and humans. BFRs including TBBPA have also been shown to be cytotoxic and potentially endocrine-disrupting to a variety of cells in culture. Furthermore, TBBPA has specifically been shown to cause disruption of Ca2+ homoeostasis within cells, which may be the underlying cause of its cytotoxicity. In this study, we have demonstrated that TBBPA is a potent non-isoform-specific inhibitor of the SERCA (sarcoplasmic/endoplasmic reticulum Ca2+-ATPase) (apparent K(i) 0.46-2.3 microM), thus we propose that TBBPA inhibition of SERCA contributes in some degree to Ca2+ signalling disruption. TBBPA binds directly to the SERCA without the need to partition into the phospholipid bilayer. From activity results and Ca2+-induced conformational results, it appears that the major effect of TBBPA is to decrease the SERCA affinity for Ca2+ (increasing the K(d) from approx. 1 microM to 30 microM in the presence of 10 microM TBBPA). Low concentrations of TBBPA can quench the tryptophan fluorescence of the SERCA and this quenching can be reversed by BHQ [2,5-di-(t-butyl)-1,4-hydroquinone] and 4-n-nonylphenol, but not thapsigargin, indicating that TBBPA and BHQ may be binding to similar regions in the SERCA.

Molecular determinants of sarco/endoplasmic reticulum calcium ATPase inhibition by hydroquinone‐based compounds

The ion transport activity of the sarco/endoplasmic reticulum calcium ATPase (SERCA) is specifically and potently inhibited by the small molecule 2,5-di-tert-butylhydroquinone (BHQ). In this study, we investigated the relative importance of the nature and position of BHQ's four substituents for enzyme inhibition by employing a combination of experimental and computational techniques. The inhibitory potencies of 21 commercially available or synthesized BHQ derivatives were determined in ATPase activity assays, and 11 compounds were found to be active. Maximum inhibitory potency was observed in compounds with two para hydroxyl groups, whereas BHQ analogues with only one hydroxyl group were still active, albeit with a reduced potency. The results also demonstrated that two alkyl groups were an absolute requirement for activity, with the most potent compounds having 2,5-substituents with four or five carbon atoms at each position. Using the program GOLD in conjunction with the ChemScore scoring function, the structures of the BHQ analogues were docked into the crystal structure of SERCA mimicking the enzyme's E(2) conformation. Analysis of the docking results indicated that inhibitor binding to SERCA was primarily mediated by a hydrogen bond between a hydroxyl group and Asp-59 and by hydrophobic interactions involving the bulky inhibitor alkyl groups. Attempts to dock BHQ into crystal structures corresponding to the E(1) conformation of the enzyme failed, because the conformational changes accompanying the E(2)/E(1) transition severely restricted the size of the binding site, suggesting that BHQ stabilizes the enzyme in its E(2) form. The potential role of Glu309 in enzyme inhibition is discussed in the context of the computational results. The docking scores correlated reasonably well with the measured inhibitory potencies and allowed the distinction between active and inactive compounds, which is a key requirement for future virtual screening of large compound databases for novel SERCA inhibitors.

Side-chain protonation and mobility in the sarcoplasmic reticulum Ca2+-ATPase: implications for proton countertransport and Ca2+ release

Protonation of acidic residues in the sarcoplasmic reticulum Ca(2+)-ATPase (SERCA 1a) was studied by multiconformation continuum electrostatic calculations in the Ca(2+)-bound state Ca(2)E1, in the Ca(2+)-free state E2(TG) with bound thapsigargin, and in the E2P (ADP-insensitive phosphoenzyme) analog state with MgF(4)(2-) E2(TG+MgF(4)(2-)). Around physiological pH, all acidic Ca(2+) ligands (Glu(309), Glu(771), Asp(800), and Glu(908)) were unprotonated in Ca(2)E1; in E2(TG) and E2(TG+MgF(4)(2-)) Glu(771), Asp(800), and Glu(908) were protonated. Glu(771) and Glu(908) had calculated pK(a) values larger than 14 in E2(TG) and E2(TG+MgF(4)(2-)), whereas Asp(800) titrated with calculated pK(a) values near 7.5. Glu(309) had very different pK(a) values in the Ca(2+)-free states: 8.4 in E2(TG+MgF(4)(2-)) and 4.7 in E2(TG) because of a different local backbone conformation. This indicates that Glu(309) can switch between a high and a low pK(a) mode, depending on the local backbone conformation. Protonated Glu(309) occupied predominantly two main, very differently orientated side-chain conformations in E2(TG+MgF(4)(2-)): one oriented inward toward the other Ca(2+) ligands and one oriented outward toward a protein channel that seems to be in contact with the cytoplasm. Upon deprotonation, Glu(309) adopted completely the outwardly orientated side-chain conformation. The contact of Glu(309) with the cytoplasm in E2(TG+MgF(4)(2-)) makes this residue unlikely to bind lumenal protons. Instead it might serve as a proton shuttle between Ca(2+)-binding site I and the cytoplasm. Glu(771), Asp(800), and Glu(908) are proposed to take part in proton countertransport.

Protonation and hydrogen bonding of Ca2+ site residues in the E2P phosphoenzyme intermediate of sarcoplasmic reticulum Ca2+-ATPase studied by a combination of infrared spectroscopy and electrostatic calculations

Protonation of the Ca(2+) ligands of the SR Ca(2+)-ATPase (SERCA1a) was studied by a combination of rapid scan FTIR spectroscopy and electrostatic calculations. With FTIR spectroscopy, we investigated the pH dependence of C=O bands of the Ca(2+)-free phosphoenzyme (E2P) and obtained direct experimental evidence for the protonation of carboxyl groups upon Ca(2+) release. At least three of the infrared signals from protonated carboxyl groups of E2P are pH dependent with pK(a) values near 8.3: a band at 1758 cm(-1) characteristic of nonhydrogen-bonded carbonyl groups, a shoulder at 1720 cm(-1), and part of a band at 1710 cm(-1), both characteristic of hydrogen-bonded carbonyl groups. The bands are thus assigned to H(+) binding residues, some of which are involved in H(+) countertransport. At pH 9, bands at 1743 and 1710 cm(-1) remain which we do not attribute to Ca(2+)/H(+) exchange. We also obtained evidence for a pH-dependent conformational change in beta-sheet or turn structures of the ATPase. With MCCE on the E2P analog E2(TG+MgF(4)(2-)), we assigned infrared bands to specific residues and analyzed whether or not the carbonyl groups of the acidic Ca(2+) ligands are hydrogen bonded. The carbonyl groups of Glu(771), Asp(800), and Glu(908) were found to be hydrogen bonded and will thus contribute to the lower wave number bands. The carbonyl group of some side-chain conformations of Asp(800) is left without a hydrogen-bonding partner; they will therefore contribute to the higher wave number band.

Cloning, characterization, and expression patterns of three sarco/endoplasmic reticulum Ca2+-ATPase isoforms from pearl oyster (Pinctada fucata)

A large amount of calcium is required for mollusk biomineralization. Sarco/endoplasmic reticulum Ca2+-ATPase (SERCA) is a well-known protein with the function of sustaining the calcium homeostasis. How does it possibly function in the process of pearl oyster biomineralization? Three SERCA isoforms, namely PSERA, PSERB, and PSERC were cloned from the pearl oyster, Pinctada fucata. The cDNAs of the three isoforms were isolated by reverse transcription-polymerase chain reaction (RT-PCR) and rapid amplification of cDNA ends. PSERA consisted of 3568 bp encoding 1007 amino acids, PSERB included 3953 bp encoding 1024 amino acids, and PSERC comprised of 3450 bp encoding 1000 amino acids. The three isoforms showed high homology (65%-87%) with SERCAs from other species. Consistent with the results from other invertebrates, Southern blot analysis revealed that the three isoforms originated from a single gene that was also related to SERCA1, SERCA2, and SERCA3 of vertebrates. The splicing mechanism of the three isoforms was similar to that of isoforms of vertebrate SERCA3. Semiquantitative RT-PCR was carried out to study the expression patterns of the three isoforms. The results showed that PSERB was ubiquitously expressed in all tested tissues and was a potential "housekeeping" SERCA isoform; PSERA was expressed in the adductor muscle and foot and was likely to be a muscle-specific isoform, and PSERC was expressed in the other tissues except the adductor muscle or foot with the highest expression levels in the gill and mantle, indicating that it was a non-muscle-specific isoform and might be involved in calcium homeostasis during pearl oyster biomineralization.

Proton paths in the sarcoplasmic reticulum Ca(2+) -ATPase

The sarcoplasmic reticulum Ca(2+)-ATPase (SERCA1a) pumps Ca(2+) and countertransport protons. Proton pathways in the Ca(2+) bound and Ca(2+)-free states are suggested based on an analysis of crystal structures to which water molecules were added. The pathways are indicated by chains of water molecules that interact favorably with the protein. In the Ca(2+) bound state Ca(2)E1, one of the proposed Ca(2+) entry paths is suggested to operate additionally or alternatively as proton pathway. In analogs of the ADP-insensitive phosphoenzyme E2P and in the Ca(2+)-free state E2, the proton path leads between transmembrane helices M5 to M8 from the lumenal side of the protein to the Ca(2+) binding residues Glu-771, Asp-800 and Glu-908. The proton path is different from suggested Ca(2+) dissociation pathways. We suggest that separate proton and Ca(2+) pathways enable rapid (partial) neutralization of the empty cation binding sites. For this reason, transient protonation of empty cation binding sites and separate pathways for different ions are advantageous for P-type ATPases in general.

Toward a general method to observe the phosphate groups of phosphoenzymes with infrared spectroscopy

A general method to study the phosphate group of phosphoenzymes with infrared difference spectroscopy by helper enzyme-induced isotope exchange was developed. This allows the selective monitoring of the phosphate P-O vibrations in large proteins, which provides detailed information on several band parameters. Here, isotopic exchange was achieved at the oxygen atoms of the catalytically important phosphate group that transiently binds to the sarcoplasmic reticulum Ca(2+)-ATPase (SERCA1a). [gamma-(18)O(3)]ATP phosphorylated the ATPase, which produced phosphoenzyme that was initially isotopically labeled. The helper enzyme adenylate kinase regenerated the substrate ATP from ADP (added or generated upon ATP hydrolysis) with different isotopic composition than used initially. With time this produced the unlabeled phosphoenzyme. The method was tested on the ADP-insensitive phosphoenzyme state of the Ca(2+)-ATPase for which the vibrational frequencies of the phosphate group are known, and it was established that the helper enzyme is effective in mediating the isotope exchange process.

Interactions of phosphate groups of ATP and Aspartyl phosphate with the sarcoplasmic reticulum Ca2+-ATPase: an FTIR study

Phosphate binding to the sarcoplasmic reticulum Ca2+-ATPase was studied by time-resolved Fourier transform infrared spectroscopy with ATP and isotopically labeled ATP ([beta-18O2, betagamma-18O]ATP and [gamma-18O3]ATP). Isotopic substitution identified several bands that can be assigned to phosphate groups of bound ATP: bands at 1260, 1207, 1145, 1110, and 1085 cm(-1) are affected by labeling of the beta-phosphate, bands likely near 1154, and 1098-1089 cm(-1) are affected by gamma-phosphate labeling. The findings indicate that the strength of interactions of beta- and gamma- phosphate with the protein are similar to those in aqueous solution. Two bands, at 1175 and 1113 cm(-1), were identified for the phosphate group of the ADP-sensitive phosphoenzyme Ca2E1P. They indicate terminal and bridging P-O bond strengths that are intermediate between those of ADP-insensitive phosphoenzyme E2P and the model compound acetyl phosphate in water. The bridging bond of Ca2E1P is weaker than for acetyl phosphate, which will facilitate phosphate transfer to ADP, but is stronger than for E2P, which will make the Ca2E1P phosphate less susceptible to attack by water.

Calcium pumps in the central nervous system

Two families of Ca2+ transport ATPases are involved in the maintenance of Ca2+ homeostasis in the nervous system, the plasma membrane Ca2+-ATPase that pumps Ca2+ to the extracellular medium and the intracellular sarco/endoplasmic reticulum Ca2+-ATPase that transports Ca2+ from the cytosol to the endoplasmic reticulum. Both types of calcium pumps show precise regulatory properties and they are localized in specific subcellular regions. In this review, we describe the functional and regulatory properties of both families of calcium pumps, their distribution in nerve cells, and their involvement in neurological disorders. The functional characterization of neuronal calcium pumps is very important in order to understand the biochemical processes involved in the maintenance of intracellular calcium in synaptic terminals.

Use of helper enzymes for ADP removal in infrared spectroscopic experiments: application to Ca2+-ATPase

Adenylate kinase (AdK) and apyrase were employed as helper enzymes to remove ADP in infrared spectroscopic experiments that study the sarcoplasmic reticulum Ca(2+)-ATPase. The infrared absorbance changes of their enzymatic reactions were characterized and used to monitor enzyme activity. AdK transforms ADP to ATP and AMP, whereas apyrase consumes ATP and ADP to generate AMP and inorganic phosphate. The benefits of using them as helper enzymes are severalfold: i), both remove ADP generated after ATP hydrolysis by ATPase, which enables repeat of ATP-release experiments several times with the same sample without interference by ADP; ii), AdK helps maintain the presence of ATP for a longer time by regenerating 50% of the initial ATP; iii), apyrase generates free P(i), which can help stabilize the ADP-insensitive phosphoenzyme (E2P); and iv), apyrase can be used to monitor ADP dissociation from transient enzyme intermediates with relatively high affinity to ADP, as shown here for ADP dissociation from the ADP-sensitive phosphoenzyme intermediate (Ca(2)E1P). The respective infrared spectra indicate that ADP dissociation relaxes the closed conformation immediately after phosphorylation partially back toward the open conformation of Ca(2)E1 but does not trigger the transition to E2P. The helper enzyme approach can be extended to study other nucleotide-dependent proteins.

Calcium transport by sarcoplasmic reticulum Ca-ATPase can be investigated on a solid-supported membrane

Sarcoplasmic reticulum (SR) native vesicles incorporating Ca-ATPase are adsorbed on a solid-supported lipid membrane (SSM). Upon adsorption, the ion pumps are chemically activated by concentration jumps of ATP and the capacitive current transients generated by SR Ca-ATPase are measured under potentiostatic conditions. The Michaelis-Menten constant, K(M), for ATP is evaluated by varying the concentration of ATP in the activating solution. This preliminary result shows that ion transport by SR Ca-ATPase can be suitably investigated by a technique based on concentration jumps on an SSM.

Normal mode-based fitting of atomic structure into electron density maps: application to sarcoplasmic reticulum Ca-ATPase

A method for the flexible docking of high-resolution atomic structures into lower resolution densities derived from electron microscopy is presented. The atomic structure is deformed by an iterative process using combinations of normal modes to obtain the best fit of the electron microscopical density. The quality of the computed structures has been evaluated by several techniques borrowed from crystallography. Two atomic structures of the SERCA1 Ca-ATPase corresponding to different conformations were used as a starting point to fit the electron density corresponding to a different conformation. The fitted models have been compared to published models obtained by rigid domain docking, and their relation to the known crystallographic structures are explored by normal mode analysis. We find that only a few number of modes contribute significantly to the transition. The associated motions involve almost exclusively rotation and translation of the cytoplasmic domains as well as displacement of cytoplasmic loops. We suggest that the movements of the cytoplasmic domains are driven by the conformational change that occurs between nonphosphorylated and phosphorylated intermediate, the latter being mimicked by the presence of vanadate at the phosphorylation site in the electron microscopy structure.

Structural basis of ion pumping by Ca2+-ATPase of the sarcoplasmic reticulum

The structures of the Ca2+-ATPase (SERCA1a) have been determined for five different states by X-ray crystallography. Detailed comparison of the structures in the Ca2+ bound form and unbound (but thapsigargin bound) form reveals that very large rearrangements of the transmembrane helices take place accompanying Ca2+ dissociation and binding and that they are mechanically linked with equally large movements of the cytoplasmic domains. The meanings of the rearrangements of the transmembrane helices and those of the cytoplasmic domains as well as the mechanistic roles of phosphorylation are now becoming clear. Furthermore, the roles of critical amino acid residues identified by extensive mutagenesis studies are becoming evident in terms of atomic structure.

Analysis of functional motions in Brownian molecular machines with an efficient block normal mode approach: myosin-II and Ca2+ -ATPase

The structural flexibilities of two molecular machines, myosin and Ca(2+)-ATPase, have been analyzed with normal mode analysis and discussed in the context of their energy conversion functions. The normal mode analysis with physical intermolecular interactions was made possible by an improved implementation of the block normal mode (BNM) approach. The BNM results clearly illustrated that the large-scale conformational transitions implicated in the functional cycles of the two motor systems can be largely captured with a small number of low-frequency normal modes. Therefore, the results support the idea that structural flexibility is an essential part of the construction principle of molecular motors through evolution. Such a feature is expected to be more prevalent in motor proteins than in simpler systems (e.g., signal transduction proteins) because in the former, large-scale conformational transitions often have to occur before the chemical events (e.g., ATP hydrolysis in myosin and ATP binding/phosphorylation in Ca(2+)-ATPase). This highlights the importance of Brownian motions associated with the protein domains that are involved in the functional transitions; in this sense, Brownian molecular machines is an appropriate description of molecular motors, although the normal mode results do not address the origin of the ratchet effect. The results also suggest that it might be more appropriate to describe functional transitions in some molecular motors as intrinsic elastic motions modulating local structural changes in the active site, which in turn gets stabilized by the subsequent chemical events, in contrast with the conventional idea of local changes somehow getting amplified into larger-scale motions. In the case of myosin, for example, we favor the idea that Brownian motions associated with the flexible converter propagates to the Switch I/II region, where the salt-bridge formation gets stabilized by ATP hydrolysis, in contrast with the textbook notion that ATP hydrolysis drives the converter motion. Another useful aspect of the BNM results is that selected low-frequency normal modes have been identified to form a set of collective coordinates that can be used to characterize the progress of a significant fraction of large-scale conformational transitions. Therefore, the present normal mode analysis has provided a stepping-stone toward more elaborate microscopic simulations for addressing critical issues in free energy conversions in molecular machines, such as the coupling and the causal relationship between collective motions and essential local changes at the catalytic active site where ATP hydrolysis occurs.

TNP-AMP binding to the sarcoplasmic reticulum Ca(2+)-ATPase studied by infrared spectroscopy

Infrared spectroscopy was used to monitor the conformational change of 2',3'-O-(2,4,6-trinitrophenyl)adenosine 5'-monophosphate (TNP-AMP) binding to the sarcoplasmic reticulum Ca(2+)-ATPase. TNP-AMP binding was observed in a competition experiment: TNP-AMP is initially bound to the ATPase but is then replaced by beta,gamma-iminoadenosine 5'-triphosphate (AMPPNP) after AMPPNP release from P(3)-1-(2-nitrophenyl)ethyl AMPPNP (caged AMPPNP). The resulting infrared difference spectra are compared to those of AMPPNP binding to the free ATPase, to obtain a difference spectrum that reflects solely TNP-AMP binding to the Ca(2+)-ATPase. TNP-AMP used as an ATP analog in the crystal structure of the sarcoplasmic reticulum Ca(2+)-ATPase was found to induce a conformational change upon binding to the ATPase. It binds with a binding mode that is different from that of AMPPNP, ATP, and other tri- and diphosphate nucleotides: TNP-AMP binding causes partially opposite and smaller conformational changes compared to ATP or AMPPNP. The conformation of the TNP-AMP ATPase complex is more similar to that of the E1Ca(2) state than to that of the E1ATPCa(2) state. Regarding the use of infrared spectroscopy as a technique for ligand binding studies, our results show that infrared spectroscopy is able to distinguish different binding modes.

Variable clinical expression of an identical mutation in the ATP7A gene for Menkes disease/occipital horn syndrome in three affected males in a single family

Two maternal half-brothers presented with huge cephalic hematoma, fatal in one. Skin morphology disclosed lack of elastic fibres. Their maternal uncle is moderately mentally handicapped and has extensive connective tissue disorders. In all these patients, an identical missense mutation in the ATP7A gene was found and confirmed Menkes' disease.

Transconformations of the SERCA1 Ca-ATPase: a normal mode study

The transport of Ca(2+) by Ca-ATPase across the sarcoplasmic reticulum membrane is accompanied by several transconformations of the protein. Relying on the already established functional importance of low-frequency modes in dynamics of proteins, we report here a normal mode analysis of the Ca(2+)-ATPase based on the crystallographic structures of the E1Ca(2) and E2TG forms. The lowest-frequency modes reveal that the N and A(+Nter) domains undergo the largest amplitude movements. The dynamical domain analysis performed with the DomainFinder program suggests that they behave as rigid bodies, unlike the highly flexible P domain. We highlight two types of movements of the transmembrane helices: i), a concerted movement around an axis perpendicular to the membrane which "twists open" the lumenal side of the protein and ii), an individual translational and rotational mobility which is of lower amplitude for the helices hosting the calcium binding sites. Among all modes calculated for E1Ca, only three are enough to describe the transition to E2TG; the associated movements involve almost exclusively the A and N domains, reflecting the closure of the cytoplasmic headpiece and high displacement of the L7-8 lumenal loop. Subsequently, we discuss the potential contribution of the remaining low-frequency normal modes to the transconformations occurring within the overall calcium transport cycle.

Time-resolved charge translocation by sarcoplasmic reticulum Ca-ATPase measured on a solid supported membrane

Sarcoplasmic reticulum vesicles were adsorbed on an octadecanethiol/phosphatidylcholine mixed bilayer anchored to a gold electrode, and the Ca-ATPase contained in the vesicles was activated by ATP concentration jumps both in the absence and in the presence of K(+) ions and at different pH values. Ca(2+) concentration jumps in the absence of ATP were also carried out. The resulting capacitive current transients were analyzed together with the charge under the transients. The relaxation time constants of the current transients were interpreted on the basis of an equivalent circuit. The current transient after ATP concentration jumps and the charge after Ca(2+) concentration jumps in the absence of ATP exhibit almost the same dependence upon the Ca(2+) concentration, with a half-saturating value of approximately 1.5 microM. The pH dependence of the charge after Ca(2+) translocation demonstrates the occurrence of one H(+) per one Ca(2+) countertransport at pH 7 by direct charge-transfer measurements. The presence of K(+) decreases the magnitude of the current transients without altering their shape; this decrease is explained by K(+) binding to the cytoplasmic side of the pump in the E(1) conformation and being released to the same side during the E(1)-E(2) transition.

Lipid-protein interactions in biological membranes: a structural perspective

Lipid molecules bound to membrane proteins are resolved in some high-resolution structures of membrane proteins. An analysis of these structures provides a framework within which to analyse the nature of lipid-protein interactions within membranes. Membrane proteins are surrounded by a shell or annulus of lipid molecules, equivalent to the solvent layer surrounding a water-soluble protein. The lipid bilayer extends right up to the membrane protein, with a uniform thickness around the protein. The surface of a membrane protein contains many shallow grooves and protrusions to which the fatty acyl chains of the surrounding lipids conform to provide tight packing into the membrane. An individual lipid molecule will remain in the annular shell around a protein for only a short period of time. Binding to the annular shell shows relatively little structural specificity. As well as the annular lipid, there is evidence for other lipid molecules bound between the transmembrane alpha-helices of the protein; these lipids are referred to as non-annular lipids. The average thickness of the hydrophobic domain of a membrane protein is about 29 A, with a few proteins having significantly smaller or greater thicknesses than the average. Hydrophobic mismatch between a membrane protein and the surrounding lipid bilayer generally leads to only small changes in membrane thickness. Possible adaptations in the protein to minimise mismatch include tilting of the helices and rotation of side chains at the ends of the helices. Packing of transmembrane alpha-helices is dependent on the chain length of the surrounding phospholipids. The function of membrane proteins is dependent on the thickness of the surrounding lipid bilayer, sometimes on the presence of specific, usually anionic, phospholipids, and sometimes on the phase of the phospholipid.

Sarco/endoplasmic-reticulum calcium ATPase SERCA1 is maintained in the endoplasmic reticulum by a retrieval signal located between residues 1 and 211

The location of sarco/endoplasmic-reticulum calcium ATPase (SERCA) retention/retrieval motifs in the sequence of the SERCA1 has been investigated by examining the subcellular location in COS-7 cells of enhanced-green-fluorescent-protein-tagged calcium-pump chimaeras. These chimaeras have been constructed from the fast-twitch SERCA1 and the plasma-membrane calcium ATPase PMCA3. The N-terminal, central and C-terminal segments of these calcium pumps were exchanged between SERCA1 and PMCA3. The segments exchanged correspond to residues 1-211, 212-711 and 712-994 of SERCA1, and residues 1-264, 265-788 and 789-1159 of PMCA3 respectively. Only chimaeras containing the N-terminal segment of SERCA1 were located in the endoplasmic reticulum (ER), whereas chimaeras containing the N-terminal segment from PMCA3 were able to escape from the ER and enter the endomembrane pathway en route for the plasma membrane. Co-localization of SERCA1 in COS-7 cells with the ER/Golgi-intermediate compartment marker ERGIC53 indicates that SERCA1 is maintained in the ER by a process of retrieval. These results indicate that the N-terminal region of SERCA1, containing transmembrane helices M1 and M2, contains an ER-retrieval signal.

A coarse-grained normal mode approach for macromolecules: an efficient implementation and application to Ca(2+)-ATPase

A block normal mode (BNM) algorithm, originally proposed by Tama et al., (Proteins Struct. Func. Genet. 41:1-7, 2000) was implemented into the simulation program CHARMM. The BNM approach projects the hessian matrix into local translation/rotation basis vectors and, therefore, dramatically reduces the size of the matrix involved in diagonalization. In the current work, by constructing the atomic hessian elements required in the projection operation on the fly, the memory requirement for the BNM approach has been significantly reduced from that of standard normal mode analysis and previous implementation of BNM. As a result, low frequency modes, which are of interest in large-scale conformational changes of large proteins or protein-nucleic acid complexes, can be readily obtained. Comparison of the BNM results with standard normal mode analysis for a number of small proteins and nucleic acids indicates that many properties dominated by low frequency motions are well reproduced by BNM; these include atomic fluctuations, the displacement covariance matrix, vibrational entropies, and involvement coefficients for conformational transitions. Preliminary application to a fairly large system, Ca(2+)-ATPase (994 residues), is described as an example. The structural flexibility of the cytoplasmic domains (especially domain N), correlated motions among residues on domain interfaces and displacement patterns for the transmembrane helices observed in the BNM results are discussed in relation to the function of Ca(2+)-ATPase. The current implementation of the BNM approach has paved the way for developing efficient sampling algorithms with molecular dynamics or Monte Carlo for studying long-time scale dynamics of macromolecules.

Molecular mechanism of copper transport in Wilson disease

Wilson disease is an autosomal recessive disorder of copper metabolism. The Wilson disease protein is a putative copper-transporting P-type ATPase, ATP7B, whose malfunction results in the toxic accumulation of copper in the liver and brain, causing the hepatic and/or neurological symptoms accompanying this disease. The cytosolic N-terminal domain (approximately 70 kDa) of this ATPase comprises six heavy metal-associated domains, each of which contains the conserved metal-binding motif GMTCXXC. The N-terminal domain (Wilson disease copper-binding domain [WCBD]) has been expressed, purified, and characterized using various techniques. The WCBD binds six atoms of copper in the +1 oxidation state competitively, and with a greater affinity than all other metals. The copper atom is coordinated by two cysteines in a distorted linear geometry. Copper binds the WCBD in a cooperative manner and induces secondary and tertiary conformation changes. Zinc binding to the WCBD has also been characterized by circular dichroism spectroscopy and shown to produce conformational changes that are completely different from those induced by copper. The phosphorylation/nucleotide-binding domain of ATP7B has also been expressed and characterized and shown to be capable of binding ATP but lacking ATPase activity. A peptide corresponding to the sixth transmembrane domain of ATP7B has been constructed and shown to undergo secondary conformational changes upon binding a single atom of copper. Finally, a chimeric protein consisting of the WCBD and truncated ZntA, a zinc-transporting ATPase lacking the N-terminal domain, has been constructed and analyzed for metal ion selectivity. These results suggest that the core determines the metal ion specificity of P-type ATPases, and the N-terminal metal-binding domain may play a regulatory role.

Structural similarities of Na,K-ATPase and SERCA, the Ca(2+)-ATPase of the sarcoplasmic reticulum

The crystal structure of SERCA1a (skeletal-muscle sarcoplasmic-reticulum/endoplasmic-reticulum Ca(2+)-ATPase) has recently been determined at 2.6 A (note 1 A = 0.1 nm) resolution [Toyoshima, Nakasako, Nomura and Ogawa (2000) Nature (London) 405, 647-655]. Other P-type ATPases are thought to share key features of the ATP hydrolysis site and a central core of transmembrane helices. Outside of these most-conserved segments, structural similarities are less certain, and predicted transmembrane topology differs between subclasses. In the present review the homologous regions of several representative P-type ATPases are aligned with the SERCA sequence and mapped on to the SERCA structure for comparison. Homology between SERCA and the Na,K-ATPase is more extensive than with any other ATPase, even PMCA, the Ca(2+)-ATPase of plasma membrane. Structural features of the Na,K-ATPase are projected on to the Ca(2+)-ATPase crystal structure to assess the likelihood that they share the same fold. Homology extends through all ten transmembrane spans, and most insertions and deletions are predicted to be at the surface. The locations of specific residues are examined, such as proteolytic cleavage sites, intramolecular cross-linking sites, and the binding sites of certain other proteins. On the whole, the similarity supports a shared fold, with some particular exceptions.