Dwarf open reading frame (DWORF) is a direct activator of the sarcoplasmic reticulum calcium pump SERCA

  1. M'Lynn E Fisher
  2. Elisa Bovo
  3. Rodrigo Aguayo-Ortiz
  4. Ellen E Cho
  5. Marsha P Pribadi
  6. Michael P Dalton
  7. Nishadh Rathod
  8. M Joanne Lemieux
  9. L Michel Espinoza-Fonseca
  10. Seth L Robia
  11. Aleksey V Zima
  12. Howard S Young  Is a corresponding author
  1. Department of Biochemistry, University of Alberta, Canada
  2. Department of Cell and Molecular Physiology, Stritch School of Medicine, Loyola University Chicago, United States
  3. Center for Arrhythmia Research, Department of Internal Medicine, Division of Cardiovascular Medicine, University of Michigan, United States

Abstract

The sarco-plasmic reticulum calcium pump (SERCA) plays a critical role in the contraction-relaxation cycle of muscle. In cardiac muscle, SERCA is regulated by the inhibitor phospholamban. A new regulator, dwarf open reading frame (DWORF), has been reported to displace phospholamban from SERCA. Here, we show that DWORF is a direct activator of SERCA, increasing its turnover rate in the absence of phospholamban. Measurement of in-cell calcium dynamics supports this observation and demonstrates that DWORF increases SERCA-dependent calcium reuptake. These functional observations reveal opposing effects of DWORF activation and phospholamban inhibition of SERCA. To gain mechanistic insight into SERCA activation, fluorescence resonance energy transfer experiments revealed that DWORF has a higher affinity for SERCA in the presence of calcium. Molecular modeling and molecular dynamics simulations provide a model for DWORF activation of SERCA, where DWORF modulates the membrane bilayer and stabilizes the conformations of SERCA that predominate during elevated cytosolic calcium.

Introduction

The sarco-endoplasmic reticulum calcium pump (SERCA) is an ion transporting ATPase that plays a critical role in intracellular calcium signaling. SERCA maintains calcium content of the sarco-endoplasmic reticulum, creating a 2000-fold gradient that can be mobilized for signaling via inositol-1,4,5-triphosphate receptor (IP3R) and ryanodine receptor (RyR) calcium channels. This signaling mechanism is essential for normal cell physiology, and disordered calcium handling underlies a diverse array of diseases including cardiomyopathies (MacLennan and Kranias, 2003), skeletal muscle disorders (Treves et al., 2017), and neurological diseases (Mekahli et al., 2011). In the heart, the cardiac-specific isoform SERCA2a maintains sarcoplasmic reticulum (SR) calcium stores, which is the major source of calcium for cardiac contractility. While SERCA calcium transport is essential for all cells, it is particularly important in cardiomyocytes where defects in SERCA activity or regulation are associated with cardiomyopathies and heart failure (Schmitt et al., 2003; Haghighi et al., 2006; Haghighi et al., 2003). This connection with heart failure has focused attention on SERCA as a possible target for therapy. An initially promising approach in animal models, increasing SERCA expression via gene delivery (Byrne et al., 2008), has proven challenging in human clinical trials (Hulot et al., 2017). Nonetheless, the rationale for treating heart failure by improving calcium transport function remains compelling and orthogonal approaches may be required.

In the heart, SERCA activity is regulated to allow dynamic calcium homeostasis, which changes in response to the need for cardiac output during rest or exertion. A small transmembrane protein, phospholamban (PLN), is the primary regulatory subunit of SERCA in cardiac muscle. PLN inhibits SERCA by decreasing its apparent calcium affinity, thereby reducing both the rate of calcium reuptake and the amount of calcium in the SR. PLN inhibition of SERCA controls the rate of relaxation and the SR calcium load for subsequent contractions, modulating both the lusitropic and inotropic properties of the heart. PLN inhibition of SERCA is relieved by phosphorylation with the main mechanism involving protein kinase A (PKA) and the β-adrenergic pathway (Kirchberber et al., 1975). Under resting conditions, PLN is an important brake that prevents SR calcium overload and the arrhythmogenic consequences (Desantiago et al., 2008). In turn, regulation of PLN by phosphorylation – regulation of the regulator – creates a dynamic SR calcium load that can respond to the sympathetic need for increased cardiac output, the so-called ‘fight-or-flight’ response. Thus, there exists a pool of SERCA pumps in the SR membrane and their combined calcium transport capacity can be finely tuned – upregulated or downregulated – depending on the requirement for cardiac contractility.

Since the initial discovery of PLN decades ago (Kirchberber et al., 1975), it has stood as the only regulatory subunit of SERCA in ventricular muscle. This changed recently with the discovery of a small transmembrane SERCA-binding protein, dwarf open reading frame (DWORF) (Nelson et al., 2016). DWORF has weak sequence similarity to PLN and its skeletal muscle homolog sarcolipin (SLN), and it appears to comprise a transmembrane peptide of unknown structure and function (Figure 1). DWORF was found to increase SERCA activity by opposing PLN function, leading to the hypothesis that DWORF is a non-inhibitory competitor that displaces PLN from the inhibitory groove of SERCA (Makarewich et al., 2018). This raised the question, why is an additional means of reversing PLN inhibition necessary? PLN inhibition can be reversed by phosphorylation (Karczewski et al., 1997; Catalucci et al., 2009), elevated calcium concentrations, and we can now add DWORF as another level of redundancy in reversing SERCA inhibition by PLN. In the present study, we show that DWORF is a direct activator of SERCA activity, in addition to its role in displacing PLN from the inhibitory groove of SERCA. To elucidate DWORF function, we compared the ability of DWORF and PLN to directly regulate SERCA function. We measured the effects of DWORF on SERCA in a membrane reconstitution system, and we measured SERCA-dependent calcium dynamics and SERCA binding by DWORF during active intracellular calcium signaling in cells. We then modeled the structure of DWORF and the SERCA-DWORF complex, and evaluated their properties using molecular dynamics (MD) simulations. The results provide insight into DWORF structure and function and a role for DWORF in the direct regulation of the SERCA calcium pump.

Sequence alignment, secondary structure prediction, and transmembrane domain prediction for phospholamban and dwarf open reading frame (DWORF).

The predicted transmembrane domains are in bold letters. The sequences are aligned between Asn34 of phospholamban and Pro15 of DWORF (Primeau et al., 2018). The helical regions of phospholamban are based on the NMR structure (PDB code 2KYV). The helical regions of DWORF are based on sequence predictions and molecular dynamics simulations in the present study.

Results

SERCA activity in the presence of DWORF

Previous studies of DWORF function focused on the apparent calcium affinity and not the maximal activity of SERCA in the presence of DWORF (Nelson et al., 2016; Makarewich et al., 2018). These studies demonstrated that DWORF does not alter the apparent calcium affinity of SERCA. A recent study claimed to show that DWORF increases the apparent calcium affinity of SERCA (Li et al., 2021), though our current data (Figure 2) and previous studies by others (Nelson et al., 2016; Makarewich et al., 2018) do not support their conclusion. To investigate if there was a direct functional effect of DWORF on the apparent calcium affinity (KCa) or maximal activity (Vmax) of SERCA, we used an in vitro membrane co-reconstitution system (e.g., Ceholski et al., 2012; Trieber et al., 2009; Trieber et al., 2005). To achieve this, DWORF was co-reconstituted into membrane vesicles with SERCA and the calcium-dependent ATPase activity was measured (Figure 2). The previously reported ATPase activity curves for SERCA in the absence and presence of PLN (Glaves et al., 2019; Gorski et al., 2013) were compared to the activity curves for SERCA in the absence and presence of DWORF. At a molar ratio of ~1–2 DWORF per SERCA, DWORF increased the turnover rate of SERCA at nearly all calcium concentrations tested (0.01–15 µM free calcium). The maximal activity (Vmax) of SERCA in the absence of DWORF increased 1.7-fold in the presence of DWORF (Vmax values were 4.1 ± 0. and 6.9 ± 0.1 µmol/min/mg, respectively). This level of SERCA activation is comparable to the small-molecule activator CDN1163 (Kang et al., 2016). There was no significant effect of DWORF on the apparent calcium affinity of SERCA (KCa values of 0.42 ± 0.03 and 0.48 ± 0.03 µM calcium, respectively). Thus, DWORF acted as a direct activator of SERCA by increasing its maximal activity (turnover rate).

ATPase activity of sarco-endoplasmic reticulum calcium pump (SERCA) co-reconstituted into membrane vesicles in the absence and presence of dwarf open reading frame (DWORF).

(A) ATPase activity of SERCA membrane vesicles in the absence (black line) and presence of DWORF (gray line). (B) Comparison to previously published ATPase activity for SERCA in the absence (black line) and presence of phospholamban (gray line). Only the fitted curves from reference (Anderson et al., 2015) are shown. (C) ATPase activity of SERCA membrane vesicles in the absence (black bar) and presence (gray bar) of the small molecule activator CDN1163 at saturating calcium concentration (3 µM).

Cellular calcium transport in the presence of DWORF

Previous studies have shown that DWORF relieves SERCA inhibition by displacing PLN from the inhibitory groove of SERCA (Nelson et al., 2016; Makarewich et al., 2018). Our data are consistent with this observation, and we add a new facet of DWORF function as a direct activator of SERCA even in the absence of PLN. To determine whether increased ATPase activity was accompanied by increased cellular calcium transport activity, we utilized a newly developed approach for measuring calcium uptake into the endoplasmic reticulum (ER) of live cells (Bovo et al., 2019). The technique used an inducible human SERCA2a stable cell line (t-Rex-293 cells) containing the calcium release channel RyR2, which allows manipulation of ER calcium load, and an ER-targeted calcium indicator R-CEPIA1er, which allows measurement of ER calcium load (Bovo et al., 2019; Bovo et al., 2016). The cell line was transiently transfected with mCer-DWORF (or mCer-PLN) to assess regulation of SERCA2a. mCer-DWORF was predominantly expressed in the ER (estimated from the overlap between the mCer-DWORF and the R-CEPIA1er signals; Figure 3A), with a distribution pattern similar to SERCA (Bovo et al., 2019). To determine SERCA2a function, the plasma membrane was permeabilized with saponin to allow control of the cytosolic-free calcium and ATP concentrations. ER calcium uptake was quantified from the increase in R-CEPIA1er fluorescence after full ER calcium depletion by caffeine followed by RyR2 inhibition with ruthenium red and tetracaine (RyR + Tetr; Figure 3B), which blocks the principal calcium leak pathway (Bovo et al., 2019). The first derivative of ER calcium uptake (d[Ca2+]ER/dt) was plotted against the corresponding [Ca2+]ER to estimate the maximum ER calcium uptake rate and the maximum ER calcium load (Figure 3C). As expected, PLN transfection decreased ER calcium uptake compared to control. The opposite effect was observed with DWORF transfection, which almost doubled SERCA calcium uptake rate over the entire range of physiological ER calcium loads. These data indicated that DWORF acts as a potent activator of SERCA2a and this regulation is direct, not requiring pre-existing inhibition of SERCA by PLN. Moreover, it was noteworthy that DWORF enhanced the maximum ER calcium load (Figure 3C, arrow). Since the thermodynamic driving force for calcium transport is the ATP/ADP ratio, we conclude that DWORF enhances SERCA catalytic efficiency and reduces the energetic cost of calcium transport (Bovo et al., 2019).

Dwarf open reading frame (DWORF) enhances sarco-endoplasmic reticulum calcium pump (SERCA)-dependent calcium dynamics in cells.

(A) Inducible human SERCA2a stable cell line (t-Rex-293 cells + Tetr) was transfected with mCer-DWORF together with the Ca2+ release channel ryanodine receptor (RyR)2 and the endoplasmic reticulum (ER)-targeted Ca2+ indicator R-CEPIA1er. MCer-DWORF expression in these cells showed a similar pattern as R-CEPIA1er. (B) The rate of [Ca2+]ER reuptake after full ER Ca2+ depletion by caffeine followed by RyR2 inhibition with ruthenium red (RR). (C) The ER Ca2+ uptake rate was plotted against the corresponding ER [Ca2+] load. DWORF expression almost doubled ER Ca2+ uptake rate through the entire range of physiological ER Ca2+ loads, which is the opposite trend seen for the SERCA inhibitor phospholamban. DWORF also increased the ER [Ca2+] load (arrow). *p<0.05 versus CTRL.

