Structure, dynamics, and ion conductance of the phospholamban pentamer

Molecular noise filtering in the β-adrenergic signaling network by phospholamban pentamers

Phospholamban (PLN) is an important regulator of cardiac calcium handling due to its ability to inhibit the calcium ATPase SERCA. β-Adrenergic stimulation reverses SERCA inhibition via PLN phosphorylation and facilitates fast calcium reuptake. PLN also forms pentamers whose physiological significance has remained elusive. Using mathematical modeling combined with biochemical and cell biological experiments, we show that pentamers regulate both the dynamics and steady-state levels of monomer phosphorylation. Substrate competition by pentamers and a feed-forward loop involving inhibitor-1 can delay monomer phosphorylation by protein kinase A (PKA), whereas cooperative pentamer dephosphorylation enables bistable PLN steady-state phosphorylation. Simulations show that phosphorylation delay and bistability act as complementary filters that reduce the effect of random fluctuations in PKA activity, thereby ensuring consistent monomer phosphorylation and SERCA activity despite noisy upstream signals. Preliminary analyses suggest that the PLN mutation R14del could impair noise filtering, offering a new perspective on how this mutation causes cardiac arrhythmias.

Perturbations of Native Membrane Protein Structure in Alkyl Phosphocholine Detergents: A Critical Assessment of NMR and Biophysical Studies

Membrane proteins perform a host of vital cellular functions. Deciphering the molecular mechanisms whereby they fulfill these functions requires detailed biophysical and structural investigations. Detergents have proven pivotal to extract the protein from its native surroundings. Yet, they provide a milieu that departs significantly from that of the biological membrane, to the extent that the structure, the dynamics, and the interactions of membrane proteins in detergents may considerably vary, as compared to the native environment. Understanding the impact of detergents on membrane proteins is, therefore, crucial to assess the biological relevance of results obtained in detergents. Here, we review the strengths and weaknesses of alkyl phosphocholines (or foscholines), the most widely used detergent in solution-NMR studies of membrane proteins. While this class of detergents is often successful for membrane protein solubilization, a growing list of examples points to destabilizing and denaturing properties, in particular for α-helical membrane proteins. Our comprehensive analysis stresses the importance of stringent controls when working with this class of detergents and when analyzing the structure and dynamics of membrane proteins in alkyl phosphocholine detergents.

Mitochondrial ADP/ATP Carrier in Dodecylphosphocholine Binds Cardiolipins with Non-native Affinity

Biophysical investigation of membrane proteins generally requires their extraction from native sources using detergents, a step that can lead, possibly irreversibly, to protein denaturation. The propensity of dodecylphosphocholine (DPC), a detergent widely utilized in NMR studies of membrane proteins, to distort their structure has been the subject of much controversy. It has been recently proposed that the binding specificity of the yeast mitochondrial ADP/ATP carrier (yAAC3) toward cardiolipins is preserved in DPC, thereby suggesting that DPC is a suitable environment in which to study membrane proteins. In this communication, we used all-atom molecular dynamics simulations to investigate the specific binding of cardiolipins to yAAC3. Our data demonstrate that the interaction interface observed in a native-like environment differs markedly from that inferred from an NMR investigation in DPC, implying that in this detergent, the protein structure is distorted. We further investigated yAAC3 solubilized in DPC and in the milder dodecylmaltoside with thermal-shift assays. The loss of thermal transition observed in DPC confirms that the protein is no longer properly folded in this environment.

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.

A computational study of ion current modulation in hVDAC3 induced by disulfide bonds

