Vertebrate membrane proteins: structure, function, and insights from biophysical approaches

Nuclear translocation of the membrane oxoeicosanoid/androgen receptor, OXER1: Possible mechanisms involved

OXER1, the receptor for the arachidonic acid metabolite 5-οxo-eicosatetraenoic acid (5-oxo-ETE), has been reported to also bind and mediate the membrane-initiated actions of androgens. Indeed, androgens antagonize the 5-oxo-ETE effects through OXER1, affecting a number of signaling pathways and inhibiting cancer cell proliferation and migration. OXER1, being a GPCR, was classically described to be localized in the plasma membrane. However, for numerous GPCRs, there is now strong evidence that they can be also found in other cellular compartments, including the nucleus. The aim of the present work was to investigate OXER1's possible localization in the nucleus and identify the mechanism(s) involved. For this purpose, we verified OXER1's nuclear presence by immunofluorescence and western blot, in whole cells and nuclei of two different prostate cancer cell lines (DU-145 and LNCaP) and in CHO cells transfected with a GFP labelled OXER1, both in untreated and OXER1 ligands' treated cells. Mutated, OXER1-tGFP expressing, CHO cells were used to verify that OXER1 agonist (5-oxo-ETE) binding is necessary for OXER1 nuclear translocation. NLS sequences were in silico identified, and a specific inhibitor, as well as, specific importins' siRNAs were also utilized to explore the mechanism involved. Moreover, we examined the role of palmitoylation in OXER1 nuclear translocation by in silico identifying possible palmitoylation sites and using a palmitoylation inhibitor. Our results clearly show that OXER1 can be localized in the nucleus, in an agonist-dependent manner, that is inhibited by androgens. We also provide evidence for two possible mechanisms for its nuclear trafficking, that involve receptor palmitoylation and importin-mediated cytoplasmic-nuclear transport. In our knowledge, it is the first time that a membrane androgen receptor is identified into the nucleus, suggesting an alternative, more direct, mode of action, involving nuclear mechanisms. Therefore, our findings provide new insights on androgen-mediated actions and androgen-lipid interactions, and reveal new possible therapeutic targets, not only for cancer, but also for other pathological conditions in which OXER1 may have an important role.

Characterizing the West Nile Virus's polyprotein from nucleotide sequence to protein structure - Computational tools

Objectives

West Nile virus (WNV) belongs to the Flaviviridae family and causes West Nile fever. The mechanism of transmission involves the culex mosquito species. Infected individuals are primarily asymptomatic, and few exhibit common symptoms. Moreover, 10 % of neuronal infection caused by this virus cause death. The proteins encoded by these genes had been uncharacterized, although understanding their function and structure is important for formulating antiviral drugs.

Methods

Herein, we used in silico approaches, including various bioinformatic tools and databases, to analyse the proteins from the WNV polyprotein individually. The characterization included GC content, physicochemical properties, conserved domains, soluble and transmembrane regions, signal localization, protein disorder, and secondary structure features and their respective 3D protein structures.

Results

Among 11 proteins, eight had >50 % GC content, eight proteins had basic pI values, three proteins were unstable under in vitro conditions, four were thermostable according to >100 AI values and some had negative GRAVY values in physicochemical analyses. All protein-conserved domains were shared among Flaviviridae family members. Five proteins were soluble and lacked transmembrane regions. Two proteins had signals for localization in the host endoplasmic reticulum. Non-structural (NS) 2A showed low protein disorder. The secondary structural features and tertiary structure models provide a valuable biochemical resource for designing selective substrates and synthetic inhibitors.

Conclusions

WNV proteins NS2A, NS2B, PM, NS3 and NS5 can be used as drug targets for the pharmacological design of lead antiviral compounds.

Sponge-phase Lipid Droplets as Synthetic Organelles: An Ultrafast Study of Hydrogen Bonding and Interfacial Environments

Bottom-up design of biomimetic organelles has gained recent attention as a route towards understanding the transition between non-living matter and life. Despite various artificial lipid membranes being developed, the specific relations between lipid structure, composition, interfacial properties, and morphology are not currently understood. Sponge-phase droplets contain dense, nonlamellar lipid bilayer networks that capture the complexities of the endoplasmic reticulum (ER), making them ideal artificial models of such organelles. Here, we combine ultrafast two-dimensional infrared (2D IR) spectroscopy and molecular dynamics simulations to investigate the interfacial H-bond networks in sponge-phase droplets composed of glycolipid and nonionic detergents. In the sponge phase, the interfacial environments are more hydrated and water molecules confined to the nanometer-scale aqueous channels in the sponge phase exhibit dynamics that are significantly slower compared to bulk water. Surfactant configurations and microscopic phase separation play a dominant role in determining membrane curvature and slow dynamics observed in the sponge phase. The studies suggest that H-bond networks within the nanometer-scale channels are disrupted not only by confinement but also by the interactions of surfactants, which extend 1-2 nm from the bilayer surface. The results provide a molecular-level description for controlling phase and morphology in the design of synthetic lipid organelles.

TMEM244 Is a Long Non-Coding RNA Necessary for CTCL Cell Growth

Transmembrane protein 244 (TMEM244) was annotated to be a member of the TMEM family, which are is a component of cell membranes and is involved in many cellular processes. To date, the expression of the TMEM244 protein has not been experimentally confirmed, and its function has not been clarified. Recently, the expression of the TMEM244 gene was acknowledged to be a diagnostic marker for Sézary syndrome, a rare cutaneous T-cell lymphoma (CTCL). In this study, we aimed to determine the role of the TMEM244 gene in CTCL cells. Two CTCL cell lines were transfected with shRNAs targeting the TMEM244 transcript. The phenotypic effect of TMEM244 knockdown was validated using green fluorescent protein (GFP) growth competition assays and AnnexinV/7AAD staining. Western blot analysis was performed to identify the TMEM244 protein. Our results indicate that TMEM244 is not a protein-coding gene but a long non-coding RNA (lncRNA) that is necessary for the growth of CTCL cells.

Mechanism of interactions between tripeptide NCW on cellular membrane using molecular dynamic simulation

Tripeptide NCW identified in Mizuhopecten yessoensis has been shown to possess in vivo antihypertensive effect. However, the poor understanding of the absorption of NCW across the membrane limits its application. In this study, we have investigated the interaction of NCW with DPPC membrane via 400 ns all-atom molecular dynamic simulation using GROMACS software. The structural variations of NCW during absorption, the location and distribution of NCW in the membrane, and the effect of NCW on the properties of membranes during simulation were analyzed to understand the dynamic behavior of NCW in DPPC membrane system. The results suggested that the structures of NCW were stable during simulation. Further, NCW could bind on the surface of the DPPC membrane and enter the hydrophobic interior of the DPPC membrane. Residue Try played an important role in the absorption of NCW across the membrane. Hydrogen bonds and hydrophobic interactions stabilized the interaction of NCW with the membrane. All the above studies analyzed the interaction mechanism between NCW and DPPC membranes at the atomic level and laid the foundation for further transmembrane studies of NCW.

