Protein contacts and ligand binding in the inward-facing model of human P-glycoprotein

Effects of Point Mutations on the Thermal Stability of the NBD1 Domain of hCFTR

Cystic fibrosis (CF) is an autosomal recessive genetic disease caused by mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) chloride channel. The first nucleotide-binding domain (NBD1) of the CFTR is considered to be a hotspot for CF-causing mutations, and some of these mutations compromise the domain's thermal stability as well as its interactions with other domains. The mechanisms by which such mutations exert their deleterious effects are important in the basic research of this complex disease as well as for the development of mutation-specific therapies. With this in mind, we studied two class-II, severe, CF-causing mutations, L467P and A559T, known to destabilize the domain by 19.3 and 10.7 °C, respectively, and to lead to a misfolded, nonfunctioning CFTR, by conducting microsecond-long molecular dynamics (MD) simulations at an elevated temperature of 410 K on L467P-NBD1 and A559T-NBD1 constructs. For comparison, similar simulations were also performed on the wild-type (WT) construct and on the 6SS-NBD1 and 2PT/M470V-NBD1 constructs, both bearing sets of stabilizing mutations that stabilize the domain by 17.5 and 8.2 °C, respectively. The resulting trajectories were analyzed using multiple metrics, leading to a good correlation between the experimental ΔTm values and the results of the simulations, as well as multiple experimental observations and results of previous modeling efforts. Specifically, our analyses point to specific regions within NBD1 that are substantially affected by the L467P and A559T mutations and, therefore, may play some role in their pathogenesis. Many of these regions are also known to be important for the proper folding and function of the full-length CFTR. Using time-dependent assignment of DSSP elements, we also found that the two mutants follow different disintegration pathways, that of L467P-NBD1 starting in region 464-471 which resides within the F1-like ATP-binding core subdomain and continues in regions 550-562 and 514-523 within the ABCα subdomain whereas that of A559T-NBD1 simultaneously starting at the 550-562 and 514-523 regions. We propose that the analyses presented in this work may pave the way toward the development of L467P and A559T-specific CF therapies and by extension to other mutation-specific therapies for CF and for other diseases involving mutations in NBDs of other proteins.

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.

Computational and artificial intelligence-based approaches for drug metabolism and transport prediction

Drug metabolism and transport, orchestrated by drug-metabolizing enzymes (DMEs) and drug transporters (DTs), are implicated in drug-drug interactions (DDIs) and adverse drug reactions (ADRs). Reliable and precise predictions of DDIs and ADRs are critical in the early stages of drug development to reduce the rate of drug candidate failure. A variety of experimental and computational technologies have been developed to predict DDIs and ADRs. Recent artificial intelligence (AI) approaches offer new opportunities for better predicting and understanding the complex processes related to drug metabolism and transport. We summarize the role of major DMEs and DTs, and provide an overview of current progress in computational approaches for the prediction of drug metabolism, transport, and DDIs, with an emphasis on AI including machine learning (ML) and deep learning (DL) modeling.

Probing the Allosteric Modulation of P-Glycoprotein: A Medicinal Chemistry Approach Toward the Identification of Noncompetitive P-Gp Inhibitors

A medicinal chemistry approach combining in silico and in vitro methodologies was performed aiming at identifying and characterizing putative allosteric drug-binding sites (aDBSs) at the interface of the transmembrane- and nucleotide-binding domains (TMD-NBD) of P-glycoprotein. Two aDBSs were identified, one in TMD1/NBD1 and another one in TMD2/NBD2, by means of in silico fragment-based molecular dynamics and characterized in terms of size, polarity, and lining residues. From a small library of thioxanthone and flavanone derivatives, experimentally described to bind at the TMD-NBD interfaces, several compounds were identified to be able to decrease the verapamil-stimulated ATPase activity. An IC₅₀ of 81 ± 6.6 μM is reported for a flavanone derivative in the ATPase assays, providing evidence for an allosteric efflux modulation in P-glycoprotein. Molecular docking and molecular dynamics gave additional insights on the binding mode on how flavanone derivatives may act as allosteric inhibitors.

Long-range communication between transmembrane- and nucleotide-binding domains does not depend on drug binding to mutant P-glycoprotein

