Identification of structural transitions in bacterial fatty acid binding proteins that permit ligand entry and exit at membranes

Fatty acid (FA) transfer proteins extract FA from membranes and sequester them to facilitate their movement through the cytosol. Detailed structural information is available for these soluble protein–FA complexes, but the structure of the protein conformation responsible for FA exchange at the membrane is unknown. Staphylococcus aureus FakB1 is a prototypical bacterial FA transfer protein that binds palmitate within a narrow, buried tunnel. Here, we define the conformational change from a “closed” FakB1 state to an “open” state that associates with the membrane and provides a path for entry and egress of the FA. Using NMR spectroscopy, we identified a conformationally flexible dynamic region in FakB1, and X-ray crystallography of FakB1 mutants captured the conformation of the open state. In addition, molecular dynamics simulations show that the new amphipathic α-helix formed in the open state inserts below the phosphate plane of the bilayer to create a diffusion channel for the hydrophobic FA tail to access the hydrocarbon core and place the carboxyl group at the phosphate layer. The membrane binding and catalytic properties of site-directed mutants were consistent with the proposed membrane docked structure predicted by our molecular dynamics simulations. Finally, the structure of the bilayer-associated conformation of FakB1 has local similarities with mammalian FA binding proteins and provides a conceptual framework for how these proteins interact with the membrane to create a diffusion channel from the FA location in the bilayer to the protein interior.     

Modularity of the hydrophobic core and evolution of functional diversity in fold A glycosyltransferases

Hydrophobic cores are fundamental structural properties of proteins typically associated with protein folding and stability; however, how the hydrophobic core shapes protein evolution and function is poorly understood. Here, we investigated the role of conserved hydrophobic cores in fold-A glycosyltransferases (GT-As), a large superfamily of enzymes that catalyze formation of glycosidic linkages between diverse donor and acceptor substrates through distinct catalytic mechanisms (inverting versus retaining). Using hidden Markov models and protein structural alignments, we identify similarities in the phosphate-binding cassette (PBC) of GT-As and unrelated nucleotide-binding proteins, such as UDP-sugar pyrophosphorylases. We demonstrate that GT-As have diverged from other nucleotide-binding proteins through structural elaboration of the PBC and its unique hydrophobic tethering to the F-helix, which harbors the catalytic base (xED-Asp). While the hydrophobic tethering is conserved across diverse GT-A fold enzymes, some families, such as B3GNT2, display variations in tethering interactions and core packing. We evaluated the structural and functional impact of these core variations through experimental mutational analysis and molecular dynamics simulations and find that some of the core mutations (T336I in B3GNT2) increase catalytic efficiency by modulating the conformational occupancy of the catalytic base between “D-in” and acceptor-accessible “D-out” conformation. Taken together, our studies support a model of evolution in which the GT-A core evolved progressively through elaboration upon an ancient PBC found in diverse nucleotide-binding proteins, and malleability of this core provided the structural framework for evolving new catalytic and substrate-binding functions in extant GT-A fold enzymes.     

Classification of α-Helical Membrane Proteins Using Predicted Helix Architectures

Despite significant methodological advances in protein structure determination high-resolution structures of membrane proteins are still rare, leaving sequence-based predictions as the only option for exploring the structural variability of membrane proteins at large scale. Here, a new structural classification approach for α-helical membrane proteins is introduced based on the similarity of predicted helix interaction patterns. Its application to proteins with known 3D structure showed that it is able to reliably detect structurally similar proteins even in the absence of any sequence similarity, reproducing the SCOP and CATH classifications with a sensitivity of 65% at a specificity of 90%. We applied the new approach to enhance our comprehensive structural classification of α-helical membrane proteins (CAMPS), which is primarily based on sequence and topology similarity, in order to find protein clusters that describe the same fold in the absence of sequence similarity. The total of 151 helix architectures were delineated for proteins with more than four transmembrane segments. Interestingly, we observed that proteins with 8 and more transmembrane helices correspond to fewer different architectures than proteins with up to 7 helices, suggesting that in large membrane proteins the evolutionary tendency to re-use already available folds is more pronounced.     