To investigate the effect of DWORF on SERCA in cellular calcium handling dynamics, we mimicked cardiac calcium handling in the heterologous cell model expressing SERCA and RyR (Figure 4, Figure 4—figure supplement 1). This system generates periodic calcium waves due to spontaneous calcium release by RyR2 followed by calcium reuptake by SERCA (Figure 4A, B). Co-expression of DWORF resulted in increased ER calcium load (Figure 4B; [Ca2+]ER) and increased frequency (Figure 4B; Tau(s)) and amplitude of spontaneous calcium waves, suggesting a significantly faster calcium reuptake due to increased SERCA activity. In contrast, PLN significantly decreased the calcium uptake rate by slowing the calcium wave recovery. Similar effects of DWORF on spontaneous calcium release were observed when calcium waves were measured as cytosolic calcium fluctuations in intact cells (Figure 4C) with a significant increase in the average amplitude and integral of RyR-mediated calcium release events in DWORF-expressing cells (Figure 4D). These results demonstrate that DWORF expression increases ER calcium load by enhancing SERCA-mediated calcium uptake. These data support the observation that DWORF directly increases the calcium-dependent ATPase activity of SERCA (Figure 2).

Figure 4 with 1 supplement see all
Effect of dwarf open reading frame (DWORF) on spontaneous calcium release from the endoplasmic reticulum (ER).

(AB) Co-expression of sarco-endoplasmic reticulum calcium pump (SERCA)2a and ryanodine receptor (RyR2) produced Ca2+ waves due to spontaneous activation of RyR2 followed by SERCA Ca2+ reuptake. DWORF increased the magnitude ([Ca2+]ER) and frequency (rate of recovery of [Ca2+]ER, Tau(s)) of spontaneous Ca2+ waves, while phospholamban significantly decreased it. Similar effects of DWORF on spontaneous calcium release were observed when Ca2+ waves were measured as cytosolic Ca2+ fluctuations in intact cells. (C) Average amplitude and frequency of RyR-mediated Ca2+ release events were significantly increased in DWORF-expressing cells. (D) The integral of RyR-mediated Ca2+ release events was significantly increased in DWORF-expressing cells.

Relative affinity of DWORF binding to SERCA

To gain insight into the mechanism of SERCA activation by DWORF, we measured the relative binding of PLN and DWORF to SERCA in a cell environment. Fluorescent protein tags were fused to SERCA, PLN, and DWORF, and intermolecular fluorescence resonance energy transfer (FRET) was quantified in cells expressing pairs of fusion proteins. This method was recently used to show that PLN and DWORF had similar apparent affinities for SERCA in intact cells (Makarewich et al., 2018; Singh et al., 2019). In addition, we previously showed that PLN affinity for SERCA is sensitive to calcium, where the affinity was decreased by high calcium in permeabilized cells or during extended calcium elevations in rapidly paced cardiomyocytes (Bidwell et al., 2011). This apparent change in affinity suggests that PLN binds more avidly to the conformations of SERCA that prevail under resting conditions when cytosolic calcium is low. The enhanced interaction with calcium-free forms of SERCA offers a mechanistic explanation for the effect of PLN on the apparent calcium affinity of SERCA. PLN binds to SERCA in a groove formed by transmembrane segments M2, M6, and M9 (Akin et al., 2013; Toyoshima et al., 2013; Winther et al., 2013). Upon calcium binding to SERCA, M2 undergoes a conformational change that closes the inhibitory groove and forms the calcium-bound E1 state of SERCA. PLN appears to act as a competitive inhibitor of calcium binding by impeding groove closure and formation of the calcium-bound E1 state.

To determine how the SERCA-DWORF regulatory complex responds to changes in calcium, we permeabilized cells expressing Cer-SERCA and YFP-PLN or YFP-DWORF in solutions mimicking low (diastolic; 0.1 µM) or high (systolic; 3 μM) intracellular calcium concentrations (Figure 5). FRET was quantified for each cell and compared to that cell's level of expression of the YFP acceptor. As previously observed (Kelly et al., 2008), FRET between SERCA-PLN was lowest for cells expressing low levels of protein and increased to a maximal level of ~25% FRET efficiency for high-expressing cells (Figure 5A). The relationship was well-described by a hyperbolic fit that yielded the apparent dissociation constant (KD) of the SERCA-PLN complex (Figure 5C). Permeabilization of cells in high calcium yielded a binding curve that was right-shifted to higher protein concentrations, suggesting a decrease in SERCA-PLN binding affinity compared to low calcium conditions (Bidwell et al., 2011). In contrast, the SERCA-DWORF regulatory complex showed the opposite response to calcium, with a small left-shift of the hyperbolic FRET versus protein concentration curve (Figure 5B). In this case, the binding curve was left-shifted to lower protein concentrations, suggesting an increase in SERCA-DWORF binding affinity at high calcium (Figure 5C). The data indicate that high calcium decreases the affinity of PLN for SERCA, yet it has the opposite effect on DWORF where it increases the affinity for SERCA. These data provide a mechanistic framework for DWORF activation of SERCA observed in vitro (Figure 2) and in cells (Figures 3 and 4).

FRET analysis of sarco-endoplasmic reticulum calcium pump–dwarf open reading frame (SERCA-DWORF) interactions.

The average acceptor sensitization FRET efficiency of cells co-transfected with mCer-SERCA2a and either (A) YFP-PLN or (B) YFP-DWORF. FRET efficiency was measured at high and low calcium concentrations to assess the relative affinity of phospholamban (PLN) and DWORF for the calcium-free and calcium-bound conformations of SERCA. (C) Hyperbolic fits to data provide quantification of the apparent dissociation constant (KD) of the SERCA-PLN or SERCA-DWORF regulatory complexes. Ca decreases the apparent affinity of PLN for SERCA, yet it increases the affinity of DWORF for SERCA.

Structure of DWORF

To gain molecular insight into the activation behavior described above, we generated two molecular models of DWORF, one as a continuous α-helix and another as an N-terminal α-helix (residues 1–13), a flexible linker (residues 14–16; centered around Pro15), and a C-terminal α-helix (residues 17–35). The models were embedded in a lipid bilayer and equilibrated using MD simulations for 4 µs (Figure 6). Both starting models quickly converged on a bent helix or helix-linker-helix structure, though the overall structure was dynamic in the simulations and varied somewhat depending on the starting model. A similar behavior has been observed with SLN (Aguayo-Ortiz et al., 2020). The MD simulations support the helix-linker-helix model, with a bent or unwound region centered on Pro15, as the more stable and probable structure of DWORF (Figure 2B). This model was maintained over 4 µs of MD simulations. To further rationalize this model, we carried out predictions of the secondary structure (Buchan and Jones, 2019; Deléage, 2017; Drozdetskiy et al., 2015; Torrisi et al., 2019) and the transmembrane region (Buchan and Jones, 2019; Tusnády and Simon, 2001; Krogh et al., 2001; Rost and Liu, 2003), as well as hydrophobic moment analysis of the DWORF peptide using only the amino acid sequence (Figure 1). The secondary structure predictions suggested that residues 9–35 are likely to be found as an α-helix, and prediction of the transmembrane domain suggested that residues 10–33 reside in the membrane.

Molecular dynamics (MD) simulations of dwarf open reading frame (DWORF) modeled as a continuous helix (A) and as a helix-linker-helix (B).

Shown are snapshots during the simulations (0, 1, 2, 3, and 4 µs), as well as the RMSD (nm), kink angle, per-residue secondary structure (coil in dark gray and helix in light gray), and tilt angles of the transmembrane helix domain (TMD) and the amphipathic helix (AH). Notice that both simulations support a helix-linker-helix structure of DWORF. The MD simulation in (B) indicates that the helix-linker-helix structure is maintained, though it is dynamic during the time course of the simulation.

To evaluate the N-terminal juxtamembrane helix in the helix-linker-helix model of DWORF and the potential ‘helix breaker’ property of Pro15, we carried out hydrophobic moment analysis of the DWORF amino acid sequence (Figures 1 and 7). Hydrophobic moment analysis revealed an amphipathic helix comprising residues 1–13, which is remarkably like the N-terminus of PLN. However, DWORF has a shorter sequence following the proline residue (Pro21 in PLN and Pro15 in DWORF), which includes only the transmembrane domain. PLN has a polar helical region (residues Gln22 to Asn30) that follows the proline residue and precedes the transmembrane domain (residues 31–51), and this feature is absent in DWORF. The presence of the amphipathic helix allows hydrophobic residues to be located on one side of the helix and hydrophilic residues on the other. In the molecular model, the hydrophobic residues are oriented towards the lipid bilayer, while the hydrophilic residues are exposed to the aqueous solvent. As a result, DWORF modulates the local lipid bilayer (Figure 7B). The amphipathic juxtamembrane region is also observed in PLN (PDB code 2KYV; Verardi et al., 2011) and the amyloid precursor protein (PDB code 2LLM; Nadezhdin et al., 2011). Thus, the formation of this juxtamembrane helix suggests that DWORF may mimic some of the structural features of PLN.

Helix-linker-helix model of dwarf open reading frame (DWORF).

(A) Hydrophobic moment analysis revealed amphipathic helices at the N-terminus of DWORF (left panel) and phospholamban (PLN) (right panel). Notice the similar distribution of polar (red) and apolar (black) residues in the helix-kink structure of DWORF and PLN N-terminal to a proline residue (Pro15 in DWORF and Pro21 in PLN). (B) Time-averaged local membrane thickness in the molecular dynamics simulation of the DWORF helix-linker-helix model. Notice how the bilayer becomes thinner around the short transmembrane helix and amphipathic helix of DWORF. The maximum and minimum widths of the lipid bilayer are indicated by the blue-white-red color-coded bar.

Molecular model of the SERCA-DWORF complex

To consider the implications of the helix-linker-helix model of DWORF, we generated structures of DWORF bound to the inhibitory groove of SERCA using protein-protein docking and MD simulations. DWORF was docked to the SERCA-PLN complex (PDB: 4KYT; Akin et al., 2013) with Pro15 of DWORF aligned with Asn34 of PLN (Figure 8). Three replicate SERCA-DWORF complexes were generated using the most populated structures of DWORF based on root-mean-square deviation (RMSD) clustering analysis of the MD simulations. During the 1 µs MD simulation times, the SERCA-DWORF complexes were maintained and the interaction remained stable (Figure 8A) with the exception of MD replicate MD2, which lost contact with DWORF for approximately 25% of the simulation time. However, the SERCA-DWORF interaction was reestablished. While the main focus of this study was not the precise molecular interactions of the complex, there are some interactions between DWORF and SERCA that appear to stabilize the complex (Figure 8B). Pro15 and Ile16 of DWORF interact with Trp107 on M2 of SERCA. At the base of the transmembrane helix of DWORF, Tyr31 packs against a short orthogonal helix at the base of M9 in the inhibitory groove of SERCA (Pro952 and Ile956). In addition to these anchoring interactions, hydrophobic interactions dominate the interface, with the inhibitory groove of SERCA lined by DWORF hydrophobic residues including Leu17, Ile20, Ile23, Val24, Ile27, Ile28, and Ile30. The MD simulations also reveal that the bilayer thickness is modulated by DWORF in the vicinity of the inhibitory groove and the calcium access funnel of SERCA (Figure 8C).

Molecular model for the interaction of sarco-endoplasmic reticulum calcium pump (SERCA) with dwarf open reading frame (DWORF).

(A) Backbone RMSD, number of contacts between SERCA and DWORF, and tilt angle of DWORF transmembrane helix domain in the three 1 µs replicates of SERCA-DWORF molecular dynamics (MD) simulations. The figure to the right shows the representative structure of the SERCA-DWORF complex from the RMSD clustering analysis of the three simulations. SERCA is colored yellow, with the nucleotide-binding domain in green, the phosphorylation domain in blue, and the actuator domain in red. DWORF is shown in gray. (B) SERCA residues located within 3.0 Å of DWORF (left panel) and DWORF residues located within 3.0 Å of SERCA (right panel) colored by the computed occupancy fraction (OF) of the three independent simulations. Only the residues that presented an OF ≥ 0.6 are shown in the figure. (C) Representation of the averaged local membrane thickness of SERCA-DWORF complex calculated with the first replicate MD1.

