Relationship between Ca2+-affinity and shielding of bulk water in the Ca2+-pump from molecular dynamics simulations

Crystal structures of Na+ ,K+ -ATPase reveal the mechanism that converts the K+ -bound form to Na+ -bound form and opens and closes the cytoplasmic gate

Na+ ,K+ -ATPase (NKA) plays a pivotal role in establishing electrochemical gradients for Na+ and K+ across the cell membrane by alternating between the E1 (showing high affinity for Na+ and low affinity for K+ ) and E2 (low affinity to Na+ and high affinity to K+ ) forms. Presented here are two crystal structures of NKA in E1·Mg2+ and E1·3Na+ states at 2.9 and 2.8 Å resolution, respectively. These two E1 structures fill a gap in our description of the NKA reaction cycle based on the atomic structures. We describe how NKA converts the K+ -bound E2·2K+ form to an E1 (E1·Mg2+ ) form, which allows high-affinity Na+ binding, eventually closing the cytoplasmic gate (in E1 ~ P·ADP·3Na+ ) after binding three Na+ , while keeping the extracellular ion pathway sealed. We now understand previously unknown functional roles for several parts of NKA and that NKA uses even the lipid bilayer for gating the ion pathway.

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.

GENESIS 1.1: A hybrid-parallel molecular dynamics simulator with enhanced sampling algorithms on multiple computational platforms

GENeralized-Ensemble SImulation System (GENESIS) is a software package for molecular dynamics (MD) simulation of biological systems. It is designed to extend limitations in system size and accessible time scale by adopting highly parallelized schemes and enhanced conformational sampling algorithms. In this new version, GENESIS 1.1, new functions and advanced algorithms have been added. The all-atom and coarse-grained potential energy functions used in AMBER and GROMACS packages now become available in addition to CHARMM energy functions. The performance of MD simulations has been greatly improved by further optimization, multiple time-step integration, and hybrid (CPU + GPU) computing. The string method and replica-exchange umbrella sampling with flexible collective variable choice are used for finding the minimum free-energy pathway and obtaining free-energy profiles for conformational changes of a macromolecule. These new features increase the usefulness and power of GENESIS for modeling and simulation in biological research. © 2017 Wiley Periodicals, Inc.

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.

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.

CHARMM Force-Fields with Modified Polyphosphate Parameters Allow Stable Simulation of the ATP-Bound Structure of Ca(2+)-ATPase

Adenosine triphosphate (ATP) is an indispensable energy source in cells. In a wide variety of biological phenomena like glycolysis, muscle contraction/relaxation, and active ion transport, chemical energy released from ATP hydrolysis is converted to mechanical forces to bring about large-scale conformational changes in proteins. Investigation of structure-function relationships in these proteins by molecular dynamics (MD) simulations requires modeling of ATP in solution and ATP bound to proteins with accurate force-field parameters. In this study, we derived new force-field parameters for the triphosphate moiety of ATP based on the high-precision quantum calculations of methyl triphosphate. We tested our new parameters on membrane-embedded sarcoplasmic reticulum Ca(2+)-ATPase and four soluble proteins. The ATP-bound structure of Ca(2+)-ATPase remains stable during MD simulations, contrary to the outcome in shorter simulations using original parameters. Similar results were obtained with the four ATP-bound soluble proteins. The new force-field parameters were also tested by investigating the range of conformations sampled during replica-exchange MD simulations of ATP in explicit water. Modified parameters allowed a much wider range of conformational sampling compared with the bias toward extended forms with original parameters. A diverse range of structures agrees with the broad distribution of ATP conformations in proteins deposited in the Protein Data Bank. These simulations suggest that the modified parameters will be useful in studies of ATP in solution and of the many ATP-utilizing proteins.

Mechanism of the αβ conformational change in F1-ATPase after ATP hydrolysis: free-energy simulations

