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

Phospholamban inhibits the cardiac calcium pump by interrupting an allosteric activation pathway

Phospholamban (PLB) is a transmembrane micropeptide that regulates the sarcoplasmic reticulum Ca2+-ATPase (SERCA) in cardiac muscle, but the physical mechanism of this regulation remains poorly understood. PLB reduces the Ca2+ sensitivity of active SERCA, increasing the Ca2+ concentration required for pump cycling. However, PLB does not decrease Ca2+ binding to SERCA when ATP is absent, suggesting PLB does not inhibit SERCA Ca2+ affinity. The prevailing explanation for these seemingly conflicting results is that PLB slows transitions in the SERCA enzymatic cycle associated with Ca2+ binding, altering transport Ca2+ dependence without actually affecting the equilibrium binding affinity of the Ca2+-coordinating sites. Here, we consider another hypothesis, that measurements of Ca2+ binding in the absence of ATP overlook important allosteric effects of nucleotide binding that increase SERCA Ca2+ binding affinity. We speculated that PLB inhibits SERCA by reversing this allostery. To test this, we used a fluorescent SERCA biosensor to quantify the Ca2+ affinity of non-cycling SERCA in the presence and absence of a non-hydrolyzable ATP-analog, AMPPCP. Nucleotide activation increased SERCA Ca2+ affinity, and this effect was reversed by co-expression of PLB. Interestingly, PLB had no effect on Ca2+ affinity in the absence of nucleotide. These results reconcile the previous conflicting observations from ATPase assays versus Ca2+ binding assays. Moreover, structural analysis of SERCA revealed a novel allosteric pathway connecting the ATP- and Ca2+-binding sites. We propose this pathway is disrupted by PLB binding. Thus, PLB reduces the equilibrium Ca2+ affinity of SERCA by interrupting allosteric activation of the pump by ATP.

Multiscale computational modeling of the effects of 2'-deoxy-ATP on cardiac muscle calcium handling

2'-Deoxy-ATP (dATP), a naturally occurring near analog of ATP, is a well-documented myosin activator that has been shown to increase contractile force, improve pump function, and enhance lusitropy in the heart. Calcium transients in cardiomyocytes with elevated levels of dATP show faster calcium decay compared with cardiomyocytes with basal levels of dATP, but the mechanisms behind this are unknown. Here, we design and utilize a multiscale computational modeling framework to test the hypothesis that dATP acts on the sarcoendoplasmic reticulum calcium-ATPase (SERCA) pump to accelerate calcium re-uptake into the sarcoplasmic reticulum during cardiac relaxation. Gaussian accelerated molecular dynamics simulations of human cardiac SERCA2A in the E1 apo, ATP-bound and dATP-bound states showed that dATP forms more stable contacts in the nucleotide binding pocket of SERCA and leads to increased closure of cytosolic domains. These structural changes ultimately lead to changes in calcium binding, which we assessed using Brownian dynamics simulations. We found that dATP increases calcium association rate constants to SERCA and that dATP binds to apo SERCA more rapidly than ATP. Using a compartmental ordinary differential equation model of human cardiomyocyte excitation-contraction coupling, we found that these increased association rate constants contributed to the accelerated rates of calcium transient decay observed experimentally. This study provides clear mechanistic evidence of enhancements in cardiac SERCA2A pump function due to interactions with dATP.

Structural basis of the conformational and functional regulation of human SERCA2b, the ubiquitous endoplasmic reticulum calcium pump

Sarco/endoplasmic reticulum Ca²⁺ ATPase 2b (SERCA2b), a member of the SERCA family, is expressed ubiquitously and transports Ca²⁺ into the sarco/endoplasmic reticulum using the energy provided by ATP binding and hydrolysis. The crystal structure of SERCA2b in its Ca²⁺ - and ATP-bound (E1∙2Ca²⁺ -ATP) state and cryo-electron microscopy (cryo-EM) structures of the protein in its E1∙2Ca²⁺ -ATP and Ca²⁺ -unbound phosphorylated (E2P) states have provided essential insights into how the overall conformation and ATPase activity of SERCA2b is regulated by the transmembrane helix 11 and the subsequent luminal extension loop, both of which are specific to this isoform. More recently, our cryo-EM analysis has revealed that SERCA2b likely adopts open and closed conformations of the cytosolic domains in the Ca²⁺ -bound but ATP-free (E1∙2Ca²⁺ ) state, and that the closed conformation represents a state immediately prior to ATP binding. This review article summarizes the unique mechanisms underlying the conformational and functional regulation of SERCA2b.

