Global marine pollutants inhibit P-glycoprotein: Environmental levels, inhibitory effects, and cocrystal structure

See allHide authors and affiliations

Science Advances  15 Apr 2016:
Vol. 2, no. 4, e1600001
DOI: 10.1126/sciadv.1600001


The world’s oceans are a global reservoir of persistent organic pollutants to which humans and other animals are exposed. Although it is well known that these pollutants are potentially hazardous to human and environmental health, their impacts remain incompletely understood. We examined how persistent organic pollutants interact with the drug efflux transporter P-glycoprotein (P-gp), an evolutionarily conserved defense protein that is essential for protection against environmental toxicants. We identified specific congeners of organochlorine pesticides, polychlorinated biphenyls, and polybrominated diphenyl ethers that inhibit mouse and human P-gp, and determined their environmental levels in yellowfin tuna from the Gulf of Mexico. In addition, we solved the cocrystal structure of P-gp bound to one of these inhibitory pollutants, PBDE (polybrominated diphenyl ether)–100, providing the first view of pollutant binding to a drug transporter. The results demonstrate the potential for specific binding and inhibition of mammalian P-gp by ubiquitous congeners of persistent organic pollutants present in fish and other foods, and argue for further consideration of transporter inhibition in the assessment of the risk of exposure to these chemicals.


Persistent organic pollutants (POPs) are hazardous, man-made chemicals that endure in the environment and bioaccumulate in animals. Their environmental persistence ensues from properties such as halogenation and hydrophobicity that slow degradation and promote partitioning into organisms. At the same time, these properties also favor POP bioaccumulation by slowing their elimination. Indeed, although all animals have numerous metabolic enzymes, conjugation systems, and transporter proteins that normally act to eliminate xenobiotics, these systems appear ineffective at limiting POP bioaccumulation.

A critical step toward understanding the persistence and organismal impacts of POPs is defining their interactions with xenobiotic elimination systems. Drug transporters are plasma membrane proteins that both limit the entry of foreign chemicals into the body and speed their clearance, and are already well studied for their roles in drug disposition (1). Previous studies have suggested that environmental chemicals can also interact with drug transporters, such as P-glycoprotein (P-gp), but that they are poorly transported, and that these interactions ultimately lead to inhibition of transporter function (210). Of concern is that this inhibition reduces the efficacy of transport, thereby sensitizing animals to toxic chemicals that would otherwise be effluxed (11, 12).

Here, we took a multilevel approach to examine transporter-pollutant interactions, from levels in the environment down to the cocrystal structure of an environmental chemical bound to the transporter. We focused on P-gp, an adenosine triphosphate (ATP)–binding cassette (ABC) transporter (1315), which plays a major role in the disposition of xenobiotics (1, 16) and which is one of the best-studied drug transporters to date. P-gp has a large binding pocket that interacts with a wide variety of structurally divergent hydrophobic molecules (1619), and binding within this large pocket can have different impacts on the transporter, from stimulation to inhibition of function (20, 21). P-gp is conserved and is typically expressed at environmental barrier tissues, such as the small intestine or gills (2224).

To identify POPs that interact with P-gp, we used robust biochemical and cellular assays of mouse and human P-gp and identified specific congeners that inhibit this transporter. Using x-ray crystallography, we validated the binding of one of these chemicals, polybrominated diphenyl ether (PBDE)–100, deep within the ligand pocket of the transporter, providing the first snapshots of P-gp bound to a pollutant. To gain insight into the environmental relevance of P-gp inhibitors, we measured their levels in yellowfin tuna (Thunnus albacares) from the Gulf of Mexico (GOM) and used these data to examine the effects of a representative POP mixture on the transport function of the human P-gp.


Persistent pollutants inhibit P-gp

A gap in knowledge exists with regard to the identity and levels of environmental modulators of P-gp. To begin to address this question, we applied a systematic screen to test ubiquitous chemicals for effects on P-gp and to identify specific congeners most relevant to human and environmental health (Fig. 1). In this screen, we used two independent assays to test 37 pollutants for inhibitory effects, and published data on the incidence of the specific congeners in humans to identify relevant compounds (Fig. 1 and Table 1). We took advantage of mouse P-gp (fig. S1), an established model of mammalian P-gp that is amenable to high-resolution x-ray crystallography (25, 26) and has a high amino acid sequence homology with human P-gp (18).

Fig. 1 Identification of P-gp–inhibiting POPs.

Thirty-seven pollutants were tested for interactions with mouse P-gp using two independent assays (see Table 1). Sixteen compounds were identified as inhibitors in both assays. We focused on 10 congeners reported in humans on the basis of the literature (86, 87) and the Fourth National Report on Human Exposure to Environmental Chemicals of the U.S. Centers for Disease Control and Prevention (88). DDD, dichlorodiphenyldichloroethane; DDE, dichlorodiphenyldichloroethylene; DDT, dichlorodiphenyltrichloroethane; PCB, polychlorinated biphenyl; PBDE, polybrominated diphenyl ether.

Table 1 Summary of the interaction kinetics of environmental chemicals with mouse P-gp.

Compounds that were inhibitors in both assays and present in humans are in boldface. Compounds that were inhibitors in both assays but not known to be detected in humans are italicized. Y, yes; N, no; HCH, hexachlorocyclohexane; DDD, dichlorodiphenyldichloroethane; DDE, dichlorodiphenyldichloroethylene; DDT, dichlorodiphenyltrichloroethane; PFOA, perfluorooctanoic acid; PFOS, perfluorooctane sulfonate; DEHP, di(2-ethylhexyl)phthalate; NI, noninteraction; DL, dioxin-like PCB congeners.

View this table:

For chemosensitization assays, we expressed mouse P-gp in drug-sensitive yeast (Saccharomyces cerevisiae). In these strains, three endogenous transporters are deleted, rendering the cells sensitive to the cytotoxic effects of doxorubicin (DOX) (27). Expression of P-gp protects the cells against DOX (fig. S2), and inhibition of recombinant P-gp by exogenous inhibitors restores DOX sensitivity, causing a dose-dependent decrease in their growth (Fig. 2A). For ATPase assays, we used purified, detergent-solubilized protein and measured the liberation of inorganic phosphate in the presence or absence of chemicals (28, 29). To determine the potency of chemicals as inhibitors, P-gp was prestimulated with verapamil (30, 31) and then incubated with each pollutant individually.

Fig. 2 P-gp–inhibiting effects of POPs.

(A) Representative images showing chemosensitization of P-gp–expressing yeast by POPs. Inhibition of mouse P-gp heterologously expressed in drug (+DOX)–sensitive yeast. Inhibition is indicated by the reduction in yeast growth with increasing POP concentration. None of the POPs were toxic to yeast in the absence of DOX (−DOX). Yeast assays were replicated three times and representative micrographs are shown. (B) Upper panel: Inhibition of verapamil-stimulated P-gp adenosine triphosphatase (ATPase) activity by POPs. Graphs show P-gp ATPase inhibition kinetics with the 10 transporter-inhibiting POPs. ATPase assays were performed with purified, recombinant mouse P-gp protein. Lower panel: Lack of P-gp ATPase activation by POPs. ATPase activation of P-gp–inhibiting organochlorine pesticides (OCPs), PBDEs, and PCBs was determined with increasing concentrations of each compound and without verapamil prestimulation. The black curves show verapamil stimulation. Points were normalized to 100 μM verapamil stimulation (gray dashed line) and represent the average ATPase activity ± SD from three to six experiments. Where not visible, error bars are smaller than symbols; R2 values were all >0.99.

From our assays, 16 chemicals that exhibited inhibitory action in both assays were identified. Of these, 10 have been reported in humans (Fig. 1 and Table 1), including DDT, DDE, DDD, dieldrin, endrin, PCB-146, PCB-170, PCB-187, PBDE-47, and PBDE-100 (Fig. 1). All 10 compounds chemosensitized P-gp–overexpressing yeast toward DOX and had IC50 (median inhibitory concentration) values of 1.1 to 72.5 μM in the ATPase assays (Fig. 2 and Table 1). We also observed that congener-specific effects were similar in the two assays, suggesting common modes of interaction with the transporter. For instance, in the cellular assay, endrin had the most pronounced chemosensitizing effect, whereas its stereoisomer, dieldrin, had a markedly lower inhibitory effect (Fig. 2A). Similarly, in ATPase assays, endrin was nearly 20 times more potent than dieldrin (Fig. 2B and Table 1).

We also examined whether the 10 compounds might stimulate the basal (that is, unstimulated) activity of P-gp (32), which is sometimes seen when P-gp is incubated with its transported substrates, and we found that none of them significantly stimulated activity (Fig. 2B). This, along with their ready bioaccumulation, suggests that POPs are not well transported; however, such a possibility is not completely ruled out because the ATPase assay does not directly measure transport. In addition, none of these 10 compounds were toxic to yeast in the absence of DOX (Fig. 2A), indicating that the observed chemosensitizing effects were not related to nonspecific effects, such as membrane perturbation.

PBDE-100 binds to a conserved region deep within the drug-binding pocket of P-gp

