Binding of modulators to mouse and human multidrug resistance P-glycoprotein. A computational study

Study Models of Drug-Drug Interactions Involving P-Glycoprotein: The Potential Benefit of P-Glycoprotein Modulation at the Kidney and Intestinal Levels

P-glycoprotein (P-gp) is a crucial membrane transporter situated on the cell's apical surface, being responsible for eliminating xenobiotics and endobiotics. P-gp modulators are compounds that can directly or indirectly affect this protein, leading to changes in its expression and function. These modulators can act as inhibitors, inducers, or activators, potentially causing drug-drug interactions (DDIs). This comprehensive review explores diverse models and techniques used to assess drug-induced P-gp modulation. We cover several approaches, including in silico, in vitro, ex vivo, and in vivo methods, with their respective strengths and limitations. Additionally, we explore the therapeutic implications of DDIs involving P-gp, with a special focus on the renal and intestinal elimination of P-gp substrates. This involves enhancing the removal of toxic substances from proximal tubular epithelial cells into the urine or increasing the transport of compounds from enterocytes into the intestinal lumen, thereby facilitating their excretion in the feces. A better understanding of these interactions, and of the distinct techniques applied for their study, will be of utmost importance for optimizing drug therapy, consequently minimizing drug-induced adverse and toxic effects.

Effect of the Force Field on Molecular Dynamics Simulations of the Multidrug Efflux Protein P-Glycoprotein

Molecular dynamics (MD) simulations have been used extensively to study P-glycoprotein (P-gp), a flexible multidrug transporter that is a key player in the development of multidrug resistance to chemotherapeutics. A substantial body of literature has grown from simulation studies that have employed various simulation conditions and parameters, including AMBER, CHARMM, OPLS, GROMOS, and coarse-grained force fields, drawing conclusions from simulations spanning hundreds of nanoseconds. Each force field is typically parametrized and validated on different data and observables, usually of small molecules and peptides; there have been few comparisons of force field performance on large protein-membrane systems. Here we compare the conformational ensembles of P-gp embedded in a POPC/cholesterol bilayer generated over 500 ns of replicate simulation with five force fields from popular biomolecular families: AMBER 99SB-ILDN, CHARMM 36, OPLS-AA/L, GROMOS 54A7, and MARTINI. We find considerable differences among the ensembles with little conformational overlap, although they correspond to similar extents to structural data obtained from electron paramagnetic resonance and cross-linking studies. Moreover, each trajectory was still sampling new conformations at a high rate after 500 ns of simulation, suggesting the need for more sampling. This work highlights the need to consider known limitations of the force field used (e.g., biases toward certain secondary structures) and the simulation itself (e.g., whether sufficient sampling has been achieved) when interpreting accumulated results of simulation studies of P-gp and other transport proteins.

Structure activity relationships and the binding mode of quinolinone-pyrimidine hybrids as reversal agents of multidrug resistance mediated by P-gp

P-gp-associated multidrug resistance is a major impediment to the success of chemotherapy. With the aim of finding non-toxic and effective P-gp inhibitors, we investigated a panel of quinolin-2-one-pyrimidine hybrids. Among the active compounds, two of them significantly increased intracellular doxorubicin and rhodamine 123 accumulation by inhibiting the efflux mediated by P-gp and restored doxorubicin toxicity at nanomolar range. Structure-activity relationships showed that the number of methoxy groups, an optimal length of the molecule in its extended conformation, and at least one flexible methylene group bridging the quinolinone to the moiety bearing the pyrimidine favored the inhibitory potency of P-gp. The best compounds showed a similar binding pattern and interactions to those of doxorubicin and tariquidar, as revealed by MD and hybrid QM/MM simulations performed with the recent experimental structure of P-gp co-crystallized with paclitaxel. Analysis of the molecular interactions stabilizing the different molecular complexes determined by MD and QTAIM showed that binding to key residues from TMH 4-7 and 12 is required for inhibition.

Computational approaches in cancer multidrug resistance research: Identification of potential biomarkers, drug targets and drug-target interactions