In this study, the impact of four P-gp mutations (G185V, G830V, F978A and ΔF335) on drug-binding and efflux-related signal-transmission mechanism was comprehensively evaluated in the presence of ligands within the drug-binding pocket (DBP), experimentally related with changes in their drug efflux profiles. The severe repacking of the transmembrane helices (TMH), induced by mutations and exacerbated by the presence of ligands, indicates that P-gp is sensitive to perturbations in the transmembrane region. Alterations on drug-binding were also observed as a consequence of the TMH repacking, but were not always correlated with alterations on ligands binding mode and/or binding affinity. Finally, and although all P-gp variants holo systems showed considerable changes in the intracellular coupling helices/nucleotide-binding domain (ICH-NBD) interactions, they seem to be primarily induced by the mutation itself rather than by the presence of ligands within the DBP. The data further suggest that the changes in drug efflux experimentally reported are mostly related with changes on drug specificity rather than effects on signal-transmission mechanism. We also hypothesize that an increase in the drug-binding affinity may also be related with the decreased drug efflux, while minor changes in binding affinities are possibly related with the increased drug efflux observed in transfected cells.Communicated by Ramaswamy H. Sarma.

Structure of ABCB1/P-Glycoprotein in the Presence of the CFTR Potentiator Ivacaftor

ABCB1/P-glycoprotein is an ATP binding cassette transporter that is involved in the clearance of xenobiotics, and it affects the disposition of many drugs in the body. Conformational flexibility of the protein within the membrane is an intrinsic part of its mechanism of action, but this has made structural studies challenging. Here, we have studied different conformations of P-glycoprotein simultaneously in the presence of ivacaftor, a known competitive inhibitor. In order to conduct this, we used high contrast cryo-electron microscopy imaging with a Volta phase plate. We associate the presence of ivacaftor with the appearance of an additional density in one of the conformational states detected. The additional density is in the central aqueous cavity and is associated with a wider separation of the two halves of the transporter in the inward-facing state. Conformational changes to the nucleotide-binding domains are also observed and may help to explain the stimulation of ATPase activity that occurs when transported substrate is bound in many ATP binding cassette transporters.

N-Substituted piperazine derivatives as potential multitarget agents acting on histamine H3 receptor and cancer resistance proteins

Taking into account that multidrug resistance (MDR) is the main cause for chemotherapeutic failure in cancer treatment, the ability of novel histamine H₃ receptor ligands to reverse the cancer MDR was evaluated, using the ABCB1 efflux pump inhibition assay in mouse MDR T-lymphoma cells. The most active compounds displayed significant cytotoxic and antiproliferative effects as well as a very potent MDR efflux pump inhibitory action, 3-5-fold stronger than that of reference inhibitor verapamil. Although these compounds possess weak antagonistic properties against histamine H₃ receptors, they are valuable pharmacological tools in the search for novel anticancer molecules. Furthermore, for the most active compounds, an insight into mechanisms of action using either, the luminescent Pgp-Glo™ Assay in vitro or docking studies to human Pgp, was performed.

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.

Improvement of conventional anti-cancer drugs as new tools against multidrug resistant tumors

Multidrug resistance (MDR) is the dominant cause of the failure of cancer chemotherapy. The design of antitumor drugs that are able to evade MDR is rapidly evolving, showing that this area of biomedical research attracts great interest in the scientific community. The current review explores promising recent approaches that have been developed with the aim of circumventing or overcoming MDR. Encouraging results have been obtained in the investigation of the MDR-modulating properties of various classes of natural compounds and their analogues. Inhibition of P-gp or downregulation of its expression have proven to be the main mechanisms by which MDR can be surmounted. The use of hybrid molecules that are able to simultaneously interact with two or more cancer cell targets is currently being explored as a means to circumvent drug resistance. This strategy is based on the design of hybrid compounds that are obtained either by merging the structural features of separate drugs, or by conjugating two drugs or pharmacophores via cleavable/non-cleavable linkers. The approach is highly promising due to the pharmacokinetic and pharmacodynamic advantages that can be achieved over the independent administration of the two individual components. However, it should be stressed that the task of obtaining successful multivalent drugs is a very challenging one. The conjugation of anticancer agents with nitric oxide (NO) donors has recently been developed, creating a particular class of hybrid that can combat tumor drug resistance. Appropriate NO donors have been shown to reverse drug resistance via nitration of ABC transporters and by interfering with a number of metabolic enzymes and signaling pathways. In fact, hybrid compounds that are produced by covalently attaching NO-donors and antitumor drugs have been shown to elicit a synergistic cytotoxic effect in a variety of drug resistant cancer cell lines. Another strategy to circumvent MDR is based on nanocarrier-mediated transport and the controlled release of chemotherapeutic drugs and P-gp inhibitors. Their pharmacokinetics are governed by the nanoparticle or polymer carrier and make use of the enhanced permeation and retention (EPR) effect, which can increase selective delivery to cancer cells. These systems are usually internalized by cancer cells via endocytosis and accumulate in endosomes and lysosomes, thus preventing rapid efflux. Other modalities to combat MDR are described in this review, including the pharmaco-modulation of acridine, which is a well-known scaffold in the development of bioactive compounds, the use of natural compounds as means to reverse MDR, and the conjugation of anticancer drugs with carriers that target specific tumor-cell components. Finally, the outstanding potential of in silico structure-based methods as a means to evaluate the ability of antitumor drugs to interact with drug transporters is also highlighted in this review. Structure-based design methods, which utilize 3D structural data of proteins and their complexes with ligands, are the most effective of the in silico methods available, as they provide a prediction regarding the interaction between transport proteins and their substrates and inhibitors. The recently resolved X-ray structure of human P-gp can help predict the interaction sites of designed compounds, providing insight into their binding mode and directing possible rational modifications to prevent them from becoming P-gp drug substrates. In summary, although major efforts were invested in the search for new tools to combat drug resistant tumors, they all require further implementation and methodological development. Further investigation and progress in the abovementioned strategies will provide significant advances in the rational combat against cancer MDR.

