Positioning of proteins in membranes: a computational approach

The contribution of surface residues to membrane binding and ligand transfer by the α-tocopherol transfer protein (α-TTP)

Previous work has shown that the α-tocopherol transfer protein (α-TTP) can bind to vesicular or immobilized phospholipid membranes. Revealing the molecular mechanisms by which α-TTP associates with membranes is thought to be critical to understanding its function and role in the secretion of tocopherol from hepatocytes into the circulation. Calculations presented in the Orientations of Proteins in Membranes database have provided a testable model for the spatial arrangement of α-TTP and other CRAL-TRIO family proteins with respect to the lipid bilayer. These calculations predicted that a hydrophobic surface mediates the interaction of α-TTP with lipid membranes. To test the validity of these predictions, we used site-directed mutagenesis and examined the substituted mutants with regard to intermembrane ligand transfer, association with lipid layers and biological activity in cultured hepatocytes. Substitution of residues in helices A8 (F165A and F169A) and A10 (I202A, V206A and M209A) decreased the rate of intermembrane ligand transfer as well as protein adsorption to phospholipid bilayers. The largest impairment was observed upon mutation of residues that are predicted to be fully immersed in the lipid bilayer in both apo (open) and holo (closed) conformations such as Phe165 and Phe169. Mutation F169A, and especially F169D, significantly impaired α-TTP-assisted secretion of α-tocopherol outside cultured hepatocytes. Mutation of selected basic residues (R192H, K211A, and K217A) had little effect on transfer rates, indicating no significant involvement of nonspecific electrostatic interactions with membranes.

The mechanism of sodium and substrate release from the binding pocket of vSGLT

Membrane co-transport proteins that use a five-helix inverted repeat motif have recently emerged as one of the largest structural classes of secondary active transporters. However, despite many structural advances there is no clear evidence of how ion and substrate transport are coupled. Here we report a comprehensive study of the sodium/galactose transporter from Vibrio parahaemolyticus (vSGLT), consisting of molecular dynamics simulations, biochemical characterization and a new crystal structure of the inward-open conformation at a resolution of 2.7 Å. Our data show that sodium exit causes a reorientation of transmembrane helix 1 that opens an inner gate required for substrate exit, and also triggers minor rigid-body movements in two sets of transmembrane helical bundles. This cascade of events, initiated by sodium release, ensures proper timing of ion and substrate release. Once set in motion, these molecular changes weaken substrate binding to the transporter and allow galactose readily to enter the intracellular space. Additionally, we identify an allosteric pathway between the sodium-binding sites, the unwound portion of transmembrane helix 1 and the substrate-binding site that is essential in the coupling of co-transport.

Interaction of calmodulin with Sec61α limits Ca2+ leakage from the endoplasmic reticulum

In eukaryotes, protein transport into the endoplasmic reticulum (ER) is facilitated by a protein-conducting channel, the Sec61 complex. The presence of large, water-filled pores with uncontrolled ion permeability, as formed by Sec61 complexes in the ER membrane, would seriously interfere with the regulated release of calcium from the ER lumen into the cytosol, an essential mechanism for intracellular signalling. We identified a calmodulin (CaM)-binding motif in the cytosolic N-terminus of mammalian Sec61α that bound CaM but not Ca2+-free apocalmodulin with nanomolar affinity and sequence specificity. In single-channel measurements, CaM potently mediated Sec61-channel closure in Ca2+-dependent manner. At the cellular level, two different CaM antagonists stimulated calcium release from the ER through Sec61 channels. However, protein transport into microsomes was not modulated by Ca2+-CaM. Molecular modelling of the ribosome/Sec61/CaM complexes supports the view that simultaneous ribosome and CaM binding to the Sec61 complex may be possible. Overall, CaM is involved in limiting Ca2+ leakage from the ER.

Open and closed conformations of two SpoIIAA-like proteins (YP_749275.1 and YP_001095227.1) provide insights into membrane association and ligand binding

The crystal structures of the proteins encoded by the YP_749275.1 and YP_001095227.1 genes from Shewanella frigidimarina and S. loihica, respectively, have been determined at 1.8 and 2.25 Å resolution, respectively. These proteins are members of a novel family of bacterial proteins that adopt the α/β SpoIIAA-like fold found in STAS and CRAL-TRIO domains. Despite sharing 54% sequence identity, these two proteins adopt distinct conformations arising from different dispositions of their α2 and α3 helices. In the `open' conformation (YP_001095227.1), these helices are 15 Å apart, leading to the creation of a deep nonpolar cavity. In the `closed' structure (YP_749275.1), the helices partially unfold and rearrange, occluding the cavity and decreasing the solvent-exposed hydrophobic surface. These two complementary structures are reminiscent of the conformational switch in CRAL-TRIO carriers of hydrophobic compounds. It is suggested that both proteins may associate with the lipid bilayer in their `open' monomeric state by inserting their amphiphilic helices, α2 and α3, into the lipid bilayer. These bacterial proteins may function as carriers of nonpolar substances or as interfacially activated enzymes.

Structural insights into maize viviparous14, a key enzyme in the biosynthesis of the phytohormone abscisic acid

The key regulatory step in the biosynthesis of abscisic acid (ABA), a hormone central to the regulation of several important processes in plants, is the oxidative cleavage of the 11,12 double bond of a 9-cis-epoxycarotenoid. The enzyme viviparous14 (VP14) performs this cleavage in maize (Zea mays), making it a target for the rational design of novel chemical agents and genetic modifications that improve plant behavior through the modulation of ABA levels. The structure of VP14, determined to 3.2-Å resolution, provides both insight into the determinants of regio- and stereospecificity of this enzyme and suggests a possible mechanism for oxidative cleavage. Furthermore, mutagenesis of the distantly related CCD1 of maize shows how the VP14 structure represents a template for all plant carotenoid cleavage dioxygenases (CCDs). In addition, the structure suggests how VP14 associates with the membrane as a way of gaining access to its membrane soluble substrate.

The membrane environment modulates self-association of the human GpA TM domain--implications for membrane protein folding and transmembrane signaling