The two main assumptions that underlie this model include (i) DWORF binds to the inhibitory groove of SERCA and (ii) Pro15 of DWORF plays a role in disrupting interactions that contribute to SERCA inhibition (Figure 9A,B). First, DWORF binding to the inhibitory groove of SERCA is supported by the existing model that DWORF activates SERCA by displacing PLN (Anderson et al., 2015). Second, we tested the importance of Pro15 in the activation of SERCA by measuring the ER calcium uptake rate and ER calcium load in the presence of a Pro15-Ala mutant of DWORF (Figure 9C). The enhanced cellular calcium dynamics observed in the presence wild-type DWORF were abolished in the presence of the Pro15-Ala missense variant of DWORF.

Molecular model for the interaction of sarco-endoplasmic reticulum calcium pump (SERCA) with dwarf open reading frame (DWORF).

(A) SERCA-PLN, SERCA-DWORF, and the overlay are shown in cartoon format. The molecular model of DWORF (Figure 6A) was superimposed on the structure of the SERCA-PLN complex (PDB ID: 4KYT) according to the sequence alignment in Figure 1. SERCA is colored yellow, with the nucleotide-binding domain in green, the phosphorylation domain in blue, and the actuator domain in red. Phospholamban (PLN) is shown in cyan and DWORF in gray. (B) Close-up view of the SERCA-PLN and SERCA-DWORF complexes. The relative positions of Leu31 and Asn34 are indicated, as well as Pro15 of DWORF, which aligns with Asn34 of PLN. (C) Endoplasmic reticulum calcium uptake rate for the Pro15-Ala (P15A) mutant of DWORF compared to SERCA alone (CTRL) and SERCA in the presence of wild-type DWORF (DWORF). *p<0.01.

Discussion

Functional effect of DWORF

The data presented here demonstrate that DWORF is a direct activator of SERCA, enhancing calcium-dependent ATPase activity of SERCA in reconstituted membranes and SERCA-dependent calcium dynamics in living cells. While DWORF displaces PLN from SERCA as previously reported (Nelson et al., 2016; Makarewich et al., 2018), we have shown that DWORF also has a direct effect on SERCA function that is opposite to the inhibitory properties of PLN. The opposing functions of DWORF and PLN represent a previously unknown regulatory axis for fine-tuning SERCA-dependent calcium homeostasis in response to demand for cardiac output. Previous studies of DWORF overexpression in mouse models with normal PLN expression level demonstrated that peak systolic calcium transient amplitude and SR calcium load were significantly increased with DWORF overexpression (Nelson et al., 2016; Makarewich et al., 2018). While these findings were attributed to DWORF displacement of PLN and relief of SERCA inhibition, these findings are also consistent with a dual mechanism that includes direct activation of SERCA by DWORF.

The membrane reconstitution system allowed us to assess the effect of DWORF on SERCA function in a controlled membrane environment. This SERCA-DWORF two-component system is well-suited for measurements of ATPase activity (Ceholski et al., 2012; Trieber et al., 2009; Trieber et al., 2005; Gorski et al., 2013; Douglas et al., 2005) and calcium transport (Smeazzetto et al., 2017) in the presence of regulatory subunits such as PLN and DWORF. PLN is a known inhibitor of SERCA in that it reduces the apparent affinity for calcium (Figure 2A). Herein, we measured SERCA ATPase activity in the presence of DWORF, which allowed assessment of both the apparent calcium affinity and maximal activity of SERCA. Previous studies of DWORF function focused on the apparent calcium affinity of SERCA (Nelson et al., 2016; Makarewich et al., 2018). At a near equimolar ratio, DWORF acted as a direct activator by increasing the maximal activity of SERCA 1.7-fold without an effect on the apparent calcium affinity (Figure 2B). The magnitude of the Vmax increase seen in the presence of DWORF (~1.7 fold) is similar to that seen in the presence of the small-molecule activator of SERCA, CDN1163 (~1.5-fold; Kang et al., 2016). Thus, the important feature of DWORF and PLN regulation in cardiac muscle is that the regulatory subunits have different functions. The primary effect of PLN is to alter the apparent calcium affinity of SERCA, and the primary effect of DWORF is to alter the maximal activity of SERCA. These distinct but complimentary functions allow for a fine level of control of SERCA-mediated calcium homeostasis.

Cellular effect of DWORF

Membrane co-reconstitution systems can provide mechanistic insight into SERCA regulation; however, it is both an advantage and disadvantage that they lack the complexity of a cellular environment. HEK293 cells have been established as a novel cardiomimetic system for evaluating SERCA-dependent calcium-handling and its regulation by PLN (Bovo et al., 2019). This approach provides the complex environment of a cell, while avoiding the confounding effects of imaging cells during contractions. Co-expression of SERCA2a and RyR2 in this cell system causes periodic calcium transients (waves) similar to those that elicit contractions in cardiomyocytes. The HEK cell model lacks endogenous SERCA regulators, and it can be stably transfected with the exogenous regulators PLN and DWORF. This enabled the observation of ER calcium load, calcium uptake rate, and changes in the magnitude and frequency of calcium waves. These are SERCA-dependent cellular responses to changes in calcium concentration and the presence of SERCA regulatory subunits. This system was used to test the effects of DWORF on SERCA-dependent cellular calcium dynamics. As was observed for PLN, DWORF co-localized with SERCA in ER membranes (Figure 3A). DWORF increased two key parameters of SERCA function, ER calcium content (Figure 3B) and ER calcium uptake rate (Figure 3C). These parameters reflect the kinetics and thermodynamics of SERCA in the absence and presence of DWORF, and the effects were opposite to what is seen for PLN.

An alternative measure of ER calcium load and SERCA activity was provided by the cardiomimetic HEK system, which generates periodic calcium waves due to spontaneous calcium release by RyR2 followed by SERCA calcium reuptake. Calcium release via RyR2 is known to depend on ER calcium load, which is determined by SERCA activity. Thus, the cytosolic calcium concentration during periodic calcium waves is an indirect measure of SERCA function in the absence and presence of DWORF (Figure 4A, B). DWORF increased the amplitude and frequency of RyR-mediated calcium release events (Figure 4C, D). The SERCA-mediated calcium reuptake rate, ER calcium load, and calcium waves were all increased in DWORF-expressing cells, strongly suggesting that DWORF directly activates SERCA and increases its catalytic efficiency. These cellular data (Figures 3 and 4) agree with the activity data from reconstituted membranes containing SERCA and DWORF (Figure 2).

This raised the question, what is the mechanism by which DWORF activates SERCA? To address this, FRET efficiency between SERCA-PLN and SERCA-DWORF was measured at high and low calcium concentrations to assess the relative affinity of PLN and DWORF for the calcium-bound and calcium-free conformations of SERCA (Figure 5A, B). The data provided relative measures of the dissociation constants (KD) of the SERCA-PLN and SERCA-DWORF complexes (Figure 5C). Previous studies have shown that PLN remains associated with SERCA even in the presence of calcium (e.g., Bidwell et al., 2011; Dong and Thomas, 2014; Mueller et al., 2004; Li et al., 2004). Consistent with these observations, elevated calcium did not abolish binding of PLN to SERCA, instead it reduced the apparent affinity of PLN for SERCA (Figure 5A, C). In contrast, calcium increased the apparent affinity of DWORF for SERCA (Figure 5B, C), revealing that DWORF prefers to interact with conformations of SERCA that prevail at high calcium. In the context of the current model of DWORF as a competitive inhibitor of PLN binding to SERCA, the data suggest that DWORF would compete more effectively when calcium is elevated, helping to relieve SERCA inhibition each time that increased transport function is required. In addition, the inverted calcium dependence of the SERCA-DWORF interaction offers an explanation for the direct activation of SERCA by DWORF observed in reconstituted membranes (Figure 2) and cellular calcium dynamics (Figures 3 and 4). While PLN inhibition involves preferential interaction and stabilization of SERCA enzymatic states that prevail at low calcium (e.g., Akin et al., 2013; Asahi et al., 2000; Jones et al., 2002), DWORF preferentially binds and stabilizes SERCA states that are populated during calcium elevations. This DWORF interaction enhances the kinetics of rate-limiting steps in the SERCA transport cycle.

Unique structure of DWORF

A critical step was to evaluate the structure of DWORF in comparison to the well-characterized SERCA regulatory subunits PLN (Akin et al., 2013; Verardi et al., 2011) and SLN (Toyoshima et al., 2013; Winther et al., 2013; Mascioni et al., 2002). Both SLN and PLN form continuous transmembrane helices with well-defined orientations in the membrane. The cytoplasmic domain of PLN is longer than that found in SLN and DWORF, and it lies along the membrane surface in the structure determined by NMR spectroscopy (PDB ID: 2KYV; Verardi et al., 2011). In the X-ray crystal structures of the SERCA-PLN (Akin et al., 2013) and SERCA-SLN (Toyoshima et al., 2013; Winther et al., 2013) complexes, PLN and SLN are found as continuous transmembrane helices (residues ~ 24–48 and ~1–31, respectively), though the cytoplasmic domain of PLN was not resolved. The helical transmembrane domains of PLN and SLN facilitate the structural interactions in the SERCA-bound states, and it is a critical feature of these peptides. In contrast, the DWORF transmembrane domain appears to be discontinuous (residues Leu17 to Ser35), with a break at Pro15 and an N-terminal helix that lies along the membrane surface (residues Met1 to Leu13).

The molecular structure of the SERCA-DWORF complex remains unknown. The current hypothesis is that DWORF binds to the inhibitory groove and displaces PLN (Nelson et al., 2016), and that DWORF and PLN have similar affinities for SERCA (Singh et al., 2019). That said, the proposed helix-linker-helix structure of DWORF (Figure 6s and 7) allows us to speculate about potential mechanisms. The break at Pro15 of DWORF occurs at a critical location for PLN and SLN inhibition of SERCA (Figures 1 and 9). Asn34 of PLN (Asn11 in SLN) is an essential residue for SERCA inhibition, and it is positioned to interact with SERCA by the continuous transmembrane helix of PLN. Thus, if DWORF binds to SERCA and replaces PLN, DWORF lacks Asn34 and N-terminal residues of PLN that contribute to SERCA inhibition (Figure 9). This model offers an explanation for why DWORF itself does not inhibit SERCA – the discontinuous transmembrane helix and the substitution of a proline residue are inconsistent with structural features known to contribute to SERCA inhibition. The model also offers an explanation for how DWORF activates SERCA. We have previously proposed that modulation of the lipid bilayer by PLN is a mechanism for enhancing SERCA maximal activity (Glaves et al., 2019). DWORF binding to the inhibitory groove of SERCA would be expected to modulate the lipid bilayer as suggested by the MD simulations (Figure 8C). Importantly, the modulation of the lipid bilayer observed in the MD simulations is more pronounced on the cytoplasmic side of the membrane where it could impact the dynamics of transmembrane helices (e.g., M1 and M2) that form the calcium access funnel. The displacement of PLN, the absence of key inhibitory interactions (two essential residues of PLN include Leu31 and Asn34Trieber et al., 2009; Kimura et al., 1997), and perturbation of the membrane bilayer provide rationales for the higher maximal activity of SERCA in the presence of DWORF.

Materials and methods

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Strain, strain background (Escherichia coli)BL21(DE3)Sigma-AldrichCMC0016Chemically competent cells
Recombinant DNA reagentpMAL c2x DWORF plasmidThis studyMBP-DWORF fusion; H.S. Young lab
Peptide, recombinant proteinHuman DWORFThis studyNCBI:
NP_001339058.1
Purified human DWORF peptide; H.S. Young lab
Chemical compound, drugCDN1163; 4-(1-methylethoxy)-N-(2-methyl-8-quinolinyl)-benzamideSigma-AldrichSML1682Previously published SERCA activator (PMID:26702054)
Chemical compound, drugC12E8; octaethylene glycol monododecyl etherNikkolBL-8SYDetergent for solubilizing SERCA
Chemical compound, drugEgg PC; L-α-phosphatidylcholine (Egg, Chicken)Avanti840051Natural lipid
Chemical compound, drugEgg PE; L-α-phosphatidylethanolamine (Egg, Chicken)Avanti840021Natural lipid
Chemical compound, drugEgg PA; L-α-phosphatidic acid (Egg, Chicken)Avanti840101Natural lipid
Chemical compound, drugIonomycin, calcium saltSigma-AldrichI3909Facilitates the transport of calcium across the plasma membrane
Chemical compound, drugRuthenium redSigma-AldrichR2751RyR inhibitor
Chemical compound, drugCaffeineSigma-AldrichC0750RyR2 agonist
Recombinant DNA reagentpEGFP_RyR2Donated by Dr. C. George (PMID:15047862)Mammalian expression construct containing GFP fused RyR2
Cell line (HEK 293)T-Rex-293 Cell LineThermo Fisher
Scientific
R71007Mammalian cell line for the stable expression of proteins of interest
Recombinant DNA reagentEYFP-PLNSingh et al., 2019Constitutively fluorescent green/yellow fluorescent protein fusion construct
Recombinant DNA reagentEYFP-DWORFSingh et al., 2019Constitutively fluorescent green/yellow fluorescent protein fusion construct
Recombinant DNA reagentmCerulean-SERCABovo et al., 2019Fluorescent protein fusion construct
Cell line
(HEK 293)
AAV-293 cellsAgilent Technologies240073HEK293 cell line optimized for AAV transfection