Novel Heat Shock Protein 90 Inhibitors Suppress P-Glycoprotein Activity and Overcome Multidrug Resistance in Cancer Cells

Heat Shock Protein 90 (Hsp90) chaperone interacts with a broad range of client proteins involved in cancerogenesis and cancer progression. However, Hsp90 inhibitors were unsuccessful as anticancer agents due to their high toxicity, lack of selectivity against cancer cells and extrusion by membrane transporters responsible for multidrug resistance (MDR) such as P-glycoprotein (P-gp). Recognizing the potential of new compounds to inhibit P-gp function and/or expression is essential in the search for effective anticancer drugs. Eleven Hsp90 inhibitors containing an isoxazolonaphtoquinone core were synthesized and evaluated in two MDR models comprised of sensitive and corresponding resistant cancer cells with P-gp overexpression (human non-small cell lung carcinoma and colorectal adenocarcinoma). We investigated the effect of Hsp90 inhibitors on cell growth inhibition, P-gp activity and P-gp expression. Structure-activity relationship analysis was performed in respect to cell growth and P-gp inhibition. Compounds 5, 7, and 9 directly interacted with P-gp and inhibited its ATPase activity. Their potential P-gp binding site was identified by molecular docking studies. In addition, these compounds downregulated P-gp expression in MDR colorectal carcinoma cells, showed good relative selectivity towards cancer cells, while compound 5 reversed resistance to doxorubicin and paclitaxel in concentration-dependent manner. Therefore, compounds 5, 7 and 9 could be promising candidates for treating cancers with P-gp overexpression.

New insights in the in vitro characterisation and molecular modelling of the P-glycoprotein inhibitory promiscuity

The presence of several binding sites for both substrates and inhibitors is yet a poorly explored thematic concerning the assessment of the drug-drug interactions risk due to interactions of multiple drugs with the human transport protein P-glycoprotein (P-gp or MDR1, gene ABCB1). In this study we measured the inhibitory behaviour of a set of known drugs towards P-gp by using three different probe substrates (digoxin, Hoechst 33,342 and rhodamine 123). A structure-based model was built to unravel the different substrates binding sites and to rationalize the cases where drugs were not inhibiting all the substrates. A separate set of experiments was used to validate the model and confirmed its suitability to either detect the substrate-dependent P-gp inhibition and to anticipate proper substrates for in vitro experiments case by case. The modelling strategy described can be applied for either design safer drugs (P-gp as antitarget) or to target specific sub-site inhibitors towards other drugs (P-gp as target).

A short cross-linker activates human P-glycoprotein missing a catalytic carboxylate