The influence of lipid bilayer properties on a defined and sequence-specific transmembrane helix-helix interaction is not well characterized yet. To study the potential impact of changing bilayer properties on a sequence-specific transmembrane helix-helix interaction, we have traced the association of fluorescent-labeled glycophorin A transmembrane peptides by fluorescence spectroscopy in model membranes with varying lipid compositions. The observed changes of the glycophorin A dimerization propensities in different lipid bilayers suggest that the lipid bilayer thickness severely influences the monomer-dimer equilibrium of this transmembrane domain, and dimerization was most efficient under hydrophobic matching conditions. Moreover, cholesterol considerably promotes self-association of transmembrane helices in model membranes by affecting the lipid acyl chain ordering. In general, the order of the lipid acyl chains appears to be an important factor involved in determining the strength and stability of transmembrane helix-helix interactions. As discussed, the described influences of membrane properties on transmembrane helix-helix interactions are highly important for understanding the mechanism of transmembrane protein folding and functioning as well as for gaining a deeper insight into the regulation of signal transduction via membrane integral proteins by bilayer properties.

Structure of the cytosolic portion of the motor protein prestin and functional role of the STAS domain in SLC26/SulP anion transporters

Prestin is the motor protein responsible for the somatic electromotility of cochlear outer hair cells and is essential for normal hearing sensitivity and frequency selectivity of mammals. Prestin is a member of mammalian solute-linked carrier 26 (SLC26) anion exchangers, a family of membrane proteins capable of transporting a wide variety of monovalent and divalent anions. SLC26 transporters play important roles in normal human physiology in different tissues, and many of them are involved in genetic diseases. SLC26 and related SulP transporters carry a hydrophobic membrane core and a C-terminal cytosolic portion that is essential in plasma membrane targeting and protein function. This C-terminal portion is mainly composed of a STAS (sulfate transporters and anti-sigma factor antagonist) domain, whose name is due to a remote but significant sequence similarity with bacterial ASA (anti-sigma factor antagonist) proteins. Here we present the crystal structure at 1.57 A resolution of the cytosolic portion of prestin, the first structure of a SulP transporter STAS domain, and its characterization in solution by heteronuclear multidimensional NMR spectroscopy. Prestin STAS significantly deviates from the related bacterial ASA proteins, especially in the N-terminal region, which-although previously considered merely as a generic linker between the domain and the last transmembrane helix-is indeed fully part of the domain. Hence, unexpectedly, our data reveal that the STAS domain starts immediately after the last transmembrane segment and lies beneath the lipid bilayer. A structure-function analysis suggests that this model can be a general template for most SLC26 and SulP anion transporters and supports the notion that STAS domains are involved in functionally important intramolecular and intermolecular interactions. Mapping of disease-associated or functionally harmful mutations on STAS structure indicates that they can be divided into two categories: those causing significant misfolding of the domain and those altering its interaction properties.

Molecular Model of Plasma PAF Acetylhydrolase-Lipoprotein Association: Insights from the Structure

Plasma platelet-activating factor acetylhydrolase (PAF-AH), also called lipoprotein-associated phospholipase A₂ (Lp-PLA₂), is a group VIIA PLA₂ enzyme that catalyzes the hydrolysis of PAF and certain oxidized phospholipids. Although the role of PAF-AH as a pro- or anti-atherosclerotic enzyme is highly debated, several studies have shown it to be an independent marker of cardiovascular diseases. In humans the majority of plasma PAF-AH is bound to LDL and a smaller portion to HDL; the majority of the enzyme being associated with small dense LDL and VHDL-1 subclasses. Several studies suggest that the anti- or pro-atherosclerotic tendency of PAF-AH might be dependent on the type of lipoprotein it is associated with. Amino acid residues in PAF-AH necessary for binding to LDL and HDL have been identified. However our understanding of the interaction of PAF-AH with LDL and HDL is still incomplete. In this review we present an overview of what is already known about the interaction of PAF-AH with lipoprotein particles, and we pose questions that are yet to be answered. The recently solved crystal structure of PAF-AH, along with functional work done by others is used as a guide to develop a model of interaction of PAF-AH with lipoprotein particles.

Contribution of charged and polar residues for the formation of the E1-E2 heterodimer from Hepatitis C Virus

The transmembrane domains of the envelope glycoprotein E1 and E2 have crucial multifunctional roles in the biogenesis of hepatitis C virus. We have performed molecular dynamics simulations to investigate a structural model of the transmembrane segments of the E1-E2 heterodimer. The simulations support the key role of the Lys370-Asp728 ion pair for mediating the E1-E2 heterodimerization. In comparison to these two residues, the simulation results also reveal the differential effect of the conserved Arg730 residue that has been observed in experimental studies. Furthermore, we discovered the formation of inter-helical hydrogen bonds via Asn367 that stabilize dimer formation. Simulations of single and double mutants further demonstrate the importance of the ion-pair and polar interactions between the interacting helix monomers. The conformation of the E1 fragment in the simulation of the E1-E2 heterodimer is in close agreement with an NMR structure of the E1 transmembrane segment. The proposed model of the E1-E2 heterodimer supports the postulated cooperative insertion of both helices by the translocon complex into the bilayer.

Structural determination and tryptophan fluorescence of heterokaryon incompatibility C2 protein (HET-C2), a fungal glycolipid transfer protein (GLTP), provide novel insights into glycolipid specificity and membrane interaction by the GLTP fold