Powerful Sequence Similarity Search Methods and In-Depth Manual Analyses Can Identify Remote Homologs in Many Apparently “Orphan” Viral Proteins

The genome sequences of new viruses often contain many “orphan” or “taxon-specific” proteins apparently lacking homologs. However, because viral proteins evolve very fast, commonly used sequence similarity detection methods such as BLAST may overlook homologs. We analyzed a data set of proteins from RNA viruses characterized as “genus specific” by BLAST. More powerful methods developed recently, such as HHblits or HHpred (available through web-based, user-friendly interfaces), could detect distant homologs of a quarter of these proteins, suggesting that these methods should be used to annotate viral genomes. In-depth manual analyses of a subset of the remaining sequences, guided by contextual information such as taxonomy, gene order, or domain cooccurrence, identified distant homologs of another third. Thus, a combination of powerful automated methods and manual analyses can uncover distant homologs of many proteins thought to be orphans. We expect these methodological results to be also applicable to cellular organisms, since they generally evolve much more slowly than RNA viruses. As an application, we reanalyzed the genome of a bee pathogen, Chronic bee paralysis virus (CBPV). We could identify homologs of most of its proteins thought to be orphans; in each case, identifying homologs provided functional clues. We discovered that CBPV encodes a domain homologous to the Alphavirus methyltransferase-guanylyltransferase; a putative membrane protein, SP24, with homologs in unrelated insect viruses and insect-transmitted plant viruses having different morphologies (cileviruses, higreviruses, blunerviruses, negeviruses); and a putative virion glycoprotein, ORF2, also found in negeviruses. SP24 and ORF2 are probably major structural components of the virions.     

Structural Insights into the Catalytic Mechanism of Synechocystis Magnesium Protoporphyrin IX O-Methyltransferase (ChlM)

Magnesium protoporphyrin IX O-methyltransferase (ChlM) catalyzes transfer of the methyl group from S-adenosylmethionine to the carboxyl group of the C13 propionate side chain of magnesium protoporphyrin IX. This reaction is the second committed step in chlorophyll biosynthesis from protoporphyrin IX. Here we report the crystal structures of ChlM from the cyanobacterium Synechocystis sp. PCC 6803 in complex with S-adenosylmethionine and S-adenosylhomocysteine at resolutions of 1.6 and 1.7 Å, respectively. The structures illustrate the molecular basis for cofactor and substrate binding and suggest that conformational changes of the two “arm” regions may modulate binding and release of substrates/products to and from the active site. Tyr-28 and His-139 were identified to play essential roles for methyl transfer reaction but are not indispensable for cofactor/substrate binding. Based on these structural and functional findings, a catalytic model is proposed.     

Enzymatic Detoxication,Conformational Selection,and the Role of Molten Globule Active Sites

The role of conformational ensembles in enzymatic reactions remains unclear. Discussion concerning “induced fit” versus “conformational selection” has, however, ignored detoxication enzymes, which exhibit catalytic promiscuity. These enzymes dominate drug metabolism and determine drug-drug interactions. The detoxication enzyme glutathione transferase A1–1 (GSTA1–1), exploits a molten globule-like active site to achieve remarkable catalytic promiscuity wherein the substrate-free conformational ensemble is broad with barrierless transitions between states. A quantitative index of catalytic promiscuity is used to compare engineered variants of GSTA1–1 and the catalytic promiscuity correlates strongly with characteristics of the thermodynamic partition function, for the substrate-free enzymes. Access to chemically disparate transition states is encoded by the substrate-free conformational ensemble. Pre-steady state catalytic data confirm an extension of the conformational selection model, wherein different substrates select different starting conformations. The kinetic liability of the conformational breadth is minimized by a smooth landscape. We propose that “local” molten globule behavior optimizes detoxication enzymes.     