One of the motive forces for F1-ATPase rotation is the conformational change of the catalytically active β subunit due to closing and opening motions caused by ATP binding and hydrolysis, respectively. The closing motion is accomplished in two steps: the hydrogen-bond network around ATP changes and then the entire structure changes via B-helix sliding, as shown in our previous study. Here, we investigated the opening motion induced by ATP hydrolysis using all-atom free-energy simulations, combining the nudged elastic band method and umbrella sampling molecular-dynamics simulations. Because hydrolysis requires residues in the α subunit, the simulations were performed with the αβ dimer. The results indicate that the large-scale opening motion is also achieved by the B-helix sliding (in the reverse direction). However, the sliding mechanism is different from that of ATP binding because sliding is triggered by separation of the hydrolysis products ADP and Pi. We also addressed several important issues: 1), the timing of the product Pi release; 2), the unresolved half-closed β structure; and 3), the ADP release mechanism. These issues are fundamental for motor function; thus, the rotational mechanism of the entire F1-ATPase is also elucidated through this αβ study. During the conformational change, conserved residues among the ATPase proteins play important roles, suggesting that the obtained mechanism may be shared with other ATPase proteins. When combined with our previous studies, these results provide a comprehensive view of the β-subunit conformational change that drives the ATPase.

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.

K+ congeners that do not compromise Na+ activation of the Na+,K+-ATPase: hydration of the ion binding cavity likely controls ion selectivity

The Na(+),K(+)-ATPase is essential for ionic homeostasis in animal cells. The dephosphoenzyme contains Na(+) selective inward facing sites, whereas the phosphoenzyme contains K(+) selective outward facing sites. Under normal physiological conditions, K(+) inhibits cytoplasmic Na(+) activation of the enzyme. Acetamidinium (Acet(+)) and formamidinium (Form(+)) have been shown to permeate the pump through the outward facing sites. Here, we show that these cations, unlike K(+), are unable to enter the inward facing sites in the dephosphorylated enzyme. Consistently, the organic cations exhibited little to no antagonism to cytoplasmic Na(+) activation. Na(+),K(+)-ATPase structures revealed a previously undescribed rotamer transition of the hydroxymethyl side chain of the absolutely conserved Thr(772) of the α-subunit. The side chain contributes its hydroxyl to Na(+) in site I in the E1 form and rotates to contribute its methyl group toward K(+) in the E2 form. Molecular dynamics simulations to the E1·AlF4 (-)·ADP·3Na(+) structure indicated that 1) bound organic cations differentially distorted the ion binding sites, 2) the hydroxymethyl of Thr(772) rotates to stabilize bound Form(+) through water molecules, and 3) the rotamer transition is mediated by water traffic into the ion binding cavity. Accordingly, dehydration induced by osmotic stress enhanced the interaction of the congeners with the outward facing sites and profoundly modified the organization of membrane domains of the α-subunit. These results assign a catalytic role for water in pump function, and shed light on a backbone-independent but a conformation-dependent switch between H-bond and dispersion contact as part of the catalytic mechanism of the Na(+),K(+)-ATPase.

Investigation of the Josephin Domain protein-protein interaction by molecular dynamics

Spinocerebellar ataxia (SCA) 3, the most common form of SCA, is a neurodegenerative rare disease characterized by polyglutamine tract expansion and self-assembly of Ataxin3 (At3) misfolded proteins into highly organized fibrillar aggregates. The At3 N-terminal Josephin Domain (JD) has been suggested as being responsible for mediating the initial phase of the At3 double-step fibrillogenesis. Several issues concerning the residues involved in the JD's aggregation and, more generally, the JD clumping mechanism have not been clarified yet. In this paper we present an investigation focusing on the JD protein-protein interaction by means of molecular modeling. Our results suggest possible aminoacids involved in JD contact together with local and non-local effects following JD dimerization. Surprisingly, JD conformational changes following the binding may involve ubiquitin binding sites and hairpin region even though they do not pertain to the JD interaction surfaces. Moreover, the JD binding event has been found to alter the hairpin open-like conformation toward a closed-like arrangement over the simulated timescale. Finally, our results suggest that the JD aggregation might be a multi-step process, with an initial fast JD-JD binding mainly driven by Arg101, followed by slower structural global rearrangements involving the exposure to the solvent of Leu84-Trp87, which might play a role in a second step of JD aggregation.

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.

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.

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

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