HET-C2 is a fungal protein that transfers glycosphingolipids between membranes and has limited sequence homology with human glycolipid transfer protein (GLTP). The human GLTP fold is unique among lipid binding/transfer proteins, defining the GLTP superfamily. Herein, GLTP fold formation by HET-C2, its glycolipid transfer specificity, and the functional role(s) of its two Trp residues have been investigated. X-ray diffraction (1.9 A) revealed a GLTP fold with all key sugar headgroup recognition residues (Asp(66), Asn(70), Lys(73), Trp(109), and His(147)) conserved and properly oriented for glycolipid binding. Far-UV CD showed secondary structure dominated by alpha-helices and a cooperative thermal unfolding transition of 49 degrees C, features consistent with a GLTP fold. Environmentally induced optical activity of Trp/Tyr/Phe (2:4:12) detected by near-UV CD was unaffected by membranes containing glycolipid but was slightly altered by membranes lacking glycolipid. Trp fluorescence was maximal at approximately 355 nm and accessible to aqueous quenchers, indicating free exposure to the aqueous milieu and consistent with surface localization of the two Trps. Interaction with membranes lacking glycolipid triggered significant decreases in Trp emission intensity but lesser than decreases induced by membranes containing glycolipid. Binding of glycolipid (confirmed by electrospray injection mass spectrometry) resulted in a blue-shifted emission wavelength maximum (approximately 6 nm) permitting determination of binding affinities. The unique positioning of Trp(208) at the HET-C2 C terminus revealed membrane-induced conformational changes that precede glycolipid uptake, whereas key differences in residues of the sugar headgroup recognition center accounted for altered glycolipid specificity and suggested evolutionary adaptation for the simpler glycosphingolipid compositions of filamentous fungi.

Critical review of general guidelines for membrane proteins model building and analysis

Membrane proteins play major roles in many biological processes such as signalling, transport, etc. They have been shown to be involved in the development of many diseases and have become important drug targets per se. The understanding of their functional properties may be facilitated if a 3D structure is available. However, in the case of membrane proteins, only a few 3D structures have been solved to date. Bioinformatics and molecular modelling approaches are thus powerful alternatives to fill the gap between the sequence and the structure. Here, a review of the most recent approaches is proposed together with guidelines on how to use them. In addition, insofar as important biological processes require conformational changes, we discuss some interesting methods aimed at exploring the dynamic behaviour of proteins in their membrane environment. The paper ends with a brief description of useful approaches for determining oligomerisation or ligand binding sites.

Topology of 4-Pyridoxic acid dehydrogenase in transformed Escherichia coli cells

The topology of 4-pyridoxic acid dehydrogenase in the Escherichia coli cell membrane was examined with transformed E. coli cells overexpressing the enzyme from Mesorhizobium loti. The recombinant enzymes with a His(6)-tag either in the N-terminal region or at the C-terminus were localized on the E. coli cell membrane like the wild-type enzyme without a His(6)-tag. The His(6)-tags were labelled with Ni-NTA AP conjugate only when the E. coli protoplast cells were broken. The membrane-bound enzyme in the intact protoplast cells was not digested by trypsin, although the one in the gently broken protoplast cells was almost totally digested. Thus, 4-pyridoxic acid dehydrogenase was an integral monotopic protein, protruding into a cytoplasm side from the bacterial membrane. The deletion or mutation of a deduced transmembrane segment in 4-pyridoxic acid dehydrogenase made it an inclusion body, and the enzyme protein was not found in the E. coli cell membrane. Thus, it was suggested that the intact deduced transmembrane segment was necessary for 4-pyridoxic acid dehydrogenase to be localized on the bacterial cell membrane.

Molecular dynamics simulation of the heterodimeric mGluR2/5HT(2A) complex. An atomistic resolution study of a potential new target in psychiatric conditions

Homo- and heterodimerization is becoming an assessed concept in G-protein coupled receptor (GPCR) pharmacology, and the notion that GPCRs may dimerize or oligomerize is allowing for a reinterpretation of some inconsistencies or anomalies and is providing medicinal chemists with potentially relevant novel molecular targets for a variety of therapeutic conditions. Recently, it has been reported that two unrelated GPCRs, namely class C metabotropic glutamate receptor type-2 (mGluR2) and class A 5HT(2A) serotoninergic receptor, can heterodimerize at the transmembrane domain level. We performed a 40 ns molecular dynamics simulation of the mGluR2/5HT(2A) heterocomplex constructed around a TM4/TM5 interface and embedded in an explicit phospholipidic bilayer surrounded by water molecules. In a separate experiment, the monomeric 5HT(2A) receptor was simulated for additional 40 ns under the same conditions. The analysis and the comparison of the two simulations allowed us to clearly identify a cross-talk between the two protomers and to put forward an effect of the heterodimerization on the shape of the binding pocket of 5HT(2A). This result provides the first molecular explanation for the reported allosteric effect of mGluR2 on 5HT(2A)-mediated response and suggests that the heterocomplex can be a more suitable target for in silico screening than the monomeric protomers.

Effect of bilayer phospholipid composition and curvature on ligand transfer by the alpha-tocopherol transfer protein

We report here our preliminary investigations on the mechanism of alpha-TTP-mediated ligand transfer as assessed using fluorescence resonance energy transfer (FRET) assays. These assays monitor the movement of the model alpha-tocopherol fluorescent derivative ((R)-2,5,7,8-tetramethyl-chroman-2-[9-(7-nitro-benzo[1,2,5]oxadiazol-4-yl amino)-nonyl]-chroman-6-ol; NBD-Toc) from protein to acceptor vesicles containing the fluorescence quencher TRITC-PE. We have found that alpha-TTP utilizes a collisional mechanism of ligand transfer requiring direct protein-membrane contact, that rates of ligand transfer are greater to more highly curved lipid vesicles, and that such rates are insensitive to the presence of anionic phospholipids in the acceptor membrane. These results point to hydrophobic features of alpha-TTP dominating the binding energy between protein and membrane.

A comparative study of the COX-1 and COX-2 isozymes bound to lipid membranes

The monotopic proteins COX-1 and -2 in dimeric form bound to lipid bilayer membranes are studied using molecular dynamics simulations within an aqueous environment. The 25-ns simulations are performed for both isozymes with arachidonic acid bound in the cyclooxygenase sites. The interactions between the enzymes and the lipids are analyzed, providing insight into the attachment mechanism of monotopic proteins to membranes. Our study reveals some key differences between the two isozymes that include the orientations at which they sit on the surface of the membranes and the depths to which they embed within the membranes. The differences in membrane association of the isozymes indicate that they may integrate distinctively with the same membrane, and/or with different membranes or their lipid components. Our results indicate that arachidonic acid can be bound in the cyclooxygenase active site in distinct catalytically competent conformations that lead to certain hydroperoxy acids; and the arachidonic acid and/or cyclooxygenase sites undergo a conformational change which makes only one subunit of each homodimer catalytically active.

