Structure-activity relationships of flavonoids as inhibitors of breast cancer resistance protein (BCRP).

Flavonoids are an interesting group of natural products ubiquitously present in human diet. Their consumption has been associated with various and differing beneficial health effects. However, several flavonoids have been reported to inhibit the breast cancer resistance protein (BCRP) encoded by the ABCG2 gene. Thus, the consumption of flavonoids with high inhibitory activity could change pharmacokinetics and drug levels of drugs that are BCRP substrates. In cancer patients receiving chemotherapy an increased intake of such flavonoids could lead to adverse effects. We investigated a structurally diverse set of flavonoids, including derivatives with a rare C-methylated structure that were isolated from plants used in traditional medicine. The flavones retusin and ayanin were found to be highly potent inhibitors of BCRP, showing only slightly less potency than Ko143, the most potent ABCG2 inhibitor known so far. The activity data were analyzed by 2D and 3D QSAR analyses and the results revealed the impact of the different substituents at the various positions of the flavonoid core on activity. Additionally, a lateral 2D QSAR analysis of data collected from the literature was performed aiming to derive more general information about the influence of distinct structural features on the inhibitory potency of flavonoids. The comparative QSAR analyses led to a consistent picture of the effects of the different substituents at various positions of the flavone backbone. The following structural features were found to contribute positively to BCRP inhibition: a hydroxyl group in position 5, double bond between position 2 and 3, and a methoxy group in position 3. The exchange of a 3-methoxy group by an OH-group acting also as a hydrogen bond donor, resulted in decrease in activity underlining the potential role of the hydrogen bond acceptor 3-OCH(3) for the interaction with BCRP.

Synthesis and biological evaluation of a small molecule library of 3rd generation multidrug resistance modulators.

The development of new modulators possessing high efficacy, low toxicity and high selectivity is a pivotal approach to overcoming P-glycoprotein (P-gp) mediated multidrug resistance (MDR) in tumour cells. In this study 39 compounds are presented which have been synthesized and pharmacologically investigated in our laboratory. Similarly to the potent 3rd generation MDR modulator tariquidar (XR9576) the compounds contain a tetrahydroisoquinoline-ethyl-phenylamine substructure that, in contrast to XR9576, is connected to a smaller hydrophobic part, thus leading to molecules of lower molecular weight. The connection between the tetrahydroisoquinoline-ethyl-phenylamine substructure and the hydrophobic part was achieved through four different types of linkers: amide, urea, amide-ether and amide-styryl. A number of structural modifications in the hydrophobic part were created. The calcein AM assay served as test system to determine the P-gp transport inhibitory potencies of the compounds. For the amide linker derivatives a structure-activity relationship analysis was performed outlining which structural modifications contributed to the inhibitory potency. The compounds containing a bicyclic hydrophobic part with a particular substituent in a specific orientation were identified as the most potent amide derivatives. Among the urea derivatives the compounds with highest inhibitory potency possessed an ortho-nitro substituent. The conformational analysis revealed that this position enables the formation of a hydrogen bond to the urea linker thus stabilizing the conformation. Regarding the amide-styryl derivatives the elongation of the amide linker seemed to be most decisive for the observed increase in activity. The most promising candidate in the whole library possess an amide-ether linker and an ortho-nitro substituent in the hydrophobic part. This compound inhibites P-gp slightly less than tariquidar and can serve as a lead structure for new potent P-gp modulators.

Functional assay and structure-activity relationships of new third-generation P-glycoprotein inhibitors.

Twenty-eight compounds, including 24 structurally related derivatives of tariquidar synthesized in our laboratory, and four XR compounds, reported by Xenova group Ltd, were investigated by the Hoechst 33342 and Calcein AM functional assays for estimation of their inhibitory effects on the transport activity of P-glycoprotein (P-gp). A high correlation between the effects obtained in both assays was observed at the substrate concentrations used. The analyses of kinetics data from experiments at different substrate concentrations revealed non-competitive inhibition in the Calcein AM assay and competitive inhibition in the Hoechst 33342 assay. The 3D structures of the compounds were further aligned on Hoechst 33342 using flexible and pharmacophore alignments. The results suggested that inhibitors could interact with the H-binding site of P-gp and this could potentially be achieved by different ways of binding. The best 3D-QSAR models, generated by CoMFA and CoMSIA, yielded an internal predictive squared correlation coefficient higher than 0.8 and included electrostatic, steric, hydrogen bond acceptor, and hydrophobic fields. Validation of the models on an external test set of 30 XR compounds gave predictive squared correlation coefficients of up to 0.66. An excellent correspondence between the experimental and modeled activities of the test compounds was observed. The models can be used for prediction and rational design of new P-gp inhibitors.

Pharmacophore model of drugs involved in P-glycoprotein multidrug resistance: explanation of structural variety (hypothesis).

