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

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.

Isolation of the Sarcoplasmic Reticulum Ca2+-ATPase from Rabbit Fast-Twitch Muscle

The sarcoendoplasmic reticulum Ca2+-ATPase (SERCA) is a membrane protein that is destabilized during purification in the absence of calcium ions. The disaccharide trehalose is a protein stabilizer that accumulates in the yeast cytoplasm when under stress. In the present work, SERCA was purified by including trehalose in the purification protocol. The purified SERCA showed high protein purity (~95%) and ATPase activity. ATP hydrolysis was dependent on the presence of Ca2+ and the enzyme kinetics showed a hyperbolic dependence on ATP (Km = 12.16 ± 2.25 μM ATP). FITC labeling showed the integrity of the ATP-binding site and the identity of the isolated enzyme as a P-type ATPase. Circular dichroism (CD) spectral changes at a wavelength of 225 nm were observed upon titration with ATP, indicating α-helical rearrangements in the nucleotide-binding domain (N-domain), which correlated with ATP affinity (Km). The presence of Ca2+ did not affect FITC labeling or the ATP-mediated structural changes at the N-domain. The use of trehalose in the SERCA purification protocol stabilized the enzyme. The isolated SERCA appears to be suitable for structural and ligand binding studies, e.g., for testing newly designed or natural inhibitors. The use of trehalose is recommended for the isolation of unstable enzymes.

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.

Structural Basis for the Function of the C-Terminal Proton Release Pathway in the Calcium Pump

The calcium pump (sarco/endoplasmic reticulum Ca²⁺-ATPase, SERCA) plays a major role in calcium homeostasis in muscle cells by clearing cytosolic Ca²⁺ during muscle relaxation. Active Ca²⁺ transport by SERCA involves the structural transition from a low-Ca²⁺ affinity E2 state toward a high-Ca²⁺ affinity E1 state of the pump. This structural transition is accompanied by the countertransport of protons to stabilize the negative charge and maintain the structural integrity of the transport sites and partially compensate for the positive charges of the two Ca²⁺ ions passing through the membrane. X-ray crystallography studies have suggested that a hydrated pore located at the C-terminal domain of SERCA serves as a conduit for proton countertransport, but the existence and function of this pathway have not yet been fully characterized. We used atomistic simulations to demonstrate that in the protonated E2 state and the absence of initially bound water molecules, the C-terminal pore becomes hydrated in the nanosecond timescale. Hydration of the C-terminal pore is accompanied by the formation of water wires that connect the transport sites with the cytosol. Water wires are known as ubiquitous proton-transport devices in biological systems, thus supporting the notion that the C-terminal domain serves as a conduit for proton release. Additional simulations showed that the release of a single proton from the transport sites induces bending of transmembrane helix M5 and the interaction between residues Arg762 and Ser915. These structural changes create a physical barrier against full hydration of the pore and prevent the formation of hydrogen-bonded water wires once proton transport has occurred through this pore. Together, these findings support the notion that the C-terminal proton release pathway is a functional element of SERCA and also provide a mechanistic model for its operation in the catalytic cycle of the pump.

Deciphering the Mechanism of Inhibition of SERCA1a by Sarcolipin Using Molecular Simulations

SERCA1a is an ATPase calcium pump that transports Ca²⁺ from the cytoplasm to the sarco/endoplasmic reticulum lumen. Sarcolipin (SLN), a transmembrane peptide, regulates the activity of SERCA1a by decreasing its Ca²⁺ transport rate, but its mechanism of action is still not well-understood. To decipher this mechanism, we have performed normal mode analysis in the all-atom model, with the SERCA1a-SLN complex, or the isolated SERCA1a, embedded in an explicit membrane. The comparison of the results allowed us to provide an explanation at the atomic level for the action of SLN that is in good agreement with experimental observations. In our analyses, the presence of SLN locally perturbs the TM6 transmembrane helix and as a consequence modifies the position of D800, one of the key metal-chelating residues. Additionally, it reduces the flexibility of the gating residues, V304, and E309 in TM4, at the entrance of the Ca²⁺ binding sites, which would decrease the affinity for Ca²⁺. Unexpectedly, SLN has also an effect on the ATP binding site more than 35 Å away, due to the straightening of TM5, a long helix considered as the spine of the protein. The straightening of TM5 modifies the structure of the P-N linker that sits above it, and which comprises the ³⁵¹DKTG³⁵⁴ conserved motif, resulting in an increase of the distance between ATP and the phosphorylation site. As a consequence, the turn-over rate could be affected. All this gives SERCA1a the propensity to go toward a Ca²⁺ low-affinity E2-like state in the presence of SLN and toward a Ca²⁺ high-affinity E1-like state in the absence of SLN. In addition to a general mechanism of inhibition of SERCA1a regulatory peptides, this study also provides an insight into the conformational transition between the E2 and E1 states.