Materials

All reagents were of the highest purity available: octaethylene glycol monododecyl ether (C12E8; Nikkol Group, Nikko Chemicals Co., Ltd., Tokyo, Japan); egg yolk phosphatidylcholine (EYPC), phosphatidylethanolamine (EYPE), and phosphatidic acid (EYPA) (Avanti Polar Lipids, Alabaster, AL); all reagents used in the coupled enzyme assay include NADH, ATP, PEP, lactate dehydrogenase, and pyruvate kinase (Sigma-Aldrich, Oakville, ON Canada).

Co-reconstitution of DWORF and SERCA

Request a detailed protocol

Recombinant human DWORF was expressed as a maltose-binding protein (MBP) fusion with a TEV cleavage site for removal of MBP. DWORF was purified by a combination of organic extraction (chloroform-isopropanol-water) and reverse-phase HPLC (Douglas et al., 2005). Purified DWORF was stored as lyophilized thin films (63 µg aliquots). SERCA1a was purified from rabbit skeletal muscle SR. For co-reconstitution, lyophilized DWORF (63 µg) was suspended in a 100 µl mixture of trifluoroethanol-water (5:1) and mixed with lipids (360 µg EYPC and 40 µg EYPA) from stock chloroform solutions. The peptide-lipid mixture was dried to a thin film under nitrogen gas and desiccated under vacuum overnight. The peptide-lipid mixture was hydrated in buffer (20 mM imidazole pH 7.0; 100 mM NaCl; 0.02% NaN3) at 50°C for 15 min, cooled to room temperature, and detergent-solubilized by the addition of C12E8 (0.2% final concentration) with vigorous vortexing. Detergent-solubilized SERCA1a was added (300 µg in a total volume of 300 µl) and the reconstitution was stirred gently at room temperature. Detergent was slowly removed by the addition of SM-2 Bio-Beads (Bio-Rad, Hercules, CA) over a 4 hr time course (final weight ratio of 25 Bio-Beads to one detergent). Following detergent removal, the reconstitution was centrifuged over a sucrose step-gradient (20–50% layers) for 1 hr at 100,000 g. The reconstituted proteoliposomes at the gradient interface were removed, flash-frozen in liquid-nitrogen, and stored at −80°C. The final molar ratio was 120 lipids to 2 DWORF to 1 SERCA. To ensure the proper ratio of ~2 DWORF per SERCA in the reconstituted proteoliposomes, quantitative SDS-PAGE was used with known quantities of SERCA and DWORF proteins as standards (Trieber et al., 2009; Trieber et al., 2005; Young et al., 2001).

ATPase activity assays of SERCA reconstitutions

View detailed protocol

ATPase activity of the co-reconstituted proteoliposomes was measured by a coupled-enzyme assay over a range of calcium concentrations from 0.1 μM to 10 μM (Trieber et al., 2009; Trieber et al., 2005). The assay has been adapted to a 96-well format utilizing Synergy 4 (BioTek Instruments) or SpectraMax M3 (Molecular Devices) microplate readers. Data points were collected at 340 nm wavelength, with a well volume of 155 μl containing 10–20 nM SERCA at 30°C (data points collected every 28–39 s for 1 hr). The reactions were initiated by the addition of proteoliposomes to the assay solution. The Vmax (maximal activity) and KCa (apparent calcium affinity) were determined based on non-linear least-squares fitting of the activity data to the Hill equation (Sigma Plot software, SPSS Inc, Chicago, IL). Errors were calculated as the standard error of the mean for a minimum of four independent reconstitutions.

In-cell calcium uptake and spontaneous calcium release methods

Request a detailed protocol

pCMV R-CEPIA1er was a gift from Dr. Masamitsu Iino (Addgene, USA). The vector encoding the human RyR2 cDNA fused to GFP at the N-terminus was kindly provided by Dr. Christopher George (University of Cardiff, UK). The vector encoding human SERCA2a cDNA was kindly provided by Dr. David Thomas (University of Minnesota, USA). The SERCA2a cDNA was cloned into the mCerluean-M1 modified plasmid (Addgene, USA), yielding SERCA2a fused to a modified Cerulean fluorescent protein (mCer) at the N-terminus. SERCA2a cDNA was also cloned into the inducible expression vector pcDNA5/FRT/TO for the generation of SERCA2a stable cell line. The sequences were all verified by single pass primer extension analysis (ACGT Inc, USA).

Generation of SERCA2a stable cell line

Request a detailed protocol

Stable inducible Flp-In T-Rex-293 cell line expressing SERCA2a was generated using the Flp-In T-REx Core Kit (Invitrogen, USA) as described (Bovo et al., 2019). Flp-In T-REx-293 cells were co-transfected with the pOG44 vector encoding the Flp recombinase and the expression vector pcDNA5/FRT/TO containing the SERCA2a cDNA. 48 hr after transfection, the growth medium was replaced with a selection medium containing 100 μg/ml hygromycin. The hygromycin-resistant cell foci were selected and expanded. Stable cell lines were cultured in high glucose Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 100 units/ml penicillin, 100 mg/ml streptomycin, and 10% fetal bovine serum (FBS) at 5% CO2 and 37°C. Expression of SERCA2a in the stable cell line was verified by western blot analysis 48 hr after induction of recombinant pump expression with 1 μg/ml tetracycline. This strategy results in SERCA2a expressed 10-fold over endogenous SERCA2b (Bovo et al., 2019), so calcium uptake is dominated by the exogenous cardiac isoform. MCer-DWORF expression in these cells showed a similar pattern as CEPIA-1er (Figure 4A), indicating that DWORF is preferentially localized in the ER membrane. Since CICR depends on RyR2 expression level, we analyzed calcium waves only in cells that had similar GFP-RyR2 signal.

Confocal microscopy

Request a detailed protocol

Changes in the luminal ER [Ca2+] ([Ca2+]ER) and cytosolic [Ca2+] ([Ca2+]Cyt) were measured with laser scanning confocal microscopy (Radiance 2000 MP, Bio-Rad, UK) equipped with a ×40 oil-immersion objective lens (N.A. = 1.3). T-Rex-293 cells expressing SERCA2a were transiently co-transfected with plasmids containing the cDNA of RyR2, DWORF, and R-CEPIA1er. Experiments were conducted 48 hr after transfection to obtain the optimal level of recombinant protein expression. The surface membrane was permeabilized with saponin (0.005%). Experiments were conducted after wash out of saponin with a solution of 100 mM K-aspartate, 15 mM KCl, 5 mM KH2PO4, 5 mM MgATP, 0.35 mM EGTA, 0.22 mM CaCl2, 0.75 mM MgCl2, 10 mM HEPES (pH 7.2), and 2% dextran (MW: 40,000). Free [Ca2+] and [Mg2+] of this solution were 200 nM and 1 mM, respectively.

[Ca2+]ER measurements

Request a detailed protocol

[Ca2+]ER was recorded as changes in fluorescence intensity of the genetically encoded ER-targeted Ca2+ sensor R-CEPIA1er (Bovo et al., 2016). R-CEPIA1er was excited with a 514 nm line of the argon laser and signal was collected at >560 nm. Line-scan images were collected at a speed of 10 ms/line. The R-CEPIA1er signal (F) was converted to [Ca2+]ER by the following formula: [Ca2+]SE = Kd × [(F - Fmin)/(Fmax - F)]. Fmax was recorded in 5 mM Ca2+ and 5 μM ionomycin and Fmin was recorded after ER Ca2+ depletion with 5 mM caffeine. The Kd (Ca2+ dissociation constant) was 564 μM (Mekahli et al., 2011). SERCA-mediated Ca2+ uptake was calculated as the first derivative of [Ca2+]ER refilling (d[Ca2+]ER/dt) after RyR2 inhibition with ruthenium red (15 μM) and tetracaine (1 mM). RyR2-independent Ca2+ leak was analyzed as the first derivative of [Ca2+]ER decline (d[Ca2+]ER/dt) after simultaneous inhibition of RyR2 and SERCA. ER Ca2+ uptake and Ca2+ leak rates were plotted as a function of [Ca2+]ER to analyze maximum ER Ca2+ uptake rate and maximum ER Ca2+ load. All 2-D images and line-scan measurements for [Ca2+]ER were analyzed with ImageJ software (NIH, USA).

Cytosolic [Ca2+] ([Ca2+]i) measurements

Request a detailed protocol

[Ca2+]i was measured in intact cells with the high-affinity Ca2+ indicator Fluo-4 (Molecular Probes/Invitrogen, Carlsbad, CA). To load the cytosol with Fluo-4, cells were incubated at room temperature with 10 μM Fluo-4 AM for 15 min in Tyrode solution (140 mM NaCl, 4 mM KCl, 0.5 mM CaCl2, 1 mM MgCl2, 10 mM glucose, 10 mM HEPES, pH 7.4), followed by a 20 min wash. Fluo-4 was excited with the 488 nm line of an argon laser and the emission signal collected at wavelengths above 515 nm. Spontaneous Ca2+ waves were measured at different [Ca2+] (0.5, 1, 2, 5, and 10 mM). In the end of each experiment, caffeine was applied to induce maximal ER Ca2+ release.

Statistics

Data are presented as mean ± standard error of the mean (SEM) of n measurements. Statistical comparisons between groups were performed with the Student's t test for unpaired data sets. Differences were considered statistically significant at p<0.05. Statistical analysis and graphical representation of averaged data were carried out on OriginPro7.5 software (OriginLab, USA).

In-cell FRET methods

Fluorescence resonance energy transfer measurements

Request a detailed protocol

Acceptor sensitization FRET was quantified as previously described (Bidwell et al., 2011). AAV 293 cells were cultured in DMEM cell culture medium supplemented with 10% FBS (ThermoScientific, Waltham, MA). Following culture, cells were transiently transfected using either MBS mammalian transfection kit (Agilent Technologies, Stratagene, La Jolla, CA) or Lipofectamine 3000 transfection kit (Invitrogen) with either EYFP-PLN or DWORF, and mCerulean (mCer)-SERCA constructs in a 1:5 molar plasmid ratio with the fluorescent protein fused via a five amino acid linker to the N-terminus (Makarewich et al., 2018; Singh et al., 2019; Autry et al., 2011; Robia et al., 2007). Cells were then plated 24 hr before imaging in four-well chambered coverglass plates coated with poly-D-lysine and imaged utilizing wide-field fluorescent microscopy. Cells were imaged on a Nikon Eclipse Ti 2 equipped with a Photometrix Prime 95B CMOS camera (Tucson, AZ) and Lumencor Spectra X (Beaverton, OR). Imaging was done in a permeabilization buffer containing 120 mM potassium aspartate, 15 mM potassium chloride, 5 mM magnesium ATP, 0.75 mM magnesium chloride, 2% dextran, 5 mM potassium phosphate, 2 mM EGTA, 20 mM HEPES, and either 0 or 1.7 mM added calcium chloride. Cells were permeabilized in 0.005% saponin buffer for 1 min and then washed twice with buffer before imaging. Data were collected using a 20× 0.75 N.A. objective using 50 ms exposure times for all channels and analyzed using in FIJI using a macro selecting cells with a CFP intensity above 150 AU above background, circularity of 0.4–1 AU, and size of 500–2500 pixels with a rolling background size of 200. FRET efficiency was calculated according to Eapp = IDA – a(IAA) – d(IDD) /(IDA – a(IAA) + (G-d)) IDD (Zal and Gascoigne, 2004), where IDA is the intensity of fluorescence of acceptor emission with donor excitation; IAA is the intensity of acceptor fluorescence with acceptor excitation; and IDD is the fluorescence intensity of donor emission with donor excitation. G represents the ratio of sensitized emission to the corresponding amount of donor recover during acceptor photobleaching and acts as a correction factor, which is constant for a given fluorophore and image conditions. Constants a and d are bleed-through constants calculated from a = IDA/ IAA for a control sample transfected with only YFP-labeled SERCA, and d = IDA/ IDD for a sample transfected with only Cer-labeled SERCA. These values were determined to be a = 0.185, d = 0.405, and G = 2.782. FRET intensity of each cell was then compared to the cell’s YFP intensity, which was used as a measure of protein expression ([PLN]). This leads to estimates of apparent dissociation constant (KD), YFP intensity at ½ FRETmax, and intrinsic FRET of the SERCA-PLN complex, which is equal to the FRETmax. Data were fit with a hyperbola of the function FRET = (FRETmax)([PLN])/(K+ [PLN]).