Rapid Identification of Sequences for Orphan Enzymes to Power Accurate Protein Annotation

The power of genome sequencing depends on the ability to understand what those genes and their proteins products actually do. The automated methods used to assign functions to putative proteins in newly sequenced organisms are limited by the size of our library of proteins with both known function and sequence. Unfortunately this library grows slowly, lagging well behind the rapid increase in novel protein sequences produced by modern genome sequencing methods. One potential source for rapidly expanding this functional library is the “back catalog” of enzymology – “orphan enzymes,” those enzymes that have been characterized and yet lack any associated sequence. There are hundreds of orphan enzymes in the Enzyme Commission (EC) database alone. In this study, we demonstrate how this orphan enzyme “back catalog” is a fertile source for rapidly advancing the state of protein annotation. Starting from three orphan enzyme samples, we applied mass-spectrometry based analysis and computational methods (including sequence similarity networks, sequence and structural alignments, and operon context analysis) to rapidly identify the specific sequence for each orphan while avoiding the most time- and labor-intensive aspects of typical sequence identifications. We then used these three new sequences to more accurately predict the catalytic function of 385 previously uncharacterized or misannotated proteins. We expect that this kind of rapid sequence identification could be efficiently applied on a larger scale to make enzymology’s “back catalog” another powerful tool to drive accurate genome annotation.     

The Crystal Structures of Yeast Get3 Suggest a Mechanism for Tail-Anchored Protein Membrane Insertion

Tail-anchored (TA) proteins represent a unique class of membrane proteins that contain a single C-terminal transmembrane helix. The post-translational insertion of the yeast TA proteins into the ER membrane requires the Golgi ER trafficking (GET) complex which contains Get1, Get2 and Get3. Get3 is an ATPase that recognizes and binds the C-terminal transmembrane domain (TMD) of the TA proteins. We have determined the crystal structures of Get3 from two yeast species, S. cerevisiae and D. hansenii, respectively. These high resolution crystal structures show that Get3 contains a nucleotide-binding domain and a “finger” domain for binding the TA protein TMD. A large hydrophobic groove on the finger domain of S. cerevisiae Get3 structure might represent the binding site for TMD of TA proteins. A hydrophobic helix from a symmetry-related Get3 molecule sits in the TMD-binding groove and mimics the TA binding scenario. Interestingly, the crystal structures of the Get3 dimers from S. cerevisiae and D. hansenii exhibit distinct conformations. The S. cerevisiae Get3 dimer structure does not contain nucleotides and maintains an “open” conformation, while the D. hansenii Get3 dimer structure binds ADP and stays in a “closed” conformation. We propose that the conformational changes to switch the Get3 between the open and closed conformations may facilitate the membrane insertions for TA proteins.     

Crystal Structures of the Human RNA Demethylase Alkbh5 Reveal Basis for Substrate Recognition

N⁶-Methylation of adenosine is the most ubiquitous and abundant modification of nucleoside in eukaryotic mRNA and long non-coding RNA. This modification plays an essential role in the regulation of mRNA translation and RNA metabolism. Recently, human AlkB homolog 5 (Alkbh5) and fat mass- and obesity-associated protein (FTO) were shown to erase this methyl modification on mRNA. Here, we report five high resolution crystal structures of the catalytic core of Alkbh5 in complex with different ligands. Compared with other AlkB proteins, Alkbh5 displays several unique structural features on top of the conserved double-stranded β-helix fold typical of this protein family. Among the unique features, a distinct “lid” region of Alkbh5 plays a vital role in substrate recognition and catalysis. An unexpected disulfide bond between Cys-230 and Cys-267 is crucial for the selective binding of Alkbh5 to single-stranded RNA/DNA by bringing a “flipping” motif toward the central β-helix fold. We generated a substrate binding model of Alkbh5 based on a demethylation activity assay of several structure-guided site-directed mutants. Crystallographic and biochemical studies using various analogs of α-ketoglutarate revealed that the active site cavity of Alkbh5 is much smaller than that of FTO and preferentially binds small molecule inhibitors. Taken together, our findings provide a structural basis for understanding the substrate recognition specificity of Alkbh5 and offer a foundation for selective drug design against AlkB members.     