Molecular dynamics simulations reveal proton transfer pathways in cytochrome C-dependent nitric oxide reductase

Nitric oxide reductases (NORs) are membrane proteins that catalyze the reduction of nitric oxide (NO) to nitrous oxide (N₂O), which is a critical step of the nitrate respiration process in denitrifying bacteria. Using the recently determined first crystal structure of the cytochrome c-dependent NOR (cNOR) [Hino T, Matsumoto Y, Nagano S, Sugimoto H, Fukumori Y, et al. (2010) Structural basis of biological N2O generation by bacterial nitric oxide reductase. Science 330: 1666–70.], we performed extensive all-atom molecular dynamics (MD) simulations of cNOR within an explicit membrane/solvent environment to fully characterize water distribution and dynamics as well as hydrogen-bonded networks inside the protein, yielding the atomic details of functionally important proton channels. Simulations reveal two possible proton transfer pathways leading from the periplasm to the active site, while no pathways from the cytoplasmic side were found, consistently with the experimental observations that cNOR is not a proton pump. One of the pathways, which was newly identified in the MD simulation, is blocked in the crystal structure and requires small structural rearrangements to allow for water channel formation. That pathway is equivalent to the functional periplasmic cavity postulated in cbb₃ oxidase, which illustrates that the two enzymes share some elements of the proton transfer mechanisms and confirms a close evolutionary relation between NORs and C-type oxidases. Several mechanisms of the critical proton transfer steps near the catalytic center are proposed.

Denitrification is an anaerobic process performed by several bacteria as an alternative to aerobic respiration. A key intermediate step is catalyzed by the nitric oxide reductase (NOR) enzyme, which is situated in the cytoplasmic membrane. Proton delivery to the catalytic site inside NOR is an important part of its functioning. In this work we use molecular dynamics simulations to describe water distribution and to identify proton transfer pathways in cNOR. Our results reveal two channels from the periplasmic side of the membrane and none from the cytoplasmic side, indicating that cNOR is not a proton pump. It is our hope that these results will provide a basis for further experimental and computational studies aimed to understand details of the NOR mechanism. Furthermore, this work sheds light on the molecular evolution of respiratory enzymes.

Identifying the sarco(endo)plasmic reticulum Ca2+ ATPase (SERCA) as a potential target for hypericin--a theoretical study

The exact cellular target for the potent anti-cancer agent hypericin has not yet been determined; this thus encourages the application of computational chemistry tools to be employed in order to provide insights that can be employed in further drug development studies. In the present study computational docking and molecular dynamics simulations are applied to investigate possible interactions between hypericin and the Ca(2+) pump SERCA as proposed in the literature. Hypericin was found to bind strongly both in pockets within the transmembrane region and in the cytosolic region of the protein, although the two studied isoforms of SERCA differ slightly in their preferred binding sites. The calculated binding energies for hypericin in the four investigated sites were of the same magnitude as for thapsigargin (TG), the most potent SERCA inhibitor, or in the range between TG and di-tert-butylhydroquinone (BHQ), which is also known to possess inhibitory activity. The hydrophobic character of hypericin indicates that the molecule initially binds in the ER membrane from which it diffuses into the transmembrane region of the protein and to binding pockets therein. The transmembrane TG and BHQ binding pockets provide suitable locations for hypericin as they allow for favourable interactions with the lipid tails that surround these. High binding energies were noted for hypericin in these pockets and are expected to constitute highly possible binding sites due to their accessibility from the ER membrane. Hypericin most likely binds to both isoforms of SERCA and acts as an inhibitor or, under light irradiation, as a singlet oxygen generator that in turn degrades the protein or induces lipid peroxidation.

Molecular mechanism of the allosteric enhancement of the umami taste sensation