Insights into Membrane Curvature Sensing and Membrane Remodeling by Intrinsically Disordered Proteins and Protein Regions

Cellular membranes are highly dynamic in shape. They can rapidly and precisely regulate their shape to perform various cellular functions. The protein’s ability to sense membrane curvature is essential in various biological events such as cell signaling and membrane trafficking. As they are bound, these curvature-sensing proteins may also change the local membrane shape by one or more curvature driving mechanisms. Established curvature-sensing/driving mechanisms rely on proteins with specific structural features such as amphipathic helices and intrinsically curved shapes. However, the recent discovery and characterization of many proteins have shattered the protein structure–function paradigm, believing that the protein functions require a unique structural feature. Typically, such structure-independent functions are carried either entirely by intrinsically disordered proteins or hybrid proteins containing disordered regions and structured domains. It is becoming more apparent that disordered proteins and regions can be potent sensors/inducers of membrane curvatures. In this article, we outline the basic features of disordered proteins and regions, the motifs in such proteins that encode the function, membrane remodeling by disordered proteins and regions, and assays that may be employed to investigate curvature sensing and generation by ordered/disordered proteins.

Identification of small molecule allosteric modulators that act as enhancers/disrupters of rhodopsin oligomerization

The elongated cilia of the outer segment of rod and cone photoreceptor cells can contain concentrations of visual pigments of up to 5 mM. The rod visual pigments, G protein–coupled receptors called rhodopsins, have a propensity to self-aggregate, a property conserved among many G protein–coupled receptors. However, the effect of rhodopsin oligomerization on G protein signaling in native cells is less clear. Here, we address this gap in knowledge by studying rod phototransduction. As the rod outer segment is known to adjust its size proportionally to overexpression or reduction of rhodopsin expression, genetic perturbation of rhodopsin cannot be used to resolve this question. Therefore, we turned to high-throughput screening of a diverse library of 50,000 small molecules and used a novel assay for the detection of rhodopsin dimerization. This screen identified nine small molecules that either disrupted or enhanced rhodopsin dimer contacts in vitro. In a subsequent cell-free binding study, we found that all nine compounds decreased intrinsic fluorescence without affecting the overall UV-visible spectrum of rhodopsin, supporting their actions as allosteric modulators. Furthermore, ex vivo electrophysiological recordings revealed that a disruptive, hit compound #7 significantly slowed down the light response kinetics of intact rods, whereas compound #1, an enhancing hit candidate, did not substantially affect the photoresponse kinetics but did cause a significant reduction in light sensitivity. This study provides a monitoring tool for future investigation of the rhodopsin signaling cascade and reports the discovery of new allosteric modulators of rhodopsin dimerization that can also alter rod photoreceptor physiology.

A multiscale approach for bridging the gap between potency, efficacy, and safety of small molecules directed at membrane proteins

Membrane proteins constitute a substantial fraction of the human proteome, thus representing a vast source of therapeutic drug targets. Indeed, newly devised technologies now allow targeting "undruggable" regions of membrane proteins to modulate protein function in the cell. Despite the advances in technology, the rapid translation of basic science discoveries into potential drug candidates targeting transmembrane protein domains remains challenging. We address this issue by harmonizing single molecule-based and ensemble-based atomistic simulations of ligand-membrane interactions with patient-derived induced pluripotent stem cell (iPSC)-based experiments to gain insights into drug delivery, cellular efficacy, and safety of molecules directed at membrane proteins. In this study, we interrogated the pharmacological activation of the cardiac Ca2+ pump (Sarcoplasmic reticulum Ca2+-ATPase, SERCA2a) in human iPSC-derived cardiac cells as a proof-of-concept model. The combined computational-experimental approach serves as a platform to explain the differences in the cell-based activity of candidates with similar functional profiles, thus streamlining the identification of drug-like candidates that directly target SERCA2a activation in human cardiac cells. Systematic cell-based studies further showed that a direct SERCA2a activator does not induce cardiotoxic pro-arrhythmogenic events in human cardiac cells, demonstrating that pharmacological stimulation of SERCA2a activity is a safe therapeutic approach targeting the heart. Overall, this novel multiscale platform encompasses organ-specific drug potency, efficacy, and safety, and opens new avenues to accelerate the bench-to-patient research aimed at designing effective therapies directed at membrane protein domains.

Characterization of Five Transmembrane Proteins: With Focus on the Tweety, Sideroflexin, and YIP1 Domain Families

Transmembrane proteins are involved in many essential cell processes such as signal transduction, transport, and protein trafficking, and hence many are implicated in different disease pathways. Further, as the structure and function of proteins are correlated, investigating a group of proteins with the same tertiary structure, i.e., the same number of transmembrane regions, may give understanding about their functional roles and potential as therapeutic targets. This analysis investigates the previously unstudied group of proteins with five transmembrane-spanning regions (5TM). More than half of the 58 proteins identified with the 5TM architecture belong to 12 families with two or more members. Interestingly, more than half the proteins in the dataset function in localization activities through movement or tethering of cell components and more than one-third are involved in transport activities, particularly in the mitochondria. Surprisingly, no receptor activity was identified within this dataset in large contrast with other TM groups. The three major 5TM families, which comprise nearly 30% of the dataset, include the tweety family, the sideroflexin family and the Yip1 domain (YIPF) family. We also analyzed the evolutionary origin of these three families. The YIPF family appears to be the most ancient with presence in bacteria and archaea, while the tweety and sideroflexin families are first found in eukaryotes. We found no evidence of common decent for these three families. About 30% of the 5TM proteins have prominent expression in the brain, liver, or testis. Importantly, 60% of these proteins are identified as cancer prognostic markers, where they are associated with clinical outcomes of various tumor types. Nearly 10% of the 5TMs are still not fully characterized and further investigation of their functional activities and expression is warranted. This study provides the first comprehensive analysis of proteins with the 5TM architecture, providing details of their unique characteristics.

Bridging protein structure, dynamics, and function using hydrogen/deuterium-exchange mass spectrometry