Effects of amino acid modifications on the permeability of the pentameric sarcolipin channel

Sarcolipin (SLN) is an important transmembrane (TM) protein encoded by long noncoding RNA. SLN is expressed in the sarcoplasmic reticulum and regulates cardiac and skeletal muscle contractions. SLN forms a pentameric hydrophobic ligand-gated ion channel. The protonation of Glu7 (protonated SLN, pSLN) and mutation of Thr18 to Ala18 (T18A) have been reported to exert a significant influence on the permeability of the channel. In this study, the altered permeability of both the pSLN and T18A pentameric channels was simulated. Combined with molecular dynamics simulation, the free-energy landscape for single ions, computational electrophysiology, diffusion coefficient, and pore geometrical characteristic analyses were performed to further understand the properties of amino acid modifications in the SLN pentameric channel. The results suggest that both the pSLN and T18A pentameric channels form stable hydrophobic ligand-gated channels. The TM voltage has a positive effect on the permeability of water molecules and ions. By using pSLN and T18A, our study provides helpful information on the pore-forming mechanism of SLN and furthers our understanding of the regulatory mechanisms underlying the permeation of ions and water molecules in the pentameric SLN channel.

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.

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

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

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.

Ca2+-based allosteric switches and shape shifting in RGLG1 VWA domain

RGLG1 is an E3 ubiquitin ligase in Arabidopsis thaliana that participates in ABA signaling and regulates apical dominance. Here, we present crystal structures of RGLG1 VWA domain, revealing two novel calcium ions binding sites (NCBS1 and NCBS2). Furthermore, the structures with guided mutagenesis in NCBS1 prove that Ca²⁺ ions play important roles in controlling conformational change of VWA, which is stabilized in open state with Ca²⁺ bound and converted to closed state after Ca²⁺ removal. This allosteric regulation mechanism is distinct from the ever reported one involving the C-terminal helix in integrin α and β I domains. The mutation of a key residue in NCBS2 do not abolish its Ca²⁺-binding potential, with no conformational change. MD simulations reveals that open state of RGLG1 VWA has higher ligand affinity than its closed state, consisting with integrin. Structural comparison of ion-free-MIDAS with Mg²⁺-MIDAS reveals that Mg²⁺ binding to MIDAS does not induce conformational change. With acquisition of first structure of plant VWA domain in both open state and closed state, we carefully analyze the conformational change and propose a totally new paradigm for its transition of open-closed states, which will be of great value for guiding future researches on VWA proteins and their important biological significance.

ANS Interacts with the Ca2+-ATPase Nucleotide Binding Site

The binding of 8-anilino-1-naphthalene sulfonate (ANS) to the nucleotide binding domain (N-domain) of the sarcoplasmic reticulum Ca²⁺-ATPase (SERCA) was studied. Molecular docking predicted two ANS binding modes (BMI and BMII) in the nucleotide binding site. The molecular interaction was confirmed as the fluorescence intensity of ANS was dramatically increased when in the presence of an engineered recombinant N-domain. Molecular dynamics simulation showed BMI (which occupies the ATP binding site) as the mode that is stable in solution. The above was confirmed by the absence of ANS fluorescence in the presence of a fluorescein isothiocyanate (FITC)-labeled N-domain. Further, the labeling of the N-domain with FITC was hindered by the presence of ANS, i.e., ANS was bound to the ATP binding site. Importantly, ANS displayed a higher affinity than ATP. In addition, ANS binding led to quenching the N-domain intrinsic fluorescence displaying a FRET pattern, which suggested the existence of a Trp-ANS FRET couple. Nonetheless, the chemical modification of the sole Trp residue with N-bromosuccinimide (NBS) discarded the existence of FRET and instead indicated structural rearrangements in the nucleotide binding site during ANS binding. Finally, Ca²⁺-ATPase kinetics in the presence of ANS showed a partial mixed-type inhibition. The Dixon plot showed the ANS-Ca²⁺-ATPase complex as catalytically active, hence supporting the existence of a functional dimeric Ca²⁺-ATPase in sarcoplasmic reticulum vesicles. ANS may be used as a molecular platform for the development of more effective inhibitors of Ca²⁺-ATPase and appears to be a new fluorescent probe for the nucleotide binding site. Graphical Abstract Molecular docking of ANS to the nucleotide binding site of Ca²⁺-ATPase. ANS fluorescence increase reveals molecular interaction.

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.