Like physics in the 19th century, biology and molecular biology in particular, has been fertilized and enhanced like few other scientific fields, by the incorporation of mathematical methods. In the last decades, a whole new scientific field, bioinformatics, has developed with an output of over 30,000 papers a year (Pubmed search using the keyword "bioinformatics"). Huge databases of mass throughput data have been established, with ArrayExpress alone containing more than 2.7 million assays (October 2019). Computational methods have become indispensable tools in molecular biology, particularly in one of the most challenging areas of cancer research, multidrug resistance (MDR). However, confronted with a plethora of different algorithms, approaches, and methods, the average researcher faces key questions: Which methods do exist? Which methods can be used to tackle the aims of a given study? Or, more generally, how do I use computational biology/bioinformatics to bolster my research? The current review is aimed at providing guidance to existing methods with relevance to MDR research. In particular, we provide an overview on: a) the identification of potential biomarkers using expression data; b) the prediction of treatment response by machine learning methods; c) the employment of network approaches to identify gene/protein regulatory networks and potential key players; d) the identification of drug-target interactions; e) the use of bipartite networks to identify multidrug targets; f) the identification of cellular subpopulations with the MDR phenotype; and, finally, g) the use of molecular modeling methods to guide and enhance drug discovery. This review shall serve as a guide through some of the basic concepts useful in MDR research. It shall give the reader some ideas about the possibilities in MDR research by using computational tools, and, finally, it shall provide a short overview of relevant literature.

An Energetically Favorable Ligand Entrance Gate of a Multidrug Transporter Revealed by Partial Nudged Elastic Band Simulations

P-glycoprotein (P-gp) is a multidrug transporter, which harnesses the chemical energy of ATP to power the efflux of diverse chemotherapeutics out of cells and thus contributes to the development of multidrug resistance (MDR) in cancer. It has been proved that the ligand-binding pocket of P-gp is located at the transmembrane domains (TMDs). However, the access of ligands into the binding pocket remains to be elucidated, which definitely hinder the development of P-gp inhibitors. Herein, the access pathways of a well-known substrate rhodamine-123 and a cyclopeptide inhibitor QZ-Leu were characterized by time-independent partial nudged elastic band (PNEB) simulations. The decreasing free energies along the PNEB-optimized access pathway indicated that TM4/6 cleft may be an energetically favorable entrance gate for ligand entry into the binding pocket of P-gp. The results can be reconciled with a range of experimental studies, further corroborating the reliability of the gate revealed by computational simulations. Our atomic level description of the ligand access pathway provides valuable insights into the gating mechanism for drug uptake and transport by P-gp and other multidrug transporters.

•

P-gp contributes to the development of multidrug resistance in cancer.

•

The entrance of drugs into P-gp binding pocket has yet to be elucidated.

•

An energetically favorable entrance gate was revealed by PNEB simulations.

•

The computational results were reconciled with the experimental data.

•

The atomic simulations provide insights into the gating mechanism of P-gp.

Analogues of the Lignan Pinoresinol as Novel Lead Compounds for P-glycoprotein (P-gp) Inhibitors

To find novel P-gp-inhibitors, a library of pregnane X receptor (PXR) ligands and the ZINC DrugsNow library were superimposed on the P-gp inhibitor (+)-pinoresinol (1) used as a query for a three-dimensional similarity search. After determining the TanimotoCombo index of similarity with 1, eight compounds from the PXR library and two ZINC compounds were selected for biological evaluation. The P-gp inhibition study showed that compounds 7, 8, and 9 successfully increased intracellular doxorubicin (DOX) accumulation in the P-gp overexpressed Lucena 1 cells from 25, 12.5, and 6.25 μM, respectively. Among a series of analogues of 9, compounds 26-30 were shown to be active, with 26 and 27 causing a significant increase in DOX accumulation from 1.56 μM and rendering Lucena 1 sensitive to DOX from 1.56 and 0.78 μM, respectively. Molecular modeling studies showed that both compounds bind to the P-gp at transmembrane helices (TMH) 4, 5, and 6, with 27 also showing contacts with TMH 3.

Discovery of a non-toxic [1,2,4]triazolo[1,5-a]pyrimidin-7-one (WS-10) that modulates ABCB1-mediated multidrug resistance (MDR)

