Predicting binding to p-glycoprotein by flexible receptor docking

The pH-dependence of efflux ratios determined with bidirectional transport assays across cellular monolayers

MDCK/Caco-2 assays serve as essential in vitro tools for evaluating membrane permeability and active transport, especially mediated by P-glycoprotein (P-gp). Despite their utility, challenges remain in quantifying active transport and using the efflux ratio (ER) to determine intrinsic values for active efflux. Such an intrinsic value for P-gp facilitated efflux necessitates knowing whether this transporter transports the neutral or ionic species of a compound. Utilising MDCK-MDR1 assays, we investigate a method for determining transporter substrate fraction preference by studying ER pH-dependence for basic, acidic and non-dissociating compounds. These results are compared with model fits based on various assumptions of transporter species preference. As an unexpected consequence of these assays, we also give evidence for an additional influx transporter at the basolateral membrane, and further extend our model to incorporate this transport. The combined influences of paracellular transport, the previously unaccounted for basolateral influx transporter, as well as potential pH effects on the transporter impedes the extraction of intrinsic values for active transport from the ER. Furthermore, we determined that using inhibitor affects the measurement of paracellular transport. While clear indications of transporter species preference remain elusive, this study enhances understanding of the MDCK system.

Model of P-Glycoprotein Ligand Binding and Validation with Efflux Substrate Matched Pairs

The blood-brain barrier (BBB) poses a significant obstacle in developing therapeutics for neurodegenerative diseases and central nervous system (CNS) disorders. P-glycoprotein (P-gp), a multidrug resistance protein, is a critical gatekeeper in the BBB and plays a role in cancer chemoresistance. This paper uses cryo-EM P-gp structures as starting points with an induced fit docking (IFD) model to evaluate 19 pairs of compounds with known P-gp efflux data. The study reveals significant differences in binding energy and sheds light on structural modifications' impact on efflux properties. In the cases examined, fluorine incorporation influences the efflux by altering the molecular conformation rather than proximal heteroatom basicity. Although there are limitations in addressing covalent interactions or when binding extends into the more flexible vestibule region of the protein, the results provide valuable insights and potential strategies to overcome P-gp efflux, contributing to the advancement of drug development for both CNS disorders and cancer therapies.

Ligand- and Structure-based Approaches for Transmembrane Transporter Modeling

The study of transporter proteins is key to understanding the mechanism behind multidrug resistance and drug-drug interactions causing severe side effects. While ATP-binding transporters are well-studied, solute carriers illustrate an understudied family with a high number of orphan proteins. To study these transporters, in silico methods can be used to shed light on the basic molecular machinery by studying protein-ligand interactions. Nowadays, computational methods are an integral part of the drug discovery and development process. In this short review, computational approaches, such as machine learning, are discussed, which try to tackle interactions between transport proteins and certain compounds to locate target proteins. Furthermore, a few cases of selected members of the ATP binding transporter and solute carrier family are covered, which are of high interest in clinical drug interaction studies, especially for regulatory agencies. The strengths and limitations of ligand-based and structure-based methods are discussed to highlight their applicability for different studies. Furthermore, the combination of multiple approaches can improve the information obtained to find crucial amino acids that explain important interactions of protein-ligand complexes in more detail. This allows the design of drug candidates with increased activity towards a target protein, which further helps to support future synthetic efforts.

An Oral 3D Printed PLGA-Tocopherol PEG Succinate Nanocomposite Hydrogel for High-Dose Methotrexate Delivery in Maintenance Chemotherapy

High-dose methotrexate (HDMTX) is one of the chemotherapeutic agents used to treat a variety of cancers in both adults and children. However, the toxicity associated with HDMTX has resulted in the spread of infections and treatment interruption. Further, poor bioavailability due to efflux pump activities mediated by P-glycoprotein has also been linked to poor therapeutic effects of methotrexate following oral administrations. D-α-Tocopheryl poly-ethylene glycol 1000 succinate (TPGS) is known to improve the bioavailability of poorly soluble drugs by inhibiting P-gp efflux activities, thus enhancing cellular uptake. Therefore, to achieve improved bioavailability for MTX, this study aimed to design and develop a novel drug delivery system employing TPGS and a biodegradable polymer, i.e., PLGA, to construct methotrexate-loaded nanoparticles fixated in alginate-gelatine 3D printable hydrogel ink to form a solid 3D printed tablet for oral delivery. The results indicated that high accuracy (>95%) of the 3D printed tablets was achieved using a 25 G needle. In vitro, drug release profiles were investigated at pH 1.2 and pH 7.4 to simulate the gastrointestinal environment. The in vitro release profile displayed a controlled and prolonged release of methotrexate over 24 h. The in silico modeling study displayed P-gp ATPase inhibition, suggesting enhanced MTX absorption from the gastrointestinal site. The 3D-printed hydrogel-based tablet has the potential to overcome the chemotherapeutic challenges that are experienced with conventional therapies.

Systematic Exploration of Passive Permeability in Tetrapeptides with Hydrogen---Bond Accepting Amino Acid Side Chains

We synthesized and experimentally tested the passive permeability of more than thirty tetrapeptides mimicking the N -terminus of the pro-apoptotic protein Smac (Second mitochondria-derived activator of caspases). Each peptide bore one or two unnatural Hydrogen Bond Acceptor-bearing Amino Acid (HBA-AA) residues, such that intramolecular hydrogen bonding with proximal backbone amide N-H donors is feasible. Passive permeability of the synthetic peptides was determined using the parallel artificial membrane permeability assay (PAMPA). Experimental permeability values were found to span three orders of magnitude, providing useful empirical guidance for the design of more permeable Smac mimetics specifically, and peptidic ligands generally.

Structure-Function Relationships in the Human P-Glycoprotein (ABCB1): Insights from Molecular Dynamics Simulations