Cell lines

Request a detailed protocol

The Flip-In T-Rex Cell Line was purchased from Thermo Fisher Scientific. The Flip-In T-Rex Cell Line is derived from 293 human embryonic kidney cells. The 293 parental cell line was obtained from the American Type Culture Collection, and every batch of cells has a Certificate of Analysis. Low passage AAV-293 cells (a type of HEK-293 cell) were obtained from Agilent Technologies. All cells were routinely screened for Mycoplasma using DAPI staining and PCR.

Molecular modeling of DWORF

Request a detailed protocol

The protein structure homology-modeling program MODELLER (Webb and Sali, 2016) was used to generate molecular models of DWORF as a continuous α-helix (residues 1–35) and a helix-linker-helix (residues 1–13 modeled as an α-helix; residues 14–16 as random coil; and residues 17–35 as an α-helix). In evaluating the amino acid sequence of DWORF, we carried out secondary structure prediction using PSIPRED (Buchan and Jones, 2019), MLRC (Deléage, 2017), JPRED v4.0 (Drozdetskiy et al., 2015), and Porter v5.0 (34). We performed transmembrane region predictions using MEMSAT-SVM (Buchan and Jones, 2019), HMMTOP (Tusnády and Simon, 2001), TMHMM v2.0 (Krogh et al., 2001), and PredictProtein (Rost and Liu, 2003). Finally, we carried out hydrophobic moment analysis to identify regions of amphipathic helices using PEPWHEEL (Di Scala et al., 2014) and HMOMENT (Rice et al., 2000). The SERCA-DWORF complex was modeled based on the SERCA-PLN structure (PDB: 4KYT; Akin et al., 2013) using PyMol v1.7 (Schrodinger, LLC, 2015) and the SwissModel server (Waterhouse et al., 2018).

Molecular dynamics simulations

Request a detailed protocol

The DWORF and SERCA-DWORF models were inserted in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer containing a total 180 and 370 lipid molecules, respectively, using the membrane builder module of CHARMM-GUI web server (Jo et al., 2008; Jo et al., 2009). The systems were solvated using a TIP3P water model with a minimum margin of 20 Å between the protein and the edges of the periodic box in the z-axis. Chloride and sodium ions were added to reach a concentration of 150 mM and neutralize the total charge of the system. MD simulations were carried out using the Amber ff19SB (Tian et al., 2020) and Lipid 17 force field topologies with the parameters implemented in Amber 18 and AmberTools package (Salomon-Ferrer et al., 2013). The systems were energy minimized and equilibrated with NVT and NPT ensembles following the six-step preparation protocol recommended by CHARMM-GUI (Lee et al., 2016). Langevin thermostat algorithm was used to maintain the temperature at 310 K and the Monte Carlos barostat to set the pressure at 1.0 bar. All bonds involving hydrogens were constrained using the SHAKE algorithm. Each DWORF model was subjected to 4 µs MD simulation, while three independent replicates of the SERCA-DWORF complex were subjected to 1 µs MD simulations. For the MD analysis, MDAnalysis python library (Michaud-Agrawal et al., 2011) was used to convert AMBER MD trajectories and coordinates to GROMACS files. We calculated the backbone RMSD, the tilt angle relative to the membrane normal of the TMD and AH (helanal library), the secondary structure (do_dssp), helix kink angle, number of contacts (mindist) and selection of residues within 3.0 Å of a query structure (select) for MD trajectories using the built-in tools of GROMACS 5.1.4 (Abraham et al., 2015) and MDAnalysis. Local membrane property analysis was performed for both systems with the g_lomepro package (Gapsys et al., 2013) implemented in GROMACS. The frames of each simulation were subjected to a backbone RMSD clustering analysis using a cutoff value of 2.0 and 2.5 Å to retrieve the most representative structures of DWORF and SERCA-DWORF, respectively. The plots were constructed with Gnuplot v5.2 (http://www.gnuplot.info) and figures were generated with PyMOL v1.7.

Data availability

All data generated or analyzed during this study are included in the manuscript.

References

    1. Kirchberber MA
    2. Tada M
    3. Katz AM
    (1975)
    Phospholamban: a regulatory protein of the cardiac sarcoplasmic reticulum
    Recent Advances in Studies on Cardiac Structure and Metabolism 5:103–115.
  1. Software
    1. Schrodinger, LLC
    (2015) The PyMol Molecular Graphics System, version 1.7
    The PyMol Molecular Graphics System.

Decision letter

  1. Mark T Nelson
    Reviewing Editor; University of Vermont, United States
  2. Richard W Aldrich
    Senior Editor; The University of Texas at Austin, United States
  3. Amreen Mughal
    Reviewer; University of Vermont, United States

In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.

Acceptance summary:

The authors found that the peptide, dwarf open reading frame (DWORF), is a direct activator of the sarcoplasmic reticulum (SR) calcium pump (SERCA), which is responsible for maintaining high calcium within the SR. DWORF is the first known peptide activator of SERCA. Furthermore, they show that DWORF has a high affinity of SERCA in the presence of calcium.

Decision letter after peer review:

[Editors’ note: the authors submitted for reconsideration following the decision after peer review. What follows is the decision letter after the first round of review.]

Thank you for submitting your work entitled "Dwarf open reading frame (DWORF) is a direct activator of the sarcoplasmic reticulum calcium pump SERCA" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and a Senior Editor. The following individual involved in review of your submission has agreed to reveal their identity: Amreen Mughal (Reviewer #1).

We are sorry to say that, after extensive consultation with the reviewers, we have decided that your work will not be considered further for publication by eLife.

It was felt that the conclusions were not supported by the data. Also previous studies including by one of the current co-authors provided clear evidence that DWORF competes with phospholamban, contrary to the current study. The general interest is also unclear for the readership of eLife without some understanding of physiological regulation of DWORF's interaction with SERCA.

Reviewer #1:

In this manuscript, Fished and Bravo et al. show that DWORF can act as a direct activator of SERCA. Building on seminal earlier work from Nelson et al. (2016), they proposed that DWORF increases SERCA dependent calcium reuptake. The findings are interesting, but in the absence of in vivo experiments, have limited physiological relevance.

The authors raised an important question, namely why an additional DWORF dependent mechanism is necessary for reversing PLN inhibition, but they have not yet answered this question.

Comments for the authors:

Introduction:

1. Authors did not describe Sarcolipin in the introduction. As they performed experiments with Sarcolipin as well, it is important that authors introduce it and provide information to the readers: why those experiments are important for this manuscript.

2. Authors raised an important question "This raised the question, why is an additional means of reversing PLN inhibition necessary?" But they did not answer this question in the manuscript. It will be worth if authors include some relevant discussion based on their results to show why having DWORF dependent SERCA activation mechanism is important.

Results:

3. Figure 4A: There is a baseline shift in Ca2+[ER] in DWORF group compared to control and PLN. This shift could cause an increase in ER Ca2+ concentration, but absolute change looks similar to other groups. Authors should clarify these results.

4. Figure 4: Please include representative cell images showing changes in Ca2+ concentration in different conditions.

5. All the experiments are performed in the absence of Phospholamban. However, phospholamban will be present in the physiological conditions, hence it is challenging to understand the physiological relevance of this work. There are no in vivo experiments performed to supports these results. Hence, authors should include this as a limitation.

Reviewer #2:

This work report the new finding that Dwarf open reading frame (DWORF) directly stimulates SERCA activity in vitro and in situ and it provides a molecular model for this action derived from FRET experiments and Molecular Dynamics simulations. The authors propose an interesting model for DWORF stimulation of SERCA, which awaits structural determination of the DWORF/SERCA2a complex: that the perturbation of the cytosolic side of the lipid bilayer may impact M1 and M2 movements near the Ca2+ access funnel and facilitate pump loading and concomitantly stimulate transport efficiency.

Strengths:

The combination of multiple techniques to demonstrate the direct interaction between DWORF and SERCA, as well as the clear stimulation of SERCA enzymatic activity (ATPase) and transport activity (Ca2+ uptake) in the presence and more importantly absence of phosopholamban convincingly validates the hypothesis that DWORF directly activates SERCA. The co-reconstitution of SERCA with DWORF clearly shows enzymatic increases to levels comparable to maxes reported in the presence of high phospholamban and the known SERCA agonist CDN1163. The functional measurements of Ca2+ transport and ER loading in the stable SERCA2a HEK-293 cell line is an excellent complement to the enzymology and provides additional confidence for direct DWORF activation. Clearly, this observation expands our knowledge of how SERCA activity is cellularly regulated and is likely an important component to our understanding of excitation-contraction coupling.

Weaknesses:

There are no apparent weaknesses with experimental results or the stated interpretation of these data. However, the discussion of the proposed model of DWORF stimulation generated by the modeling leaves a few unaddressed questions, the most significant being the following. The molecular modeling suggests that DWORF, like PLN binds in the inhibitory groove formed by M2, M6, and M9. Indeed, the proposed mechanism of PLN inhibition of SERCA is by competing for the enzyme conformation which binds Ca2+; specifically, the binding of Ca2+ causes a shift in M2 that closes the inhibitory groove such that the binding of Ca2+ and PLN are mutually exclusive. It is then puzzling how Ca2+-binding facilitates DWORF binding, if indeed it is within the same inhibitory groove.

Comments for the authors:

1) The only grammatical errors identified is in the Methods on page 12, line 6. "Data was fit with…." Should be "Data were fit with…."

2) The authors talk about structural differences between PLN and DWORF and the reader can make their own inferences as to how these minor differences might explain how DWORF binding is stabilized by Ca2+ binding, where PLN is destabilized. Nevertheless, this is something the authors may consider discussing a bit further as they certainly must have pondered the dichotomy. Along these lines, a few of the authors published a recent paper suggesting that SERCA dimerization plays a role in pump stimulation (Biophysical J. 119(7), 2020). It may be worth discussing how DWORF may contribute (or not) to the quaternary structure of SERCA and its physiologic function.

3) A few years ago a couple of the current authors published similar work on DWORF stimulation of SERCA (Makarewich et al., eLife 2018;7:e38319). In that report, they were pretty certain the DWORF mode of action was via competition with phospholamban. Indeed, the reviewers asked them to address the possibility that DWORF might directly modulate SERCA activity and they performed further experiments and concluded that DWORF only worked via PLN competition. Here is Makarewich et al. rebuttal to the reviewer's inquiry, "We have performed suggested experiments 'a' and 'c' and have included these results as Figure 2—figure supplement 3. Consistent with our previous publication on DWORF (Nelson et al., 2016), we found that co-expression of DWORF with unregulated SERCA did not lead to a direct or additive stimulation of SERCA activity (Figure 2—figure supplement 3B and C). Instead, we found that the ability of DWORF to simulate SERCA activity lies in its capacity to compete PLN off of SERCA and reduce the inhibitory effects of PLN on SERCA. This competitive displacement of PLN from SERCA can be very clearly seen in immunoprecipitation experiments shown in Figure 2—figure supplement 3A." Thus, this apparent mechanistic reversal deserves some explanation beyond "Compared to previous reports (12),…" stated on page 3. Specifically, in the 2018 eLife paper, Figure 2, supplement Figure 3B and 3C (purple), the direct stimulation of SERCA by DWORF is not apparent. What was missing in those previous experiments that prevented the observation that seems clear in the current manuscript? This information may be useful to other investigators studying SERCA regulation.

Reviewer #3:

In this study, Fisher et al. investigated the interaction of the sarco-endoplasmic reticulum calcium pump (SERCA) and the SERCA-binding protein small dwarf open reading frame (DWORF). Previous studies have shown that DWORF does not directly activate the SERCA pump and does not change the apparent affinity of SERCA for Ca+2. The current consensus in the field is that DWORF regulates the relative affinity of SERCA for Ca2+ by relieving the inhibition of PLN. In contrast to previous studies, the authors claim that DWORF directly regulates SERCA activity and calcium reuptake that is independent of its known role in relieving the inhibition of PLN. However, this major claim is not supported by the data presented in the manuscript.

