Structural adaptations of proteins to different biological membranes

Archaeal lipids

The major archaeal membrane glycerolipids are distinguished from those of bacteria and eukaryotes by the contrasting stereochemistry of their glycerol backbones, and by the use of ether-linked isoprenoid-based alkyl chains rather than ester-linked fatty acyl chains for their hydrophobic moieties. These fascinating compounds play important roles in the extremophile lifestyles of many species, but are also present in the growing numbers of recently discovered mesophilic archaea. The past decade has witnessed significant advances in our understanding of archaea in general and their lipids in particular. Much of the new information has come from the ability to screen large microbial populations via environmental metagenomics, which has revolutionised our understanding of the extent of archaeal biodiversity that is coupled with a strict conservation of their membrane lipid compositions. Significant additional progress has come from new culturing and analytical techniques that are gradually enabling archaeal physiology and biochemistry to be studied in real time. These studies are beginning to shed light on the much-discussed and still-controversial process of eukaryogenesis, which probably involved both bacterial and archaeal progenitors. Puzzlingly, although eukaryotes retain many attributes of their putative archaeal ancestors, their lipid compositions only reflect their bacterial progenitors. Finally, elucidation of archaeal lipids and their metabolic pathways have revealed potentially interesting applications that have opened up new frontiers for biotechnological exploitation of these organisms. This review is concerned with the analysis, structure, function, evolution and biotechnology of archaeal lipids and their associated metabolic pathways.

Supramolecular organization and dynamics of mannosylated phosphatidylinositol lipids in the mycobacterial plasma membrane

Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis (TB), a disease that claims ~1.6 million lives annually. The current treatment regime is long and expensive, and missed doses contribute to drug resistance. Therefore, development of new anti-TB drugs remains one of the highest public health priorities. Mtb has evolved a complex cell envelope that represents a formidable barrier to antibiotics. The Mtb cell envelop consists of four distinct layers enriched for Mtb specific lipids and glycans. Although the outer membrane, comprised of mycolic acid esters, has been extensively studied, less is known about the plasma membrane, which also plays a critical role in impacting antibiotic efficacy. The Mtb plasma membrane has a unique lipid composition, with mannosylated phosphatidylinositol lipids (phosphatidyl-myoinositol mannosides, PIMs) comprising more than 50% of the lipids. However, the role of PIMs in the structure and function of the membrane remains elusive. Here, we used multiscale molecular dynamics (MD) simulations to understand the structure-function relationship of the PIM lipid family and decipher how they self-organize to shape the biophysical properties of mycobacterial plasma membranes. We assess both symmetric and asymmetric assemblies of the Mtb plasma membrane and compare this with residue distributions of Mtb integral membrane protein structures. To further validate the model, we tested known anti-TB drugs and demonstrated that our models agree with experimental results. Thus, our work sheds new light on the organization of the mycobacterial plasma membrane. This paves the way for future studies on antibiotic development and understanding Mtb membrane protein function.

Hysteretic Pressure Dependence of Ca2+ Binding in LH1 Bacterial Membrane Chromoproteins

Much of the thermodynamic parameter values that support life are set by the properties of proteins. While the denaturing effects of pressure and temperature on proteins are well documented, their precise structural nature is rarely revealed. This work investigates the destabilization of multiple Ca²⁺ binding sites in the cyclic LH1 light-harvesting membrane chromoprotein complexes from two Ca-containing sulfur purple bacteria by hydrostatic high-pressure perturbation spectroscopy. The native (Ca-saturated) and denatured (Ca-depleted) phases of these complexes are well distinguishable by much-shifted bacteriochlorophyll a exciton absorption bands serving as innate optical probes in this study. The pressure-induced denaturation of the complexes related to the failure of the protein Ca-binding pockets and the concomitant breakage of hydrogen bonds between the pigment chromophores and protein environment were found cooperative, involving all or most of the Ca²⁺ binding sites, but irreversible. The strong hysteresis observed in the spectral and kinetic characteristics of phase transitions along the compression and decompression pathways implies asymmetry in the relevant free energy landscapes and activation free energy distributions. A phase transition pressure equal to about 1.9 kbar was evaluated for the complexes from Thiorhodovibrio strain 970 from the pressure dependence of biphasic kinetics observed in the minutes to 100 h time range.

Length matters: Functional flip of the short TatA transmembrane helix

The twin arginine translocase (Tat) exports folded proteins across bacterial membranes. The putative pore-forming or membrane-weakening component (TatA_d in B. subtilis) is anchored to the lipid bilayer via an unusually short transmembrane α-helix (TMH), with less than 16 residues. Its tilt angle in different membranes was analyzed under hydrophobic mismatch conditions, using synchrotron radiation circular dichroism and solid-state NMR. Positive mismatch (introduced either by reconstitution in short-chain lipids or by extending the hydrophobic TMH length) increased the helix tilt of the TMH as expected. Negative mismatch (introduced either by reconstitution in long-chain lipids or by shortening the TMH), on the other hand, led to protein aggregation. These data suggest that the TMH of TatA is just about long enough for stable membrane insertion. At the same time, its short length is a crucial factor for successful translocation, as demonstrated here in native membrane vesicles using an in vitro translocation assay. Furthermore, when reconstituted in model membranes with negative spontaneous curvature, the TMH was found to be aligned parallel to the membrane surface. This intrinsic ability of TatA to flip out of the membrane core thus seems to play a key role in its membrane-destabilizing effect during Tat-dependent translocation.

Using Alphafold2 to Predict the Structure of the Gp5/M Dimer of Porcine Respiratory and Reproductive Syndrome Virus

Porcine reproductive and respiratory syndrome virus is a positive-stranded RNA virus of the family Arteriviridae. The Gp5/M dimer, the major component of the viral envelope, is required for virus budding and is an antibody target. We used alphafold2, an artificial-intelligence-based system, to predict a credible structure of Gp5/M. The short disulfide-linked ectodomains lie flat on the membrane, with the exception of the erected N-terminal helix of Gp5, which contains the antibody epitopes and a hypervariable region with a changing number of carbohydrates. The core of the dimer consists of six curved and tilted transmembrane helices, and three are from each protein. The third transmembrane regions extend into the cytoplasm as amphiphilic helices containing the acylation sites. The endodomains of Gp5 and M are composed of seven β-strands from each protein, which interact via β-strand seven. The area under the membrane forms an open cavity with a positive surface charge. The M and Orf3a proteins of coronaviruses have a similar structure, suggesting that all four proteins are derived from the same ancestral gene. Orf3a, like Gp5/M, is acylated at membrane-proximal cysteines. The role of Gp5/M during virus replication is discussed, in particular the mechanisms of virus budding and models of antibody-dependent virus neutralization.