Multidrug resistance (MDR) has been shown to reduce the effectiveness of chemotherapy. Strategies to overcoming MDR have been widely explored in the last decades, leading to a generation of numerous small molecules targeting ABC and MRP transporters. Among the ABC family, ABCB1 plays key roles in the development of drug resistance and is the most well studied. In this work, we report the discovery of non-toxic [1,2,4]triazolo[1,5-a]pyrimidin-7-one (WS-10) from our structurally diverse in-house compound collection that selectively modulates ABCB1-mediated multidrug resistance. WS-10 enhanced the intracellular accumulation of paclitaxel in SW620/Ad300 cells, but did not affect the expression of ABCB1 Protein and ABCB1 localization. The cellular thermal shift assay (CETSA) showed that WS-10 was able to bind to ABCB1, which could be responsible for the reversal effect of WS-10 toward paclitaxel and doxorubicin in SW620/Ad300 cells. Docking simulations were performed to show the possible binding modes of WS-10 within ABCB1 transporter. To conclude, WS-10 could be used as a template for designing new ABCB1 modulators to overcome ABCB1-mediated multidrug resistance.

Structural and functional aspects of P-glycoprotein and its inhibitors

P-glycoprotein (P-gp) is a member of ATP-binding cassette (ABC) superfamily which extrudes chemotherapeutic agents out of the cell. Suppression of this efflux activity has been the subject of numerous attempts to develop P-gp inhibitors. The aim of this review is to present up-to-date information on the structural and functional aspects of P-gp and its known inhibitors. The data presented also provide some information on drug discovery approaches for candidate P-gp inhibitors. Nucleotide-binding domains (NBDs) and drug-binding domains (DBDs) have been extensively studied to gain more information about P-gp inhibition and it looks that the ATPase activity of this pump has been the most attractive target for designing inhibitors. Hydrophobic and π-π (aromatic) interactions between P-gp binding domains and inhibitors are dominant intermolecular forces that have been reported in many studies using different methods. Many synthetic and natural products have been found to possess inhibitory or modulatory effects on drug transporter proteins. Log P value is an important factor in studying these inhibitors and has a crucial role on absorption, distribution, metabolism, and excretion (ADME) properties of candidate P-gp inhibitors.

In vitro and in silico evaluation of P-glycoprotein inhibition through 99m Tc-methoxyisobutylisonitrile uptake

P-glycoprotein (P-gp) is a multidrug resistance (MDR) transporter with unknown structural details. This macromolecule is normally responsible for extruding xenobiotics from normal cells. Overexpression of P-gp in tumor cells is a major obstacle in cancer chemotherapy. In this study, human 3D model of P-gp was built by homology modeling based on mouse P-gp crystallographic structure and stabilized through 1 ns molecular dynamics (MD) simulation. Stabilized human P-gp structure was used for flexible docking of 80 drugs into the putative active site of P-gp. Accordingly, digoxin, itraconazole, risperidone, ketoconazole, prazosin, verapamil, cyclosporine A, and ranitidine were selected for further in vitro assay. Subsequently, cell-based P-gp inhibition assay was performed on Caco-2 cells while ^99m Tc-methoxyisobutylisonitrile (MIBI) was used as a P-gp efflux substrate for calculating IC50 values. Results of the ^99m Tc-MIBI uptake in drug-treated Caco-2 cells were in agreement with the previously reported activities. This study for the first time described the relation between molecular dynamics and flexible docking with cellular experiments using ^99m Tc-MIBI radiotracer for evaluation of potencies of P-gp inhibitors. Finally, results showed that our radiotracer-cell-based assay is an accurate and fast screening tool for detecting P-gp inhibitors and non-inhibitors in drug development process.

The reliability of molecular dynamics simulations of the multidrug transporter P-glycoprotein in a membrane environment

Despite decades of research, the mechanism of action of the ABC multidrug transporter P-glycoprotein (P-gp) remains elusive. Due to experimental limitations, many researchers have turned to molecular dynamics simulation studies in order to investigate different aspects of P-gp function. However, such studies are challenging and caution is required when interpreting the results. P-gp is highly flexible and the time scale on which it can be simulated is limited. There is also uncertainty regarding the accuracy of the various crystal structures available, let alone the structure of the protein in a physiologically relevant environment. In this study, three alternative structural models of mouse P-gp (3G5U, 4KSB, 4M1M), all resolved to 3.8 Å, were used to initiate sets of simulations of P-gp in a membrane environment in order to determine: a) the sensitivity of the results to differences in the starting configuration; and b) the extent to which converged results could be expected on the times scales commonly simulated for this system. The simulations suggest that the arrangement of the nucleotide binding domains (NBDs) observed in the crystal structures is not stable in a membrane environment. In all simulations, the NBDs rapidly associated (within 10 ns) and changes within the transmembrane helices were observed. The secondary structure within the transmembrane domain was best preserved in the 4M1M model under the simulation conditions used. However, the extent to which replicate simulations diverged on a 100 to 200 ns timescale meant that it was not possible to draw definitive conclusions as to which structure overall was most stable, or to obtain converged and reliable results for any of the properties examined. The work brings into question the reliability of conclusions made in regard to the nature of specific interactions inferred from previous simulation studies on this system involving similar sampling times. It also highlights the need to demonstrate the statistical significance of any results obtained in simulations of large flexible proteins, especially where the initial structure is uncertain.