Retrograde Analysis of Calcium Signaling by CaMPARI2 Shows Cytosolic Calcium in Chondrocytes Is Unaffected by Parabolic Flights

Calcium (Ca2+) elevation is an essential secondary messenger in many cellular processes, including disease progression and adaptation to external stimuli, e.g., gravitational load. Therefore, mapping and quantifying Ca2+ signaling with a high spatiotemporal resolution is a key challenge. However, particularly on microgravity platforms, experiment time is limited, allowing only a small number of replicates. Furthermore, experiment hardware is exposed to changes in gravity levels, causing experimental artifacts unless appropriately controlled. We introduce a new experimental setup based on the fluorescent Ca2+ reporter CaMPARI2, onboard LED arrays, and subsequent microscopic analysis on the ground. This setup allows for higher throughput and accuracy due to its retrograde nature. The excellent performance of CaMPARI2 was demonstrated with human chondrocytes during the 75th ESA parabolic flight campaign. CaMPARI2 revealed a strong Ca2+ response triggered by histamine but was not affected by the alternating gravitational load of a parabolic flight.

Cryo-EM analysis provides new mechanistic insight into ATP binding to Ca2+ -ATPase SERCA2b

Sarco/endoplasmic reticulum Ca²⁺ -ATPase (SERCA) 2b is a ubiquitous SERCA family member that conducts Ca²⁺ uptake from the cytosol to the ER. Herein, we present a 3.3 Å resolution cryo-electron microscopy (cryo-EM) structure of human SERCA2b in the E1·2Ca²⁺ state, revealing a new conformation for Ca²⁺ -bound SERCA2b with a much closer arrangement of cytosolic domains than in the previously reported crystal structure of Ca²⁺ -bound SERCA1a. Multiple conformations generated by 3D classification of cryo-EM maps reflect the intrinsically dynamic nature of the cytosolic domains in this state. Notably, ATP binding residues of SERCA2b in the E1·2Ca²⁺ state are located at similar positions to those in the E1·2Ca²⁺ -ATP state; hence, the cryo-EM structure likely represents a preformed state immediately prior to ATP binding. Consistently, a SERCA2b mutant with an interdomain disulfide bridge that locks the closed cytosolic domain arrangement displayed significant autophosphorylation activity in the presence of Ca²⁺ . We propose a novel mechanism of ATP binding to SERCA2b.

Structural Changes beyond the EF-Hand Contribute to Apparent Calcium Binding Affinities: Insights from Parvalbumins

Members of the parvalbumin (PV) family of calcium (Ca2+) binding proteins (CBPs) share a relatively high level of sequence similarity. However, their Ca2+ affinities and selectivities against competing ions like Mg2+ can widely vary. We conducted molecular dynamics simulations of several α-parvalbumin (αPV) constructs with micromolar to nanomolar Ca2+ affinities to identify structural and dynamic features that contribute to their binding of ions. Specifically, we examined a D94S/G98E construct with a lower Ca2+ affinity (≈-18 kcal/mol) relative to the wild type (WT) (≈-22 kcal/mol) and an S55D/E59D variant with enhanced affinity (≈-24 kcal/mol). Additionally, we also examined the binding of Mg2+ to these isoforms, which is much weaker than Ca2+. We used mean spherical approximation (MSA) theory to evaluate ion binding thermodynamics within the proteins' EF-hand domains to account for the impact of ions' finite sizes and the surrounding electrolyte composition. While the MSA scores differentiated Mg2+ from Ca2+, they did not indicate that Ca2+ binding affinities at the binding loop differed between the PV isoforms. Instead, molecular mechanics generalized Born surface area (MM/GBSA) approximation energies, which we used to quantify the thermodynamic cost of structural rearrangement of the proteins upon binding ions, indicated that S55D/E59D αPV favored Ca2+ binding by -20 kcal/mol relative to WT versus 30 kcal/mol for D94S/G98E αPV. Meanwhile, Mg2+ binding was favored for the S55D/E59D αPV and D94S/G98E αPV variants by -18.32 and -1.65 kcal/mol, respectively. These energies implicate significant contributions to ion binding beyond oxygen coordination at the binding loop, which stemmed from changes in α-helicity, β-sheet character, and hydrogen bonding. Hence, Ca2+ affinity and selectivity against Mg2+ are emergent properties stemming from both local effects within the proteins' ion binding sites as well as non-local contributions elsewhere. Our findings broaden our understanding of the molecular bases governing αPV ion binding that are likely shared by members of the broad family of CBPs.