Comments for the authors:

1. In Figure 2b it is shown that DWORF directly increases SERCA activity. Is this effect dose depended? Does the molar ratio of DWORF/SERCA affects the Vmax? It is also important to show the regulation of SERCA activity by DWORF in the presence and absence of PLN, phosphorylated PLN and SLN.

2. Related to the previous comment, there is a leftward shift of the Ca2+ affinity curve in the presence of DWORF (Figure 2b), suggesting that the affinity for Ca2+ is increased. Is this significant? How do the authors explain this effect?

3. In Figure 2a the data shows that PLN increases the Vmax of SERCA. This is puzzling as PLN is known to regulate the affinity os SERCA for Ca2+ but does not affect the Vmax. Is there an explanation for this observation? Also, the dose response curve and error bars should be shown on the graph.

4. The dose response curve of the small molecule CDN1163 should be shown in Figure 2C.

5. The authors used the skeletal muscle SERCA isoform SERCA1b in the reconstituted vesicles system and cardiac isoform, Serca2a, in the cell-based experiments. Are there any particular reasons why different isoforms were used in different experiments? Does DWORF regulate the activity of all SERCA isoforms?

6. Figure 3. How do the authors explain the discrepancies between the findings of the current study and the previous studies? For example, using a similar cellular model of DWORF overexpression (Cos7 cells co-transfected with SERCA2a) it has been reported that cells expressing SERCA2a and DWORF alone do not exhibit enhanced SERCA activity indicating that DWORF itself does not biochemically activate SERCA. Additional experiments are needed to clary this discrepancy. What is the effect of DWORF on SERCA2a in the presence of PLN and/or SLN in the cellular model used here? Is the Ca2+ uptake kinetics affected by the stoichiometry of the overexpressed proteins? As the levels of the overexpressed proteins are not titrated, it is likely that these observations represent artifacts associated with protein overexpression.

7. Figure 4a. The calcium transient frequency is significantly increased in the presence of DWORF. This suggest an increased RYR activity in the presence of DOWRF. Is there an indirect/direct effect of the RyR2 overexpression on the observed SERCA Ca2+ kinetics? Is RyR2 necessary for the regulation of SERCA by DWORF?

8. Figure 4b. Ideally the Ca2+ reuptake kinetics, eg Tau, should be measured under electrically stimulated conditions as the frequency of the Ca2+ release will affected the kinetics. It is likely that the increase of Tau by DWORF reflect the increase of the oscillation frequency.

https://doi.org/10.7554/eLife.65545.sa1

Author response

[Editors’ note: The authors appealed the original decision. What follows is the authors’ response to the first round of review.]

Reviewer #1:

[…] Comments for the authors:

Introduction:

1. Authors did not describe Sarcolipin in the introduction. As they performed experiments with Sarcolipin as well, it is important that authors introduce it and provide information to the readers: why those experiments are important for this manuscript.

In our manuscript, we presented previously published data for SERCA-SLN and SEERCA-PLN for comparison to SERCA-DWORF. However, after further consideration we decided the comparison to SERCA-SLN is not physiologically relevant and it has been removed (Figures 1 and 2 have been modified).

2. Authors raised an important question "This raised the question, why is an additional means of reversing PLN inhibition necessary?" But they did not answer this question in the manuscript. It will be worth if authors include some relevant discussion based on their results to show why having DWORF dependent SERCA activation mechanism is important.

The answer to this question is that DWORF is more than just an additional means of reversing PLN inhibition. DWORF is a direct activator of SERCA, which increases the dynamic physiological range of calcium that can be mobilized in response to demand for cardiac output. This question appears in the Introduction (on page 3 of our manuscript) and we have modified this paragraph to provide an immediate answer and preview of our results. The statement now reads: “This raised the question, why is an additional means of reversing PLN inhibition necessary? PLN inhibition can be reversed by phosphorylation (13, 14), elevated calcium concentrations, and we can now add DWORF as another level of redundancy in reversing SERCA inhibition by PLN. In the present study, we show that DWORF is a direct activator of SERCA activity, in addition to its role in displacing PLN from the inhibitory groove of SERCA.”

Results:

3. Figure 4A: There is a baseline shift in Ca2+[ER] in DWORF group compared to control and PLN. This shift could cause an increase in ER Ca2+ concentration, but absolute change looks similar to other groups. Authors should clarify these results.

We will be happy to clarify these results. The baseline shift is caused by an increase in ER calcium concentration, which is due to increased SERCA activity with DWORF expression. This can be seen in the upper panel of Figure 4B. DWORF causes both a baseline shift (increase in ER calcium load) and an absolute change (increase in the amplitude of calcium release from the ER). Changes in these two parameters are directly related to changes in SERCA function and they underscore the functional significance of the novel mechanisms described in our manuscript. We have clarified the description of these results to indicate that ER calcium load increases with DWORF co-expression (page 4). The text now reads:

“Coexpression of DWORF resulted in increased ER calcium load (Figure 4B; [Ca2+]ER) and increased frequency (Figure 4B; Tau(s)) and amplitude of spontaneous calcium waves, suggesting a significantly faster calcium reuptake due to increased SERCA activity.”

4. Figure 4: Please include representative cell images showing changes in Ca2+ concentration in different conditions.

We have added representative cell image line scans showing changes in calcium concentration as Supplemental Figure 1. Further experimental details for these methods, as well as example cell images and line scans, have been previously published and can be found in Bovo et al., Am J Physiol Heart Circ Physiol (2019) 316:H1323-H1331.

5. All the experiments are performed in the absence of Phospholamban. However, phospholamban will be present in the physiological conditions, hence it is challenging to understand the physiological relevance of this work. There are no in vivo experiments performed to supports these results. Hence, authors should include this as a limitation.

The in vivo experiments have been done by others for DWORF in the presence of phospholamban (1,2). The original paper describing DWORF looked at WT, DWORF overexpression, and DWORF knockout mice (2). All mouse models had normal phospholamban expression levels. Peak systolic calcium transient amplitude and SR calcium load were significantly increased in the DWORF overexpression adult cardiomyocytes. These data are completely consistent with the novel observations presented in our manuscript.

As requested by the reviewer, we have elaborated on these points in the Discussion section of our revised manuscript (page 7). The following statements have been added:

“Previous studies of DWORF overexpression in mouse models with normal PLN expression level demonstrated that peak systolic calcium transient amplitude and SR calcium load were significantly increased with DWORF overexpression (11, 12). While these findings were attributed to DWORF displacement of PLN and relief of SERCA inhibition, these findings are also consistent with a dual mechanism that includes direct activation of SERCA by DWORF.”

Reviewer #2:

[…] There are no apparent weaknesses with experimental results or the stated interpretation of these data. However, the discussion of the proposed model of DWORF stimulation generated by the modeling leaves a few unaddressed questions, the most significant being the following. The molecular modeling suggests that DWORF, like PLN binds in the inhibitory groove formed by M2, M6, and M9. Indeed, the proposed mechanism of PLN inhibition of SERCA is by competing for the enzyme conformation which binds Ca2+; specifically, the binding of Ca2+ causes a shift in M2 that closes the inhibitory groove such that the binding of Ca2+ and PLN are mutually exclusive. It is then puzzling how Ca2+-binding facilitates DWORF binding, if indeed it is within the same inhibitory groove.

The reviewer’s assessment that “the binding of Ca2+ and PLN are mutually exclusive” is not supported by current evidence. Early models for the SERCA-PLN interaction suggested that calcium relieves SERCA inhibition by dissociating the complex. However, there is now extensive evidence indicating that PLN remains bound to SERCA in the presence of calcium (e.g. (3-6)). In addition, the crystal structure of the SERCA-PLN complex clearly shows that SERCA is in an E1-like state similar to the calcium-bound state of SERCA (7-9). Importantly, the FRET data shown in Figure 5 of our manuscript clearly shows that PLN remains associated with SERCA in the presence of calcium, though the binding affinity is lower.

We have clarified this point in the Discussion section of our manuscript (page 8). The statement now reads:

“Previous studies have shown that PLN remains associated with SERCA even in the presence of calcium (e.g. (24, 42-44)). […] In contrast, calcium increased the apparent affinity of DWORF for SERCA (Figure 5B and C), revealing that DWORF prefers to interact with conformations of SERCA that prevail at high calcium.”

Comments for the authors:

1) The only grammatical errors identified is in the Methods on page 12, line 6. "Data was fit with…." Should be "Data were fit with…."

This has been corrected.

2) The authors talk about structural differences between PLN and DWORF and the reader can make their own inferences as to how these minor differences might explain how DWORF binding is stabilized by Ca2+ binding, where PLN is destabilized. Nevertheless, this is something the authors may consider discussing a bit further as they certainly must have pondered the dichotomy. Along these lines, a few of the authors published a recent paper suggesting that SERCA dimerization plays a role in pump stimulation (Biophysical J. 119(7), 2020). It may be worth discussing how DWORF may contribute (or not) to the quaternary structure of SERCA and its physiologic function.

The reviewer raises a very interesting question. We do not have direct experimental evidence of a relationship between DWORF/PLN regulation and SERCA dimerization; as far as we know these are distinct structure/function mechanisms. However, both dimerization and DWORF/PLN modulate specific steps in the enzymatic cycle so there may be an undiscovered mechanistic connection. While this is beyond the scope of the present manuscript, a future detailed kinetic study may correlate SERCA dimer conformational coupling and SERCA regulation by DWORF/PLN.

3) A few years ago a couple of the current authors published similar work on DWORF stimulation of SERCA (Makarewich et al., eLife 2018;7:e38319). In that report, they were pretty certain the DWORF mode of action was via competition with phospholamban. Indeed, the reviewers asked them to address the possibility that DWORF might directly modulate SERCA activity and they performed further experiments and concluded that DWORF only worked via PLN competition. Here is Makarewich et al. rebuttal to the reviewer's inquiry, "We have performed suggested experiments 'a' and 'c' and have included these results as Figure 2—figure supplement 3. Consistent with our previous publication on DWORF (Nelson et al., 2016), we found that co-expression of DWORF with unregulated SERCA did not lead to a direct or additive stimulation of SERCA activity (Figure 2—figure supplement 3B and C). Instead, we found that the ability of DWORF to simulate SERCA activity lies in its capacity to compete PLN off of SERCA and reduce the inhibitory effects of PLN on SERCA. This competitive displacement of PLN from SERCA can be very clearly seen in immunoprecipitation experiments shown in Figure 2—figure supplement 3A." Thus, this apparent mechanistic reversal deserves some explanation beyond "Compared to previous reports (12),…" stated on page 3. Specifically, in the 2018 eLife paper, Figure 2, supplement Figure 3B and 3C (purple), the direct stimulation of SERCA by DWORF is not apparent. What was missing in those previous experiments that prevented the observation that seems clear in the current manuscript? This information may be useful to other investigators studying SERCA regulation.

We thank the reviewer for raising this concern, as the answer to this question is quite simple. The previous studies of DWORF, including the original paper (2) and the mentioned eLife paper (1), focused on the apparent affinity of SERCA for calcium (KCa), which is not significantly changed by DWORF. They did not examine the maximal activity (Vmax) of SERCA in these studies, which is the parameter modulated by DWORF. It is difficult to measure Vmax from mouse heart homogenates, but this parameter is accessible in the well-defined biochemical preparations used in the present study. In summary, we do not consider our present results to be in conflict with the previous papers (1,2), as we did not observe an effect of DWORF on the KCa of SERCA. Our results provide a new facet of SERCA regulation by DWORF and build on the important contributions of previous studies (1,2).

To clarify this point, we have added statements to the Results section of our manuscript (page 3). The statements read:

“Previous studies of DWORF function focused on the apparent calcium affinity and not the maximal activity of SERCA in the presence of DWORF (11, 12).” And “Previous studies have shown that DWORF relieves SERCA inhibition by displacing PLN from the inhibitory groove of SERCA (11, 12). Our data are consistent with this observation, and we add a new facet of DWORF function as a direct activator of SERCA even in the absence of PLN.”