The human VDAC channel exists in three isoforms characterized by high sequence homology and structural similarity. Yet the function and mode of action of hVDAC3 are still elusive. The presence of six surface cysteines exposed to the oxidizing environment of the mitochondrial inter-membrane space suggests the possible establishment of intramolecular disulfide bonds. Two natural candidates for disulfide bridge formation are Cys2 and Cys8 that, located on the flexible N-terminal domain, can easily come in contact. A third potentially important residue is Cys122 that is close to Cys2 in the homology model of VDAC3. Here we analyzed the impact of SS bonds through molecular dynamics simulations of derivatives of hVDAC3 (dubbed SS-2-8, SS-2-122, SS-8-122) including a single disulfide bond. Simulations showed that in SS-8-122, the fragment 1-7 crosses the top part of the barrel partially occluding the pore and causing a 20% drop of conductance. In order to identify other potential channel-occluding disulfide bonds, we used a set of neural networks and structural bioinformatics algorithms, after filtering with the steric constraints imposed by the 3D-structure. We identified other three species, namely SS-8-65, SS-2-36 and SS-8-36. While the conductance of SS-8-65 and SS-2-36 is about 30% lower than that of the species without disulfide bonds, the conductance of SS-8-36 was 40-50% lower. The results show how VDAC3 is able to modulate its pore size and current by exploiting the mobility of the N-terminal and forming, upon external stimuli, disulfide bridges with cysteine residues located on the barrel and exposed to the inter-membrane space.

Simulations of outer membrane channels and their permeability

Channels in the outer membrane of Gram-negative bacteria provide essential pathways for the controlled and unidirectional transport of ions, nutrients and metabolites into the cell. At the same time the outer membrane serves as a physical barrier for the penetration of noxious substances such as antibiotics into the bacteria. Most antibiotics have to pass through these membrane channels to either reach cytoplasmic bound targets or to further cross the hydrophobic inner membrane. Considering the pharmaceutical significance of antibiotics, understanding the functional role and mechanism of these channels is of fundamental importance in developing strategies to design new drugs with enhanced permeation abilities. Due to the biological complexity of membrane channels and experimental limitations, computer simulations have proven to be a powerful tool to investigate the structure, dynamics and interactions of membrane channels. Considerable progress has been made in computer simulations of membrane channels during the last decade. The goal of this review is to provide an overview of the computational techniques and their roles in modeling the transport across outer membrane channels. A special emphasis is put on all-atom molecular dynamics simulations employed to better understand the transport of molecules. Moreover, recent molecular simulations of ion, substrate and antibiotics translocation through membrane pores are briefly summarized. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Structure-function relation of phospholamban: modulation of channel activity as a potential regulator of SERCA activity

Phospholamban (PLN) is a small integral membrane protein, which binds and inhibits in a yet unknown fashion the Ca²⁺-ATPase (SERCA) in the sarcoplasmic reticulum. When reconstituted in planar lipid bilayers PLN exhibits ion channel activity with a low unitary conductance. From the effect of non-electrolyte polymers on this unitary conductance we estimate a narrow pore with a diameter of ca. 2.2 Å for this channel. This value is similar to that reported for the central pore in the structure of the PLN pentamer. Hence the PLN pentamer, which is in equilibrium with the monomer, is the most likely channel forming structure. Reconstituted PLN mutants, which either stabilize (K27A and R9C) or destabilize (I47A) the PLN pentamer and also phosphorylated PLN still generate the same unitary conductance of the wt/non-phosphorylated PLN. However the open probability of the phosphorylated PLN and of the R9C mutant is significantly lower than that of the respective wt/non-phosphorylated control. In the context of data on PLN/SERCA interaction and on Ca²⁺ accumulation in the sarcoplasmic reticulum the present results are consistent with the view that PLN channel activity could participate in the balancing of charge during Ca²⁺ uptake. A reduced total conductance of the K⁺ transporting PLN by phosphorylation or by the R9C mutation may stimulate Ca²⁺ uptake in the same way as an inhibition of K⁺ channels in the SR membrane. The R9C-PLN mutation, a putative cause of dilated cardiomyopathy, might hence affect SERCA activity also via its inherent low open probability.

Cation selectivity is a conserved feature in the OccD subfamily of Pseudomonas aeruginosa