To further investigate these interactions, we next sought to visualize the binding of a pollutant with P-gp. Until recently, little has been known about the mechanisms by which environmental chemicals might interact with P-gp, and an unresolved question has been whether their interactions involve specific binding to the ligand-binding domain or nonspecific interactions with the transporter. Although solving structures of large transmembrane proteins bound to inhibitors is a challenge because of the resolution range of ~3 to 4 Å, this limitation can be circumvented by the use of compounds that have an anomalous x-ray scatter signal from atoms such as selenium, mercury, or bromine. For this reason, we chose to cocrystallize and solve the structure of mouse P-gp bound to the penta-brominated compound PBDE-100. The resulting crystal diffracted to 3.5 Å resolution (table S1). The PBDE-100 binding site was located on one-half of the transporter’s pseudo-symmetric dimer structure (Fig. 3, A and B). Fifteen residues from transmembrane segments 5, 6, 7, 8, and 12 mediated hydrophobic interactions with the diphenyl core of PBDE-100 (Fig. 3, C and D), burying more than 90% of its solvent-accessible surface. Eleven of these residues (Y303, Y306, A307, F310, F331, F724, I727, V731, S752, F755, and S975) were unique to this structure and not shared with those previously identified in structures of P-gp crystallized with a series of selenium-labeled cyclopeptide ligands (26). Alternative side-chain conformations were observed for all interacting residues upon PBDE-100 binding (fig. S3).

Fig. 3 POP interactions at the substrate-binding site of mouse P-gp.

(A) Structure of mouse P-gp cocrystallized with PBDE-100. (B) Location of PBDE-100 at a distinct binding site in the internal cavity of P-gp, viewed from the intracellular side. TM, transmembrane. (C) 2mFoDFc electron density (where m is the figure of merit and D is the Sigma-A weighting factor) for PBDE-100 (blue; contour level of 1.2σ) and anomalous difference density peaks (purple; contour level of 3.5σ). (D) Stereo view of the binding pocket, with key residues important for the interaction with the diphenyl backbone of PBDE-100 shown as sticks. (E) Conserved binding site for PBDE-100. Top: Side chains found to interact with PBDE-100 are shown in blue (conserved in human and mouse) or green (not conserved). These residues are Y303, Y306, A307, F310, F331, Q721, F724, S725, I727, F728, V731, S752, F755, S975, and F979. Bottom: Amino acid sequence alignment of mouse and human P-gp highlighting the 15 interacting residues with PBDE-100 in TM5, TM6, TM7, TM8, and TM12.

To compare PBDE-100 binding sites of human and mouse P-gp, we aligned the interacting regions identified in our cocrystal structure (Fig. 3E). In addition, we examined potential conservation or divergence of this region in human, mouse, zebrafish, and sea urchin P-gp (fig. S4). These comparisons revealed a high degree of similarity in PBDE-100 binding residues, with 11 of the 15 residues being identical in vertebrates and with 13 being identical in humans and mice (Fig. 3E). Nine of these residues were conserved in sea urchins (fig. S4), which diverge from humans at the base of the deuterostome lineage (33). Collectively, these results suggest evolutionary conservation of the newly identified PBDE-100 binding sites.

Environmentally relevant inhibitor mixtures impair the transport function of human P-gp

The chemicals that we found to inhibit P-gp are ubiquitous, legacy pollutants that have been widely reported in humans and wildlife. Humans are likely to be exposed to these compounds through consumption of contaminated food, raising questions on the environmental levels and impacts of these inhibitor mixtures. Previous studies have demonstrated that marine environments are major global repositories of POPs (34) and that fish can have high levels of these compounds as compared to other foods (35, 36). Thus, to determine the environmental levels of the 10 inhibitors, we measured them in wild-caught yellowfin tuna. We focused on tuna because they are among the most widely consumed fish in the world (37). More than 1 million metric tons of yellowfin tuna are caught annually, accounting for 27% of the global tuna catch (38).

White dorsal muscle fillets of eight fish (Fig. 4A and table S2) were collected and analyzed for levels of OCPs, PCBs, PBDEs, and perfluorochemicals (PFCs) (fig. S5). Nine of the 10 inhibitory POPs were detected in these tuna, and the mean lipid-normalized concentration of the nine inhibitors was 1.01 μM, accounting for approximately 23% of total POPs (Fig. 4, B and C, and table S2). PCBs accounted for 79% of total POPs (fig. S5 and table S2) and 44% of inhibitors (Fig. 4C and table S2).

Fig. 4 Levels of P-gp inhibitors in yellowfin tuna (T. albacares).

(A) Sampling site for the eight yellowfin tunas (T. albacares) caught in the GOM. The inset shows a yellowfin tuna with the sampled dorsal muscle tissue marked in red. (B) Lipid-normalized concentrations of the total POPs and the 10 P-gp inhibitors. The red-filled circles represent the minimum and maximum values. The white diamonds represent the mean value. The horizontal lines represent the 50th percentile, and the boxes represent the 25th and 75th percentiles. (C) Range of concentrations of nine inhibitory POPs measured in yellowfin tuna muscle from the GOM.