Interaction of monotopic membrane enzymes with a lipid bilayer: a coarse-grained MD simulation study

Monotopic membrane proteins bind tightly to cell membranes but do not generally span the lipid bilayer. Their interactions with lipid bilayers may be studied via coarse-grained molecular dynamics (CG-MD) simulations. Understanding such interactions is important as monotopic enzymes frequently act on hydrophobic substrates, while X-ray structures rarely provide direct information about their interactions with membranes. CG-MD self-assembly simulations enable prediction of the orientation and depth of insertion into a lipid bilayer of a monotopic protein, and also of the interactions of individual protein residues with lipid molecules. The CG-MD method has been evaluated via comparison with extended (>30 ns) atomistic simulations of monoamine oxidase, revealing good agreement between the results of coarse-grained and atomistic simulations. CG-MD simulations have been applied to a set of 11 monotopic proteins for which three-dimensional structures are available. These proteins may be divided into two groups on the basis of the results of the simulations. One group consists of those proteins which are inserted into the lipid bilayer to a limited extent, interacting mainly at the phospholipid-water interface. The second group consists of those which are inserted more deeply into the bilayer. Those monotopic proteins which are inserted more deeply cause significant local perturbation of bilayer properties such as bilayer thickness. Deeper insertion seems to correlate with a greater number of basic residues in the "foot" whereby a monotopic protein interacts with the membrane.

Determining the orientation of protegrin-1 in DLPC bilayers using an implicit solvent-membrane model

Continuum models that describe the effects of solvent and biological membrane molecules on the structure and behavior of antimicrobial peptides, holds a promise to improve our understanding of the mechanisms of antimicrobial action of these peptides. In such methods, a lipid bilayer model membrane is implicitly represented by multiple layers of relatively low dielectric constant embedded in a high dielectric aqueous solvent, while an antimicrobial peptide is accounted for by a dielectric cavity with fixed partial charge at the center of each one of its atoms. In the present work, we investigate the ability of continuum approaches to predict the most probable orientation of the β-hairpin antimicrobial peptide Protegrin-1 (PG-1) in DLPC lipid bilayers by calculating the difference in the transfer free energy from an aqueous environment to a membrane-water environment for multiple orientations. The transfer free energy is computed as a sum of two terms; polar/electrostatic and non-polar. They both include energetic and entropic contributions to the free energy. We numerically solve the Poisson-Boltzmann equation to calculate the electrostatic contribution to the transfer free energy, while the non-polar contribution to the free energy is approximated using a linear solvent accessible surface area relationships. The most probable orientation of PG-1 is that with the lowest relative transfer free energy. Our simulation results indicate that PG-1 assumes an oblique orientation in DLPC lipid bilayers. The predicted most favorable orientation was with a tilt angle of 19°, which is in qualitative agreement with the experimentally observed orientations derived from solid-state NMR data.

Finding and characterizing tunnels in macromolecules with application to ion channels and pores

We describe a new algorithm, CHUNNEL, to automatically find, characterize, and display tunnels or pores in proteins. The correctness and accuracy of the algorithm is verified on a constructed set of proteins and used to analyze large sets of real proteins. The verification set contains proteins with artificially created pores of known path and width profile. The previous benchmark algorithm, HOLE, is compared with the new algorithm. Results show that the major advantage of the new algorithm is that it can successfully find and characterize tunnels with no a priori guidance or clues about the location of the tunnel mouth, and it will successfully find multiple tunnels if present. CHUNNEL can also be used in conjunction with HOLE, with the former used toprime HOLE and the latter to track and characterize the pores. Analysis was conducted on families of membrane protein structures culled from the Protein Data Bank as well as on a set of transmembrane proteins with predicted membrane-aqueous phase interfaces, yielding the first completely automated examination of tunnels through membrane proteins, including tunnels that exit in the membrane bilayer.

Lipid-dependent membrane protein topogenesis

The topology of polytopic membrane proteins is determined by topogenic sequences in the protein, protein-translocon interactions, and interactions during folding within the protein and between the protein and the lipid environment. Orientation of transmembrane domains is dependent on membrane phospholipid composition during initial assembly as well as on changes in lipid composition postassembly. The membrane translocation potential of negative amino acids working in opposition to the positive-inside rule is largely dampened by the normal presence of phosphatidylethanolamine, thus explaining the dominance of positive residues as retention signals. Phosphatidylethanolamine provides the appropriate charge density that permits the membrane surface to maintain a charge balance between membrane translocation and retention signals and also allows the presence of negative residues in the cytoplasmic face of proteins for other purposes.

Crystal structure of human plasma platelet-activating factor acetylhydrolase: structural implication to lipoprotein binding and catalysis

Human plasma platelet-activating factor (PAF) acetylhydrolase functions by reducing PAF levels as a general anti-inflammatory scavenger and is linked to anaphylactic shock, asthma, and allergic reactions. The enzyme has also been implicated in hydrolytic activities of other pro-inflammatory agents, such as sn-2 oxidatively fragmented phospholipids. This plasma enzyme is tightly bound to low and high density lipoprotein particles and is also referred to as lipoprotein-associated phospholipase A2. The crystal structure of this enzyme has been solved from x-ray diffraction data collected to a resolution of 1.5 angstroms. It has a classic lipase alpha/beta-hydrolase fold, and it contains a catalytic triad of Ser273, His351, and Asp296. Two clusters of hydrophobic residues define the probable interface-binding region, and a prediction is given of how the enzyme is bound to lipoproteins. Additionally, an acidic patch of 10 carboxylate residues and a neighboring basic patch of three residues are suggested to play a role in high density lipoprotein/low density lipoprotein partitioning. A crystal structure is also presented of PAF acetylhydrolase reacted with the organophosphate compound paraoxon via its active site Ser273. The resulting diethyl phosphoryl complex was used to model the tetrahedral intermediate of the substrate PAF to the active site. The model of interface binding begins to explain the known specificity of lipoprotein-bound substrates and how the active site can be both close to the hydrophobic-hydrophilic interface and at the same time be accessible to the aqueous phase.