Dynamics-Driven Allostery Underlies Ca2+-Mediated Release of SERCA Inhibition by Phospholamban

Sarcoplasmic reticulum (SR) Ca2+-ATPase (SERCA) and phospholamban (PLB) are essential for intracellular Ca2+ transport in myocytes. Ca2+-dependent activation of SERCA-PLB provides a control function that regulates cytosolic and SR Ca2+ levels. Although experimental and computational studies alone have led to a greater insight into SERCA-PLB regulation, the structural mechanisms for Ca2+ binding reversing inhibition of the complex remain poorly understood. Therefore, we have performed atomistic simulations totaling 32.7 μs and cell-based intramolecular fluorescence resonance energy transfer (FRET) experiments to determine structural changes of PLB-bound SERCA in response to binding of a single Ca2+ ion. Complementary MD simulations and FRET experiments showed that open-to-closed transitions in the structure of the headpiece underlie PLB inhibition of SERCA, and binding of a single Ca2+ ion is sufficient to shift the protein population toward a structurally closed structure of the complex. Closure is accompanied by functional interactions between the N-domain β5-β6 loop and the A-domain and the displacement of the catalytic phosphorylation domain toward a competent structure. We propose that reversal of SERCA-PLB inhibition is achieved by stringing together its controlling modules (A-domain and loop Nβ5-β6) with catalytic elements (P-domain) to regulate function during intracellular Ca2+ signaling. We conclude that binding of a single Ca2+ is a critical mediator of allosteric signaling that dictates structural changes and motions that relieve SERCA inhibition by PLB. Understanding allosteric regulation is of paramount importance to guide therapeutic modulation of SERCA and other evolutionarily related ion-motive ATPases.

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

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

Endoplasmic Reticulum Calcium Pumps and Tumor Cell Differentiation

Endoplasmic reticulum (ER) calcium homeostasis plays an essential role in cellular calcium signaling, intra-ER protein chaperoning and maturation, as well as in the interaction of the ER with other organelles. Calcium is accumulated in the ER by sarco/endoplasmic reticulum calcium ATPases (SERCA enzymes) that generate by active, ATP-dependent transport, a several thousand-fold calcium ion concentration gradient between the cytosol (low nanomolar) and the ER lumen (high micromolar). SERCA enzymes are coded by three genes that by alternative splicing give rise to several isoforms, which can display isoform-specific calcium transport characteristics. SERCA expression levels and isoenzyme composition vary according to cell type, and this constitutes a mechanism whereby ER calcium homeostasis is adapted to the signaling and metabolic needs of the cell, depending on its phenotype, its state of activation and differentiation. As reviewed here, in several normal epithelial cell types including bronchial, mammary, gastric, colonic and choroid plexus epithelium, as well as in mature cells of hematopoietic origin such as pumps are simultaneously expressed, whereas in corresponding tumors and leukemias SERCA3 expression is selectively down-regulated. SERCA3 expression is restored during the pharmacologically induced differentiation of various cancer and leukemia cell types. SERCA3 is a useful marker for the study of cell differentiation, and the loss of SERCA3 expression constitutes a previously unrecognized example of the remodeling of calcium homeostasis in tumors.

Structural dynamics of P-type ATPase ion pumps