P-Glycoprotein (P-gp) is a transmembrane protein belonging to the ATP binding cassette superfamily of transporters, and it is a xenobiotic efflux pump that limits intracellular drug accumulation by pumping compounds out of cells. P-gp contributes to a reduction in toxicity, and has broad substrate specificity. It is involved in the failure of many cancer and antiviral chemotherapies due to the phenomenon of multidrug resistance (MDR), in which the membrane transporter removes chemotherapeutic drugs from target cells. Understanding the details of the ligand–P-gp interaction is therefore critical for the development of drugs that can overcome the MDR phenomenon, for the early identification of P-gp substrates that will help us to obtain a more effective prediction of toxicity, and for the subsequent outdesign of substrate properties if needed. In this work, a series of molecular dynamics (MD) simulations of human P-gp (hP-gp) in an explicit membrane-and-water environment were performed to investigate the effects of binding different compounds on the conformational dynamics of P-gp. The results revealed significant differences in the behaviour of P-gp in the presence of active and non-active compounds within the binding pocket, as different patterns of movement were identified that could be correlated with conformational changes leading to the activation of the translocation mechanism. The predicted ligand–P-gp interactions are in good agreement with the available experimental data, as well as the estimation of the binding-free energies of the studied complexes, demonstrating the validity of the results derived from the MD simulations.

Triterpenoids with modified A-ring as modulators of P-gp-dependent drug-resistance in cancer cells

Semi-synthetic A-cycle modified triterpenic derivatives with A-cycle condensed with a heterocyclic fragment (compound 1) and fragmented A-ring (compound 2) were tested for cytotoxicity against several tumor cell cultures and doxorubicin (Dox)-resistant cell lines. The equal cytotoxicity of the tested compounds to the parental tumor cell lines (HBL-100, K562) and their resistant subclones (HBL-100/Dox, K562/i-S9) was revealed. The overexpression of ABCB1 (MDR1) gene and P-glycoprotein (P-gp) was confirmed for both resistant subclones of tumor cells. Compounds 1 and 2 were shown to inhibit the ABC-transporter gene expression (MDR1, MRP, MVP, and BCRP) and the transport of well-known P-gp substrate Rhodamine 123 from resistant cells. The docking of triterpenoids 1 and 2 into the drug binding site of P-gp revealed a similarity between the conformation of the tested triterpenoids and that of classical inhibitor verapamil, thus assuming these compounds to be more likely the inhibitors than the substrates of P-gp. Any tested triterpenic derivatives, when combined at non-toxic concentrations with doxorubicin, improved cytotoxic effect of the therapeutic drug against resistant subclones of tumor cells.

Molecular docking study of lignanamides from Cannabis sativa against P-glycoprotein

P-glycoprotein (P-gp), which was first identified in cancer cells, is an ATP-dependent efflux transporter that expels a wide variety of cytotoxic compounds out of cells. This transporter can decrease the bioavailability of therapeutic drugs by preventing their sufficient intracellular accumulation. Over expression of P-gp in cancer cells lead to multidrug resistance (MDR) phenotype that is one of the main reasons for the failure of chemotherapy. Hence, P-gp inhibition is a favorable method to reverse MDR. In this study, the lignanamides from Cannabis sativa were docked against P-gp to recognize potential binding affinities of these phytochemicals. Tariquidar and zosuquidar, two well-known P-gp inhibitors, were selected as the control ligands. It was observed that cannabisin M and cannabisin N exhibited higher binding affinities (- 10.2 kcal/mol) to drug-binding pocket of P-gp when compared with tariquidar and zosuquidar that showed binding affinities of - 10.1 and - 9.6 kcal/mol, respectively. Based on these findings, cannabisin M and cannabisin N could be good drug candidates against P-gp.

Theoretical insights on helix repacking as the origin of P-glycoprotein promiscuity

P-glycoprotein (P-gp, ABCB1) overexpression is, currently, one of the most important multidrug resistance (MDR) mechanisms in tumor cells. Thus, modulating drug efflux by P-gp has become one of the most promising approaches to overcome MDR in cancer. Yet, more insights on the molecular basis of drug specificity and efflux-related signal transmission mechanism between the transmembrane domains (TMDs) and the nucleotide binding domains (NBDs) are needed to develop molecules with higher selectivity and efficacy. Starting from a murine P-gp crystallographic structure at the inward-facing conformation (PDB ID: 4Q9H), we evaluated the structural quality of the herein generated human P-gp homology model. This initial human P-gp model, in the presence of the “linker” and inserted in a suitable lipid bilayer, was refined through molecular dynamics simulations and thoroughly validated. The best human P-gp model was further used to study the effect of four single-point mutations located at the TMDs, experimentally related with changes in substrate specificity and drug-stimulated ATPase activity. Remarkably, each P-gp mutation is able to induce transmembrane α-helices (TMHs) repacking, affecting the drug-binding pocket volume and the drug-binding sites properties (e.g. volume, shape and polarity) finally compromising drug binding at the substrate binding sites. Furthermore, intracellular coupling helices (ICH) also play an important role since changes in the TMHs rearrangement are shown to have an impact in residue interactions at the ICH-NBD interfaces, suggesting that identified TMHs repacking affect TMD-NBD contacts and interfere with signal transmission from the TMDs to the NBDs.

Search for ABCB1 Modulators Among 2-Amine-5-Arylideneimidazolones as a New Perspective to Overcome Cancer Multidrug Resistance

Multidrug resistance (MDR) is a severe problem in the treatment of cancer with overexpression of glycoprotein P (Pgp, ABCB1) as a reason for chemotherapy failure. A series of 14 novel 5-arylideneimidazolone derivatives containing the morpholine moiety, with respect to two different topologies (groups A and B), were designed and obtained in a three- or four-step synthesis, involving the Dimroth rearrangement. The new compounds were tested for their inhibition of the ABCB1 efflux pump in both sensitive (parental (PAR)) and ABCB1-overexpressing (MDR) T-lymphoma cancer cells in a rhodamine 123 accumulation assay. Their cytotoxic and antiproliferative effects were investigated by a thiazolyl blue tetrazolium bromide (MTT) assay. For active compounds, an insight into the mechanisms of action using either the luminescent Pgp-Glo™ Assay in vitro or docking studies to human Pgp was performed. The safety profile in vitro was examined. Structure–activity relationship (SAR) analysis was discussed. The most active compounds, representing both 2-substituted- (11) and Dimroth-rearranged 3-substituted (18) imidazolone topologies, displayed 1.38–1.46 fold stronger efflux pump inhibiting effects than reference verapamil and were significantly safer than doxorubicin in cell-based toxicity assays in the HEK-293 cell line. Results of mechanistic studies indicate that active imidazolones are substrates with increasing Pgp ATPase activity, and their dye-efflux inhibition via competitive action on the Pgp verapamil binding site was predicted in silico.