P-glycoprotein (P-gp) is an ATP-dependent drug pump that protects us from toxic agents and confers multidrug resistance. It has a tweezer-like structure with each arm consisting of a transmembrane domain (TMD) and a nucleotide-binding domain (NBD). Drug substrates bind to sites within the TMDs to activate ATPase activity by promoting a tweezer-like closing of the gap between the NBDs. The catalytic carboxylates may be critical for NBD movements because the E556Q(NBD1) or E1201Q(NBD2) mutation inhibited drug-stimulated ATPase activity. If the catalytic carboxylates were components of the mechanism to bring the NBDs together, then we predicted that insertion of a flexible cross-linker between the arms would increase ATPase activity of the mutants. We found that cross-linking (between L175C(TMD1) and N820C(TMD2)) with a short flexible cross-linker (7.8Å maximum) restored high levels of drug-stimulated ATPase activity of the E556Q or E1201Q mutants. Cross-linking with a longer cross-linker (22Å maximum) however, did not restore activity. Cross-linking could not rescue all ATPase deficient mutants. For example, cross-linking L175C/N820C with short or long cross-linkers did not activate the H-loop mutants H587A or H1232A or the Walker A K433M or K1076M mutants. The results suggest that the E556 and E1201 catalytic carboxylates are part of a spring-like mechanism that is required to facilitate movements between the open and closed conformations of P-gp during ATP hydrolysis.

Homology modelling of human P-glycoprotein

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

Toward Understanding the Catalytic Mechanism of Human Paraoxonase 1: Site-Specific Mutagenesis at Position 192

Human paraoxonase 1 (h-PON1) is a serum enzyme that can hydrolyze a variety of substrates. The enzyme exhibits anti-inflammatory, anti-oxidative, anti-atherogenic, anti-diabetic, anti-microbial and organophosphate-hydrolyzing activities. Thus, h-PON1 is a strong candidate for the development of therapeutic intervention against a variety conditions in human. However, the crystal structure of h-PON1 is not solved and the molecular details of how the enzyme hydrolyzes different substrates are not clear yet. Understanding the catalytic mechanism(s) of h-PON1 is important in developing the enzyme for therapeutic use. Literature suggests that R/Q polymorphism at position 192 in h-PON1 dramatically modulates the substrate specificity of the enzyme. In order to understand the role of the amino acid residue at position 192 of h-PON1 in its various hydrolytic activities, site-specific mutagenesis at position 192 was done in this study. The mutant enzymes were produced using Escherichia coli expression system and their hydrolytic activities were compared against a panel of substrates. Molecular dynamics simulation studies were employed on selected recombinant h-PON1 (rh-PON1) mutants to understand the effect of amino acid substitutions at position 192 on the structural features of the active site of the enzyme. Our results suggest that, depending on the type of substrate, presence of a particular amino acid residue at position 192 differentially alters the micro-environment of the active site of the enzyme resulting in the engagement of different subsets of amino acid residues in the binding and the processing of substrates. The result advances our understanding of the catalytic mechanism of h-PON1.

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

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

Do adsorbed drugs onto P-glycoprotein influence its efflux capability?

The membrane biophysical aspects by which multidrug resistance (MDR) relate to the ABC transporter function still remain largely unknown. Notwithstanding the central role that efflux pumps like P-glycoprotein have in MDR onset, experimental studies classified additionally the lipid micro-environment where P-gp is inserted as a determinant for the increased efflux capability demonstrated in MDR cell lines. Recently, a nonlinear model for drug-membrane interactions showed that, upon drug adsorption, long-range mechanical alterations are predicted to affect the P-gp ATPase function at external drug concentrations of ∼10-100 μM. However, our results also show that drug adsorption may also occur at P-gp nucleotide-binding domains where conformational changes drive the efflux cycle. Thus, we assessed the effect of drug adsorption to both protein-water and lipid-water interfaces by means of molecular dynamics simulations. The results show that free energies of adsorption are lower for modulators in both lipid/water and protein/water interfaces. Important differences in drug-protein interactions, protein dynamics and membrane biophysical characteristics were observed between the different classes. Therefore, we hypothesize that drug adsorption to the protein and lipid-water interface accounts for a complex network of events that affect the ability of transporters to efflux drugs.

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.

Three- and four-class classification models for P-glycoprotein inhibitors using counter-propagation neural networks

P-glycoprotein (P-gp) is an ATP binding cassette (ABC) transporter that helps to protect several certain human organs from xenobiotic exposure. This efflux pump is also responsible for multi-drug resistance (MDR), an issue of the chemotherapy approach in the fight against cancer. Therefore, the discovery of P-gp inhibitors is considered one of the most popular strategies to reverse MDR in tumour cells and to improve therapeutic efficacy of commonly used cytotoxic drugs. Until now, several generations of P-gp inhibitors have been developed but they have largely failed in preclinical and clinical studies due to lack of selectivity, poor solubility and severe pharmacokinetic interactions. In this study, three models (SION, SIO, SIN) to classify specific 'true' P-gp inhibitors as well as three other models (CPBN, CPB1, CPN) to distinguish between P-gp inhibitors, CYP 3A inhibitors and co-inhibitors of these proteins with rather high accuracy values for the test set and the external set were generated based on counter-propagation neural networks (CPG-NN). Such three and four-class classification models helped provide more information about the bioactivities of compounds not only on one target (P-gp), but also on a combination of multiple targets (P-gp, CYP 3A).

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.

