Suppressor mutations in the transmembrane segments of P-glycoprotein promote maturation of processing mutants and disrupt a subset of drug-binding sites

In Silico and In Vitro Identification of P-Glycoprotein Inhibitors from a Library of 375 Phytochemicals

Cancer therapy with clinically established anticancer drugs is frequently hampered by the development of drug resistance of tumors and severe side effects in normal organs and tissues. The demand for powerful, but less toxic, drugs is high. Phytochemicals represent an important reservoir for drug development and frequently exert less toxicity than synthetic drugs. Bioinformatics can accelerate and simplify the highly complex, time-consuming, and expensive drug development process. Here, we analyzed 375 phytochemicals using virtual screenings, molecular docking, and in silico toxicity predictions. Based on these in silico studies, six candidate compounds were further investigated in vitro. Resazurin assays were performed to determine the growth-inhibitory effects towards wild-type CCRF-CEM leukemia cells and their multidrug-resistant, P-glycoprotein (P-gp)-overexpressing subline, CEM/ADR5000. Flow cytometry was used to measure the potential to measure P-gp-mediated doxorubicin transport. Bidwillon A, neobavaisoflavone, coptisine, and z-guggulsterone all showed growth-inhibitory effects and moderate P-gp inhibition, whereas miltirone and chamazulene strongly inhibited tumor cell growth and strongly increased intracellular doxorubicin uptake. Bidwillon A and miltirone were selected for molecular docking to wildtype and mutated P-gp forms in closed and open conformations. The P-gp homology models harbored clinically relevant mutations, i.e., six single missense mutations (F336Y, A718C, Q725A, F728A, M949C, Y953C), three double mutations (Y310A-F728A; F343C-V982C; Y953A-F978A), or one quadruple mutation (Y307C-F728A-Y953A-F978A). The mutants did not show major differences in binding energies compared to wildtypes. Closed P-gp forms generally showed higher binding affinities than open ones. Closed conformations might stabilize the binding, thereby leading to higher binding affinities, while open conformations may favor the release of compounds into the extracellular space. In conclusion, this study described the capability of selected phytochemicals to overcome multidrug resistance.

Extended-ensemble docking to probe dynamic variation of ligand binding sites during large-scale structural changes of proteins

Proteins can sample a broad landscape as they undergo conformational transition between different functional states. At the same time, as key players in almost all cellular processes, proteins are important drug targets. Considering the different conformational states of a protein is therefore central for a successful drug-design strategy. Here we introduce a novel docking protocol, termed extended-ensemble docking, pertaining to proteins that undergo large-scale (global) conformational changes during their function. In its application to multidrug ABC-transporter P-glycoprotein (Pgp), extensive non-equilibrium molecular dynamics simulations employing system-specific collective variables are first used to describe the transition cycle of the transporter. An extended set of conformations (extended ensemble) representing the full transition cycle between the inward- and the outward-facing states is then used to seed high-throughput docking calculations of known substrates, non-substrates, and modulators of the transporter. Large differences are predicted in the binding affinities to different conformations, with compounds showing stronger binding affinities to intermediate conformations compared to the starting crystal structure. Hierarchical clustering of the binding modes shows all ligands preferably bind to the large central cavity of the protein, formed at the apex of the transmembrane domain (TMD), whereas only small binding populations are observed in the previously described R and H sites present within the individual TMD leaflets. Based on the results, the central cavity is further divided into two major subsites, first preferably binding smaller substrates and high-affinity inhibitors, whereas the second one shows preference for larger substrates and low-affinity modulators. These central subsites along with the low-affinity interaction sites present within the individual TMD leaflets may respectively correspond to the proposed high- and low-affinity binding sites in Pgp. We propose further an optimization strategy for developing more potent inhibitors of Pgp, based on increasing its specificity to the extended ensemble of the protein, instead of using a single protein structure, as well as its selectivity for the high-affinity binding site. In contrast to earlier in silico studies using single static structures of Pgp, our results show better agreement with experimental studies, pointing to the importance of incorporating the global conformational flexibility of proteins in future drug-discovery endeavors.

Preclinical studies of the triazolo[1,5-a]pyrimidine derivative WS-716 as a highly potent, specific and orally active P-glycoprotein (P-gp) inhibitor

Multidrug resistance (MDR) is the main cause of clinical treatment failure and poor prognosis in cancer. Targeting P-glycoprotein (P-gp) has been regarded as an effective strategy to overcome MDR. In this work, we reported our preclinical studies of the triazolo[1,5-a]pyrimidine-based compound WS-716 as a highly potent, specific, and orally active P-gp inhibitor. Through direct binding to P-gp, WS-716 inhibited efflux function of P-gp and specifically reversed P-gp-mediated MDR to paclitaxel (PTX) in multiple resistant cell lines, without changing its expression or subcellular localization. WS-716 and PTX synergistically inhibited formation of colony and 3D spheroid, induced apoptosis and cell cycle arrest at G2/M phase in resistant SW620/Ad300 cells. In addition, WS-716 displayed minimal effect on the drug-metabolizing enzyme cytochrome P4503A4 (CYP3A4). Importantly, WS-716 increased sensitivity of both pre-clinically and clinically derived MDR tumors to PTX in vivo with the T/C value of 29.7% in patient-derived xenograft (PDX) models. Relative to PTX treatment alone, combination of WS-716 and PTX caused no obvious adverse reactions. Taken together, our preclinical studies revealed therapeutic promise of WS-716 against MDR cancer, the promising data warrant its further development for cancer therapy.

Identification and characterisation of putative drug binding sites in human ATP-binding cassette B5 (ABCB5) transporter

The human ATP-binding cassette B5 (ABCB5) transporter, a member of the ABC transporter superfamily, is linked to chemoresistance in tumour cells by drug effluxion. However, little is known about its structure and drug-binding sites. In this study, we generated an atomistic model of the full-length human ABCB5 transporter with the highest quality using the X-ray crystal structure of mouse ABCB1 (Pgp1), a close homologue of ABCB5 and a well-studied member of the ABC family. Molecular dynamics simulations were used to validate the atomistic model of ABCB5 and characterise its structural properties in model cell membranes. Molecular docking simulations of known ABCB5 substrates such as taxanes, anthracyclines, camptothecin and etoposide were then used to identify at least three putative binding sites for chemotherapeutic drugs transported by ABCB5. The location of these three binding sites is predicted to overlap with the corresponding binding sites in Pgp1. These findings will serve as the basis for future in vitro studies to validate the nature of the identified substrate-binding sites in the full-length ABCB5 transporter.