Note that during the review of our manuscript by eLife, a publication appeared in the Journal of Biological Chemistry (Li, Yuen, Stroik, Kleinboehl, Cornea, and Thomas). The authors claim to show that DWORF activates SERCA by altering the calcium affinity of SERCA. These data are in conflict with previous studies (1,2), as well as the data in our manuscript. We have addressed this inconsistency in the revised version of our manuscript. We have added the following statement to the Results section of our manuscript (page 3):

“A recent study claimed to show that DWORF increases the apparent calcium affinity of SERCA (15), though our current data (Figure 2) and previous studies by others (11, 12) do not support their conclusion.”

The publication of the Li et al. study further increases the urgency for dissemination of our data and manuscript, which clearly show that DWORF does not alter the KCa of SERCA and instead it increases the Vmax of SERCA. Note that the manuscript deposited in BioRxiv by the Veglia lab (https://www.biorxiv.org/content/10.1101/2021.05.05.442831v1) also shows that DWORF increases the Vmax of SERCA.

Reviewer #3:

In this study, Fisher et al. investigated the interaction of the sarco-endoplasmic reticulum calcium pump (SERCA) and the SERCA-binding protein small dwarf open reading frame (DWORF). Previous studies have shown that DWORF does not directly activate the SERCA pump and does not change the apparent affinity of SERCA for Ca+2. The current consensus in the field is that DWORF regulates the relative affinity of SERCA for Ca2+ by relieving the inhibition of PLN. In contrast to previous studies, the authors claim that DWORF directly regulates SERCA activity and calcium reuptake that is independent of its known role in relieving the inhibition of PLN. However, this major claim is not supported by the data presented in the manuscript.

We would like to point out that our data are completely consistent with the current model that DWORF competes with PLN for SERCA. The previous studies of the DWORF functional effect on SERCA only examined the effect on the KCa of SERCA. Our results agree with their observation; we also see no effect on KCa. Importantly, those previous studies did not look at the maximal activity (Vmax) of SERCA, thus our discovery of a direct effect on Vmax by DWORF are not in conflict with previous reports (1,2). In addition, the in-cell studies are unique and have not been performed by others to the best of our knowledge. The in-cell studies are also completely consistent with previous reports (1,2).

This same concern was raised by reviewer #2, please refer to question 3 above and our response. Our results confirm the previous observation that DWORF does not alter SERCA calcium affinity (2,3) and we have added two new mechanistic insights: (i) DWORF and PLN bind to SERCA in a calcium-dependent manner and (ii) DWORF directly increases SERCA maximal activity (Vmax). In support of these proposed novel mechanisms, we provide biochemical measurements of SERCA ATPase activity and unique, quantitative in-cell measurements of SERCA calcium transport kinetics.

Comments for the authors:

1. In Figure 2b it is shown that DWORF directly increases SERCA activity. Is this effect dose depended? Does the molar ratio of DWORF/SERCA affects the Vmax? It is also important to show the regulation of SERCA activity by DWORF in the presence and absence of PLN, phosphorylated PLN and SLN.

Indeed, we have done a dose response of DWORF in the presence of SERCA. At equimolar ratios of SERCA-DWORF, DWORF directly activates SERCA as presented. Based on the available literature, this condition best represents the physiological situation. Under extreme, non-physiological conditions of excess DWORF, we have observed modest SERCA inhibition manifested as a decrease in calcium affinity and a decrease in maximal activity. We have chosen not to include this latter data as it is certainly non-physiological, and it will confuse readers. The revised manuscript emphasizes the importance of controlling this parameter. The following statement has been added to the Methods section (page 10):

“To ensure the proper ratio of ~2 DWORF per SERCA in the reconstituted proteoliposomes, quantitative SDS-PAGE was used with known quantities of SERCA and DWORF proteins as standards (17, 18, 50).”

For the second point, our goal was to elucidate a direct effect of DWORF on SERCA and our results clearly support this novel mechanism. It is generally accepted that DWORF displaces PLN and our goal was not to reiterate the conclusions of previous reports (1,2).

2. Related to the previous comment, there is a leftward shift of the Ca2+ affinity curve in the presence of DWORF (Figure 2b), suggesting that the affinity for Ca2+ is increased. Is this significant? How do the authors explain this effect?

There is a small rightward shift (not a leftward shift) of the calcium affinity curve in the presence of DWORF (in Figure 2b), but this is not statistically significant.

3. In Figure 2a the data shows that PLN increases the Vmax of SERCA. This is puzzling as PLN is known to regulate the affinity os SERCA for Ca2+ but does not affect the Vmax. Is there an explanation for this observation? Also, the dose response curve and error bars should be shown on the graph.

SERCA and PLN are present together in cardiac SR membranes at molar ratios ranging from 1:1 to 1:4. At the higher molar ratios (1:4), PLN inhibits SERCA at low calcium concentrations (effect on KCa) and it increases the maximal activity of SERCA at higher calcium concentrations. The Young lab has previously studied this effect (10-15). Figure 2a showed fitted curves from previously published data and it was presented for comparison with the new SERCA-DWORF data. In retrospect, the effect of PLN on SERCA at different molar ratios is not important for this study and is confusing for the reader. We have modified the figure to show the effect of PLN on SERCA at an equimolar 1:1 ratio, which will better compare to the data for equimolar ratio of DWORF and SERCA.

4. The dose response curve of the small molecule CDN1163 should be shown in Figure 2C.

The CDN1163 compound has been thoroughly studied and published by others (16). We used this compound in our system as a control for comparison with the activation of SERCA by DWORF.

5. The authors used the skeletal muscle SERCA isoform SERCA1b in the reconstituted vesicles system and cardiac isoform, Serca2a, in the cell-based experiments. Are there any particular reasons why different isoforms were used in different experiments? Does DWORF regulate the activity of all SERCA isoforms?

The reconstitution of SERCA for quantitative ATPase measurements requires a lot of protein, which is readily available from rabbit skeletal muscle. Because it expresses high levels of the SERCA1a isoform, it is the standard source of enzyme in the field. The two isoforms (SERCA1a and 2a) are 84% identical, and their molecular structures have been determined to be identical (17). We have previously shown that PLN regulates SERCA1a and SERCA2a in a similar manner (18). Most important, DWORF has been shown to co-immunoprecipitate with all SERCA isoforms (2).

6. Figure 3. How do the authors explain the discrepancies between the findings of the current study and the previous studies? For example, using a similar cellular model of DWORF overexpression (Cos7 cells co-transfected with SERCA2a) it has been reported that cells expressing SERCA2a and DWORF alone do not exhibit enhanced SERCA activity indicating that DWORF itself does not biochemically activate SERCA. Additional experiments are needed to clary this discrepancy. What is the effect of DWORF on SERCA2a in the presence of PLN and/or SLN in the cellular model used here? Is the Ca2+ uptake kinetics affected by the stoichiometry of the overexpressed proteins? As the levels of the overexpressed proteins are not titrated, it is likely that these observations represent artifacts associated with protein overexpression.

Please see our responses to similar questions above. In our revision, we have explained more clearly that there is no discrepancy. In a previous study, COS7 cells co-transfected with SERCA2a and DWORF were used for co-immunoprecipitation studies and measurement of the calcium affinity of SERCA (2). To quote this study “We found that coexpression of DWORF alone with SERCA2a did not change the apparent affinity of SERCA for Ca2+, but it relieved the inhibition by PLN in a dose-dependent manner (Figure S15).” Thus, these studies focused on the KCa of SERCA and not the Vmax. Our results are completely consistent with these prior studies by the Olson laboratory (2).

7. Figure 4a. The calcium transient frequency is significantly increased in the presence of DWORF. This suggest an increased RYR activity in the presence of DOWRF. Is there an indirect/direct effect of the RyR2 overexpression on the observed SERCA Ca2+ kinetics? Is RyR2 necessary for the regulation of SERCA by DWORF?

Our data shows that DWORF increases the ER calcium content, which is dependent on the presence of SERCA. As the ER calcium content increases, the frequency of spontaneous calcium release mediated by RyR2 also increases. A similar effect is observed in cardiomyocytes. Calcium transient (calcium wave) frequency is increased in cells expressing DWORF (Figure 4C). This effect is mediated by increased SERCA activity and the resulting increase in ER calcium load. Indeed, it is well known that an increase in SR calcium load during SERCA activation (adrenergic activation) increases the frequency of spontaneous calcium transients in myocytes.

We can rule out an effect of DWORF on RyR, since it has been shown by many labs (including (19)) that RyR activation would only transiently increase calcium transient frequency. As SR calcium gets depleted during RyR activation, calcium transients would be reduced to control levels. Sustained increases in spontaneous calcium transient frequency can only be achieved by increased SERCA activity.

Moreover, in vitro ATPase measurements show that RyR is not required for SERCA activation by DWORF. We express RyR for the live cell calcium uptake assay as a way to manipulate ER calcium load (by producing spontaneous CICR or by activating/inhibiting RyR), so SERCA function can be studied throughout a wide range of ER calcium loads.

8. Figure 4b. Ideally the Ca2+ reuptake kinetics, eg Tau, should be measured under electrically stimulated conditions as the frequency of the Ca2+ release will affected the kinetics. It is likely that the increase of Tau by DWORF reflect the increase of the oscillation frequency.

The HEK cells expressing SERCA, RyR, DWORF and PLN are a useful incell cardiomimetic system; however, it is not possible to do electrical stimulation of HEK cells. Calcium release frequency only affects uptake kinetics insofar as it may affect the ER calcium load, a parameter that is carefully quantified here. An important aspect of our in-cell experiments is the calibration of the fluorescent ER calcium indicator at the end of each experiment. ER content calibration allows us to specifically determine uptake rate as a function of ER calcium load. We also control the cytosolic calcium concentration (buffered to 100 nM), and we quantify ER leak rate independent of SERCA and RyR. Our unique approach yields quantitative measurements of maximal uptake rate and maximal ER calcium load as intrinsic features of SERCA function.

References:

1. Makarewich, C. A., Munir, A. Z., Schiattarella, G. G., Bezprozvannaya, S., Raguimova, O. N., Cho, E. E., Vidal, A. H., Robia, S. L., Bassel-Duby, R., and Olson, E. N. (2018) The DWORF micropeptide enhances contractility and prevents heart failure in a mouse model of dilated cardiomyopathy. eLife 7

2. Nelson, B. R., Makarewich, C. A., Anderson, D. M., Winders, B. R., Troupes, C. D., Wu, F., Reese,

A. L., McAnally, J. R., Chen, X., Kavalali, E. T., Cannon, S. C., Houser, S. R., Bassel-Duby, R., and Olson, E. N. (2016) A peptide encoded by a transcript annotated as long noncoding RNA enhances SERCA activity in muscle. Science 351, 271-275

3. Bidwell, P., Blackwell, D. J., Hou, Z., Zima, A. V., and Robia, S. L. (2011) Phospholamban binds with differential affinity to calcium pump conformers. J Biol Chem 286, 35044-35050

4. Dong, X., and Thomas, D. D. (2014) Time-resolved FRET reveals the structural mechanism of SERCA-PLB regulation. Biochem Biophys Res Commun 449, 196-201

5. Mueller, B., Karim, C. B., Negrashov, I. V., Kutchai, H., and Thomas, D. D. (2004) Direct detection of phospholamban and sarcoplasmic reticulum Ca-ATPase interaction in membranes using fluorescence resonance energy transfer. Biochemistry 43, 8754-8765

6. Li, J., Bigelow, D. J., and Squier, T. C. (2004) Conformational changes within the cytosolic portion of phospholamban upon release of Ca-ATPase inhibition. Biochemistry 43, 3870-3879

7. Akin, B. L., Hurley, T. D., Chen, Z., and Jones, L. R. (2013) The structural basis for phospholamban inhibition of the calcium pump in sarcoplasmic reticulum. J Biol Chem 288, 30181-30191

8. Toyoshima, C., Iwasawa, S., Ogawa, H., Hirata, A., Tsueda, J., and Inesi, G. (2013) Crystal structures of the calcium pump and sarcolipin in the Mg2+-bound E1 state. Nature 495, 260-264

9. Winther, A. M., Bublitz, M., Karlsen, J. L., Moller, J. V., Hansen, J. B., Nissen, P., and BuchPedersen, M. J. (2013) The sarcolipin-bound calcium pump stabilizes calcium sites exposed to the cytoplasm. Nature 495, 265-269

10. Glaves, J. P., Primeau, J. O., Gorski, P. A., Espinoza-Fonseca, L. M., Lemieux, M. J., and Young, H.

S. (2020) Interaction of a Sarcolipin Pentamer and Monomer with the Sarcoplasmic Reticulum Calcium Pump, SERCA. Biophys J 118, 518-531

11. Glaves, J. P., Primeau, J. O., Espinoza-Fonseca, L. M., Lemieux, M. J., and Young, H. S. (2019) The Phospholamban Pentamer Alters Function of the Sarcoplasmic Reticulum Calcium Pump SERCA. Biophys J 116, 633-647

12. Ceholski, D. K., Trieber, C. A., and Young, H. S. (2012) Hydrophobic imbalance in the cytoplasmic domain of phospholamban is a determinant for lethal dilated cardiomyopathy. J Biol Chem 287, 16521-16529

13. Glaves, J. P., Trieber, C. A., Ceholski, D. K., Stokes, D. L., and Young, H. S. (2011) Phosphorylation and mutation of phospholamban alter physical interactions with the sarcoplasmic reticulum calcium pump. J Mol Biol 405, 707-723

14. Trieber, C. A., Afara, M., and Young, H. S. (2009) Effects of phospholamban transmembrane mutants on the calcium affinity, maximal activity, and cooperativity of the sarcoplasmic reticulum calcium pump. Biochemistry 48, 9287-9296

15. Trieber, C. A., Douglas, J. L., Afara, M., and Young, H. S. (2005) The effects of mutation on the regulatory properties of phospholamban in co-reconstituted membranes. Biochemistry 44, 3289-3297

16. Kang, S., Dahl, R., Hsieh, W., Shin, A., Zsebo, K. M., Buettner, C., Hajjar, R. J., and Lebeche, D.

(2016) Small Molecular Allosteric Activator of the Sarco/Endoplasmic Reticulum Ca2+-ATPase

(SERCA) Attenuates Diabetes and Metabolic Disorders. J Biol Chem 291, 5185-5198

17. Sitsel, A., De Raeymaecker, J., Drachmann, N. D., Derua, R., Smaardijk, S., Andersen, J. L., Vandecaetsbeek, I., Chen, J., De Maeyer, M., Waelkens, E., Olesen, C., Vangheluwe, P., and

Nissen, P. (2019) Structures of the heart specific SERCA2a Ca(2+)-ATPase. EMBO J 38

18. Gorski, P. A., Trieber, C. A., Ashrafi, G., and Young, H. S. (2015) Regulation of the sarcoplasmic reticulum calcium pump by divergent phospholamban isoforms in zebrafish. J Biol Chem 290, 6777-6788

19. Trafford, A. W., Diaz, M. E., Sibbring, G. C., and Eisner, D. A. (2000) Modulation of CICR has no maintained effect on systolic Ca2+: simultaneous measurements of sarcoplasmic reticulum and sarcolemmal Ca2+ fluxes in rat ventricular myocytes. J Physiol 522 Pt 2, 259-270

https://doi.org/10.7554/eLife.65545.sa2

Article and author information

Author details

  1. M'Lynn E Fisher

    Department of Biochemistry, University of Alberta, Edmonton, Canada
    Contribution
    Conceptualization, Data curation, Formal analysis, Investigation, Methodology
    Competing interests
    No competing interests declared
  2. Elisa Bovo

    Department of Cell and Molecular Physiology, Stritch School of Medicine, Loyola University Chicago, Maywood, United States
    Contribution
    Formal analysis, Investigation, Methodology, Writing - review and editing
    Contributed equally with
    Rodrigo Aguayo-Ortiz
    Competing interests
    No competing interests declared
  3. Rodrigo Aguayo-Ortiz

    Center for Arrhythmia Research, Department of Internal Medicine, Division of Cardiovascular Medicine, University of Michigan, Ann Arbor, United States
    Contribution
    Formal analysis, Investigation, Methodology, Writing - review and editing
    Contributed equally with
    Elisa Bovo
    Competing interests
    No competing interests declared
  4. Ellen E Cho

    Department of Cell and Molecular Physiology, Stritch School of Medicine, Loyola University Chicago, Maywood, United States
    Contribution
    Formal analysis, Investigation
    Competing interests
    No competing interests declared
  5. Marsha P Pribadi

    Department of Cell and Molecular Physiology, Stritch School of Medicine, Loyola University Chicago, Maywood, United States
    Contribution
    Formal analysis, Investigation, Methodology
    Competing interests
    No competing interests declared
  6. Michael P Dalton

    Department of Cell and Molecular Physiology, Stritch School of Medicine, Loyola University Chicago, Maywood, United States
    Contribution
    Formal analysis, Investigation, Methodology
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-5296-5099
  7. Nishadh Rathod

    Department of Biochemistry, University of Alberta, Edmonton, Canada
    Contribution
    Investigation, Methodology
    Competing interests
    No competing interests declared
  8. M Joanne Lemieux

    Department of Biochemistry, University of Alberta, Edmonton, Canada
    Contribution
    Resources, Methodology, Writing - review and editing
    Competing interests
    Reviewing editor, eLife
  9. L Michel Espinoza-Fonseca

    Center for Arrhythmia Research, Department of Internal Medicine, Division of Cardiovascular Medicine, University of Michigan, Ann Arbor, United States
    Contribution
    Formal analysis, Funding acquisition, Investigation, Methodology, Writing - review and editing
    Competing interests
    No competing interests declared
  10. Seth L Robia

    Department of Cell and Molecular Physiology, Stritch School of Medicine, Loyola University Chicago, Maywood, United States
    Contribution
    Conceptualization, Data curation, Formal analysis, Funding acquisition, Methodology, Writing - review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1193-9510
  11. Aleksey V Zima

    Department of Cell and Molecular Physiology, Stritch School of Medicine, Loyola University Chicago, Maywood, United States
    Contribution
    Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Writing - review and editing
    Competing interests
    No competing interests declared
  12. Howard S Young

    Department of Biochemistry, University of Alberta, Edmonton, Canada
    Contribution
    Conceptualization, Resources, Data curation, Formal analysis, Supervision, Funding acquisition, Investigation, Visualization, Methodology, Writing - original draft, Project administration, Writing - review and editing
    For correspondence
    hyoung@ualberta.ca
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5990-8422

Funding

National Institutes of Health (R01HL092321)

  • Seth L Robia
  • Howard S Young

National Institutes of Health (R01HL143816)

  • Seth L Robia
  • Howard S Young

National Institutes of Health (R01GM120142)

  • L Michel Espinoza-Fonseca

National Institutes of Health (R01HL148068)

  • L Michel Espinoza-Fonseca

National Institutes of Health (R01HL130231)

  • Aleksey V Zima

Natural Sciences and Engineering Research Council of Canada (RGPIN-2016-06478)

  • M Joanne Lemieux

Heart and Stroke Foundation of Canada

  • Howard S Young

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

This research was supported in part through computational resources and services provided by Advanced Research Computing at the University of Michigan, Ann Arbor, and Compute Canada (http://www.computecanada.ca). Funding: This work was supported in part by grants from the National Institutes of Health (R01HL092321 and R01HL143816 to HSY and SLR; R01GM120142 and R01HL148068 to LMEF; and R01HL130231 to AVZ), from the Heart and Stroke Foundation of Canada (to HSY), and the Natural Sciences and Engineering Research Council of Canada (RGPIN-2016-06478 to MJL).

Senior Editor

  1. Richard W Aldrich, The University of Texas at Austin, United States

Reviewing Editor

  1. Mark T Nelson, University of Vermont, United States

Reviewer

  1. Amreen Mughal, University of Vermont, United States

Version history

  1. Received: December 7, 2020
  2. Accepted: June 1, 2021
  3. Accepted Manuscript published: June 2, 2021 (version 1)
  4. Version of Record published: June 14, 2021 (version 2)

Copyright

© 2021, Fisher et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,999
    Page views
  • 274
    Downloads
  • 19
    Citations

Article citation count generated by polling the highest count across the following sources: PubMed Central, Crossref, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. M'Lynn E Fisher
  2. Elisa Bovo
  3. Rodrigo Aguayo-Ortiz
  4. Ellen E Cho
  5. Marsha P Pribadi
  6. Michael P Dalton
  7. Nishadh Rathod
  8. M Joanne Lemieux
  9. L Michel Espinoza-Fonseca
  10. Seth L Robia
  11. Aleksey V Zima
  12. Howard S Young
(2021)
Dwarf open reading frame (DWORF) is a direct activator of the sarcoplasmic reticulum calcium pump SERCA
eLife 10:e65545.
https://doi.org/10.7554/eLife.65545

Further reading

    1. Biochemistry and Chemical Biology
    2. Cell Biology
    Maria Körner, Susanne R Meyer ... Alexander Buchberger
    Research Article Updated

    The ATPase p97 (also known as VCP, Cdc48) has crucial functions in a variety of important cellular processes such as protein quality control, organellar homeostasis, and DNA damage repair, and its de-regulation is linked to neuromuscular diseases and cancer. p97 is tightly controlled by numerous regulatory cofactors, but the full range and function of the p97–cofactor network is unknown. Here, we identify the hitherto uncharacterized FAM104 proteins as a conserved family of p97 interactors. The two human family members VCP nuclear cofactor family member 1 and 2 (VCF1/2) bind p97 directly via a novel, alpha-helical motif and associate with p97-UFD1-NPL4 and p97-UBXN2B complexes in cells. VCF1/2 localize to the nucleus and promote the nuclear import of p97. Loss of VCF1/2 results in reduced nuclear p97 levels, slow growth, and hypersensitivity to chemical inhibition of p97 in the absence and presence of DNA damage, suggesting that FAM104 proteins are critical regulators of nuclear p97 functions.

    1. Biochemistry and Chemical Biology
    2. Epidemiology and Global Health
    Takashi Sasaki, Yoshinori Nishimoto ... Yasumichi Arai
    Research Article

    Background: High levels of circulating adiponectin are associated with increased insulin sensitivity, low prevalence of diabetes, and low body mass index (BMI); however, high levels of circulating adiponectin are also associated with increased mortality in the 60-70 age group. In this study, we aimed to clarify factors associated with circulating high-molecular-weight (cHMW) adiponectin levels and their association with mortality in the very old (85-89 years old) and centenarians.

    Methods: The study included 812 (women: 84.4%) for centenarians and 1,498 (women: 51.7%) for the very old. The genomic DNA sequence data were obtained by whole genome sequencing or DNA microarray-imputation methods. LASSO and multivariate regression analyses were used to evaluate cHMW adiponectin characteristics and associated factors. All-cause mortality was analyzed in three quantile groups of cHMW adiponectin levels using Cox regression.

    Results: The cHMW adiponectin levels were increased significantly beyond 100 years of age, were negatively associated with diabetes prevalence, and were associated with SNVs in CDH13 (p = 2.21 × 10-22) and ADIPOQ (p = 5.72 × 10-7). Multivariate regression analysis revealed that genetic variants, BMI, and high-density lipoprotein cholesterol (HDLC) were the main factors associated with cHMW adiponectin levels in the very old, whereas the BMI showed no association in centenarians. The hazard ratios for all-cause mortality in the intermediate and high cHMW adiponectin groups in very old men were significantly higher rather than those for all-cause mortality in the low level cHMW adiponectin group, even after adjustment with BMI. In contrast, the hazard ratios for all-cause mortality were significantly higher for high cHMW adiponectin groups in very old women, but were not significant after adjustment with BMI.

    Conclusions: cHMW adiponectin levels increased with age until centenarians, and the contribution of known major factors associated with cHMW adiponectin levels, including BMI and HDLC, varies with age, suggesting that its physiological significance also varies with age in the oldest old.

    Funding: This study was supported by grants from the Ministry of Health, Welfare, and Labour for the Scientific Research Projects for Longevity; a Grant-in-Aid for Scientific Research (No 21590775, 24590898, 15KT0009, 18H03055, 20K20409, 20K07792, 23H03337) from the Japan Society for the Promotion of Science; Keio University Global Research Institute (KGRI), Kanagawa Institute of Industrial Science and Technology (KISTEC), Japan Science and Technology Agency (JST) Research Complex Program 'Tonomachi Research Complex' Wellbeing Research Campus: Creating new values through technological and social innovation (JP15667051), the Program for an Integrated Database of Clinical and Genomic Information from the Japan Agency for Medical Research and Development (No. 16kk0205009h001, 17jm0210051h0001, 19dk0207045h0001); the medical-welfare-food-agriculture collaborative consortium project from the Japan Ministry of Agriculture, Forestry, and Fisheries; and the Biobank Japan Program from the Ministry of Education, Culture, Sports, and Technology.