Much of our understanding of protein structure and mechanistic function has been derived from static high-resolution structures. As structural biology has continued to evolve it has become clear that high-resolution structures alone are unable to fully capture the mechanistic basis for protein structure and function in solution. Recently Hydrogen/Deuterium-exchange Mass Spectrometry (HDX-MS) has developed into a powerful and versatile tool for structural biologists that provides novel insights into protein structure and function. HDX-MS enables direct monitoring of a protein's structural fluctuations and conformational changes under native conditions in solution even as it is carrying out its functions. In this review, we focus on the use of HDX-MS to monitor these dynamic changes in proteins. We examine how HDX-MS has been applied to study protein structure and function in systems ranging from large, complex assemblies to intrinsically disordered proteins, and we discuss its use in probing conformational changes during protein folding and catalytic function. STATEMENT FOR A BROAD AUDIENCE: The biophysical and structural characterization of proteins provides novel insight into their functionalities. Protein motions, ranging from small scale local fluctuations to larger concerted structural rearrangements, often determine protein function. Hydrogen/Deuterium-exchange Mass Spectrometry (HDX-MS) has proven a powerful biophysical tool capable of probing changes in protein structure and dynamic protein motions that are often invisible to most other techniques.

Structure retention of silica gel-encapsulated bacteriorhodopsin in purple membrane and in lipid nanodiscs

The integral membrane protein, bacteriorhodopsin (BR) was encapsulated in sol-gel derived porous silica gel monoliths in native purple membrane (BR-PM) and synthetic lipid nanodisc (BR nanodisc) environments. BR nanodiscs were synthesized by solubilizing purple membrane in discoidal phospholipid bilayer stabilized by amphipathic Styrene-Maleic Acid (SMA) copolymer. UV-vis absorbance spectroscopy and dynamic-light scattering indicated the formation of BR monomers solubilized in lipid nanodiscs 10.2 ± 0.7 nm in average diameter. Fluorescence and absorbance spectroscopic techniques were utilized to probe conformational, environmental, and rotational changes associated with the tryptophan residues and the covalently-bound retinal moiety of BR upon entrapment in the silica matrix. We show that the immobilized BR in both membrane environments retained its bound retinal cofactor and the ability of the cofactor to undergo conformational changes upon light illumination necessary for BR's activity as a proton transporter. For purple membrane fragments, the results indicated that the local pH in the pores around BR after encapsulation was important for its stability at temperatures higher than 50 °C. Under the same buffering conditions, retinal was released from silica-encapsulated BR-PM and BR nanodiscs beginning at 80 °C (without a conformational change) and 50 °C (with a conformational change), respectively, reflecting differences in protein-protein (trimeric vs. monomeric) and protein-lipid interactions.

Shaping membranes with disordered proteins

Membrane proteins control and shape membrane trafficking processes. The role of protein structure in shaping cellular membranes is well established. However, a significant fraction of membrane proteins is disordered or contains long disordered regions. It becomes more and more clear that these disordered regions contribute to the function of membrane proteins. While the fold of a structured protein is essential for its function, being disordered seems to be a crucial feature of membrane bound intrinsically disordered proteins and protein regions. Here we outline the motifs that encode function in disordered proteins and discuss how these functional motifs enable disordered proteins to modulate membrane properties. These changes in membrane properties facilitate and regulate membrane trafficking processes which are highly abundant in eukaryotes.

Inter-residue interactions in alpha-helical transmembrane proteins

MOTIVATION

The number of available membrane protein structures has markedly increased in the last years and, in parallel, the reliability of the methods to detect transmembrane (TM) segments. In the present report, we characterized inter-residue interactions in α-helical membrane proteins using a dataset of 3462 TM helices from 430 proteins. This is by far the largest analysis published to date.

RESULTS

Our analysis of residue-residue interactions in TM segments of membrane proteins shows that almost all interactions involve aliphatic residues and Phe. There is lack of polar-polar, polar-charged and charged-charged interactions except for those between Thr or Ser sidechains and the backbone carbonyl of aliphatic and Phe residues. The results are discussed in the context of the preferences of amino acids to be in the protein core or exposed to the lipid bilayer and to occupy specific positions along the TM segment. Comparison to datasets of β-barrel membrane proteins and of α-helical globular proteins unveils the specific patterns of interactions and residue composition characteristic of α-helical membrane proteins that are the clue to understanding their structure.

AVAILABILITY AND IMPLEMENTATION

Results data and datasets used are available at http://lmc.uab.cat/TMalphaDB/interactions.php.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Site-Specific Peptide Probes Detect Buried Water in a Lipid Membrane

Transmembrane peptides contain polar residues in the interior of the membrane, which may alter the electrostatic environment and favor hydration in the otherwise nonpolar environment of the membrane core. Here, we demonstrate a general, nonperturbative strategy to probe hydration of the peptide backbone at specific depths within the bilayer using a combination of site-specific isotope labels, ultrafast two-dimensional infrared spectroscopy, and spectral modeling based on molecular dynamics simulations. Our results show that the amphiphilic pH-low insertion peptide supports a highly heterogeneous environment, with significant backbone hydration of nonpolar residues neighboring charged residues. For example, a leucine residue located as far as 1 nm into the hydrophobic bulk reports hydrogen-bonded populations as high as ∼20%. These findings indicate that the polar nature of these residues may facilitate the transport of water molecules into the hydrophobic core of the membrane.

ProLIF - quantitative integrin protein-protein interactions and synergistic membrane effects on proteoliposomes

Integrin transmembrane receptors control a wide range of biological interactions by triggering the assembly of large multiprotein complexes at their cytoplasmic interface. Diverse methods have been used to investigate interactions between integrins and intracellular proteins, and predominantly include peptide-based pulldowns and biochemical immuno-isolations from detergent-solubilised cell lysates. However, quantitative methods to probe integrin-protein interactions in a more biologically relevant context where the integrin is embedded within a lipid bilayer have been lacking. Here, we describe 'protein-liposome interactions by flow cytometry' (denoted ProLIF), a technique to reconstitute recombinant integrin transmembrane domains (TMDs) and cytoplasmic tail (CT) fragments in liposomes as individual subunits or as αβ heterodimers and, via flow cytometry, allow rapid and quantitative measurement of protein interactions with these membrane-embedded integrins. Importantly, the assay can analyse binding of fluorescent proteins directly from cell lysates without further purification steps. Moreover, the effect of membrane composition, such as PI(4,5)P₂ incorporation, on protein recruitment to the integrin CTs can be analysed. ProLIF requires no specific instrumentation and can be applied to measure a broad range of membrane-dependent protein-protein interactions with the potential for high-throughput/multiplex analyses.This article has associated First Person interviews with the first authors of the paper (see doi: 10.1242/jcs.223644 and doi: 10.1242/jcs.223719).

Protein-enriched outer membrane vesicles as a native platform for outer membrane protein studies