On the basis of these data, a pollutant mixture was formulated, reflecting relative levels of inhibitors found in tuna (Fig. 5A). We used this mixture to examine effects on the transport function of human P-gp. In this assay, inside-out membrane vesicles of mammalian cells overexpressing human P-gp were generated, and effects of POPs were assayed by measurement of the accumulation of the P-gp substrate N-methylquinidine (NMQ) into the vesicles (fig. S6A). In the presence of the POP mixture, NMQ transport was inhibited with an IC50 of 28.7 μM and an IC10 of 7.1 μM (Fig. 5B). We also verified the inhibition of NMQ transport by the individual compounds (fig. S6B). The effects of the POP mixture on mouse P-gp ATPase activity were similar to those seen on human P-gp, with an IC50 of 25 μM and an IC10 of 5.7 μM (Fig. 5C). Collectively, these results indicate that POPs identified in our assays are likely to exhibit similar inhibitory effects on human P-gp.

Fig. 5 An environmentally relevant POP mixture inhibits the transport function of human and mouse P-gp.

(A) Relative ratio of the mean concentrations of P-gp inhibitors found in yellowfin tuna. (B) Inhibition of human P-gp by the POP mixture. Points represent the average percentage of NMQ uptake ± SD relative to the control from nine different experiments and with increasing concentration of the POP mixture. (C) Inhibition of verapamil-stimulated ATPase activity of mouse P-gp by the POP mixture. Shown is the respective dose-response curve as ATPase activity relative to 100 μM verapamil stimulation. The ATPase activity of the purified protein was measured in the presence of increasing concentrations of the POP mixture on the basis of the relative concentration of nine inhibitory POPs identified in this study. All data were fitted using a Hill function [y = v1 + (v2v1) * xn/(kn + xn)]. The R2 value was >0.99.


Here, we show the direct binding of a pollutant to the ligand-binding pocket of P-gp and the potential for POPs to inhibit this transporter. POPs are ubiquitous contaminants and P-gp is a key protein for xenobiotic elimination in all animals. P-gp, along with cytochrome P450, is typically expressed on apical plasma membranes at sites of xenobiotic uptake (39, 40), such as the intestine, and plays a key role in first-pass elimination of dietary toxins. Thus, the inhibition of P-gp by POPs could represent a mechanism by which this critical cellular defense becomes compromised.

A notable finding of this study was that the levels of transporter inhibitors were sometimes high, raising concerns about transporter inhibition after consumption of highly contaminated foods and about exposure to inhibitors in vulnerable populations of humans and wildlife. High exposure could result from consumption of individual fish with elevated levels of inhibitors. For instance, in this study, although the mean tuna inhibitor concentration was 1.1 μM, two of our eight fish had 2.8 and 2.9 μM POP inhibitors and 11.9 and 12.6 μM total POPs (table S2). This is notable considering that we measured levels in lean cuts of tuna and that more fatty cuts or other types of fish (41) could have even higher levels.

Vulnerable populations could include human neonates, who can have limited intestinal xenobiotic metabolism activity in the first several months of life (42, 43) and, in certain places, can be exposed to high levels of POPs through breast milk. Indeed, in areas where DDT was used for malaria control, high concentrations (>28 μM) were reported in milk fat (44, 45). In wildlife, P-gp inhibitors may pose similar dangers to vulnerable life history stages and geographic populations. For instance, in the GOM, there are high levels of oil hydrocarbons (46, 47) and P-gp could be important for protection against their toxic effects (4851). Inhibition of P-gp by POPs could be especially problematic for early life history stages because the developing fish heart is sensitive to polycyclic aromatic hydrocarbons (52, 53), and high levels of POPs could be maternally loaded through the yolk (54, 55).

Although our study sheds light on the potential effects of specific POPs on P-gp, the mechanisms of their inhibitory interactions remain incompletely understood. The compounds that we identified are hydrophobic halogenated compounds (table S3), with relatively high potential for passive movement from the environment into tissues (56). Although we did not directly determine the binding constants of these POPs to P-gp, the cocrystal structure demonstrates that a representative compound from our study can form intimate interactions with the ligand-binding domain of the transporter. The high passive flux of POPs into cells and their specific binding properties could be factors that reduce the effectiveness of P-gp in eliminating POPs (5759). In addition, the conformational changes of the protein that occur upon ligand binding (fig. S3) might also play some as yet unknown role in the observed inhibition by POPs (60, 61).

Finally, it will be important for future studies to consider how metabolism could alter the inhibitory potential and/or transport of pollutants. For instance, at sites of excretion, metabolism pathways, including those involved in conjugation and oxidation, might mitigate some of the inhibitory effects of pollutants and alter their potential for transport. Indeed, in mice, PBDE-47 is conjugated to carrier proteins that improve its efflux by P-gp into urine (62). In addition, other ABC transporters, such as ABCC1, ABCC2, and ABCG2, show a broad tissue distribution and can serve as an additional line of defense by expelling modified and/or unmodified pollutants (16, 24, 40, 63, 64).