Probing the effects of nonannular lipid binding on the stability of the calcium pump SERCA

The calcium pump SERCA is a transmembrane protein that is critical for calcium transport in cells. SERCA resides in an environment made up largely by the lipid bilayer, so lipids play a central role on its stability and function. Studies have provided insights into the effects of annular and bulk lipids on SERCA activation, but the role of a nonannular lipid site in the E2 intermediate state remains elusive. Here, we have performed microsecond molecular dynamics simulations to probe the effects of nonannular lipid binding on the stability and structural dynamics of the E2 state of SERCA. We found that the structural integrity and stability of the E2 state is independent of nonannular lipid binding, and that occupancy of a lipid molecule at this site does not modulate destabilization of the E2 state, a step required to initiate the transition toward the competent E1 state. We also found that binding of the nonannular lipid does not induce direct allosteric control of the intrinsic functional dynamics the E2 state. We conclude that nonannular lipid binding is not necessary for the stability of the E2 state, but we speculate that it becomes functionally significant during the E2-to-E1 transition of the pump.

Proton Countertransport and Coupled Gating in the Sarcoplasmic Reticulum Calcium Pump

The calcium pump of the sarcoplasmic reticulum (SERCA) is an ATP-driven active transporter of Ca²⁺ ions that functions via an "alternating-access" cycle mechanism. In each cycle, SERCA transports two Ca²⁺ ions toward the lumen of the sarcoplasmic reticulum and two to three protons to the cytoplasm. How the latter conformational transition is coupled to cytoplasmic release of protons remains poorly understood. The present computational study shows how the mechanism of proton countertransport is coupled to the alternating access gating process in SERCA. Molecular dynamics simulation trajectories are generated starting from a series of configurations taken along the E2 to E1 transition pathway determined by the string method with swarms-of-trajectories. Simulations of different protonation configurations at the binding sites reveal how deprotonation events affect the opening of the cytoplasmic gate. The results show that there is a strong coupling between the chronological order of deprotonation, the entry of water molecules into the TM region, and the opening of the cytoplasmic gate. Deprotonation of E309 and E771 is sequential with E309 being the first to lose the proton. The deprotonation promotes the opening of the cytoplasmic gate but leads to a productive gating transition only if it occurs after the transmembrane domain has reached an intermediate conformation. Deprotonation of E309 and E771 is unproductive when it occurs too early because it causes the re-opening of the luminal gate.

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.

The Ca2+-ATPase pump facilitates bidirectional proton transport across the sarco/endoplasmic reticulum

Ca²⁺ transport across the sarco/endoplasmic reticulum (SR) plays an essential role in intracellular Ca²⁺ homeostasis, signalling, cell differentiation and muscle contractility. During SR Ca²⁺ uptake and release, proton fluxes are required to balance the charge deficit generated by the exchange of Ca²⁺ and other ions across the SR. During Ca²⁺ uptake by the SR Ca²⁺-ATPase (SERCA), two protons are countertransported from the SR lumen to the cytosol, thus partially compensating for the charge moved by Ca²⁺ transport. Studies have shown that protons are also transported from the cytosol to the lumen during Ca²⁺ release, but a transporter that facilitates proton transport into the SR lumen has not been described. In this article we propose that SERCA forms pores that facilitate bidirectional proton transport across the SR. We describe the location and structure of water-filled pores in SERCA that form cytosolic and luminal pathways for protons to cross the SR membrane. Based on this structural information, we suggest mechanistic models for proton translocation to the cytosol during active Ca²⁺ transport, and into the SR lumen during SERCA inhibition by endogenous regulatory proteins. Finally, we discuss the physiological consequences of SERCA-mediated bidirectional proton transport across the SR membrane of muscle and non-muscle cells.

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

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

Conformational Transitions and Alternating-Access Mechanism in the Sarcoplasmic Reticulum Calcium Pump