Most studies characterizing the folding, structure, and function of membrane proteins rely on solubilized or reconstituted samples. Whereas solubilized membrane proteins lack the functionally important lipid membrane, reconstitution embeds them into artificial lipid bilayers, which lack characteristic features of cellular membranes including lipid diversity, composition and asymmetry. Here, we utilize outer membrane vesicles (OMVs) released from Escherichia coli to study outer membrane proteins (Omps) in the native membrane environment. Enriched in the native membrane of the OMV we characterize the assembly, folding, and structure of OmpG, FhuA, Tsx, and BamA. Comparing Omps in OMVs to those reconstituted into artificial lipid membranes, we observe different unfolding pathways for some Omps. This observation highlights the importance of the native membrane environment to maintain the native structure and function relationship of Omps. Our fast and easy approach paves the way for functional and structural studies of Omps in the native membrane.

Johannes Thoma et al. overexpress outer membrane proteins (Omps) in Escherichia coli and collect the expelled outer membrane vesicles (OMVs) to study Omp assembly, folding and structure. They find that Omps in OMVs show different unfolding pathways compared to Omps reconstituted in artificial lipid membranes.

Genome-wide comparative analysis of putative Pth11-related G protein-coupled receptors in fungi belonging to Pezizomycotina

Background

G-protein coupled receptors (GPCRs) are the largest family of transmembrane receptors in fungi, where they play important roles in signal transduction. Among them, the Pth11-related GPCRs form a large and divergent protein family, and are only found in fungi in Pezizomycotina. However, the evolutionary process and potential functions of Pth11-related GPCRs remain largely unknown.

Results

Twenty genomes of fungi in Pezizomycotina covering different nutritional strategies were mined for putative Pth11-related GPCRs. Phytopathogens encode much more putative Pth11-related GPCRs than symbionts, saprophytes, or entomopathogens. Based on the phylogenetic tree, these GPCRs can be divided into nine clades, with each clade containing fungi in different taxonomic orders. Instead of fungi from the same order, those fungi with similar nutritional strategies were inclined to share orthologs of putative Pth11-related GPCRs. Most of the CFEM domain-containing Pth11-related GPCRs, which were only included in two clades, were detected in phytopathogens. Furthermore, many putative Pth11-related GPCR genes of phytopathogens were upregulated during invasive plant infection, but downregulated under biotic stress. The expressions of putative Pth11-related GPCR genes of saprophytes and entomopathogens could be affected by nutrient conditions, especially the carbon source. The gene expressions revealed that Pth11-related GPCRs could respond to biotic/abiotic stress and invasive plant infection with different expression patterns.

Conclusion

Our results indicated that the Pth11-related GPCRs existed before the diversification of Pezizomycotina and have been gained and/or lost several times during the evolutionary process. Tandem duplications and trophic variations have been important factors in this evolution.

Electronic supplementary material

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

Functional Contributions of Positive Charges in the Pore-Lining Helix 3 of the Bordetella pertussis CyaA-Hemolysin to Hemolytic Activity and Ion-Channel Opening

The Bordetella pertussis CyaA-hemolysin (CyaA-Hly) domain was previously demonstrated to be an important determinant for hemolysis against target erythrocytes and ion-channel formation in planar lipid bilayers (PLBs). Here, net-charge variations in the pore-lining helix of thirteen related RTX cytolysins including CyaA-Hly were revealed by amino acid sequence alignments, reflecting their different degrees of hemolytic activity. To analyze possible functional effects of net-charge alterations on hemolytic activity and channel formation of CyaA-Hly, specific mutations were made at Gln⁵⁷⁴ or Glu⁵⁸¹ in its pore-lining α3 of which both residues are highly conserved Lys in the three highly active RTX cytolysins (i.e., Escherichia coli α-hemolysin, Actinobacillus pleuropneumoniae toxin, and Aggregatibacter actinomycetemcomitans leukotoxin). All six constructed CyaA-Hly mutants that were over-expressed in E. coli as 126 kDa His-tagged soluble proteins were successfully purified via immobilized Ni²⁺-affinity chromatography. Both positive-charge substitutions (Q574K, Q574R, E581K, E581R) and negative-charge elimination (E581Q) appeared to increase the kinetics of toxin-induced hemolysis while the substitution with a negatively-charged side-chain (Q574E) completely abolished its hemolytic activity. When incorporated into PLBs under symmetrical conditions (1.0 M KCl, pH 7.4), all five mutant toxins with the increased hemolytic activity produced clearly-resolved single channels with higher open probability and longer lifetime than the wild-type toxin, albeit with a half decrease in their maximum conductance. Molecular dynamics simulations for 50 ns of a trimeric CyaA-Hly pore model comprising three α2-loop-α3 transmembrane hairpins revealed a significant role of the positive charge at both target positions in the structural stability and enlarged diameter of the simulated pore. Altogether, our present data have disclosed functional contributions of positively-charged side-chains substituted at positions Gln⁵⁷⁴ and Glu⁵⁸¹ in the pore-lining α3 to the enhanced hemolytic activity and ion-channel opening of CyaA-Hly that actually mimics the highly-active RTX (repeat-in-toxin) cytolysins.

Novel signaling of dynorphin at κ-opioid receptor/bradykinin B2 receptor heterodimers

The κ-opioid receptor (KOR) and bradykinin B2 receptor (B2R) are involved in a variety of important physiological processes and share many similar characteristics in terms of their distribution and functions in the nervous system. We first demonstrated the endogenous expression of KOR and B2R in human SH-SY5Y cells and their co-localization on the membrane of human embryonic kidney 293 (HEK293) cells. Bioluminescence and fluorescence resonance energy transfer and the proximity ligation assay were exploited to demonstrate the formation of functional KOR and B2R heteromers in transfected cells. KOR/B2R heteromers triggered dynorphin A (1-13)-induced Gαs/protein kinase A signaling pathway activity, including upregulation of intracellular cAMP levels and cAMP-response element luciferase reporter activity, resulting in increased cAMP-response element-binding protein (CREB) phosphorylation, which could be dampened by the protein kinase A (PKA) inhibitor H89. This indicated that the co-existence of KOR and B2R is critical for CREB phosphorylation. In addition, dynorphin A (1-13) induced a significantly higher rate of proliferation in HEK293-KOR/B2R and human SH-SY5Y cells than in the control group. These results indicate that KOR can form a heterodimer with B2R and this leads to increased protein kinase A activity by the CREB signaling pathway, leading to a significant increase in cell proliferation. The nature of this signaling pathway has significant implications for the role of dynorphin in the regulation of neuroprotective effects.

Topology based identification and comprehensive classification of four-transmembrane helix containing proteins (4TMs) in the human genome

BACKGROUND