Comparative Molecular Dynamics Simulation Studies of Realistic Eukaryotic, Prokaryotic, and Archaeal Membranes

We present a comparative all-atom molecular dynamics simulation study of 18 biomembrane systems with lipid compositions corresponding to eukaryotic, bacterial, and archaebacterial membranes together with three single-component lipid bilayers. A total of 105 lipid types used in this study include diverse sterols and glycerol-based lipids with acyl chains of various lengths, unsaturation degrees, and branched or cyclic moieties. Our comparative analysis provides deeper insight into the influences of sterols and lipid unsaturation on the structural and mechanical properties of these biomembranes, including water permeation into the membrane hydrocarbon core. For sterol-containing membranes, sterol fraction is correlated with the membrane thickness, the area compressibility modulus, and lipid order but anticorrelated with the area per lipid and sterol tilt angles. Similarly, for all 18 biomembranes, lipid order is correlated with the membrane thickness and area compressibility modulus. Sterols and lipid unsaturation produce opposite effects on membrane thickness, but only sterols influence water permeation into the membrane. All membrane systems are accessible for public use in CHARMM-GUI Archive. They can be used as templates to expedite future modeling of realistic cell membranes with transmembrane and peripheral membrane proteins to study their structure, dynamics, molecular interactions, and function in a nativelike membrane environment.

Spatial arrangement of proteins in planar and curved membranes by PPM 3.0

Cellular protrusions, invaginations, and many intracellular organelles have strongly curved membrane regions. Transmembrane and peripheral membrane proteins that induce, sense, or stabilize such regions cannot be properly fitted into a single flat bilayer. To treat such proteins, we developed a new method and a web tool, PPM 3.0, for positioning proteins in curved or planar, single or multiple membranes. This method determines the energetically optimal spatial position, the hydrophobic thickness, and the radius of intrinsic curvature of a membrane-deforming protein structure by arranging it in a single or several sphere-shaped or planar membrane sections. In addition, it can define the lipid-embedded regions of a protein that simultaneously spans several membranes or determine the optimal position of a peptide in a spherical micelle. The PPM 3.0 web server operates with 17 types of biological membranes and 4 types of artificial bilayers. It is publicly available at https://opm.phar.umich.edu/ppm_server3. PPM 3.0 was applied to identify and characterize arrangements in membranes of 128 proteins with a significant intrinsic curvature, such as BAR domains, annexins, Piezo, and MscS mechanosensitive channels, cation-chloride cotransporters, as well as mitochondrial ATP synthases, calcium uniporters, and TOM complexes. These proteins form large complexes that are mainly localized in mitochondria, plasma membranes, and endosomes. Structures of bacterial drug efflux pumps, AcrAB-TolC, MexAB-OrpM, and MacAB-TolC, were positioned in both membranes of the bacterial cell envelop, while structures of multimeric gap-junction channels were arranged in two opposed cellular membranes.

Hydrostatic High-Pressure-Induced Denaturation of LH2 Membrane Proteins

The denaturation of globular proteins by high pressure is frequently associated with the release of internal voids and/or the exposure of the hydrophobic protein interior to a polar aqueous solvent. Similar evidence with respect to membrane proteins is not available. Here, we investigate the impact of hydrostatic pressures reaching 12 kbar on light-harvesting 2 integral membrane complexes of purple photosynthetic bacteria using two types of innate chromophores in separate strategic locations: bacteriochlorophyll-a in the hydrophobic interior and tryptophan at both protein-solvent interfacial gateways to internal voids. The complexes from mutant Rhodobacter sphaeroides with low resilience against pressure were considered in parallel with the naturally robust complexes of Thermochromatium tepidum. In the former case, a firm correlation was established between the abrupt blue shift of the bacteriochlorophyll-a exciton absorption, a known indicator of the breakage of tertiary structure pigment-protein hydrogen bonds, and the quenching of tryptophan fluorescence, a supposed result of further protein solvation. No such effects were observed in the reference complex. While these data may be naively taken as supporting evidence of the governing role of hydration, the analysis of atomistic model structures of the complexes confirmed the critical part of the structure in the pressure-induced denaturation of the membrane proteins studied.

Identification and assessment of cardiolipin interactions with E. coli inner membrane proteins

Integral membrane proteins are localized and/or regulated by lipids present in the surrounding bilayer. While bacteria have relatively simple membranes, there is ample evidence that many bacterial proteins bind to specific lipids, especially the anionic lipid cardiolipin. Here, we apply molecular dynamics simulations to assess lipid binding to 42 different Escherichia coli inner membrane proteins. Our data reveal an asymmetry between the membrane leaflets, with increased anionic lipid binding to the inner leaflet regions of the proteins, particularly for cardiolipin. From our simulations, we identify >700 independent cardiolipin binding sites, allowing us to identify the molecular basis of a prototypical cardiolipin binding site, which we validate against structures of bacterial proteins bound to cardiolipin. This allows us to construct a set of metrics for defining a high-affinity cardiolipin binding site on bacterial membrane proteins, paving the way for a heuristic approach to defining other protein-lipid interactions.

The Det.Belt Server: A Tool to Visualize and Estimate Amphipathic Solvent Belts around Membrane Proteins

Detergents wrap around membrane proteins to form a belt covering the hydrophobic part of the protein serving for membrane insertion and interaction with lipids. The number of detergent monomers forming this belt is usually unknown to investigators, unless dedicated detergent quantification is undertaken, which for many projects is difficult to setup. Yet, having an approximate knowledge of the amount of detergent forming the belt is extremely useful, to better grasp the protein of interest in interaction with its direct environment rather than picturing the membrane protein "naked". We created the Det.Belt server to dress up membrane proteins and represent in 3D the bulk made by detergent molecules wrapping in a belt. Many detergents are included in a database, allowing investigators to screen in silico the effect of different detergents around their membrane protein. The input number of detergents is changeable with fast recomputation of the belt for interactive usage. Metrics representing the belt are readily available together with scripts to render quality 3D images for publication. The Det.Belt server is a tool for biochemists to better grasp their sample.

