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

Emerging Diversity in Lipid-Protein Interactions

Membrane lipids interact with proteins in a variety of ways, ranging from providing a stable membrane environment for proteins to being embedded in to detailed roles in complicated and well-regulated protein functions. Experimental and computational advances are converging in a rapidly expanding research area of lipid-protein interactions. Experimentally, the database of high-resolution membrane protein structures is growing, as are capabilities to identify the complex lipid composition of different membranes, to probe the challenging time and length scales of lipid-protein interactions, and to link lipid-protein interactions to protein function in a variety of proteins. Computationally, more accurate membrane models and more powerful computers now enable a detailed look at lipid-protein interactions and increasing overlap with experimental observations for validation and joint interpretation of simulation and experiment. Here we review papers that use computational approaches to study detailed lipid-protein interactions, together with brief experimental and physiological contexts, aiming at comprehensive coverage of simulation papers in the last five years. Overall, a complex picture of lipid-protein interactions emerges, through a range of mechanisms including modulation of the physical properties of the lipid environment, detailed chemical interactions between lipids and proteins, and key functional roles of very specific lipids binding to well-defined binding sites on proteins. Computationally, despite important limitations, molecular dynamics simulations with current computer power and theoretical models are now in an excellent position to answer detailed questions about lipid-protein interactions.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

Distinct pH dependencies of Na+/K+ selectivity at the two faces of Na,K-ATPase

The sodium pump (Na,K-ATPase) in animal cells is vital for actively maintaining ATP hydrolysis-powered Na⁺ and K⁺ electrochemical gradients across the cell membrane. These ion gradients drive co- and countertransport and are critical for establishing the membrane potential. It has been an enigma how Na,K-ATPase discriminates between Na⁺ and K⁺, despite the pumped ion on each side being at a lower concentration than the other ion. Recent crystal structures of analogs of the intermediate conformations E2·Pi·2K⁺ and Na⁺-bound E1∼P·ADP suggest that the dimensions of the respective binding sites in Na,K-ATPase are crucial in determining its selectivity. Here, we found that the selectivity at each membrane face is pH-dependent and that this dependence is unique for each face. Most notable was a strong increase in the specific affinity for K⁺ at the extracellular face (i.e. E2 conformation) as the pH is lowered from 7.5 to 5. We also observed a smaller increase in affinity for K⁺ on the cytoplasmic side (E1 conformation), which reduced the selectivity for Na⁺ Theoretical analysis of the pK_a values of ion-coordinating acidic amino acid residues suggested that the face-specific pH dependences and Na⁺/K⁺ selectivities may arise from the protonation or ionization of key residues. The increase in K⁺ selectivity at low pH on the cytoplasmic face, for instance, appeared to be associated with Asp⁸⁰⁸ protonation. We conclude that changes in the ionization state of coordinating residues in Na,K-ATPase could contribute to altering face-specific ion selectivity.

How does ytterbium chloride interact with DMPC bilayers? A computational and experimental study

Lanthanide salts have been studied for many years, primarily in Nuclear Magnetic Resonance (NMR) experiments of mixed lipid–protein systems and more recently to study lipid flip-flop in model membrane systems.

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.

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.

On the thermodynamic efficiency of Ca²⁺-ATPase molecular machines

Experimental studies have shown that the activity of the reconstituted molecular pump Ca(2+)-ATPase strongly depends on the thickness of the supporting bilayer. It is thus expected that the bilayer structure will have an impact on the thermodynamic efficiency of this nanomachine. Here, we introduce a nonequilibrium-thermodynamics theoretical approach to estimate the thermodynamic efficiency of the Ca(2+)-ATPase from analysis of available experimental data about ATP hydrolysis and Ca(2+) transport. We find that the entropy production, i.e., the heat released to the surroundings under working conditions, is approximately constant for bilayers containing phospholipids with hydrocarbon chains of 18-22 carbon atoms. Our estimates for the heat released during the pump operation agree with results obtained from separate calorimetric experiments on the Ca(2+)-ATPase derived from sarcoplasmic reticulum. We show that the thermodynamic efficiency of the reconstituted Ca(2+)-ATPase reaches a maximum for bilayer thicknesses corresponding to maximum activity. Surprisingly, the estimated thermodynamic efficiency is very low, ∼12%. We discuss the significance of this result as representative of the efficiency of other nanomachines, and we address the influence of the experimental set-up on such a low efficiency. Overall, our approach provides a general route to estimate thermodynamic efficiencies and heat dissipation in experimental studies of nanomachines.

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.

Flexible P-type ATPases interacting with the membrane

Cation pumps and lipid flippases of the P-type ATPase family maintain electrochemical gradients and asymmetric lipid distributions across membranes, and offer significant insight of protein:membrane interactions. The sarcoplasmic reticulum Ca(2+)-ATPase features flexible and adaptive interactions with the surrounding membrane, while the Na(+),K(+)-ATPase complex is modulated by membrane components and a role for the γ-subunit as a stabilizer of a specific lipid interaction with the α-subunit has been proposed. The first crystal structure of a heavy-metal transporting ATPase shows a markedly amphipathic helix at the cytoplasmic membrane surface, highlighting this structure as a general motif of all P-type ATPases although with specialization to different membranes. Residues of central importance for the lipid flippase activity of the P4-type ATPase subfamily have been pinpointed by mutational studies, but the transport pathway and mechanism remain unknown.