Membrane proteins are key components in a large spectrum of diverse functions and thus account for the major proportion of the drug-targeted portion of the genome. From a structural perspective, the α-helical transmembrane proteins can be categorized into major groups based on the number of transmembrane helices and these groups are often associated with specific functions. When compared to the well-characterized seven-transmembrane containing proteins (7TM), other TM groups are less explored and in particular the 4TM group. In this study, we identify the complete 4TM complement from the latest release of the human genome and assess the 4TM structure group as a whole. We functionally characterize this dataset and evaluate the resulting groups and ubiquitous functions, and furthermore describe disease and drug target involvement.

RESULTS

We classified 373 proteins, which represents ~7 % of the human membrane proteome, and includes 69 more proteins than our previous estimate. We have characterized the 4TM dataset based on functional, structural, and/or evolutionary similarities. Proteins that are involved in transport activity constitute 37 % of the dataset, 23 % are receptor-related, and 13 % have enzymatic functions. Intriguingly, proteins involved in transport are more than double the 15 % of transporters in the entire human membrane proteome, which might suggest that the 4TM topological architecture is more favored for transporting molecules over other functions. Moreover, we found an interesting exception to the ubiquitous intracellular N- and C-termini localization that is found throughout the entire membrane proteome and 4TM dataset in the neurotransmitter gated ion channel families. Overall, we estimate that 58 % of the dataset has a known association to disease conditions with 19 % of the genes possibly involved in different types of cancer.

CONCLUSIONS

We provide here the most robust and updated classification of the 4TM complement of the human genome as a platform to further understand the characteristics of 4TM functions and to explore pharmacological opportunities.

Recognition and binding of voltage-dependent anion channel-1 with ATP and NADH by spectroscopic analysis and molecular docking

Recognition and binding of voltage-dependent anion channel-1 with ATP and NADH by spectroscopic analysis and molecular docking.

Disruption of Rhodopsin Dimerization with Synthetic Peptides Targeting an Interaction Interface

Although homo- and heterodimerizations of G protein-coupled receptors (GPCRs) are well documented, GPCR monomers are able to assemble in different ways, thus causing variations in the interactive interface between receptor monomers among different GPCRs. Moreover, the functional consequences of this phenomenon, which remain to be clarified, could be specific for different GPCRs. Synthetic peptides derived from transmembrane (TM) domains can interact with a full-length GPCR, blocking dimer formation and affecting its function. Here we used peptides corresponding to TM helices of bovine rhodopsin (Rho) to investigate the Rho dimer interface and functional consequences of its disruption. Incubation of Rho with TM1, TM2, TM4, and TM5 peptides in rod outer segment (ROS) membranes shifted the resulting detergent-solubilized protein migration through a gel filtration column toward smaller molecular masses with a reduced propensity for dimer formation in a cross-linking reaction. Binding of these TM peptides to Rho was characterized by both mass spectrometry and a label-free assay from which dissociation constants were calculated. A BRET (bioluminescence resonance energy transfer) assay revealed that the physical interaction between Rho molecules expressed in membranes of living cells was blocked by the same four TM peptides identified in our in vitro experiments. Although disruption of the Rho dimer/oligomer had no effect on the rates of G protein activation, binding of Gt to the activated receptor stabilized the dimer. However, TM peptide-induced disruption of dimer/oligomer decreased receptor stability, suggesting that Rho supramolecular organization could be essential for ROS stabilization and receptor trafficking.

Nanoscale monitoring of drug actions on cell membrane using atomic force microscopy

Knowledge of the nanoscale changes that take place in individual cells in response to a drug is useful for understanding the drug action. However, due to the lack of adequate techniques, such knowledge was scarce until the advent of atomic force microscopy (AFM), which is a multifunctional tool for investigating cellular behavior with nanometer resolution under near-physiological conditions. In the past decade, researchers have applied AFM to monitor the morphological and mechanical dynamics of individual cells following drug stimulation, yielding considerable novel insight into how the drug molecules affect an individual cell at the nanoscale. In this article we summarize the representative applications of AFM in characterization of drug actions on cell membrane, including topographic imaging, elasticity measurements, molecular interaction quantification, native membrane protein imaging and manipulation, etc. The challenges that are hampering the further development of AFM for studies of cellular activities are aslo discussed.

Single-molecule force spectroscopy of membrane proteins from membranes freely spanning across nanoscopic pores

Single-molecule force spectroscopy (SMFS) provides detailed insight into the mechanical (un)folding pathways and structural stability of membrane proteins. So far, SMFS could only be applied to membrane proteins embedded in native or synthetic membranes adsorbed to solid supports. This adsorption causes experimental limitations and raises the question to what extent the support influences the results obtained by SMFS. Therefore, we introduce here SMFS from native purple membrane freely spanning across nanopores. We show that correct analysis of the SMFS data requires extending the worm-like chain model, which describes the mechanical stretching of a polypeptide, by the cubic extension model, which describes the bending of a purple membrane exposed to mechanical stress. This new experimental and theoretical approach allows to characterize the stepwise (un)folding of the membrane protein bacteriorhodopsin and to assign the stability of single and grouped secondary structures. The (un)folding and stability of bacteriorhodopsin shows no significant difference between freely spanning and directly supported purple membranes. Importantly, the novel experimental SMFS setup opens an avenue to characterize any protein from freely spanning cellular or synthetic membranes.

Conformational rearrangements in the transmembrane domain of CNGA1 channels revealed by single-molecule force spectroscopy

Cyclic nucleotide-gated (CNG) channels are activated by binding of cyclic nucleotides. Although structural studies have identified the channel pore and selectivity filter, conformation changes associated with gating remain poorly understood. Here we combine single-molecule force spectroscopy (SMFS) with mutagenesis, bioinformatics and electrophysiology to study conformational changes associated with gating. By expressing functional channels with SMFS fingerprints in Xenopus laevis oocytes, we were able to investigate gating of CNGA1 in a physiological-like membrane. Force spectra determined that the S4 transmembrane domain is mechanically coupled to S5 in the open state, but S3 in the closed state. We also show there are multiple pathways for the unfolding of the transmembrane domains, probably caused by a different degree of α-helix folding. This approach demonstrates that CNG transmembrane domains have dynamic structure and establishes SMFS as a tool for probing conformational change in ion channels.

Cyclic nucleotide gated channels are activated after binding cyclic nucleotides. Here, using single molecule force spectroscopy, the authors reveal that cyclic nucleotide binding causes conformational changes and tighter coupling of the S4 helix to the pore forming domain.

Crystallization of proteins from crude bovine rod outer segments