Genetic Code-guided Protein Synthesis and Folding in Escherichia coli

Universal genetic codes are degenerated with 61 codons specifying 20 amino acids, thus creating synonymous codons for a single amino acid. Synonymous codons have been shown to affect protein properties in a given organism. To address this issue and explore how Escherichia coli selects its “codon-preferred” DNA template(s) for synthesis of proteins with required properties, we have designed synonymous codon libraries based on an antibody (scFv) sequence and carried out bacterial expression and screening for variants with altered properties. As a result, 342 codon variants have been identified, differing significantly in protein solubility and functionality while retaining the identical original amino acid sequence. The soluble expression level varied from completely insoluble aggregates to a soluble yield of ∼2.5 mg/liter, whereas the antigen-binding activity changed from no binding at all to a binding affinity of > 10⁻⁸ m. Not only does our work demonstrate the involvement of genetic codes in regulating protein synthesis and folding but it also provides a novel screening strategy for producing improved proteins without the need to substitute amino acids.     

Structome Analysis of Virulent Mycobacterium tuberculosis,Which Survives with Only 700 Ribosomes per 0.1 fl of Cytoplasm

We previously reported the exquisite preservation of the ultrastructures of virulent Mycobacterium tuberculosis cells processed through cryofixation and rapid freeze substitution. Here, we report the “structome” analysis (i.e., the quantitative three-dimensional structural analysis of a whole cell at the electron microscopic level) of virulent M. tuberculosis using serial ultrathin sections prepared after cryofixation and rapid freeze substitution and analyzed by transmission electron microscopy. Five M. tuberculosis cells, which were contained in the serial ultrathin cross sections encompassing from one end to the other, were cut into 24, 36, 69, 55, and 63 serial ultrathin sections, respectively. On average, the cells were 2.71 ± 1.05 μm in length, and the average diameter of the cell was 0.345 ± 0.029 μm. The outer membrane and plasma membrane surface areas were 3.04 ± 1.33 μm² and 2.67 ± 1.19 μm², respectively. The cell, outer membrane, periplasm, plasma membrane, and cytoplasm volumes were 0.293 ± 0.113 fl (= μm³), 0.006 ± 0.003 fl, 0.060 ± 0.021 fl, 0.019 ± 0.008 fl, and 0.210 ± 0.091 fl, respectively. The average total ribosome number was 1,672 ± 568, and the ribosome density was 716.5 ± 171.4/0.1 fl. This is the first report of a structome analysis of M. tuberculosis cells prepared as serial ultrathin sections following cryofixation and rapid freeze substitution and examined by transmission electron microscopy. These data are based on the direct measurement and enumeration of exquisitely preserved single-cell structures in transmission electron microscopy images rather than calculations or assumptions from indirect biochemical or molecular biological data. In addition, these data may explain the slow growth of M. tuberculosis and enhance understanding of the structural properties related to the expression of antigenicity, acid-fastness, and the mechanism of drug resistance, particularly in regard to the ratio of target to drug concentrations.     

Functional Insights into Human HMG-CoA Lyase from Structures of Acyl-CoA-containing Ternary Complexes