The fifth taste quality, umami, arises from binding of glutamate to the umami receptor T1R1/T1R3. The umami taste is enhanced several-fold upon addition of free nucleotides such as guanosine-5'-monophosphate (GMP) to glutamate-containing food. GMP may operate via binding to the ligand-binding domain of the T1R1 part of the umami receptor at an allosteric site. Using molecular dynamics simulations, we show that GMP can stabilize the closed (active) state of T1R1 by binding to the outer vestibule of the so-called Venus flytrap domain of the receptor. The transition between the closed and open conformations was accessed in the simulations. Using principal component analysis, we show that the dynamics of the Venus flytrap domain along the hinge-bending motion that activates signaling is dampened significantly upon binding of glutamate, and further slows down upon binding of GMP at an allosteric site, thus suggesting a molecular mechanism of cooperativity between GMP and glutamate.

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

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

Endoplasmic reticulum calcium pumps and cancer cell differentiation

The endoplasmic reticulum (ER) is a major intracellular calcium storage pool and a multifunctional organelle that accomplishes several calcium-dependent functions involved in many homeostatic and signaling mechanisms. Calcium is accumulated in the ER by Sarco/Endoplasmic Reticulum Calcium ATPase (SERCA)-type calcium pumps. SERCA activity can determine ER calcium content available for intra-ER functions and for calcium release into the cytosol, and can shape the spatiotemporal characteristics of calcium signals. SERCA function therefore constitutes an important nodal point in the regulation of cellular calcium homeostasis and signaling, and can exert important effects on cell growth, differentiation and survival. In several cell types such as cells of hematopoietic origin, mammary, gastric and colonic epithelium, SERCA2 and SERCA3-type calcium pumps are simultaneously expressed, and SERCA3 expression levels undergo significant changes during cell differentiation, activation or immortalization. In addition, SERCA3 expression is decreased or lost in several tumor types when compared to the corresponding normal tissue. These observations indicate that ER calcium homeostasis is remodeled during cell differentiation, and may present defects due to decreased SERCA3 expression in tumors. Modulation of the state of differentiation of the ER reflected by SERCA3 expression constitutes an interesting new aspect of cell differentiation and tumor biology.

Molecular dynamics simulations of the Ca2+-pump: a structural analysis

We report large scale molecular dynamics computer simulations, ∼100 ns, of the ion pump Ca(2+)-ATPase immersed in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer. The structure simulated here, E1, one of the several conformations resolved using X-ray diffraction techniques, hosts two Ca(2+)-ions in the hydrophobic domain. Our results indicate that protonated residues lead to stronger ion-residue interactions, supporting previous conclusions regarding the sensitivity of the Ca(2+) behaviour to the protonated state of the amino acid binding sites. We also investigate how the protein perturbs the bilayer structure. We show that the POPC bilayer is ∼12% thinner than the pure bilayer, near the protein surface. This perturbation decays exponentially with the distance from the protein with a characteristic decay length of 0.8 nm. We find that the projected area per lipid also decreases near the protein. Using an analytical model we show that this change in the area is only apparent and it can be explained by considering the local curvature of the membrane. Our results indicate that the real area per lipid near the protein is not significantly modified with respect to the pure bilayer result. Further our results indicate that the local deformation of the membrane around the protein might be compatible with the enhanced protein activity observed in experiments over a narrow range of membrane thicknesses.

Lipid recognition propensities of amino acids in membrane proteins from atomic resolution data

Background

Protein-lipid interactions play essential roles in the conformational stability and biological functions of membrane proteins. However, few of the previous computational studies have taken into account the atomic details of protein-lipid interactions explicitly.

Results

To gain an insight into the molecular mechanisms of the recognition of lipid molecules by membrane proteins, we investigated amino acid propensities in membrane proteins for interacting with the head and tail groups of lipid molecules. We observed a common pattern of lipid tail-amino acid interactions in two different data sources, crystal structures and molecular dynamics simulations. These interactions are largely explained by general lipophilicity, whereas the preferences for lipid head groups vary among individual proteins. We also found that membrane and water-soluble proteins utilize essentially an identical set of amino acids for interacting with lipid head and tail groups.