Structure-function relationships in ABCG2: insights from molecular dynamics simulations and molecular docking studies

Efflux pumps of the ATP-binding cassette transporters superfamily (ABC transporters) are frequently involved in the multidrug-resistance (MDR) phenomenon in cancer cells. Herein, we describe a new atomistic model for the MDR-related ABCG2 efflux pump, also named breast cancer resistance protein (BCRP), based on the recently published crystallographic structure of the ABCG5/G8 heterodimer sterol transporter, a member of the ABCG family involved in cholesterol homeostasis. By means of molecular dynamics simulations and molecular docking, a far-reaching characterization of the ABCG2 homodimer was obtained. The role of important residues and motifs in the structural stability of the transporter was comprehensively studied and was found to be in good agreement with the available experimental data published in literature. Moreover, structural motifs potentially involved in signal transmission were identified, along with two symmetrical drug-binding sites that are herein described for the first time, in a rational attempt to better understand how drug binding and recognition occurs in ABCG2 homodimeric transporters.

Mechanism Underlying the Reversal of Drug Resistance in P-Glycoprotein-Expressing Leukemia Cells by Pinoresinol and the Study of a Derivative

P-glycoprotein (P-gp) is a membrane protein associated with multidrug resistance (MDR) due to its key role in mediating the traffic of chemotherapeutic drugs outside cancer cells, leading to a cellular response that hinders efforts toward successful therapy. With the aim of finding agents that circumvent the MDR phenotype mediated by P-gp, 15 compounds isolated from native and naturalized plants of Argentina were screened. Among these, the non-cytotoxic lignan (±) pinoresinol successfully restored sensitivity to doxorubicin from 7 μM in the P-gp overexpressed human myelogenous leukemia cells, Lucena 1. This resistance-reversing effect was confirmed by competitively increasing the intracellular doxorubicin accumulation and by significantly inhibiting the efflux of doxorubicin and, to a lesser extent, that of rhodamine 123. The activity obtained was similar to that observed with verapamil. No such results were observed in the sensitive parental K562 cell line. To gain deeper insight into the mode of action of pinoresinol, its effect on P-gp function and expression was examined. The docking simulations indicated that the lignan bound to P-gp at the apex of the V-shaped transmembrane cavity, involving transmembrane helices 4, 5, and 6, and partially overlapped the binding region of tariquidar, which was used as a positive control. These results would shed some light on the nature of its interaction with P-gp at molecular level and merit further mechanistic and kinetic studies. In addition, it showed a maximum 29% activation of ATP hydrolysis and antagonized verapamil-stimulated ATPase activity with an IC₅₀ of 20.9 μM. On the other hand, pinoresinol decreased the presence of P-gp in the cell surface. Derivatives of pinoresinol with improved activity were identified by docking studies. The most promising one, the non-cytotoxic 1-acetoxypinoresinol, caused a reversion of doxorubicin resistance from 0.11 μM and thus higher activity than the lead compound. It also caused a significant increase in doxorubicin accumulation. Results were similar to those observed with verapamil. The results obtained positioned these compounds as potential candidates for effective agents to overcome P-gp-mediated MDR, leading to better outcomes for leukemia chemotherapy.

Biophysical Approaches Facilitate Computational Drug Discovery for ATP-Binding Cassette Proteins