A general pharmacophore model of P-glycoprotein (P-gp) drugs is proposed that is based on a highly diverse data set and relates to the verapamil binding site of the protein. It is derived from structurally different drugs using the program GASP. The pharmacophore model consists of two hydrophobic points, three hydrogen bond (HB) acceptor points, and one HB donor point. Pharmacophore patterns of various drugs are obtained, and different binding modes are presumed for some of them. It is concluded that the binding affinity of the drugs depends on the number of the pharmacophore points simultaneously involved in the interaction with P-gp. On the basis of the obtained results, a hypothesis is proposed to explain the broad structural variety of the P-gp substrates and inhibitors: (i) the verapamil binding site of P-gp has several points that can participate in hydrophobic and HB interactions; (ii) different drugs can interact with different receptor points in different binding modes.

Structure-function relationships of multidrug resistance P-glycoprotein.

The direct structure-function relationships of P-glycoprotein (P-gp) are presently unknown. In this paper two P-gp models are described: a homology model based on the Escherichia coli MsbA lipid transporter and a model based on the cross-linking results of Loo and Clarke. The pharmacophore pattern for the H-site (Hoechst 33342) is derived and binding sites on the transmembrane domains TM5 and TM11 are identified. Binding sites of rhodamines are also proposed on TM6 and TM12 in accordance with the published data. Location of the binding sites is opposite in both models, suggesting that TMs undergo rotation exposing the substrate bound from the membrane to the pore. It has been concluded that the models derived represent two different functional states of P-gp corresponding to nucleotide-free and nucleotide-bound P-gp. A qualitative correspondence to the P-gp crystallographic structure at 20 A resolution is found. A hypothesis is proposed about rearrangement of TMs upon state transition.

HAGE, the helicase antigen as a biomarker for breast cancer prognosis (WO2013144616).

INTRODUCTION: Breast cancer (BC) is the most common cancer in women and it ranks second as a cause of cancer death in women (after lung cancer). The receptor-based diagnosis of BC tumors allows application of more individual therapies. Depending on the status of the receptors and other risk markers, like tumor size and lymph node status, patients are assigned to risk classes. Invention of new biomarkers that could improve diagnosis and prognosis of BC patients is thus of an increased need. AREAS COVERED: The invention estimates the possibility of using the amount of expression of helicase antigen (HAGE) in samples of BC tissue as a biomarker for screening and prognosis. A total of 1650 BC patients were tested for HAGE expression and analyzed for well-characterized prognostic and predictive factors. HAGE expression was found to correlate significantly with aggressive clinicopathological features. A total of 443 triple negative patients were analyzed for therapeutic treatment received and survival. The HAGE expression was identified as a predictor of response to anthracycline treatment. The result was confirmed by comparison with an independent validation group. EXPERT OPINION: Тo identify patients who could be treated, a more detailed analysis taking into account, the full distribution spectrum of HAGE expression is desirable.

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

The primary aim of this work was to analyze the contacts between residues in the nucleotide binding domains (NBDs) and at the interface between the transmembrane domains (TMDs) and the NBDs in the inward-open homology model of human P-glycoprotein (P-gp). The analysis revealed communication nets through hydrogen bonding in the NBD and at the NBD-TMD interface of each half involving residues from the adenosine triphosphate (ATP) motifs and the coupling helices of the intracellular loops. Similar networks have been identified in P-gp conformations generated by molecular dynamics simulation. Differences have been recorded in the networking between both halves of P-gp. Many of the residue contacts have also been observed in the X-ray crystal structures of other ATP binding cassette (ABC) transporters, which confirms their validity. Next, possible binding pockets involving residues of importance for the TMD-NBD communication were identified. By studying these pockets, binding sites were suggested for rhodamine 123 (R-site) and prazosin (regulatory site) at the NBD-TMD interface that agreed with the experimental data on their location. Additionally, one more R-site in the protein cavity was proposed, in accordance with the available biochemical data. Together with the previously suggested Hoechst 33342 site (H-site), all sites were interpreted with respect to their effects on the protein ATPase activity, in correspondence with the experimental observations. Several residues involved in key contacts in the P-gp NBDs were proposed for further targeted mutagenesis experiments.

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.

Structure-activity relationships of tariquidar analogs as multidrug resistance modulators.

The review summarizes the most recent achievements in structure-activity relationship (SAR) studies of tariquidar and its analogs. Tariquidar is one of the most promising representatives of the third generation of multidrug resistance (MDR) modulators created so far. This fact determines the strong interest of different research groups in the development of tariquidar-like structures as selective inhibitors of MDR transporters in resistant human cancer cells. After the discovery of tariquidar, a number of analogs have been synthesized and pharmacologically tested, thus supplying good data for comprehensive analyses of their structure-activity relationships. In the review, the structural and pharmacological data of newly synthesized tariquidar-like compounds are first presented. Next, the main achievements in the SAR studies are described focusing on two main transport proteins: P-glycoprotein and breast cancer resistance protein. The reported results are discussed from the point of view of their significance and importance for future directions in the rational design of effective MDR modulators.

Comparison of the inward- and outward-open homology models and ligand binding of human P-glycoprotein.