Thermodynamics-Based Molecular Modeling of α-Helices in Membranes and Micelles

The Folding of Membrane-Associated Peptides (FMAP) method was developed for modeling α-helix formation by linear peptides in micelles and lipid bilayers. FMAP 2.0 identifies locations of α-helices in the amino acid sequence, generates their three-dimensional models in planar bilayers or spherical micelles, and estimates their thermodynamic stabilities and tilt angles, depending on temperature and pH. The method was tested for 723 peptides (926 data points) experimentally studied in different environments and for 170 single-pass transmembrane (TM) proteins with available crystal structures. FMAP 2.0 detected more than 95% of experimentally observed α-helices with an average error in helix end determination of around 2, 3, 4, and 5 residues per helix for peptides in water, micelles, bilayers, and TM proteins, respectively. Helical and nonhelical residue states were predicted with an accuracy from 0.86 to 0.96, and the Matthews correlation coefficient was from 0.64 to 0.88 depending on the environment. Experimental micelle- and membrane-binding energies and tilt angles of peptides were reproduced with a root-mean-square deviation of around 2 kcal/mol and 7°, respectively. The TM and non-TM states of hydrophobic and pH-triggered α-helical peptides in various lipid bilayers were reproduced in more than 95% of cases. The FMAP 2.0 web server (https://membranome.org/fmap) is publicly available to explore the structural polymorphism of antimicrobial, cell-penetrating, fusion, and other membrane-binding peptides, which is important for understanding the mechanisms of their biological activities.

Biological Functions and Identification of Novel Biomarker Expressed on the Surface of Breast Cancer-Derived Cancer Stem Cells via Proteomic Analysis

Breast cancer is one of the most common life-threatening malignancies and the top cause of cancer deaths in women. Although many conventional therapies exist for its treatment, breast cancer still has many handicaps to overcome. Cancer stem cells (CSCs) are a well-known cause of tumor recurrences due to the ability of CSCs for self-renewal and differentiation into cell subpopulations, similar to stem cells. To fully treat breast cancer, a strategy for the treatment of both cancer cells and CSCs is required. However, current strategies for the eradication of CSCs are non-specific and have low efficacy. Therefore, surface biomarkers to selectively treat CSCs need to be developed. Here, 34 out of 641 surface biomarkers on CSCs were identified by proteomic analysis between the human breast adenocarcinoma cell line MCF-7 and MCF-7-derived CSCs. Among them, carcinoembryonic antigen-related cell adhesion molecules 6 (CEACAM6 or CD66c), a member of the CEA family, was selected as a novel biomarker on the CSC surface. This biomarker was then experimentally validated and evaluated for use as a CSC-specific marker. Its biological effects were assessed by treating breast cancer stem cells (BCSCs) with short hairpin (sh)-RNA under oxidative cellular conditions. This study is the first to evaluate the biological function of CD66c as a novel biomarker on the surface of CSCs. This marker is available as a moiety for use in the development of targeted therapeutic agents against CSCs.

MemSTATS: A Benchmark Set of Membrane Protein Symmetries and Pseudosymmetries

In membrane proteins, symmetry and pseudosymmetry often have functional or evolutionary implications. However, available symmetry detection methods have not been tested systematically on this class of proteins because of the lack of an appropriate benchmark set. Here we present MemSTATS, a publicly available benchmark set of both quaternary- and internal-symmetries in membrane protein structures. The symmetries are described in terms of order, repeated elements, and orientation of the axis with respect to the membrane plane. Moreover, using MemSTATS, we compare the performance of four widely used symmetry detection algorithms and highlight specific challenges and areas for improvement in the future.

Cardiolipin-Dependent Properties of Model Mitochondrial Membranes from Molecular Simulations

Cardiolipin is an anionic lipid found in the mitochondrial membranes of eukaryotes ranging from unicellular microorganisms to metazoans. This unique lipid contributes to various mitochondrial functions, including metabolism, mitochondrial membrane fusion and/or fission dynamics, and apoptosis. However, differences in cardiolipin content between the two mitochondrial membranes, as well as dynamic fluctuations in cardiolipin content in response to stimuli and cellular signaling events, raise questions about how cardiolipin concentration affects mitochondrial membrane structure and dynamics. Although cardiolipin’s structural and dynamic roles have been extensively studied in binary mixtures with other phospholipids, the biophysical properties of cardiolipin in higher number lipid mixtures are still not well resolved. Here, we used molecular dynamics simulations to investigate the cardiolipin-dependent properties of ternary lipid bilayer systems that mimic the major components of mitochondrial membranes. We found that changes to cardiolipin concentration only resulted in minor changes to bilayer structural features but that the lipid diffusion was significantly affected by those alterations. We also found that cardiolipin position along the bilayer surfaces correlated to negative curvature deflections, consistent with the induction of negative curvature stress in the membrane monolayers. This work contributes to a foundational understanding of the role of cardiolipin in altering the properties in ternary lipid mixtures composed of the major mitochondrial phospholipids, providing much-needed insights to help understand how cardiolipin concentration modulates the biophysical properties of mitochondrial membranes.

PerMM: A Web Tool and Database for Analysis of Passive Membrane Permeability and Translocation Pathways of Bioactive Molecules

The PerMM web server and database were developed for quantitative analysis and visualization of passive translocation of bioactive molecules across lipid membranes. The server is the first physics-based web tool that calculates membrane binding energies and permeability coefficients of diverse molecules through artificial and natural membranes (phospholipid bilayers, PAMPA-DS, blood-brain barrier, and Caco-2/MDCK cell membranes). It also visualizes the transmembrane translocation pathway as a sequence of translational and rotational positions of a permeant as it moves across the lipid bilayer, along with the corresponding changes in solvation energy. The server can be applied for prediction of permeability coefficients of compounds with diverse chemical scaffolds to facilitate selection and optimization of potential drug leads. The complementary PerMM database allows comparison of computationally and experimentally determined permeability coefficients for more than 500 compounds in different membrane systems. The website and database are freely accessible at https://permm.phar.umich.edu/ .

Physics-Based Method for Modeling Passive Membrane Permeability and Translocation Pathways of Bioactive Molecules

Assessment of permeability is a critical step in the drug development process for selection of drug candidates with favorable ADME properties. We have developed a novel physics-based method for fast computational modeling of passive permeation of diverse classes of molecules across lipid membranes. The method is based on heterogeneous solubility-diffusion theory and operates with all-atom 3D structures of solutes and the anisotropic solvent model of the lipid bilayer characterized by transbilayer profiles of dielectric and hydrogen bonding capacity parameters. The optimal translocation pathway of a solute is determined by moving an ensemble of representative conformations of the molecule through the dioleoyl-phosphatidylcholine (DOPC) bilayer and optimizing their rotational orientations in every point of the transmembrane trajectory. The method calculates (1) the membrane-bound state of the solute molecule; (2) free energy profile of the solute along the permeation pathway; and (3) the permeability coefficient obtained by integration over the transbilayer energy profile and assuming a constant size-dependent diffusivity along the membrane normal. The accuracy of the predictions was evaluated against experimental permeability coefficients measured in pure lipid membranes (for 78 compounds, R² was 0.88 and rmse was 1.15 log units), PAMPA-DS (for 280 compounds, R² was 0.75 and rmse was 1.59 log units), BBB (for 182 compounds, R² was 0.69 and rmse was 0.87 log units), and Caco-2/MDCK assays (for 165 compounds, R² was 0.52 and rmse was 0.89 log units).

Exploring the impact of proteins on the line tension of a phase-separating ternary lipid mixture

The separation of lipid mixtures into thermodynamically stable phase-separated domains is dependent on lipid composition, temperature, and system size. Using molecular dynamics simulations, the line tension between thermodynamically stable lipid domains formed from ternary mixtures of di-C16:0 PC:di-C18:2 PC:cholesterol at 40:40:20 mol. % ratio was investigated via two theoretical approaches. The line tension was found to be 3.1 ± 0.2 pN by capillary wave theory and 4.7 ± 3.7 pN by pressure tensor anisotropy approaches for coarse-grained models based on the Martini force field. Using an all-atom model of the lipid membrane based on the CHARMM36 force field, the line tension was found to be 3.6 ± 0.9 pN using capillary wave theory and 1.8 ± 2.2 pN using pressure anisotropy approaches. The discrepancy between estimates of the line tension based on capillary wave theory and pressure tensor anisotropy methods is discussed. Inclusion of protein in Martini membrane lipid mixtures was found to reduce the line tension by 25%-35% as calculated by the capillary wave theory approach. To further understand and predict the behavior of proteins in phase-separated membranes, we have formulated an analytical Flory-Huggins model and parameterized it against the simulation results. Taken together these results suggest a general role for proteins in reducing the thermodynamic cost associated with domain formation in lipid mixtures and quantifies the thermodynamic driving force promoting the association of proteins to domain interfaces.

Polyphosphoinositide-Binding Domains: Insights from Peripheral Membrane and Lipid-Transfer Proteins

Within eukaryotic cells, biochemical reactions need to be organized on the surface of membrane compartments that use distinct lipid constituents to dynamically modulate the functions of integral proteins or influence the selective recruitment of peripheral membrane effectors. As a result of these complex interactions, a variety of human pathologies can be traced back to improper communication between proteins and membrane surfaces; either due to mutations that directly alter protein structure or as a result of changes in membrane lipid composition. Among the known structural lipids found in cellular membranes, phosphatidylinositol (PtdIns) is unique in that it also serves as the membrane-anchored precursor of low-abundance regulatory lipids, the polyphosphoinositides (PPIn), which have restricted distributions within specific subcellular compartments. The ability of PPIn lipids to function as signaling platforms relies on both non-specific electrostatic interactions and the selective stereospecific recognition of PPIn headgroups by specialized protein folds. In this chapter, we will attempt to summarize the structural diversity of modular PPIn-interacting domains that facilitate the reversible recruitment and conformational regulation of peripheral membrane proteins. Outside of protein folds capable of capturing PPIn headgroups at the membrane interface, recent studies detailing the selective binding and bilayer extraction of PPIn species by unique functional domains within specific families of lipid-transfer proteins will also be highlighted. Overall, this overview will help to outline the fundamental physiochemical mechanisms that facilitate localized interactions between PPIn lipids and the wide-variety of PPIn-binding proteins that are essential for the coordinate regulation of cellular metabolism and membrane dynamics.

Au courant computation of the PDB to audit diffraction anisotropy of soluble and membrane proteins

This data article makes available the informed computation of the whole Protein Data Bank (PDB) to investigate diffraction anisotropy on a large scale and to perform statistics. This data has been investigated in detail in "X-ray diffraction reveals the intrinsic difference in the physical properties of membrane and soluble proteins" [1]. Diffraction anisotropy is traditionally associated with absence of contacts in-between macromolecules within the crystals in a given direction of space. There are however many case that do not follow this empirical rule. To investigate and sort out this discrepancy, we computed diffraction anisotropy for every entry of the PDB, and put them in context of relevant metrics to compare X-ray diffraction in reciprocal space to the crystal packing in real space. These metrics were either extracted from PDB files when available (resolution, space groups, cell parameters, solvent content), or calculated using standard procedures (anisotropy, crystal contacts, presence of ligands). More specifically, we separated entries to compare soluble vs membrane proteins, and further separated the later in subcategories according to their insertion in the membrane, function, or type of crystallization (Type I vs Type II crystal packing). This informed database is being made available to investigators in the raw and curated formats that can be re-used for further downstream studies. This dataset is useful to test ideas and to ascertain hypothesis based on statistical analysis.

Physics of HIV

This review summarizes over a decade of investigations into how membrane-binding proteins from the HIV-1 virus interact with lipid membrane mimics various HIV and host T-cell membranes. The goal of the work was to characterize at the molecular level both the elastic and structural changes that occur due to HIV protein/membrane interactions, which could lead to new drugs to thwart the HIV virus. The main technique used to study these interactions is diffuse X-ray scattering, which yields the bending modulus, K_C, as well as structural parameters such as membrane thickness, area/lipid and position of HIV peptides (parts of HIV proteins) in the membrane. Our methods also yield information about lipid chain order or disorder caused by the peptides. This review focuses on three stages of the HIV-1 life cycle: 1) infection, 2) Tat membrane transport, and 3) budding. In the infection stage, our lab studied three different parts of HIV-1 gp41 (glycoprotein 41 fusion protein): 1) FP23, the N-terminal 23 amino acids that interact non-specifically with the T-cell host membrane to cause fusion of two membranes, and its trimer version, 2) CRAC (cholesterol recognition amino acid consensus sequence), on the MPER (membrane proximal external region) near the membrane-spanning domain, and 3) LLP2 (lentiviral lytic peptide 2) on the CTT (cytoplasmic C-terminal tail). For Tat transport, we used membrane mimics of the T-cell nuclear membrane as well as simpler models that varied charge and negative curvature. For membrane budding, we varied the myristoylation of the MA₃₁ peptide as well as the negatively charged lipid. These studies show that HIV peptides with different roles in the HIV life cycle affect differently the relevant membrane mimics. In addition, the membrane lipid composition plays an important role in the peptides' effects.