Although membrane proteins represent most therapeutically relevant drug targets, the availability of atomic resolution structures for this class of proteins has been limited. Structural characterization has been hampered by the biophysical nature of these polytopic transporters, receptors, and channels, and recent innovations to in vitro techniques aim to mitigate these challenges. One such class of membrane proteins, the ATP-binding cassette (ABC) superfamily, are broadly expressed throughout the human body, required for normal physiology and disease-causing when mutated, yet lacks sufficient structural representation in the Protein Data Bank. However, recent improvements to biophysical techniques (e.g., cryo-electron microscopy) have allowed for previously “hard-to-study” ABC proteins to be characterized at high resolution, providing insight into molecular mechanisms-of-action as well as revealing novel druggable sites for therapy design. These new advances provide ample opportunity for computational methods (e.g., virtual screening, molecular dynamics simulations, and structure-based drug design) to catalyze the discovery of novel small molecule therapeutics that can be easily translated from computer to bench and subsequently to the patient's bedside. In this review, we explore the utility of recent advances in biophysical methods coupled with well-established in silico techniques towards drug development for diseases caused by dysfunctional ABC proteins.

Computational modelling of efflux pumps and their inhibitors

Antimicrobial resistance is based on the multifarious strategies that bacteria adopt to face antibiotic therapies, making it a key public health concern of our era. Among these strategies, efflux pumps (EPs) contribute significantly to increase the levels and profiles of resistance by expelling a broad range of unrelated compounds - buying time for the organisms to develop specific resistance. In Gram-negative bacteria, many of these chromosomally encoded transporters form multicomponent 'pumps' that span both inner and outer membranes and are driven energetically by a primary or secondary transporter component.One of the strategies to reinvigorate the efficacy of antimicrobials is by joint administration with EP inhibitors (EPI), which either block the substrate binding and/or hinder any of the transport-dependent steps of the pump. In this review, we provide an overview of multidrug-resistance EPs, their inhibition strategies and the relevant findings from the various computational simulation studies reported to date with respect to deciphering the mechanism of action of inhibitors with the purpose of improving their rational design.

Reversal of Multidrug Resistance in Mouse Lymphoma Cells by Extracts and Flavonoids from Pistacia integerrima

Phytochemical investigation of Pistacia integerrima has highlighted isolation of two known compounds naringenin (1) and dihydrokaempferol (2). A crude extract and these isolated compounds were here evaluated for their effects on reversion of multidrug resistance (MDR) mediated by P-glycoprotein (P-gp). The multidrug resistance P-glycoprotein is a target for chemotherapeutic drugs from cancer cells. In the present study rhodamine- 123 exclusion screening test on human mdr1 gene transfected mouse gene transfected L5178 and L5178Y mouse T-cell lymphoma cells showed excellent MDR reversing effects in a dose dependent manner. In-silico molecular docking investigations demonstrated a common binding site for Rhodamine123, and compounds naringenin and dihydrokaempferol. Our results showed that the relative docking energies estimated by docking softwares were in satisfactory correlation with the experimental activities. Preliminary interaction profile of P-gp docked complexes were also analysed in order to understand the nature of binding modes of these compounds. Our computational investigation suggested that the compounds interactions with the hydrophobic pocket of P-gp are mainly related to the inhibitory activity. Moreover this study s a platform for the discovery of novel natural compounds from herbal origin, as inhibitor molecules against the P-glycoprotein for the treatment of cancer.

Homology modelling of human P-glycoprotein

P-glycoprotein (P-gp) is an ATP-binding cassette transporter that exports a huge range of compounds out of cells and is thus one of the key proteins in conferring multi-drug resistance in cancer. Understanding how it achieves such a broad specificity and the series of conformational changes that allow export to occur form major, on-going, research objectives around the world. Much of our knowledge to date has been derived from mutagenesis and assay data. However, in recent years, there has also been great progress in structural biology and although the structure of human P-gp has not yet been solved, there are now a handful of related structures on which homology models can be built to aid in the interpretation of the vast amount of experimental data that currently exists. Many models for P-gp have been built with this aim, but the situation is complicated by the apparent flexibility of the system and by the fact that although many potential templates exist, there is large variation in the conformational state in which they have been crystallized. In this review, we summarize how homology modelling has been used in the past, how models are typically selected and finally illustrate how MD simulations can be used as a means to give more confidence about models that have been generated via this approach.

Structural and dynamic perspectives on the promiscuous transport activity of P-glycoprotein