An homology model of human P-glycoprotein, based on the X-ray structure of the recently resolved mouse P-glycoprotein, is presented. The model corresponds to the inward-facing conformation competent for drug binding. From the model, the residues involved in the protein-binding cavity are identified and compared with those in the outward-facing conformation of human P-glycoprotein developed previously based on the Sav1866 structure. A detailed analysis of the interactions of the cyclic peptides QZ59-RRR and QZ59-SSS is presented in both the X-ray structures of mouse P-glycoprotein and the human P-glycoprotein model generated by ligand docking. The results confirm the functional role of transmembrane domains TM4, TM6, TM10 and TM12 as entrance gates to the protein cavity, and also imply differences in their functions. The analysis of the cavities in both models suggests that the ligands remain bound to the same residues during the transition from the inward- to the outward-facing conformations. The analysis of the ligand-protein interactions in the X-ray complexes shows differences in the residues involved, as well as in the specific interactions performed by the same ligand within the same protein. This observation is supported by docking of the QZ59 ligands into human P-glycoprotein, thus aiding in the understanding of the complex behavior of P-glycoprotein substrates and inhibitors. The results confirm the possibility for multispecific drug interactions of the protein, and are important for elucidating the P-glycoprotein function and ligand interactions.

Combined pharmacophore modeling, docking, and 3D QSAR studies of ABCB1 and ABCC1 transporter inhibitors.

Quinazolinones, indolo- and pyrrolopyrimidines with inhibitory effects toward ABCB1 (P-gp) and ABCC1 (MRP1) transporters were studied by pharmacophore modeling, docking, and 3D QSAR to describe the binding preferences of the proteins. The pharmacophore overlays between dual and/or highly selective inhibitors point to binding sites of different topology and physiochemical properties for MRP1 and P-gp. Docking of selective inhibitors into the P-gp binding cavity by the use of a structural model based on the recently resolved P-gp structure confirms the P-gp pharmacophore features identified, and reveals the interactions of some functional groups and atoms in the structures with particular protein residues. The 3D QSAR analysis of the dual-effect inhibitors allows satisfactory prediction of the selectivity index of the compounds and outlines electrostatics as most important for selectivity. The results from the combined modeling approach complement each other and could improve our understanding of the protein-ligand interactions involved, and could aid in the development of highly selective and potent inhibitors of P-gp and MRP1.

Identification of putative binding sites of P-glycoprotein based on its homology model.

A homology model of P-glycoprotein based on the crystal structure of the multidrug transporter Sav1866 is developed, incorporated into a membrane environment, and optimized. The resulting model is analyzed in relation to the functional state and potential binding sites. The comparison of modeled distances to distances reported in experimental studies between particular residues suggests that the model corresponds most closely to the first ATP hydrolysis step of the protein transport cycle. Comparison to the protein 3D structure confirms this suggestion. Using SiteID and Site Finder programs three membrane related binding regions are identified: a region at the interface between the membrane and cytosol and two regions located in the transmembrane domains. The regions contain binding pockets of different size, orientation, and amino acids. A binding pocket located inside the membrane cavity is also identified. The pockets are analyzed in relation to amino acids shown experimentally to influence the protein function. The results suggest that the protein has multiple binding sites and may bind and/or release substrates in multiple pathways.

New functional assay of P-glycoprotein activity using Hoechst 33342.

In this study we describe a simplified, HTS-capable functional assay for the multidrug resistance (MDR) transporter P-glycoprotein (P-gp) based on its substrate Hoechst 33342. The physicochemical properties of Hoechst 33342 and the enormous milieu dependency of its fluorescence intensity allowed performing the assay in a homogeneous manner. This new assay served as an effective tool to estimate the potency of 10 well recognized P-gp substrates and modulators. Further, the potency of these compounds was also estimated in the calcein AM assay. The Hoechst 33342 and calcein AM assays yielded significantly comparable results for all compounds tested. Principal component analysis (PCA) applied to literature data on inhibition of P-gp activity and our results obtained in the Hoechst 33342 and calcein AM assay indicated similarity of compared functional transport assays. However, no correlation could be detected between these functional assays and the ATPase activity assay.

Structure-activity relationships of a series of tariquidar analogs as multidrug resistance modulators.

Tariquidar (XR9576) analogs, modulators of cancer multidrug resistance (MDR), were subjected to QSAR and 3D-QSAR analyses. The structural features contributing to anti-MDR activity were identified by the Free-Wilson analysis and pharmacophore search using Hoechst 33342 as a template. 3D-QSAR CoMFA and CoMSIA models were derived and tested. The best models yielded an external predictivity of 0.66-0.75 squared correlation coefficient and outlined HB-acceptor, steric, and hydrophobic fields as the most important 3D properties. On the basis of the QSAR and 3D-QSAR analyses it was suggested that the strong inhibitory potency of the compounds studied is related to the presence of a bulky aromatic ring system with a 3rd positioned heteroatom toward the anthranilamide nucleus in the opposite end of the tetrahydroquinoline group. The results can help in directing the rational design of new generations of potent P-glycoprotein MDR modulators.