Potential synergic mechanism of Wutou-Gancao herb-pair by inhibiting efflux transporter P-glycoprotein

Wutou-Gancao herb-pair is extensively used to attenuate the toxicity and enhance the efficacy of aconite. In this study, potential synergic mechanism of the herb pair was investigated by utilizing multiple approaches. In silico and in vitro Caco-2 cell models were applied to study the potential binding mode of bioactive ingredients existing in liquorice with P-glycoprotein (P-gp), as well as the inhibition effects on P-gp. Additionally, anti-inflammatory activity of aconitine (AC) combined with active ingredients of liquorice, as well as pharmacokinetic patterns of AC after co-administration was investigated. Anti-inflammatory effect of AC (1 mg/kg) in rats was enhanced in combination with bioactive ingredients of liquorice (10 mg/kg). In the meanwhile, the exposure of AC in vivo was altered, in terms of Cmax and AUC. For instance, the Cmax and AUC were increased to 1.9 and 1.3 folds, respectively, when used in combination with liquiritigenin. The in silico study revealed the potential binding mode with outward facing conformation of P-gp. The resulting data obtained from transport of rhodamine-123 (Rh-123) across Caco-2 cell monolayer further indicated that the function of P-gp was inhibited by chemicals in liquorice. The synergic effect was therefore proposed to be attributed to inhibition of P-gp by liquorice since AC has been demonstrated to be the substrate of P-gp. The resuls revealed that potential synergic mechanism of Wutou-Gancao herb-pair by inhibiting function of key efflux transporter P-gp to enhance the exposure of AC in systematic circulation, and further the anti-inflammatory effect, which helps clarify the compatibility rationale of these two herbs.

Multiclass Classifier for P-Glycoprotein Substrates, Inhibitors, and Non-Active Compounds

P-glycoprotein (P-gp) is a transmembrane protein that actively transports a wide variety of chemically diverse compounds out of the cell. It is highly associated with the ADMET (absorption, distribution, metabolism, excretion and toxicity) properties of drugs/drug candidates and contributes to decreasing toxicity by eliminating compounds from cells, thereby preventing intracellular accumulation. Therefore, in the drug discovery and toxicological assessment process it is advisable to pay attention to whether a compound under development could be transported by P-gp or not. In this study, an in silico multiclass classification model capable of predicting the probability of a compound to interact with P-gp was developed using a counter-propagation artificial neural network (CP ANN) based on a set of 2D molecular descriptors, as well as an extensive dataset of 2512 compounds (1178 P-gp inhibitors, 477 P-gp substrates and 857 P-gp non-active compounds). The model provided a good classification performance, producing non error rate (NER) values of 0.93 for the training set and 0.85 for the test set, while the average precision (AvPr) was 0.93 for the training set and 0.87 for the test set. An external validation set of 385 compounds was used to challenge the model’s performance. On the external validation set the NER and AvPr values were 0.70 for both indices. We believe that this in silico classifier could be effectively used as a reliable virtual screening tool for identifying potential P-gp ligands.

Different structures of berberine and five other protoberberine alkaloids that affect P-glycoprotein-mediated efflux capacity

Berberine, berberrubine, thalifendine, demethyleneberberine, jatrorrhizine, and columbamine are six natural protoberberine alkaloid (PA) compounds that display extensive pharmacological properties and share the same protoberberine molecular skeleton with only slight substitution differences. The oral delivery of most PAs is hindered by their poor bioavailability, which is largely caused by P-glycoprotein (P-gp)-mediated drug efflux. Meanwhile, P-gp undergoes large-scale conformational changes (from an inward-facing to an outward-facing state) when transporting substrates, and these changes might strongly affect the P-gp-binding specificity. To confirm whether these six compounds are substrates of P-gp, to investigate the differences in efflux capacity caused by their trivial structural differences and to reveal the key to increasing their binding affinity to P-gp, we conducted a series of in vivo, in vitro, and in silico assays. Here, we first confirmed that all six compounds were substrates of P-gp by comparing the drug concentrations in wild-type and P-gp-knockout mice in vivo. The efflux capacity (net efflux) ranked as berberrubine > berberine > columbamine ~ jatrorrhizine > thalifendine > demethyleneberberine based on in vitro transport studies in Caco-2 monolayers. Using molecular dynamics simulation and molecular docking techniques, we determined the transport pathways of the six compounds and their binding affinities to P-gp. The results suggested that at the early binding stage, different hydrophobic and electrostatic interactions collectively differentiate the binding affinities of the compounds to P-gp, whereas electrostatic interactions are the main determinant at the late release stage. In addition to hydrophobic interactions, hydrogen bonds play an important role in discriminating the binding affinities.

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.

Theoretical Prediction of the Complex P-Glycoprotein Substrate Efflux Based on the Novel Hierarchical Support Vector Regression Scheme

P-glycoprotein (P-gp), a membrane-bound transporter, can eliminate xenobiotics by transporting them out of the cells or blood⁻brain barrier (BBB) at the expense of ATP hydrolysis. Thus, P-gp mediated efflux plays a pivotal role in altering the absorption and disposition of a wide range of substrates. Nevertheless, the mechanism of P-gp substrate efflux is rather complex since it can take place through active transport and passive permeability in addition to multiple P-gp substrate binding sites. A nonlinear quantitative structure⁻activity relationship (QSAR) model was developed in this study using the novel machine learning-based hierarchical support vector regression (HSVR) scheme to explore the perplexing relationships between descriptors and efflux ratio. The predictions by HSVR were found to be in good agreement with the observed values for the molecules in the training set (n = 50, r² = 0.96, qCV2 = 0.94, RMSE = 0.10, s = 0.10) and test set (n = 13, q² = 0.80⁻0.87, RMSE = 0.21, s = 0.22). When subjected to a variety of statistical validations, the developed HSVR model consistently met the most stringent criteria. A mock test also asserted the predictivity of HSVR. Consequently, this HSVR model can be adopted to facilitate drug discovery and development.

Coarse-grained molecular dynamics simulations reveal lipid access pathways in P-glycoprotein