HMG-CoA lyase (HMGCL) is crucial to ketogenesis, and inherited human mutations are potentially lethal. Detailed understanding of the HMGCL reaction mechanism and the molecular basis for correlating human mutations with enzyme deficiency have been limited by the lack of structural information for enzyme liganded to an acyl-CoA substrate or inhibitor. Crystal structures of ternary complexes of WT HMGCL with the competitive inhibitor 3-hydroxyglutaryl-CoA and of the catalytically deficient HMGCL R41M mutant with substrate HMG-CoA have been determined to 2.4 and 2.2 Å, respectively. Comparison of these β/α-barrel structures with those of unliganded HMGCL and R41M reveals substantial differences for Mg²⁺ coordination and positioning of the flexible loop containing the conserved HMGCL “signature” sequence. In the R41M-Mg²⁺-substrate ternary complex, loop residue Cys²⁶⁶ (implicated in active-site function by mechanistic and mutagenesis observations) is more closely juxtaposed to the catalytic site than in the case of unliganded enzyme or the WT enzyme-Mg²⁺-3-hydroxyglutaryl-CoA inhibitor complex. In both ternary complexes, the S-stereoisomer of substrate or inhibitor is specifically bound, in accord with the observed Mg²⁺ liganding of both C3 hydroxyl and C5 carboxyl oxygens. In addition to His²³³ and His²³⁵ imidazoles, other Mg²⁺ ligands are the Asp⁴² carboxyl oxygen and an ordered water molecule. This water, positioned between Asp⁴² and the C3 hydroxyl of bound substrate/inhibitor, may function as a proton shuttle. The observed interaction of Arg⁴¹ with the acyl-CoA C1 carbonyl oxygen explains the effects of Arg⁴¹ mutation on reaction product enolization and explains why human Arg⁴¹ mutations cause drastic enzyme deficiency.     

Functional Analysis of a Rickettsial OmpA Homology Domain of Shigella flexneri IcsA

Shigella flexneri is a gram-negative bacterium that causes diarrhea and dysentery by invasion and spread through the colonic epithelium. Bacteria spread by assembling actin and other cytoskeletal proteins of the host into “actin tails” at the bacterial pole; actin tail assembly provides the force required to move bacteria through the cell cytoplasm and into adjacent cells. The 120-kDa S. flexneri outer membrane protein IcsA is essential for actin assembly. IcsA is anchored in the outer membrane by a carboxy-terminal domain (the β domain), such that the amino-terminal 706 amino acid residues (the α domain) are exposed on the exterior of the bacillus. The α domain is therefore likely to contain the domains that are important to interactions with host factors. We identify and characterize a domain of IcsA within the α domain that bears significant sequence similarity to two repeated domains of rickettsial OmpA, which has been implicated in rickettsial actin tail formation. Strains of S. flexneri and Escherichia coli that carry derivatives of IcsA containing deletions within this domain display loss of actin recruitment and increased accessibility to IcsA-specific antibody on the surface of intracytoplasmic bacteria. However, site-directed mutagenesis of charged residues within this domain results in actin assembly that is indistinguishable from that of the wild type, and in vitro competition of a polypeptide of this domain fused to glutathione S-transferase did not alter the motility of the wild-type construct. Taken together, our data suggest that the rickettsial homology domain of IcsA is required for the proper conformation of IcsA and that its disruption leads to loss of interactions of other IcsA domains within the amino terminus with host cytoskeletal proteins.     

A New Way to Degrade Heme: THE MYCOBACTERIUM TUBERCULOSIS ENZYME MhuD CATALYZES HEME DEGRADATION WITHOUT GENERATING CO*?