P-type ATPases transport ions across biological membranes against concentration gradients and are essential for all cells. They use the energy from ATP hydrolysis to propel large intramolecular movements, which drive vectorial transport of ions. Tight coordination of the motions of the pump is required to couple the two spatially distant processes of ion binding and ATP hydrolysis. Here, we review our current understanding of the structural dynamics of P-type ATPases, focusing primarily on Ca2+ pumps. We integrate different types of information that report on structural dynamics, primarily time-resolved fluorescence experiments including single-molecule Förster resonance energy transfer and molecular dynamics simulations, and interpret them in the framework provided by the numerous crystal structures of sarco/endoplasmic reticulum Ca2+-ATPase. We discuss the challenges in characterizing the dynamics of membrane pumps, and the likely impact of new technologies on the field.

Thermodynamics of Cation Binding to the Sarcoendoplasmic Reticulum Calcium ATPase Pump and Impacts on Enzyme Function

Sarcoendoplasmic reticulum Ca²⁺-ATPase (SERCA) is a transmembrane pump that plays an important role in transporting calcium into the sarcoplasmic reticulum (SR). While calcium (Ca²⁺) binds SERCA with micromolar affinity, magnesium (Mg²⁺) and potassium (K⁺) also compete with Ca²⁺ binding. However, the molecular bases for these competing ions' influence on the SERCA function and the selectivity of the pump for Ca²⁺ are not well-established. We therefore used in silico methods to resolve molecular determinants of cation binding in the canonical site I and II Ca²⁺ binding sites via (1) triplicate molecular dynamics (MD) simulations of Mg²⁺, Ca²⁺, and K⁺-bound SERCA, (2) mean spherical approximation (MSA) theory to score the affinity and selectivity of cation binding to the MD-resolved structures, and (3) state models of SERCA turnover informed from MSA-derived affinity data. Our key findings are that (a) coordination at sites I and II is optimized for Ca²⁺ and to a lesser extent for Mg²⁺ and K⁺, as determined by MD-derived cation-amino acid oxygen and bound water configurations, (b) the impaired coordination and high desolvation cost for Mg²⁺ precludes favorable Mg²⁺ binding relative to Ca²⁺, while K⁺ has limited capacity to bind site I, and (c) Mg²⁺ most likely acts as inhibitor and K⁺ as intermediate in SERCA's reaction cycle, based on a best-fit state model of SERCA turnover. These findings provide a quantitative basis for SERCA function that leverages molecular-scale thermodynamic data and rationalizes enzyme activity across broad ranges of K⁺, Ca²⁺, and Mg²⁺ concentrations.

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

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

Prediction and validation of protein intermediate states from structurally rich ensembles and coarse-grained simulations

Protein conformational changes are at the heart of cell functions, from signalling to ion transport. However, the transient nature of the intermediates along transition pathways hampers their experimental detection, making the underlying mechanisms elusive. Here we retrieve dynamic information on the actual transition routes from principal component analysis (PCA) of structurally-rich ensembles and, in combination with coarse-grained simulations, explore the conformational landscapes of five well-studied proteins. Modelling them as elastic networks in a hybrid elastic-network Brownian dynamics simulation (eBDIMS), we generate trajectories connecting stable end-states that spontaneously sample the crystallographic motions, predicting the structures of known intermediates along the paths. We also show that the explored non-linear routes can delimit the lowest energy passages between end-states sampled by atomistic molecular dynamics. The integrative methodology presented here provides a powerful framework to extract and expand dynamic pathway information from the Protein Data Bank, as well as to validate sampling methods in general.

Protein conformational changes are key to a wide range of cellular functions but remain difficult to access experimentally. Here the authors describe eBDIMS, a novel approach to predict intermediates observed in structural transition pathways from experimental ensembles.

Understanding Ion Binding Affinity and Selectivity in β-Parvalbumin Using Molecular Dynamics and Mean Spherical Approximation Theory