P-glycoprotein (P-gp) exports a broad range of dissimilar compounds, including drugs, lipids, and lipid-like molecules. Because of its substrate promiscuity, P-gp is a key player in the development of cancer multidrug resistance. Although P-gp is one of the most studied ABC transporters, the mechanism by which its substrates access the cavity remains unclear. In this study, we perform coarse-grained molecular dynamics simulations to explore possible lipid access pathways in the inward-facing conformation of P-gp embedded in bilayers of different lipid compositions. In the inward-facing orientation, only lipids from the lower leaflet access the cavity of the transporter. We identify positively charged residues at the portals of P-gp that favor lipid entrance to the cavity, as well as lipid-binding sites at the portals and within the cavity, which is in good agreement with previous experimental studies. This work includes several examples of lipid pathways for phosphatidylcholine and phosphatidylethanolamine lipids that help elucidate the molecular mechanism of lipid binding in P-gp.

Selected anthraquinones as potential free radical scavengers and P-glycoprotein inhibitors

In this article, we estimated the scavenger capacity of six selected anthraquinones toward free radicals and their efficacy as inhibitors of P-glycoproteins.

Computational modeling in glioblastoma: from the prediction of blood-brain barrier permeability to the simulation of tumor behavior

The integrated in silico-in vitro-in vivo approaches have fostered the development of new treatment strategies for glioblastoma patients and improved diagnosis, establishing the bridge between biochemical research and clinical practice. These approaches have provided new insights on the identification of bioactive compounds and on the complex mechanisms underlying the interactions among glioblastoma cells, and the tumor microenvironment. This review focuses on the key advances pertaining to computational modeling in glioblastoma, including predictive data on drug permeability across the blood-brain barrier, tumor growth and treatment responses. Structure- and ligand-based methods have been widely adopted, enabling the study of dynamic and evolutionary aspects of glioblastoma. Their potential applications as predictive tools and the advantages over other well-known methodologies are outlined. Challenges regarding in silico approaches for predicting tumor properties are also discussed.

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.

Emodin reverses leukemia multidrug resistance by competitive inhibition and downregulation of P-glycoprotein

Development of multidrug resistance (MDR) is a continuous clinical challenge partially due to the overexpression of P-glycoprotein (P-gp) for chronic myelogenous leukemia (CML) patients. Herein, we evaluated the inhibitory potency of emodin, a natural anthraquinone derivative isolated from Rheum palmatum L, on P-gp in P-gp positive K562/ADM cells. Competition experiments combined with molecular docking analysis were utilized to investigate the binding modes between emodin and binding sites of P-gp. Emodin reversed adriamycin resistance in K562/ADM cells accompanied with the decrease of P-gp protein expression, further increasing the uptake of rhodamine123 in both K562/ADM and Caco-2 cells, indicating the inhibition of P-gp efflux function. Moreover, when incubated with emodin under different conditions where P-gp was inhibited, K562/ADM cells displayed increasing intracellular uptake of emodin, suggesting that emodin may be the potential substrate of P-gp. Importantly, rhodamine 123 could increase the K_intrinsic (K_i) value of emodin linearly, whereas, verapamil could not, implying that emodin competitively bound to the R site of P-gp and noncompetition existed between emodin and verapamil at the M site, in a good accordance with the results of molecular docking that emodin bound to the R site of P-gp with higher affinity. Based on our results, we suggest that emodin might be used to modulate P-gp function and expression.

Evolutionarily Conserved Roles for Blood-Brain Barrier Xenobiotic Transporters in Endogenous Steroid Partitioning and Behavior

Central nervous system (CNS) chemical protection depends upon discrete control of small-molecule access by the blood-brain barrier (BBB). Curiously, some drugs cause CNS side-effects despite negligible transit past the BBB. To investigate this phenomenon, we asked whether the highly BBB-enriched drug efflux transporter MDR1 has dual functions in controlling drug and endogenous molecule CNS homeostasis. If this is true, then brain-impermeable drugs could induce behavioral changes by affecting brain levels of endogenous molecules. Using computational, genetic, and pharmacologic approaches across diverse organisms, we demonstrate that BBB-localized efflux transporters are critical for regulating brain levels of endogenous steroids and steroid-regulated behaviors (sleep in Drosophila and anxiety in mice). Furthermore, we show that MDR1-interacting drugs are associated with anxiety-related behaviors in humans. We propose a general mechanism for common behavioral side effects of prescription drugs: pharmacologically challenging BBB efflux transporters disrupts brain levels of endogenous substrates and implicates the BBB in behavioral regulation.

Schiff base derived from thiosemicarbazone and anthracene showed high potential in overcoming multidrug resistance in vitro with low drug resistance index

Multidrug resistance (MDR) is a huge obstacle in cancer chemotherapeutics. Overcoming MDR is a great challenge for anticancer drug discovery. Here, DNA binding and cytotoxicity of Schiff base L1 and L2 were explored to assess their efficiency in fighting cancer and overcoming the MDR. L1 and L2 could treat extremely chemoresistant MCF-7/ADR cell as drug-sensitive cell, with drug resistance index (DRI) <2.13, showing high potential in overcoming the MDR. The apoptotic ratio induced by L1 and L2 was low for both MCF-7 and MCF-7/ADR cells. L1 and L2 induced an impairment of cell cycle progression of MCF-7 and MCF-7/ADR cell lines and suppressed cell growth by perturbing progress through the G0/G1 phase, with L2 causing more profound effect, which might account for lower drug resistance after L2 treatment. The molecular docking revealed weak interaction between L1/L2 and P-glycoprotein (P-gp), the most important drug efflux pump and intracellular Rhodamine 123 accumulation indicated that the activity of P-gp was not inhibited by L1 and L2. Combined with the cellular uptake results, it implied that L1 and L2 could bypass P-gp efflux to exert anticancer activity.

Molecular properties associated with transporter-mediated drug disposition

Membrane transporters play a key role in the absorption, distribution, clearance, elimination, and transport of drugs. Understanding the drug properties and structure activity relationships (SAR) for affinity to membrane transporters is critical to optimize clearance and pharmacokinetics during drug design. To facilitate the early identification of clearance mechanism, a framework named the extended clearance classification system (ECCS) was recently introduced. Using in vitro and physicochemical properties that are readily available in early drug discovery, ECCS has been successfully applied to identify major clearance mechanism and to implicate the role of membrane transporters in determining pharmacokinetics. While the crystal structures for most of the drug transporters are currently not available, ligand-based modeling approaches that use information obtained from the structure and molecular properties of the ligands have been applied to associate the drug-related properties and transporter-mediated disposition. The approach allows prospective prediction of transporter both substrate and/or inhibitor affinity and build quantitative structure-activity relationship (QSAR) to enable early optimization of pharmacokinetics, tissue distribution and drug-drug interaction risk. Drug design applications can be further improved through uncovering transporter protein crystal structure and generation of quality data to refine and develop viable predictive models.

Structure based classification for bile salt export pump (BSEP) inhibitors using comparative structural modeling of human BSEP

The bile salt export pump (BSEP) actively transports conjugated monovalent bile acids from the hepatocytes into the bile. This facilitates the formation of micelles and promotes digestion and absorption of dietary fat. Inhibition of BSEP leads to decreased bile flow and accumulation of cytotoxic bile salts in the liver. A number of compounds have been identified to interact with BSEP, which results in drug-induced cholestasis or liver injury. Therefore, in silico approaches for flagging compounds as potential BSEP inhibitors would be of high value in the early stage of the drug discovery pipeline. Up to now, due to the lack of a high-resolution X-ray structure of BSEP, in silico based identification of BSEP inhibitors focused on ligand-based approaches. In this study, we provide a homology model for BSEP, developed using the corrected mouse P-glycoprotein structure (PDB ID: 4M1M). Subsequently, the model was used for docking-based classification of a set of 1212 compounds (405 BSEP inhibitors, 807 non-inhibitors). Using the scoring function ChemScore, a prediction accuracy of 81% on the training set and 73% on two external test sets could be obtained. In addition, the applicability domain of the models was assessed based on Euclidean distance. Further, analysis of the protein-ligand interaction fingerprints revealed certain functional group-amino acid residue interactions that could play a key role for ligand binding. Though ligand-based models, due to their high speed and accuracy, remain the method of choice for classification of BSEP inhibitors, structure-assisted docking models demonstrate reasonably good prediction accuracies while additionally providing information about putative protein-ligand interactions.

A stronger reversal effect of the combination of dasatinib and menadione on P-gp-mediated multidrug resistance in human leukemia K562/Adr cell line

Multidrug resistance (MDR) leads to poor efficiency of chemotherapy.

QSAR Model of Unbound Brain-to-Plasma Partition Coefficient, Kp,uu,brain: Incorporating P-glycoprotein Efflux as a Variable

We report development and prospective validation of a QSAR model of the unbound brain-to-plasma partition coefficient, K_p,uu,brain, based on the in-house data set of ∼1000 compounds. We discuss effects of experimental variability, explore the applicability of both regression and classification approaches, and evaluate a novel, model-within-a-model approach of including P-glycoprotein efflux prediction as an additional variable. When tested on an independent test set of 91 internal compounds, incorporation of P-glycoprotein efflux information significantly improves the model performance resulting in an R² of 0.53, RMSE of 0.57, Spearman's Rho correlation coefficient of 0.73, and qualitative prediction accuracy of 0.8 (kappa = 0.6). In addition to improving the performance, one of the key advantages of this approach is the larger chemical space coverage provided indirectly through incorporation of the in vitro, higher throughput data set that is 4 times larger than the in vivo data set.

Machine learning-, rule- and pharmacophore-based classification on the inhibition of P-glycoprotein and NorA

The efflux pumps P-glycoprotein (P-gp) in humans and NorA in Staphylococcus aureus are of great interest for medicinal chemists because of their important roles in multidrug resistance (MDR). The high polyspecificity as well as the unavailability of high-resolution X-ray crystal structures of these transmembrane proteins lead us to combining ligand-based approaches, which in the case of this study were machine learning, perceptual mapping and pharmacophore modelling. For P-gp inhibitory activity, individual models were developed using different machine learning algorithms and subsequently combined into an ensemble model which showed a good discrimination between inhibitors and noninhibitors (acctrain-diverse = 84%; accinternal-test = 92% and accexternal-test = 100%). For ligand promiscuity between P-gp and NorA, perceptual maps and pharmacophore models were generated for the detection of rules and features. Based on these in silico tools, hit compounds for reversing MDR were discovered from the in-house and DrugBank databases through virtual screening in an attempt to restore drug sensitivity in cancer cells and bacteria.

Designing the angiogenic inhibitor for brain tumor via disruption of VEGF and IL17A expression

Glioblastoma multiforme is a highly malignant, heterogenic, and drug resistant tumor. The blood-brain barrier (BBB), systemic cytotoxicity, and limited specificity are the main obstacles in designing brain tumor drugs. In this study a computational approach was used to design brain tumor drugs that could downregulate VEGF and IL17A in glioblastoma multiforme type four. Computational screening tools were used to evaluate potential candidates for antiangiogenic activity, target binding, BBB permeability, and ADME physicochemical properties. Additionally, in vitro cytotoxicity, migration, invasion, tube formation, apoptosis, ROS and ELISA assays were conducted for molecule 6 that was deemed most likely to succeed. The efflux ratio of membrane permeability and calculated docking scores of permeability to glycoproteins (P-gps) were used to determine the BBB permeability of the molecules. The results showed BBB permeation for molecule 6, with the predicted efficiency of 0.55kcal/mol and binding affinity of -37kj/mol corresponding to an experimental efflux ratio of 0.625 and predicted -15kj/mol of binding affinity for P-gps. Molecule 6 significantly affected the angiogenesis pathways by 2-fold downregulation of IL17A and VEGF through inactivation of active sites of HSP90 (predicted binding: -37kj/mol, predicted efficiency: 0.55kcal/mol) and p23 (predicted binding: 12kj/mol, predicted efficiency: 0.17kcal/mol) chaperon proteins. Additionally, molecule 6 activated the 17.38% relative fold of ROS level at 18.3μg/mL and upregulated the caspase which lead the potential synergistic apoptosis through the antiangiogenic activity of molecule 6 and thereby the highly efficacious anticancer upshot. The results indicate that the binding of the molecules to the therapeutic target is not essential to produce a lethal effect on cancer cells of the brain and that antiangiogenic efficiency is much more important.

Microscopic Characterization of Membrane Transporter Function by In Silico Modeling and Simulation