MhuD is an oxygen-dependent heme-degrading enzyme from Mycobacterium tuberculosis with high sequence similarity (∼45%) to Staphylococcus aureus IsdG and IsdI. Spectroscopic and mutagenesis studies indicate that the catalytically active 1:1 heme-MhuD complex has an active site structure similar to those of IsdG and IsdI, including the nonplanarity (ruffling) of the heme group bound to the enzyme. Distinct from the canonical heme degradation, we have found that the MhuD catalysis does not generate CO. Product analyses by electrospray ionization-MS and NMR show that MhuD cleaves heme at the α-meso position but retains the meso-carbon atom at the cleavage site, which is removed by canonical heme oxygenases. The novel tetrapyrrole product of MhuD, termed “mycobilin,” has an aldehyde group at the cleavage site and a carbonyl group at either the β-meso or the δ-meso position. Consequently, MhuD catalysis does not involve verdoheme, the key intermediate of ring cleavage by canonical heme oxygenase enzymes. Ruffled heme is apparently responsible for the heme degradation mechanism unique to MhuD. In addition, MhuD heme degradation without CO liberation is biologically significant as one of the signals of M. tuberculosis transition to dormancy is mediated by the production of host CO.     

Glutamate 90 at the Luminal Ion Gate of Sarcoplasmic Reticulum Ca2+-ATPase Is Critical for Ca2+ Binding on Both Sides of the Membrane

The roles of Ser⁷², Glu⁹⁰, and Lys²⁹⁷ at the luminal ends of transmembrane helices M1, M2, and M4 of sarcoplasmic reticulum Ca²⁺-ATPase were examined by transient and steady-state kinetic analysis of mutants. The dependence on the luminal Ca²⁺ concentration of phosphorylation by P_i (“Ca²⁺ gradient-dependent E2P formation”) showed a reduction of the apparent affinity for luminal Ca²⁺ in mutants with alanine or leucine replacement of Glu⁹⁰, whereas arginine replacement of Glu⁹⁰ or Ser⁷² allowed E2P formation from P_i even at luminal Ca²⁺ concentrations much too small to support phosphorylation in wild type. The latter mutants further displayed a blocked dephosphorylation of E2P and an increased rate of conversion of the ADP-sensitive E1P phosphoenzyme intermediate to ADP-insensitive E2P as well as insensitivity of the E2·BeF₃⁻ complex to luminal Ca²⁺. Altogether, these findings, supported by structural modeling, indicate that the E2P intermediate is stabilized in the mutants with arginine replacement of Glu⁹⁰ or Ser⁷², because the positive charge of the arginine side chain mimics Ca²⁺ occupying a luminally exposed low affinity Ca²⁺ site of E2P, thus identifying an essential locus (a “leaving site”) on the luminal Ca²⁺ exit pathway. Mutants with alanine or leucine replacement of Glu⁹⁰ further displayed a marked slowing of the Ca²⁺ binding transition as well as slowing of the dissociation of Ca²⁺ from Ca₂E1 back toward the cytoplasm, thus demonstrating that Glu⁹⁰ is also critical for the function of the cytoplasmically exposed Ca²⁺ sites on the opposite side of the membrane relative to where Glu⁹⁰ is located.     

Bacterial-like PPP protein phosphatases: Novel sequence alterations in pathogenic eukaryotes and peculiar features of bacterial sequence similarity

Reversible phosphorylation is a widespread modification affecting the great majority of eukaryotic cellular proteins, and whose effects influence nearly every cellular function. Protein phosphatases are increasingly recognized as exquisitely regulated contributors to these changes. The PPP (phosphoprotein phosphatase) family comprises enzymes, which catalyze dephosphorylation at serine and threonine residues. Nearly a decade ago, “bacterial-like” enzymes were recognized with similarity to proteins from various bacterial sources: SLPs (Shewanella-like phosphatases), RLPHs (Rhizobiales-like phosphatases), and ALPHs (ApaH-like phosphatases). A recent article from our laboratory appearing in Plant Physiology characterizes their extensive organismal distribution, abundance in plant species, predicted subcellular localization, motif organization, and sequence evolution. One salient observation is the distinct evolutionary trajectory followed by SLP genes and proteins in photosynthetic eukaryotes vs. animal and plant pathogens derived from photosynthetic ancestors. We present here a closer look at sequence data that emphasizes the distinctiveness of pathogen SLP proteins and that suggests that they might represent novel drug targets. A second observation in our original report was the high degree of similarity between the bacterial-like PPPs of eukaryotes and closely related proteins of the “eukaryotic-like” phyla Myxococcales and Planctomycetes. We here reflect on the possible implications of these observations and their importance for future research.     