To flip or not to flip: lipid-protein charge interactions are a determinant of final membrane protein topology

The molecular details of how lipids influence final topological organization of membrane proteins are not well understood. Here, we present evidence that final topology is influenced by lipid-protein interactions most likely outside of the translocon. The N-terminal half of Escherichia coli lactose permease (LacY) is inverted with respect to the C-terminal half and the membrane bilayer when assembled in mutants lacking phosphatidylethanolamine and containing only negatively charged phospholipids. We demonstrate that inversion is dependent on interactions between the net charge of the cytoplasmic surface of the N-terminal bundle and the negative charge density of the membrane bilayer surface. A transmembrane domain, acting as a molecular hinge between the two halves of the protein, must also exit from the membrane for inversion to occur. Phosphatidylethanolamine dampens the translocation potential of negative residues in favor of the cytoplasmic retention potential of positive residues, thus explaining the dominance of positive over negative amino acids as co- or post-translational topological determinants.

Beta-barrel proteins that reside in the Escherichia coli outer membrane in vivo demonstrate varied folding behavior in vitro

Little is known about the dynamic process of membrane protein folding, and few models exist to explore it. In this study we doubled the number of Escherichia coli outer membrane proteins (OMPs) for which folding into lipid bilayers has been systematically investigated. We cloned, expressed, and folded nine OMPs: outer membrane protein X (OmpX), OmpW, OmpA, the crcA gene product (PagP), OmpT, outer membrane phospholipase A (OmpLa), the fadl gene product (FadL), the yaet gene product (Omp85), and OmpF. These proteins fold into the same bilayer in vivo and share a transmembrane beta-barrel motif but vary in sequence and barrel size. We quantified the ability of these OMPs to fold into a matrix of bilayer environments. Several trends emerged from these experiments: higher pH values, thinner bilayers, and increased bilayer curvature promote folding of all OMPs. Increasing the incubation temperature promoted folding of several OMPs but inhibited folding of others. We discovered that OMPs do not have the same ability to fold into any single bilayer environment. This suggests that although environmental factors influence folding, OMPs also have intrinsic qualities that profoundly modulate their folding. To rationalize the differences in folding efficiency, we performed kinetic and thermal denaturation experiments, the results of which demonstrated that OMPs employ different strategies to achieve the observed folding efficiency.

Remote homology detection of integral membrane proteins using conserved sequence features

Compared with globular proteins, transmembrane proteins are surrounded by a more intricate environment and, consequently, amino acid composition varies between the different compartments. Existing algorithms for homology detection are generally developed with globular proteins in mind and may not be optimal to detect distant homology between transmembrane proteins. Here, we introduce a new profile-profile based alignment method for remote homology detection of transmembrane proteins in a hidden Markov model framework that takes advantage of the sequence constraints placed by the hydrophobic interior of the membrane. We expect that, for distant membrane protein homologs, even if the sequences have diverged too far to be recognized, the hydrophobicity pattern and the transmembrane topology are better conserved. By using this information in parallel with sequence information, we show that both sensitivity and specificity can be substantially improved for remote homology detection in two independent test sets. In addition, we show that alignment quality can be improved for the most distant homologs in a public dataset of membrane protein structures. Applying the method to the Pfam domain database, we are able to suggest new putative evolutionary relationships for a few relatively uncharacterized protein domain families, of which several are confirmed by other methods. The method is called Searcher for Homology Relationships of Integral Membrane Proteins (SHRIMP) and is available for download at http://www.sbc.su.se/shrimp/.

The interaction of phospholipase A2 with a phospholipid bilayer: coarse-grained molecular dynamics simulations

A number of membrane-active enzymes act in a complex environment formed by the interface between a lipid bilayer and bulk water. Although x-ray diffraction studies yield structures of isolated enzyme molecules, a detailed characterization of their interactions with the interface requires a measure of how deeply such a membrane-associated protein penetrates into a lipid bilayer. Here, we apply coarse-grained (CG) molecular dynamics (MD) simulations to probe the interaction of porcine pancreatic phospholipase A2 (PLA2) with a lipid bilayer containing palmitoyl-oleoyl-phosphatidyl choline and palmitoyl-oleoyl-phosphatidyl glycerol molecules. We also used a configuration from a CG-MD trajectory to initiate two atomistic (AT) MD simulations. The results of the CG and AT simulations are evaluated by comparison with available experimental data. The membrane-binding surface of PLA2 consists of a patch of hydrophobic residues surrounded by polar and basic residues. We show this proposed footprint interacts preferentially with the anionic headgroups of the palmitoyl-oleoyl-phosphatidyl glycerol molecules. Thus, both electrostatic and hydrophobic interactions determine the location of PLA2 relative to the bilayer. From a general perspective, this study demonstrates that CG-MD simulations may be used to reveal the orientation and location of a membrane-surface-bound protein relative to a lipid bilayer, which may subsequently be refined by AT-MD simulations to probe more detailed interactions.

Orientation of the Escherichia coli outer membrane protein OmpX in phospholipid bilayer membranes determined by solid-State NMR

The solid-state NMR orientation-dependent frequencies measured for membrane proteins in macroscopically oriented lipid bilayers provide precise orientation restraints for structure determination in membranes. Here we show that this information can also be used to supplement crystallographic structural data to establish the orientation of a membrane protein in the membrane. This is achieved by incorporating a few orientation restraints, measured for the Escherichia coli outer membrane protein OmpX in magnetically oriented lipid bilayers (bicelles), in a simulated annealing calculation with the coordinates of the OmpX crystal structure. The (1)H-(15)N dipolar couplings measured for the seven Phe residues of OmpX in oriented bilayers can be assigned by back-calculation of the NMR spectrum from the crystal structure and are sufficient to establish the three-dimensional orientation of the protein in the membrane, while the (15)N chemical shifts provide a measure of cross-validation for the analysis. In C14 lipid bilayers, OmpX adopts a transmembrane orientation with a 7 degrees tilt of its beta-barrel axis relative to the membrane normal, matching the hydrophobic thickness of the barrel with that of the membrane.