Ion pumps are integral membrane proteins responsible for transporting ions against concentration gradients across biological membranes. Sarco/endoplasmic reticulum Ca2+-ATPase (SERCA), a member of the P-type ATPases family, transports two calcium ions per hydrolyzed ATP molecule via an "alternating-access" mechanism. High-resolution crystallographic structures provide invaluable insight on the structural mechanism of the ion pumping process. However, to understand the molecular details of how ATP hydrolysis is coupled to calcium transport, it is necessary to gain knowledge about the conformational transition pathways connecting the crystallographically resolved conformations. Large-scale transitions in SERCA occur at time-scales beyond the current reach of unbiased molecular dynamics simulations. Here, we overcome this challenge by employing the string method, which represents a transition pathway as a chainofstates linking two conformational endpoints. Using a multiscale methodology, we have determined all-atom transition pathways for three main conformational transitions responsible for the alternating-access mechanism. The present pathways provide a clear chronology and ordering of the key events underlying the active transport of calcium ions by SERCA. Important conclusions are that the conformational transition that leads to occlusion with bound ATP and calcium is highly concerted and cooperative, the phosphorylation of Asp351 causes areorganization of the cytoplasmic domains that subsequently drives the opening of the luminal gate, and thereclosing of luminal gate induces a shift in the cytoplasmic domains that subsequently enables the dephosphorylation of Asp351-P. Formation of transient residue-residue contacts along the conformational transitions predicted by the computations provide an experimental route to test the general validity of the computational pathways.

Sarcolipin Promotes Uncoupling of the SERCA Ca2+ Pump by Inducing a Structural Rearrangement in the Energy-Transduction Domain

We have performed microsecond (μs) molecular dynamics simulation (MDS) to identify structural mechanisms for sarcolipin (SLN) uncoupling of Ca2+ transport from ATP hydrolysis for the sarcoplasmic reticulum Ca2+-ATPase (SERCA). SLN regulates muscle metabolism and energy expenditure to provide resistance against diet-induced obesity and extreme cold. MDS demonstrated that the cytosolic domain of SLN induces a salt bridge-mediated structural rearrangement in the energy-transduction domain of SERCA. We propose that this structural change uncouples SERCA by perturbing Ca2+ occlusion at residue E309 in transport site II, thus facilitating Ca2+ backflux to the cytosol. Our results have important implications for designing muscle-based therapies for human obesity.

Optimization of Analytical Potentials for Coarse-Grained Biopolymer Models

The increasing trend in the recent literature on coarse grained (CG) models testifies their impact in the study of complex systems. However, the CG model landscape is variegated: even considering a given resolution level, the force fields are very heterogeneous and optimized with very different parametrization procedures. Along the road for standardization of CG models for biopolymers, here we describe a strategy to aid building and optimization of statistics based analytical force fields and its implementation in the software package AsParaGS (Assisted Parameterization platform for coarse Grained modelS). Our method is based on the use and optimization of analytical potentials, optimized by targeting internal variables statistical distributions by means of the combination of different algorithms (i.e., relative entropy driven stochastic exploration of the parameter space and iterative Boltzmann inversion). This allows designing a custom model that endows the force field terms with a physically sound meaning. Furthermore, the level of transferability and accuracy can be tuned through the choice of statistical data set composition. The method-illustrated by means of applications to helical polypeptides-also involves the analysis of two and three variable distributions, and allows handling issues related to the FF term correlations. AsParaGS is interfaced with general-purpose molecular dynamics codes and currently implements the "minimalist" subclass of CG models (i.e., one bead per amino acid, Cα based). Extensions to nucleic acids and different levels of coarse graining are in the course.

Engineering a Prototypic P-type ATPase Listeria monocytogenes Ca(2+)-ATPase 1 for Single-Molecule FRET Studies

Approximately 30% of the ATP generated in the living cell is utilized by P-type ATPase primary active transporters to generate and maintain electrochemical gradients across biological membranes. P-type ATPases undergo large conformational changes during their functional cycle to couple ATP hydrolysis in the cytoplasmic domains to ion transport across the membrane. The Ca(2+)-ATPase from Listeria monocytogenes, LMCA1, was found to be a suitable model of P-type ATPases and was engineered to facilitate single-molecule FRET studies of transport-related structural changes. Mutational analyses of the endogenous cysteine residues in LMCA1 were performed to reduce background labeling without compromising activity. Pairs of cysteines were introduced into the optimized low-reactivity background, and labeled with maleimide derivatives of Cy3 and Cy5 resulting in site-specifically double-labeled protein with moderate activity. Ensemble and confocal single-molecule FRET studies revealed changes in FRET distribution related to structural changes during the transport cycle, consistent with those observed by X-ray crystallography for the sarco/endoplasmic reticulum Ca(2+) ATPase (SERCA). Notably, the cytosolic headpiece of LMCA1 was found to be distinctly more compact in the E1 state than in the E2 state. Thus, the established experimental system should allow future real-time FRET studies of the structural dynamics of LMCA1 as a representative P-type ATPase.

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.

Transmembrane dynamics of the Thr-5 phosphorylated sarcolipin pentameric channel

Sarcolipin (SLN), an important membrane protein expressed in the sarcoplasmic reticulum (SR), regulates muscle contractions in cardiac and skeletal muscle. The phosphorylation at amino acid Thr5 of the SLN protein modulates the amount of Ca(2+) that passes through the SR. Using molecular dynamics simulation, we evaluated the phosphorylation at Thr5 of pentameric SLN (phospho-SLN) channel's energy barrier and pore characteristics by calculating the potential of mean force (PMF) along the channel pore and determining the diffusion coefficient. The results indicate that pentameric phospho-SLN promotes penetration of monovalent and divalent ions through the channel. The analysis of PMF, pore radius and diffusion coefficient indicates that Leu21 is the hydrophobic gate of the pentameric SLN channel. In the channel, water molecules near the Leu21 pore demonstrated a clear hydrated-dehydrated transition; however, the mutation of Leu21 to an Alanine (L21A) destroyed the hydrated-dehydrated transitions. These water-dynamic behaviors and PMF confirm that Leu21 is the key residue that regulates the ion permeability of the pentameric SLN channel. These results provide the structural-basis insights and molecular-dynamic information that are needed to understand the regulatory mechanisms of ion permeability in the pentameric SLN channel.

Pathogenic mutation R959W alters recognition dynamics of dysferlin inner DysF domain

Dysferlin, a 220 kD protein, plays a major role in regulating plasma membrane repair in muscle cells. Mutations in the dysferlin inner DysF domain are known to cause different types of muscular dystrophy, including limb-girdle muscular dystrophy type 2B (LGMD2B) and Miyoshi myopathy (MM). Replacement of arginine in position 959 by tryptophan has been frequently associated with both LGMD2B and MM, but the molecular mechanisms by which this mutation alters dysferlin function remain unknown. In this study, we have used protein binding site predictions and microsecond molecular dynamics (MD) simulations to determine the effects pathogenic mutation R959W on the structural dynamics of dysferlin inner DysF domain. Analysis of 2 μs long MD trajectories revealed that mutation R959W does not induce local destabilization, unfolding or misfolding of the domain. We used a binding site predictor to discover a protein-binding site (residues T958-I966 and E1031-H1037) that resembles pincers in shape. Cartesian principal component analysis and interresidue distance distributions of the wild-type domain showed that the predicted protein-binding site undergoes a pincer motion, and populates two structural states, open and closed. We found that mutation R959W inhibits the pincer motion of the protein-binding site and completely shifts the equilibrium toward the open state. These differences in the structural dynamics of the predicted binding site suggest that mutation R959W alters recognition dynamics of the inner DysF domain. Based on these findings and on previous experimental studies, we propose a novel role for the inner DysF domain in muscle membrane repair through recruitment of dysferlin to plasma membrane. In conclusion, these findings have important implications for our understanding of the structural aspects of muscular dystrophies in atomic-level resolution.

Atomic-level mechanisms for phospholamban regulation of the calcium pump

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

Biomembranes in atomistic and coarse-grained simulations

The architecture of biological membranes is tightly coupled to the localization, organization, and function of membrane proteins. The organelle-specific distribution of lipids allows for the formation of functional microdomains (also called rafts) that facilitate the segregation and aggregation of membrane proteins and thus shape their function. Molecular dynamics simulations enable to directly access the formation, structure, and dynamics of membrane microdomains at the molecular scale and the specific interactions among lipids and proteins on timescales from picoseconds to microseconds. This review focuses on the latest developments of biomembrane force fields for both atomistic and coarse-grained molecular dynamics (MD) simulations, and the different levels of coarsening of biomolecular structures. It also briefly introduces scale-bridging methods applicable to biomembrane studies, and highlights selected recent applications.