X-ray diffraction reveals the intrinsic difference in the physical properties of membrane and soluble proteins

Membrane proteins are distinguished from soluble proteins by their insertion into biological membranes. This insertion is achieved via a noticeable arrangement of hydrophobic amino acids that are exposed at the surface of the protein, and renders the interaction with the aliphatic tails of lipids more energetically favorable. This important difference between these two categories of proteins is the source of the need for a specific handling of membrane proteins, which transpired in the creation of new tools for their recombinant expression, purification and even crystallization. Following this line, we show here that crystals of membrane proteins display systematically higher diffraction anisotropy than those of soluble proteins. This phenomenon dramatically hampers structure solution and refinement, and has a strong impact on the quality of electron-density maps. A farther search for origins of this phenomenon showed that the type of crystallization, and thus the crystal packing, has no impact on anisotropy, nor does the nature or function of the membrane protein. Membrane proteins fully embedded within the membrane display equal anisotropy compared to the ones with extra membranous domains or fusions with soluble proteins. Overall, these results overturn common beliefs and call for a specific handling of their diffraction data.

Sequences flanking the transmembrane segments facilitate mitochondrial localization and membrane fusion by mitofusin

Mitochondria constantly connect through membrane fusion. The merging of the outer mitochondrial membrane requires mitofusin (MFN) proteins. MFN is a membrane-anchored GTPase, but whether it is sufficient to achieve fusion, and if so how, is largely unknown. We have taken advantage of a similar GTPase named atlastin (ATL), which mediates fusion of the endoplasmic reticulum (ER), as its mechanism is better understood. Domain swapping experiments show that MFN is capable of fusing membranes, even on the ER. The C-terminal tail of MFN contains an amphipathic helix that promotes fusion. MFN is properly inserted into the mitochondrial membrane with the help of the helix and neighboring hydrophobic residues. These findings provide insight into how mitochondria fuse.