To achieve the uptake of small, water-soluble nutrients, Pseudomonas aeruginosa, a pathogenic Gram-negative bacterium, employs substrate-specific channels located within its outer membrane. In this paper, we present a detailed description of the single-channel characteristics of six members of the outer membrane carboxylate channel D (OccD) subfamily. Recent structural studies showed that the OccD proteins share common features, such as a closely related, monomeric, 18-stranded β-barrel conformation and large extracellular loops, which are folded back into the channel lumen. Here, we report that the OccD proteins displayed single-channel activity with a unitary conductance covering an unusually broad range, between 20 and 670pS, as well as a diverse gating dynamics. Interestingly, we found that cation selectivity is a conserved trait among all members of the OccD subfamily, bringing a new distinction between the members of the OccD subfamily and the anion-selective OccK channels. Conserved cation selectivity of the OccD channels is in accord with an increased specificity and selectivity of these proteins for positively charged, carboxylate-containing substrates.

Constant electric field simulations of the membrane potential illustrated with simple systems

Advances in modern computational methods and technology make it possible to carry out extensive molecular dynamics simulations of complex membrane proteins based on detailed atomic models. The ultimate goal of such detailed simulations is to produce trajectories in which the behavior of the system is as realistic as possible. A critical aspect that requires consideration in the case of biological membrane systems is the existence of a net electric potential difference across the membrane. For meaningful computations, it is important to have well validated methodologies for incorporating the latter in molecular dynamics simulations. A widely used treatment of the membrane potential in molecular dynamics consists of applying an external uniform electric field E perpendicular to the membrane. The field acts on all charged particles throughout the simulated system, and the resulting applied membrane potential V is equal to the applied electric field times the length of the periodic cell in the direction perpendicular to the membrane. A series of test simulations based on simple membrane-slab models are carried out to clarify the consequences of the applied field. These illustrative tests demonstrate that the constant-field method is a simple and valid approach for accounting for the membrane potential in molecular dynamics studies of biomolecular systems. This article is part of a Special Issue entitled: Membrane protein structure and function.

Structural topology of phospholamban pentamer in lipid bilayers by a hybrid solution and solid-state NMR method

Phospholamban (PLN) is a type II membrane protein that inhibits the sarcoplasmic reticulum Ca(2+)-ATPase (SERCA), thereby regulating calcium homeostasis in cardiac muscle. In membranes, PLN forms pentamers that have been proposed to function either as a storage for active monomers or as ion channels. Here, we report the T-state structure of pentameric PLN solved by a hybrid solution and solid-state NMR method. In lipid bilayers, PLN adopts a pinwheel topology with a narrow hydrophobic pore, which excludes ion transport. In the T state, the cytoplasmic amphipathic helices (domains Ia) are absorbed into the lipid bilayer with the transmembrane domains arranged in a left-handed coiled-coil configuration, crossing the bilayer with a tilt angle of approximately 11° with respect to the membrane normal. The tilt angle difference between the monomer and pentamer is approximately 13°, showing that intramembrane helix-helix association forces dominate over the hydrophobic mismatch, driving the overall topology of the transmembrane assembly. Our data reveal that both topology and function of PLN are shaped by the interactions with lipids, which fine-tune the regulation of SERCA.

An allosteric mechanism inferred from molecular dynamics simulations on phospholamban pentamer in lipid membranes

Phospholamban functions as a regulator of Ca(2+) concentration of cardiac muscle cells by triggering the bioactivity of sarcoplasmic reticulum Ca(2+)-ATPase. In order to understand its dynamic mechanism in the environment of bilayer surroundings, we performed long time-scale molecular dynamic simulations based on the high-resolution NMR structure of phospholamban pentamer. It was observed from the molecular dynamics trajectory analyses that the conformational transitions between the "bellflower" and "pinwheel" modes were detected for phospholamban. Particularly, the two modes became quite similar to each other after phospholamban was phosphorylated at Ser16. Based on these findings, an allosteric mechanism was proposed to elucidate the dynamic process of phospholamban interacting with Ca(2+)-ATPase.