Energetics of hydrophobic matching in lipid-protein interactions

Lipid chain length modulates the activity of transmembrane proteins by mismatch between the hydrophobic span of the protein and that of the lipid membrane. Relative binding affinities of lipids with different chain lengths are used to estimate the excess free energy of lipid-protein interaction that arises from hydrophobic mismatch. For a wide range of integral proteins and peptides, the energy cost is much less than the elastic penalty of fully stretching or compressing the lipid chains to achieve complete hydrophobic matching. The chain length dependences of the free energies of lipid association are described by a model that combines elastic chain extension with a free energy term that depends linearly on the extent of residual mismatch. The excess free energy densities involved lie in the region of 0.5-2.0 k(B)T x nm(-2). Values of this size could arise from exposure of hydrophobic groups to polar portions of the lipid or protein, but not directly to water, or alternatively from changes in tilt of the transmembrane helices that are energetically comparable to those activating mechanosensitive channels. The influence of hydrophobic mismatch on dimerization of transmembrane helices and their transfer between lipid vesicles, and on shifts in chain-melting transitions of lipid bilayers by incorporated proteins, is analyzed by using the same thermodynamic model. Segmental order parameters confirm that elastic lipid chain distortions are insufficient to compensate fully for the mismatch, but the dependence on chain length with tryptophan-anchored peptides requires that the free energy density of hydrophobic mismatch should increase with increasing extent of mismatch.

Incorporation of outer membrane protein OmpG in lipid membranes: protein-lipid interactions and beta-barrel orientation

OmpG is an intermediate size, monomeric, outer membrane protein from Escherichia coli, with n beta = 14 beta-strands. It has a large pore that is amenable to modification by protein engineering. The stoichiometry ( N b = 20) and selectivity ( K r = 0.7-1.2) of lipid-protein interaction with OmpG incorporated in dimyristoyl phosphatidylcholine bilayer membranes was determined with various 14-position spin-labeled lipids by using EPR spectroscopy. The limited selectivity for different lipid species is consistent with the disposition of charged residues in the protein. The conformation and orientation (beta-strand tilt and beta-barrel order parameters) of OmpG in disaturated phosphatidylcholines of odd and even chain lengths from C(12:0) to C(17:0) was determined from polarized infrared spectroscopy of the amide I and amide II bands. A discontinuity in the protein orientation (deduced from the beta-barrel order parameters) is observed at the point of hydrophobic matching of the protein with lipid chain length. Compared with smaller (OmpA; n beta = 8) and larger (FhuA; n beta = 22) monomeric E. coli outer membrane proteins, the stoichiometry of motionally restricted lipids increases linearly with the number of beta-strands, the tilt (beta approximately 44 degrees ) of the beta-strands is comparable for the three proteins, and the order parameter of the beta-barrel increases regularly with n beta. These systematic features of the integration of monomeric beta-barrel proteins in lipid membranes could be useful for characterizing outer membrane proteins of unknown structure.

In silico local structure approach: a case study on outer membrane proteins

The detection of Outer Membrane Proteins (OMP) in whole genomes is an actual question, their sequence characteristics have thus been intensively studied. This class of protein displays a common beta-barrel architecture, formed by adjacent antiparallel strands. However, due to the lack of available structures, few structural studies have been made on this class of proteins. Here we propose a novel OMP local structure investigation, based on a structural alphabet approach, i.e., the decomposition of 3D structures using a library of four-residue protein fragments. The optimal decomposition of structures using hidden Markov model results in a specific structural alphabet of 20 fragments, six of them dedicated to the decomposition of beta-strands. This optimal alphabet, called SA20-OMP, is analyzed in details, in terms of local structures and transitions between fragments. It highlights a particular and strong organization of beta-strands as series of regular canonical structural fragments. The comparison with alphabets learned on globular structures indicates that the internal organization of OMP structures is more constrained than in globular structures. The analysis of OMP structures using SA20-OMP reveals some recurrent structural patterns. The preferred location of fragments in the distinct regions of the membrane is investigated. The study of pairwise specificity of fragments reveals that some contacts between structural fragments in beta-sheets are clearly favored whereas others are avoided. This contact specificity is stronger in OMP than in globular structures. Moreover, SA20-OMP also captured sequential information. This can be integrated in a scoring function for structural model ranking with very promising results.

Coarse-grained MD simulations of membrane protein-bilayer self-assembly

Complete determination of a membrane protein structure requires knowledge of the protein position within the lipid bilayer. As the number of determined structures of membrane proteins increases so does the need for computational methods which predict their position in the lipid bilayer. Here we present a coarse-grained molecular dynamics approach to lipid bilayer self-assembly around membrane proteins. We demonstrate that this method can be used to predict accurately the protein position in the bilayer for membrane proteins with a range of different sizes and architectures.

Coarse-grained simulation: a high-throughput computational approach to membrane proteins