Mitochondria constantly divide and fuse. Homotypic fusion of the outer mitochondrial membranes requires the mitofusin (MFN) proteins, a family of dynamin-like GTPases. MFNs are anchored in the membrane by transmembrane (TM) segments, exposing both the N-terminal GTPase domain and the C-terminal tail (CT) to the cytosol. This arrangement is very similar to that of the atlastin (ATL) GTPases, which mediate fusion of endoplasmic reticulum (ER) membranes. We engineered various MFN-ATL chimeras to gain mechanistic insight into MFN-mediated fusion. When MFN1 is localized to the ER by TM swapping with ATL1, it functions in the maintenance of ER morphology and fusion. In addition, an amphipathic helix in the CT of MFN1 is exchangeable with that of ATL1 and critical for mitochondrial localization of MFN1. Furthermore, hydrophobic residues N-terminal to the TM segments of MFN1 play a role in membrane targeting but not fusion. Our findings provide important insight into MFN-mediated membrane fusion.

Charged residues next to transmembrane regions revisited: "Positive-inside rule" is complemented by the "negative inside depletion/outside enrichment rule"

Background

Transmembrane helices (TMHs) frequently occur amongst protein architectures as means for proteins to attach to or embed into biological membranes. Physical constraints such as the membrane’s hydrophobicity and electrostatic potential apply uniform requirements to TMHs and their flanking regions; consequently, they are mirrored in their sequence patterns (in addition to TMHs being a span of generally hydrophobic residues) on top of variations enforced by the specific protein’s biological functions.

Results

With statistics derived from a large body of protein sequences, we demonstrate that, in addition to the positive charge preference at the cytoplasmic inside (positive-inside rule), negatively charged residues preferentially occur or are even enriched at the non-cytoplasmic flank or, at least, they are suppressed at the cytoplasmic flank (negative-not-inside/negative-outside (NNI/NO) rule). As negative residues are generally rare within or near TMHs, the statistical significance is sensitive with regard to details of TMH alignment and residue frequency normalisation and also to dataset size; therefore, this trend was obscured in previous work. We observe variations amongst taxa as well as for organelles along the secretory pathway. The effect is most pronounced for TMHs from single-pass transmembrane (bitopic) proteins compared to those with multiple TMHs (polytopic proteins) and especially for the class of simple TMHs that evolved for the sole role as membrane anchors.

Conclusions

The charged-residue flank bias is only one of the TMH sequence features with a role in the anchorage mechanisms, others apparently being the leucine intra-helix propensity skew towards the cytoplasmic side, tryptophan flanking as well as the cysteine and tyrosine inside preference. These observations will stimulate new prediction methods for TMHs and protein topology from a sequence as well as new engineering designs for artificial membrane proteins.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0404-4) contains supplementary material, which is available to authorized users.

Comprehensive analysis of the numbers, lengths and amino acid compositions of transmembrane helices in prokaryotic, eukaryotic and viral integral membrane proteins of high-resolution structure