In vivo imaging of multidrug resistance using a third generation MDR1 inhibitor

Cellular up-regulation of multidrug resistance protein 1 (MDR1) is a common cause for resistance to chemotherapy; development of third generation MDR1 inhibitors-several of which contain a common 6,7-dimethoxy-2-phenethyl-1,2,3,4-tetrahydroisoquinoline substructure-is underway. Efficacy of these agents has been difficult to ascertain, partly due to a lack of pharmacokinetic reporters for quantifying inhibitor localization and transport dynamics. Some of the recent third generation inhibitors have a pendant heterocycle, for example, a chromone moiety, which we hypothesized could be converted to a fluorophore. Following synthesis and teasing of a small set of analogues, we identified one lead compound that can be used as a cellular imaging agent that exhibits structural similarity and behavior akin to the latest generation of MDR1 inhibitors.

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.

Mutations in intracellular loops 1 and 3 lead to misfolding of human P-glycoprotein (ABCB1) that can be rescued by cyclosporine A, which reduces its association with chaperone Hsp70

P-glycoprotein (P-gp) is an ATP binding cassette transporter that effluxes a variety of structurally diverse compounds including anticancer drugs. Computational models of human P-gp in the apo- and nucleotide-bound conformation show that the adenine group of ATP forms hydrogen bonds with the conserved Asp-164 and Asp-805 in intracellular loops 1 and 3, respectively, which are located at the interface between the nucleotide binding domains and transmembrane domains. We investigated the role of Asp-164 and Asp-805 residues by substituting them with cysteine in a cysteine-less background. It was observed that the D164C/D805C mutant, when expressed in HeLa cells, led to misprocessing of P-gp, which thus failed to transport the drug substrates. The misfolded protein could be rescued to the cell surface by growing the cells at a lower temperature (27 °C) or by treatment with substrates (cyclosporine A, FK506), modulators (tariquidar), or small corrector molecules. We also show that short term (4-6 h) treatment with 15 μM cyclosporine A or FK506 rescues the pre-formed immature protein trapped in the endoplasmic reticulum in an immunophilin-independent pathway. The intracellularly trapped misprocessed protein associates more with chaperone Hsp70, and the treatment with cyclosporine A reduces the association of mutant P-gp, thus allowing it to be trafficked to the cell surface. The function of rescued cell surface mutant P-gp is similar to that of wild-type protein. These data demonstrate that the Asp-164 and Asp-805 residues are not important for ATP binding, as proposed earlier, but are critical for proper folding and maturation of a functional transporter.

Interactions of the multidrug resistance modulators tariquidar and elacridar and their analogues with P-glycoprotein

Tariquidar and elacridar are among the most potent inhibitors of the multidrug resistance transporter P-glycoprotein (P-gp), but how they interact with the protein is yet unknown. In this work, we describe a possible way in which these inhibitors interact with P-gp. We rely on structure-activity relationship analysis of a small group of tariquidar and elacridar analogues that was purposefully selected, designed, and tested. Structural modifications of the compounds relate to the presence or absence of functional groups in the tariquidar and elacridar scaffolds. The activity of the compounds was evaluated by their effects on the accumulation of P-gp substrates rhodamine 123 and Hoechst 33342 in resistant tumor cells. The data allow estimation of the ability of the compounds to interact with the experimentally proposed R- and H-sites to which rhodamine 123 and Hoechst 33342 bind, respectively. Using an inward-facing homology model of human P-gp based on the crystallographic structure of mouse P-gp, we demonstrate that these binding sites may overlap with the binding sites of the QZ59 ligands co-crystallized with mouse P-gp. Based on this SAR analysis, and using flexible alignment and docking, we propose possible binding modes for tariquidar and elacridar. Our results suggest the possibility for the studied compounds to bind to sites that coincide or overlap with the binding sites of rhodamine 123 and Hoechst 33342. These results contribute to further understanding of structure-function relationships of P-gp and can help in the design of selective and potent P-gp inhibitors with potential clinical use.