The idea that environmental pollutants could perturb P-gp was first postulated nearly 20 years ago (65); however, the precise mechanisms and environmental relevance have largely remained a mystery. Here, we show the inhibition of the transporter by ubiquitous pollutants found in fish, and likely in other foods. Transporter interactions are already among the criteria evaluated in assessing drug availability (1, 66), and our results argue for consideration of these interactions in risk assessment of environmental chemicals.

An important step taken by this study was solving the cocrystal structure of P-gp bound to a pollutant. The structure revealed a high degree of conservation in the PBDE-100 binding residues and provides a useful tool to help explain and potentially predict both similarities and differences in compound selectivity across species (67). These interactions could also help experimentally define binding sites that are responsible for the adaptive evolution of P-gp into anthropogenic pollutants (68) and could prove applicable to the design of chemicals (6971) with better potential for elimination by transporters. Indeed, with other transporter structures almost within reach, transporter interactions could soon represent an exciting new avenue for green chemistry.



All OCPs, cyclosporine A (CsA), verapamil, DOX, Triton X-100, pesticides, and dimethyl sulfoxide (DMSO) were purchased from Sigma. All PCB and PBDE congeners were purchased from AccuStandard. With the exception of verapamil, which was dissolved in water, all stock solutions were prepared in DMSO and diluted to the final concentrations in reaction buffer (ATPase assay) or sterile water (yeast assay). The yeast nitrogen base without amino acids was purchased from Fluka, and the amino acid dropout supplement (−Leu) was obtained from Clontech. The final DMSO concentration in the ATPase and yeast cell assays did not exceed 2 and 0.5%, respectively.

Expression of mouse P-gp in DOX-sensitive yeast

Mouse P-gp (MmABCB1a) was cloned into Bam HI/Hind III restriction sites of the constitutive yeast expression vector p415GPD (Ampr LEU2 ARS/CEN) (72). The hypersensitized MATa S. cerevisiae strain lacking three ABC transporters (Δpdr5 Δsnq2 Δyor1) (27) was transformed with p415GPD plasmids containing MmABCB1a using the standard lithium acetate/polyethylene glycol method (73). The empty p415GPD plasmid was transformed as a control. Transformed cells were selected on synthetic complete dropout (SC-Leu) plates. A single colony from each transformation was picked and cultured overnight in SC-Leu media at 30°C for the spot assays. Overnight cultures were washed in sterile water and their concentrations were adjusted to an OD600 (optical density at 600 nm) of 0.1. To test the functional expression of mouse ABCB1 in yeast cells, fivefold serial dilutions of transformed cells were spotted on SC-Leu and SC-Leu + 10 μM DOX plates and incubated at 30°C for 36 and 48 hours, respectively. For the chemosensitization assay, 15 μl of the culture of p415GPD and p415GPD-MmABCB1a transformed cells was soaked in indicated concentrations of pollutants and spotted on SC-Leu and SC-Leu + 10 μM DOX plates. SC-Leu plates were incubated at 30°C for 40 hours and SC-Leu + DOX plates were incubated at 30°C for 48 to 72 hours.

Expression and purification of mouse P-gp protein

To generate purified mouse P-gp, we overexpressed both genes in Pichia pastoris and purified the protein using combined affinity tag and size exclusion chromatography (SEC). The expression and purification of mouse P-gp in P. pastoris were described previously (18, 25, 26). Briefly, P. pastoris transformed with yeast codon–optimized mouse P-gp (mouse ABCB1a, GenBank JF834158) harboring a C-terminal 6×His-tag was grown as 10-liter cultures in a Bioflow 415 bioreactor (New Brunswick Scientific), induced using slow methanol induction, and harvested. Cells were lysed at 40,000 psi by a single pass through a constant cell disrupter (TS-Series; Constant Systems Inc.), cell debris was separated by centrifugation at 12,500g, and membranes were isolated by centrifugation at 38,400g. Membranes were then solubilized and P-gp was purified using a Ni–nitrilotriacetic acid Superflow resin (Qiagen) via fast protein liquid chromatography (AEKTA, GE Life Sciences). The protein was concentrated (Centricon YM-100; Millipore), ultracentrifuged, and subjected to SEC (Superdex 200 16/60; GE Healthcare). The protein concentration was determined using the Micro BCA Protein Assay Kit (Pierce).

To evaluate the purity of the recombinant P-gp, 5 μg was separated by electrophoresis in a 7.5% SDS–polyacrylamide gel electrophoresis gel and subjected to immunoblotting using a 6×His epitope tag antibody (fig. S1A) (Pierce). The proteins were transferred to a polyvinylidene difluoride membrane by wet electroblotting (tank transfer) for 30 min at 100 V in transfer buffer [25 mM tris-base, 150 mM glycine, and 20% (v/v) methanol (pH 7.4)]. After the transfer, the membranes were blocked in 5% skim milk powder/1× TBST [0.1 M tris-base, 150 mM NaCl, and 0.05% Tween 20 (pH 7.4)] overnight at 23°C. After three washes with 1× TBST, the membranes were incubated for 2 hours with a 6×His epitope tag monoclonal antibody (mouse) in 5% skim milk powder/1× TBST at a ratio of 1:2000 (Pierce). Goat anti-mouse immunoglobulin G–horseradish peroxidase was used as secondary antibody in 5% skim milk powder/1× TBST at a ratio of 1:5000 for 1 hour at 23°C (Pierce). The proteins were visualized using the SuperSignal West Pico Chemiluminescent Substrate Kit (Pierce).