We report a comprehensive analysis of the numbers, lengths and amino acid compositions of transmembrane helices in 235 high-resolution structures of integral membrane proteins. The properties of 1551 transmembrane helices in the structures were compared with those obtained by analysis of the same amino acid sequences using topology prediction tools. Explanations for the 81 (5.2%) missing or additional transmembrane helices in the prediction results were identified. Main reasons for missing transmembrane helices were mis-identification of N-terminal signal peptides, breaks in α-helix conformation or charged residues in the middle of transmembrane helices and transmembrane helices with unusual amino acid composition. The main reason for additional transmembrane helices was mis-identification of amphipathic helices, extramembrane helices or hairpin re-entrant loops. Transmembrane helix length had an overall median of 24 residues and an average of 24.9 ± 7.0 residues and the most common length was 23 residues. The overall content of residues in transmembrane helices as a percentage of the full proteins had a median of 56.8% and an average of 55.7 ± 16.0%. Amino acid composition was analysed for the full proteins, transmembrane helices and extramembrane regions. Individual proteins or types of proteins with transmembrane helices containing extremes in contents of individual amino acids or combinations of amino acids with similar physicochemical properties were identified and linked to structure and/or function. In addition to overall median and average values, all results were analysed for proteins originating from different types of organism (prokaryotic, eukaryotic, viral) and for subgroups of receptors, channels, transporters and others.

Ligand and membrane-binding behavior of the phosphatidylinositol transfer proteins PITPα and PITPβ

Phosphatidylinositol transfer proteins (PITPs) are believed to be lipid transfer proteins because of their ability to transfer either phosphatidylinositol (PI) or phosphatidylcholine (PC) between membrane compartments, in vitro. However, the detailed mechanism of this transfer process is not fully established. To further understand the transfer mechanism of PITPs we examined the interaction of PITPs with membranes using dual polarization interferometry (DPI), which measures protein binding affinity on a flat immobilized lipid surface. In addition, a fluorescence resonance energy transfer (FRET)-based assay was also employed to monitor how quickly PITPs transfer their ligands to lipid vesicles. DPI analysis revealed that PITPβ had a higher affinity to membranes compared with PITPα. Furthermore, the FRET-based transfer assay revealed that PITPβ has a higher ligand transfer rate compared with PITPα. However, both PITPα and PITPβ demonstrated a preference for highly curved membrane surfaces during ligand transfer. In other words, ligand transfer rate was higher when the accepting vesicles were highly curved.

Studies on the interaction of 2-amino-3-hydroxy-anthraquinone with surfactant micelles reveal its nucleation in human MDA-MB-231 breast adinocarcinoma cells

Structural, spectroscopic and biochemical studies on 2-amino-3-hydroxy-anthraquinone (AQ) were carried out.

Membrane Protein Production in E. coli for Applications in Drug Discovery

Producing high quality purified membrane proteins for structure-based drug design and biophysical assays compatible with typical timelines in drug discovery is a significant challenge. Escherichia coli has been an expression host of the utmost importance for soluble proteins and has applications for membrane proteins as well. However, membrane protein overexpression in E. coli may lead to toxicity and low yields of functional product. Here, we review the challenges encountered with heterologous overproduction of α-helical membrane proteins in E. coli and a range of strategies to overcome them. A detailed protocol is also provided for expression and screening of membrane proteins in E. coli using a His-specific fluorescent probe and fluorescent size-exclusion chromatography.

Open Questions on the Origin of Eukaryotes

Despite recent progress, the origin of the eukaryotic cell remains enigmatic. It is now known that the last eukaryotic common ancestor was complex and that endosymbiosis played a crucial role in eukaryogenesis at least via the acquisition of the alphaproteobacterial ancestor of mitochondria. However, the nature of the mitochondrial host is controversial, although the recent discovery of an archaeal lineage phylogenetically close to eukaryotes reinforces models proposing archaea-derived hosts. We argue that, in addition to improved phylogenomic analyses with more comprehensive taxon sampling to pinpoint the closest prokaryotic relatives of eukaryotes, determining plausible mechanisms and selective forces at the origin of key eukaryotic features, such as the nucleus or the bacterial-like eukaryotic membrane system, is essential to constrain existing models.

Bacterial gene import and mesophilic adaptation in archaea

It is widely believed that the archaeal ancestor was hyperthermophilic, but during archaeal evolution, several lineages - including haloarchaea and their sister methanogens, the Thaumarchaeota, and the uncultured Marine Group II and Marine Group III Euryarchaeota (MGII/III) - independently adapted to lower temperatures. Recent phylogenomic studies suggest that the ancestors of these lineages were recipients of massive horizontal gene transfer from bacteria. Many of the acquired genes, which are often involved in metabolism and cell envelope biogenesis, were convergently acquired by distant mesophilic archaea. In this Opinion article, we explore the intriguing hypothesis that the import of these bacterial genes was crucial for the adaptation of archaea to mesophilic lifestyles.

MemProtMD: Automated Insertion of Membrane Protein Structures into Explicit Lipid Membranes

There has been exponential growth in the number of membrane protein structures determined. Nevertheless, these structures are usually resolved in the absence of their lipid environment. Coarse-grained molecular dynamics (CGMD) simulations enable insertion of membrane proteins into explicit models of lipid bilayers. We have automated the CGMD methodology, enabling membrane protein structures to be identified upon their release into the PDB and embedded into a membrane. The simulations are analyzed for protein-lipid interactions, identifying lipid binding sites, and revealing local bilayer deformations plus molecular access pathways within the membrane. The coarse-grained models of membrane protein/bilayer complexes are transformed to atomistic resolution for further analysis and simulation. Using this automated simulation pipeline, we have analyzed a number of recently determined membrane protein structures to predict their locations within a membrane, their lipid/protein interactions, and the functional implications of an enhanced understanding of the local membrane environment of each protein.

•

A simulation pipeline for predicting the location of a membrane protein in a bilayer

•

A protocol for identifying novel membrane protein structures in the PDB

•

Analysis of lipid binding sites and local bilayer deformation by membrane proteins

•

Functional implications from enhanced understanding of local membrane environments

Plant membrane assays with cytokinin receptors underpin the unique role of free cytokinin bases as biologically active ligands