Sarcolipin and phospholamban inhibit the calcium pump by populating a similar metal ion-free intermediate state

We have performed microsecond molecular dynamics (MD) simulations and protein pKa calculations of the muscle calcium pump (sarcoplasmic reticulum Ca(2+)-ATPase, SERCA) in complex with sarcolipin (SLN) to determine the mechanism by which SLN inhibits SERCA. SLN and its close analog phospholamban (PLN) are membrane proteins that regulate SERCA by inhibiting Ca(2+) transport in skeletal and cardiac muscle. Although SLN and PLB binding to SERCA have different functional outcomes on the coupling efficiency of SERCA, both proteins decrease the apparent Ca(2+) affinity of the pump, suggesting that SLN and PLB inhibit SERCA by using a similar mechanism. Recently, MD simulations showed that PLB inhibits SERCA by populating a metal ion-free, partially-protonated E1 state of the pump, E1· [Formula: see text] . X-ray crystallography studies at 40-80 mM Mg(2+) have proposed that SLN-bound SERCA populates E1·Mg(2+), an intermediate with Mg(2+) bound near transport site I. To test this proposed mode of SLN regulation, we performed a 0.5-μs MD simulation of E1·Mg(2+)-SLN in a solution containing 100 mM K(+) and 3 mM Mg(2+), with calculation of domain dynamics in the cytosolic headpiece and side-chain ionization and occupancy in the transport sites. We found that SLN increases the distance between residues E771 and D800, thereby rendering E1·Mg(2+) incapable of producing a competent Ca(2+) transport site I. Following removal of Mg(2+,) a 2-μs MD simulation of Mg(2+)-free SERCA-SLN showed that Mg(2+) does not re-bind to the transport sites, indicating that SERCA-SLN does not populate E1·Mg(2+) at physiological conditions. Instead, protein pKa calculations indicate that SLN stabilizes a metal ion-free SERCA state (E1· [Formula: see text] ) protonated at residue E771, but ionized at E309 and D800. We conclude that both SLN and PLB inhibit SERCA by populating a similar metal ion-free intermediate state. We propose that (i) this partially-protonated intermediate serves as the consensus mechanism for SERCA inhibition by other members of the SERCA regulatory subunit family including myoregulin and sarcolamban, and (ii) this consensus mechanism is utilized to regulate Ca(2+) transport in skeletal and cardiac muscle, with important implications for therapeutic approaches to muscle dystrophy and heart failure.

The N Terminus of Sarcolipin Plays an Important Role in Uncoupling Sarco-endoplasmic Reticulum Ca2+-ATPase (SERCA) ATP Hydrolysis from Ca2+ Transport

The sarcoendoplasmic reticulum Ca(2+)-ATPase (SERCA) is responsible for intracellular Ca(2+) homeostasis. SERCA activity in muscle can be regulated by phospholamban (PLB), an affinity modulator, and sarcolipin (SLN), an uncoupler. Although PLB gets dislodged from Ca(2+)-bound SERCA, SLN continues to bind SERCA throughout its kinetic cycle and promotes uncoupling of Ca(2+) transport from ATP hydrolysis. To determine the structural regions of SLN that mediate uncoupling of SERCA, we employed mutagenesis and generated chimeras of PLB and SLN. In this study we demonstrate that deletion of SLN N-terminal residues (2)ERSTQ leads to loss of the uncoupling function even though the truncated peptide can target and constitutively bind SERCA. Furthermore, molecular dynamics simulations of SLN and SERCA interaction showed a rearrangement of SERCA residues that is altered when the SLN N terminus is deleted. Interestingly, transfer of the PLB cytosolic domain to the SLN transmembrane (TM) and luminal tail causes the chimeric protein to lose SLN-like function. Further introduction of the PLB TM region into this chimera resulted in conversion to full PLB-like function. We also found that swapping PLB N and C termini with those from SLN caused the resulting chimera to acquire SLN-like function. Swapping the C terminus alone was not sufficient for this conversion. These results suggest that domains can be switched between SLN and PLB without losing the ability to regulate SERCA activity; however, the resulting chimeras acquire functions different from the parent molecules. Importantly, our studies highlight that the N termini of SLN and PLB influence their respective unique functions.

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.