The yeast codon–optimized mouse P-gp used in this study has three N-glycosylation sites converted into glutamine (N83Q, N87Q, and N90Q) and a C-terminal His6-tag fusion. The calculated molecular mass of the resulting protein is ~142 kD, which was confirmed by matrix-assisted laser desorption/ionization mass spectrometry (MALDI MS) (fig. S1B). MALDI MS analysis was used with a time-of-flight mass detector on a Voyager Mass Spectrometer DE-STR (Applied Biosystems) and an effective ion path length of 2 m in the positive ion reflector mode. Sinapinic acid in 50% acetonitrile and 0.1% trifluoroacetic acid was used as matrix solution. Samples (1 mg/ml) were diluted 1:20 with the matrix solution and 1 μl was spotted onto the MALDI sample target plates and air-dried. Spectra were obtained in the mass range between 5000 and 200,000 daltons with 256 laser shots per spectrum. Accelerating voltage was 25,000 V with a grid voltage of 93% and a guide wire percentage of 0.3%. Spectra were recorded in delayed extraction with a delay time of 700 ns. Internal calibration was performed using bovine serum albumin standard (Sigma-Aldrich) with a calculated molecular mass of 66.5 kD. All data were analyzed using Voyager Data Explorer (Applied Biosystems) and plotted using Origin 7.0 (OriginLab).

ATPase activity of purified mouse P-gp

To test for molecular interactions of pollutants with the P-gp transporters, we optimized an ATPase assay based on the malachite green method (28, 74). This assay is well established (29, 75), requiring only 1 to 2 μg of total protein per reaction for measurements of the ATPase activity of efflux transporters. Figure S1C shows the respective dose-response curves of mouse P-gp with the model ATPase stimulator verapamil and the model inhibitor CsA. Activation of mouse P-gp ATPase using verapamil resulted in a half maximal stimulation concentration of 9.4 ± 0.7 μM, which is in agreement with that found previously (76, 77). Verapamil-stimulated (100 μM) ATPase activity was inhibited by CsA with an IC50 of 1.3 ± 0.1 μM. Briefly, we used 2 μg of purified, solubilized mouse P-gp and added the protein to the wells of a chilled 96-well plate containing 60 μl of ATP-free reaction buffer [10 mM MgSO4, 0.05% (w/v) DDM, 1 mM tris(2-carboxyethyl)phosphine (TCEP), and Escherichia coli polar extract lipids (0.1 mg/ml) in 50 mM tris-Cl buffer (pH 7.5)] with serial dilutions of verapamil (control activator). To test inhibition, we used serial dilutions of CsA (control inhibitor) or pollutant compounds plus 100 μM verapamil. Then, 60 μl of ATP solution [5 mM Na-ATP, 10 mM MgSO4, 0.05% (w/v) DDM, 1 mM TCEP, and E. coli polar extract lipids (0.1 mg/ml) in 50 mM tris-Cl buffer (pH 7.5)] was added, mixed, and incubated for 3 min on ice. After incubation, the reaction mixtures in the 96-well polymerase chain reaction plate were transferred to a thermocycler and the reaction was started with the following cycling parameters: 3 s at 4°C, 5 min at 37°C, 15 s at 80°C (heat inactivation), and hold at 4°C. ATPase reactions (30 μl) were transferred to a 96-well enzyme-linked immunosorbent assay plate and the liberated inorganic phosphate was measured by adding 150 μl of an activated color development solution [17 mg of malachite green in 3.75 ml of Milli-Q H2O and 0.525 g of ammonium molybdate tetrahydrate in 12.5 ml of 4N HCl, activated with 0.02% (v/v) Triton X-100] in each sample well. The absorbance of each sample was immediately measured at 600 nm in a microplate reader (Spectramax M2) to minimize acid-catalyzed ATP hydrolysis in the color development solution (pH ~1 to 2). Control samples containing buffer and DMSO (CsA, pollutants) or H2O (verapamil) without added P-gp protein were subtracted as background values. Inorganic phosphate standards (KH2PO4) from 0.125 to 2 nmol served as controls. The ATPase activity data are given as mean ± SEM from three to six different measurements, and representative experiments are shown. For the calculation of IC50 values, the data were fitted to a Hill function: y = v1 + (v2v1) * xn/(kn + xn), where v1 and v2 are the initial and final reaction velocities, respectively, n is the Hill coefficient or the cooperativity of the dependence on x, and k is the concentration of the inhibitor (IC50) or the stimulator (EC50) that corresponds to 50% of the maximum (Michaelis constant). All calculations were performed using Origin software.