Parvalbumin (PV) is a globular calcium (Ca(2+))-selective protein expressed in a variety of biological tissues. Our computational studies of the rat β-parvalbumin (β-PV) isoform seek to elucidate the molecular thermodynamics of Ca(2+) versus magnesium (Mg(2+)) binding at the protein's two EF-hand motifs. Specifically, we have utilized molecular dynamics (MD) simulations and a mean-field electrolyte model (mean spherical approximation (MSA) theory) to delineate how the EF-hand scaffold controls the "local" thermodynamics of Ca(2+) binding selectivity over Mg(2+). Our MD simulations provide the probability density of metal-chelating oxygens within the EF-hand scaffolds for both Ca(2+) and Mg(2+), as well the conformational strain induced by Mg(2+) relative to Ca(2+) binding. MSA theory utilizes the binding domain oxygen and charge distributions to predict the chemical potential of ion binding, as well as their corresponding concentrations within the binding domain. We find that the electrostatic and steric contributions toward ion binding were similar for Mg(2+) and Ca(2+), yet the latter was 5.5 kcal/mol lower in enthalpy when internal strain within the EF hand was considered. We therefore speculate that beyond differences in dehydration energies for the Ca(2+) versus Mg(2+), strain induced in the β-PV EF hand by cation binding significantly contributes to the nearly 10,000-fold difference in binding affinity reported in the literature. We further complemented our analyses of local factors governing cation binding selectivity with whole-protein (global) contributions, such as interhelical residue-residue contacts and solvent exposure of hydrophobic surface. These contributions were found to be comparable for both Ca(2+)- and Mg(2+)-bound β-PV, which may implicate local factors, EF-hand strain, and dehydration, in providing the primary means of selectivity. We anticipate these methods could be used to estimate metal binding thermodynamics across a broad range of PV sequence homologues and EF-hand-containing, Ca(2+) binding proteins.

A structural mechanism for calcium transporter headpiece closure

To characterize the conformational dynamics of sarcoplasmic reticulum (SR) calcium pump (SERCA) we performed molecular dynamics simulations beginning with several different high-resolution structures. We quantified differences in structural disorder and dynamics for an open conformation of SERCA versus closed structures and observed that dynamic motions of SERCA cytoplasmic domains decreased with decreasing domain–domain separation distance. The results are useful for interpretation of recent intramolecular Förster resonance energy transfer (FRET) distance measurements obtained for SERCA fused to fluorescent protein tags. Those previous physical measurements revealed several discrete structural substates and suggested open conformations of SERCA are more dynamic than compact conformations. The present simulations support this hypothesis and provide additional details of SERCA molecular mechanisms. Specifically, all-atoms simulations revealed large-scale translational and rotational motions of the SERCA N-domain relative to the A- and P-domains during the transition from an open to a closed headpiece conformation over the course of a 400 ns trajectory. The open-to-closed structural transition was accompanied by a disorder-to-order transition mediated by an initial interaction of an N-domain loop (Nβ5-β6, residues 426–436) with residues 133–139 of the A-domain. Mutation of three negatively charged N-domain loop residues abolished the disorder-to-order transition and prevented the initial domain–domain interaction and subsequent closure of the cytoplasmic headpiece. Coarse-grained molecular dynamics simulations were in harmony with all-atoms simulations and physical measurements and revealed a close communication between fluorescent protein tags and the domain to which they were fused. The data indicate that previous intramolecular FRET distance measurements report SERCA structure changes with high fidelity and suggest a structural mechanism that facilitates the closure of the SERCA cytoplasmic headpiece.

POVME 2.0: An Enhanced Tool for Determining Pocket Shape and Volume Characteristics

Analysis of macromolecular/small-molecule binding pockets can provide important insights into molecular recognition and receptor dynamics. Since its release in 2011, the POVME (POcket Volume MEasurer) algorithm has been widely adopted as a simple-to-use tool for measuring and characterizing pocket volumes and shapes. We here present POVME 2.0, which is an order of magnitude faster, has improved accuracy, includes a graphical user interface, and can produce volumetric density maps for improved pocket analysis. To demonstrate the utility of the algorithm, we use it to analyze the binding pocket of RNA editing ligase 1 from the unicellular parasite Trypanosoma brucei, the etiological agent of African sleeping sickness. The POVME analysis characterizes the full dynamics of a potentially druggable transient binding pocket and so may guide future antitrypanosomal drug-discovery efforts. We are hopeful that this new version will be a useful tool for the computational- and medicinal-chemist community.

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

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