Cytokinin receptors play a key role in cytokinin-dependent processes regulating plant growth, development, and adaptation; therefore, the functional properties of these receptors are of great importance. Previously the properties of cytokinin receptors were investigated in heterologous assay systems using unicellular microorganisms, mainly bacteria, expressing receptor proteins. However, within microorganisms receptors reside in an alien environment that might distort the receptor properties. Therefore, a new assay system has been developed allowing studies of individual receptors within plant membranes (i.e. closer to their natural environment). The main ligand-binding characteristics of receptors from Arabidopsis [AHK2, AHK3, and AHK4] and maize (ZmHK1) were refined in this new system, and the properties of full-length Arabidopsis receptor AHK2 were characterized for the first time. Ligand specificity profiles of receptors towards cytokinin bases were comparable with the profiles retrieved in bacterial assay systems. In contrast, cytokinin-9-ribosides displayed a strongly reduced affinity for receptors in the plant assay system, indicating that ribosides as the common transport form of cytokinins have no or very weak cytokinin activity. This underpins the central role of free bases as the sole biologically active cytokinin compounds. According to molecular modelling and docking studies, N (9)-ribosylation alters the bonding pattern in cytokinin-receptor interaction and prevents β6-β7 loop movement important for tight hormone binding. A common feature of all receptors was a greatly reduced ligand binding at low (5.0-5.5) pH. The particularly high sensitivity of ZmHK1 to pH changes leads to the suggestion that some cytokinin receptors may play an additional role as pH sensors in the lumen of the endoplasmic reticulum.

Water penetration profile at the protein-lipid interface in Na,K-ATPase membranes

The affinity of ionized fatty acids for the Na,K-ATPase is used to determine the transmembrane profile of water penetration at the protein-lipid interface. The standardized intensity of the electron spin echo envelope modulation (ESEEM) from (2)H-hyperfine interaction with D2O is determined for stearic acid, n-SASL, spin-labeled systematically at the C-n atoms throughout the chain. In both native Na,K-ATPase membranes from shark salt gland and bilayers of the extracted membrane lipids, the D2O-ESEEM intensities of fully charged n-SASL decrease progressively with position down the fatty acid chain toward the terminal methyl group. Whereas the D2O intensities decrease sharply at the n = 9 position in the lipid bilayers, a much broader transition region in the range n = 6 to 10 is found with Na,K-ATPase membranes. Correction for the bilayer population in the membranes yields the intrinsic D2O-intensity profile at the protein-lipid interface. For positions at either end of the chains, the D2O concentrations at the protein interface are greater than in the lipid bilayer, and the positional profile is much broader. This reveals the higher polarity, and consequently higher intramembrane water concentration, at the protein-lipid interface. In particular, there is a significant water concentration adjacent to the protein at the membrane midplane, unlike the situation in the bilayer regions of this cholesterol-rich membrane. Experiments with protonated fatty acid and phosphatidylcholine spin labels, both of which have a considerably lower affinity for the Na,K-ATPase, confirm these results.

Quantitative and Systems-Based Approaches for Deciphering Bacterial Membrane Interactome and Gene Function

High-throughput genomic and proteomic methods provide a concise description of the molecular constituents of a cell, whereas systems biology strives to understand the way these components function as a whole. Recent developments, such as genome editing technologies and protein epitope-tagging coupled with high-sensitivity mass-spectrometry, allow systemic studies to be performed at an unprecedented scale. Available methods can be successfully applied to various goals, both expanding fundamental knowledge and solving applied problems. In this review, we discuss the present state and future of bacterial cell envelope interactomics, with a specific focus on host-pathogen interactions and drug target discovery. Both experimental and computational methods will be outlined together with examples of their practical implementation.

Interfacial partitioning of a loop hinge residue contributes to diacylglycerol affinity of conserved region 1 domains

Conventional and novel isoenzymes of PKC are activated by the membrane-embedded second messenger diacylglycerol (DAG) through its interactions with the C1 regulatory domain. The affinity of C1 domains to DAG varies considerably among PKCs. To gain insight into the origin of differential DAG affinities, we conducted high-resolution NMR studies of C1B domain from PKCδ (C1Bδ) and its W252Y variant. The W252Y mutation was previously shown to render C1Bδ less responsive to DAG (Dries, D. R., Gallegos, L. L., and Newton, A. C. (2007) A single residue in the C1 domain sensitizes novel protein kinase C isoforms to cellular diacylglycerol production. J. Biol. Chem. 282, 826-830) and thereby emulate the behavior of C1B domains from conventional PKCs that have a conserved Tyr at the equivalent position. Our data revealed that W252Y mutation did not perturb the conformation of C1Bδ in solution but significantly reduced its propensity to partition into a membrane-mimicking environment in the absence of DAG. Using detergent micelles doped with a paramagnetic lipid, we determined that both the residue identity at position 252 and complexation with diacylglycerol influence the geometry of C1Bδ-micelle interactions. In addition, we identified the C-terminal helix α1 of C1Bδ as an interaction site with the head groups of phosphatidylserine, a known activator of PKCδ. Taken together, our studies (i) reveal the identities of C1Bδ residues involved in interactions with membrane-mimicking environment, DAG, and phosphatidylserine, as well as the affinities associated with each event and (ii) suggest that the initial ligand-independent membrane recruitment of C1B domains, which is greatly facilitated by the interfacial partitioning of Trp-252, is responsible, at least in part, for the differential DAG affinities.

Life at the border: adaptation of proteins to anisotropic membrane environment

This review discusses main features of transmembrane (TM) proteins which distinguish them from water-soluble proteins and allow their adaptation to the anisotropic membrane environment. We overview the structural limitations on membrane protein architecture, spatial arrangement of proteins in membranes and their intrinsic hydrophobic thickness, co-translational and post-translational folding and insertion into lipid bilayers, topogenesis, high propensity to form oligomers, and large-scale conformational transitions during membrane insertion and transport function. Special attention is paid to the polarity of TM protein surfaces described by profiles of dipolarity/polarizability and hydrogen-bonding capacity parameters that match polarity of the lipid environment. Analysis of distributions of Trp resides on surfaces of TM proteins from different biological membranes indicates that interfacial membrane regions with preferential accumulation of Trp indole rings correspond to the outer part of the lipid acyl chain region-between double bonds and carbonyl groups of lipids. These "midpolar" regions are not always symmetric in proteins from natural membranes. We also examined the hydrophobic effect that drives insertion of proteins into lipid bilayer and different free energy contributions to TM protein stability, including attractive van der Waals forces and hydrogen bonds, side-chain conformational entropy, the hydrophobic mismatch, membrane deformations, and specific protein-lipid binding.

E. coli outer membrane and interactions with OmpLA