Solution NMR Structure and Functional Analysis of the Integral Membrane Protein YgaP from Escherichia coli

The solution NMR structure of the α-helical integral membrane protein YgaP from Escherichia coli in mixed 1,2-diheptanoyl-sn-glycerol-3-phosphocholine/1-myristoyl-2-hydroxy-sn-glycero-3-phospho-(1′-rac-glycerol) micelles is presented. In these micelles, YgaP forms a homodimer with the two transmembrane helices being the dimer interface, whereas the N-terminal cytoplasmic domain includes a rhodanese-fold in accordance to its sequence homology to the rhodanese family of sulfurtransferases. The enzymatic sulfur transfer activity of full-length YgaP as well as of the N-terminal rhodanese domain only was investigated performing a series of titrations with sodium thiosulfate and potassium cyanide monitored by NMR and EPR. The data indicate the thiosulfate concentration-dependent addition of several sulfur atoms to the catalytic Cys-63, which process can be reversed by the addition of potassium cyanide. The catalytic reaction induces thereby conformational changes within the rhodanese domain, as well as on the transmembrane α-helices of YgaP. These results provide insights into a potential mechanism of YgaP during the catalytic thiosulfate activity in vivo.     

Phospho‐dependent signaling during the general stress response by the atypical response regulator and ClpXP adaptor RssB

In the model organism Escherichia coli and related species, the general stress response relies on tight regulation of the intracellular levels of the promoter specificity subunit RpoS. RpoS turnover is exclusively dependent on RssB, a two‐domain response regulator that functions as an adaptor that delivers RpoS to ClpXP for proteolysis. Here, we report crystal structures of the receiver domain of RssB both in its unphosphorylated form and bound to the phosphomimic BeF₃ ⁻. Surprisingly, we find only modest differences between these two structures, suggesting that truncating RssB may partially activate the receiver domain to a “meta‐active” state. Our structural and sequence analysis points to RssB proteins not conforming to either the Y–T coupling scheme for signaling seen in prototypical response regulators, such as CheY, or to the signaling model of the less understood FATGUY proteins.     

Breast Cancer-Derived Microparticles Display Tissue Selectivity in the Transfer of Resistance Proteins to Cells

Microparticles (MPs) play a vital role in cell communication by facilitating the horizontal transfer of cargo between cells. Recently, we described a novel “non-genetic” mechanism for the acquisition of multidrug resistance (MDR) in cancer cells by intercellular transfer of functional P-gp, via MPs. MDR is caused by the overexpression of the efflux transporters P-glycoprotein (P-gp) and Multidrug Resistance-Associated Protein 1 (MRP1). These transporters efflux anticancer drugs from resistant cancer cells and maintain sublethal intracellular drug concentrations. By conducting MP transfer experiments, we show that MPs derived from DX breast cancer cells selectively transfer P-gp to malignant MCF-7 breast cells only, in contrast to VLB₁₀₀ leukaemic cell-derived MPs that transfer P-gp and MRP1 to both malignant and non-malignant cells. The observed transfer selectivity is not the result of membrane restrictions for intercellular exchange, limitations in MP binding to recipient cells or the differential expression of the cytoskeletal protein, Ezrin. CD44 (isoform 10) was found to be selectively present on the breast cancer-derived MPs and not on leukaemic MPs and may contribute to the observed selective transfer of P-gp to malignant breast cells observed. Using the MCF-7 murine tumour xenograft model we demonstrated the stable transfer of P-gp by MPs in vivo, which was found to localize to the tumour core as early as 24 hours post MP exposure and to remain stable for at least 2 weeks. These findings demonstrate a remarkable capacity by MPs to disseminate a stable resistant trait in the absence of any selective pressure.