Membrane transporters mediate one of the most fundamental processes in biology. They are the main gatekeepers controlling active traffic of materials in a highly selective and regulated manner between different cellular compartments demarcated by biological membranes. At the heart of the mechanism of membrane transporters lie protein conformational changes of diverse forms and magnitudes, which closely mediate critical aspects of the transport process, most importantly the coordinated motions of remotely located gating elements and their tight coupling to chemical processes such as binding, unbinding and translocation of transported substrate and cotransported ions, ATP binding and hydrolysis, and other molecular events fueling uphill transport of the cargo. An increasing number of functional studies have established the active participation of lipids and other components of biological membranes in the function of transporters and other membrane proteins, often acting as major signaling and regulating elements. Understanding the mechanistic details of these molecular processes require methods that offer high spatial and temporal resolutions. Computational modeling and simulations technologies empowered by advanced sampling and free energy calculations have reached a sufficiently mature state to become an indispensable component of mechanistic studies of membrane transporters in their natural environment of the membrane. In this article, we provide an overview of a number of major computational protocols and techniques commonly used in membrane transporter modeling and simulation studies. The article also includes practical hints on effective use of these methods, critical perspectives on their strengths and weak points, and examples of their successful applications to membrane transporters, selected from the research performed in our own laboratory.

Enzastaurin inhibits ABCB1-mediated drug efflux independently of effects on protein kinase C signalling and the cellular p53 status

The PKCβ inhibitor enzastaurin was tested in parental neuroblastoma and rhabdomyosarcoma cell lines, their vincristine-resistant sub-lines, primary neuroblastoma cells, ABCB1-transduced, ABCG2-transduced, and p53-depleted cells. Enzastaurin IC₅₀s ranged from 3.3 to 9.5 μM in cell lines and primary cells independently of the ABCB1, ABCG2, or p53 status. Enzastaurin 0.3125 μM interfered with ABCB1-mediated drug transport. PKCα and PKCβ may phosphorylate and activate ABCB1 under the control of p53. However, enzastaurin exerted similar effects on ABCB1 in the presence or absence of functional p53. Also, enzastaurin inhibited PKC signalling only in concentrations ≥ 1.25 μM. The investigated cell lines did not express PKCβ. PKCα depletion reduced PKC signalling but did not affect ABCB1 activity. Intracellular levels of the fluorescent ABCB1 substrate rhodamine 123 rapidly decreased after wash-out of extracellular enzastaurin, and enzastaurin induced ABCB1 ATPase activity resembling the ABCB1 substrate verapamil. Computational docking experiments detected a direct interaction of enzastaurin and ABCB1. These data suggest that enzastaurin directly interferes with ABCB1 function. Enzastaurin further inhibited ABCG2-mediated drug transport but by a different mechanism since it reduced ABCG2 ATPase activity. These findings are important for the further development of therapies combining enzastaurin with ABC transporter substrates.

Essential oils from Inula japonica and Angelicae dahuricae enhance sensitivity of MCF-7/ADR breast cancer cells to doxorubicin via multiple mechanisms

ETHNOPHARMACOLOGICAL RELEVANCE

Angelicae dahurica (Hoffm.) Benth. & Hook.f.ex Franch. & Sav combined with Pueraria and Gastrodia elata Bl. combined with Inula japonica Thunb. are widely used in herb-pairs of traditional chinese medicine. Previous studies have shown that Angelicae dahuricae essential oil (ADO) enhanced puerarin internalization into ABCB1-overexpressed Caco-2 cells. These findings suggest the possibility that essential oils may enhance the absorption via certain mechanisms related to ABCB1 and reverse multidrug resistance (MDR).

AIM OF THE STUDY

ADO and essential oils from Inula japonica (IJO) may reverse ABCB1-mediated MDR, but this ability has not been investigated in detail in the well-established cancer cell lines. In this study, the underlying molecular mechanisms were further investigated to examine how IJO and ADO reverse MDR in the resistant human breast cancer cell line of MCF-7/ADR. Also this work may help uncover the conceivable compatibility mechanisms of above herb-pairs involved in ABCB1.

MATERIALS AND METHODS

The MDR human breast cancer MCF-7/ADR cells were treated with IJO, its sesquiterpene component isoalantolactone (ISO) or ADOat non- cytotoxic concentrations. The MDR ability was examined by measuring the sensitivity to doxorubicin (DOX), DOX accumulation and efflux, ABCB1 ATPase activity, ABCB1 expression, membrane fluidity, and stability and localization of lipid rafts and caveolae. Finally, the molecular modeling was performed to postulate how ISO interacts with ABCB1.

RESULTS

Treating MCF-7/ADR cells with IJ oil, ISO or AD oil reversed MDR 2- to 3-fold, without affecting the sensitivity of the non-MDR parental cell line. Mechanistic studies showed that these oils down-regulated mRNA and protein expression of ABCB1, and reduced the stability of lipid rafts in the cell membrane, which has previously been shown to reduce ABCB1-mediated transport. On the other hand, IJO, ISO and ADO did not inhibit ABCB1 ATPase activity, and fluorescence polarization experiments showed that low concentrations of the oils did not appear to alter membrane fluidity, unlike some MDR-reversing agents, ISO showed a higher docking score than verapamil but lower than dofequidar and tariquidar.

CONCLUSIONS

Our results suggest that IJO, ISO and ADO could reverse MDR by down-regulating ABCB1 expression and reducing lipid raft stability. These findings may be useful for developing safer and effective MDR reversal agents and also help find out the compatibility mechanisms.

Human-Mouse Chimeras with Normal Expression and Function Reveal That Major Domain Swapping Is Tolerated by P-Glycoprotein (ABCB1)

The efflux transporter P-glycoprotein (P-gp) plays a vital role in the transport of molecules across cell membranes and has been shown to interact with a panoply of functionally and structurally unrelated compounds. How human P-gp interacts with this large number of drugs has not been well understood, although structural flexibility has been implicated. To gain insight into this transporter's broad substrate specificity and to assess its ability to accommodate a variety of molecular and structural changes, we generated human-mouse P-gp chimeras by the exchange of homologous transmembrane and nucleotide-binding domains. High-level expression of these chimeras by BacMam- and baculovirus-mediated transduction in mammalian (HeLa) and insect cells, respectively, was achieved. There were no detectable differences between wild-type and chimeric P-gp in terms of cell surface expression, ability to efflux the P-gp substrates rhodamine 123, calcein-AM, and JC-1, or to be inhibited by the substrate cyclosporine A and the inhibitors tariquidar and elacridar. Additionally, expression of chimeric P-gp was able to confer a paclitaxel-resistant phenotype to HeLa cells characteristic of P-gp-mediated drug resistance. P-gp ATPase assays and photo-cross-linking with [(125)I]iodoarylazidoprazosin confirmed that transport and biochemical properties of P-gp chimeras were similar to those of wild-type P-gp, although differences in drug binding were detected when human and mouse transmembrane domains were combined. Overall, chimeras with one or two mouse P-gp domains were deemed functionally equivalent to human wild-type P-gp, demonstrating the ability of human P-gp to tolerate major structural changes.