Obtaining protein crystals suitable for X-ray diffraction studies comprises the greatest challenge in the determination of protein crystal structures, especially for membrane proteins and protein complexes. Although high purity has been broadly accepted as one of the most significant requirements for protein crystallization, a recent study of the Escherichia coli proteome showed that many proteins have an inherent propensity to crystallize and do not require a highly homogeneous sample (Totir et al., 2012). As exemplified by RPE65 (Kiser, Golczak, Lodowski, Chance, & Palczewski, 2009), there also are cases of mammalian proteins crystallized from less purified samples. To test whether this phenomenon can be applied more broadly to the study of proteins from higher organisms, we investigated the protein crystallization profile of bovine rod outer segment (ROS) crude extracts. Interestingly, multiple protein crystals readily formed from such extracts, some of them diffracting to high resolution that allowed structural determination. A total of seven proteins were crystallized, one of which was a membrane protein. Successful crystallization of proteins from heterogeneous ROS extracts demonstrates that many mammalian proteins also have an intrinsic propensity to crystallize from complex biological mixtures. By providing an alternative approach to heterologous expression to achieve crystallization, this strategy could be useful for proteins and complexes that are difficult to purify or obtain by recombinant techniques.

Folding and stability of the aquaglyceroporin GlpF: Implications for human aqua(glycero)porin diseases

Aquaporins are highly selective polytopic transmembrane channel proteins that facilitate the permeation of water across cellular membranes in a large diversity of organisms. Defects in aquaporin function are associated with common diseases, such as nephrogenic diabetes insipidus, congenital cataract and certain types of cancer. In general, aquaporins have a highly conserved structure; from prokaryotes to humans. The conserved structure, together with structural dynamics and the structural framework for substrate selectivity is discussed. The folding pathway of aquaporins has been a topic of several studies in recent years. These studies revealed that a conserved protein structure can be reached by following different folding pathways. Based on the available data, we suggest a complex folding pathway for aquaporins, starting from the insertion of individual helices up to the formation of the tetrameric aquaporin structure. The consequences of some known mutations in human aquaporin-encoding genes, which most likely affect the folding and stability of human aquaporins, are discussed.

Atomic force microscopy on plasma membranes from Xenopus laevis oocytes containing human aquaporin 4

Atomic force microscopy (AFM) is a unique tool for imaging membrane proteins in near-native environment (embedded in a membrane and in buffer solution) at ~1 nm spatial resolution. It has been most successful on membrane proteins reconstituted in 2D crystals and on some specialized and densely packed native membranes. Here, we report on AFM imaging of purified plasma membranes from Xenopus laevis oocytes, a commonly used system for the heterologous expression of membrane proteins. Isoform M23 of human aquaporin 4 (AQP4-M23) was expressed in the X. laevis oocytes following their injection with AQP4-M23 cRNA. AQP4-M23 expression and incorporation in the plasma membrane were confirmed by the changes in oocyte volume in response to applied osmotic gradients. Oocyte plasma membranes were then purified by ultracentrifugation on a discontinuous sucrose gradient, and the presence of AQP4-M23 proteins in the purified membranes was established by Western blotting analysis. Compared with membranes without over-expressed AQP4-M23, the membranes from AQP4-M23 cRNA injected oocytes showed clusters of structures with lateral size of about 10 nm in the AFM topography images, with a tendency to a fourfold symmetry as may be expected for higher-order arrays of AQP4-M23. In addition, but only infrequently, AQP4-M23 tetramers could be resolved in 2D arrays on top of the plasma membrane, in good quantitative agreement with transmission electron microscopy analysis and the current model of AQP4. Our results show the potential and the difficulties of AFM studies on cloned membrane proteins in native eukaryotic membranes.

Structure-based drug design for G protein-coupled receptors

Our understanding of the structural biology of G protein-coupled receptors has undergone a transformation over the past 5 years. New protein-ligand complexes are described almost monthly in high profile journals. Appreciation of how small molecules and natural ligands bind to their receptors has the potential to impact enormously how medicinal chemists approach this major class of receptor targets. An outline of the key topics in this field and some recent examples of structure- and fragment-based drug design are described. A table is presented with example views of each G protein-coupled receptor for which there is a published X-ray structure, including interactions with small molecule antagonists, partial and full agonists. The possible implications of these new data for drug design are discussed.

G protein-coupled receptor 56 regulates matrix production and motility of lung fibroblasts

Idiopathic pulmonary fibrosis is a chronic, progressive, and fatal fibrotic lung disease with a poor prognosis, but no effective treatment is available. G protein-coupled receptor 56 (GPR56) plays a role in cell adhesion and tumor progression, but its function in fibrogenesis has not been explored. In this in vitro study, we found that GPR56 in IPF fibroblasts was lower than in normal fibroblasts. GPR56 regulated the production of fibronectin and type I collagen, and also changed the migratory and invasive capacity of lung fibroblasts. However, it was not sufficient to activate some classic markers of fibroblast and myofibroblast, such as α-smooth muscle actin and fibroblast specific protein 1. These findings demonstrate that reduced expression of GPR56 in lung fibroblasts may be an important link with pulmonary fibrosis, playing a role in regulating some important fibroblast functions.

From atomic structures to neuronal functions of g protein-coupled receptors

G protein-coupled receptors (GPCRs) are essential mediators of signal transduction, neurotransmission, ion channel regulation, and other cellular events. GPCRs are activated by diverse stimuli, including light, enzymatic processing of their N-termini, and binding of proteins, peptides, or small molecules such as neurotransmitters. GPCR dysfunction caused by receptor mutations and environmental challenges contributes to many neurological diseases. Moreover, modern genetic technology has helped identify a rich array of mono- and multigenic defects in humans and animal models that connect such receptor dysfunction with disease affecting neuronal function. The visual system is especially suited to investigate GPCR structure and function because advanced imaging techniques permit structural studies of photoreceptor neurons at both macro and molecular levels that, together with biochemical and physiological assessment in animal models, provide a more complete understanding of GPCR signaling.

Calcium-dependent ligand binding and G-protein signaling of family B GPCR parathyroid hormone 1 receptor purified in nanodiscs

GPCRs mediate intracellular signaling upon external stimuli, making them ideal drug targets. However, little is known about their activation mechanisms due to the difficulty in purification. Here, we introduce a method to purify GPCRs in nanodiscs, which incorporates GPCRs into lipid bilayers immediately after membrane solubilization, followed by single-step purification. Using this approach, we purified a family B GPCR, parathyroid hormone 1 receptor (PTH1R), which regulates calcium and phosphate homeostasis and is a drug target for osteoporosis. We demonstrated that the purified PTH1R in nanodiscs can bind to PTH(1-34) and activate G protein. We also observed that Ca(2+) is a weak agonist of PTH1R, and Ca(2+) in millimolar concentration can switch PTH(1-34) from an inverse agonist to an agonist. Hence, our results show that nanodiscs are a viable vehicle for GPCR purification, enabling studies of GPCRs under precise experimental conditions without interference from other cellular or membrane components.