Conclusions

We showed that the lipophilicity of amino acid residues determines the amino acid preferences for lipid tail groups in both membrane and water-soluble proteins, suggesting that tightly-bound lipid molecules and lipids in the annular shell interact with membrane proteins in a similar manner. In contrast, interactions between lipid head groups and amino acids showed a more variable pattern, apparently constrained by each protein's specific molecular function.

Protonation states of important acidic residues in the central Ca²⁺ ion binding sites of the Ca²⁺-ATPase: a molecular modeling study

The P-type ATPases are responsible for the transport of cations across cell membranes. The sarco(endo)plasmic reticulum Ca²⁺-ATPase (SERCA) transports two Ca²⁺ ions from the cytoplasm to the lumen of the sarco(endo)plasmic reticulum and countertransports two or three protons per catalytic cycle. Two binding sites for Ca²⁺ ions have been located via protein crystallography, including four acidic amino acid residues that are essential to the ion coordination. In this study, we present molecular dynamics (MD) simulations examining the protonation states of these amino acid residues in a Ca²⁺-free conformation of SERCA. Such knowledge will be important for an improved understanding of atomistic details of the transport mechanism of protons and Ca²⁺ ions. Eight combinations of the protonation of four central acidic residues, Glu309, Glu771, Asp800, and Glu908, are tested from 10 ns MD simulations with respect to protein stability and ability to maintain a structure similar to the crystal structure. The trajectories for the most prospective combinations of protonation states were elongated to 50 ns and subjected to more detailed analysis, including prediction of pK(a) values of the four acidic residues over the trajectories. From the simulations we find that the combination leaving only Asp800 as charged is most likely. The results are compared to available experimental data and explain the observed destabilization upon full deprotonation, resulting in the entry of cytoplasmic K⁺ ions into the Ca²⁺ binding sites during the simulation in which Ca²⁺ ions are absent. Furthermore, a hypothesis for the exchange of protons from the central binding cavity is proposed.

Atomic-level characterization of the activation mechanism of SERCA by calcium

We have performed molecular dynamics (MD) simulations to elucidate, in atomic detail, the mechanism by which the sarcoplasmic reticulum Ca²⁺-ATPase (SERCA) is activated by Ca²⁺. Crystal structures suggest that activation of SERCA occurs when the cytoplasmic head-piece, in an open (E1) conformation stabilized by Ca²⁺, undergoes a large-scale open-to-closed (E1 to E2) transition that is induced by ATP binding. However, spectroscopic measurements in solution suggest that these structural states (E1 and E2) are not tightly coupled to biochemical states (defined by bound ligands); the closed E2 state predominates even in the absence of ATP, in both the presence and absence of Ca²⁺. How is this loose coupling consistent with the high efficiency of energy transduction in the Ca²⁺-ATPase? To provide insight into this question, we performed long (500 ns) all-atom MD simulations starting from the open crystal structure, including a lipid bilayer and water. In both the presence and absence of Ca²⁺, we observed a large-scale open-to-closed conformational transition within 400 ns, supporting the weak coupling between structural and biochemical states. However, upon closer inspection, it is clear that Ca²⁺ is necessary and sufficient for SERCA to reach the precise geometrical arrangement necessary for activation of ATP hydrolysis. Contrary to suggestions from crystal structures, but in agreement with solution spectroscopy, the presence of ATP is not required for this activating transition. Principal component analysis showed that Ca²⁺ reshapes the free energy landscape of SERCA to create a path between the open conformation and the activated closed conformation. Thus the malleability of the free energy landscape is essential for SERCA efficiency, ensuring that ATP hydrolysis is tightly coupled to Ca²⁺ transport. These results demonstrate the importance of real-time dynamics in the formation of catalytically competent conformations of SERCA, with broad implications for understanding enzymatic catalysis in atomic detail.