Hepatic drug transporters: the journey so far

INTRODUCTION

The key role of transporter biology in both the manifestation and treatment of disease is now firmly established. Experiences of sub-optimal drug exposure due to drug-transporter interplay have supported incorporation of studies aimed at understanding the interactions between compounds and drug transporters much earlier in drug discovery. While drug transporters can impact the most pivotal pharmacokinetic parameter with respect to human dose and exposure projections, clearance, at a renal or hepatobiliary level, the latter will form the focus of this perspective.

AREAS COVERED

A synopsis of guidelines on which transporters to study together with an overview of the currently available toolkit is presented. A perspective on when to conduct studies with various hepatic transporters is also provided together with structural "alerts" which should prompt early investigation.

EXPERT OPINION

Great progress has been made in individual laboratories and via consortia to understand the role of drug transporters in disease, drug disposition, drug-drug interactions and toxicity. A systematic analysis of the value posed by the available approaches and an inter-lab comparison now seems warranted. The emerging ability to use physico-chemical properties to guide future screening cascades promises to revolutionise the efficiency of early drug discovery.

Drug-protein hydrogen bonds govern the inhibition of the ATP hydrolysis of the multidrug transporter P-glycoprotein

P-glycoprotein (P-gp) is a member of the ATP-binding cassette transporter superfamily. This multidrug transporter utilizes energy from ATP hydrolysis for the efflux of a variety of hydrophobic and amphipathic compounds including anticancer drugs. Most of the substrates and modulators of P-gp stimulate its basal ATPase activity, although some inhibit it. The molecular mechanisms that are in play in either case are unknown. In this report, mutagenesis and molecular modeling studies of P-gp led to the identification of a pair of phenylalanine-tyrosine structural motifs in the transmembrane region that mediate the inhibition of ATP hydrolysis by certain drugs (zosuquidar, elacridar and tariquidar), with high affinity (IC50's ranging from 10 to 30nM). Upon mutation of any of these residues, drugs that inhibit the ATPase activity of P-gp switch to stimulation of the activity. Molecular modeling revealed that the phenylalanine residues F978 and F728 interact with tyrosine residues Y953 and Y310, respectively, in an edge-to-face conformation, which orients the tyrosines in such a way that they establish hydrogen-bond contacts with the inhibitor. Biochemical investigations along with transport studies in intact cells showed that the inhibitors bind at a high affinity site to produce inhibition of ATP hydrolysis and transport function. Upon mutation, they bind at lower affinity sites, stimulating ATP hydrolysis and only poorly inhibiting transport. These results also reveal that screening chemical compounds for their ability to inhibit the basal ATP hydrolysis can be a reliable tool to identify modulators with high affinity for P-gp.

Nitric oxide releasing acridone carboxamide derivatives as reverters of doxorubicin resistance in MCF7/Dx cancer cells

A series of nitric oxide donating acridone derivatives are synthesized and evaluated for in vitro cytotoxic activity against different sensitive and resistant cancer cell lines MCF7/Wt, MCF7/Mr (BCRP overexpression) and MCF7/Dx (P-gp expression). The results showed that NO-donating acridones are potent against both the sensitive and resistant cells. Structure activity relationship indicate that the nitric oxide donating moiety connected through a butyl chain at N(10) position as well as morpholino moiety linkage through an amide bridge on the acridone ring system at C-2 position, are required to exert a good cytotoxic effect. Further, good correlations were observed when cytotoxic properties were compared with in vitro nitric oxide release rate, nitric oxide donating group potentiated the cytotoxic effect of the acridone derivatives. Exogenous release of nitric oxide by NO donating acridones enhanced the accumulation of doxorubicin in MCF7/Dx cell lines when it was coadministered with doxorubicin, which inhibited the efflux process of doxorubicin. In summary, a nitric oxide donating group can potentiate the anti-MDR property of acridones.

Computational modeling to predict the functions and impact of drug transporters

Transport proteins are important mediators of cellular drug influx and efflux and play crucial roles in drug distribution, disposition and clearance. Drug-drug interactions have increasingly been found to occur at the transporter level and, hence, computational tools for studying drug-transporter interactions have gained in interest. In this short review, we present the most important transport proteins for drug influx and efflux. Computational tools for predicting and understanding the substrate and inhibitor interactions with these membrane-bound proteins are discussed. We have primarily focused on ligand-based and structure-based modeling, for which the state-of-the-art and future challenges are also discussed.

In silico methods for predicting drug-drug interactions with cytochrome P-450s, transporters and beyond

Drug-drug interactions (DDIs) are associated with severe adverse effects that may lead to the patient requiring alternative therapeutics and could ultimately lead to drug withdrawal from the market if they are severe. To prevent the occurrence of DDI in the clinic, experimental systems to evaluate drug interaction have been integrated into the various stages of the drug discovery and development process. A large body of knowledge about DDI has also accumulated through these studies and pharmacovigillence systems. Much of this work to date has focused on the drug metabolizing enzymes such as cytochrome P-450s as well as drug transporters, ion channels and occasionally other proteins. This combined knowledge provides a foundation for a hypothesis-driven in silico approach, using either cheminformatics or physiologically based pharmacokinetics (PK) modeling methods to assess DDI potential. Here we review recent advances in these approaches with emphasis on hypothesis-driven mechanistic models for important protein targets involved in PK-based DDI. Recent efforts with other informatics approaches to detect DDI are highlighted. Besides DDI, we also briefly introduce drug interactions with other substances, such as Traditional Chinese Medicines to illustrate how in silico modeling can be useful in this domain. We also summarize valuable data sources and web-based tools that are available for DDI prediction. We finally explore the challenges we see faced by in silico approaches for predicting DDI and propose future directions to make these computational models more reliable, accurate, and publically accessible.

Prediction of drug-ABC-transporter interaction--Recent advances and future challenges

With the discovery of P-glycoprotein (P-gp), it became evident that ABC-transporters play a vital role in bioavailability and toxicity of drugs. They prevent intracellular accumulation of toxic compounds, which renders them a major defense mechanism against xenotoxic compounds. Their expression in cells of all major barriers (intestine, blood–brain barrier, blood–placenta barrier) as well as in metabolic organs (liver, kidney) also explains their influence on the ADMET properties of drugs and drug candidates. Thus, in silico models for the prediction of the probability of a compound to interact with P-gp or analogous transporters are of high value in the early phase of the drug discovery process. Within this review, we highlight recent developments in the area, with a special focus on the molecular basis of drug–transporter interaction. In addition, with the recent availability of X-ray structures of several ABC-transporters, also structure-based design methods have been applied and will be addressed.

Predicting the extent of metabolism using in vitro permeability rate measurements and in silico permeability rate predictions

The Biopharmaceutics Drug Disposition Classification System (BDDCS) can be utilized to predict drug disposition, including interactions with other drugs and transporter or metabolizing enzyme effects based on the extent of metabolism and solubility of a drug. However, defining the extent of metabolism relies upon clinical data. Drugs exhibiting high passive intestinal permeability rates are extensively metabolized. Therefore, we aimed to determine if in vitro measures of permeability rate or in silico permeability rate predictions could predict the extent of metabolism, to determine a reference compound representing the permeability rate above which compounds would be expected to be extensively metabolized, and to predict the major route of elimination of compounds in a two-tier approach utilizing permeability rate and a previously published model predicting the major route of elimination of parent drug. Twenty-two in vitro permeability rate measurement data sets in Caco-2 and MDCK cell lines and PAMPA were collected from the literature, while in silico permeability rate predictions were calculated using ADMET Predictor or VolSurf+. The potential for permeability rate to differentiate between extensively and poorly metabolized compounds was analyzed with receiver operating characteristic curves. Compounds that yielded the highest sensitivity-specificity average were selected as permeability rate reference standards. The major route of elimination of poorly permeable drugs was predicted by our previously published model, and the accuracies and predictive values were calculated. The areas under the receiver operating curves were >0.90 for in vitro measures of permeability rate and >0.80 for the VolSurf+ model of permeability rate, indicating they were able to predict the extent of metabolism of compounds. Labetalol and zidovudine predicted greater than 80% of extensively metabolized drugs correctly and greater than 80% of poorly metabolized drugs correctly in Caco-2 and MDCK, respectively, while theophylline predicted greater than 80% of extensively and poorly metabolized drugs correctly in PAMPA. A two-tier approach predicting elimination route predicts 72 ± 9%, 49 ± 10%, and 66 ± 7% of extensively metabolized, biliarily eliminated, and renally eliminated parent drugs correctly when the permeability rate is predicted in silico and 74 ± 7%, 85 ± 2%, and 73 ± 8% of extensively metabolized, biliarily eliminated, and renally eliminated parent drugs correctly when the permeability rate is determined in vitro.

A high throughput flow cytometric assay platform targeting transporter inhibition

This review highlights the concepts, recent applications and limitations of High Throughput Screening (HTS) flow cytometry-based efflux inhibitory assays. This platform has been employed in mammalian and yeast efflux systems leading to the identification of small molecules with transporter inhibitory capabilities. This technology offers the possibility of substrate multiplexing and may promote novel strategies targeting microbial efflux systems. This platform can generate a comprehensive dataset that may support efforts to map the interface between chemistry and transporter biology in a variety of pathogenic systems.

Molecular basis of the polyspecificity of P-glycoprotein (ABCB1): recent biochemical and structural studies

ABCB1 (P-glycoprotein/P-gp) is an ATP-binding cassette transporter well known for its association with multidrug resistance in cancer cells. Powered by the hydrolysis of ATP, it effluxes structurally diverse compounds. In this chapter, we discuss current views on the molecular basis of the substrate polyspecificity of P-gp. One of the features that accounts for this property is the structural flexibility observed in P-gp. Several X-ray crystal structures of mouse P-gp have been published recently in the absence of nucleotide, with and without bound inhibitors. All the structures are in an inward-facing conformation exhibiting different degrees of domain separation, thus revealing a highly flexible protein. Biochemical and biophysical studies also demonstrate this flexibility in mouse as well as human P-gp. Site-directed mutagenesis has revealed the existence of multiple transport-active binding sites in P-gp for a single substrate. Thus, drugs can bind at either primary or secondary sites. Biochemical, molecular modeling, and structure-activity relationship studies suggest a large, common drug-binding pocket with overlapping sites for different substrates. We propose that in addition to the structural flexibility, the molecular or chemical flexibility also contributes to the binding of substrates to multiple sites forming the basis of polyspecificity.

Molecular insight into conformational transmission of human P-glycoprotein

P-glycoprotein (P-gp), a kind of ATP-binding cassette transporter, can export candidates through a channel at the two transmembrane domains (TMDs) across the cell membranes using the energy released from ATP hydrolysis at the two nucleotide-binding domains (NBDs). Considerable evidence has indicated that human P-gp undergoes large-scale conformational changes to export a wide variety of anti-cancer drugs out of the cancer cells. However, molecular mechanism of the conformational transmission of human P-gp from the NBDs to the TMDs is still unclear. Herein, targeted molecular dynamics simulations were performed to explore the atomic detail of the conformational transmission of human P-gp. It is confirmed that the conformational transition from the inward- to outward-facing is initiated by the movement of the NBDs. It is found that the two NBDs move both on the two directions (x and y). The movement on the x direction leads to the closure of the NBDs, while the movement on the y direction adjusts the conformations of the NBDs to form the correct ATP binding pockets. Six key segments (KSs) protruding from the TMDs to interact with the NBDs are identified. The relative movement of the KSs along the y axis driven by the NBDs can be transmitted through α-helices to the rest of the TMDs, rendering the TMDs to open towards periplasm in the outward-facing conformation. Twenty eight key residue pairs are identified to participate in the interaction network that contributes to the conformational transmission from the NBDs to the TMDs of human P-gp. In addition, 9 key residues in each NBD are also identified. The studies have thus provided clear insight into the conformational transmission from the NBDs to the TMDs in human P-gp.