Reductive methylation and crystallization of mouse P-gp

We added freshly made borane dimethylamine (50 mM) and formaldehyde (100 mM) to purified mouse P-gp (~1 to 2 mg). The mixture was incubated (2 hours at 4°C), and the reaction was quenched with glycine, washed via concentration in SEC buffer [20 mM tris-HCl (pH 8.0), 20 mM NaCl, 0.01% LMNG (lauryl maltose neopentyl glycol), 0.01% sodium cholate, 0.2 mM TCEP, and 0.5 mM EDTA], and treated with 0.2 to 0.5 mM PBDE-100 dissolved in DMSO for an overnight incubation in the dark (12 to 16 hours at 4°C). The sample was then diluted 1:10 in SEC buffer, concentrated to ~12 to 15 mg/ml, and used for crystal trials. P-gp crystals were grown in 24-well Cryschem plates (Hampton Research) at a protein concentration of ~12 to 15 mg/ml using 4-μl sitting drops at a ratio of 1:1 protein/mother liquor [0.1 M Hepes, 50 mM lithium sulfate, 10 mM EDTA, and 24 to 29.5% (w/v) polyethylene glycol (PEG) 600 (pH 7.0 to 8.4)]. Crystals grown at 4°C typically appeared after 1 to 3 days and continued to grow to full size in approximately 2 weeks. Collected crystals were first cryoprotected by soaking in 0.1 M Hepes (at a pH identical to the crystal growth condition), 50 mM lithium sulfate, 10 mM EDTA, and 32% PEG 600. Collected crystals were typically ~650 μm × 400 μm × 300 μm.

X-ray data collection, structure determination, and refinement of the mouse P-gp/PBDE-100 cocrystal structure