Proteomic and bioinformatic analysis of membrane proteome in type 2 diabetic mouse liver

Type 2 diabetes mellitus (T2DM) is the most prevalent and serious metabolic disease affecting people worldwide. T2DM results from insulin resistance of the liver, muscle, and adipose tissue. In this study, we used proteomic and bioinformatic methodologies to identify novel hepatic membrane proteins that are related to the development of hepatic insulin resistance, steatosis, and T2DM. Using FT-ICR MS, we identified 95 significantly differentially expressed proteins in the membrane fraction of normal and T2DM db/db mouse liver. These proteins are primarily involved in energy metabolism pathways, molecular transport, and cellular signaling, and many of them have not previously been reported in diabetic studies. Bioinformatic analysis revealed that 16 proteins may be related to the regulation of insulin signaling in the liver. In addition, six proteins are associated with energy stress-induced, nine proteins with inflammatory stress-induced, and 14 proteins with endoplasmic reticulum stress-induced hepatic insulin resistance. Moreover, we identified 19 proteins that may regulate hepatic insulin resistance in a c-Jun amino-terminal kinase-dependent manner. In addition, three proteins, 14-3-3 protein beta (YWHAB), Slc2a4 (GLUT4), and Dlg4 (PSD-95), are discovered by comprehensive bioinformatic analysis, which have correlations with several proteins identified by proteomics approach. The newly identified proteins in T2DM should provide additional insight into the development and pathophysiology of hepatic steatosis and insulin resistance, and they may serve as useful diagnostic markers and/or therapeutic targets for these diseases.

High-resolution imaging of 2D outer membrane protein F crystals by atomic force microscopy

In this chapter the methodological bases are provided to achieve subnanometer resolution on two-dimensional (2D) membrane protein crystals by atomic force microscopy (AFM). This is outlined in detail with the example of AFM studies of the outer membrane protein F (OmpF) from the bacterium Escherichia coli (E. coli). We describe in detail the high-resolution imaging of 2D OmpF crystals in aqueous solution and under near-physiological conditions. The topographs of OmpF, and stylus effects and artifacts encountered when imaging by AFM are discussed.

3D imaging and quantitative analysis of small solubilized membrane proteins and their complexes by transmission electron microscopy

Inherently unstable, detergent-solubilized membrane protein complexes can often not be crystallized. For complexes that have a mass of >300 kDa, cryo-electron microscopy (EM) allows their three-dimensional (3D) structure to be assessed to a resolution that makes secondary structure elements visible in the best case. However, many interesting complexes exist whose mass is below 300 kDa and thus need alternative approaches. Two methods are reviewed: (i) Mass measurement in a scanning transmission electron microscope, which has provided important information on the stoichiometry of membrane protein complexes. This technique is applicable to particulate, filamentous and sheet-like structures. (ii) 3D-EM of negatively stained samples, which determines the molecular envelope of small membrane protein complexes. Staining and dehydration artifacts may corrupt the quality of the 3D map. Staining conditions thus need to be optimized. 3D maps of plant aquaporin SoPIP2;1 tetramers solubilized in different detergents illustrate that the flattening artifact can be partially prevented and that the detergent itself contributes significantly. Another example discussed is the complex of G protein-coupled receptor rhodopsin with its cognate G protein transducin.

G protein-coupled receptors--recent advances

The years 2000 and 2007 witnessed milestones in current understanding of G protein-coupled receptor (GPCR) structural biology. In 2000 the first GPCR, bovine rhodopsin, was crystallized and the structure was solved, while in 2007 the structure of β(2)-adrenergic receptor, the first GPCR with diffusible ligands, was determined owing to advances in microcrystallization and an insertion of the fast-folding lysozyme into the receptor. In parallel with those crystallographic studies, the biological and biochemical characterization of GPCRs has advanced considerably because those receptors are molecular targets for many of currently used drugs. Therefore, the mechanisms of activation and signal transduction to the cell interior deduced from known GPCRs structures are of the highest importance for drug discovery. These proteins are the most diversified membrane receptors encoded by hundreds of genes in our genome. They participate in processes responsible for vision, smell, taste and neuronal transmission in response to photons or binding of ions, hormones, peptides, chemokines and other factors. Although the GPCRs share a common seven-transmembrane α-helical bundle structure their binding sites can accommodate thousands of different ligands. The ligands, including agonists, antagonists or inverse agonists change the structure of the receptor. With bound agonists they can form a complex with a suitable G protein, be phosphorylated by kinases or bind arrestin. The discovered signaling cascades invoked by arrestin independently of G proteins makes the GPCR activating scheme more complex such that a ligand acting as an antagonist for G protein signaling can also act as an agonist in arrestin-dependent signaling. Additionally, the existence of multiple ligand-dependent partial activation states as well as dimerization of GPCRs result in a 'microprocessor-like' action of these receptors rather than an 'on-off' switch as was commonly believed only a decade ago.

Another Piece of the Membrane Puzzle: Extending Slipids Further

To be able to model complex biological membranes in a more realistic manner, the force field Slipids (Stockholm lipids) has been extended to include parameters for sphingomyelin (SM), phosphatidylglycerol (PG), phosphatidylserine (PS) lipids, and cholesterol. Since the parametrization scheme was faithful to the scheme used in previous editions of Slipids, all parameters are consistent and fully compatible. The results of careful validation of a number of key structural properties for one and two component lipid bilayers are in excellent agreement with experiments. Potentials of mean force for transferring water across binary mixtures of lipids and cholesterol were also computed in order to compare water permeability rates to experiments. In agreement with experimental and simulation studies, it was found that the permeability and partitioning of water is affected by cholesterol in lipid bilayers made of saturated lipids to the largest extent. With the extensions of Slipids presented here, it is now possible to study complex systems containing many different lipids and proteins in a fully atomistic resolution in the isothermic-isobaric (NPT) ensemble, which is the proper ensemble for membrane simulations.

The role of tryptophan side chains in membrane protein anchoring and hydrophobic mismatch

Tryptophan (Trp) is abundant in membrane proteins, preferentially residing near the lipid-water interface where it is thought to play a significant anchoring role. Using a total of 3 μs of molecular dynamics simulations for a library of hydrophobic WALP-like peptides, a long poly-Leu α-helix, and the methyl-indole analog, we explore the thermodynamics of the Trp movement in membranes that governs the stability and orientation of transmembrane protein segments. We examine the dominant hydrogen-bonding interactions between the Trp and lipid carbonyl and phosphate moieties, cation-π interactions to lipid choline moieties, and elucidate the contributions to the thermodynamics that serve to localize the Trp, by ~4 kcal/mol, near the membrane glycerol backbone region. We show a striking similarity between the free energy to move an isolated Trp side chain to that found from a wide range of WALP peptides, suggesting that the location of this side chain is nearly independent of the host transmembrane segment. Our calculations provide quantitative measures that explain Trp's role as a modulator of responses to hydrophobic mismatch, providing a deeper understanding of how lipid composition may control a range of membrane active peptides and proteins.