Effect of ABC transporter expression and mutational status on survival rates of cancer patients

ATP-binding cassette (ABC) transporters mediate multidrug resistance in cancer. In contrast to DNA single nucleotide polymorphisms in normal tissues, the role of mutations in tumors is unknown. Furthermore, the significance of their expression for prediction of chemoresistance and survival prognosis is still under debate. We investigated 18 tumors by RNA-sequencing. The mutation rate varied from 27,507 to 300885. In ABCB1, three hotspots with novel mutations were in transmembrane domains 3, 8, and 9. We also mined the cBioPortal database with 11,814 patients from 23 different tumor entities. We performed Kaplan-Meier survival analyses to investigate the effect of ABC transporter expression on survival rates of cancer patients. Novel mutations were also found in ABCA2, ABCA3, ABCB2, ABCB5, ABCC1-6, and ABCG2. Mining the cBioPortal database with 11,814 patients from 23 different tumor entities validated our results. Missense and in-frame mutations led to altered binding of anticancer drugs in molecular docking approaches. The ABCB1 nonsense mutation Q856* led to a truncated P-glycoprotein, which may sensitize tumors to anticancer drugs. The search for ABC transporter nonsense mutations represents a novel approach for precision medicine.. Low ABCB1 mRNA expression correlated with significantly longer survival in ovarian or kidney cancer and thymoma. In cancers of breast, kidney or lung, ABC transporter expression correlated with different tumor stages and human populations as further parameters to refine strategies for more individualized chemotherapy.

ABCB transporters in a leaf beetle respond to sequestered plant toxins