X-ray diffraction data were collected at 100 K at the Canadian Light Source (08ID-1). X-ray fluorescence spectra were collected to verify the presence of bromine in P-gp crystals soaked with PBDE-100, and subsequent multiwavelength anomalous dispersion scans were conducted around the Br-K edge to maximize their anomalous signal contribution during data collection. All diffraction data were indexed and integrated with MOSFLM (78), processed with AIMLESS (79), and truncated with TRUNCATE within the CCP4 suite of programs (80). Phases were determined using the recently improved model of P-gp (26) as a search model in molecular replacement in Phaser (81). The resultant model underwent rigid-body and restrained positional refinement, with hydrogens applied in their riding positions, using PHENIX.REFINE (82) against a maximum likelihood target function with grouped B factors and secondary structure restraints. Rounds of refinement were interspersed with manual inspection and correction against Sigma-A–weighted electron density maps in Coot (83), and improvements to model geometry and stereochemistry were monitored using MolProbity (84). The side chains of residues proximal to PBDE-100 were modeled during the final rounds of refinement to avoid biasing their placement. Ligand description dictionaries were calculated using PHENIX.ELBOW (82), and the crystallographic position of PBDE-100 was validated using the anomalous scattering from the bromine atoms (Fig. 3C). The refined structure was judged to have excellent geometry as determined by MolProbity (84). The resulting refinement statistics are listed in table S1. Figures displayed in this paper were prepared using PyMOL ( Atomic coordinates and structure factors for the cocrystal structure presented here were deposited with the Protein Data Bank (accession code 4XWK). Sequence alignments were performed using ClustalW2 Version 2.1 [European Molecular Biology Laboratory–European Bioinformatics Institute (EMBL-EBI)], and molecular interfaces were examined using PDBePISA (EMBL-EBI).

Tuna collection and pollutant data analysis

Yellowfin tuna (T. albacares) were collected from the GOM (28°29′N, 89°27′W). Tuna muscles were sent to AXYS (AXYS Analytical Services Ltd.) for analysis. Pollutants were measured according to the following EPA (U.S. Environmental Protection Agency) methods: for OCPs, 608, 625, 1625, 8081, and 8270; for PBDEs, 1614 and 625; and for PCBs, 1668 and 8270. For PFCs, methods developed by AXYS were used (AXYS MLA-041, AXYS MLA-042, AXYS MLA-043, and AXYS MLA-060). Analyses of pollutant data were performed according to previous studies (35, 85). Where pollutant values were below the detection limit, the values were treated as nondetectable (ND). Most blank measurements were at or below the detection limits and thus not subtracted from samples. In runs where blank values were above detection limits, only samples with values more than twice the blank value were reported (after subtraction of the blank value). The remaining values were treated as ND. The “ND” values were assigned zero when calculating the molar concentrations of POPs in the lipid fraction of the muscle tissue and for further statistical analysis, so as not to overestimate the level of pollutant in samples.

All tuna pollutant data were converted into micromolar lipid-normalized concentrations. Because the lipid composition of the tuna tissue samples was unknown, we assumed a lipid mass density of 1 for volume conversion (that is, 1 g of lipid is equal to 1 ml). To meet assumptions of normality and homogeneity of variance for statistical analyses, all data were log10-transformed. Zero values were excluded from the analyses. To avoid overestimating pollutant levels, one of nine fish initially collected from the GOM was treated as an outlier because of its high levels (more than two times the interquartile range) that did not appear to represent the concentration signature for the group. All data were analyzed and plotted using Origin software.

Vesicle inhibition assays and data analysis

Inhibition of the efflux of probe substrates into membrane vesicles expressing P-gp was carried out according to the manufacturer’s instructions, with some modifications (Sigma-Aldrich). Briefly, vesicle membrane suspensions were added to a 96-well plate stored on ice. Incubation media containing the single POP congener, the POP mixture, or verapamil as positive control were added to the plate and incubated for 15 min. Substrate solutions containing either Mg-ATP or Mg-AMP (adenosine monophosphate) and probe substrate (NMQ) were added to the plate for the designated time (3 min). The incubation was ended by the addition of the ice-cold wash mix. The sample solution was transferred to a filter plate and washed five times with wash mix. Plates were allowed to dry at room temperature for approximately 1 hour, after which a 50:50 methanol/water solution containing the internal standard was added to the filter plate wells and incubated for 15 min before it was transferred to a 96-well analytical plate. For each inhibitor, three assays were performed in duplicate. IC50 values were determined from the average percent inhibition values from each experiment ± SE of the measurement.

Analyst Instrument Control and Data Processing Software (AB SCIEX, version 1.6.1) was used to analyze the unlabeled probe substrate in P-gp vesicle assays (NMQ) for data collection and integration, which were then processed with Microsoft Excel 2007 (Microsoft). Calibration standards were used to calculate concentration on the basis of analyte/internal standard peak-area ratios with Analyst Instrument Control and Data Processing Software (AB SCIEX, version 1.6.1). A Shimadzu API 4000 mass spectrometer in positive mode (4500 V) was used in tandem with a Waters Atlantis (dC18, 5 μm, 100 × 2.1 mm) column (at 40°C) and a Luna C8 guard column (4.0 × 2.0 mm) for separation in a mobile phase of 0.2% formic acid in water and 0.2% formic acid in methanol and an injection volume of 1 μl. The mobile phase flow rate was 0.6 ml/min and the mass transitions used to identify NMQ were m/z (mass/charge ratio) = 339.1 and 339.3. Deuterated NMQ (d3-NMQ) was used as an internal standard and identified with mass transitions of m/z = 342.1 and 342.3.


Supplementary material for this article is available at

fig. S1. ATPase activity of purified, recombinant mouse P-gp.

fig. S2. Functional expression of mouse P-gp in yeast cells.

fig. S3. Comparison of PBDE-100 interacting residues.

fig. S4. Amino acid sequence alignment of human, mouse, zebrafish, and sea urchin P-gp.

fig. S5. Levels of POPs in wild-caught yellowfin tuna.

fig. S6. Individual POP congeners from the mixture inhibiting human P-gp.

table S1. Data collection and refinement statistics of the mouse P-gp/PBDE-100 cocrystal structure.

table S2. Metadata on yellowfin tuna specimens used in this study.

table S3. Physical and chemical properties of the 10 POP inhibitors.

References (8994)

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.


Acknowledgments: We thank G. Moy for technical assistance, R. Doshi and I. L. Urbatsch for helpful discussions on ATPase assays, and R. J. Ferl and J. D. Rine for the yeast strain and the p415GPD plasmid. We thank B. Zgliczynski for providing the GOM map. We thank D. Epel, V. D. Vacquier, and J. Stegeman for discussions of earlier versions of the manuscript. Funding: This work was supported by the program on Oceans and Human Health through NIH grant ES021985 and NSF grant 1314480 to A.H. and G.C., by a WAITT Foundation grant to A.H. and S.S., and by a University of California San Diego Academic Senate Grant to A.H. A.P.M. was supported by a National Health and Medical Research Council CJ Martin Postdoctoral Research Fellowship. Author contributions: S.C.T.N. designed the research, conducted the assays, analyzed the data, and co-authored the manuscript. S.D.R. and A.P.M. conducted x-ray and protein purification work, and edited the manuscript. T.G. designed and performed the yeast assays, and edited the manuscript. L.T.B. collected tuna samples and analyzed the data. L.M.V. performed vesicle assays and analyzed the data. C.C. conducted protein purification work. G.L., S.S., and G.C. designed and supervised the research, and edited the manuscript. A.H. analyzed the data, designed the research, co-authored the manuscript, and supervised the project. Competing interests: The authors declare that they have no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors. The GOM map file was created using ArcGIS Desktop v10.3.1 and is licensed under the Educational Site License.
View Abstract

Stay Connected to Science Advances

Navigate This Article