The Role of Extramembranous Cytoplasmic Termini in Assembly and Stability of the Tetrameric K+-Channel KcsA

Membrane-active alcohol 2,2,2-trifluoroethanol has been proven to be an attractive tool in the investigation of the intrinsic stability of integral membrane protein complexes by taking K+-channel KcsA as a suitable and representative ion channel. In the present study, the roles of both cytoplasmic N and C termini in channel assembly and stability of KcsA were determined. The N terminus (1–18 residues) slightly increased tetramer stability via electrostatic interactions in the presence of 30 mol.% acidic phosphatidylglycerol (PG) in phosphatidylcholine lipid bilayer. Furthermore, the N terminus was found to be potentially required for efficient channel (re)assembly. In contrast, truncation of the C terminus (125–160 residues) greatly facilitated channel reversibility from either a partially or a completely unfolded state, and this domain was substantially involved in stabilizing the tetramer in either the presence or absence of PG in lipid bilayer. These studies provide new insights into how extramembranous parts play their crucial roles in the assembly and stability of integral membrane protein complexes

Conformational Changes and Slow Dynamics through Microsecond Polarized Atomistic Molecular Simulation of an Integral Kv1.2 Ion Channel

Structure and dynamics of voltage-gated ion channels, in particular the motion of the S4 helix, is a highly interesting and hotly debated topic in current membrane protein research. It has critical implications for insertion and stabilization of membrane proteins as well as for finding how transitions occur in membrane proteins—not to mention numerous applications in drug design. Here, we present a full 1 µs atomic-detail molecular dynamics simulation of an integral Kv1.2 ion channel, comprising 120,000 atoms. By applying 0.052 V/nm of hyperpolarization, we observe structural rearrangements, including up to 120° rotation of the S4 segment, changes in hydrogen-bonding patterns, but only low amounts of translation. A smaller rotation (∼35°) of the extracellular end of all S4 segments is present also in a reference 0.5 µs simulation without applied field, which indicates that the crystal structure might be slightly different from the natural state of the voltage sensor. The conformation change upon hyperpolarization is closely coupled to an increase in 310 helix contents in S4, starting from the intracellular side. This could support a model for transition from the crystal structure where the hyperpolarization destabilizes S4–lipid hydrogen bonds, which leads to the helix rotating to keep the arginine side chains away from the hydrophobic phase, and the driving force for final relaxation by downward translation is partly entropic, which would explain the slow process. The coordinates of the transmembrane part of the simulated channel actually stay closer to the recently determined higher-resolution Kv1.2 chimera channel than the starting structure for the entire second half of the simulation (0.5–1 µs). Together with lipids binding in matching positions and significant thinning of the membrane also observed in experiments, this provides additional support for the predictive power of microsecond-scale membrane protein simulations

A neuroscientist's guide to lipidomics

Nerve cells mould the lipid fabric of their membranes to ease vesicle fusion, regulate ion fluxes and create specialized microenvironments that contribute to cellular communication. The chemical diversity of membrane lipids controls protein traffic, facilitates recognition between cells and leads to the production of hundreds of molecules that carry information both within and across cells. With so many roles, it is no wonder that lipids make up half of the human brain in dry weight. The objective of neural lipidomics is to understand how these molecules work together; this difficult task will greatly benefit from technical advances that might enable the testing of emerging hypotheses

MD simulations of spontaneous membrane protein/detergent micelle formation.

The in vitro study of membrane proteins for the purpose of physicochemical analysis or structure determination often relies upon successful reconstitution into detergent micelles. Moreover, a number of biological processes such as membrane protein folding and transport rely on lipid interactions which may resemble the micellar environment. Little is known about the structures of these micelles or the processes which lead to their formation. We therefore present two 50 ns all-atom molecular dynamics simulations of spontaneous dodecylphosphocholine micelle formation around representatives of the two major families of membrane proteins, a small beta-barrel protein, OmpA, and a model alpha-helical protein, glycophorin A. Despite differences in protein architecture, we highlight common mechanistic pathways in micelle formation, which are consistent with experimental studies. We characterize the exponential kinetics of detergent-protein adsorption and suggest a simple model which may explain the aggregation process. We also compare the results with 25 and 50 ns simulations of preformed micelles containing the same proteins. We confirm that the end structures of the self-assembled micelles are similar to those from their preformed counterparts, with each micelle presenting a bilayerlike environment to the enclosed protein

Molecular dynamics simulations of GlpF in a micelle vs in a bilayer: conformational dynamics of a membrane protein as a function of environment.

Octyl glucoside (OG) is a detergent widely employed in structural and functional studies of membrane proteins. To better understand the nature of protein-OG interactions, molecular dynamics simulations (duration 10 ns) have been used to explore an alpha-helical membrane protein, GlpF, in OG micelles and in DMPC bilayers. Greater conformational drift of the extramembraneous protein loops, from the initial X-ray structure, is seen for the GlpF-OG simulations than for the GlpF-DMPC simulation. The mobility of the transmembrane alpha-helices is approximately 1.3x higher in the GlpF-OG than the GlpF-DMPC simulations. The detergent is seen to form an irregular torus around the protein. The presence of the protein leads to a small perturbation in the behavior of the alkyl chains in the OG micelle, namely an approximately 15% increase in the trans-gauche(-)-gauche(+) transition time. Aromatic side chains (Trp, Tyr) and basic side chains (Arg, Lys) play an important role in both protein-detergent (OG) and protein-lipid (DMPC) interactions

Lipid/protein interactions and the membrane/water interfacial region.

Despite the importance of lipid/protein interactions in the folding, assembly, stability, and function of membrane proteins, information at an atomic level on how such proteins interact with the lipids that surround them remains sparse. The dynamics and flexible nature of the protein/bilayer interaction make it difficult to study, for example, by crystallographic means. However, based on recent progress in molecular simulations of membranes it is possible to address this problem computationally. This communication reports one of the first attempts to use multiple ns molecular simulations to establish a qualitative picture of the intermolecular interactions between the lipids of a bilayer and two topologically different membrane proteins for which a high resolution (2 A or better) X-ray structure is available

Molecular simulations and lipid-protein interactions: potassium channels and other membrane proteins.

Molecular dynamics simulations may be used to probe the interactions of membrane proteins with lipids and with detergents at atomic resolution. Examples of such simulations for ion channels and for bacterial outer membrane proteins are described. Comparison of simulations of KcsA (an alpha-helical bundle) and OmpA (a beta-barrel) reveals the importance of two classes of side chains in stabilizing interactions with the head groups of lipid molecules: (i) tryptophan and tyrosine; and (ii) arginine and lysine. Arginine residues interacting with lipid phosphate groups play an important role in stabilizing the voltage-sensor domain of the KvAP channel within a bilayer. Simulations of the bacterial potassium channel KcsA reveal specific interactions of phosphatidylglycerol with an acidic lipid-binding site at the interface between adjacent protein monomers. A combination of molecular modelling and simulation reveals a potential phosphatidylinositol 4,5-bisphosphate-binding site on the surface of Kir6.2