An Extension and Further Validation of an All-Atomistic Force Field for Biological Membranes

Biological membranes are versatile in composition and host intriguing molecular processes. In order to be able to study these systems, an accurate model Hamiltonian or force field (FF) is a necessity. Here, we report the results of our extension of earlier developed all-atomistic FF parameters for fully saturated phospholipids that complements an earlier parameter set for saturated phosphatidylcholine lipids (J. Phys. Chem. B, 2012, 116, 3164-3179). The FF, coined Slipids (Stockholm lipids), now also includes parameters for unsaturated phosphatidylcholine and phosphatidylethanolamine lipids, e.g., POPC, DOPC, SOPC, POPE, and DOPE. As the extended set of parameters is derived with the same philosophy as previously applied, the resulting FF has been developed in a fully consistent manner. The capabilities of Slipids are demonstrated by performing long simulations without applying any surface tension and using the correct isothermal-isobaric (NPT) ensemble for a range of temperatures and carefully comparing a number of properties with experimental findings. Results show that several structural properties are very well reproduced, such as scattering form factors, NMR order parameters, thicknesses, and area per lipid. Thermal dependencies of different thicknesses and area per lipid are reproduced as well. Lipid diffusion is systematically slightly underestimated, whereas the normalized lipid diffusion follows the experimental trends. This is believed to be due to the lack of collective movement in the relatively small bilayer patches used. Furthermore, the compatibility with amino acid FFs from the AMBER family is tested in explicit transmembrane complexes of the WALP23 peptide with DLPC and DOPC bilayers, and this shows that Slipids can be used to study more complex and biologically relevant systems.

Reconstitution of glucosylceramide flip-flop across endoplasmic reticulum: implications for mechanism of glycosphingolipid biosynthesis

Most glycosphingolipids are synthesized by the sequential addition of monosaccharides to glucosylceramide (GlcCer) in the lumen of the Golgi apparatus. Because GlcCer is synthesized on the cytoplasmic face of Golgi membranes, it must be flipped to the non-cytoplasmic face by a lipid flippase in order to nucleate glycosphingolipid synthesis. Halter et al. (Halter, D., Neumann, S., van Dijk, S. M., Wolthoorn, J., de Mazière, A. M., Vieira, O. V., Mattjus, P., Klumperman, J., van Meer, G., and Sprong, H. (2007) Pre- and post-Golgi translocation of glucosylceramide in glycosphingolipid synthesis. J. Cell Biol. 179, 101-115) proposed that this essential flipping step is accomplished via a complex trafficking itinerary; GlcCer is moved from the cytoplasmic face of the Golgi to the endoplasmic reticulum (ER) by FAPP2, a cytoplasmic lipid transfer protein, flipped across the ER membrane, then delivered to the lumen of the Golgi complex by vesicular transport. We now report biochemical reconstitution studies to analyze GlcCer flipping at the ER. Using proteoliposomes reconstituted from Triton X-100-solubilized rat liver ER membrane proteins, we demonstrate rapid (t(½) < 20 s), ATP-independent flip-flop of N-(6-((7-nitro-2-1,3-benzoxadiazol-4-yl)amino)hexanoyl)-D-glucosyl-β1-1'-sphingosine, a fluorescent GlcCer analog. Further studies involving protein modification, biochemical fractionation, and analyses of flip-flop in proteoliposomes reconstituted with ER membrane proteins from yeast indicate that GlcCer translocation is facilitated by well characterized ER phospholipid flippases that remain to be identified at the molecular level. By reason of their abundance and membrane bending activity, we considered that the ER reticulons and the related Yop1 protein could function as phospholipid-GlcCer flippases. Direct tests showed that these proteins have no flippase activity.

Convergence of Free Energy Profile of Coumarin in Lipid Bilayer

Atomistic molecular dynamics (MD) simulations of druglike molecules embedded in lipid bilayers are of considerable interest as models for drug penetration and positioning in biological membranes. Here we analyze partitioning of coumarin in dioleoylphosphatidylcholine (DOPC) bilayer, based on both multiple, unbiased 3 μs MD simulations (total length) and free energy profiles along the bilayer normal calculated by biased MD simulations (∼7 μs in total). The convergences in time of free energy profiles calculated by both umbrella sampling and z-constraint techniques are thoroughly analyzed. Two sets of starting structures are also considered, one from unbiased MD simulation and the other from "pulling" coumarin along the bilayer normal. The structures obtained by pulling simulation contain water defects on the lipid bilayer surface, while those acquired from unbiased simulation have no membrane defects. The free energy profiles converge more rapidly when starting frames from unbiased simulations are used. In addition, z-constraint simulation leads to more rapid convergence than umbrella sampling, due to quicker relaxation of membrane defects. Furthermore, we show that the choice of RESP, PRODRG, or Mulliken charges considerably affects the resulting free energy profile of our model drug along the bilayer normal. We recommend using z-constraint biased MD simulations based on starting geometries acquired from unbiased MD simulations for efficient calculation of convergent free energy profiles of druglike molecules along bilayer normals. The calculation of free energy profile should start with an unbiased simulation, though the polar molecules might need a slow pulling afterward. Results obtained with the recommended simulation protocol agree well with available experimental data for two coumarin derivatives.

Emerging methodologies to investigate lipid-protein interactions

Cellular membranes are composed of hundreds of different lipids, ion channels, receptors and scaffolding complexes that act as signalling and trafficking platforms for processes fundamental to life. Cellular signalling and membrane trafficking are often regulated by peripheral proteins, which reversibly interact with lipid molecules in highly regulated spatial and temporal fashions. In most cases, one or more modular lipid-binding domain(s) mediate recruitment of peripheral proteins to specific cellular membranes. These domains, of which more than 10 have been identified since 1989, harbour structurally selective lipid-binding sites. Traditional in vitro and in vivo studies have elucidated how these domains coordinate their cognate lipids and thus how the parent proteins associate with membranes. Cellular activities of peripheral proteins and subsequent physiological processes depend upon lipid binding affinities and selectivity. Thus, the development of novel sensitive and quantitative tools is essential in furthering our understanding of the function and regulation of these proteins. As this field expands into new areas such as computational biology, cellular lipid mapping, single molecule imaging, and lipidomics, there is an urgent need to integrate technologies to detail the molecular architecture and mechanisms of lipid signalling. This review surveys emerging cellular and in vitro approaches for studying protein-lipid interactions and provides perspective on how integration of methodologies directs the future development of the field.