The outer membrane of Gram-negative bacteria is a unique asymmetric lipid bilayer composed of phospholipids (PLs) in the inner leaflet and lipopolysaccharides (LPSs) in the outer leaflet. Its function as a selective barrier is crucial for the survival of bacteria in many distinct environments, and it also renders Gram-negative bacteria more resistant to antibiotics than their Gram-positive counterparts. Here, we report the structural properties of a model of the Escherichia coli outer membrane and its interaction with outer membrane phospholipase A (OmpLA) utilizing molecular dynamics simulations. Our results reveal that given the lipid composition used here, the hydrophobic thickness of the outer membrane is ∼3 Å thinner than the corresponding PL bilayer, mainly because of the thinner LPS leaflet. Further thinning in the vicinity of OmpLA is observed due to hydrophobic matching. The particular shape of the OmpLA barrel induces various interactions between LPS and PL leaflets, resulting in asymmetric thinning around the protein. The interaction between OmpLA extracellular loops and LPS (headgroups and core oligosaccharides) stabilizes the loop conformation with reduced dynamics, which leads to secondary structure variation and loop displacement compared to that in a DLPC bilayer. In addition, we demonstrate that the LPS/PL ratios in asymmetric bilayers can be reliably estimated by the per-lipid surface area of each lipid type, and there is no statistical difference in the overall membrane structure for the outer membranes with one more or less LPS in the outer leaflet, although individual lipid properties vary slightly.

Solid-state NMR studies of full-length BamA in lipid bilayers suggest limited overall POTRA mobility

The outer membrane protein BamA is the key player in β-barrel assembly in Gram-negative bacteria. Despite the availability of high-resolution crystal structures, the dynamic behavior of the transmembrane domain and the large periplasmic extension consisting of five POTRA (POlypeptide-TRansport-Associated) domains remains unclear. We demonstrate reconstitution of full-length BamA in proteoliposomes at low lipid-to-protein ratio, leading to high sensitivity and resolution in solid-state NMR (ssNMR) experiments. We detect POTRA domains in ssNMR experiments probing rigid protein segments in our preparations. These results suggest that the periplasmic region of BamA is firmly attached to the β-barrel and does not experience fast global motion around the angle between POTRA 2 and 3. We show that this behavior holds at lower protein concentrations and elevated temperatures. Chemical shift variations observed after reconstitution in lipids with different chain lengths and saturation levels are compatible with conformational plasticity of BamA's transmembrane domain. Electron microscopy of the ssNMR samples shows that BamA can cause local disruptions of the lipid bilayer in proteoliposomes. The observed interplay between protein-protein and protein-lipid interactions may be critical for BamA-mediated insertion of substrates into the outer membrane.

X-ray structure, thermodynamics, elastic properties and MD simulations of cardiolipin/dimyristoylphosphatidylcholine mixed membranes

Cardiolipins (CLs) are important biologically for their unique role in biomembranes that couple phosphorylation and electron transport like bacterial plasma membranes, chromatophores, chloroplasts and mitochondria. CLs are often tightly coupled to proteins involved in oxidative phosphorylation. The first step in understanding the interaction of CL with proteins is to obtain the pure CL structure, and the structure of mixtures of CL with other lipids. In this work we use a variety of techniques to characterize the fluid phase structure, material properties and thermodynamics of mixtures of dimyristoylphosphatidylcholine (DMPC) with tetramyristoylcardiolipin (TMCL), both with 14-carbon chains, at several mole percentages. X-ray diffuse scattering was used to determine structure, including bilayer thickness and area/lipid, the bending modulus, KC, and SXray, a measure of chain orientational order. Our results reveal that TMCL thickens DMPC bilayers at all mole percentages, with a total increase of ∼6 Å in pure TMCL, and increases AL from 64 Å(2) (DMPC at 35 °C) to 109 Å(2) (TMCL at 50 °C). KC increases by ∼50%, indicating that TMCL stiffens DMPC membranes. TMCL also orders DMPC chains by a factor of ∼2 for pure TMCL. Coarse grain molecular dynamics simulations confirm the experimental thickening of 2 Å for 20mol% TMCL and locate the TMCL headgroups near the glycerol-carbonyl region of DMPC; i.e., they are sequestered below the DMPC phosphocholine headgroup. Our results suggest that TMCL plays a role similar to cholesterol in that it thickens and stiffens DMPC membranes, orders chains, and is positioned under the umbrella of the PC headgroup. CL may be necessary for hydrophobic matching to inner mitochondrial membrane proteins. Differential scanning calorimetry, SXray and CGMD simulations all suggest that TMCL does not form domains within the DMPC bilayers. We also determined the gel phase structure of TMCL, which surprisingly displays diffuse X-ray scattering, like a fluid phase lipid. AL=40.8 Å(2) for the ½TMCL gel phase, smaller than the DMPC gel phase with AL=47.2 Å(2), but similar to AL of DLPE=41 Å(2), consistent with untilted chains in gel phase TMCL.

Single tryptophan and tyrosine comparisons in the N-terminal and C-terminal interface regions of transmembrane GWALP peptides

Hydrophobic membrane-spanning helices often are flanked by interfacial aromatic or charged residues. In this paper, we compare the consequences of single Trp → Tyr substitutions at each interface for the properties of a defined transmembrane helix in the absence of charged residues. The choice of molecular framework is critical for these single-residue experiments because the presence of "too many" aromatic residues (more than one at either membrane-water interface) introduces excess dynamic averaging of solid state NMR observables. To this end, we compare the outcomes when changing W(5) or W(19), or both of them, to tyrosine in the well-characterized transmembrane peptide acetyl-GGALW(5)(LA)6LW(19)LAGA-amide ("GWALP23"). By means of solid-state (2)H and (15)N NMR experiments, we find that Y(19)GW(5)ALP23 displays similar magnitudes of peptide helix tilt as Y(5)GW(19)ALP23 and responds similarly to changes in bilayer thickness, from DLPC to DMPC to DOPC. The presence of Y(19) changes the azimuthal rotation angle ρ (about the helix axis) to a similar extent as Y(5), but in the opposite direction. When tyrosines are substituted for both tryptophans to yield GY(5,19)ALP23, the helix tilt angle is again of comparable magnitude, and furthermore, the preferred azimuthal rotation angle ρ is relatively unchanged from that of GW(5,19)ALP23. The extent of dynamic averaging increases marginally when Tyr replaces Trp. Yet, importantly, all members of the peptide family having single Tyr or Trp residues near each interface exhibit only moderate and not highly extensive dynamic averaging. The results provide important benchmarks for evaluating conformational and dynamic control of membrane protein function.