Phytophagous insects can tolerate and detoxify toxic compounds present in their host plants and have evolved intricate adaptations to this end. Some insects even sequester the toxins for their defence. This necessitates specific mechanisms, especially carrier proteins that regulate uptake and transport to specific storage sites or protect sensitive tissues from noxious compounds. We identified three ATP-binding cassette subfamily B (ABCB) transporters from the transcriptome of the cardenolide-sequestering leaf beetle Chrysochus auratus and analysed their functional role in the sequestration process. These were heterologously expressed and tested for their ability to interact with various potential substrates: verapamil (standard ABCB substrate), the cardenolides digoxin (commonly used), cymarin (present in the species's host plant) and calotropin (present in the ancestral host plants). Verapamil stimulated all three ABCBs and each was activated by at least one cardenolide, however, they differed as to which they were activated by. While the expression of the most versatile transporter fits with a protective role in the blood-brain barrier, the one specific for cymarin shows an extreme abundance in the elytra, coinciding with the location of the defensive glands. Our data thus suggest a key role of ABCBs in the transport network needed for cardenolide sequestration.

Medicinal chemistry strategies to discover P-glycoprotein inhibitors: An update

The presence of multidrug resistance (MDR) in malignant tumors is one of the primary causes of treatment failure in cancer chemotherapy. The overexpression of the ATP binding cassette (ABC) transporter, P-glycoprotein (P-gp), which significantly increases the efflux of certain anticancer drugs from tumor cells, produces MDR. Therefore, inhibition of P-gp may represent a viable therapeutic strategy to overcome cancer MDR. Over the past 4 decades, many compounds with P-gp inhibitory efficacy (referred to as first- and second-generation P-gp inhibitors) have been identified or synthesized. However, these compounds were not successful in clinical trials due to a lack of efficacy and/or untoward toxicity. Subsequently, third- and fourth-generation P-gp inhibitors were developed but dedicated clinical trials did not indicate a significant therapeutic effect. In recent years, an extraordinary array of highly potent, selective, and low-toxicity P-gp inhibitors have been reported. Herein, we provide a comprehensive review of the synthetic and natural products that have specific inhibitory activity on P-gp drug efflux as well as promising chemosensitizing efficacy in MDR cancer cells. The present review focuses primarily on the structural features, design strategies, and structure-activity relationships (SAR) of these compounds.

Evidence for the critical role of transmembrane helices 1 and 7 in substrate transport by human P-glycoprotein (ABCB1)

P-glycoprotein (P-gp) is an ABC transporter that exports many amphipathic or hydrophobic compounds, including chemically and functionally dissimilar anticancer drugs, from cells. To understand the role of transmembrane helices (TMH) 1 and 7 in drug-binding and transport, we selected six residues from both TMH1 (V53, I59, I60, L65, M68 and F72) and TMH7 (V713, I719, I720, Q725, F728 and F732); and substituted them with alanine by gene synthesis to generate a variant termed "TMH1,7 mutant P-gp". The expression and function of TMH1,7 mutant P-gp with twelve mutations was characterized using the BacMam baculovirus-HeLa cell expression system. The expression and conformation of TMH1,7 mutant P-gp was not altered by the introduction of the twelve mutations, as confirmed by using the human P-gp-specific antibodies UIC2, MRK16 and 4E3. We tested 25 fluorescently-labeled substrates and found that only three substrates, NBD-cyclosporine A, Rhod-2-AM and X-Rhod-1-AM were transported by the TMH1,7 mutant. The basal ATPase activity of TMH1,7 mutant P-gp was lower (40-50%) compared to wild-type (WT) P-gp, despite similar level of expression. Although most of the substrates modulate ATPase activity of P-gp, the activity of TMH1,7 mutant transporter was not significantly modulated by any of the tested substrates. Docking of selected substrates in homology models showed comparable docking scores for the TMH1,7 mutant and WT P-gp, although the binding conformations were different. Both the ATPase assay and in silico docking analyses suggest that the interactions with residues in the drug-binding pocket are altered as a consequence of the mutations. We demonstrate that it is possible to generate a variant of P-gp with a loss of broad substrate specificity and propose that TMH1 and TMH7 play a critical role in the drug efflux function of this multidrug transporter.

Targeting P-glycoprotein: Investigation of piperine analogs for overcoming drug resistance in cancer

P-glycoprotein (P-gp) is a drug transporter that effluxes chemotherapeutic drugs and is implicated in the development of resistance of cancer cells to chemotherapeutic drugs. To date, no drug has been approved to inhibit P-gp and restore chemotherapy efficacy. Moreover, majority of the reported inhibitors have high molecular weight and complex structures, making it difficult to understand the basic structural requirement for P-gp inhibition. In this study, two structurally simple, low molecular weight piperine analogs Pip1 and Pip2 were designed and found to better interact with P-gp than piperine in silico. A one step, acid-amine coupling reaction between piperic acid and 6,7-dimethoxytetrahydroisoquinoline or 2-(3,4-dimethoxyphenyl)ethylamine afforded Pip1 and Pip2, respectively. In vitro testing in drug resistant P-gp overexpressing KB (cervical) and SW480 (colon) cancer cells showed that both analogs, when co-administered with vincristine, colchicine or paclitaxel were able to reverse the resistance. Moreover, accumulation of P-gp substrate (rhodamine 123) in the resistant cells, a result of alteration of the P-gp efflux, was also observed. These investigations suggest that the natural product analog - Pip1 ((2E,4E)-5-(benzo[d][1,3]dioxol-5-yl)-1-(6,7-dimethoxy-3,4-dihydroisoquinolin-2(1 H)-yl)penta-2,4-dien-1-one) - is superior to piperine and could inhibit P-gp function. Further studies are required to explore the full potential of Pip1 in treating drug resistant cancer.

Biochemical studies on the structure-function relationship of major drug transporters in the ATP-binding cassette family and solute carrier family

Human drug transporters often play key roles in determining drug accumulation within cells. Their activities are often directly related to therapeutic efficacy, drug toxicity as well as drug-drug interactions. However, the progress for interpretation of their crystal structures is relatively slow. Hence, conventional biochemical studies together with computer modeling became useful manners to reveal essential structures of these membrane proteins. Over the years, quite a few structure-function relationship information had been obtained for members of the two major transporter families: the ATP-binding cassette family and the solute carrier family. Critical structural features of drug transporters include transmembrane domains, post-translational modification sites and domains for cell surface assembly and protein-protein interactions. Alterations at these important sites may affect protein stability, trafficking to the plasma membrane and/or ability of transporters to interact with substrates.

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.

Multiple Drug Transport Pathways through Human P-Glycoprotein

P-Glycoprotein (P-gp) is a plasma membrane efflux pump that is commonly associated with therapy resistances in cancers and infectious diseases. P-gp can lower the intracellular concentrations of many drugs to subtherapeutic levels by translocating them out of the cell. Because of the broad range of substrates transported by P-gp, overexpression of P-gp causes multidrug resistance. We reported previously on dynamic transitions of P-gp as it moved through conformations based on crystal structures of homologous ABCB1 proteins using in silico targeted molecular dynamics techniques. We expanded these studies here by docking transport substrates to drug binding sites of P-gp in conformations open to the cytoplasm, followed by cycling the pump through conformations that opened to the extracellular space. We observed reproducible transport of two substrates, daunorubicin and verapamil, by an average of 11-12 Å through the plane of the membrane as P-gp progressed through a catalytic cycle. Methylpyrophosphate, a ligand that should not be transported by P-gp, did not show this movement through P-gp. Drug binding to either of two subsites on P-gp appeared to determine the initial pathway used for drug movement through the membrane. The specific side-chain interactions with drugs within each pathway seemed to be, at least in part, stochastic. The docking and transport properties of a P-gp inhibitor, tariquidar, were also studied. A mechanism of inhibition by tariquidar that involves stabilization of an outward open conformation with tariquidar bound in intracellular loops or at the drug binding domain of P-gp is presented.

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.

The Transmission Interfaces Contribute Asymmetrically to the Assembly and Activity of Human P-glycoprotein

P-glycoprotein (P-gp; ABCB1) is an ABC drug pump that protects us from toxic compounds. It is clinically important because it confers multidrug resistance. The homologous halves of P-gp each contain a transmembrane (TM) domain (TMD) with 6 TM segments followed by a nucleotide-binding domain (NBD). The drug- and ATP-binding sites reside at the interface between the TMDs and NBDs, respectively. Each NBD is connected to the TMDs by a transmission interface involving a pair of intracellular loops (ICLs) that form ball-and-socket joints. P-gp is different from CFTR (ABCC7) in that deleting NBD2 causes misprocessing of only P-gp. Therefore, NBD2 might be critical for stabilizing ICLs 2 and 3 that form a tetrahelix bundle at the NBD2 interface. Here we report that the NBD1 and NBD2 transmission interfaces in P-gp are asymmetric. Point mutations to 25 of 60 ICL2/ICL3 residues at the NBD2 transmission interface severely reduced P-gp assembly while changes to the equivalent residues in ICL1/ICL4 at the NBD1 interface had little effect. The hydrophobic nature at the transmission interfaces was also different. Mutation of Phe-1086 or Tyr-1087 to arginine at the NBD2 socket blocked activity or assembly while the equivalent mutations at the NBD1 socket had only modest effects. The results suggest that the NBD transmission interfaces are asymmetric. In contrast to the ICL2/3-NBD2 interface, the ICL1/4-NBD1 transmission interface is more hydrophilic and insensitive to mutations. Therefore the ICL2/3-NBD2 transmission interface forms a precise hydrophobic connection that acts as a linchpin for assembly and trafficking of P-gp.

Convallatoxin: a new P-glycoprotein substrate

Digitalis-like compounds (DLCs), such as digoxin and digitoxin that are derived from digitalis species, are currently used to treat heart failure and atrial fibrillation, but have a narrow therapeutic index. Drug-drug interactions at the transporter level are frequent causes of DLCs toxicity. P-glycoprotein (P-gp, ABCB1) is the primary transporter of digoxin and its inhibitors influence pharmacokinetics and disposition of digoxin in the human body; however, the involvement of P-gp in the disposition of other DLCs is currently unknown. In present study, the transport of fourteen DLCs by human P-gp was studied using membrane vesicles originating from human embryonic kidney (HEK293) cells overexpressing P-gp. DLCs were quantified by liquid chromatography-mass spectrometry (LC-MS). The Lily of the Valley toxin, convallatoxin, was identified as a P-gp substrate (Km: 1.1±0.2 mM) in the vesicular assay. Transport of convallatoxin by P-gp was confirmed in rat in vivo, in which co-administration with the P-gp inhibitor elacridar, resulted in increased concentrations in brain and kidney cortex. To address the interaction of convallatoxin with P-gp on a molecular level, the effect of nine alanine mutations was compared with the substrate N-methyl quinidine (NMQ). Phe343 appeared to be more important for transport of NMQ than convallatoxin, while Val982 was particularly relevant for convallatoxin transport. We identified convallatoxin as a new P-gp substrate and recognized Val982 as an important amino acid involved in its transport. These results contribute to a better understanding of the interaction of DLCs with P-gp.

The cystic fibrosis V232D mutation inhibits CFTR maturation by disrupting a hydrophobic pocket rather than formation of aberrant interhelical hydrogen bonds

Processing mutations that inhibit folding and trafficking of CFTR are the main cause of cystic fibrosis. Repair of CFTR mutants requires an understanding of the mechanisms of misfolding caused by processing mutations. Previous studies on helix-loop-helix fragments of the V232D processing mutation suggested that its mechanism was to lock transmembrane (TM) segments 3 and 4 together by a non-native hydrogen bond (Asp232(TM4)/Gln207(TM3)). Here, we performed mutational analysis to test for Asp232/Gln207 interactions in full-length CFTR. The rationale was that a V232N mutation should mimic V232D and a V232D/Q207A mutant should mature if the processing defect was caused by hydrogen bonds. We report that only Val232 mutations to charged amino acids severely blocked CFTR maturation. The V232N mutation did not mimic V232D as V232N showed 40% maturation compared to 2% for V232D. Mutation of Val232 to large nonpolar residues (Leu, Phe) had little effect. The Q207L mutation did not rescue V232D because Q207L showed about 50% maturation in the presence of corrector VX-809 while V232D/Q207A could no longer be rescued. These results suggest that V232D inhibits maturation by disrupting a hydrophobic pocket between TM segments rather than forming a non-native hydrogen bond. Disulfide cross-linking analysis of cysteines W356C(TM6) and W1145C(TM12) suggest that the V232D mutation inhibits maturation by trapping CFTR as a partially folded intermediate. Since correctors can efficiently rescue V232D CFTR, the results suggest that hydrophilic processing mutations facing a hydrophobic pocket are good candidates for rescue with pharmacological chaperones.

Locking intracellular helices 2 and 3 together inactivates human P-glycoprotein

The P-glycoprotein (P-gp) drug pump (ABCB1) has two transmembrane domains and two nucleotide-binding domains (NBDs). Coupling of the drug-binding sites in the transmembrane domains to the NBDs occurs through interaction of the intracellular helices (IHs) with residues in the NBDs (IH1/IH4/NBD1 and IH2/IH3/NBD2). We showed previously that cross-linking of cysteines in IH3 and IH1 with a short cross-linker mimicked drug binding as it activated P-gp ATPase activity. Here we show that residue A259C(IH2) could be directly cross-linked to W803C(IH3). Cross-linking was inhibited by the presence of ATP and adenosine 5'-(β,γ-imino)triphosphate but not by ADP. Cross-linking of mutant A259C/W803C inhibited its verapamil-stimulated ATPase activity mutant, but activity was restored after addition of dithiothreitol. Because these residues are close to the ball-and-socket joint A266C(IH2)/Phe(1086)(NBD2), we mutated the adjacent Tyr(1087)(NBD2) close to IH3. Mutants Y1087A and Y1087L, but not Y1087F, were misprocessed, and all inhibited ATPase activity. Mutation of hydrophobic residues (F793A, L797A, L814A, and L818A) flanking IH3 also inhibited maturation. The results suggest that these residues, together with Trp(803) and Phe(804), form a large hydrophobic pocket. The results show that there is an important hydrophobic network at the IH2/IH3/NBD2 transmission interface that is critical for folding and activity of P-gp.

Refined structures of mouse P-glycoprotein

The recently determined C. elegans P-glycoprotein (Pgp) structure revealed significant deviations compared to the original mouse Pgp structure, which suggested possible misinterpretations in the latter model. To address this concern, we generated an experimental electron density map from single-wavelength anomalous dispersion phasing of an original mouse Pgp dataset to 3.8 Å resolution. The map exhibited significantly more detail compared to the original MAD map and revealed several regions of the structure that required de novo model building. The improved drug-free structure was refined to 3.8 Å resolution with a 9.4 and 8.1% decrease in R_work and R_free, respectively, (R_work = 21.2%, R_free = 26.6%) and a significant improvement in protein geometry. The improved mouse Pgp model contains ∼95% of residues in the favorable Ramachandran region compared to only 57% for the original model. The registry of six transmembrane helices was corrected, revealing amino acid residues involved in drug binding that were previously unrecognized. Registry shifts (rotations and translations) for three transmembrane (TM)4 and TM5 and the addition of three N-terminal residues were necessary, and were validated with new mercury labeling and anomalous Fourier density. The corrected position of TM4, which forms the frame of a portal for drug entry, had backbone atoms shifted >6 Å from their original positions. The drug translocation pathway of mouse Pgp is 96% identical to human Pgp and is enriched in aromatic residues that likely play a collective role in allowing a high degree of polyspecific substrate recognition.

Drug rescue distinguishes between different structural models of human P-glycoprotein

There is no high-resolution crystal structure of the human P-glycoprotein (P-gp) drug pump. Homology models of human P-gp based on the crystal structures of mouse or Caenorhabditis elegans P-gps show large differences in the orientation of transmembrane segment 5 (TM5). TM5 is one of the most important transmembrane segments involved in drug–substrate interactions. Drug rescue of P-gp processing mutants containing an arginine at each position in TM5 was used to identify positions facing the lipid or internal aqueous chamber. Only the model based on the C. elegans P-gp structure was compatible with the drug rescue results.

Insight into pleiotropic drug resistance ATP-binding cassette pump drug transport through mutagenesis of Cdr1p transmembrane domains

The fungal ATP-binding cassette (ABC) transporter Cdr1 protein (Cdr1p), responsible for clinically significant drug resistance, is composed of two transmembrane domains (TMDs) and two nucleotide binding domains (NBDs). We have probed the nature of the drug binding pocket by performing systematic mutagenesis of the primary sequences of the 12 transmembrane segments (TMSs) found in the TMDs. All mutated proteins were expressed equally well and localized properly at the plasma membrane in the heterologous host Saccharomyces cerevisiae, but some variants differed significantly in efflux activity, substrate specificity, and coupled ATPase activity. Replacement of the majority of the amino acid residues with alanine or glycine yielded neutral mutations, but about 42% of the variants lost resistance to drug efflux substrates completely or selectively. A predicted three-dimensional homology model shows that all the TMSs, apart from TMS4 and TMS10, interact directly with the drug-binding cavity in both the open and closed Cdr1p conformations. However, TMS4 and TMS10 mutations can also induce total or selective drug susceptibility. Functional data and homology modeling assisted identification of critical amino acids within a drug-binding cavity that, upon mutation, abolished resistance to all drugs tested singly or in combinations. The open and closed Cdr1p models enabled the identification of amino acid residues that bordered a drug-binding cavity dominated by hydrophobic residues. The disposition of TMD residues with differential effects on drug binding and transport are consistent with a large polyspecific drug binding pocket in this yeast multidrug transporter.

Human P-glycoprotein contains a greasy ball-and-socket joint at the second transmission interface

The P-glycoprotein drug pump protects us from toxins. Drug-binding sites in the transmembrane (TM) domains (TMDs) are connected to the nucleotide-binding domains (NBDs) by intracellular helices (IHs). TMD-NBD cross-talk is a key step in the transport mechanism because drug binding stimulates ATP hydrolysis followed by drug efflux. Here, we tested whether the IHs are critical for maturation and TMD-NBD coupling by characterizing the effects of mutations to the IH1 and IH2 interfaces. Although IH1 mutations had little effect, most mutations at the IH2-NBD2 interface inhibited maturation or activity. For example, the F1086A mutation at the IH2-NBD2 interface abolished drug-stimulated ATPase activity. The mutant F1086A, however, retained the ability to bind ATP and drug substrates. The mutant was defective in mediating ATP-dependent conformational changes in the TMDs because binding of ATP no longer promoted cross-linking between cysteines located at the extracellular ends of TM segments 6 and 12. Replacement of Phe-1086 (in NBD2) with hydrophobic but not charged residues yielded active mutants. The activity of the F1086A mutant could be restored when the nearby residue Ala-266 (in IH2) was replaced with aromatic residues. These results suggest that Ala-266/Phe-1086 lies in a hydrophobic IH2-NBD2 "ball-and-socket" joint.

A salt bridge in intracellular loop 2 is essential for folding of human p-glycoprotein

There is no high-resolution structure of the human P-glycoprotein (P-gp, ABCB1) drug pump. Homology models based on the crystal structures of mouse and Caenorhabditis elegans P-gps show extensive contacts between intracellular loop 2 (ICL2, in the first transmembrane domain) and the second nucleotide-binding domain. Human P-gp modeled on these P-gp structures yields different ICL2 structures. Only the model based on the C. elegans P-gp structure predicts the presence of a salt bridge. We show that the Glu256–Arg276 salt bridge was critical for P-gp folding.

Interaction of digitalis-like compounds with p-glycoprotein

Digitalis-like compounds (DLCs), or cardiac glycosides, are produced and sequestered by certain plants and animals as a protective mechanism against herbivores or predators. Currently, the DLCs digoxin and digitoxin are used in the treatment of cardiac congestion and some types of cardiac arrhythmia, despite a very narrow therapeutic index. P-glycoprotein (P-gp; ABCB1) is the only known ATP-dependent efflux transporter that handles digoxin as a substrate. Ten alanine mutants of human P-gp drug-binding amino acids-Leu(65), Ile(306), Phe(336), Ile(340), Phe(343), Phe(728), Phe(942), Thr(945), Leu(975), and Val(982)-were generated and expressed in HEK293 cells with a mammalian baculovirus system. The uptake of [(3)H]-N-methyl-quinidine (NMQ), the P-gp substrate in vesicular transport assays, was determined. The mutations I306A, F343A, F728A, T945A, and L975A abolished NMQ transport activity of P-gp. For the other mutants, the apparent affinities for six DLCs (cymarin, digitoxin, digoxin, peruvoside, proscillaridin A, and strophanthidol) were determined. The affinities of digoxin, proscillaridin A, peruvoside, and cymarin for mutants F336A and I340A were decreased two- to fourfold compared with wild type, whereas that of digitoxin and strophanthidol did not change. In addition, the presence of a hydroxyl group at position 12β seems to reduce the apparent affinity when the side chain of Phe(336) and Phe(942) is absent. Our results showed that a δ-lactone ring and a sugar moiety at 3β of the steroid body are favorable for DLC binding to P-gp. Moreover, DLC inhibition is increased by hydroxyl groups at positions 5β and 19, whereas inhibition is decreased by those at positions 1β, 11α, 12β, and 16β. The understanding of the P-gp-DLC interaction improves our insight into DLCs toxicity and might enhance the replacement of digoxin with other DLCs that have less adverse drug effects.

Crystal structure of the multidrug transporter P-glycoprotein from Caenorhabditis elegans

P-glycoprotein (P-gp) is an ATP-binding cassette transporter that confers multidrug resistance in cancer cells. It also affects the absorption, distribution and clearance of cancer-unrelated drugs and xenobiotics. For these reasons, the structure and function of P-gp have been studied extensively for decades. Here we present biochemical characterization of P-gp from Caenorhabditis elegans and its crystal structure at a resolution of 3.4 ångströms. We find that the apparent affinities of P-gp for anticancer drugs actinomycin D and paclitaxel are approximately 4,000 and 100 times higher, respectively, in the membrane bilayer than in detergent. This affinity enhancement highlights the importance of membrane partitioning when a drug accesses the transporter in the membrane. Furthermore, the transporter in the crystal structure opens its drug pathway at the level of the membrane's inner leaflet. In the helices flanking the opening to the membrane, we observe extended loops that may mediate drug binding, function as hinges to gate the pathway or both. We also find that the interface between the transmembrane and nucleotide-binding domains, which couples ATP hydrolysis to transport, contains a ball-and-socket joint and salt bridges similar to the ATP-binding cassette importers, suggesting that ATP-binding cassette exporters and importers may use similar mechanisms to achieve alternating access for transport. Finally, a model of human P-gp derived from the structure of C. elegans P-gp not only is compatible with decades of biochemical analysis, but also helps to explain perplexing functional data regarding the Phe335Ala mutant. These results increase our understanding of the structure and function of this important molecule.

The ATPase activity of the P-glycoprotein drug pump is highly activated when the N-terminal and central regions of the nucleotide-binding domains are linked closely together

The P-glycoprotein (P-gp, ABCB1) drug pump protects us from toxic compounds and confers multidrug resistance. Each of the homologous halves of P-gp is composed of a transmembrane domain (TMD) with 6 TM segments followed by a nucleotide-binding domain (NBD). The predicted drug- and ATP-binding sites reside at the interface between the TMDs and NBDs, respectively. Crystal structures and EM projection images suggest that the two halves of P-gp are separated by a central cavity that closes upon binding of nucleotide. Binding of drug substrates may induce further structural rearrangements because they stimulate ATPase activity. Here, we used disulfide cross-linking with short (8 Å) or long (22 Å) cross-linkers to identify domain-domain interactions that activate ATPase activity. It was found that cross-linking of cysteines that lie close to the LSGGQ (P517C) and Walker A (I1050C) sites of NBD1 and NBD2, respectively, as well as the cytoplasmic extensions of TM segments 3 (D177C or L175C) and 9 (N820C) with a short cross-linker activated ATPase activity over 10-fold. A pyrylium compound that inhibits ATPase activity blocked cross-linking at these sites. Cross-linking between the NBDs was not inhibited by tariquidar, a drug transport inhibitor that stimulates P-gp ATPase activity but is not transported. Cross-linking between extracellular cysteines (T333C/L975C) predicted to lock P-gp into a conformation that prevents close NBD association inhibited ATPase activity. The results suggest that trapping P-gp in a conformation in which the NBDs are closely associated likely mimics the structural rearrangements caused by binding of drug substrates that stimulate ATPase activity.

Exhaustive sampling of docking poses reveals binding hypotheses for propafenone type inhibitors of P-glycoprotein

Overexpression of the xenotoxin transporter P-glycoprotein (P-gp) represents one major reason for the development of multidrug resistance (MDR), leading to the failure of antibiotic and cancer therapies. Inhibitors of P-gp have thus been advocated as promising candidates for overcoming the problem of MDR. However, due to lack of a high-resolution structure the concrete mode of interaction of both substrates and inhibitors is still not known. Therefore, structure-based design studies have to rely on protein homology models. In order to identify binding hypotheses for propafenone-type P-gp inhibitors, five different propafenone derivatives with known structure-activity relationship (SAR) pattern were docked into homology models of the apo and the nucleotide-bound conformation of the transporter. To circumvent the uncertainty of scoring functions, we exhaustively sampled the pose space and analyzed the poses by combining information retrieved from SAR studies with common scaffold clustering. The results suggest propafenone binding at the transmembrane helices 5, 6, 7 and 8 in both models, with the amino acid residue Y307 playing a crucial role. The identified binding site in the non-energized state is overlapping with, but not identical to, known binding areas of cyclic P-gp inhibitors and verapamil. These findings support the idea of several small binding sites forming one large binding cavity. Furthermore, the binding hypotheses for both catalytic states were analyzed and showed only small differences in their protein-ligand interaction fingerprints, which indicates only small movements of the ligand during the catalytic cycle.

Exhaustive sampling of docking poses reveals binding hypotheses for propafenone type inhibitors of P-glycoprotein

Overexpression of the xenotoxin transporter P-glycoprotein (P-gp) represents one major reason for the development of multidrug resistance (MDR), leading to the failure of antibiotic and cancer therapies. Inhibitors of P-gp have thus been advocated as promising candidates for overcoming the problem of MDR. However, due to lack of a high-resolution structure the concrete mode of interaction of both substrates and inhibitors is still not known. Therefore, structure-based design studies have to rely on protein homology models. In order to identify binding hypotheses for propafenone-type P-gp inhibitors, five different propafenone derivatives with known structure-activity relationship (SAR) pattern were docked into homology models of the apo and the nucleotide-bound conformation of the transporter. To circumvent the uncertainty of scoring functions, we exhaustively sampled the pose space and analyzed the poses by combining information retrieved from SAR studies with common scaffold clustering. The results suggest propafenone binding at the transmembrane helices 5, 6, 7 and 8 in both models, with the amino acid residue Y307 playing a crucial role. The identified binding site in the non-energized state is overlapping with, but not identical to, known binding areas of cyclic P-gp inhibitors and verapamil. These findings support the idea of several small binding sites forming one large binding cavity. Furthermore, the binding hypotheses for both catalytic states were analyzed and showed only small differences in their protein-ligand interaction fingerprints, which indicates only small movements of the ligand during the catalytic cycle.

Intragenic suppressing mutations correct the folding and intracellular traffic of misfolded mutants of Yor1p, a eukaryotic drug transporter

ATP-binding cassette (ABC) transporters play pivotal physiological roles in substrate transport across membranes, and defective assembly of these proteins can cause severe disease associated with improper drug or ion flux. The yeast protein Yor1p is a useful model to study the biogenesis of ABC transporters; deletion of a phenylalanine residue in the first nucleotide-binding domain (NBD1) causes misassembly and retention in the endoplasmic reticulum (ER) of the resulting protein Yor1p-ΔF670, similar to the predominant disease-causing allele in humans, CFTR-ΔF508. Here we describe two novel Yor1p mutants, G278R and I1084P, which fail to assemble and traffic similar to Yor1p-ΔF670. These mutations are located in the two intracellular loops (ICLs) that interface directly with NBD1, and thus disrupt a functionally important structural module. We isolated 2 second-site mutations, F270S and R1168M, which partially correct the folding injuries associated with the G278R, I1084P, and ΔF670 mutants and reinstate their trafficking. The position of both corrective mutations at the cytoplasmic face of a transmembrane helix suggests that they restore biogenesis by influencing the behavior of the transmembrane domains rather than by direct restoration of the ICL1-ICL4-NBD1 structural module. Given the conserved topology of many ABC transporters, our findings provide new understanding of functionally important inter-domain interactions and suggest new potential avenues for correcting folding defects caused by abrogation of those domain interfaces.

A promiscuous conformational switch in the secondary multidrug transporter MdfA

Multidrug (Mdr) transporters are membrane proteins that actively export structurally dissimilar drugs from the cell, thereby rendering the cell resistant to toxic compounds. Similar to substrate-specific transporters, Mdr transporters also undergo substrate-induced conformational changes. However, the mechanism by which a variety of dissimilar substrates are able to induce similar transport-compatible conformational responses in a single transporter remains unclear. To address this major aspect of Mdr transport, we studied the conformational behavior of the Escherichia coli Mdr transporter MdfA. Our results show that indeed, different substrates induce similar conformational changes in the transporter. Intriguingly, in addition, we observed that compounds other than substrates are able to confer similar conformational changes when covalently attached at the putative Mdr recognition pocket of MdfA. Taken together, the results suggest that the Mdr-binding pocket of MdfA is conformationally sensitive. We speculate that the same conformational switch that usually drives active transport is triggered promiscuously by merely occupying the Mdr-binding site.

Identification of residues in the drug translocation pathway of the human multidrug resistance P-glycoprotein by arginine mutagenesis

P-glycoprotein (P-gp, ATP-binding cassette B1) is a drug pump that extracts toxic drug substrates from the plasma membrane and catalyzes their ATP-dependent efflux. To map the residues in the drug translocation pathway, we performed arginine-scanning mutagenesis on all transmembrane (TM) segments (total = 237 residues) of a P-gp processing mutant (G251V) defective in folding (15% maturation efficiency) (glycosylation state used to monitor folding). The rationale was that arginines introduced into the drug-binding sites would mimic drug rescue and enhance maturation of wild-type or processing mutants of P-gp. It was found that 38 of the 89 mutants that matured had enhanced maturation. Enhancer mutations were found in 11 of the 12 TM segments with the largest number found in TMs 6 and 12 (seven in each), TMs that are critical for P-gp-drug substrate interactions. Modeling of the TM segments showed that the enhancer arginines were found on the hydrophilic face, whereas inhibitory arginines were located on a hydrophobic face that may be in contact with the lipid bilayer. It was found that many of the enhancer arginines caused large alterations in P-gp-drug interactions in ATPase assays. For example, mutants A302R (TM5), L339R (TM6), G872R (TM10), F942R (TM11), Q946R (TM11), V982R (TM12), and S993R (TM12) reduced the apparent affinity for verapamil by approximately 10-fold, whereas the F336R (TM6) and M986R (TM12) mutations caused at least a 10-fold increase in apparent affinity for rhodamine B. The results suggest that P-gp contains a large aqueous-filled drug translocation pathway with multiple drug-binding sites that can accommodate the bulky arginine side chains to promote folding of the protein.

Molecular models of human P-glycoprotein in two different catalytic states

Background

P-glycoprotein belongs to the family of ATP-binding cassette proteins which hydrolyze ATP to catalyse the translocation of their substrates through membranes. This protein extrudes a large range of components out of cells, especially therapeutic agents causing a phenomenon known as multidrug resistance. Because of its clinical interest, its activity and transport function have been largely characterized by various biochemical studies. In the absence of a high-resolution structure of P-glycoprotein, homology modeling is a useful tool to help interpretation of experimental data and potentially guide experimental studies.

Results

We present here three-dimensional models of two different catalytic states of P-glycoprotein that were developed based on the crystal structures of two bacterial multidrug transporters. Our models are supported by a large body of biochemical data. Measured inter-residue distances correlate well with distances derived from cross-linking data. The nucleotide-free model features a large cavity detected in the protein core into which ligands of different size were successfully docked. The locations of docked ligands compare favorably with those suggested by drug binding site mutants.

Conclusion

Our models can interpret the effects of several mutants in the nucleotide-binding domains (NBDs), within the transmembrane domains (TMDs) or at the NBD:TMD interface. The docking results suggest that the protein has multiple binding sites in agreement with experimental evidence. The nucleotide-bound models are exploited to propose different pathways of signal transmission upon ATP binding/hydrolysis which could lead to the elaboration of conformational changes needed for substrate translocation. We identified a cluster of aromatic residues located at the interface between the NBD and the TMD in opposite halves of the molecule which may contribute to this signal transmission.

Our models may characterize different steps in the catalytic cycle and may be important tools to understand the structure-function relationship of P-glycoprotein.

Processing mutations disrupt interactions between the nucleotide binding and transmembrane domains of P-glycoprotein and the cystic fibrosis transmembrane conductance regulator (CFTR)

P-glycoprotein (P-gp, ABCB1) is an ATP-dependent drug pump. Each of its two homologous halves contains a transmembrane domain (TMD) that has six transmembrane (TM) segments and a nucleotide-binding domain (NBD). Determining how the two halves interact may provide insight into the folding of P-gp as the drug-binding pocket and nucleotide-binding sites are predicted to be at the interface between the two halves. Here, we present evidence for NBD1-TMD2 and NBD2-TMD1 interactions. We also show that TMD-NBD interactions in immature and mature P-gp can be affected by the presence of a processing mutation. We found that the NBD-TMD mutants L443C(NBD1)/S909C(TMD2) and A266C(TMD1)/F1086C(NBD2) could be cross-linked at 0 degrees C with oxidant (copper phenanthroline). Cross-linking was inhibited by vanadate-trapping of nucleotide. The presence of a processing mutation (G268V/L443C(NBD1)/S909C(TMD2); L1260A/A266C(TMD1)/F1086C(NBD2)) resulted in the synthesis of the immature (150 kDa) protein as the major product and the mutants could not be cross-linked with copper phenanthroline. Expression of the processing mutants in the presence of a pharmacological chaperone (cyclosporin A), however, resulted in the expression of mature (170 kDa) protein at the cell surface that could be cross-linked. Similarly, CFTR mutants A274C(TMD1)/L1260C(NBD2) and V510C(NBD1)/A1067C(TMD2) could be cross-linked at 0 degrees C with copper phenanthroline. Introduction of DeltaF508 mutation in these mutants, however, resulted in the synthesis of immature CFTR that could not be cross-linked. These results suggest that establishment of NBD interactions with the opposite TMD is a key step in folding of ABC transporters.

Chalcogenopyrylium dyes as inhibitors/modulators of P-glycoprotein in multidrug-resistant cells

A series of chalcogenopyrylium dyes were evaluated as modulators/inhibitors of P-glycoprotein (Pgp). Their ability to inhibit verapamil (VER)-dependent ATPase activity (IC(50) values) in lipid-activated, mouse Cys-less mdr3 Pgp was determined. Their ability to promote calcein-AM (CAM) uptake in MDCKII-MDR1 cells and their capacity to be transported by Pgp in monolayers of MDCKII-MDR1 cells were also evaluated. The chalcogenopyrylium dyes promoted CAM uptake with values of EC(50) between 5 x 10(-6) and 3.5 x 10(-5)M and 7 of the 9 dyes examined in transport studies were substrates for Pgp with efflux ratios (P(BA/AB)) between 14 and 390. Binding of three compounds (1-S, 3-S, and 4-S) to Pgp was also assessed by fluorescence. These three thiopyrylium dyes showed increased fluorescence upon binding to Pgp, giving apparent binding constants, K(app), on the order of 10(-7) to 10(-6)M. Compound 8-Te was particularly intriguing since it appeared to influence Pgp at low micromolar concentrations as evidenced by its influence on VER-stimulated ATPase activity (IC(50) of 1.2 x 10(-6)M), CAM uptake (EC(50) of 5.4 x 10(-6)M), as well as [(3)H]-vinblastine transport by Pgp in cells (IC(50) of 4.3 x 10(-6)M) and within inside-out membrane vesicles (IC(50) of 9.6 x 10(-6)M). Yet, Pgp did not influence the distribution of 8-Te in MDCKII-MDR1 monolayers suggesting that 8-Te may bind to an allosteric site.

Correctors promote folding of the CFTR in the endoplasmic reticulum

Cystic fibrosis (CF) is most commonly caused by deletion of a residue (DeltaF508) in the CFTR (cystic fibrosis transmembrane conductance regulator) protein. The misfolded mutant protein is retained in the ER (endoplasmic reticulum) and is not trafficked to the cell surface (misprocessed mutant). Corrector molecules such as corr-2b or corr-4a are small molecules that increase the amount of functional CFTR at the cell surface. Correctors may function by stabilizing CFTR at the cell surface or by promoting folding in the ER. To test whether correctors promoted folding of CFTR in the ER, we constructed double-cysteine CFTR mutants that would be retained in the ER and only undergo cross-linking when the protein folds into a native structure. The mature form, but not the immature forms, of M348C(TM6)/T1142C(TM12) (where TM is transmembrane segment), T351C(TM6)/T1142C(TM12) and W356C(TM6)/W1145C(TM12) mutants were efficiently cross-linked. Mutations to the COPII (coatamer protein II) exit motif (Y(563)KDAD(567)) were then made in the cross-linkable cysteine mutants to prevent the mutant proteins from leaving the ER. Membranes were prepared from the mutants expressed in the absence or presence of correctors and subjected to disulfide cross-linking analysis. The presence of correctors promoted folding of the mutants as the efficiency of cross-linking increased from approx. 2-5% to 22-35%. The results suggest that correctors interact with CFTR in the ER to promote folding of the protein into a native structure.

Arginines in the first transmembrane segment promote maturation of a P-glycoprotein processing mutant by hydrogen bond interactions with tyrosines in transmembrane segment 11

A key goal is to correct defective folding of mutant ATP binding cassette (ABC) transporters, as they cause diseases such as cystic fibrosis. P-glycoprotein (ABCB1) is a useful model system because introduction of an arginine at position 65 of the first transmembrane (TM) segment could repair folding defects. To determine the mechanism of arginine rescue, we first tested the effects of introducing arginines at other positions in TM1 (residues 52-72) of a P-glycoprotein processing mutant (G251V) that is defective in folding and trafficking to the cell surface (20% maturation efficiency). We found that arginines introduced into one face of the TM1 helix (positions 52, 55, 56, 59, 60, 62, 63, 66, and 67) inhibited maturation, whereas arginines on the opposite face of the helix promoted (positions 64, 65, 68, and 71) or had little effect (positions 61, and 69) on maturation. Arginines at positions 61, 64, 65, and 68 appeared to lie close to the drug binding sites as they reduced the apparent affinity for drug substrates such as vinblastine and verapamil. Therefore, arginines that promoted maturation may face an aqueous drug translocation pathway, whereas those that inhibited maturation may face the lipid bilayer. The highest maturation efficiencies (60-85%) were observed with the Arg-65 and Arg-68 mutants. Mutations that removed hydrogen bond acceptors (Y950F/Y950A or Y953F/Y953A) in TM11 predicted to lie close to Arg-65 or Arg-68 inhibited maturation but did not affect maturation of the G251V parent. Therefore, arginine may rescue defective folding by promoting packing of the TM segments through hydrogen bond interactions.