Computational modeling of subcellular transport and signaling

Numerous signaling processes in the cell are controlled in microdomains that are defined by cellular structures ranging from nm to μm in size. Recent improvements in microscopy enable the resolution and reconstruction of these micro domains, while new computational methods provide the means to elucidate their functional roles. Collectively these tools allow for a biophysical understanding of the cellular environment and its pathological progression in disease. Here we review recent advancements in microscopy, and subcellular modeling on the basis of reconstructed geometries, with a special focus on signaling microdomains that are important for the excitation contraction coupling in cardiac myocytes.

Phosphorylated phospholamban stabilizes a compact conformation of the cardiac calcium-ATPase

The sarcoendoplasmic reticulum calcium ATPase (SERCA) plays a key role in cardiac calcium handling and is considered a high-value target for the treatment of heart failure. SERCA undergoes conformational changes as it harnesses the chemical energy of ATP for active transport. X-ray crystallography has provided insight into SERCA structural substates, but it is not known how well these static snapshots describe in vivo conformational dynamics. The goals of this work were to quantify the direction and magnitude of SERCA motions as the pump performs work in live cardiac myocytes, and to identify structural determinants of SERCA regulation by phospholamban. We measured intramolecular fluorescence resonance energy transfer (FRET) between fluorescent proteins fused to SERCA cytoplasmic domains. We detected four discrete structural substates for SERCA expressed in cardiac muscle cells. The relative populations of these discrete states oscillated with electrical pacing. Low FRET states were most populated in low Ca (diastole), and were indicative of an open, disordered structure for SERCA in the E2 (Ca-free) enzymatic substate. High FRET states increased with Ca (systole), suggesting rigidly closed conformations for the E1 (Ca-bound) enzymatic substates. Notably, a special compact E1 state was observed after treatment with β-adrenergic agonist or with coexpression of phosphomimetic mutants of phospholamban. The data suggest that SERCA calcium binding induces the pump to undergo a transition from an open, dynamic conformation to a closed, ordered structure. Phosphorylated phospholamban stabilizes a unique conformation of SERCA that is characterized by a compact architecture.

New crystal structures of PII-type ATPases: excitement continues

P-type ATPases are ATP-powered ion pumps, classified into five subfamilies (PI-PV). Of these, PII-type ATPases, including Ca2+-ATPase, Na+,K+-ATPase and gastric H+,K+-ATPase, among others, have been the most intensively studied. Best understood structurally and biochemically is Ca2+-ATPase from sarcoplasmic reticulum of fast twitch skeletal muscle (sarco(endo)plasmic reticulum Ca2+-ATPase 1a, SERCA1a). Since publication of the first crystal structure in 2000, it has continuously been a source of excitement, as crystal structures for new reaction intermediates always show large structural changes. Crystal structures now exist for most of the reaction intermediates, almost covering the entire reaction cycle. This year the crystal structure of a missing link, the E1·Mg2+ state, finally appeared, bringing another surprise: bound sarcolipin (SLN). The current status of two other important PII-type ATPases, Na+,K+-ATPase and H+,K+-ATPase, is also briefly described.

Finite Element Estimation of Protein-Ligand Association Rates with Post-Encounter Effects: Applications to Calcium binding in Troponin C and SERCA

We introduce a computational pipeline and suite of software tools for the approximation of diffusion-limited binding based on a recently developed theoretical framework. Our approach handles molecular geometries generated from high-resolution structural data and can account for active sites buried within the protein or behind gating mechanisms. Using tools from the FEniCS library and the APBS solver, we implement a numerical code for our method and study two Ca(2+)-binding proteins: Troponin C and the Sarcoplasmic Reticulum Ca(2+) ATPase (SERCA). We find that a combination of diffusional encounter and internal 'buried channel' descriptions provide superior descriptions of association rates, improving estimates by orders of magnitude.