The multidrug transporter P-glycoprotein (P-gp) is expressed in the blood-brain barrier endothelium where it effluxes a range of drug substrates, preventing their accumulation within the brain. P-gp has been studied extensively for 40 years because of its crucial role in the absorption, distribution, metabolism and elimination of a range of pharmaceutical compounds. Despite this, many aspects of the structure-function mechanism of P-gp are unresolved. Here we review the emerging role of molecular dynamics simulation techniques in our understanding of the membrane-embedded conformation of P-gp. We discuss its conformational plasticity in the presence and absence of ATP, and recent efforts to characterize the drug binding sites and uptake pathways.

Mapping the Binding Site of the Inhibitor Tariquidar That Stabilizes the First Transmembrane Domain of P-glycoprotein

ABC (ATP-binding cassette) transporters are clinically important because drug pumps like P-glycoprotein (P-gp, ABCB1) confer multidrug resistance and mutant ABC proteins are responsible for many protein-folding diseases such as cystic fibrosis. Identification of the tariquidar-binding site has been the subject of intensive molecular modeling studies because it is the most potent inhibitor and corrector of P-gp. Tariquidar is a unique P-gp inhibitor because it locks the pump in a conformation that blocks drug efflux but activates ATPase activity. In silico docking studies have identified several potential tariquidar-binding sites. Here, we show through cross-linking studies that tariquidar most likely binds to sites within the transmembrane (TM) segments located in one wing or at the interface between the two wings (12 TM segments form 2 divergent wings). We then introduced arginine residues at all positions in the 12 TM segments (223 mutants) of P-gp. The rationale was that a charged residue in the drug-binding pocket would disrupt hydrophobic interaction with tariquidar and inhibit its ability to rescue processing mutants or stimulate ATPase activity. Arginines introduced at 30 positions significantly inhibited tariquidar rescue of a processing mutant and activation of ATPase activity. The results suggest that tariquidar binds to a site within the drug-binding pocket at the interface between the TM segments of both structural wings. Tariquidar differed from other drug substrates, however, as it stabilized the first TM domain. Stabilization of the first TM domain appears to be a key mechanism for high efficiency rescue of ABC processing mutants that cause disease.

Inhibition of P-glycoprotein Gene Expression and Function Enhances Triptolide-induced Hepatotoxicity in Mice

Triptolide (TP) is the major active principle of Tripterygium wilfordii Hook f. and very effective in treatment of autoimmune diseases. However, TP induced hepatotoxicity limited its clinical applications. Our previous study found that TP was a substrate of P-glycoprotein and its hepatobiliary clearance was markedly affected by P-gp modulation in sandwich-cultured rat hepatocytes. In this study, small interfering RNA (siRNA) and specific inhibitor tariquidar were used to investigate the impact of P-gp down regulation on TP-induced hepatotoxicity. The results showed that when the function of P-gp was inhibited by mdr1a-1 siRNA or tariquidar, the systemic and hepatic exposures of TP were significantly increased. The aggravated hepatotoxicity was evidenced with the remarkably lifted levels of serum biomarkers (ALT and AST) and pathological changes in liver. The other toxicological indicators (MDA, SOD and Bcl-2/Bax) were also significantly changed by P-gp inhibition. The data analysis showed that the increase of TP exposure in mice was quantitatively correlated to the enhanced hepatotoxicity, and the hepatic exposure was more relevant to the toxicity. P-gp mediated clearance played a significant role in TP detoxification. The risk of herb-drug interaction likely occurs when TP is concomitant with P-gp inhibitors or substrates in clinic.

Identification of Possible Binding Sites for Morphine and Nicardipine on the Multidrug Transporter P-Glycoprotein Using Umbrella Sampling Techniques

The multidrug transporter P-glycoprotein (P-gp) is central to the development of multidrug resistance in cancer. While residues essential for transport and binding have been identified, the location, composition, and specificity of potential drug binding sites are uncertain. Here molecular dynamics simulations are used to calculate the free energy profile for the binding of morphine and nicardipine to P-gp. We show that morphine and nicardipine primarily interact with key residues implicated in binding and transport from mutational studies, binding at different but overlapping sites within the transmembrane pore. Their permeation pathways were distinct but involved overlapping sets of residues. The results indicate that the binding location and permeation pathways of morphine and nicardipine are not well separated and cannot be considered as unique. This has important implications for our understanding of substrate uptake and transport by P-gp. Our results are independent of the choice of starting structure and consistent with a range of experimental studies.