An understanding of the interactions of membrane proteins with a lipid bilayer environment is central to relating their structure to their function and stability. A high-throughput approach to prediction of membrane protein interactions with a lipid bilayer based on coarse-grained Molecular Dynamics simulations is described. This method has been used to develop a database of CG simulations (coarse-grained simulations) of membrane proteins (http://sbcb.bioch.ox.ac.uk/cgdb). Comparison of CG simulations and AT simulations (atomistic simulations) of lactose permease reveals good agreement between the two methods in terms of predicted lipid headgroup contacts. Both CG and AT simulations predict considerable local bilayer deformation by the voltage sensor domain of the potassium channel KvAP.

Implicit membrane models for membrane protein simulation

Implicit models of membrane environments offer computational advantages in simulations of membrane-interacting proteins and peptides. Such methods are especially useful for studies of long time scale processes, such as folding and aggregation, or very large complexes that are otherwise intractable with explicit lipid environments. Implicit models replace explicit solute-solvent interactions with a mean-field approach. In the most physical models, continuum dielectric electrostatics is combined with empirical formulations for the nonpolar components of the free energy of solvation. The practical use of a number of implicit membrane models ranging from the empirical IMM1 method to generalized Born-based methods with two-dielectric and multidielectric representations of biological membrane characteristics is presented.

Straightening and sequential buckling of the pore-lining helices define the gating cycle of MscS

We describe a mechanism connecting the adaptive behavior of the bacterial mechanosensitive channel MscS to the flexibility of the pore-lining helix TM3. Simulated expansion of the channel structure revealed straightening of a characteristic kink near Gly113 in the open state; return to the closed state produced an alternative kink at Gly121. Patch-clamp experiments showed that higher helical propensity introduced by a G113A mutation prevented inactivation. A similar mutation, G121A, kinetically impeded both closure and inactivation. Duplicating the glycines at each of these sites to increase flexibility produced directly opposite effects. The severely toxic G113A G121A mutation resulted in channels that could not inactivate or close with the release of tension. These data suggest that the open MscS features straight TM3 helices, which act as collapsible 'struts'. Closure and desensitization rely on buckling at Gly121, whereas the crystal-like kink at Gly113 is a feature of the inactivated state.

Characterization of the resting MscS: modeling and analysis of the closed bacterial mechanosensitive channel of small conductance

Channels from the MscS family are adaptive tension-activated osmolyte release valves that regulate turgor in prokaryotes and volume in plant chloroplasts. The crystal structure of Escherichia coli MscS has provided a starting point for detailed descriptions of its mechanism. However, solved in the absence of the lipid bilayer, this structure may deviate from a native conformation. In this study, we utilized molecular dynamics simulations and a new iterative extrapolated-motion protocol to pack the splayed peripheral TM1 and TM2 transmembrane helices along the central TM3 shaft. This modification restored the tension transmission route between the membrane and the channel gate. We also modeled the structure of the 26-amino acid N-terminal segments that were unresolved in the crystals. The resulting compact conformation, which we believe approximates the closed resting state of MscS, matches the hydrophobic thickness of the lipid bilayer with arginines 46, 54, and 74 facing the polar lipid headgroups. The pore-lining helices in this resting state feature alternative kinks near the conserved G121 instead of the G113 kinks observed in the crystal structure and the transmembrane barrel remains stable in extended molecular dynamics simulations. Further analysis of the dynamics of the pore constriction revealed several moderately asymmetric and largely dehydrated states. Biochemical and patch-clamp experiments with engineered double-cysteine mutants demonstrated cross-linking between predicted adjacent residue pairs, which formed either spontaneously or under moderate oxidation. The L72C-V99C bridge linking more peripheral TM2 to TM3 caused a shift of channel activation to higher pressures. TM3 to TM3 cross-links through the A84C-T93C, S95C-I97C, and A106C-G108C cysteine pairs were shown to lock MscS in a nonconductive state. Normal channel activity in these mutants could be recovered upon disulfide reduction with dithiothreitol. These results confirmed our modeling predictions of a closed MscS channel featuring a TM3 barrel that largely resembles the crystal conformation though with more tightly packed peripheral helices. From this closed-resting conformation, the TM3 helices must expand to allow for channel opening.

TOPDB: topology data bank of transmembrane proteins

The Topology Data Bank of Transmembrane Proteins (TOPDB) is the most complete and comprehensive collection of transmembrane protein datasets containing experimentally derived topology information currently available. It contains information gathered from the literature and from public databases available on the internet for more than a thousand transmembrane proteins. TOPDB collects details of various experiments that were carried out to learn about the topology of particular transmembrane proteins. In addition to experimental data from the literature, an extensive collection of structural data was also compiled from PDB and from PDBTM. Because topology information is often incomplete, for each protein in the database the most probable topology that is consistent with the collected experimental constraints was also calculated using the HMMTOP transmembrane topology prediction algorithm. Each record in TOPDB also contains information on the given protein sequence, name, organism and cross references to various other databases. The web interface of TOPDB includes tools for searching, relational querying and data browsing as well as for visualization. TOPDB is designed to bridge the gap between the number of transmembrane proteins available in sequence databases and the publicly accessible topology information of experimentally or computationally studied transmembrane proteins. TOPDB is available at http://topdb.enzim.hu.

Role of Aromatic Side Chains in the Folding and Thermodynamic Stability of Integral Membrane Proteins

Aromatic residues are frequently found in helical and beta-barrel integral membrane proteins enriched at the membrane-water interface. Although the importance of these residues in membrane protein folding has been rationalized by thermodynamic partition measurements using peptide model systems, their contribution to the stability of bona fide membrane proteins has never been demonstrated. Here, we have investigated the contribution of interfacial aromatic residues to the thermodynamic stability of the beta-barrel outer membrane protein OmpA from Escherichia coli in lipid bilayers by performing extensive mutagenesis and equilibrium folding experiments. Isolated interfacial tryptophanes contribute -2.0 kcal/mol, isolated interfacial tyrosines contribute -2.6 kcal/mol, and isolated interfacial phenylalanines contribute -1.0 kcal/mol to the stability of this protein. These values agree well with the prediction from the Wimley-White interfacial hydrophobicity scale, except for tyrosine residues, which contribute more than has been expected from the peptide models. Double mutant cycle analysis reveals that interactions between aromatic side chains become significant when their centroids are separated by less than 6 A but are nearly insignificant above 7 A. Aromatic-aromatic side chain interactions are on the order of -1.0 to -1.4 kcal/mol and do not appear to depend on the type of aromatic residue. These results suggest that the clustering of aromatic side chains at membrane interfaces provides an additional heretofore not yet recognized driving force for the folding and stability of integral membrane proteins.

The role of hydrophobic interactions in positioning of peripheral proteins in membranes

BACKGROUND

Three-dimensional (3D) structures of numerous peripheral membrane proteins have been determined. Biological activity, stability, and conformations of these proteins depend on their spatial positions with respect to the lipid bilayer. However, these positions are usually undetermined.

RESULTS

We report the first large-scale computational study of monotopic/peripheral proteins with known 3D structures. The optimal translational and rotational positions of 476 proteins are determined by minimizing energy of protein transfer from water to the lipid bilayer, which is approximated by a hydrocarbon slab with a decadiene-like polarity and interfacial regions characterized by water-permeation profiles. Predicted membrane-binding sites, protein tilt angles and membrane penetration depths are consistent with spin-labeling, chemical modification, fluorescence, NMR, mutagenesis, and other experimental studies of 53 peripheral proteins and peptides. Experimental membrane binding affinities of peripheral proteins were reproduced in cases that did not involve a helix-coil transition, specific binding of lipids, or a predominantly electrostatic association. Coordinates of all examined peripheral proteins and peptides with the calculated hydrophobic membrane boundaries, subcellular localization, topology, structural classification, and experimental references are available through the Orientations of Proteins in Membranes (OPM) database.

CONCLUSION

Positions of diverse peripheral proteins and peptides in the lipid bilayer can be accurately predicted using their 3D structures that represent a proper membrane-bound conformation and oligomeric state, and have membrane binding elements present. The success of the implicit solvation model suggests that hydrophobic interactions are usually sufficient to determine the spatial position of a protein in the membrane, even when electrostatic interactions or specific binding of lipids are substantial. Our results demonstrate that most peripheral proteins not only interact with the membrane surface, but penetrate through the interfacial region and reach the hydrocarbon interior, which is consistent with published experimental studies.

Carnitine palmitoyltransferase 2: analysis of membrane association and complex structure with a substrate analog

The mitochondrial membrane-associated carnitine palmitoyltransferase system is a validated target for the treatment of type 2 diabetes mellitus. To further facilitate structure-based drug discovery, we determined the crystal structure of rat CPT-2 (rCPT-2) in complex with the substrate analogue palmitoyl-aminocarnitine at 1.8A resolution. Biochemical analyses revealed a strong effect of this compound on rCPT-2 activity and stability. Using a computational approach we examined the membrane association of rCPT-2. The protein interacts with the membrane as a functional monomer and the calculations confirm the presence of a membrane association domain that consists of layers of hydrophobic and positively charged residues.

Glycolipid transfer proteins

Glycolipid transfer proteins (GLTPs) are small (24 kDa), soluble, ubiquitous proteins characterized by their ability to accelerate the intermembrane transfer of glycolipids in vitro. GLTP specificity encompasses both sphingoid- and glycerol-based glycolipids, but with a strict requirement that the initial sugar residue be beta-linked to the hydrophobic lipid backbone. The 3D architecture of GLTP reveals liganded structures with unique lipid-binding modes. The biochemical properties of GLTP action at the membrane surface have been studied rather comprehensively, but the biological role of GLTP remains enigmatic. What is clear is that GLTP differs distinctly from other known glycolipid-binding proteins, such as nonspecific lipid transfer proteins, lysosomal sphingolipid activator proteins, lectins, lung surfactant proteins as well as other lipid-binding/transfer proteins. Based on the unique conformational architecture that targets GLTP to membranes and enables glycolipid binding, GLTP is now considered the prototypical and founding member of a new protein superfamily in eukaryotes.

On the derivation of propensity scales for predicting exposed transmembrane residues of helical membrane proteins

Helical membrane proteins (HMPs) play a crucial role in diverse physiological processes. Given the difficulty in determining their structures by experimental techniques, it is desired to develop computational methods for predicting the burial status of transmembrane residues. Deriving a propensity scale for the 20 amino acids to be exposed to the lipid bilayer from known structures is central to developing such methods. A fundamental problem in this regard is what would be the optimal way of deriving propensity scales. Here, we show that this problem can be reformulated such that an optimal scale is straightforwardly obtained in an analytical fashion. The derived scale favorably compares with others in terms of both algorithmic optimality and practical prediction accuracy. It also allows interesting insights into the structural organization of HMPs. Furthermore, the presented approach can be applied to other bioinformatics problems of HMPs, too. All the data sets and programs used in the study and detailed primary results are available upon request.

TOAC spin labels in the backbone of alamethicin: EPR studies in lipid membranes

Alamethicin is a 19-amino-acid residue hydrophobic peptide that produces voltage-dependent ion channels in membranes. Analogues of the Glu(OMe)(7,18,19) variant of alamethicin F50/5 that are rigidly spin-labeled in the peptide backbone have been synthesized by replacing residue 1, 8, or 16 with 2,2,6,6-tetramethyl-piperidine-1-oxyl-4-amino-4-carboxyl (TOAC), a helicogenic nitroxyl amino acid. Conventional electron paramagnetic resonance spectra are used to determine the insertion and orientation of the TOAC(n) alamethicins in fluid lipid bilayer membranes of dimyristoyl phosphatidylcholine. Isotropic (14)N-hyperfine couplings indicate that TOAC(8) and TOAC(16) are situated in the hydrophobic core of the membrane, whereas the TOAC(1) label resides closer to the membrane surface. Anisotropic hyperfine splittings show that alamethicin is highly ordered in the fluid membranes. Experiments with aligned membranes demonstrate that the principal diffusion axis lies close to the membrane normal, corresponding to a transmembrane orientation. Combination of data from the three spin-labeled positions yields both the dynamic order parameter of the peptide backbone and the intramolecular orientations of the TOAC groups. The latter are compared with x-ray diffraction results from alamethicin crystals. Saturation transfer electron paramagnetic resonance, which is sensitive to microsecond rotational motion, reveals that overall rotation of alamethicin is fast in fluid membranes, with effective correlation times <30 ns. Thus, alamethicin does not form large stable aggregates in fluid membranes, and ionic conductance must arise from transient or voltage-induced associations.