Abstract
P-glycoprotein (Pgp) is known for its dichotomous roles as both a safeguarding efflux transporter against xenobiotics and as a catalyst for multidrug resistance. Given the susceptibility of numerous therapeutic compounds to Pgp-mediated resistance, compliance with Food and Drug Administration (FDA) guidelines mandates an in-depth in vitro transport assay during drug development. This study introduces an innovative transport assay that aligns with these regulatory imperatives but also addresses limitations in the currently established techniques. Using Pgp-reconstituted liposomes and employing surface plasmon resonance (SPR), this study developed a distinct method of measuring the relative transport rates of Pgp substrates in a controlled microenvironment. The Pgp substrates selected for this study—quinidine, methadone, and desipramine—resulted in transport ratios that corroborate with trends previously observed. To assess the kinetics of Pgp-mediated transport, the results were analyzed by fitting the data to both currently proposed Pgp substrate translocation models—the vacuum cleaner and flippase models. While the resulting kinetic analysis in this study lends support predominantly to the vacuum cleaner model, this study most notably developed a novel method of assessing Pgp-mediated transport rates and real-time kinetics using surface plasmon resonance.
1 Introduction
P-glycoprotein (Pgp), a transmembrane transporter, functions as a protective mechanism within the human biological system by expelling xenobiotics from tissues. However, Pgp’s involvement in multidrug resistance in cancer and other diseases is well-documented. Most small molecule therapeutics are Pgp substrates, leading to their efflux by Pgp at therapeutic doses and thereby complicating treatment. Efforts to counter this challenge involve designing non-Pgp substrate small molecules. A fundamental step in this endeavor is measuring the transport rate of Pgp for specific compounds and determining their status as Pgp substrates. Given the critical role of assessing Pgp substrate characteristics in novel compound development, the Food and Drug Administration (FDA) issued a guideline in January 2020 named “In vitro drug interaction studies-cytochrome P450 enzyme-and transporter-mediated drug interactions guidance for industry”. This guideline outlines acceptable in vitro systems and methodologies, including the use of Caco-2 cells, commercial or in-house vesicles, knock-out/down cells, and transfected cells, for investigating Pgp transport rates ().
Presently, three widely used methods adhere to the outlined guidelines for quantifying Pgp efflux. The initial approach involves a vesicle transport assay utilizing inside-out vesicles (; ; ; ; ). This method gauges the substance’s permeation into vesicles when exposed to Pgp and a Pgp inhibitor. A quantifiable substrate, often radioactively or fluorescently labeled, is employed (; ; ). Advanced quantification may also involve mass spectrometry (), although both quantification techniques follow the assay endpoint and might not precisely capture real-time Pgp transport. The second method entails Pgp overexpression in systems like Caco-2 or MDCK cells (; ; ). Dubbed the drug resistance/accumulation test, it observes how chemicals penetrate and impact cell viability in the presence of Pgp or its inhibitor. A discernible split in the cell viability curve between Pgp inhibitor absence and presence designates a Pgp substrate (; ). However, this method’s drawback lies in its use of a mammalian cell system containing diverse transporters, yielding cumulative effects rather than strict Pgp specificity. The third method for Pgp efflux measurement employs a transcellular system utilizing Pgp overexpressing cells (; ). A cell monolayer is established on one side of the transwell system, assessing compound transport across to the other side. Like the previous method, this approach may encounter interference and cross-reactivity from other cellular transporters.
In addition to efflux rate, understanding the transport kinetics of a specific compound through Pgp is vital for comprehending its transport mechanism, ultimately aiding in improved drug design to circumvent Pgp. However, a more comprehensive grasp of Pgp’s transport process kinetics is required. Existing studies on Pgp kinetics predominantly relate to substrate-coupled ATP hydrolysis and employ Michaelis-Menten analysis (; ). While valuable for comparing how diverse substrates influence ATP hydrolysis in Pgp, this analysis typically yields substrate Km and Vmax values but falls short of revealing the complete transport mechanism. Furthermore, certain biochemical and biophysical investigations delve into the kinetics of substrate partitioning into the membrane preceding Pgp interaction (; ). A comprehensive examination of the full transport process kinetics, encompassing substrate-membrane and substrate-Pgp interactions, still merits deeper investigation.
This study aimed to employ surface plasmon resonance (SPR) to gauge the efflux rate of Pgp and elucidate the transport kinetics of Pgp with distinct substrates. The SPR-based transport assay aligns with FDA guidelines for in vitro transporter investigations and addresses limitations of current methods. Additionally, the transport kinetics analysis via SPR offers insights into both substrate-membrane and substrate-Pgp interactions, thereby enhancing our understanding of Pgp’s transport mechanism.
2 Materials and methods
2.1 Materials
Quinidine was sourced from ThermoFisher Scientific (Waltham, MA), while methadone hydrochloride came from the U. S. Pharmacopeia (Rockville, MD). Rhodamine B, desipramine hydrochloride, and ADP sodium were procured from Sigma-Aldrich (St. Louis, MO). Cholesterol, tris-HCl, and disodium ATP were obtained from Amresco (Solon, OH). Dithiothreitol (DTT) was purchased from Gold Biotechnology (Olivette, MO), and ethylene glycol tetraacetic acid (EGTA) and imidazole were sourced from Alfa Aesar (Tewksbury, MA). Escherichia (E.) coli total extract lipid and E. coli polar extract lipid were acquired from Avanti Polar Lipids Inc. (Alabaster, AL). The detergent n-dodecyl-β-D-maltoside (DDM) was procured from MilliporeSigma (Formerly, EMD Millipore Corporation) (Burlington, MA), and 4-(2-hydroxyethyl)−1−piperazineethanesulfonic acid (HEPES) and acrylamide from MilliporeSigma (Formerly, Calbiochem) (Burlington, MA). Sodium orthovanadate (Na3VO4) was sourced from Enzo Life Sciences (Farmingdale, NY). All other chemicals were acquired from ThermoFisher Scientific (Waltham, MA).
2.2 Expression, purification, and reconstitution of Pgp
The his-tagged wild-type mouse Pgp (Abcb1a, MDR3) was purified from Pichia (P.) pastoris using affinity chromatography with nickel-nitrilotriacetic acid (Ni-NTA) (Thermo Fisher Scientific) followed by diethylaminoethyl cellulose (DEAE) resin (Thermo Fisher Scientific) as previously described (; ). Briefly, Pgp solubilized by DDM was reconstituted into 200 nm liposomes using a previously established procedure so that the orientation of the transporter was inside-out with nucleotide-binding domains projected outside of the liposome (; ). Liposomes were prepared using 80% wt/vol E. coli Polar Lipid Extract (Avanti Polar Lipids) and 20% wt/vol cholesterol. A thin lipid film was created by evaporating the mixture of lipid extract and cholesterol dissolved in chloroform using a rotorvap system for 1 h. Rehydration of the dried lipid film in an aqueous phase (0.1 mM EGTA and 50 mM Tris-HCl) was followed by ten cycles of freeze-thaw in liquid nitrogen, generating liposomes of various sizes. The unilamellar liposomes with an average size of 200 nm were created by extruding the multilamellar liposome eleven times through a 400 nm filter with a LIPEX extruder. Bigger liposome size might lower SPR signal and decrease the signal sensitivity because analytes used are relatively small compared to ligand. Detergent-solubilized Pgp was dialyzed in HEPES buffer (20 mM HEPES, 100 mM NaCl, 5 mM MgCl2, 2 mM DTT, pH 7.4) for 2 h to remove DDM and glycerol. Dialyzed-protein was incubated with the extruded liposome with a lipid-to-protein ratio of 6.25 μM Pgp (mg ml−1 lipid)−1 for an hour and dialyzed for another 2 h to make proteoliposomes. The concentration of protein reconstituted in the liposomes was determined using a DC Protein Assay Kit II (Bio-Rad, Hercules, CA).
2.3 Vesicle transport assay
The ability of Pgp-embedded liposomes to transport pharmaceutical agents is conducted using a dialysis tube in Chifflet buffer (150 mM NH4Cl, 5 mM MgSO4, 0.02% wt/vol NaN3, 50 mM Tris-HCl, pH 7.4) (). Briefly, 1 mg/mL of liposome (control) or proteoliposome suspension in Chifflet buffer was transferred into narrow dialysis tubing (3.5 kDa MWCO), which was then placed into 50 mL centrifuge tubes containing 20 mL of Chifflet buffer added with 2 mg of rhodamine B and with or without 5 mM of ATP. Then, the tubes were incubated at 37 °C for up to 2 h. At time points 0, 30 min, 60 min, and 120 min, 200 μL of the buffer outside the dialysis tubing was collected, and 200 μL of fresh buffer was added to the tubes to maintain a total of 20 mL release media. The collected samples had their absorbances tested at 550 nm using UV-vis spectroscopy. The study was carried out in triplicate, and samples were evaluated through a standard curve created successive 2-fold dilutions of free rhodamine B (R2 = 0.9987) to calculate the ratio of absorbance to rhodamine B concentration.
After calculating rhodamine B concentrations using absorbance through the standard curve, the transport efficiency was calculated using the following equation. Briefly, concentration (I) is the initial concentration of rhodamine B in the buffer, while concentration (S) is the concentration of rhodamine B in the buffer at each time point. Reduced values indicate the drug being transported into the proteoliposome.
2.4 ATPase assay
ATPase assays were carried out to measure the ATP hydrolysis of Pgp in the presence of three different substrates, including desipramine, methadone, and quinidine. The Chifflet method, which detects the formation of inorganic phosphate following ATP hydrolysis, was used as described previously (; ; ). The ATPase activity of the three drugs was measured in the presence of 50 nM Pgp reconstituted in liposomes in Chifflet buffer.
The ATPase activity of desipramine, methadone, and quinidine were fit by nonlinear regression using Igor Pro 6.2 software (Wavemetrics, Tigard, OK) as described previously (; ; ). The ATP hydrolysis kinetics of desipramine and methadone fit with the monophasic equation derived from the Michaelis-Menten equation, where v is the ATP hydrolysis rate, Vmax is the maximum ATP hydrolysis rate at saturating drug, [L] is the ligand concentration, Km is the Michaelis-Menten constant, and vbasal is the basal ATPase activity in the absence of the drug.
Unlike the ATP hydrolysis kinetics of methadone and desipramine, quinidine established a biphasic ATPase activity profile. Substrate inhibition was used for biphasic kinetics, where is the inhibitory constant.
2.5 Surface plasmon resonance
Surface plasmon resonance is typically used to measure the real-time kinetics of the interaction between two molecules. The classic SPR monitors kinetic events by measuring the change in refractive angle of the incident light upon changes on the sensor surface (). The sensor in the classic SPR system is normally made of gold, which exhibits a broad plasmonic field. The SPR system we used in this study is a localized SPR (LSPR) from Nicoya Lifescience (Canada) instead of a classic SPR. Localized SPR monitors the change in absorbance wavelength at the local point where binding events occur instead of the change in the refractive angle (; ). The sensor of the LSPR system is made of nanogold, which has a smaller plasmonic field than a gold film and is therefore sensitive to localized changes instead of bulk changes. LSPR is known to have similar sensitivity for biomolecular binding events to SPR but is less susceptible to drift due to bulk changes ().
The SPR experiment was carried out with the three drugs, desipramine, methadone, and quinidine, in the presence of 5 mM ATP, 5 mM ADP, or no nucleotide as analytes shown in Figure 1. Analyte concentrations used were 500, 250, 125, 62.5, and 31.25 . Pgp reconstituted in liposomes was used as the ligand and were immobilized onto the liposome or NTA sensor at a saturated level, which typically takes two injections of the ligand at 20 μM of the protein. The Liposome sensor has a hydrophobic chain to capture the bilayer portion of the proteoliposome, while the NTA sensor captures the His tag on Pgp. The running buffer contained 100 mM Tris, 150 mM NaCl, 5 mM MgCl2, 0.5% DMSO, and pH 7.4. The running buffer contained the same amount of nucleotide (ATP, ADP, or no nucleotide) as in the analyte solution. The association phase was set to 1,200 s and the dissociation phase was set to 800 s. All the experiments were performed at a speed of 5 . Regeneration of the liposome sensor was carried out using 20 mM CHAPs, while the regeneration of the NTA sensor was carried out using 300 mM EDTA and 250 mM imidazole. Double referencing was carried out and subtracted from the data before further analysis.
FIGURE 1
2.6 Data fitting using COPASI
COPASI is the software used previously by other studies to fit complex kinetics data. SPR data with multiple reactions and compartments were also successfully modeled using COPASI (; ). Therefore, we used this software to fit the transport data obtained from the SPR experiment. The data were fit with two models, the vacuum cleaner and the flippase model. Each model consists of five different compartments, including bulk solvent, the sensor’s surface, the liposome’s outer leaflet, the liposome’s inner leaflet, and the liposome’s inner space. There are six reactions used in each model to represent the whole transport process. Specifically, in the vacuum cleaner model, the following reaction was used: (): analyte diffused from the bulk to the surface of the sensor, (), analyte diffused from outside to outer leaflet of the proteoliposome, (), analyte bound to Pgp, (), Pgp transported analyte into the inner space of the proteoliposome, (), analyte diffused from the outer leaflet to the inner leaflet of the proteoliposome, (), analyte diffused from the inner leaflet to the inner space of the proteoliposome. The flippase model used similar reactions except for reaction (), in which Pgp transported the analyte to the inner leaflet instead of the inner space of the proteoliposome.
The global quantities were defined as the SPR output signal, represented by the following equation:
Where is the relative unit in SPR sensorgram, is the overall scaling, is the baseline constant, , , , and represent the time-dependent evolution of the analytes incorporated into the inner space of the proteoliposome, the outer leaflet of the proteoliposome, the Pgp, and the inner leaflet of the proteoliposome, accordingly. A time course kinetics simulation of the SPR signal was carried out with intervals of 50,000 and an interval size of 0.06 s for a duration of 3,000 s using the deterministic method. Following the simulation, parameter estimation was used to fit the experimental SPR data. All available convergence methods were implemented, but the best result was obtained with an evolutionary programming algorithm. Data with different concentrations of the analytes were fit globally with the kinetics for each part of Equation 4, varied freely. The scaling factor and baseline constant were fitted locally for each concentration. As indicated in the previous AFM study, the amount of Pgp reconstituted into the liposome has 80% inside-out orientation (). Only Pgp with the nucleotide binding domains projecting out of the proteoliposome has access to ATP and can efflux substrates. Therefore, a restraint was applied to the number of Pgp involved in active transport during simulation.
3 Results
3.1 Vesicle transport assay as a standard transport method showed that the Pgp substrate was transported inside the proteoliposome by Pgp
Utilizing rhodamine B, a well-established Pgp substrate (; ), the vesicle transport assay was conducted. The transport efficiency of rhodamine B was depicted in Figure 2, illustrating rates at 10, 30, 60, and 120 min into the assay. Rhodamine B is recognized for its Pgp-mediated membrane transport facilitated by ATP (). Our findings demonstrated a substantial increase in rhodamine B accumulation within proteoliposome in the presence of ATP compared to its absence. Notably, within 10 min, rhodamine B uptake without ATP was around 5%, whereas with 5 mM ATP, it escalated to 10%, with statistical significance. Over 120 min, ATP-assisted transport raised vesicular rhodamine B content to 17%, while this remained at 12% without ATP. Rhodamine uptakes in proteoliposome in the absence of ATP is due to passive diffusion. These outcomes highlight successful Pgp-mediated transport of rhodamine B within proteoliposome with ATP support, validating the functionality of the employed inside-out vesicle system. ATPase of methadone, quinidine, and desipramine confirms the activity of the transporter.
FIGURE 2
To further assess proteoliposome activity, an ATPase assay was conducted using three additional Pgp substrates: methadone, quinidine, and desipramine. Figure 3 depicts the ATP hydrolysis kinetics curve for each substrate across varying concentrations. The ATPase data reveals that, in the presence of these drugs, Pgp’s basal ATP hydrolysis rate measured 482 ± 50 nmolmin−1mg−1, consistent with prior investigations (; ). The presence of sodium orthovanadate, a non-competitive ATP inhibitor, resulted in a decline of Pgp’s ATPase activity to approximately 100 nmolmin−1mg−1, corroborating findings from previous studies (; ).
FIGURE 3
Methadone incited monophasic ATP hydrolysis activation in Pgp, revealing kinetics with a peak of approximately 800 nmolmin−1mg−1 at saturating methadone levels. Curve fitting yielded and values of 849 ± 45 nmolmin−1mg−1 and 28 ± 10 μM, respectively, consistent with a prior methadone-Pgp study (
In contrast to methadone and desipramine, quinidine induced Pgp’s ATP hydrolysis in a biphasic manner, implying the presence of at least two quinidine binding sites on Pgp, as observed with other substrates like verapamil, lidocaine, and loperamide. The kinetics demonstrated a peak at approximately 900 nmolmin−1mg−1 at a quinidine concentration of 15.6 μM. Fitting the biphasic kinetics curve to Equation 3 yielded , , and values of 626 ± 44 nmolmin−1mg−1, 2.9 ± 1 μM, and 131 ± 34 μM, respectively. Notably, the and values obtained in this study were consistent with those reported in a prior investigation involving quinidine and Pgp (
3.2 SPR signal amplitude correlated accordingly to the transport rate of each Pgp substrate
The SPR signal amplitude corresponds to alterations in mass on the sensor surface (
Experiments involving three distinct Pgp substrates consistently demonstrated that the SPR signal amplitude in the presence of 5 mM ATP significantly surpassed the amplitude in the absence of nucleotides or with ADP (Figure 4). The ratio between SPR amplitude in the presence of ATP and that without nucleotides is an indicator of each substrate’s relative transport rate. Table 1 presents the ratios for quinidine, methadone, and desipramine as 1.64 ± 0.12, 1.31 ± 0.1, and 1.30 ± 0.08, respectively. These ratios align with prior investigations that assessed Pgp transport rates of these three compounds using a cell-based system. Mahar et al. and Tournier et al. reported Pgp efflux rates of 2 for desipramine and methadone (
FIGURE 4

The SPR signals of desipramine, methadone, and quinidine in the presence of 5 mM ATP, no nucleotide, and 5 mM ADP. The association phase is 1200s, and the dissociation phase is 800s, both injected at a 5 μL/min speed. The SPR sensorgrams for each Pgp substrate comprised 5 different substrate concentrations: 500, 250, 125, 62.5, and 31.25 μM. (A) is the SPR signal of desipramine with 5 mM ATP, (B) is desipramine without nucleotide, and (C) is desipramine with 5 mM ADP. (D) is the SPR signal of methadone with 5 mM ATP, (E) is methadone without nucleotide, and (F) is methadone with 5 mM ADP. (G) is the SPR signal of quinidine with 5 mM ATP, (H) is quinidine without nucleotide, and (I) is quinidine with 5 mM ADP. Dashed line represents start of association phase and dissociation phase. Figure shows the representative curve of each experiment, the number of replicates for each experiment is n=5.
TABLE 1
| Literature | Experimental | Methods used to measure Pgp efflux ratio in literature | |
|---|---|---|---|
| Methadone | 1.9 (MDR-MDCK) 2.4 (Caco-2) | 1.31 ± 0.1 | Transwell system with Caco-2 or MDR1-MDCK cells |
| Quinidine | 10 | 1.64 ± 0.12 | Transwell system with Caco-2 or MDR1-MDCK cells |
| Desipramine | 2 | 1.30 ± 0.08 | Drug accumulation assay with MDR1-MDCKII cells |
The relative transport rate of methadone, desipramine, and quinidine was measured using SPR and extracted from the literature. n=5. Error bar is ±SD.
3.3 The kinetics of SPR transport of desipramine supported the vacuum cleaner model for the Pgp transport mechanism
Recognized for its capability to track real-time binding kinetics, SPR serves as a means to not only gauge Pgp’s relative transport rate across diverse compounds but also to delve into the transport kinetics of each compound (
FIGURE 5

COPASI fitting of SPR signal in the presence of 5 mM ATP at 500 μM of the substrates with two different models: vacuum cleaner and flippase. The black line is the experimental SPR data. The red line is the fitted line. The dashed line is the residual of the fit. (A) is the fitting for desipramine with vacuum cleaner model, (B) is the fitting for desipramine with flippase model, (C) is the fitting for methadone with vacuum cleaner model, (D) is the fitting for methadone with flippase model, (E) is the fitting for quinidine with vacuum cleaner model, (F) is the fitting for quinidine with flippase model. Figure shows the representative curve of each experiment, the number of replicates for each experiment is n=5.
The kinetics of each reaction encompassed in the transport process of the three drugs were also assessed, as outlined in Table 2. The values of kon and koff for the reaction involving substrate diffusion from the exterior to the outer leaflet exhibited remarkable similarity for desipramine (within the 10−6 M−1s−1 and 10−6 s−1 range), suggesting a relatively balanced rate for this bidirectional process in desipramine. Conversely, in the case of methadone, the on-rate (5.07 x 10−5 M−1s−1) for diffusion from the exterior into the outer leaflet slightly exceeded the off-rate (8.34 x 10−6 s−1), implying a slight bias against methadone’s diffusion into the outer leaflet. For quinidine, the kon for diffusion into the outer leaflet of the proteoliposome was 4.62 x 10−8 M−1s−1, while the koff was 1.20 x 10−5 s−1, indicating a greater likelihood for quinidine to stay in the exterior than to diffuse into the outer leaflet of the bilayer.
TABLE 2
| Desipramine | Methadone | Quinidine | ||||
|---|---|---|---|---|---|---|
| Kon (M−1s−1) | Koff (s−1) | Kon (M−1s−1) | Koff (s−1) | Kon (M−1s−1) | Koff (s−1) | |
| Diffuse from outside to outer leaflet | 1.15 x 10−6 | 4.1 x 10−6 | 5.07 x 10−5 | 8.34 x 10−6 | 4.62 x 10−8 | 1.20 x 10−5 |
| Binding to Pgp | 8.9 x 10−3 | 129 | 1.61 x 10−6 | 6.39 x 105 | 6.69 x 10−4 | 823 |
| Pgp transport to inside the proteoliposome | 3.57 x 10−5 | 6.44 x 10−8 | 2.75 x 10−4 | 1.89 x 10−4 | 1.6 x 10−3 | 7.91 x 10−8 |
| Diffuse from outer leaflet to inner leaflet | 4.00 x 10−5 | 2.68 | 1.83 x 10−7 | 8.81 x 10−7 | 0.01 | 164 |
| Diffuse from inner leaflet to inside the proteoliposome | 1.10 x 10−5 | 1.09 | 0 | 493 | 1.1 x 10−3 | 0.13 |
The kinetics from the SPR curve fitted to the vacuum cleaner model for three Pgp substrates: desipramine, methadone, and quinidine. n=5.
In the reactions involving desipramine’s diffusion from the outer leaflet to the inner leaflet and the departure of substrate from the outer leaflet to bind to Pgp, the kon was significantly lower than the koff. This observation suggests a preference for the desipramine-outer leaflet complex over the desipramine-inner leaflet and desipramine-Pgp complex. Similar trends were observed with both methadone and quinidine, indicating that all three drugs displayed a preference for complex formation with the outer leaflet rather than with Pgp or the inner leaflet of the bilayer.
Remarkably, the reactions in which Pgp facilitated the transport of substrates inside the vesicle exhibited significantly smaller koff values (6.44 x 10−8 s−1 for desipramine, 1.89 x 10−4 s−1 for methadone, and 7.91 x 10−8 for quinidine) compared to the corresponding kon values (3.57 x 10−5 M−1s−1 for desipramine, 2.75 x 10−4 M−1s−1 for methadone, and 1.6 x 10−3 M−1s−1 for quinidine), indicating a pronounced preference for substrate transportation into the proteoliposome interior. In the final reaction, involving substrates diffusing from the inner leaflet to the inner space of the proteoliposome, the kon rates were significantly lower than the koff rates for all three drugs. This observation suggests that all three substrates preferred the hydrophobic inner leaflet over the hydrophilic core of the proteoliposome.
4 Discussion
Current methods, such as vesicle transport assay, drug resistance/accumulation assay, and transwell assay, lack the capability to monitor real-time substrate transport and isolate the specific transport rate mediated solely by Pgp from the combined effects of multiple transporters present in the cell membrane. To address these limitations and adhere to FDA guidelines for investigating efflux transporters (
The relative transport rates of methadone and quinidine were assessed under conditions of ATP, ADP, and no nucleotide. The signal amplitudes differed significantly between the presence of 5 mM ADP and the absence of nucleotide. Additionally, the signal amplitude with ADP was notably lower than without nucleotide, suggesting distinct Pgp functionality with ADP compared to no nucleotide. This discrepancy is likely the result of ADP inducing a specific Pgp conformation that hinders substrate expulsion, as previous research proposed that an inward-facing Pgp-ADP conformation is less conducive to efficient efflux than the outward-facing conformation (
In addition to assessing Pgp transport rate, the SPR technique also offers valuable insights into the real-time kinetics of Pgp transport. This not only enhances the comprehension of Pgp-mediated transport but also lends support to existing models elucidating the mechanism behind Pgp-mediated transport. The two models under consideration in this study are the vacuum cleaner and flippase models, both of which propose distinct pathways through which Pgp expels substrates. The vacuum cleaner model suggests that Pgp effluxes substrates from the inner leaflet of the lipid bilayer to the extracellular space, while the flippase model proposes expulsion from the inner leaflet to the outer leaflet of the bilayer (
Both of these models find support in both biochemical and structural investigations. The vacuum cleaner model finds substantiation in the crystal structures of hPgp in an outward-facing conformation and mPgp in an inward-facing conformation (
Conversely, there exists evidence supporting the substrate flippase model. Investigations have indicated that Pgp can effectively transport a range of lipids, encompassing phosphatidylcholine, phosphatidylethanolamine, and sphingomyelin (
The kinetics data derived from the SPR transport assay for the three compounds, namely, desipramine, methadone, and quinidine, provided intriguing insights into the behavior of the Pgp transport system. Notably, all three compounds exhibited a clear preference for binding to the outer leaflet of the membrane compared to interacting with Pgp or diffusing into the inner leaflet. This observation suggests that the outer leaflet could play a significant role in driving the passive diffusion of chemicals into cells. However, once binding with Pgp occurred, the substrates were actively expelled into the core of the proteoliposome, indicating that Pgp possessed the ability to overcome the passive diffusion of these drugs.
The SPR transport method developed in this study is constrained by the localized SPR technology employed. This approach hinges on alterations in sensor absorbance when interactions transpire on the sensor surface. Consequently, it is unable to detect molecules with absorbance akin to that of the sensor, positioned around the 500 nm wavelength. Nonetheless, conventional SPR, reliant on changes in reflectance angle rather than absorbance, can identify these molecules. Furthermore, the study’s employment of mouse Pgp reconstituted in liposomes, as opposed to human Pgp, presents another limitation. Consequently, the delineated kinetics are applicable solely to mouse Pgp. However, given the substantial 80% sequence similarity between mouse and human Pgp, insights gleaned from investigations involving mouse Pgp still offer valuable understanding into the functioning of human Pgp.
Undoubtedly, this study has pioneered a novel approach for quantifying both Pgp efflux rate and Pgp efflux kinetics through the utilization of SPR. The primary objective was to align with the FDA guidelines while simultaneously addressing the limitations inherent in prevailing in vitro transport techniques. Most notably, this study stands as the pioneering effort to dissect the intricate kinetics of Pgp transport into distinct and detailed stages.
Statements
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Author contributions
PN: Writing–original draft, Writing–review and editing. SC: Writing–review and editing. AK: Writing–review and editing. AR: Writing–review and editing. AGR: Writing–original draft. MPX: Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Institutes of Health (R01CA204846, AGR and R01DK099596, MPX).
Acknowledgments
We give thanks to Ina L. Urbatsch of Texas Tech University Health Sciences Center for her generous gift of P. pastoris with the wild-type mouse Pgp transporter gene.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Abbreviations
AMPPNP, adenosine 5’-(β,γ-imido)triophsophate; CHAPS, 3-((3-cholamidopropyl) dimethylammonio)−1−propanesulfonate; COPASI, complex pathway simulator; DDM, n-dodecyl-β-D-maltoside; DEAE, diethylaminoethyl cellulose; FDA, Food and Drug Administration; KPi, potassium phosphate; NATA, N-acetyl-L-tryptophanamide; NBD, nucleotide-binding domain; Ni-NTA, nickel-nitrilotriacetic acid; Pgp, P-glycoprotein; Pi, inorganic phosphate; SPR, surface plasmon resonance.
References
1
AbulrobA. n. G.GumbletonM. (1999). Transport of phosphatidylcholine in MDR3-negative epithelial cell lines via drug-induced MDR1 P-glycoprotein. Biochem. Biophysical Res. Commun.262 (1), 121–126. 10.1006/bbrc.1999.1120
2
AmbudkarS. V.LelongI. H.ZhangJ.CardarelliC. (1998). “[36] Purification and reconstitution of human P-glycoprotein,” in Methods in enzymology (Academic Press), 292, 492–504. (ABC Transporters: Biochemical, Cellular, and Molecular Aspects; . Available from: https://www.sciencedirect.com/science/article/pii/S0076687998920389 (Accessed October 20, 2022).
3
BaiJ.SwartzD. J.ProtasevichI. I.BrouilletteC. G.HarrellP. M.HildebrandtE.et al (2011). A gene optimization strategy that enhances production of fully functional P-glycoprotein in Pichia pastoris. PLOS ONE6 (8), e22577. 10.1371/journal.pone.0022577
4
BorgniaM. J.EytanG. D.AssarafY. G. (1996). Competition of hydrophobic peptides, cytotoxic drugs, and chemosensitizers on a common P-glycoprotein pharmacophore as revealed by its ATPase activity. J. Biol. Chem.271 (6), 3163–3171. 10.1074/jbc.271.6.3163
5
BoschI.Dunussi-JoannopoulosK.WuR. L.FurlongS. T.CroopJ. (1997). Phosphatidylcholine and phosphatidylethanolamine behave as substrates of the human MDR1 P-glycoprotein. Biochemistry36 (19), 5685–5694. 10.1021/bi962728r
6
BoultonD. W.DeVaneC. L.ListonH. L.MarkowitzJ. S. (2002). In vitro P-glycoprotein affinity for atypical and conventional antipsychotics. Life Sci.71 (2), 163–169. 10.1016/s0024-3205(02)01680-6
7
BrumL. (2020). Vitro drug interaction studies - cytochrome P450 enzyme - and transporter - mediated drug interactions guidance for industry.
8
ClayA. T.SharomF. J. (2013). Lipid bilayer properties control membrane partitioning, binding, and transport of P-glycoprotein substrates. Biochemistry52 (2), 343–354. 10.1021/bi301532c
9
CollierJ. H. (2012). Driving a wedge into acute allergic responses. Sci. Transl. Med.4 (163), 163ec219. 10.1126/scitranslmed.3005421
10
CuiS.QiaoJ.XiongM. P. (2022). Antibacterial and biofilm-eradicating activities of pH-responsive vesicles against Pseudomonas aeruginosa. Mol. Pharm.19 (7), 2406–2417. 10.1021/acs.molpharmaceut.2c00165
11
DoanK. M. M.HumphreysJ.WebsterL.WringS.ShampineL.Serabjit-SinghC.et al (2002). Polli JWPassive permeability and P-glycoprotein-mediated efflux differentiate central nervous system (CNS) and non-CNS marketed drugs. J Pharmacol Exp Ther 303:1029-1037. J. Pharmacol. Exp. Ther.303, 1029–1037. 10.1124/jpet.102.039255
12
GameiroM.SilvaR.Rocha-PereiraC.CarmoH.CarvalhoF.BastosM. de L.et al (2017). Cellular models and in vitro assays for the screening of modulators of P-gp, MRP1 and BCRP. Molecules22 (4), E600. 10.3390/molecules22040600
13
GibbsM. E.WiltL. A.LedwitchK. V.RobertsA. G. (2018). A conformationally gated model of methadone and loperamide transport by P-glycoprotein. J. Pharm. Sci.107 (7), 1937–1947. 10.1016/j.xphs.2018.02.019
14
GlavinasH.MéhnD.JaniM.OosterhuisB.Herédi-SzabóK.KrajcsiP. (2008). Utilization of membrane vesicle preparations to study drug-ABC transporter interactions. Expert Opin. Drug Metabolism Toxicol.4 (6), 721–732. 10.1517/17425255.4.6.721
15
HelvoortA. v.GiudiciM. L.ThielemansM.MeerG. v. (1997). Transport of sphingomyelin to the cell surface is inhibited by brefeldin A and in mitosis, where C6-NBD-sphingomyelin is translocated across the plasma membrane by a multidrug transporter activity. J. Cell Sci.110 (Pt 1), 75–83. 10.1242/jcs.110.1.75
16
HutterE.FendlerJ. H. (2004). Exploitation of localized surface plasmon resonance. Adv. Mater.16 (19), 1685–1706. 10.1002/adma.200400271
17
KimB.EggelA.TarchevskayaS. S.VogelM.PrinzH.JardetzkyT. S. (2012). Accelerated disassembly of IgE:receptor complexes by a disruptive macromolecular inhibitor. Nature491 (7425), 613–617. 10.1038/nature11546
18
KimY.ChenJ. (2018). Molecular structure of human P-glycoprotein in the ATP-bound, outward-facing conformation. 6.
19
LedwitchK. V.GibbsM. E.BarnesR. W.RobertsA. G. (2016). Cooperativity between verapamil and ATP bound to the efflux transporter P-glycoprotein. Biochem. Pharmacol.118, 96–108. 10.1016/j.bcp.2016.08.013
20
Lerner-MarmaroshN.GimiK.UrbatschI. L.GrosP.SeniorA. E. (1999). Large scale purification of detergent-soluble P-glycoprotein fromPichia pastoris cells and characterization of nucleotide binding properties of wild-type, walker A, and walker B mutant proteins. J. Biol. Chem.274 (49), 34711–34718. 10.1074/jbc.274.49.34711
21
LiJ.JaimesK. F.AllerS. G. (2014). Refined structures of mouse P-glycoprotein. Protein Sci.23 (1), 34–46. 10.1002/pro.2387
22
LoetchutinatC.SaengkhaeC.Marbeuf-GueyeC.Garnier-SuillerotA. (2003). New insights into the P-glycoprotein-mediated effluxes of rhodamines. Eur. J. Biochem.270 (3), 476–485. 10.1046/j.1432-1033.2003.03403.x
23
MatsumotoT.KaifuchiN.MizuharaY.WarabiE.WatanabeJ. (2018). Use of a Caco-2 permeability assay to evaluate the effects of several Kampo medicines on the drug transporter P-glycoprotein. J. Nat. Med.72 (4), 897–904. 10.1007/s11418-018-1222-x
24
MatsunagaT.KoseE.YasudaS.IseH.IkedaU.OhmoriS. (2006). Determination of p-glycoprotein ATPase activity using luciferase. Biol. Pharm. Bull.29 (3), 560–564. 10.1248/bpb.29.560
25
MyszkaD. G. (1997). Kinetic analysis of macromolecular interactions using surface plasmon resonance biosensors. Curr. Opin. Biotechnol.8 (1), 50–57. 10.1016/s0958-1669(97)80157-7
26
NguyenP. H.SigdelK. P.SchaeferK. G.MensahG. a. K.KingG. M.RobertsA. G. (2020). The effects of anthracycline drugs on the conformational distribution of mouse P-glycoprotein explains their transport rate differences. Biochem. Pharmacol.174, 113813. 10.1016/j.bcp.2020.113813
27
OhashiR.WatanabeR.EsakiT.TaniguchiT.Torimoto-KatoriN.WatanabeT.et al (2019). Development of simplified in vitro P-glycoprotein substrate assay and in silico prediction models to evaluate transport potential of P-glycoprotein. Mol. Pharm.16 (5), 1851–1863. 10.1021/acs.molpharmaceut.8b01143
28
OmoteH.Al-ShawiM. K. (2002). A novel electron paramagnetic resonance approach to determine the mechanism of drug transport by P-glycoprotein. J. Biol. Chem.277 (47), 45688–45694. 10.1074/jbc.m206479200
29
O’ShannessyD. J. (1994). Determination of kinetic rate and equilibrium binding constants for macromolecular interactions: a critique of the surface plasmon resonance literature. Curr. Opin. Biotechnol.5 (1), 65–71. 10.1016/s0958-1669(05)80072-2
30
PatilA. G.D’SouzaR.DixitN.DamreA. (2011). Validation of quinidine as a probe substrate for the in vitro P-gp inhibition assay in Caco-2 cell monolayer. Eur. J. Drug Metab. Pharmacokinet.36 (3), 115–119. 10.1007/s13318-011-0046-9
31
PohlA.LageH.MüllerP.PomorskiT.HerrmannA. (2002). Transport of phosphatidylserine via MDR1 (multidrug resistance 1)P-glycoprotein in a human gastric carcinoma cell line. Biochem. J.365 (Pt 1), 259–268. 10.1042/bj20011880
32
RomsickiY.SharomF. J. (1999). The membrane lipid environment modulates drug interactions with the P-glycoprotein multidrug transporter. Biochemistry38 (21), 6887–6896. 10.1021/bi990064q
33
RomsickiY.SharomF. J. (2001). Phospholipid flippase activity of the reconstituted P-glycoprotein multidrug transporter. Biochemistry40 (23), 6937–6947. 10.1021/bi0024456
34
SarverJ. G.KlisW. A.ByersJ. P.ErhardtP. W. (2002). Microplate screening of the differential effects of test agents on Hoechst 33342, rhodamine 123, and rhodamine 6G accumulation in breast cancer cells that overexpress P-glycoprotein. SLAS Discov.7 (1), 29–34. 10.1177/108705710200700105
35
SharomF. J. (2014). Complex interplay between the P-glycoprotein multidrug efflux pump and the membrane: its role in modulating protein function. Front. Oncol.4, 41. 10.3389/fonc.2014.00041
36
SharomF. J.LiuR.RomsickiY.LuP. (1999). Insights into the structure and substrate interactions of the P-glycoprotein multidrug transporter from spectroscopic studies. Biochimica Biophysica Acta (BBA) - Biomembr.1461 (2), 327–345. 10.1016/s0005-2736(99)00166-2
37
SigdelK. P.WiltL. A.MarshB. P.RobertsA. G.KingG. M. (2018). The conformation and dynamics of P-glycoprotein in a lipid bilayer investigated by atomic force microscopy. Biochem. Pharmacol.156, 302–311. 10.1016/j.bcp.2018.08.017
38
SilvaR.PalmeiraA.CarmoH.BarbosaD. J.GameiroM.GomesA.et al (2015a). P-glycoprotein induction in Caco-2 cells by newly synthetized thioxanthones prevents paraquat cytotoxicity. Arch. Toxicol.89 (10), 1783–1800. 10.1007/s00204-014-1333-4
39
SilvaR.Vilas-BoasV.CarmoH.Dinis-OliveiraR. J.CarvalhoF.de Lourdes BastosM.et al (2015b). Modulation of P-glycoprotein efflux pump: induction and activation as a therapeutic strategy. Pharmacol. Ther.149, 1–123. 10.1016/j.pharmthera.2014.11.013
40
TakadaT.YamamotoT.MatsuoH.TanJ. K.OoyamaK.SakiyamaM.et al (2018). Identification of ABCG2 as an exporter of uremic toxin indoxyl sulfate in mice and as a crucial factor influencing CKD progression. Sci. Rep.8 (1), 11147. 10.1038/s41598-018-29208-w
41
TournierN.DeclèvesX.SaubaméaB.ScherrmannJ. M.CisterninoS. (2011). Opioid transport by ATP-binding cassette transporters at the blood-brain barrier: implications for neuropsychopharmacology. Curr. Pharm. Des.17 (26), 2829–2842. 10.2174/138161211797440203
42
ToyodaY.TakadaT.SuzukiH. (2019). Inhibitors of human ABCG2: from technical background to recent updates with clinical implications. Front. Pharmacol.10, 208. 10.3389/fphar.2019.00208
43
van HelvoortA.SmithA. J.SprongH.FritzscheI.SchinkelA. H.BorstP.et al (1996). MDR1 P-glycoprotein is a lipid translocase of broad specificity, while MDR3 P-glycoprotein specifically translocates phosphatidylcholine. Cell87 (3), 507–517. 10.1016/s0092-8674(00)81370-7
44
VerhalenB.DastvanR.ThangapandianS.PeskovaY.KoteicheH. A.NakamotoR. K.et al (2017). Energy transduction and alternating access of the mammalian ABC transporter P-glycoprotein. Nature543 (7647), 738–741. 10.1038/nature21414
45
WilletsK. A.Van DuyneR. P. (2007). Localized surface plasmon resonance spectroscopy and sensing. Annu. Rev. Phys. Chem.58 (1), 267–297. 10.1146/annurev.physchem.58.032806.104607
46
WiltL. A.NguyenD.RobertsA. G. (2017). Insights into the molecular mechanism of triptan transport by P-glycoprotein. J. Pharm. Sci.106 (6), 1670–1679. 10.1016/j.xphs.2017.02.032
Summary
Keywords
p-glycoprotein, surface plasmon resonace, efflux kinetics, efflux analysis, transport assay
Citation
Nguyen PH, Cui S, Kozarich AM, Rautio A, Roberts AG and Xiong MP (2024) Utilizing surface plasmon resonance as a novel method for monitoring in-vitro P-glycoprotein efflux. Front. Biophys. 2:1367511. doi: 10.3389/frbis.2024.1367511
Received
08 January 2024
Accepted
26 February 2024
Published
08 March 2024
Volume
2 - 2024
Edited by
Olga Vinogradova, University of Connecticut, United States
Reviewed by
Parkson Lee-Gau Chong, Temple University, United States
Nicholas James, University of Hawaii at Manoa, United States
Updates

Check for updates
Copyright
© 2024 Nguyen, Cui, Kozarich, Rautio, Roberts and Xiong.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: May P. Xiong, mpxiong@uga.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.