Dataset Information

Molecular dynamics simulation of lipid reorientation at bilayer edges.

ABSTRACT: Understanding cellular membrane processes is critical for the study of events such as viral entry, neurotransmitter exocytosis, and immune activation. Supported lipid bilayers are commonly used to model these membrane processes experimentally. Despite the relative simplicity of such a system, many important structural and dynamic parameters are not experimentally observable with current techniques. Computational approaches allow the development of a high-resolution model of bilayer processes. We have performed molecular dynamics simulations of dimyristoylphosphatidylcholine (DMPC) bilayers to model the creation of bilayer gaps-a common process in bilayer patterning-and to analyze their structure and dynamics. We propose a model for gap formation in which the bilayer edges form metastable micelle-like structures on a nanosecond timescale. Molecules near edges structurally resemble lipids in ungapped bilayers but undergo small-scale motions more rapidly. These data suggest that lipids may undergo rapid local rearrangements during membrane fusion, facilitating the formation of fusion intermediates thought key to the infection cycle of viruses such as influenza, Ebola, and HIV.

SUBMITTER: Kasson PM

PROVIDER: S-EPMC1304275 | biostudies-literature | 2004 Jun

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Molecular dynamics simulation of lipid reorientation at bilayer edges.

Kasson Peter M PM Pande Vijay S VS

Biophysical journal 20040601 6

Understanding cellular membrane processes is critical for the study of events such as viral entry, neurotransmitter exocytosis, and immune activation. Supported lipid bilayers are commonly used to model these membrane processes experimentally. Despite the relative simplicity of such a system, many important structural and dynamic parameters are not experimentally observable with current techniques. Computational approaches allow the development of a high-resolution model of bilayer processes. We ...[more]

PMID: 15189870

Similar Datasets

Project description:Transient receptor potential vanilloid subtype 1 (TRPV1) is a heat-sensitive ion channel also involved in pain sensation, and is the receptor for capsaicin, the active ingredient of hot chili peppers. The recent structures of TRPV1 revealed putative ligand density within the S1 to S4 voltage-sensor-like domain of the protein. However, questions remain regarding the dynamic role of the lipid bilayer in ligand binding to TRPV1. Molecular dynamics simulations were used to explore behavior of capsaicin in a 1-palmitoyl-2-oleoyl phosphatidylcholine bilayer and with the target S1-S4 transmembrane helices of TRPV1. Equilibrium simulations reveal a preferred interfacial localization for capsaicin. We also observed a capsaicin molecule flipping from the extracellular to the intracellular leaflet, and subsequently able to access the intracellular TRPV1 binding site. Calculation of the potential of mean force (i.e., free energy profile) of capsaicin along the bilayer normal confirms that it prefers an interfacial localization. The free energy profile indicates that there is a nontrivial but surmountable barrier to the flipping of capsaicin between opposing leaflets of the bilayer. Molecular dynamics of the S1-S4 transmembrane helices of the TRPV1 in a lipid bilayer confirm that Y511, known to be crucial to capsaicin binding, has a distribution along the bilayer normal similar to that of the aromatic group of capsaicin. Simulations were conducted of the TRPV1 S1-S4 transmembrane helices in the presence of capsaicin placed in the aqueous phase, in the lipid, or docked to the protein. No stable interaction between ligand and protein was seen for simulations initiated with capsaicin in the bilayer. However, interactions were seen between TRPV1 and capsaicin starting from the cytosolic aqueous phase, and capsaicin remained stable in the majority of simulations from the docked pose. We discuss the significance of capsaicin flipping from the extracellular to the intracellular leaflet and mechanisms of binding site access by capsaicin.

Project description:All atom molecular dynamics simulations of the 18-residue beta-hairpin antimicrobial peptide protegrin-1 (PG-1, RGGRLCYCRRRFCVCVGR-NH(2)) in a fully hydrated dilauroylphosphatidylcholine (DLPC) lipid bilayer have been implemented. The goal of the reported work is to investigate the structure of the peptide in a membrane environment (previously solved only in solution [R.L. Fahrner, T. Dieckmann, S.S.L. Harwig, R.I. Lehrer, D. Eisenberg, J. Feigon, Solution structure of protegrin-1, a broad-spectrum antimicrobial peptide from porcine leukocytes. Chemistry and Biology, 3 (1996) 543-550]), and to delineate specific peptide-membrane interactions which are responsible for the peptide's membrane binding properties. A novel, previously unknown, "kick" shaped conformation of the peptide was detected, where a bend at the C-terminal beta-strand of the peptide caused the peptide backbone at residues 16-18 to extend perpendicular to the beta-hairpin plane. This bend was driven by a highly persistent hydrogen-bond between the polar peptide side-chain of TYR7 and the unshielded backbone carbonyl oxygen atom of GLY17. The H-bond formation relieves the unfavorable free energy of insertion of polar groups into the hydrophobic membrane core. PG-1 was anchored to the membrane by strong electrostatic binding of the protonated N-terminus of the peptide to the lipid head group phosphate anions. The orientation of the peptide in the membrane, and its influence on bilayer structural and dynamic properties are in excellent agreement with solid state NMR measurements [S. Yamaguchi, T. Hong, A. Waring, R.I. Lehrer, M. Hong, Solid-State NMR Investigations of Peptide-Lipid Interaction and Orientation of a b-Sheet Antimicrobial Peptide, Protegrin, Biochemistry, 41 (2002) 9852-9862]. Importantly, two simulations which started from different initial orientations of the peptide converged to the same final equilibrium orientation of the peptide relative to the bilayer. The kick-shaped conformation was observed only in one of the two simulations.

Project description:Spin-labeled lipids are commonly used as fluorescence quenchers in studies of membrane penetration of dye-labeled proteins and peptides using depth-dependent quenching. Accurate calculations of depth of the fluorophore rely on the use of several spin labels placed in the membrane at various positions. The depth of the quenchers (spin probes) has to be determined independently; however, experimental determination of transverse distributions of spin probe depths is difficult. In this Article, we use molecular dynamics (MD) simulations to study the membrane behavior and depth distributions of spin-labeled phospholipids in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) bilayer. To probe different depths within the bilayer, a series containing five Doxyl-labeled lipids (n-Doxyl PC) has been studied, in which a spin moiety was covalently attached to nth carbon atoms (where n = 5, 7, 10, 12, and 14) of the sn-2 stearoyl chain of the host phospholipid. Our results demonstrate that the chain-attached spin labels are broadly distributed across the model membrane and their environment is characterized by a high degree of mobility and structural heterogeneity. Despite the high thermal disorder, the depth distributions of the Doxyl labels were found to correlate well with their attachment positions, indicating that the distribution of the spin label within the model membrane is dictated by the depth of the nth lipid carbon atom and not by intrinsic properties of the label. In contrast, a much broader and heterogeneous distribution was observed for a headgroup-attached Tempo spin label of Tempo-PC lipids. MD simulations reveal that, due to the hydrophobic nature, a Tempo moiety favors partitioning from the headgroup region deeper into the membrane. Depending on the concentration of Tempo-PC lipids, the probable depth of the Tempo moiety could span a range from 14.4 to 18.2 Å from the membrane center. Comparison of the MD-estimated immersion depths of Tempo and n-Doxyl labels with their suggested experimental depth positions allows us to review critically the possible sources of error in depth-dependent fluorescence quenching studies.

Dataset Information

Molecular dynamics simulation of lipid reorientation at bilayer edges.

Publications

Molecular dynamics simulation of lipid reorientation at bilayer edges.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets